Haplotype inference based on Hidden Markov Models in the QTL-MAS 2010 multi-generational dataset

<p>Abstract</p> <p>Background</p> <p>We have previously demonstrated an approach for efficient computation of genotype probabilities, and more generally probabilities of allele inheritance in inbred as well as outbred populations. That work also included an extension for haplotype inference, or phasing, using Hidden Markov Models. Computational phasing of multi-thousand marker datasets has not become common as of yet. In this communication, we further investigate the method presented earlier for such problems, in a multi-generational dataset simulated for QTL detection.</p> <p>Results</p> <p>When analyzing the dataset simulated for the 14th QTLMAS workshop, the phasing produced showed zero deviations compared to original simulated phase in the founder generation. In total, 99.93% of all markers were correctly phased. 97.68% of the individuals were correct in all markers over all 5 simulated chromosomes. Results were produced over a weekend on a small computational cluster. The specific algorithmic adaptations needed for the Markov model training approach in order to reach convergence are described.</p> <p>Conclusions</p> <p>Our method provides efficient, near-perfect haplotype inference allowing the determination of completely phased genomes in dense pedigrees. These developments are of special value for applications where marker alleles are not corresponding directly to QTL alleles, thus necessitating tracking of allele origin, and in complex multi-generational crosses. The cnF2freq codebase, which is in a current state of active development, is available under a BSD-style license.</p

Isoforms of U1-70k control subunit dynamics in the human spliceosomal U1 snRNP

Most human protein-encoding genes contain multiple exons that are spliced together, frequently in alternative arrangements, by the spliceosome. It is established that U1 snRNP is an essential component of the spliceosome, in human consisting of RNA and ten proteins, several of which are post- translationally modified and exist as multiple isoforms. Unresolved and challenging to investigate are the effects of these post translational modifications on the dynamics, interactions and stability of the particle. Using mass spectrometry we investigate the composition and dynamics of the native human U1 snRNP and compare native and recombinant complexes to isolate the effects of various subunits and isoforms on the overall stability. Our data reveal differential incorporation of four protein isoforms and dynamic interactions of subunits U1-A, U1-C and Sm-B/B’. Results also show that unstructured post- ranslationally modified C-terminal tails are responsible for the dynamics of Sm-B/B’ and U1-C and that their interactions with the Sm core are controlled by binding to different U1-70k isoforms and their phosphorylation status in vivo. These results therefore provide the important functional link between proteomics and structure as well as insight into the dynamic quaternary structure of the native U1 snRNP important for its function.This work was funded by: BBSRC (OVM), BBSRC and EPSRC (HH and NM), EU Prospects (HH), European Science Foundation (NM), the Royal Society (CVR), and fellowship from JSPS and HFSP (YM and DAPK respectively)

Quantitative trait loci conferring grain mineral nutrient concentrations in durum wheat 3 wild emmer wheat RIL population

Mineral nutrient malnutrition, and particularly deficiency in zinc and iron, afflicts over 3 billion people worldwide. Wild emmer wheat, Triticum turgidum ssp. dicoccoides, genepool harbors a rich allelic repertoire for mineral nutrients in the grain. The genetic and physiological basis of grain protein, micronutrients (zinc, iron, copper and manganese) and macronutrients (calcium, magnesium, potassium, phosphorus and sulfur) concentration was studied in tetraploid wheat population of 152 recombinant inbred lines (RILs), derived from a cross between durum wheat (cv. Langdon) and wild emmer (accession G18-16). Wide genetic variation was found among the RILs for all grain minerals, with considerable transgressive effect. A total of 82 QTLs were mapped for 10 minerals with LOD score range of 3.2–16.7. Most QTLs were in favor of the wild allele (50 QTLs). Fourteen pairs of QTLs for the same trait were mapped to seemingly homoeologous positions, reflecting synteny between the A and B genomes. Significant positive correlation was found between grain protein concentration (GPC), Zn, Fe and Cu, which was supported by significant overlap between the respective QTLs, suggesting common physiological and/or genetic factors controlling the concentrations of these mineral nutrients. Few genomic regions (chromosomes 2A, 5A, 6B and 7A) were found to harbor clusters of QTLs for GPC and other nutrients. These identified QTLs may facilitate the use of wild alleles for improving grain nutritional quality of elite wheat cultivars, especially in terms of protein, Zn and Fe

Antibodies to citrullinated proteins and differences in clinical progression of rheumatoid arthritis

Antibodies to citrullinated proteins (anti-cyclic-citrullinated peptide [anti-CCP] antibodies) are highly specific for rheumatoid arthritis (RA) and precede the onset of disease symptoms, indicating a pathogenetic role for these antibodies in RA. We recently showed that distinct genetic risk factors are associated with either anti-CCP-positive disease or anti-CCP-negative disease. These data are important as they indicate that distinct pathogenic mechanisms are underlying anti-CCP-positive disease or anti-CCP-negative disease. Likewise, these observations raise the question of whether anti-CCP-positive RA and anti-CCP-negative RA are clinically different disease entities. We therefore investigated whether RA patients with anti-CCP antibodies have a different clinical presentation and disease course compared with patients without these autoantibodies. In a cohort of 454 incident patients with RA, 228 patients were anti-CCP-positive and 226 patients were anti-CCP-negative. The early symptoms, tender and swollen joint count, and C-reactive protein level at inclusion, as well as the swollen joint count and radiological destruction during 4 years of follow-up, were compared for the two groups. There were no differences in morning stiffness, type, location and distribution of early symptoms, patients' rated disease activity and C-reactive protein at inclusion between RA patients with and without anti-CCP antibodies. The mean tender and swollen joint count for the different joints at inclusion was similar. At follow-up, patients with anti-CCP antibodies had more swollen joints and more severe radiological destruction. Nevertheless, the distribution of affected joints, for swelling, bone erosions and joint space narrowing, was similar. In conclusion, the phenotype of RA patients with or without anti-CCP antibodies is similar with respect to clinical presentation but differs with respect to disease course

Tissue Harvester with Functional Valve (THFV): Shidham's device for reproducibly higher specimen yield by fine needle aspiration biopsy with easy to perform steps

BACKGROUND: Fine needle aspiration biopsy (FNAB) cytology has been a highly effective methodology for tissue diagnosis and for various ancillary studies including molecular tests. In addition to other benefits, FNAB predominantly retrieves the diagnostic loosely cohesive cells in the lesion as compared to the adjacent supporting stroma with relatively higher cohesiveness. However, FNAB procedure performed with currently available resources is highly skill dependent with inter-performer variability, which compromises its full potential as a diagnostic tool. In this study we report a device overcoming these limitations. METHODS: 'Tissue Harvester with Functional Valve' (THFV) was evaluated as part of a phase 1 National Institute of Health (NIH) research grant under Small Business Technology Transfer (STTR) Program. Working prototypes of the device were prepared. Each of the four cytopathologists with previous cytopathology fellowship training and experience in performing FNAB evaluated 5 THFV and 5 hypodermic needles resulting in 40 specimens (20 with THFV, 20 with hypodermic needles). A piece of fresh cattle liver stuffed in latex glove was used as the specimen. Based on these results a finished design was finalized. RESULTS: The smears and cell blocks prepared from the specimens obtained by THFV were superior in terms of cellularity to specimens obtained with hypodermic needles. The tissuecrit of specimens obtained with THFV ranged from 70 to 100 μl (mean 87, SD 10), compared to 17 to 30 μl (mean 24, SD 4) with conventional hypodermic needles (p < .0001, Student t-test). The technical ease [on a scale of 1 (easy) to 5 (difficult)] with THFV ranged from 1 to 2 as compared to 2 to 3 with hypodermic needles. CONCLUSION: The specimen yield with the new THFV was significantly higher when compared to hypodermic needles. Also, the FNAB procedure with THFV was relatively easier in comparison with hypodermic needles. The final version of Shidham's THFV device would improve the FNAB specimen yield by eliminating the skill factor. The increased specimen yield by this device would also facilitate wider application of FNAB specimens for various ancillary tests, including molecular tests

Strong interface-induced spin-orbit coupling in graphene on WS2

Interfacial interactions allow the electronic properties of graphene to be modified, as recently demonstrated by the appearance of satellite Dirac cones in the band structure of graphene on hexagonal boron nitride (hBN) substrates. Ongoing research strives to explore interfacial interactions in a broader class of materials in order to engineer targeted electronic properties. Here we show that at an interface with a tungsten disulfide (WS2) substrate, the strength of the spin-orbit interaction (SOI) in graphene is very strongly enhanced. The induced SOI leads to a pronounced low-temperature weak anti-localization (WAL) effect, from which we determine the spin-relaxation time. We find that spin-relaxation time in graphene is two-to-three orders of magnitude smaller on WS2 than on SiO2 or hBN, and that it is comparable to the intervalley scattering time. To interpret our findings we have performed first-principle electronic structure calculations, which both confirm that carriers in graphene-on-WS2 experience a strong SOI and allow us to extract a spin-dependent low-energy effective Hamiltonian. Our analysis further shows that the use of WS2 substrates opens a possible new route to access topological states of matter in graphene-based systems.Comment: Originally submitted version in compliance with editorial guidelines. Final version with expanded discussion of the relation between theory and experiments to be published in Nature Communication

Diagnosis of prostate cancer by detection of minichromosome maintenance 5 protein in urine sediments

Background: The accuracy of prostate-specific antigen (PSA) testing in prostate cancer detection is constrained by low sensitivity and specificity. Dysregulated expression of minichromosome maintenance (Mcm) 2–7 proteins is an early event in epithelial multistep carcinogenesis and thus MCM proteins represent powerful cancer diagnostic markers. In this study we investigate Mcm5 as a urinary biomarker for prostate cancer detection. Methods: Urine was obtained from 88 men with prostate cancer and from two control groups negative for malignancy. A strictly normal cohort included 28 men with complete, normal investigations, no urinary calculi and serum PSA <2 ng ml–1. An expanded control cohort comprised 331 men with a benign final diagnosis, regardless of PSA level. Urine was collected before and after prostate massage in the cancer patient cohort. An immunofluorometric assay was used to measure Mcm5 levels in urine sediments. Results: The Mcm5 test detected prostate cancer with 82% sensitivity (confidence interval (CI)= 72–89%) and with a specificity ranging from 73 (CI=68–78%) to 93% (CI=76–99%). Prostate massage led to increased Mcm5 signals compared with pre-massage samples (median 3440 (interquartile range (IQR) 2280 to 5220) vs 2360 (IQR <1800 to 4360); P=0.009), and was associated with significantly increased diagnostic sensitivity (82 vs 60%; P=0.012). Conclusions: Urinary Mcm5 detection seems to be a simple, accurate and noninvasive method for identifying patients with prostate cancer. Large-scale prospective trials are now required to evaluate this test in diagnosis and screening

Conducting High-Value Secondary Dataset Analysis: An Introductory Guide and Resources

Secondary analyses of large datasets provide a mechanism for researchers to address high impact questions that would otherwise be prohibitively expensive and time-consuming to study. This paper presents a guide to assist investigators interested in conducting secondary data analysis, including advice on the process of successful secondary data analysis as well as a brief summary of high-value datasets and online resources for researchers, including the SGIM dataset compendium (www.sgim.org/go/datasets). The same basic research principles that apply to primary data analysis apply to secondary data analysis, including the development of a clear and clinically relevant research question, study sample, appropriate measures, and a thoughtful analytic approach. A real-world case description illustrates key steps: (1) define your research topic and question; (2) select a dataset; (3) get to know your dataset; and (4) structure your analysis and presentation of findings in a way that is clinically meaningful. Secondary dataset analysis is a well-established methodology. Secondary analysis is particularly valuable for junior investigators, who have limited time and resources to demonstrate expertise and productivity

Gene-Boosted Assembly of a Novel Bacterial Genome from Very Short Reads

Recent improvements in technology have made DNA sequencing dramatically faster and more efficient than ever before. The new technologies produce highly accurate sequences, but one drawback is that the most efficient technology produces the shortest read lengths. Short-read sequencing has been applied successfully to resequence the human genome and those of other species but not to whole-genome sequencing of novel organisms. Here we describe the sequencing and assembly of a novel clinical isolate of Pseudomonas aeruginosa, strain PAb1, using very short read technology. From 8,627,900 reads, each 33 nucleotides in length, we assembled the genome into one scaffold of 76 ordered contiguous sequences containing 6,290,005 nucleotides, including one contig spanning 512,638 nucleotides, plus an additional 436 unordered contigs containing 416,897 nucleotides. Our method includes a novel gene-boosting algorithm that uses amino acid sequences from predicted proteins to build a better assembly. This study demonstrates the feasibility of very short read sequencing for the sequencing of bacterial genomes, particularly those for which a related species has been sequenced previously, and expands the potential application of this new technology to most known prokaryotic species

Observation of associated near-side and away-side long-range correlations in √sNN=5.02 TeV proton-lead collisions with the ATLAS detector

Two-particle correlations in relative azimuthal angle (Δϕ) and pseudorapidity (Δη) are measured in √sNN=5.02  TeV p+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using approximately 1  μb-1 of data as a function of transverse momentum (pT) and the transverse energy (ΣETPb) summed over 3.1<η<4.9 in the direction of the Pb beam. The correlation function, constructed from charged particles, exhibits a long-range (2<|Δη|<5) “near-side” (Δϕ∼0) correlation that grows rapidly with increasing ΣETPb. A long-range “away-side” (Δϕ∼π) correlation, obtained by subtracting the expected contributions from recoiling dijets and other sources estimated using events with small ΣETPb, is found to match the near-side correlation in magnitude, shape (in Δη and Δϕ) and ΣETPb dependence. The resultant Δϕ correlation is approximately symmetric about π/2, and is consistent with a dominant cos⁡2Δϕ modulation for all ΣETPb ranges and particle pT