Ribosome profiling reveals post-transcriptional buffering of divergent gene expression in yeast - PubMed (original) (raw)

Ribosome profiling reveals post-transcriptional buffering of divergent gene expression in yeast

C Joel McManus et al. Genome Res. 2014 Mar.

Abstract

Understanding the patterns and causes of phenotypic divergence is a central goal in evolutionary biology. Much work has shown that mRNA abundance is highly variable between closely related species. However, the extent and mechanisms of post-transcriptional gene regulatory evolution are largely unknown. Here we used ribosome profiling to compare transcript abundance and translation efficiency in two closely related yeast species (S. cerevisiae and S. paradoxus). By comparing translation regulatory divergence to interspecies differences in mRNA sequence features, we show that differences in transcript leaders and codon bias substantially contribute to divergent translation. Globally, we find that translation regulatory divergence often buffers species differences in mRNA abundance, such that ribosome occupancy is more conserved than transcript abundance. We used allele-specific ribosome profiling in interspecies hybrids to compare the relative contributions of cis- and trans-regulatory divergence to species differences in mRNA abundance and translation efficiency. The mode of gene regulatory divergence differs for these processes, as trans-regulatory changes play a greater role in divergent mRNA abundance than in divergent translation efficiency. Strikingly, most genes with aberrant transcript abundance in F1 hybrids (either over- or underexpressed compared to both parent species) did not exhibit aberrant ribosome occupancy. Our results show that interspecies differences in translation contribute substantially to the evolution of gene expression. Compensatory differences in transcript abundance and translation efficiency may increase the robustness of gene regulation.

PubMed Disclaimer

Figures

Figure 1.

Figure 1.

Overview of allele-specific ribosome profiling to measure divergence in ribosome occupancy, mRNA abundance, and translational efficiency. (A) Ribosome profiling was performed on log-phase cultures of S. cerevisiae, S. paradoxus, and their F1 hybrid. Ribosome protected fragments (RPFs) were purified and cloned into Illumina high-throughput sequencing libraries (left). Poly-adenylated mRNA sequencing libraries were prepared in parallel (right). (B) Sequence reads from each sample were aligned to both species genomes. Allele-specific reads were identified by comparing genomic alignments and mapped to corresponding regions of orthologous ORFs. (C) Comparison of separate and allele-specific coverage in ribosome profiling experiments. IGV browser tracks showing normalized coverage of RPF and mRNA sequence reads from S. cerevisiae (blue) and S. paradoxus (magenta) over the COX6 gene (YHR051W). Measurements of interspecies differences in read coverage are equivalent for separate species (upper) and allele-specific alignments (“mock hybrid,” lower) (see also Supplemental Fig. S1).

Figure 2.

Figure 2.

Comparing regulatory divergence of ribosome occupancy, mRNA abundance, and translation efficiency. (A–C) Scatter plots compare the normalized average number of sequence reads for S. cerevisiae (_x_-axis) and S. paradoxus (_y_-axis). Genes with statistically significant differences in read counts (FDR < 5%, minimum 1.5-fold difference) are plotted as open circles with black edges. Translation efficiency is defined here as the number of ribosome protected fragment reads (RPF) divided by the number of mRNA-seq reads covering an ORF.

Figure 3.

Figure 3.

Translation regulatory divergence buffers interspecies differences in mRNA abundance. (A) Cartoon depicting buffering (left) and amplification (right). mRNA are shown as blue lines, and ribosomes are shown as black circles. Buffered genes have divergent mRNA abundance, with less divergent ribosome occupancy such that protein production is more conserved. In contrast, amplified genes have divergent mRNA abundance and even more divergent ribosome occupancy. (B) Scatterplot comparing divergent translation efficiency (_y_-axis, RPF/mRNA) with divergent mRNA abundance (_x_-axis). Buffered and amplified genes are plotted in red and blue, respectively. The negative correlation between mRNA abundance and translation efficiency suggests a genome-wide trend toward buffering. (C) Example of RPF and RNA-seq coverage over FPK1, a buffered gene. IGV browser tracks showing normalized coverage for S. cerevisiae (blue) and S. paradoxus (red).

Figure 4.

Figure 4.

Contributions of _cis_- and _trans_-regulatory divergence in mRNA abundance, ribosome occupancy, and translation efficiency. (A) Bar plot shows the number of genes affected by significant regulatory divergence in _cis-_acting sequences (C) and _trans_-acting factors (T). _Trans_-acting factors contribute most to mRNA abundance. (B) Box plot showing the fraction of regulatory divergence attributable to differences in _cis_-regulatory elements (%cis). Ribosome occupancy and translation efficiency both have higher %cis than mRNA abundance. (C) Venn diagram showing the overlap of genes with _cis_-regulatory divergence in mRNA abundance, ribosome occupancy, and translation efficiency. (D) Box plot depicting %cis for genes with buffering (Buf) and amplifying (Amp) regulatory divergence of mRNA abundance (left) and translation efficiency (right). Asterisks indicate results of Wilcoxon rank-sum test comparisons: (*) P < 0.001; (**) P < 0.0005.

Figure 5.

Figure 5.

Inheritance of gene expression in F1 hybrids of S. cerevisiae and S. paradoxus. (A) Hypothetical patterns of gene expression in S. paradoxus (red), S. cerevisiae (blue), and F1 hybrid yeast (purple), depicting six modes of gene expression inheritance. (B) Scatterplots showing the difference between mRNA expression levels (left) in the F1 hybrid and S. cerevisiae (_x_-axis) and S. paradoxus (_y_-axis). The difference between ribosome occupancy levels (RPF) in the F1 hybrid and parental species is shown on the right. Bar plots show the number of genes in each inheritance category from ribosome occupancy data. Genes with overdominant (overexpressed) and underdominant (underexpressed) inheritance occur much less frequently when considering ribosome occupancy, as compared to misexpression in mRNA abundance.

Similar articles

Cited by

References

    1. Altman DG, Bland JM 2003. Interaction revisited: The difference between two estimates. BMJ 326: 219. - PMC - PubMed
    1. Angov E 2011. Codon usage: Nature's roadmap to expression and folding of proteins. Biotechnol J 6: 650–659 - PMC - PubMed
    1. Benjamini Y, Hochberg Y 1995. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol 57: 289–300
    1. Brar GA, Yassour M, Friedman N, Regev A, Ingolia NT, Weissman JS 2011. High-resolution view of the yeast meiotic program revealed by ribosome profiling. Science 335: 552–557 - PMC - PubMed
    1. Brar GA, Rouskin S, McGeachy AM, Ingolia NT, Weissman JS 2012. The ribosome profiling strategy for monitoring translation in vivo by deep sequencing of ribosome-protected mRNA fragments. Nat Protoc 7: 1534–1550 - PMC - PubMed

Publication types

MeSH terms

Substances

LinkOut - more resources