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.

Most organisms encode multiple, distinct copies of rRNA genes, rendering the composition of the ribosome pool intrinsically heterogeneous. Here, Kurylo et al. show that nutrient limitation in E. coli upregulates the expression of ribosomes bearing conserved sequence variation in 16S rRNA that can regulate gene expression and phenotype.

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.

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.

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.

Induction of cytochrome P450 family 1 mRNAs and activities in a cell line from the frog Xenopus laevis

Cytochrome P450 family 1 (CYP1) includes four subfamilies of enzymes: CYP1A, CYP1B, CYP1C, and CYP1D. In many vertebrates, CYP1A, 1B, and 1C expression is induced by agonists of the aryl hydrocarbon receptor, including toxic contaminants such as chlorinated dioxins, coplanar chlorinated biphenyls, and polynuclear aromatic hydrocarbons. Assessed at the level of mRNA, protein, or enzyme activity, CYP1s (especially CYP1As) represent potent and popular biomarkers of contaminant exposure in aquatic vertebrates. Alkylated resorufins are synthetic substrates used to detect, quantify, and describe catalytic activities of cytochrome P450s. The ability to oxidize specific resorufin-based substrates can distinguish the catalytic activities of individual CYP1s. Xenopus laevis, the African clawed frog, is the most widely employed amphibian model in aquatic toxicology, yet the number, inducibility, and activities of CYP1s have not been systematically characterized in this species. Here we report the cloning of cDNAs encoding two new CYP1 family members, X. laevis CYP1B and CYP1C, along with an integrated assessment of the induction of alkyloxyuresorufin-O-dealkylase (AROD) activities and mRNA expression of four known X. laevis CYP1s: CYP1A6, CYP1A7, CYP1B, and CYP1C. Using XLK-WG, an X. laevis kidney epithelial cell line, we determined that EROD (ethoxyresorufin substrate) and MROD (methoxyresorufin) were both induced 3000- to 5000-fold following 2,3,7,8 tetrachlorodibenzo-p-dioxin (TCDD) exposure up to 250 nM, while BROD (benzyloxyresorufin) and PROD (pentyloxyresorufin) activity was not detectable regardless of TCDD treatment. TCDD induced CYP1A6 and CYP1A7 mRNAs by 2-3 orders of magnitude, while CYP1B and CYP1C were unchanged. The more potent AHR agonist, FICZ (6-formylindolo[3,2-b]carbazole), induced CYP1B up to 10-fold at concentrations between 0.1 and 250 nM, while CYP1C induction was less than 3-fold. CYP1B mRNA showed the highest constitutive mRNA expression, 5- to 75-fold greater than the other CYP1 transcripts. Taken together, these results suggest that CYP1A6 and CYP1A7 perform the bulk of EROD and MROD activities we observed in these cells. The ability of each X. laevis CYP1 to catalyze oxidation of individual resorufin substrates remains to be determined. Correlating CYP1 mRNA and induced AROD activity is a significant step toward clarifying the biochemical meaning of these biomarkers and the roles of CYP1 enzymes in X. laevis. The cell culture approach represents an important complement to the long standing use of frog embryos and tadpoles in toxicological studies, providing a well suited model system for determining the molecular mechanisms underlying the regulation of these important biomarkers of contaminant exposure.