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.     

REVIEW—CONTRASTING MITOCHONDRIAL GENOME ORGANIZATIONS AND SEQUENCE AFFILIATIONS AMONG GREEN ALGAE: POTENTIAL FACTORS, MECHANISMS, AND EVOLUTIONARY SCENARIOS

The three green algal mitochondrial genomes completely sequenced to date — those of Chlamydomonas reinhardtii Dangeard, Chlamydomonas eugametos Gerloff, and Prototheca wickerhamii Soneda & Tubaki — revealed very different mitochondrial genome organizations and sequence affiliations. The Chlamydomonas genomes resemble the ciliate / fungal / animal counterparts, and the Prototheca genome resembles land plant homologues. This review points out that all the green algal mitochondrial genomes examined to date resemble either the Chlamydomonas or the Prototheca mitochondrial genome; the Chlamydomonas- like mitochondrial genomes are small and have a reduced gene content (no ribosomal protein or 5S rRNA genes and only a few protein-coding and tRNA genes) and fragmented and scrambled rRNA coding regions, whereas the Prototheca- like mitochondrial genomes are larger and have a larger set of protein-coding genes (including ribosomal protein genes), more tRNA genes, and 5S rRNA and conventional continuous small-subunit (SSU) and large-subunit (LSU) rRNA coding regions. It appears, therefore, that the differences previously observed between the mitochondrial genomes of C. reinhardtii and P. wickerhamii extend to the two green algal mitochondrial lineages to which they belong and are significant enough to raise questions about the causes and mechanisms responsible for such contrasting evolutionary strategies among green algae. This review suggests an integrative approach in explaining the occurrence of distinct evolutionary strategies and apparent phylogenetic affiliations among the known green algal mitochondrial lineages. The observed differences could be the result of distinct genetic potentials differentiated during the previous evolutionary history of the flagellate ancestors and / or of subsequent changes in habitat and life history of the more advanced green algal lineages.     

The mitochondrial genome of the legume Vigna radiata and the analysis of recombination across short mitochondrial repeats

The mitochondrial genomes of seed plants are exceptionally fluid in size, structure, and sequence content, with the accumulation and activity of repetitive sequences underlying much of this variation. We report the first fully sequenced mitochondrial genome of a legume, Vigna radiata (mung bean), and show that despite its unexceptional size (401,262 nt), the genome is unusually depauperate in repetitive DNA and "promiscuous" sequences from the chloroplast and nuclear genomes. Although Vigna lacks the large, recombinationally active repeats typical of most other seed plants, a PCR survey of its modest repertoire of short (38-297 nt) repeats nevertheless revealed evidence for recombination across all of them. A set of novel control assays showed, however, that these results could instead reflect, in part or entirely, artifacts of PCR-mediated recombination. Consequently, we recommend that other methods, especially high-depth genome sequencing, be used instead of PCR to infer patterns of plant mitochondrial recombination. The average-sized but repeat- and feature-poor mitochondrial genome of Vigna makes it ever more difficult to generalize about the factors shaping the size and sequence content of plant mitochondrial genomes.     

Seventeen new complete mtDNA sequences reveal extensive mitochondrial genome evolution within the Demospongiae

Two major transitions in animal evolution--the origins of multicellularity and bilaterality--correlate with major changes in mitochondrial DNA (mtDNA) organization. Demosponges, the largest class in the phylum Porifera, underwent only the first of these transitions and their mitochondrial genomes display a peculiar combination of ancestral and animal-specific features. To get an insight into the evolution of mitochondrial genomes within the Demospongiae, we determined 17 new mtDNA sequences from this group and analyzing them with five previously published sequences. Our analysis revealed that all demosponge mtDNAs are 16- to 25-kbp circular molecules, containing 13-15 protein genes, 2 rRNA genes, and 2-27 tRNA genes. All but four pairs of sampled genomes had unique gene orders, with the number of shared gene boundaries ranging from 1 to 41. Although most demosponge species displayed low rates of mitochondrial sequence evolution, a significant acceleration in evolutionary rates occurred in the G1 group (orders Dendroceratida, Dictyoceratida, and Verticillitida). Large variation in mtDNA organization was also observed within the G0 group (order Homosclerophorida) including gene rearrangements, loss of tRNA genes, and the presence of two introns in Plakortis angulospiculatus. While introns are rare in modern-day demosponge mtDNA, we inferred that at least one intron was present in cox1 of the common ancestor of all demosponges. Our study uncovered an extensive mitochondrial genomic diversity within the Demospongiae. Although all sampled mitochondrial genomes retained some ancestral features, including a minimally modified genetic code, conserved structures of tRNA genes, and presence of multiple non-coding regions, they vary considerably in their size, gene content, gene order, and the rates of sequence evolution. Some of the changes in demosponge mtDNA, such as the loss of tRNA genes and the appearance of hairpin-containing repetitive elements, occurred in parallel in several lineages and suggest general trends in demosponge mtDNA evolution.     

Characterization of long and short repetitive sequences in the sea urchin genome

Long and short repetitive sequences were purified from the DNA of Paracentrotus lividus under conditions designed to optimize the yield of complete, end to end sequences. Double-stranded long repeat DNA prepared in this manner ranged in length from approximately 3000 to 15 000 nucleotide pairs with average sizes of approximately 6000 base pairs. In the electron microscope, long repeat DNA was observed to possess continuous sequences that often appeared to be terminated by one or more loops and/or fold backs. Long repeat DNA sequences, resheared to 300 base pairs, were found to have an average melting point identical to that for sheared native DNA. Thus, the reassociated duplexes of long repetitive DNA seem to possess very few mismatched base pairs. Reassociation kinetic analyses indicate that the majority of the long repeat sequences are reiterated only 4--7 times per haploid amount of DNA. Melt-reassociation analyses of short repetitive DNA, at several criteria, support the previously held concept that these sequences belong the sets or families of sequences which are inexact copies of one another. Our studies also support hypotheses suggesting that short repetitive sequences belong to families which may have arisen via distinct salttatory events. The relationships between long and short repetitive DNA sequences are considered with respect to widely held concepts of their sequence organization, evolution, and possible functions within eucaryotic genomes. A model for the possible organization of short repeats within long repetitive DNA sequences is also presented.     

RNA-based phylogenetic methods: application to mammalian mitochondrial RNA sequences

The PHASE software package allows phylogenetic tree construction with a number of evolutionary models designed specifically for use with RNA sequences that have conserved secondary structure. Evolution in the paired regions of RNAs occurs via compensatory substitutions, hence changes on either side of a pair are correlated. Accounting for this correlation is important for phylogenetic inference because it affects the likelihood calculation. In the present study we use the complete set of tRNA and rRNA sequences from 69 complete mammalian mitochondrial genomes. The likelihood calculation uses two evolutionary models simultaneously for different parts of the sequence: a paired-site model for the paired sites and a single-site model for the unpaired sites. We use Bayesian phylogenetic methods and a Markov chain Monte Carlo algorithm is used to obtain the most probable trees and posterior probabilities of clades. The results are well resolved for almost all the important branches on the mammalian tree. They support the arrangement of mammalian orders within the four supra-ordinal clades that have been identified by studies of much larger data sets mainly comprising nuclear genes. Groups such as the hedgehogs and the murid rodents, which have been problematic in previous studies with mitochondrial proteins, appear in their expected position with the other members of their order. Our choice of genes and evolutionary model appears to be more reliable and less subject to biases caused by variation in base composition than previous studies with mitochondrial genomes.     

Microsatellite analysis in organelle genomes of Chlorophyta

Simple Sequence Repeats (SSRs) or microsatellites constitute a significant portion of genomes however; their significance in organellar genomes has not been completely understood. The availability of organelle genome sequences allows us to understand the organization of SSRs in their genic and intergenic regions. In the present work, SSRs were identified and categorized in 14 mitochondrial and 22 chloroplast genomes of algal species belonging to Chlorophyta. Based on the study, it was observed that number of SSRs in non-coding region were more as compared to coding region and frequency of mononucleotides repeats were highest followed by dinucleotides in both mitochondrial and chloroplast genomes. It was also observed that maximum number of SSRs was found in genes encoding for beta subunit of RNA polymerase in chloroplast genomes and NADH dehydrogenase in mitochondrial genomes. This is the first and original report on whole genomes sequence analysis of organellar genomes of green algae.     

Small, repetitive DNAs contribute significantly to the expanded mitochondrial genome of cucumber

Closely related cucurbit species possess eightfold differences in the sizes of their mitochondrial genomes. We cloned mitochondrial DNA (mtDNA) fragments showing strong hybridization signals to cucumber mtDNA and little or no signal to watermelon mtDNA. The cucumber mtDNA clones carried short (30-53 bp), repetitive DNA motifs that were often degenerate, overlapping, and showed no homology to any sequences currently in the databases. On the basis of dot-blot hybridizations, seven repetitive DNA motifs accounted for >13% (194 kb) of the cucumber mitochondrial genome, equaling >50% of the size of the Arabidopsis mitochondrial genome. Sequence analysis of 136 kb of cucumber mtDNA revealed only 11.2% with significant homology to previously characterized mitochondrial sequences, 2.4% to chloroplast DNA, and 15% to the seven repetitive DNA motifs. The remaining 71.4% of the sequence was unique to the cucumber mitochondrial genome. There was <4% sequence colinearity surrounding the watermelon and cucumber atp9 coding regions, and the much smaller watermelon mitochondrial genome possessed no significant amounts of cucumber repetitive DNAs. Our results demonstrate that the expanded cucumber mitochondrial genome is in part due to extensive duplication of short repetitive sequences, possibly by recombination and/or replication slippage.     

In silico analysis of SSRs in mitochondrial genomes of plants

Simple sequence repeats (SSRs) or microsatellites constitute a countable portion of genomes. However, the significance of SSRs in organelle genomes has not been completely understood. The availability of organelle genome sequences allows us to understand the organization of SSRs in their genic and intergenic regions. In the current study we surveyed the patterns of SSRs in mitochondrial genomes of different taxa of plants. A total of 16 mitochondrial genomes, from algae to angiosperms, have been considered to analyze the pattern of simple sequence repeats present in them. Based on study, the mononucleotide repeats of A/T were found to be more prevalent in mitochondrial genomes over other repeat types. The dinucleotides repeats, TA/AT, were the second most numerous, whereas tri-, tetra-, and pentanucleotide repeats were in less number and present in intronic or intergenic portions only. Mononucleotide repeats prevailed in protein-coding exonic portions of all organisms. These results indicates that microsatellite pattern in mitochondrial genomes is different from nuclear genomes and also focuses on organization and diversity at SSR locuses in mitochondrial genomes. This is the novel report of microsatellite polymorphism in plant mitochondrion on whole genome level.     

Comparison of base composition and codon usage in insect mitochondrial genomes

Insects, the most biodiverse taxonomic group, have high AT content in their mitochondrial genomes. Although codon usage tends to be AT-rich, base composition and codon usage of mitochondrial genomes may vary among taxa. Thus, we compare base composition and codon usage patterns of 49 insect mitochondrial genomes. For protein coding genes, AT content is as high as 80% in the Hymenoptera and Lepidoptera and as low as 72% in the Orthopotera. The AT content is high at positions 1 and 3, but A content is low at position 2. A close correlation occurs between codon usage and tRNA abundance in nuclear genomes. Optimal codons can pair well with the antr codons of the most abundant tRNAs. One tRNA gene translates a synonymous codon family in vertebrate mitochondrial genomes and these tRNA anticodons can pair with optimal codons. However, optimal codons cannot pair with anticodons in mtDNA ofCochiiomyia hominivorax (Dipteral: CaLliphoridae). Ten optimal codons cannot pair with tRNA anticodons in all 49 insect mitochondrial genomes; non-optimal codon-anticodon usage is common and codon usage is not influenced by tRNA abundance.     

Repetitive DNA sequence homologies and amplifications in South American cricetid rodents

The extent of conservation of repetitive DNA sequences in five species of South American cricetid rodents, belonging to three tribes and genera, exhibiting clear karyological differences has been investigated. Dot blot evaluation of interspecific hybridization percentages underscored sharp discrepancies in some reciprocal assays, suggesting sequence amplification events might have taken place since the evolutionary divergence of the species studied. Southern blot analysis, employing the most intense restriction bands of three species' DNAs as radiolabeled probes, allowed to confirm the sequence amplification events initially suspected. Although homologous repetitive sequences are conserved in the genomes of the spiecies investigated, differential organization of such sequences does exist, as judged by the species-specific Southern blot patterns.     

Complete mitochondrial genome sequences of three bats species and whole genome mitochondrial analyses reveal patterns of codon bias and lend support to a basal split in Chiroptera

Order Chiroptera is a unique group of mammals whose members have attained self-powered flight as their main mode of locomotion. Much speculation persists regarding bat evolution; however, lack of sufficient molecular data hampers evolutionary and conservation studies. Of ~ 1200 species, complete mitochondrial genome sequences are available for only eleven. Additional sequences should be generated if we are to resolve many questions concerning these fascinating mammals. Herein, we describe the complete mitochondrial genomes of three bats: Corynorhinus rafinesquii, Lasiurus borealis and Artibeus lituratus. We also compare the currently available mitochondrial genomes and analyze codon usage in Chiroptera. C. rafinesquii, L. borealis and A. lituratus mitochondrial genomes are 16438 bp, 17048 bp and 16709 bp, respectively. Genome organization and gene arrangements are similar to other bats. Phylogenetic analyses using complete mitochondrial genome sequences support previously established phylogenetic relationships and suggest utility in future studies focusing on the evolutionary aspects of these species. Comprehensive analyses of available bat mitochondrial genomes reveal distinct nucleotide patterns and synonymous codon preferences corresponding to different chiropteran families. These patterns suggest that mutational and selection forces are acting to different extents within Chiroptera and shape their mitochondrial genomes.     

Identification and distribution of sequences having similarity to mitochondrial plasmids in mitochondrial genomes of filamentous fungi

Mitochondrial plasmids are autonomously replicating genetic elements commonly associated with fungal and plant species. Analysis of several plant and fungal mitochondrial genomes has revealed regions that show significant homology to mitochondrial plasmids, suggesting that plasmids have had a long-term association with their mitochondrial hosts. To assess the degree to which plasmids have invaded fungal mitochondrial genomes, BLAST search parameters were modified to identify plasmid sequences within highly AT-rich mtDNAs, and output data were parsed by E value, score, and sequence complexity. High scoring hits were evaluated for the presence of shared repetitive elements and location within plasmids and mtDNAs. Our searches revealed multiple sites of sequence similarity to four distinct plasmids in the wild-type mtDNA of Neurospora crassa, which collectively comprise more than 2% of the mitochondrial genome. Regions of plasmid similarity were not restricted to plasmids known to be associated with senescence, indicating that all mt plasmids can potentially integrate into mitochondrial DNA. Unexpectedly, plasmid-related sequences were found to be clustered in regions that have disproportionately low numbers of PstI palindromic sequences, suggesting that these repetitive elements may play a role in eliminating foreign DNA. A separate class of GC-rich palindromes was identified that appear to be mobile, as indicated by their occurrence within regions of plasmid homology. Sites of sequence similarity to mitochondrial plasmids were also detected in other filamentous fungi, but to a lesser degree. The tools developed here will be useful in assessing the contribution plasmids have made to mitochondrial function and in understanding the co-evolution of mitochondrial plasmids and their hosts.Electronic Supplementary Material Supplementary material is available for this article at     

The complete mitochondrial genome analysis of the tiger (<Emphasis Type="Italic">Panthera tigris</Emphasis>)

The complete mitochondrial genomes of five tiger samples from three subspecies (P. t. sumatrae, P. t. altica, and P. t. tigris) were successfully obtained by using 26 specifically designed Panthera-specific primer sets. The genome organization and gene arrangement of the five tiger samples were similar to each other; however polymorphic tandem repeat sequences were observed in the control region (CR). This led to a difference in the genome lengths obtained from these five samples with an average size of 16,994 bp for the five tiger mitochondrial genomes. The nucleotide base composition was on average as follows: A, 31.8%; T, 27.0%; C, 26.6%; G, 14.6% and exhibited compositional asymmetry. Most of tiger mitochondrial genome characteristics are similar to those of other common vertebrate species; however, some distinctive features were observed in the CR. First, the repetitive sequence 2 (RS 2) contained two repeat units of 80 bp and the first 15 bp of what would be the third repeat motif. The repetitive sequence 3 (RS 3) contained 47–50 repeat motifs of a shorter 8 bp (ACGTAYAC)_n. Second, length heteroplasmy polycystosine (poly-C) stretches was observed at the end of the HV I locus in all tiger samples.     

DNA, chromosomes, and in situ hybridization.

In situ hybridization is a powerful and unique technique that correlates molecular information of a DNA sequence with its physical location along chromosomes and genomes. It thus provides valuable information about physical map position of sequences and often is the only means to determine abundance and distribution of repetitive sequences making up the majority of most genomes. Repeated DNA sequences, composed of units of a few to a thousand base pairs in size, occur in blocks (tandem or satellite repeats) or are dispersed (including transposable elements) throughout the genome. They are often the most variable components of a genome, often being species and, occasionally, chromosome specific. Their variability arises through amplification, diversification and dispersion, as well as homogenization and loss; there is a remarkable correlation of molecular sequence features with chromosomal organization including the length of repeat units, their higher order structures, chromosomal locations, and dispersion mechanisms. Our understanding of the structure, function, organization, and evolution of genomes and their evolving repetitive components enabled many new cytogenetic applications to both medicine and agriculture, particularly in diagnosis and plant breeding.     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

DNA sequence organization in the genomes of three related millet plant species

Summary A major portion of the genomes of three millet species, namely, barn yard millet, fox tail millet and little millet has been shown to consist of interspersed repeat and single copy DNA sequences. The interspersed repetitive DNA sequences are both short (0.15–1.0 kilo base pairs, 62–64% and long (>1.5 kilo base pairs, 36–38%) in barn yard millet and little millet while in fox tail millet, only long interspersed repeats (>1.5 kilo base pairs) are present. The length of the interspersed single copy DNA sequences varies in the range of 1.6–2.6 kilo base pairs in all the three species. The repetitive duplexes isolated after renaturation of 1.5 kilo base pairs and 20 kilo base pairs long DNA fragments exhibit a high thermal stability with Tms either equal to or greater than the corresponding native DNAs. The S1 nuclease resistant repetitive DNA duplexes also are thermally stable and reveal the presence of only 1–2% sequence divergence.The present data on the modes of sequence arrangement in millets substantiates the proposed trend in plants, namely, plants with 1C nuclear DNA content of less than 5 picograms have diverse patterns of sequence organization while those with 1C nuclear DNA content greater than 5 picograms have predominantly a short period interspersion pattern.Abbreviations kbp kilobase pairs NCL Communication No. 3606.     

Sequence and organization of the complete mitochondrial genomes of spotted halibut (Verasper variegatus) and barfin flounder (Verasper moseri).

In this work, the mitochondrial genomes for spotted halibut (Verasper variegatus) and barfin flounder (Verasper moseri) were completely sequenced. The entire mitochondrial genome sequences of the spotted halibut and barfin flounder were 17,273 and 17,588 bp in length, respectively. The organization of the two mitochondrial genomes was similar to those reported from other fish mitochondrial genomes containing 37 genes (2 rRNAs, 22 tRNAs and 13 protein-coding genes) and two non-coding regions (control region (CR) and WANCY region). In the CR, the termination associated sequence (ETAS), six central conserved block (CSB-A,B,C,D,E,F), three conserved sequence blocks (CSB1-3) and a region of 61-bp tandem repeat cluster at the end of CSB-3 were identified by similarity comparison with fishes and other vertebrates. The tandem repeat sequences show polymorphism among the different individuals of the two species. The complete mitochondrial genomes of spotted halibut and barfin flounder should be useful for evolutionary studies of flatfishes and other vertebrate species.     

五种啮齿类动物mtDNA蛋白编码基因序列变化的比较

目的利用本实验室测定的中国地鼠、金黄地鼠和GenBank中田鼠、小鼠、大鼠的线粒体全基因组序列,比较分析五种啮齿类动物的mtDNA蛋白编码基因序列的变异,探讨其分子进化关系。方法将五种动物各自的13个蛋白编码基因分别连接成一个序列,用DNAstar-EditSeq分析软件计算每个序列的碱基长度和组成,计算蛋白编码基因的碱基和氨基酸的差异。以人为外群,基于连接在一起的13个蛋白编码基因的氨基酸序列,用MEGA4.0软件通过最大简约性法(MP)和非加权成对平均数法(UPGMA)构建进化树。结果在五种啮齿动物的13个蛋白基因序列中,A、T、C、G碱基的平均含量为32.4%、29.6%、26.2%和11.9%,中国地鼠mtDNA各蛋白编码序列以及其编码的氨基酸序列与其他物种相比,与金黄地鼠的相应序列差异最小,与大鼠mtDNA各蛋白编码序列以及其编码的氨基酸序列差异较大。分子进化树也显示中国地鼠和金黄地鼠的亲缘关系最近,与小鼠、大鼠存在的差异相对大。结论五种动物的碱基组成的百分比中显示G的相对缺乏,相互之间的进化关系与传统的分类地位基本吻合。