Comprehensive comparative analysis of strand-specific RNA sequencing methods

Strand-specific, massively parallel cDNA sequencing (RNA-seq) is a powerful tool for transcript discovery, genome annotation and expression profiling. There are multiple published methods for strand-specific RNA-seq, but no consensus exists as to how to choose between them. Here we developed a comprehensive computational pipeline to compare library quality metrics from any RNA-seq method. Using the well-annotated Saccharomyces cerevisiae transcriptome as a benchmark, we compared seven library-construction protocols, including both published and our own methods. We found marked differences in strand specificity, library complexity, evenness and continuity of coverage, agreement with known annotations and accuracy for expression profiling. Weighing each method's performance and ease, we identified the dUTP second-strand marking and the Illumina RNA ligation methods as the leading protocols, with the former benefitting from the current availability of paired-end sequencing. Our analysis provides a comprehensive benchmark, and our computational pipeline is applicable for assessment of future protocols in other organisms.Howard Hughes Medical InstituteUnited States-Israel Binational Science Foundatio

Systematic discovery of unannotated genes in 11 yeast species using a database of orthologous genomic segments

<p>Abstract</p> <p>Background</p> <p>In standard BLAST searches, no information other than the sequences of the query and the database entries is considered. However, in situations where two genes from different species have only borderline similarity in a BLAST search, the discovery that the genes are located within a region of conserved gene order (synteny) can provide additional evidence that they are orthologs. Thus, for interpreting borderline search results, it would be useful to know whether the syntenic context of a database hit is similar to that of the query. This principle has often been used in investigations of particular genes or genomic regions, but to our knowledge it has never been implemented systematically.</p> <p>Results</p> <p>We made use of the synteny information contained in the Yeast Gene Order Browser database for 11 yeast species to carry out a systematic search for protein-coding genes that were overlooked in the original annotations of one or more yeast genomes but which are syntenic with their orthologs. Such genes tend to have been overlooked because they are short, highly divergent, or contain introns. The key features of our software - called SearchDOGS - are that the database entries are classified into sets of genomic segments that are already known to be orthologous, and that very weak BLAST hits are retained for further analysis if their genomic location is similar to that of the query. Using SearchDOGS we identified 595 additional protein-coding genes among the 11 yeast species, including two new genes in <it>Saccharomyces cerevisiae</it>. We found additional genes for the mating pheromone a-factor in six species including <it>Kluyveromyces lactis</it>.</p> <p>Conclusions</p> <p>SearchDOGS has proven highly successful for identifying overlooked genes in the yeast genomes. We anticipate that our approach can be adapted for study of further groups of species, such as bacterial genomes. More generally, the concept of doing sequence similarity searches against databases to which external information has been added may prove useful in other settings.</p

A pilot study to understand feasibility and acceptability of stool and cord blood sample collection for a large-scale longitudinal birth cohort

BACKGROUND: Few data are available to guide biological sample collection around the time of birth for large-scale birth cohorts. We are designing a large UK birth cohort to investigate the role of infection and the developing immune system in determining future health and disease. We undertook a pilot to develop methodology for the main study, gain practical experience of collecting samples, and understand the acceptability of sample collection to women in late pregnancy. METHODS: Between February-July 2014, we piloted the feasibility and acceptability of collecting maternal stool, baby stool and cord blood samples from participants recruited at prolonged pregnancy and planned pre-labour caesarean section clinics at University College London Hospital. Participating women were asked to complete acceptability questionnaires. RESULTS: Overall, 265 women were approached and 171 (65%) participated, with ≥1 sample collected from 113 women or their baby (66%). Women had a mean age of 34 years, were primarily of white ethnicity (130/166, 78%), and half were nulliparous (86/169, 51%). Women undergoing planned pre-labour caesarean section were more likely than those who delivered vaginally to provide ≥1 sample (98% vs 54%), but less likely to provide maternal stool (10% vs 43%). Pre-sample questionnaires were completed by 110/171 women (64%). Most women reported feeling comfortable with samples being collected from their baby (<10% uncomfortable), but were less comfortable about their own stool (19% uncomfortable) or a vaginal swab (24% uncomfortable). CONCLUSIONS: It is possible to collect a range of biological samples from women around the time of delivery, and this was acceptable for most women. These data inform study design and protocol development for large-scale birth cohorts

Establishing What Constitutes a Healthy Human Gut Microbiome: State of the Science, Regulatory Considerations, and Future Directions.

On December 17, 2018, the North American branch of the International Life Sciences Institute (ILSI North America) convened a workshop "Can We Begin to Define a Healthy Gut Microbiome Through Quantifiable Characteristics?" with &gt;40 invited academic, government, and industry experts in Washington, DC. The workshop objectives were to 1) develop a collective expert assessment of the state of the evidence on the human gut microbiome and associated human health benefits, 2) see if there was sufficient evidence to establish measurable gut microbiome characteristics that could serve as indicators of "health," 3) identify short- and long-term research needs to fully characterize healthy gut microbiome-host relationships, and 4) publish the findings. Conclusions were as follows: 1) mechanistic links of specific changes in gut microbiome structure with function or markers of human health are not yet established; 2) it is not established if dysbiosis is a cause, consequence, or both of changes in human gut epithelial function and disease; 3) microbiome communities are highly individualized, show a high degree of interindividual variation to perturbation, and tend to be stable over years; 4) the complexity of microbiome-host interactions requires a comprehensive, multidisciplinary research agenda to elucidate relationships between gut microbiome and host health; 5) biomarkers and/or surrogate indicators of host function and pathogenic processes based on the microbiome need to be determined and validated, along with normal ranges, using approaches similar to those used to establish biomarkers and/or surrogate indicators based on host metabolic phenotypes; 6) future studies measuring responses to an exposure or intervention need to combine validated microbiome-related biomarkers and/or surrogate indicators with multiomics characterization of the microbiome; and 7) because static genetic sampling misses important short- and long-term microbiome-related dynamic changes to host health, future studies must be powered to account for inter- and intraindividual variation and should use repeated measures within individuals

Antibiotic use and the risk of rheumatoid arthritis: a population-based case-control study

Background: Antibiotic-induced disturbances of the human microbiota have been implicated in the development of chronic autoimmune conditions. This study aimed to assess whether antibiotic use is associated with the onset of rheumatoid arthritis (RA). Methods: A nested case-control study was conducted utilising data from the primary care Clinical Practice Research Datalink (CPRD). Patients with an incident diagnosis of RA were identified (1995–2017). Each case was matched on age, gender, and general practice to ≥ 5 controls without RA. Conditional logistic regression was used to examine previous antibiotic prescriptions and RA onset after controlling for confounding factors. Results: We identified 22,677 cases of RA, matched to 90,013 controls, with a median follow-up of 10 years before RA diagnosis. The odds of developing RA were 60% higher in those exposed to antibiotics than in those not exposed (OR 1.60; 95% CI 1.51–1.68). A dose- or frequency-dependent association was observed between the number of previous antibiotic prescriptions and RA. All classes of antibiotics were associated with higher odds of RA, with bactericidal antibiotics carrying higher risk than bacteriostatic (45% vs. 31%). Those with antibiotic-treated upper respiratory tract (URT) infections were more likely to be RA cases. However, this was not observed for URT infections not treated with antibiotics. Antifungal (OR = 1.27; 95% CI 1.20–1.35) and antiviral (OR = 1.19; 95% CI 1.14–1.24) prescriptions were also associated with increased odds of RA. Conclusion: Antibiotic prescriptions are associated with a higher risk of RA. This may be due to microbiota disturbances or underlying infections driving risk. Further research is needed to explore these mechanisms

The Vital Role of Social Workers in Community Partnerships: The Alliance for Gay, Lesbian, Bisexual, Transgender and Questioning Youth

The account of The Alliance for Gay, Lesbian, Bisexual, Transgender, and Questioning (GLBTQ) Youth formation offers a model for developing com- munity-based partnerships. Based in a major urban area, this university-community collaboration was spearheaded by social workers who were responsible for its original conceptualization, for generating community support, and for eventual staffing, administration, direct service provision, and program evaluation design. This article presents the strategic development and evolution of this community- based service partnership, highlighting the roles of schools of social work, academics, and social work students in concert with community funders, practitioners and youth, in responding to the needs of a vulnerable population

Genomic variation and strain-specific functional adaptation in the human gut microbiome during early life

The human gut microbiome matures towards the adult composition during the first years of life and is implicated in early immune development. Here, we investigate the effects of microbial genomic diversity on gut microbiome development using integrated early childhood data sets collected in the DIABIMMUNE study in Finland, Estonia and Russian Karelia. We show that gut microbial diversity is associated with household location and linear growth of children. Single nucleotide polymorphism- and metagenomic assembly-based strain tracking revealed large and highly dynamic microbial pangenomes, especially in the genus Bacteroides, in which we identified evidence of variability deriving from Bacteroides-targeting bacteriophages. Our analyses revealed functional consequences of strain diversity; only 10% of Finnish infants harboured Bifidobacterium longum subsp. infantis, a subspecies specialized in human milk metabolism, whereas Russian infants commonly maintained a probiotic Bifidobacterium bifidum strain in infancy. Groups of bacteria contributing to diverse, characterized metabolic pathways converged to highly subject-specific configurations over the first two years of life. This longitudinal study extends the current view of early gut microbial community assembly based on strain-level genomic variation.Peer reviewe

Kinetic CRAC uncovers a role for Nab3 in determining gene expression profiles during stress

RNA-binding proteins play a key role in shaping gene expression profiles during stress, however, little is known about the dynamic nature of these interactions and how this influences the kinetics of gene expression. To address this, we developed kinetic cross-linking and analysis of cDNAs (\u3c7CRAC), an ultraviolet cross-linking method that enabled us to quantitatively measure the dynamics of protein\u2013RNA interactions in vivo on a minute time-scale. Here, using \u3c7CRAC we measure the global RNA-binding dynamics of the yeast transcription termination factor Nab3 in response to glucose starvation. These measurements reveal rapid changes in protein\u2013RNA interactions within 1\u2009min following stress imposition. Changes in Nab3 binding are largely independent of alterations in transcription rate during the early stages of stress response, indicating orthogonal transcriptional control mechanisms. We also uncover a function for Nab3 in dampening expression of stress-responsive genes. \u3c7CRAC has the potential to greatly enhance our understanding of in vivo dynamics of protein\u2013RNA interactions

Noisy Splicing Drives mRNA Isoform Diversity in Human Cells

While the majority of multiexonic human genes show some evidence of alternative splicing, it is unclear what fraction of observed splice forms is functionally relevant. In this study, we examine the extent of alternative splicing in human cells using deep RNA sequencing and de novo identification of splice junctions. We demonstrate the existence of a large class of low abundance isoforms, encompassing approximately 150,000 previously unannotated splice junctions in our data. Newly-identified splice sites show little evidence of evolutionary conservation, suggesting that the majority are due to erroneous splice site choice. We show that sequence motifs involved in the recognition of exons are enriched in the vicinity of unconserved splice sites. We estimate that the average intron has a splicing error rate of approximately 0.7% and show that introns in highly expressed genes are spliced more accurately, likely due to their shorter length. These results implicate noisy splicing as an important property of genome evolution

Spliced Leader Trapping Reveals Widespread Alternative Splicing Patterns in the Highly Dynamic Transcriptome of Trypanosoma brucei

Trans-splicing of leader sequences onto the 5′ends of mRNAs is a widespread phenomenon in protozoa, nematodes and some chordates. Using parallel sequencing we have developed a method to simultaneously map 5′splice sites and analyze the corresponding gene expression profile, that we term spliced leader trapping (SLT). The method can be applied to any organism with a sequenced genome and trans-splicing of a conserved leader sequence. We analyzed the expression profiles and splicing patterns of bloodstream and insect forms of the parasite Trypanosoma brucei. We detected the 5′ splice sites of 85% of the annotated protein-coding genes and, contrary to previous reports, found up to 40% of transcripts to be differentially expressed. Furthermore, we discovered more than 2500 alternative splicing events, many of which appear to be stage-regulated. Based on our findings we hypothesize that alternatively spliced transcripts present a new means of regulating gene expression and could potentially contribute to protein diversity in the parasite. The entire dataset can be accessed online at TriTrypDB or through: http://splicer.unibe.ch/