A clinically validated human capillary blood transcriptome test for global systems biology studies.

To prevent and treat chronic diseases, including cancer, a global application of systems biology is needed. We report here a whole blood transcriptome test that needs only 50 µl of capillary (fingerprick) blood. This test is suitable for global applications because the samples are preserved at ambient temperature for up to 4 weeks and the RNA preservative inactivates all pathogens, enabling safe transportation. Both the laboratory and bioinformatic steps are automated and performed in a clinical lab, which minimizes batch effects and creates unbiased datasets. Given its clinical testing performance and accessibility to traditionally underrepresented and diverse populations, this test offers a unique ability to reveal molecular mechanisms of disease and enable longitudinal, population-scale studies.

Correction to: Implications of sequence variation on the evolution of rRNA.

The article Implications of sequence variation on the evolution of rRNA, written by Matthew M. Parks, Chad M. Kurylo, Jake E. Batchelder, C. Theresa Vincent and Scott C. Blanchard, was originally published electronically on the publisher's internet portal (currently SpringerLink).

Ribosome biogenesis during cell cycle arrest fuels EMT in development and disease.

Ribosome biogenesis is a canonical hallmark of cell growth and proliferation. Here we show that execution of Epithelial-to-Mesenchymal Transition (EMT), a migratory cellular program associated with development and tumor metastasis, is fueled by upregulation of ribosome biogenesis during G1/S arrest. This unexpected EMT feature is independent of species and initiating signal, and is accompanied by release of the repressive nucleolar chromatin remodeling complex (NoRC) from rDNA, together with recruitment of the EMT-driving transcription factor Snai1 (Snail1), RNA Polymerase I (Pol I) and the Upstream Binding Factor (UBF). EMT-associated ribosome biogenesis is also coincident with increased nucleolar recruitment of Rictor, an essential component of the EMT-promoting mammalian target of rapamycin complex 2 (mTORC2). Inhibition of rRNA synthesis in vivo differentiates primary tumors to a benign, Estrogen Receptor-alpha (ERα) positive, Rictor-negative phenotype and reduces metastasis. These findings implicate the EMT-associated ribosome biogenesis program with cellular plasticity, de-differentiation, cancer progression and metastatic disease.

Implications of sequence variation on the evolution of rRNA.

The evolution of the multi-copy family of ribosomal RNA (rRNA) genes is unique in regard to its genetics and genome evolution. Paradoxically, rRNA genes are highly homogenized within and between individuals, yet they are globally distinct between species. Here, we discuss the implications for models of rRNA gene evolution in light of our recent discoveries that ribosomes bearing rRNA sequence variants can affect gene expression and physiology and that intra-individual rRNA alleles exhibit both context- and tissue-specific expression.

ADAM10 controls the differentiation of the coronary arterial endothelium.

The coronary vasculature is crucial for normal heart function, yet much remains to be learned about its development, especially the maturation of coronary arterial endothelium. Here, we show that endothelial inactivation of ADAM10, a key regulator of Notch signaling, leads to defects in coronary arterial differentiation, as evidenced by dysregulated genes related to Notch signaling and arterial identity. Moreover, transcriptome analysis indicated reduced EGFR signaling in A10ΔEC coronary endothelium. Further analysis revealed that A10ΔEC mice have enlarged dysfunctional hearts with abnormal myocardial compaction, and increased expression of venous and immature endothelium markers. These findings provide the first evidence for a potential role for endothelial ADAM10 in cardioprotective homeostatic EGFR signaling and implicate ADAM10/Notch signaling in coronary arterial cell specification, which is vital for normal heart development and function. The ADAM10/Notch signaling pathway thus emerges as a potential therapeutic target for improving the regenerative capacity and maturation of the coronary vasculature.

Endogenous rRNA Sequence Variation Can Regulate Stress Response Gene Expression and Phenotype.

Prevailing dogma holds that ribosomes are uniform in composition and function. Here, we show that nutrient limitation-induced stress in E. coli changes the relative expression of rDNA operons to alter the rRNA composition within the actively translating ribosome pool. The most upregulated operon encodes the unique 16S rRNA, rrsH, distinguished by conserved sequence variation within the small ribosomal subunit. rrsH-bearing ribosomes affect the expression of functionally coherent gene sets and alter the levels of the RpoS sigma factor, the master regulator of the general stress response. These impacts are associated with phenotypic changes in antibiotic sensitivity, biofilm formation, and cell motility and are regulated by stress response proteins, RelA and RelE, as well as the metabolic enzyme and virulence-associated protein, AdhE. These findings establish that endogenously encoded, naturally occurring rRNA sequence variation can modulate ribosome function, central aspects of gene expression regulation, and cellular physiology.

Using controls to limit false discovery in the era of big data.

BACKGROUND: Procedures for controlling the false discovery rate (FDR) are widely applied as a solution to the multiple comparisons problem of high-dimensional statistics. Current FDR-controlling procedures require accurately calculated p-values and rely on extrapolation into the unknown and unobserved tails of the null distribution. Both of these intermediate steps are challenging and can compromise the reliability of the results. RESULTS: We present a general method for controlling the FDR that capitalizes on the large amount of control data often found in big data studies to avoid these frequently problematic intermediate steps. The method utilizes control data to empirically construct the distribution of the test statistic under the null hypothesis and directly compares this distribution to the empirical distribution of the test data. By not relying on p-values, our control data-based empirical FDR procedure more closely follows the foundational principles of the scientific method: that inference is drawn by comparing test data to control data. The method is demonstrated through application to a problem in structural genomics. CONCLUSIONS: The method described here provides a general statistical framework for controlling the FDR that is specifically tailored for the big data setting. By relying on empirically constructed distributions and control data, it forgoes potentially problematic modeling steps and extrapolation into the unknown tails of the null distribution. This procedure is broadly applicable insofar as controlled experiments or internal negative controls are available, as is increasingly common in the big data setting.

Variant ribosomal RNA alleles are conserved and exhibit tissue-specific expression.

The ribosome, the integration point for protein synthesis in the cell, is conventionally considered a homogeneous molecular assembly that only passively contributes to gene expression. Yet, epigenetic features of the ribosomal DNA (rDNA) operon and changes in the ribosome's molecular composition have been associated with disease phenotypes, suggesting that the ribosome itself may possess inherent regulatory capacity. Analyzing whole-genome sequencing data from the 1000 Genomes Project and the Mouse Genomes Project, we find that rDNA copy number varies widely across individuals, and we identify pervasive intra- and interindividual nucleotide variation in the 5S, 5.8S, 18S, and 28S ribosomal RNA (rRNA) genes of both human and mouse. Conserved rRNA sequence heterogeneities map to functional centers of the assembled ribosome, variant rRNA alleles exhibit tissue-specific expression, and ribosomes bearing variant rRNA alleles are present in the actively translating ribosome pool. These findings provide a critical framework for exploring the possibility that the expression of genomically encoded variant rRNA alleles gives rise to physically and functionally heterogeneous ribosomes that contribute to mammalian physiology and human disease.

An exact test for comparing a fixed quantitative property between gene sets.

Motivation: A significant difference in the distribution of a feature between two gene sets can provide insight into function or regulation. This statistical setting differs from much of hypothesis testing theory because the genome is often considered to be effectively fixed, finite and entirely known in commonly studied organisms, such as human. The Mann-Whitney U test is commonly employed in this scenario despite the assumptions of the test not being met, leading to unreliable and generally underpowered results. Permutation tests are also commonly employed for this purpose, but are computationally burdensome and are not tractable for obtaining small P values or for multiple comparisons. Results: We present an exact test for the null hypothesis that gene set membership is independent of the quantitative gene feature of interest. We derive an analytic expression for the randomization distribution of the median of the quantitative feature under the null hypothesis. Efficient implementation permits calculation of precise P values of arbitrary magnitude and makes thousands of simultaneous tests of transcriptome-sized gene sets computationally tractable. The flexibility of the hypothesis testing framework presented permits extension to a variety of related tests commonly found in genomics. The exact test is used to identify signatures of translation control and protein function in the human genome. Availability and implementation: The exact test presented here is implemented in R in the package kpmt available on CRAN. Contact: map2085@med.cornell.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Wnt5a Signals through DVL1 to Repress Ribosomal DNA Transcription by RNA Polymerase I.

Ribosome biogenesis is essential for cell growth and proliferation and is commonly elevated in cancer. Accordingly, numerous oncogene and tumor suppressor signaling pathways target rRNA synthesis. In breast cancer, non-canonical Wnt signaling by Wnt5a has been reported to antagonize tumor growth. Here, we show that Wnt5a rapidly represses rDNA gene transcription in breast cancer cells and generates a chromatin state with reduced transcription of rDNA by RNA polymerase I (Pol I). These effects were specifically dependent on Dishevelled1 (DVL1), which accumulates in nucleolar organizer regions (NORs) and binds to rDNA regions of the chromosome. Upon DVL1 binding, the Pol I transcription activator and deacetylase Sirtuin 7 (SIRT7) releases from rDNA loci, concomitant with disassembly of Pol I transcription machinery at the rDNA promoter. These findings reveal that Wnt5a signals through DVL1 to suppress rRNA transcription. This provides a novel mechanism for how Wnt5a exerts tumor suppressive effects and why disruption of Wnt5a signaling enhances mammary tumor growth in vivo.

Genome Sequence and Analysis of Escherichia coli MRE600, a Colicinogenic, Nonmotile Strain that Lacks RNase I and the Type I Methyltransferase, EcoKI.

Escherichia coli strain MRE600 was originally identified for its low RNase I activity and has therefore been widely adopted by the biomedical research community as a preferred source for the expression and purification of transfer RNAs and ribosomes. Despite its widespread use, surprisingly little information about its genome or genetic content exists. Here, we present the first de novo assembly and description of the MRE600 genome and epigenome. To provide context to these studies of MRE600, we include comparative analyses with E. coli K-12 MG1655 (K12). Pacific Biosciences Single Molecule, Real-Time sequencing reads were assembled into one large chromosome (4.83 Mb) and three smaller plasmids (89.1, 56.9, and 7.1 kb). Interestingly, the 7.1-kb plasmid possesses genes encoding a colicin E1 protein and its associated immunity protein. The MRE600 genome has a G + C content of 50.8% and contains a total of 5,181 genes, including 4,913 protein-encoding genes and 268 RNA genes. We identified 41,469 modified DNA bases (0.83% of total) and found that MRE600 lacks the gene for type I methyltransferase, EcoKI. Phylogenetic, taxonomic, and genetic analyses demonstrate that MRE600 is a divergent E. coli strain that displays features of the closely related genus, Shigella. Nevertheless, comparative analyses between MRE600 and E. coli K12 show that these two strains exhibit nearly identical ribosomal proteins, ribosomal RNAs, and highly homologous tRNA species. Substantiating prior suggestions that MRE600 lacks RNase I activity, the RNase I-encoding gene, rna, contains a single premature stop codon early in its open-reading frame.

Functional Dynamics within the Human Ribosome Regulate the Rate of Active Protein Synthesis.

The regulation of protein synthesis contributes to gene expression in both normal physiology and disease, yet kinetic investigations of the human translation mechanism are currently lacking. Using single-molecule fluorescence imaging methods, we have quantified the nature and timing of structural processes in human ribosomes during single-turnover and processive translation reactions. These measurements reveal that functional complexes exhibit dynamic behaviors and thermodynamic stabilities distinct from those observed for bacterial systems. Structurally defined sub-states of pre- and post-translocation complexes were sensitive to specific inhibitors of the eukaryotic ribosome, demonstrating the utility of this platform to probe drug mechanism. The application of three-color single-molecule fluorescence resonance energy transfer (smFRET) methods further revealed a long-distance allosteric coupling between distal tRNA binding sites within ribosomes bearing three tRNAs, which contributed to the rate of processive translation.

Detecting non-allelic homologous recombination from high-throughput sequencing data.

Non-allelic homologous recombination (NAHR) is a common mechanism for generating genome rearrangements and is implicated in numerous genetic disorders, but its detection in high-throughput sequencing data poses a serious challenge. We present a probabilistic model of NAHR and demonstrate its ability to find NAHR in low-coverage sequencing data from 44 individuals. We identify NAHR-mediated deletions or duplications in 109 of 324 potential NAHR loci in at least one of the individuals. These calls segregate by ancestry, are more common in closely spaced repeats, often result in duplicated genes or pseudogenes, and affect highly studied genes such as GBA and CYP2E1.