BSG Isoform 2 (ENST00000353555) Is a Prognostic Biomarker in Predicting Unfavorable Overall Survival of Patients With Liver Hepatocellular Carcinoma.

BACKGROUND: CD147, encoded by the BSGgene, has complex transcripts that encode proteins of different lengths. Total BSG transcription is a prognostic biomarker for patients with liver cancer. This study tried to analyze the expression profile and prognostic significance of BSG transcripts in liver cancer. MATERIALS AND METHODS: RNA sequencing data from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) project, survival data from TCGA, and protein expression data from the Human Protein Atlas were systematically analyzed. RESULTS: Among the four protein-coding transcripts of BSG, ENST00000353555 encoding basigin-2 is the dominant transcript isoform. It might be an independent prognostic biomarker for unfavorable overall survival in patients with liver cancer (HR: 1.404, 95% CI: 1.1224-1.754, p = 0.003). CONCLUSIONS: ENST00000353555 might be a prognostic biomarker linking unfavorable overall survival in liver cancer patients.

A Bayesian Change Point Model for Dynamic Alternative Transcription Start Site Usage During Cellular Differentiation.

ABSTRACT An alternative transcription start site (ATSS) is a major driving force for increasing the complexity of transcripts in human tissues. As a transcriptional regulatory mechanism, ATSS has biological significance. Many studies have confirmed that ATSS plays an important role in diseases and cell development and differentiation. However, exploration of its dynamic mechanisms remains insufficient. Identifying ATSS change points during cell differentiation is critical for elucidating potential dynamic mechanisms. For relative ATSS usage as percentage data, the existing methods lack sensitivity to detect the change point for ATSS longitudinal data. In addition, some methods have strict requirements for data distribution and cannot be applied to deal with this problem. In this study, the Bayesian change point detection model was first constructed using reparameterization techniques for two parameters of a beta distribution for the percentage data type, and the posterior distributions of parameters and change points were obtained using Markov Chain Monte Carlo (MCMC) sampling. With comprehensive simulation studies, the performance of the Bayesian change point detection model is found to be consistently powerful and robust across most scenarios with different sample sizes and beta distributions. Second, differential ATSS events in the real data, whose change points were identified using our method, were clustered according to their change points. Last, for each change point, pathway and transcription factor motif analyses were performed on its differential ATSS events. The results of our analyses demonstrated the effectiveness of the Bayesian change point detection model and provided biological insights into cell differentiation.

ASAPA: a bioinformatic pipeline based on Iso-Seq that identifies the links among alternative splicing, alternative transcription initiation and alternative polyadenylation.

BACKGROUND: Although the events associated with alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI) can be identified by many approaches based on isoform sequencing (Iso-Seq), these analyses are generally independent of each other and the links between these events are still rarely mentioned. However, an interdependency analysis can be achieved because the transcriptional start site, splice sites and polyA site could be simultaneously included in a long, full-length read from Iso-Seq. RESULTS: We create ASAPA pipeline that enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data. We tested this pipeline using Arabidopsis data and found some interesting results: some adjacent introns tend to be simultaneously spliced or retained; coupling between AS and ATI or APA is limited to the initial or terminal intron; and ATI and APA are potentially linked in some special cases. CONCLUSION: Our pipeline enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data, which is conducive to a better understanding of transcription landscape generation.

BORIS/CTCFL epigenetically reprograms clustered CTCF binding sites into alternative transcriptional start sites.

BACKGROUND: Pervasive usage of alternative promoters leads to the deregulation of gene expression in carcinogenesis and may drive the emergence of new genes in spermatogenesis. However, little is known regarding the mechanisms underpinning the activation of alternative promoters. RESULTS: Here we describe how alternative cancer-testis-specific transcription is activated. We show that intergenic and intronic CTCF binding sites, which are transcriptionally inert in normal somatic cells, could be epigenetically reprogrammed into active de novo promoters in germ and cancer cells. BORIS/CTCFL, the testis-specific paralog of the ubiquitously expressed CTCF, triggers the epigenetic reprogramming of CTCF sites into units of active transcription. BORIS binding initiates the recruitment of the chromatin remodeling factor, SRCAP, followed by the replacement of H2A histone with H2A.Z, resulting in a more relaxed chromatin state in the nucleosomes flanking the CTCF binding sites. The relaxation of chromatin around CTCF binding sites facilitates the recruitment of multiple additional transcription factors, thereby activating transcription from a given binding site. We demonstrate that the epigenetically reprogrammed CTCF binding sites can drive the expression of cancer-testis genes, long noncoding RNAs, retro-pseudogenes, and dormant transposable elements. CONCLUSIONS: Thus, BORIS functions as a transcription factor that epigenetically reprograms clustered CTCF binding sites into transcriptional start sites, promoting transcription from alternative promoters in both germ cells and cancer cells.

Identification of novel TMEM231 gene splice variants and pathological findings in a fetus with Meckel Syndrome.

Background: Meckel Syndrome (MKS, OMIM #249000) is a rare and fatal autosomal recessive ciliopathy with high clinical and genetic heterogeneity. MKS shows complex allelism with other related ciliopathies such as Joubert Syndrome (JBTS, OMIM #213300). In MKS, the formation and function of the primary cilium is defective, resulting in a multisystem disorder including occipital encephalocele, polycystic kidneys, postaxial polydactyly, liver fibrosis, central nervous system malformations and genital anomalies. This study aimed to analyze the genotype of MKS patients and investigate the correlation between genotype and phenotype. Methods: A nonconsanguineous couple who conceived four times with a fetus affected by multiorgan dysfunction and intrauterine fetal death was studied. Whole exome sequencing (WES) was performed in the proband to identify the potentially pathogenic variant. Sanger sequencing was performed in family members. In silico tools were used to analyse the pathogenicity of the identified variants. cDNA TA-cloning sequencing was performed to validate the effects of intronic variants on mRNA splicing. Quantitative real-time PCR was performed to investigate the effect of the variants on gene expression. Immunofluorescence was performed to observe pathological changes of the primary cilium in kidney tissue from the proband. Results: Two splice site variants of TMEM231 (NM_001077418.2, c.583-1G>C and c.583-2_588delinsTCCTCCC) were identified in the proband, and the two variants have not been previously reported. The parents were confirmed as carriers. The two variants were predicted to be pathogenic by in silico tools and were classified as pathogenic/likely pathogenic variants according to the American College of Medical Genetics and Genomics guideline. cDNA TA cloning analysis showed that both splice site variants caused a deletion of exon 5. RT-PCR revealed that the expression of TMEM231 was significantly decreased and immunofluorescence showed that the primary cilium was almost absent in the proband's kidney tissue. Conclusion: We reported the clinical, genetic, molecular and histochemical characterisation of a family affected by MKS. Our findings not only extended the mutation spectrum of the TMEM231 gene, but also revealed for the first time the pathological aetiology of primary cilia in humans and provide a basis for genetic counselling of the parents to their offspring.

The single-cell landscape of alternative transcription start sites of diabetic retina.

More than half of mammalian protein-coding genes have multiple transcription start sites. Alternative transcription start site (TSS) modulate mRNA stability, localization, and translation efficiency on post-transcription level, and even generate novel protein isoforms. However, differential TSS usage among cell types in healthy and diabetic retina remains poorly characterized. In this study, by using 5'-tag-based single-cell RNA sequencing, we identified cell type-specific alternative TSS events and key transcription factors for each of retinal cell types. We observed that lengthening of 5'- UTRs in retinal cell types are enriched for multiple RNA binding protein binding sites, including splicing regulators Rbfox1/2/3 and Nova1. Furthermore, by comparing TSS expression between healthy and diabetic retina, we identified elevated apoptosis signal in Müller glia and microglia, which can be served as a putative early sign of diabetic retinopathy. By measuring 5'UTR isoforms in retinal single-cell dataset, our work provides a comprehensive panorama of alternative TSS and its potential consequence related to post-transcriptional regulation. We anticipate our assay can not only provide insights into cellular heterogeneity driven by transcriptional initiation, but also open up the perspectives for identification of novel diagnostic indexes for diabetic retinopathy.

The alternative transcription and expression characterization of Dmc1 in autotetraploid Carassius auratus.

Established autotetraploids often have a highly stable meiosis with high fertility compared with neo-autotetraploids. The autotetraploid Carassius auratus (4n = 200, RRRR) (4nRR), which stemmed from whole-genome duplication of Carassius auratus red var. (2n = 100, RR) (RCC), produces diploid gametes with an adopted diploid-like chromosome pairing in meiosis and maintains the formation of autotetraploid lineages. In this study, we focused on Dmc1, a meiosis-specific recombinase during the prophase of meiosis I, and elaborated on the genetic variation, alternative transcription, expression characterization, and epigenetic modification of Dmc1 in RCC and 4nRR. Two original Dmc1 from RCC were identified in 4nRR, and two duplicated Dmc1 differences in genetic composition were observed in 4nRR. Furthermore, we only noticed that one original and one duplicated Dmc1 were expressed in RCC and 4nRR, respectively. However, both possessed identical gene expression profiles, differential expression of sexual dimorphism, and hypomethylation levels. These results indicated that the specific expression of duplicated Dmc1 may be involve in the progression of meiosis of the diploid-like chromosome pairing in autotetraploid Carassius auratus. Herein, the findings significantly increase knowledge of meiosis of autopolyploid fish and provide meaningful insights into genetic breeding in polyploidy fish.

Genome-Wide Mapping of 5' Isoforms with 5'-Seq.

The transcriptome is far more complex than previously assumed. Transcripts from the same gene can differ in terms of transcription start site, transcription end site, or pattern of splicing, and growing evidence supports the functional importance of these distinct transcript isoforms. Easily identifying these isoforms experimentally via library construction and high-throughput sequencing is crucial. Current library construction methods for identifying transcription start sites (5' transcript isoforms) involve large number of steps and (expensive) reagents, utilization of cDNA intermediates for adapter ligation, and are less suitable for studying low-abundance isoforms. Here, I describe a quick protocol for the generation of sequencing libraries to define capped 5' isoforms (5'-Seq) of various abundances in yeast and suggest a 5' isoform data analysis pipeline. The protocol relies on the utilization of a dephosphorylation-decapping method (oligo-capping) to generate a sequencing library from mRNA fragments and is a simplification of previously published 5' isoform protocols in terms of the handling steps, time, and cost. This method is exemplified using Saccharomyces cerevisiae mRNA, but it can be applied to various cellular conditions to study the effects of 5' transcript isoforms on transcriptional and/or translational regulation. © 2023 Wiley Periodicals LLC. Basic Protocol: Construction of a DNA sequencing library from capped 5' isoforms Support Protocol: Sequencing data analysis.

Alternative Transcription Start Site Usage and Functional Implications in Pathogenic Fungi.

Pathogenic fungi require delicate gene regulation mechanisms to adapt to diverse living environments and escape host immune systems. Recent advances in sequencing technology have exposed the complexity of the fungal genome, thus allowing the gradual disentanglement of multiple layers of gene expression control. Alternative transcription start site (aTSS) usage, previously reported to be prominent in mammals and to play important roles in physiopathology, is also present in fungi to fine-tune gene expression. Depending on the alteration in their sequences, RNA isoforms arising from aTSSs acquire different characteristics that significantly alter their stability and translational capacity as well as the properties and biologic functions of the resulting proteins. Disrupted control of aTSS usage has been reported to severely impair growth, virulence, and the infectious capacity of pathogenic fungi. Here, we discuss principle concepts, mechanisms, and the functional implication of aTSS usage in fungi.

New start sites for plant immune gene transcription.

Appropriate spatiotemporal regulation of immune gene expression triggered by pathogen-associated molecular patterns (PAMPs) is critical for defense against pathogens. Thieffry et al. recently found that PAMP-triggered widespread alternative transcription initiation and an immediate transcription factor wave are elements of immune gene reprogramming.

Identifying genes with conserved splicing structure and orthologous isoforms in human, mouse and dog.

BACKGROUND: In eukaryote transcriptomes, a significant amount of transcript diversity comes from genes' capacity to generate different transcripts through alternative splicing. Identifying orthologous alternative transcripts across multiple species is of particular interest for genome annotators. However, there is no formal definition of transcript orthology based on the splicing structure conservation. Likewise there is no public dataset benchmark providing groups of orthologous transcripts sharing a conserved splicing structure. RESULTS: We introduced a formal definition of splicing structure orthology and we predicted transcript orthologs in human, mouse and dog. Applying a selective strategy, we analyzed 2,167 genes and their 18,109 known transcripts and identified a set of 253 gene orthologs that shared a conserved splicing structure in all three species. We predicted 6,861 transcript CDSs (coding sequence), mainly for dog, an emergent model species. Each predicted transcript was an ortholog of a known transcript: both share the same CDS splicing structure. Evidence for the existence of the predicted CDSs was found in external data. CONCLUSIONS: We generated a dataset of 253 gene triplets, structurally conserved and sharing all their CDSs in human, mouse and dog, which correspond to 879 triplets of spliced CDS orthologs. We have released the dataset both as an SQL database and as tabulated files. The data consists of the 879 CDS orthology groups with their detailed splicing structures, and the predicted CDSs, associated with their experimental evidence. The 6,861 predicted CDSs are provided in GTF files. Our data may contribute to compare highly conserved genes across three species, for comparative transcriptomics at the isoform level, or for benchmarking splice aligners and methods focusing on the identification of splicing orthologs. The data is available at https://data-access.cesgo.org/index.php/s/V97GXxOS66NqTkZ .

A novel computational framework for genome-scale alternative transcription units prediction.

Alternative transcription units (ATUs) are dynamically encoded under different conditions and display overlapping patterns (sharing one or more genes) under a specific condition in bacterial genomes. Genome-scale identification of ATUs is essential for studying the emergence of human diseases caused by bacterial organisms. However, it is unrealistic to identify all ATUs using experimental techniques because of the complexity and dynamic nature of ATUs. Here, we present the first-of-its-kind computational framework, named SeqATU, for genome-scale ATU prediction based on next-generation RNA-Seq data. The framework utilizes a convex quadratic programming model to seek an optimum expression combination of all of the to-be-identified ATUs. The predicted ATUs in Escherichia coli reached a precision of 0.77/0.74 and a recall of 0.75/0.76 in the two RNA-Sequencing datasets compared with the benchmarked ATUs from third-generation RNA-Seq data. In addition, the proportion of 5'- or 3'-end genes of the predicted ATUs, having documented transcription factor binding sites and transcription termination sites, was three times greater than that of no 5'- or 3'-end genes. We further evaluated the predicted ATUs by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes functional enrichment analyses. The results suggested that gene pairs frequently encoded in the same ATUs are more functionally related than those that can belong to two distinct ATUs. Overall, these results demonstrated the high reliability of predicted ATUs. We expect that the new insights derived by SeqATU will not only improve the understanding of the transcription mechanism of bacteria but also guide the reconstruction of a genome-scale transcriptional regulatory network.

Differential alternative RNA splicing and transcription events between tumors from African American and White patients in The Cancer Genome Atlas.

Individuals of African ancestry suffer disproportionally from higher incidence, aggressiveness, and mortality for particular cancers. This disparity likely results from an interplay among differences in multiple determinants of health, including differences in tumor biology. We used The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA aggregate expression datasets and identified differential alternative RNA splicing and transcription events (ARS/T) in cancers between self-identified African American (AA) and White (W) patients. We found that retained intron events were enriched among race-related ARS/T. In addition, on average, 12% of the most highly ranked race-related ARS/T overlapped between any two analyzed cancers. Moreover, the genes undergoing race-related ARS/T functioned in cancer-promoting pathways, and a number of race-related ARS/T were associated with patient survival. We built a web-application, CanSplice, to mine genomic datasets by self-identified race. The race-related targets have the potential to aid in the development of new biomarkers and therapeutics to mitigate cancer disparity.

Correlation of novel ALKATI with ALK immunohistochemistry and clinical outcomes in metastatic melanoma.

AIMS: Recently, a novel isoform of anaplastic lymphoma kinase, with alternative transcription initiation (ALKATI ), has been described in melanoma and is susceptible to targeted ALK-inhibitor therapy. Clinical outcomes of patients with ALKATI mutated melanoma as well as correlation with immunohistochemical (IHC) methods have not yet been described. METHODS AND RESULTS: Clinicopathological characteristics were abstracted for 324 patients with metastatic melanoma (MM). IHC, fluorescence in-situ hybridisation and RNA-based digital molecular analysis assays were performed on archival tissue from 173 stage III and 192 stage IV tumours. ALKATI was identified in 12.7 and 4.8% stage III and IV tumours, respectively. Discrete presentations of the ALKATI are seen: isolated ALKATI (n = 20) and mixed ALKATI (combined ALKATI and ALKWT ; n = 7). Isolated ALKWT expression (n = 4) was seen with no ALK fusions. Stage III patients showed improved survival with ALKATI expression compared to those with ALKWT or no expression [5-year survival 80, 95% confidence interval (CI) = 57-100% versus 43%, 95% CI = 34-55%, P = 0.013]. Clinicopathological characteristics were not statistically significant. Strong diffuse cytoplasmic staining of ALK IHC (n = 12) has a sensitivity of 52.2%, specificity 100%, PPV of 100% and NPV of 92.5% of detecting isolated ALKATI . CONCLUSION: Presence of ALKATI is a good prognostic indicator in MM. ALK IHC and digital molecular analysis can be incorporated into MM evaluation to identify patients with ALKATI for targeted therapy.

Alternative transcription start sites of the enolase-encoding gene enoA are stringently used in glycolytic/gluconeogenic conditions in Aspergillus oryzae.

Gene expression using alternative transcription start sites (TSSs) is an important transcriptional regulatory mechanism for environmental responses in eukaryotes. Here, we identify two alternative TSSs in the enolase-encoding gene (enoA) in Aspergillus oryzae, an industrially important filamentous fungus. TSS use in enoA is strictly dependent on the difference in glycolytic and gluconeogenic carbon sources. Transcription from the upstream TSS (uTSS) or downstream TSS (dTSS) predominantly occurs under gluconeogenic or glycolytic conditions, respectively. In addition to enoA, most glycolytic genes involved in reversible reactions possess alternative TSSs. The fbaA gene, which encodes fructose-bisphosphate aldolase, also shows stringent alternative TSS selection, similar to enoA. Alignment of promoter sequences of enolase-encoding genes in Aspergillus predicted two conserved regions that contain a putative cis-element required for enoA transcription from each TSS. However, uTSS-mediated transcription of the acuN gene, an enoA ortholog in Aspergillus nidulans, is not strictly dependent on carbon source, unlike enoA. Furthermore, enoA transcript levels in glycolytic conditions are higher than in gluconeogenic conditions. Conversely, acuN is more highly transcribed in gluconeogenic conditions. This suggests that the stringent usage of alternative TSSs and higher transcription in glycolytic conditions in enoA may reflect that the A. oryzae evolutionary genetic background was domesticated by exclusive growth in starch-rich environments. These findings provide novel insights into the complexity and diversity of transcriptional regulation of glycolytic/gluconeogenic genes among Aspergilli.

Tissue and cancer-specific expression of DIEXF is epigenetically mediated by an Alu repeat.

Alu repeats constitute a major fraction of human genome and for a small subset of them a role in gene regulation has been described. The number of studies focused on the functional characterization of particular Alu elements is very limited. Most Alu elements are DNA methylated and then assumed to lie in repressed chromatin domains. We hypothesize that Alu elements with low or variable DNA methylation are candidates for a functional role. In a genome-wide study in normal and cancer tissues, we pinpointed an Alu repeat (AluSq2) with differential methylation located upstream of the promoter region of the DIEXF gene. DIEXF encodes a highly conserved factor essential for the development of zebrafish digestive tract. To characterize the contribution of the Alu element to the regulation of DIEXF we analysed the epigenetic landscapes of the gene promoter and flanking regions in different cell types and cancers. Alternate epigenetic profiles (DNA methylation and histone modifications) of the AluSq2 element were associated with DIEXF transcript diversity as well as protein levels, while the epigenetic profile of the CpG island associated with the DIEXF promoter remained unchanged. These results suggest that AluSq2 might directly contribute to the regulation of DIEXF transcription and protein expression. Moreover, AluSq2 was DNA hypomethylated in different cancer types, pointing out its putative contribution to DIEXF alteration in cancer and its potential as tumoural biomarker.

N-Terminal Truncated Myb with New Transcriptional Activity Produced Through Use of an Alternative MYB Promoter in Salivary Gland Adenoid Cystic Carcinoma.

Adenoid cystic carcinoma (ACC) is an aggressive salivary gland tumor that frequently displays perineural invasion and is often associated with translocations or overexpression of the MYB oncogene. Detailed analyses of MYB transcripts from ACC patient samples revealed that ACC tumors utilize an alternative MYB promoter, which is rarely used in normal cells or other tumor types. The alternative promoter transcripts produce N-terminally truncated Myb proteins lacking a highly conserved and phosphorylated domain, which includes the pS11 epitope that is frequently used to detect Myb proteins. In RNA-seq assays, Myb isoforms lacking the N-terminal domain displayed unique transcriptional activities, regulating many genes differently than full-length Myb. Thus, a regulatory pathway unique to ACC activates the alternative MYB promoter, leading to the production of a truncated Myb protein with altered transcriptional activities. This could provide new therapeutic opportunities for ACC patients.

Fusion Transcripts of Adjacent Genes: New Insights into the World of Human Complex Transcripts in Cancer.

The awareness of genome complexity brought a radical approach to the study of transcriptome, opening eyes to single RNAs generated from two or more adjacent genes according to the present consensus. This kind of transcript was thought to originate only from chromosomal rearrangements, but the discovery of readthrough transcription opens the doors to a new world of fusion RNAs. In the last years many possible intergenic cis-splicing mechanisms have been proposed, unveiling the origins of transcripts that contain some exons of both the upstream and downstream genes. In some cases, alternative mechanisms, such as trans-splicing and transcriptional slippage, have been proposed. Five databases, containing validated and predicted Fusion Transcripts of Adjacent Genes (FuTAGs), are available for the scientific community. A comparative analysis revealed that two of them contain the majority of the results. A complete analysis of the more widely characterized FuTAGs is provided in this review, including their expression pattern in normal tissues and in cancer. Gene structure, intergenic splicing patterns and exon junction sequences have been determined and here reported for well-characterized FuTAGs. The available functional data and the possible roles in cancer progression are discussed.

mountainClimber Identifies Alternative Transcription Start and Polyadenylation Sites in RNA-Seq.

Alternative transcription start (ATS) and alternative polyadenylation (APA) create alternative RNA isoforms and modulate many aspects of RNA expression and protein production. However, ATS and APA remain difficult to detect in RNA sequencing (RNA-seq). Here, we developed mountainClimber, a de novo cumulative-sum-based approach to identify ATS and APA as change points. Unlike many existing methods, mountainClimber runs on a single sample and identifies multiple ATS or APA sites anywhere in the transcript. We analyzed 2,342 GTEx samples (36 tissues, 215 individuals) and found that tissue type is the predominant driver of transcript end variations. 75% and 65% of genes exhibited differential APA and ATS across tissues, respectively. In particular, testis displayed longer 5' untranslated regions (UTRs) and shorter 3' UTRs, often in genes related to testis-specific biology. Overall, we report the largest study of transcript ends across human tissues to our knowledge. mountainClimber is available at github.com/gxiaolab/mountainClimber.

Mitochondrial cysteinyl-tRNA synthetase is expressed via alternative transcriptional initiation regulated by energy metabolism in yeast cells.

Eukaryotes typically utilize two distinct aminoacyl-tRNA synthetase isoforms, one for cytosolic and one for mitochondrial protein synthesis. However, the genome of budding yeast (Saccharomyces cerevisiae) contains only one cysteinyl-tRNA synthetase gene (YNL247W, also known as CRS1). In this study, we report that CRS1 encodes both cytosolic and mitochondrial isoforms. The 5' complementary DNA end method and GFP reporter gene analyses indicated that yeast CRS1 expression yields two classes of mRNAs through alternative transcription starts: a long mRNA containing a mitochondrial targeting sequence and a short mRNA lacking this targeting sequence. We found that the mitochondrial Crs1 is the product of translation from the first initiation AUG codon on the long mRNA, whereas the cytosolic Crs1 is produced from the second in-frame AUG codon on the short mRNA. Genetic analysis and a ChIP assay revealed that the transcription factor heme activator protein (Hap) complex, which is involved in mitochondrial biogenesis, determines the transcription start sites of the CRS1 gene. We also noted that Hap complex-dependent initiation is regulated according to the needs of mitochondrial energy production. The results of our study indicate energy-dependent initiation of alternative transcription of CRS1 that results in production of two Crs1 isoforms, a finding that suggests Crs1's potential involvement in mitochondrial energy metabolism in yeast.