Population structure of mitochondrial genomes in Saccharomyces cerevisiae

Background

Rigorous study of mitochondrial functions and cell biology in the budding yeast, Saccharomyces cerevisiae has advanced our understanding of mitochondrial genetics. This yeast is now a powerful model for population genetics, owing to large genetic diversity and highly structured populations among wild isolates. Comparative mitochondrial genomic analyses between yeast species have revealed broad evolutionary changes in genome organization and architecture. A fine-scale view of recent evolutionary changes within S. cerevisiae has not been possible due to low numbers of complete mitochondrial sequences.

Results

To address challenges of sequencing AT-rich and repetitive mitochondrial DNAs (mtDNAs), we sequenced two divergent S. cerevisiae mtDNAs using a single-molecule sequencing platform (PacBio RS). Using de novo assemblies, we generated highly accurate complete mtDNA sequences. These mtDNA sequences were compared with 98 additional mtDNA sequences gathered from various published collections. Phylogenies based on mitochondrial coding sequences and intron profiles revealed that intraspecific diversity in mitochondrial genomes generally recapitulated the population structure of nuclear genomes. Analysis of intergenic sequence indicated a recent expansion of mobile elements in certain populations. Additionally, our analyses revealed that certain populations lacked introns previously believed conserved throughout the species, as well as the presence of introns never before reported in S. cerevisiae.

Conclusions

Our results revealed that the extensive variation in S. cerevisiae mtDNAs is often population specific, thus offering a window into the recent evolutionary processes shaping these genomes. In addition, we offer an effective strategy for sequencing these challenging AT-rich mitochondrial genomes for small scale projects.

Electronic supplementary material

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

Regulatory interactions between mitochondrial genes: Exon and intron phenotypes observed in vivo can be expressed in vitro

In Saccharomyces cerevisiae the mitochondrial gene responsible for the specification of apocytochrome b (cob-box) is believed to consist of both coding and intervening sequences. Mutations in the latter give rise to pleiotropic phenotypes in vivo, lacking not only cytochrome b but also subunit I of cytochrome oxidase, and producing sets of novel polypeptides. The experiments described here have examined 15 different mit^? mutants in this region and demonstrate that these results are faithfully reproduced by isolated mitochondria in vitro. This inference also applies to other types of mutational lesions in coding segments of the cob-box gene and of the gene oxi3, responsible for the specification of subunit I.     

Genetic and physical analysis of the mitochondrial gene for subunit II of yeast cytochrome c oxidase.

The mitochondrial genetic locus oxi 1 contains the structural gene for subunit II of Cytochrome c oxidase. In this study, the oxi 1 locus, or at least a major portion of it, has been localized to a 2·4 kb2 HpaII fragment of mitochondrial DNA, by examining the mtDNA of oxi 1 mutants, and rho^? yeast strains that selectively retained in amplified form, this region of the mitochondria) genome. The 2·4 kb fragment is missing from the mtDNA of an oxi 1 locus deletion mutant, but is present in the mtDNAs retained by two rho^? strains that genetically recombine with all 16 oxi 1 mutants tested, to produce respiring progeny. Two other rho^? strains, that retained different but overlapping portions of the oxi 1 locus as determined genetically, contained mtDNAs consisting of “cloned” segments derived from within the 2·4 kb fragment: these rho^? mtDNAs hybridized only to the 2·4 kb HpaII fragment of wild-type mtDNA and could not be cleaved with HpaII. Furthermore, these two rho^? mtDNAs were found to correspond to sequences from opposite sides of the 2·4 kb fragment that overlap for 100 to 300 base-pairs near the middle of the fragment. Thus, five oxi 1 mutations that recombine with both of these rho^? strains could be further localized to this relatively short region of overlap. One such mutation, of particular interest because it produces an altered form of subunit II, was shown to lie on a 75-base-pair fragment that maps in this region of the overlap. The 75-base-pair fragment from the mutant migrates slightly faster during electrophoresis than the corresponding wild-type fragment. In contrast, the mobility of the fragment from a spontaneous revertant was indistinguishable from wild type.     

Evidence for homologous recombination between repeated sequences containing 18S and 5S ribosomal RNA genes in wheat mitochondrial DNA

Closely linked genes for 18S and 5S rRNAs have been located on four different cloned SalI restriction fragments of wheat (Triticum aestivum L.) mitochondrial DNA. Restriction analysis has revealed that in each of the cloned fragments, the 18S and 5S rRNA genes are contained within the same basic structural unit, R, which is at least 4 kbp long. This unit is flanked by sequences designated u (0.8 kbp), v (13.7 kbp), w (0.7 kbp), and y (1.4 kbp), in the orientations v-R-w, v-R-y, u-R-w, and u-R-y in the four different SalI fragments. We conclude that 18S + 5S rRNA genes are located at several distinct sites in the wheat mitochondrial genome, and suggest that reciprocal intra- and/or intermolecular recombination between such repeated sequences could promote extensive genomic rearrangement and thus contribute to the physical heterogeneity that is a hallmark of most plant mitochondrial DNAs.     

Larger rearranged mitochondrial genomes in Dekkera/Brettanomyces yeasts are more closely related than smaller genomes with a conserved gene order

Summary Mitochondrial genomes from yeasts in the Dekkera/Brettanomyces/Eeniella group vary in size from 28 to 101 kb. Mapping of genes has shown that the three smallest genomes, of 28–42 kb, have the same gene order, whereas the three larger mitochondrial DNAs of 57–101 kb are rearranged relative to the smaller molecules and between themselves. To examine the relationships between these genomes, a phylogenetic tree has been constructed by sequence comparison of the mitochondrialencoded cytochrome oxidase subunit gene (COX2) from the six species. Contrary to expectation, the tree shows that the larger rearranged genomes are more closely related than the smaller mtDNAs. This result indicates that the gene order of the smaller mtDNAs (28–42 kb) is ancestral and that larger mtDNA molecules (57–101 kb) are more prone to rearrangement than smaller forms.Offprint requests to: G.D. Clark-Walker     

The complete mitochondrial genomes of two brown algae (Laminariales, Phaeophyceae) and phylogenetic analysis within Laminaria

Two complete mitochondrial genome sequences for Laminaria longissima (=Saccharina longissima) and Laminaria hyperborea are reported in this study. They had circular mapping organization with slight difference in size (37,628 and 37,976 bp, respectively) and contained almost the same set of mitochondrial genes, including the genes for three rRNAs (23S, 16S, and 5S), 25 tRNAs, 35 known mitochondrial proteins, and three to four Open Reading Frame genes (ORFs). Both mitochondrial genomes exhibited typical gene content and organization of Laminaria mtDNAs except for the existence of ORF157 genes being located between rRNA large subunit gene 5 (rpl5) and ORF129-139 in L. hyperborea as found in that of Laminaria digitata. The phylogenetic analysis based on mitochondrial genomes supported the hypothesis of the split of the genus Laminaria, and the result of this study provided important information on the molecular evolution of Laminaria.     

Nuclear genotype affects mitochondrial genome organization of CMS-S maize

Summary A WF9 strain of maize with the RD subtype of the S male-sterile cytoplasm (CMS-S) was converted to the inbred M825 nuclear background by recurrent backcrossing. The organization of the mitochondrial genomes of the F₁ and succeeding backcross progenies was analyzed and compared with the progenitor RD-WF9 using probes derived from the S1 and S2 mitochondrial episomes, and probes containing the genes for cytochrome c oxidase subunit I (coxI), cytochrome c oxidase subunit II (coxII) and apocytochrome b (cob). Changes in mitochondrial DNA (mtDNA) organization were observed for S1-, S2-, and coxI-homologous sequences that involve loss of homologous restriction enzyme fragments present in the RD-WF9 progenitor. With the coxI probe, the loss of certain fragments was accompanied by the appearance of a fragment not detectable in the progenitor. The changes observed indicate the effect of the nuclear genome on the differential replication of specific mitochondrial subgenomic entities.     

A comparison of three fission yeast mitochondrial genomes

The fission yeasts are members of the fungal order Schizosaccharomycetales, a candidate deep-diverging group within Ascomycota. Although a great deal of molecular information is available from Schizosaccharomyces pombe, a model eukaryote, very little is available from other members of this group. In order to better characterize mitochondrial genome evolution in this fungal lineage, the mitochondrial DNA (mtDNA) of two additional fission yeasts, Schizosaccharomyces octosporus and Schizosaccharomyces japonicus var. japonicus, was sequenced. Whereas the mtDNA of S.pombe is only 19 431 bp, the mtDNA of S.octosporus is 44 227 bp, and that of S.japonicus var. japonicus is over 80 kb. The size variation of these mtDNAs is due largely to non-coding regions. The gene content in the latter two mtDNAs is almost identical to that of the completely sequenced S.pombe mtDNA, which encodes 25 tRNA species, the large and small mitochondrial ribosomal RNAs (rnl and rns), the RNA component of mitochondrial RNaseP (rnpB), mitochondrial small subunit ribosomal protein 3 (rps3), cytochrome oxidase subunits 1, 2 and 3 (cox1, cox2 and cox3) and ATP-synthase subunits 6, 8 and 9 (atp6, atp8 and atp9). However, trnI2(cau) (C modified to lysidine) is absent in the S.octosporus mtDNA, as are corresponding ATA codons in its protein-coding genes, and rps3 and rnpB are not found in the mtDNA of S.japonicus var. japonicus. The mtDNA of S.octosporus contains five double hairpin elements, the first report of these elements in an ascomycete. This study provides further evidence in favor of the mobility of these elements, and supports their role in mitochondrial genome rearrangement. The results of our phylogenetic analysis support the monophyly of the Schizosaccharomycetales, but question their grouping within the Archiascomycota.     

Mitochondrial DNA rearrangements in Pennisetum associated with reversion from cytoplasmic male sterility to fertility

Endonuclease restriction fragment patterns of Pennisetum americanum L. mitochondrial DNAs (mtDNAs) from a cytoplasmic male-sterile (CMS-A1), fertile revertants and a normal fertile cytoplasm were variable, while chloroplast DNA from those lines lacked variation. Comparisons between mtDNAs of CMS-A1 (parental) and fertile revertant lines revealed the presence of a unique 4.7 kbp PstI fragment in the sterile line that was not detected in any of the revertant lines. A 9.7 kbp PstI fragment was found in all of the revertants, but not in the CMS-A1. Neither of those fragments was found in the normal cytoplasm mtDNA. Hybridization studies revealed two sets of multiple homologies: 1) the 4.7 kbp fragment had homology with a 10.9 kbp and a 13.6 kbp fragment; and 2) the 9.7 kbp fragment was homologous with the 13.6 kbp fragment. The presence of those two repeated mitochondrial sequences on the altered fragments suggests that they may be involved in the recombinational associated events with reversion from CMS to fertility in P. americanum.Florida Agricultural Experiment Station Journal Series No.7797.     

Rearrangements in sugar beet mitochondrial DNA induced by cell suspension,callus cultures and regeneration

Structural alterations in mitochondrial DNAs (mtDNAs) from a plant of a sterile sugar beet line, callus derived from it, suspension-cultured cells and plants regenerated from the callus were studied. BamHI restriction analysis revealed that structural alterations between the mtDNAs of the callus and the control plant had occurred. Multiple rearrangements were also demonstrated in the mtDNA from the suspension culture, of which some were similar to those appearing in the callus, and others had arisen de novo. Rearrangements were also identified by means of blot hybridization of BamHI-digested mtDNA from suspension-cultured cells with the genes encoding subunit II of cytochrome oxidase (cox II) and subunit 1 of NADH-dehydrogenase (Nd1). No alterations were observed in the mitochondrial genome of the callus and regenerants. The location of the genes for the -subunit of F1-ATPase (atpA) and apocytochrome b (cob) in the mtDNA remained unchanged.Our salient finding was of a plant with an altered mitochondrial genome as judged by EcoRI and BamHI restriction analysis. This exceptional plant had retained the sterile phenotype like all of the other regenerants and the parent. The set of plasmid-like molecules of mtDNA remained the same as that in the control plant and in all of the regenerants, callus and suspension-cultured cells. The only type of plasmid-like molecule found in all of the DNAs was the 1.6-kbp minicircle, which is a feature of sterile cytoplasms. These structural changes in mtDNA were obviously a consequence of somaclonal variation during the in vitro cultivation of the sugar beet cells.     

The Largest Mitochondrial Translation Product Copurifying with the Mitochondrial Adenosine Triphosphatase of Saccharomyces cerevisiae Is Not a Subunit of the Enzyme Complex

Mitochondrial adenosine triphosphatase isolated from a double mutant of Saccharomyces cerevisiae lacking cytochrome b apoprotein and subunit II of cytochrome oxidase does not contain the mitochondrial translation product (approximate molecular weight, 32,000) previously suggested to be a subunit of the enzyme complex.     

Physical and genetic organization of Petite and Grande yeast mitochondrial DNA. III. High resolution melting and reassociation studies

Mitochondrial DNAs from 13 petite mutants have been analyzed by means of high-resolution melting and reassociation in solution, in an attempt to relate their physical chemical properties to the mitochondrial genotype, which displays various combinations of genetic marker deletions. The kinetic complexities of the petite mtDNAs were found to range from $\text{1}\text{3}$ down to $\text{1}\text{500}$ of that of the grande mtDNA; the loss in sequential complexity undergone by petite mtDNAs parallels the loss in mitochondrial genetic complexity. Melting profiles of petite mtDNAs can be resolved into well-defined peaks. Some of them are specific to the marker genes. The genotypic specificity increases as the sequential complexity of mtDNA decreases. The mtDNA region conferring resistance to erythromycin could in this way be shown to be characterized by two melting peaks, at 72 °C and 75 °C. The results are interpreted in terms of selective enrichment in gene-specific sequences.     

Nitrobacter winogradskyi cytochrome c oxidase genes are organized in a repeated gene cluster

Cytochrome c oxidase (EC 1.9.3.1) is one of the components of the electron transport chain by which Nitrobacter, a facultative lithoautotrophic bacterium, recovers energy from nitrite oxidation. The genes encoding the two catalytic core subunits of the enzyme were isolated from a Nitrobacter winogradskyi gene library. Sequencing of one of the 14 cloned DNA segments revealed that the subunit genes are side by side in an operon-like cluster. Remarkably the cluster appears to be present in at least two copies per genome. It extends over a 5–6 kb length including, besides the catalytic core subunit genes, other cytochrome oxidase related genes, especially a heme O synthase gene. Noteworthy is the new kind of gene order identified within the cluster. Deduced sequences for the cytochrome oxidase subunits and for the heme O synthase look closest to their counterparts in other α-subdivision Proteobacteria, particularly the Rhizobiaceae. This confirms the phylogenetic relationships established only upon 16S rRNA data. Furthermore, interesting similarities exist between N. winogradskyi and mitochondrial cytochrome oxidase subunits while the heme O synthase sequence gives some new insights about the other similar published α-subdivision proteobacterial sequences.     

Characterization of the complete mitochondrial genome of Bombyx mori strain H9 (Lepidoptera: Bombycidae)

The complete mitochondrial genome (mitogenome) of Bombyx mori strain H9 (Lepidoptera: Bombycidae) is 15,670 base pairs (bp) in length, encoding 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The nucleotide composition of the genome is highly A + T biased, accounting for 81.31%, with a slightly positive AT skewness (0.059). The arrangement of 13 PCGs is similar to that of other sequenced lepidopterans. All the PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is proposed by the TTAG sequence as observed in other lepidopterans. Unlike the other PCGs, the cox1 and cytochrome c oxidase subunit 2 (cox2) genes have incomplete stop codons consisting of just a T. All tRNAs have typical structures of insect mitochondrial tRNAs, which is different from other sequenced lepidopterans. The structure of A + T-rich region is similar to that of other sequenced lepidopterans, including non-repetitive sequences, the ATAGA binding domain, a 18 bp poly-T stretch and a poly-A element upstream of transfer RNA M (trnM) gene. Phylogenetic analysis shows that the domesticated silkmoth B. mori originated from the Chinese Bombyx mandarina.     

Cytochrome oxidase subunit II gene in mitochondria of Oenothera has no intron

The cytochrome oxidase subunit II gene has been localized in the mitochondrial genome of Oenothera berteriana and the nucleotide sequence has been determined. The coding sequence contains 777 bp and, unlike the corresponding gene in Zea mays, is not interrupted by an intron. No TGA codon is found within the open reading frame. The codon CGG, as in the maize gene, is used in place of tryptophan codons of corresponding genes in other organisms. At position 742 in the Oenothera sequence the TGG of maize is changed into a CGG codon, where Trp is conserved as the amino acid in other organisms. Homologous sequences occur more than once in the mitochondrial genome as several mitochondrial DNA species hybridize with DNA probes of the cytochrome oxidase subunit II gene.     

The structure of the small mitochondrial DNA of Kluyveromyces thermotolerans is likely to reflect the ancestral gene order in fungi

Mapping the 23-kb circular mitochondrial DNA from the yeast Kluyveromyces thermotolerans has shown that only one change occurs in the gene order in comparison to the 19-kb mtDNA of Candida (Torulopsis) glabrata. Sequence analysis of the mitochondrially encoded cytochrome oxidase subunit 2 gene reveals that despite their conserved gene order, the two small genomes are more distantly related than larger mtDNA molecules with multiple rearrangements. This result supports a previous observation that larger mitochondrial genomes are more prone to rearrange than smaller forms and suggests that the architecture of the two small molecules is likely to represent the structure of an ancestor.Correspondence to: G.D. Clark-Walker 0592