Meiosis-induced alterations in transcript architecture and noncoding RNA expression in S. cerevisiae

Meiosis in budding yeast

Meiosis is a specialized cell division program that is essential for sexual reproduction. The two meiotic divisions reduce chromosome number by half, typically generating haploid genomes that are packaged into gametes. To achieve this ploidy reduction, meiosis relies on highly unusual chromosomal processes including the pairing of homologous chromosomes, assembly of the synaptonemal complex, programmed formation of DNA breaks followed by their processing into crossovers, and the segregation of homologous chromosomes during the first meiotic division. These processes are embedded in a carefully orchestrated cell differentiation program with multiple interdependencies between DNA metabolism, chromosome morphogenesis, and waves of gene expression that together ensure the correct number of chromosomes is delivered to the next generation. Studies in the budding yeast Saccharomyces cerevisiae have established essentially all fundamental paradigms of meiosis-specific chromosome metabolism and have uncovered components and molecular mechanisms that underlie these conserved processes. Here, we provide an overview of all stages of meiosis in this key model system and highlight how basic mechanisms of genome stability, chromosome architecture, and cell cycle control have been adapted to achieve the unique outcome of meiosis.

Gametogenesis: Exploring an Endogenous Rejuvenation Program to Understand Cellular Aging and Quality Control

Gametogenesis is a conserved developmental program whereby a diploid progenitor cell differentiates into haploid gametes, the precursors for sexually reproducing organisms. In addition to ploidy reduction and extensive organelle remodeling, gametogenesis naturally rejuvenates the ensuing gametes, leading to resetting of life span. Excitingly, ectopic expression of the gametogenesis-specific transcription factor Ndt80 is sufficient to extend life span in mitotically dividing budding yeast, suggesting that meiotic rejuvenation pathways can be repurposed outside of their natural context. In this review, we highlight recent studies of gametogenesis that provide emerging insight into natural quality control, organelle remodeling, and rejuvenation strategies that exist within a cell. These include selective inheritance, programmed degradation, and de novo synthesis, all of which are governed by the meiotic gene expression program entailing many forms of noncanonical gene regulation. Finally, we highlight critical questions that remain in the field and provide perspective on the implications of gametogenesis research on human health span.

RNA-Mediated Regulation of Meiosis in Budding Yeast

Cells change their physiological state in response to environmental cues. In the absence of nutrients, unicellular fungi such as budding yeast exit mitotic proliferation and enter the meiotic cycle, leading to the production of haploid cells that are encased within spore walls. These cell state transitions are orchestrated in a developmentally coordinated manner. Execution of the meiotic cell cycle program in budding yeast, Saccharomyces cerevisiae, is regulated by the key transcription factor, Ime1. Recent developments have uncovered the role of non-coding RNA in the regulation of Ime1 and meiosis. In this review, we summarize the role of ncRNA-mediated and RNA homeostasis-based processes in the regulation of meiosis in Saccharomyces cerevisiae.

Meiotic cDNA libraries reveal gene truncations and mitochondrial proteins important for competitive fitness in Saccharomyces cerevisiae

Gametogenesis is an evolutionarily conserved developmental program whereby a diploid progenitor cell undergoes meiosis and cellular remodeling to differentiate into haploid gametes, the precursors for sexual reproduction. Even in the simple eukaryotic organism Saccharomyces cerevisiae, the meiotic transcriptome is very rich and complex, thereby necessitating new tools for functional studies. Here, we report the construction of 5 stage-specific, inducible complementary DNA libraries from meiotic cells that represent over 84% of the genes found in the budding yeast genome. We employed computational strategies to detect endogenous meiotic transcript isoforms as well as library-specific gene truncations. Furthermore, we developed a robust screening pipeline to test the effect of each complementary DNA on competitive fitness. Our multiday proof-of-principle time course revealed 877 complementary DNAs that were detrimental for competitive fitness when overexpressed. The list included mitochondrial proteins that cause dose-dependent disruption of cellular respiration as well as library-specific gene truncations that expose a dominant negative effect on competitive growth. Together, these high-quality complementary DNA libraries provide an important tool for systematically identifying meiotic genes, transcript isoforms, and protein domains that are important for a specific biological function.

Long undecoded transcript isoform (LUTI) detection in meiotic budding yeast by direct RNA and transcript leader sequencing

LUTIs (Long Undecoded Transcript Isoforms) are 5'-extended and poorly translated mRNAs that can downregulate transcription from promoters more proximal to a gene's coding sequence (CDS). In this protocol, polyA RNA is extracted from budding yeast cells undergoing highly synchronized meiosis. Using a combination of long-read direct RNA sequencing and transcript leader sequencing (TL-seq), meiosis-specific LUTIs are systematically identified. Following identification, TL-seq is used to quantify the abundance of both LUTI and the more canonical gene-proximal (PROX) transcripts. For complete details on the use and execution of this protocol, please refer to Tresenrider et al. (2021).

Translation Initiation Site Profiling Reveals Widespread Synthesis of Non-AUG-Initiated Protein Isoforms in Yeast

Genomic analyses in budding yeast have helped define the foundational principles of eukaryotic gene expression. However, in the absence of empirical methods for defining coding regions, these analyses have historically excluded specific classes of possible coding regions, such as those initiating at non-AUG start codons. Here, we applied an experimental approach to globally annotate translation initiation sites in yeast and identified 149 genes with alternative N-terminally extended protein isoforms initiating from near-cognate codons upstream of annotated AUG start codons. These isoforms are produced in concert with canonical isoforms and translated with high specificity, resulting from initiation at only a small subset of possible start codons. The non-AUG initiation driving their production is enriched during meiosis and induced by low eIF5A, which is seen in this context. These findings reveal widespread production of non-canonical protein isoforms and unexpected complexity to the rules by which even a simple eukaryotic genome is decoded.

Developmentally regulated internal transcription initiation during meiosis in budding yeast

Sporulation of budding yeast is a developmental process in which cells undergo meiosis to generate stress-resistant progeny. The dynamic nature of the budding yeast meiotic transcriptome has been well established by a number of genome-wide studies. Here we develop an analysis pipeline to systematically identify novel transcription start sites that reside internal to a gene. Application of this pipeline to data from a synchronized meiotic time course reveals over 40 genes that display specific internal initiations in mid-sporulation. Consistent with the time of induction, motif analysis on upstream sequences of these internal transcription start sites reveals a significant enrichment for the binding site of Ndt80, the transcriptional activator of middle sporulation genes. Further examination of one gene, MRK1, demonstrates the Ndt80 binding site is necessary for internal initiation and results in the expression of an N-terminally truncated protein isoform. When the MRK1 paralog RIM11 is downregulated, the MRK1 internal transcript promotes efficient sporulation, indicating functional significance of the internal initiation. Our findings suggest internal transcriptional initiation to be a dynamic, regulated process with potential functional impacts on development.

Transcription of a 5' extended mRNA isoform directs dynamic chromatin changes and interference of a downstream promoter

Cell differentiation programs require dynamic regulation of gene expression. During meiotic prophase in Saccharomyces cerevisiae, expression of the kinetochore complex subunit Ndc80 is downregulated by a 5' extended long undecoded NDC80 transcript isoform. Here we demonstrate a transcriptional interference mechanism that is responsible for inhibiting expression of the coding NDC80 mRNA isoform. Transcription from a distal NDC80 promoter directs Set1-dependent histone H3K4 dimethylation and Set2-dependent H3K36 trimethylation to establish a repressive chromatin state in the downstream canonical NDC80 promoter. As a consequence, NDC80 expression is repressed during meiotic prophase. The transcriptional mechanism described here is rapidly reversible, adaptable to fine-tune gene expression, and relies on Set2 and the Set3 histone deacetylase complex. Thus, expression of a 5' extended mRNA isoform causes transcriptional interference at the downstream promoter. We demonstrate that this is an effective mechanism to promote dynamic changes in gene expression during cell differentiation.

Ndt80 activates the meiotic ORC1 transcript isoform and SMA2 via a bi-directional middle sporulation element in Saccharomyces cerevisiae

The origin of replication complex subunit ORC1 is important for DNA replication. The gene is known to encode a meiotic transcript isoform (mORC1) with an extended 5'-untranslated region (5'-UTR), which was predicted to inhibit protein translation. However, the regulatory mechanism that controls the mORC1 transcript isoform is unknown and no molecular biological evidence for a role of mORC1 in negatively regulating Orc1 protein during gametogenesis is available. By interpreting RNA profiling data obtained with growing and sporulating diploid cells, mitotic haploid cells, and a starving diploid control strain, we determined that mORC1 is a middle meiotic transcript isoform. Regulatory motif predictions and genetic experiments reveal that the activator Ndt80 and its middle sporulation element (MSE) target motif are required for the full induction of mORC1 and the divergently transcribed meiotic SMA2 locus. Furthermore, we find that the MSE-binding negative regulator Sum1 represses both mORC1 and SMA2 during mitotic growth. Finally, we demonstrate that an MSE deletion strain, which cannot induce mORC1, contains abnormally high Orc1 levels during post-meiotic stages of gametogenesis. Our results reveal the regulatory mechanism that controls mORC1, highlighting a novel developmental stage-specific role for the MSE element in bi-directional mORC1/SMA2 gene activation, and correlating mORC1 induction with declining Orc1 protein levels. Because eukaryotic genes frequently encode multiple transcripts possessing 5'-UTRs of variable length, our results are likely relevant for gene expression during development and disease in higher eukaryotes.

SUB1 Plays a Negative Role during Starvation Induced Sporulation Program in Saccharomyces cerevisiae

Saccharomyces cerevisiae Sub1 is involved in several cellular processes such as, transcription initiation, elongation, mRNA processing and DNA repair. It has also been reported to provide cellular resistance during conditions of oxidative DNA damage and osmotic stress. Here, we report a novel role of SUB1 during starvation stress-induced sporulation, which leads to meiosis and spore formation in diploid yeast cells. Deletion of SUB1 gene significantly increased sporulation efficiency as compared to the wild-type cells in S288c genetic background. Whereas, the sporulation functions of the sub1(Y66A) missense mutant were similar to Sub1. SUB1 transcript and protein levels are downregulated during sporulation, in highly synchronized and sporulation proficient wild-type SK1 cells. The changes in Sub1 levels during sporulation cascade correlate with the induction of middle sporulation gene expression. Deletion of SUB1 increased middle sporulation gene transcript levels with no effect on their induction kinetics. In wild-type cells, Sub1 associates with chromatin at these loci in a temporal pattern that correlates with their enhanced gene expression seen in sub1Δ cells. We show that SUB1 genetically interacts with HOS2, which led us to speculate that Sub1 might function with Set3 repressor complex during sporulation. Positive Cofactor 4, human homolog of Sub1, complemented the sub1Δ sporulation phenotype, suggesting conservation of function. Taken together, our results suggest that SUB1 acts as a negative regulator of sporulation.

Global alterations of the transcriptional landscape during yeast growth and development in the absence of Ume6-dependent chromatin modification

Chromatin modification enzymes are important regulators of gene expression and some are evolutionarily conserved from yeast to human. Saccharomyces cerevisiae is a major model organism for genome-wide studies that aim at the identification of target genes under the control of conserved epigenetic regulators. Ume6 interacts with the upstream repressor site 1 (URS1) and represses transcription by recruiting both the conserved histone deacetylase Rpd3 (through the co-repressor Sin3) and the chromatin-remodeling factor Isw2. Cells lacking Ume6 are defective in growth, stress response, and meiotic development. RNA profiling studies and in vivo protein-DNA binding assays identified mRNAs or transcript isoforms that are directly repressed by Ume6 in mitosis. However, a comprehensive understanding of the transcriptional alterations, which underlie the complex ume6Δ mutant phenotype during fermentation, respiration, or sporulation, is lacking. We report the protein-coding transcriptome of a diploid MAT a/α wild-type and ume6/ume6 mutant strains cultured in rich media with glucose or acetate as a carbon source, or sporulation-inducing medium. We distinguished direct from indirect effects on mRNA levels by combining GeneChip data with URS1 motif predictions and published high-throughput in vivo Ume6-DNA binding data. To gain insight into the molecular interactions between successive waves of Ume6-dependent meiotic genes, we integrated expression data with information on protein networks. Our work identifies novel Ume6 repressed genes during growth and development and reveals a strong effect of the carbon source on the derepression pattern of transcripts in growing and developmentally arrested ume6/ume6 mutant cells. Since yeast is a useful model organism for chromatin-mediated effects on gene expression, our results provide a rich source for further genetic and molecular biological work on the regulation of cell growth and cell differentiation in eukaryotes.

The histone deacetylase Rpd3/Sin3/Ume6 complex represses an acetate-inducible isoform of VTH2 in fermenting budding yeast cells

The tripartite Rpd3/Sin3/Ume6 complex represses meiotic isoforms during mitosis. We asked if it also controls starvation-induced isoforms. We report that VTH1/VTH2 encode acetate-inducible isoforms with extended 5'-regions overlapping antisense long non-coding RNAs. Rpd3 and Ume6 repress the long isoform of VTH2 during fermentation. Cells metabolising glucose contain Vth2, while the protein is undetectable in acetate and during sporulation. VTH2 is a useful model locus to study mechanisms implicating promoter directionality, lncRNA transcription and post-transcriptional control of gene expression via 5'-UTRs. Since mammalian genes encode transcript isoforms and Rpd3 is conserved, our findings are relevant for gene expression in higher eukaryotes.

The conserved histone deacetylase Rpd3 and its DNA binding subunit Ume6 control dynamic transcript architecture during mitotic growth and meiotic development

It was recently reported that the sizes of many mRNAs change when budding yeast cells exit mitosis and enter the meiotic differentiation pathway. These differences were attributed to length variations of their untranslated regions. The function of UTRs in protein translation is well established. However, the mechanism controlling the expression of distinct transcript isoforms during mitotic growth and meiotic development is unknown. In this study, we order developmentally regulated transcript isoforms according to their expression at specific stages during meiosis and gametogenesis, as compared to vegetative growth and starvation. We employ regulatory motif prediction, in vivo protein-DNA binding assays, genetic analyses and monitoring of epigenetic amino acid modification patterns to identify a novel role for Rpd3 and Ume6, two components of a histone deacetylase complex already known to repress early meiosis-specific genes in dividing cells, in mitotic repression of meiosis-specific transcript isoforms. Our findings classify developmental stage-specific early, middle and late meiotic transcript isoforms, and they point to a novel HDAC-dependent control mechanism for flexible transcript architecture during cell growth and differentiation. Since Rpd3 is highly conserved and ubiquitously expressed in many tissues, our results are likely relevant for development and disease in higher eukaryotes.

Developmental stage dependent metabolic regulation during meiotic differentiation in budding yeast

Background

The meiotic developmental pathway in yeast enables both differentiation of vegetative cells into haploid spores that ensure long-term survival, and recombination of the parental DNA to create genetic diversity. Despite the importance of proper metabolic regulation for the supply of building blocks and energy, little is known about the reprogramming of central metabolic pathways in meiotically differentiating cells during passage through successive developmental stages.

Results

Metabolic regulation during meiotic differentiation in budding yeast was analyzed by integrating information on genome-wide transcriptional activity, 26 enzymatic activities in the central metabolism, the dynamics of 67 metabolites, and a metabolic flux analysis at mid-stage meiosis. Analyses of mutants arresting sporulation at defined stages demonstrated that metabolic reprogramming is tightly controlled by the progression through the developmental pathway. The correlation between transcript levels and enzymatic activities in the central metabolism varies significantly in a developmental stage-dependent manner. The complete loss of phosphofructokinase activity at mid-stage meiosis enables a unique setup of the glycolytic pathway which facilitates carbon flux repartitioning into synthesis of spore wall precursors during the co-assimilation of glycogen and acetate. The need for correct homeostasis of purine nucleotides during the meiotic differentiation was demonstrated by the sporulation defect of the AMP deaminase mutant amd1, which exhibited hyper-accumulation of ATP accompanied by depletion of guanosine nucleotides.

Conclusions

Our systems-level analysis shows that reprogramming of the central metabolism during the meiotic differentiation is controlled at different hierarchical levels to meet the metabolic and energetic needs at successive developmental stages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-014-0060-x) contains supplementary material, which is available to authorized users.

SPO24 is a transcriptionally dynamic, small ORF-encoding locus required for efficient sporulation in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, meiosis and sporulation are highly regulated responses that are driven in part by changes in RNA expression. Alternative mRNA forms with extended 5′ UTRs are atypical in S. cerevisiae, and 5′ extensions with upstream open reading frames (uORFs) are even more unusual. Here we characterize the gene YPR036W-A, now renamed SPO24, which encodes a very small (67-amino-acid) protein. This gene gives rise to two mRNA forms: a shorter form throughout meiosis and a longer, 5′-extended form in mid-late meiosis. The latter form includes a uORF for a 14-amino-acid peptide (Spo24^u14). Deletion of the downstream ORF (dORF) leads to sporulation defects and the appearance of pseudohyphae-like projections. Experiments with luciferase reporters indicate that the uORF does not downregulate dORF translation. The protein encoded by the dORF (Spo24^d67) localizes to the prospore membrane and is differentially phosphorylated during meiosis. Transcription of the 5′-extended mRNA in mid-meiosis depends upon the presence of two middle sporulation elements (MSEs). Removal of the MSEs severely inhibits the mid-meiotic appearance of the 5′-extended mRNA and limits the ability of plasmid-borne SPO24 to rescue the sporulation defect of a spo24Δ mutant, suggesting that the 5′-extended mRNA is functionally important. These results reveal Spo24^d67 as a sporulation-related factor that is encoded by a transcriptionally dynamic, uORF-containing locus.

The splicing fate of plant SPO11 genes

Toward the global understanding of plant meiosis, it seems to be essential to decipher why all as yet sequenced plants need or at least encode for two different meiotic SPO11 genes. This is in contrast to mammals and fungi, where only one SPO11 is present. Both SPO11 in Arabidopsis thaliana are essential for the initiation of double strand breaks (DSBs) during the meiotic prophase. In nearly all eukaryotic organisms DSB induction during prophase I by SPO11 leads to meiotic DSB repair, thereby ensuring the formation of a necessary number of crossovers (CO) as physical connections between the homologous chromosomes. We aim to investigate the specific functions and evolution of both SPO11 genes in land plants. Therefore, we identified and cloned the respective orthologous genes from Brassica rapa, Carica papaya, Oryza sativa, and Physcomitrella patens. In parallel we determined the full length cDNA sequences of SPO11-1 and -2 from all of these plants by RT-PCR. During these experiments we observed that the analyzed plants exhibit a pattern of alternative splicing products of both SPO11 mRNAs. Such an aberrant splicing has previously been described for Arabidopsis and therefore seems to be conserved throughout evolution. Most of the splicing forms of SPO11-1 and -2 seem to be non-functional as they either showed intron retention (IR) or shortened exons. However, the positional distribution and number of alternative splicing events vary strongly between the different plants. The cDNAs showed in most cases premature termination codons (PTCs) due to frameshift. Nevertheless, in some cases we found alternatively spliced but functional cDNAs. These findings let us suggest that alternative splicing of SPO11 depends on the respective gene sequence and on the plant species. Therefore, this conserved mechanism might play a role concerning regulation of SPO11.

Extensive transcript diversity and novel upstream open reading frame regulation in yeast

To understand the diversity of transcripts in yeast (Saccharomyces cerevisiae) we analyzed the transcriptional landscapes for cells grown under 18 different environmental conditions. Each sample was analyzed using RNA-sequencing, and a total of 670,446,084 uniquely mapped reads and 377,263 poly-adenylated end tags were produced. Consistent with previous studies, we find that the majority of yeast genes are expressed under one or more different conditions. By directly comparing the 5′ and 3′ ends of the transcribed regions, we find extensive differences in transcript ends across many conditions, especially those of stationary phase, growth in grape juice, and salt stimulation, suggesting differential choice of transcription start and stop sites is pervasive in yeast. Relative to the exponential growth condition (i.e., YPAD), transcripts differing at the 5′ ends and 3′ ends are predicted to differ in their annotated start codon in 21 genes and their annotated stop codon in 63 genes. Many (431) upstream open reading frames (uORFs) are found in alternate 5′ ends and are significantly enriched in transcripts produced during the salt response. Mutational analysis of five genes with uORFs revealed that two sets of uORFs increase the expression of a reporter construct, indicating a role in activation which had not been reported previously, whereas two other uORFs decreased expression. In addition, RNA binding protein motifs are statistically enriched for alternate ends under many conditions. Overall, these results demonstrate enormous diversity of transcript ends, and that this heterogeneity is regulated under different environmental conditions. Moreover, transcript end diversity has important biological implications for the regulation of gene expression. In addition, our data also serve as a valuable resource for the scientific community.