Ordered processing of the polygenic transcripts from a mitochondrial tRNA gene cluster in K. lactis

The transcriptome of Candida albicans mitochondria and the evolution of organellar transcription units in yeasts

Background

Yeasts show remarkable variation in the organization of their mitochondrial genomes, yet there is little experimental data on organellar gene expression outside few model species. Candida albicans is interesting as a human pathogen, and as a representative of a clade that is distant from the model yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe. Unlike them, it encodes seven Complex I subunits in its mtDNA. No experimental data regarding organellar expression were available prior to this study.

Methods

We used high-throughput RNA sequencing and traditional RNA biology techniques to study the mitochondrial transcriptome of C. albicans strains BWP17 and SN148.

Results

The 14 protein-coding genes, two ribosomal RNA genes, and 24 tRNA genes are expressed as eight primary polycistronic transcription units. We also found transcriptional activity in the noncoding regions, and antisense transcripts that could be a part of a regulatory mechanism. The promoter sequence is a variant of the nonanucleotide identified in other yeast mtDNAs, but some of the active promoters show significant departures from the consensus. The primary transcripts are processed by a tRNA punctuation mechanism into the monocistronic and bicistronic mature RNAs. The steady state levels of various mature transcripts exhibit large differences that are a result of posttranscriptional regulation. Transcriptome analysis allowed to precisely annotate the positions of introns in the RNL (2), COB (2) and COX1 (4) genes, as well as to refine the annotation of tRNAs and rRNAs. Comparative study of the mitochondrial genome organization in various Candida species indicates that they undergo shuffling in blocks usually containing 2–3 genes, and that their arrangement in primary transcripts is not conserved. tRNA genes with their associated promoters, as well as GC-rich sequence elements play an important role in these evolutionary events.

Conclusions

The main evolutionary force shaping the mitochondrial genomes of yeasts is the frequent recombination, constantly breaking apart and joining genes into novel primary transcription units. The mitochondrial transcription units are constantly rearranged in evolution shaping the features of gene expression, such as the presence of secondary promoter sites that are inactive, or act as “booster” promoters, simplified transcriptional regulation and reliance on posttranscriptional mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2078-z) contains supplementary material, which is available to authorized users.

Complete sequence and analysis of the mitochondrial genome of the methylotrophic yeast Hansenula polymorpha DL-1

We determined the complete nucleotide sequence of the 41 719 bp mitochondrial genome of the methylotrophic yeast Hansenula polymorpha strain DL-1. It contains genes for three subunits of cytochrome oxidase (cox1, cox2 and cox3), three subunits of ATP synthase (atp6, atp8 and atp9), seven subunits of NADH dehydrogenase (nad1-6 and nad4L), apocytochrome b (cob), four endonuclease/maturase homologs, a ribosomal protein (rps3), large and small rRNAs and a complete set of tRNAs. The structural genes are organized in two major transcriptional units. Phylogenetic, gene content and gene order analyses revealed the close phylogenetic relationship between H. polymorpha and Brettanomyces custersianus, and support the assignment of strain DL-1 to a separate genus rather than including it in the polyphyletic genus Pichia.

RNA maturation in mitochondria of S. cerevisiae and S. pombe

Although the gene content is rather conserved, the genomes in mitochondria of yeasts vary dramatically in size [Clark-Walker, G.D., Evans, R.J., Hoeben, P., McArthur, C.R., 1985. Basis of diversity in yeast mitochondrial DNAs. In: Quagliariello, E.C., Palmieri, F., Saccone, C., Kroon, A.M. (Eds.). Achievements and Perspectives of Mitochondrial Research 2. Science Publishers, Amsterdam, pp. 71-78] and in the number of transcription units. Since the fidelity and processivity of the mitochondrial single-subunit phage-like RNA polymerase present in yeast mitochondria are certainly limited, one might speculate that the density of transcription initiation sites on the mitochondrial genomes is one of the factors influencing the genome size. In an effort to find common features among the apparent idiosyncrasies of Saccharomyces cerevisiae (with its extremely large mtDNA) and Schizosaccharomyces pombe (with its extremely small mitochondrial genome), the aim of this review is to compare recent data about transcription and generation of 5' and 3' ends of mature RNA transcripts in S. cerevisiae and in S. pombe. Both organisms are two attractive model systems enabling investigation of various aspects of mitochondrial genetics.

Transcription and RNA-processing in fission yeast mitochondria

We systematically examined transcription and RNA-processing in mitochondria of the petite-negative fission yeast Schizosaccharomyces pombe. Two presumptive transcription initiation sites at opposite positions on the circular-mapping mtDNA were confirmed by in vitro capping of primary transcripts with guanylyl-transferase. The major promoter (Pma) is located adjacent to the 5'-end of the rnl gene, and a second, minor promoter (Pmi) upstream from cox3. The primary 5'-termini of the mature rnl and cox3 transcripts remain unmodified. A third predicted accessory transcription initiation site is within the group IIA1 intron of the cob gene (cobI1). The consensus promoter motif of S. pombe closely resembles the nonanucleotide promoter motifs of various yeast mtDNAs. We further characterized all mRNAs and the two ribosomal RNAs by Northern hybridization, and precisely mapped their 5'- and 3'-ends. The mRNAs have leader sequences with a length of 38 up to 220 nt and, in most instances, are created by removal of tRNAs from large precursor RNAs. Like cox2 and rnl, cox1 and cox3 are not separated by tRNA genes; instead, transcription initiation from the promoters upstream from rnl and cox3 compensates for the lack of tRNA-mediated 5'-processing. The 3'-termini of mRNAs and of SSU rRNA are processed at distinct, C-rich motifs that are located at a variable distance (1-15 nt) downstream from mRNA and SSU-rRNA coding regions. The accuracy of RNA-processing at these sites is sequence-dependent. Similar 3'-RNA-processing motifs are present in species of the genus Schizosaccharomyces, but not in budding yeasts that have functionally analogous A+T-rich dodecamer processing signals.

Highly diverged homologs of Saccharomyces cerevisiae mitochondrial mRNA-specific translational activators have orthologous functions in other budding yeasts

Translation of mitochondrially coded mRNAs in Saccharomyces cerevisiae depends on membrane-bound mRNA-specific activator proteins, whose targets lie in the mRNA 5'-untranslated leaders (5'-UTLs). In at least some cases, the activators function to localize translation of hydrophobic proteins on the inner membrane and are rate limiting for gene expression. We searched unsuccessfully in divergent budding yeasts for orthologs of the COX2- and COX3-specific translational activator genes, PET111, PET54, PET122, and PET494, by direct complementation. However, by screening for complementation of mutations in genes adjacent to the PET genes in S. cerevisiae, we obtained chromosomal segments containing highly diverged homologs of PET111 and PET122 from Saccharomyces kluyveri and of PET111 from Kluyveromyces lactis. All three of these genes failed to function in S. cerevisiae. We also found that the 5'-UTLs of the COX2 and COX3 mRNAs of S. kluyveri and K. lactis have little similarity to each other or to those of S. cerevisiae. To determine whether the PET111 and PET122 homologs carry out orthologous functions, we deleted them from the S. kluyveri genome and deleted PET111 from the K. lactis genome. The pet111 mutations in both species prevented COX2 translation, and the S. kluyveri pet122 mutation prevented COX3 translation. Thus, while the sequences of these translational activator proteins and their 5'-UTL targets are highly diverged, their mRNA-specific functions are orthologous.

Transcript mapping and processing of mitochondrial RNA in the chlorophyte alga Prototheca wickerhamii

The detailed transcript map of the circular 55328 bp mitochondrial (mt) genome from the colourless chlorophycean alga Prototheca wickerhamii has been determined. On each half of this genome the genes are encoded only on one DNA strand, forming transcriptional units comprising variable numbers of genes. With the exception of four genes coding for ribosomal proteins, transcripts of the three rRNA genes and all protein-coding genes have been detected by both northern analysis and primer extension experiments. Polycistronic transcripts of protein coding and tRNA genes were verified by northern analyses, primer extension and RNAse mapping experiments. The 5' and 3' ends of different RNA species are often located in close proximity to putative stem-loop structures and some 5' termini of mRNAs coincide with the 3' end of tRNAs located immediately upstream. Transcript mapping in a putative promoter region revealed two different possible transcription initiation sites; no significant sequence homology to putative mt promoters from higher plants could be found. In addition, two out of three group I introns residing in the cox1 gene were found to be self-splicing in vitro under reaction conditions developed for related mt introns from a filamentous fungus. Mitochondrial gene expression of P. wickerhamii and of filamentous fungi has several features in common, such as intron splicing and the processing of longer polycistronic transcripts. The similarities in RNA maturation between higher-plant and P. wickerhamii mitochondria are less pronounced, since plants rarely use tRNAs as processing signals for their relatively short mitochondrial co-transcripts.

tRNA genes transcribed from the plastid-like DNA of Plasmodium falciparum

Besides their mitochondrial genome, malarial parasites contain a second organellar DNA. This 35 kb circular molecule has a number of features reminiscent of plastid DNAs. Sequence analysis shows that along with other genes the circle codes for 25 different tRNAs all of which are transcribed. Six of the tRNAs have some unusual features, and one has an intron, the only one found so far on the circle. Comparison of codon and anticodon usage indicates that the 25 tRNAs are sufficient to decode all the protein genes present on the circle. The maintenance of such a parsimonious but complete translation system is further evidence for the functionality of the circle.