Size and structural variations of mitochondrial DNA

Mitogenomic architecture and evolution of the soil ciliates Colpoda

Colpoda are cosmopolitan unicellular eukaryotes primarily inhabiting soil and benefiting plant growth, but they remain one of the least understood taxa in genetics and genomics within the realm of ciliated protozoa. Here, we investigate the architecture of de novo assembled mitogenomes of six Colpoda species, using long-read sequencing and involving 36 newly isolated natural strains in total. The mitogenome sizes span from 43 to 63 kbp and typically contain 28-33 protein-coding genes. They possess a linear structure with variable telomeres and central repeats, with one Colpoda elliotti strain isolated from Tibet harboring the longest telomeres among all studied ciliates. Phylogenomic analyses reveal that Colpoda species started to diverge more than 326 million years ago, eventually evolving into two distinct groups. Collinearity analyses also reveal significant genomic divergences and a lack of long collinear blocks. One of the most notable features is the exceptionally high level of gene rearrangements between mitochondrial genomes of different Colpoda species, dominated by gene loss events. Population-level mitogenomic analysis on natural strains also demonstrates high sequence divergence, regardless of geographic distance, but the gene order remains highly conserved within species, offering a new species identification criterion for Colpoda species. Furthermore, we identified underlying heteroplasmic sites in the majority of strains of three Colpoda species, albeit without a discernible recombination signal to account for this heteroplasmy. This comprehensive study systematically unveils the mitogenomic structure and evolution of these ancient and ecologically significant Colpoda ciliates, thus laying the groundwork for a deeper understanding of the evolution of unicellular eukaryotes.IMPORTANCEColpoda, one of the most widespread ciliated protozoa in soil, are poorly understood in regard to their genetics and evolution. Our research revealed extreme mitochondrial gene rearrangements dominated by gene loss events, potentially leading to the streamlining of Colpoda mitogenomes. Surprisingly, while interspecific rearrangements abound, our population-level mitogenomic study revealed a conserved gene order within species, offering a potential new identification criterion. Phylogenomic analysis traced their lineage over 326 million years, revealing two distinct groups. Substantial genomic divergence might be associated with the lack of extended collinear blocks and relaxed purifying selection. This study systematically reveals Colpoda ciliate mitogenome structures and evolution, providing insights into the survival and evolution of these vital soil microorganisms.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

Characterization and Comparative Analyses of Mitochondrial Genomes in Single-Celled Eukaryotes to Shed Light on the Diversity and Evolution of Linear Molecular Architecture

Determination and comparisons of complete mitochondrial genomes (mitogenomes) are important to understand the origin and evolution of mitochondria. Mitogenomes of unicellular protists are particularly informative in this regard because they are gene-rich and display high structural diversity. Ciliates are a highly diverse assemblage of protists and their mitogenomes (linear structure with high A+T content in general) were amongst the first from protists to be characterized and have provided important insights into mitogenome evolution. Here, we report novel mitogenome sequences from three representatives (Strombidium sp., Strombidium cf. sulcatum, and Halteria grandinella) in two dominant ciliate lineages. Comparative and phylogenetic analyses of newly sequenced and previously published ciliate mitogenomes were performed and revealed a number of important insights. We found that the mitogenomes of these three species are linear molecules capped with telomeric repeats that differ greatly among known species. The genomes studied here are highly syntenic, but larger in size and more gene-rich than those of other groups. They also all share an AT-rich tandem repeat region which may serve as the replication origin and modulate initiation of bidirectional transcription. More generally we identified a split version of ccmf, a cytochrome c maturation-related gene that might be a derived character uniting taxa in the subclasses Hypotrichia and Euplotia. Finally, our mitogenome comparisons and phylogenetic analyses support to reclassify Halteria grandinella from the subclass Oligotrichia to the subclass Hypotrichia. These results add to the growing literature on the unique features of ciliate mitogenomes, shedding light on the diversity and evolution of their linear molecular architecture.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

Incorporating mitogenome sequencing into integrative taxonomy: The multidisciplinary redescription of the ciliate Thuricola similis (Peritrichia, Vaginicolidae) provides new insights into the evolutionary relationships among Oligohymenophorea subclasses

The evolutionary relationships among Oligohymenophorea subclasses are under debate as the phylogenomic analysis using a large dataset of nuclear coding genes is significantly different to the 18S rDNA phylogeny, and it is unfortunately not stable within and across different published studies. In addition to nuclear genes, the faster-evolving mitochondrial genes have also shown the ability to solve phylogenetic problems in many ciliated taxa. However, due to the paucity of mitochondrial data, the corresponding work is scarce, let alone the phylogenomic analysis based on mitochondrial gene dataset. In this work, we presented the characterization on Thuricola similis Bock, 1963, a loricate peritrich (Oligohymenophorea), incorporating mitogenome sequencing into integrative taxonomy. As the first mitogenome for the subclass Peritrichia, it is linear, 38,802 bp long, and contains two rRNAs, 12 tRNAs, and 43 open reading frames (ORFs). As a peculiarity, it includes a central repeated region composed of tandemly repeated A-T rich units working as a bi-transcriptional start. Moreover, taking this opportunity, the phylogenomic analyses based on a set of mitochondrial genes were also performed, revealing that T. similis, as a representative of Peritrichia subclass, branches basally to other three Oligohymenophorea subclasses, namely Hymenostomatia, Peniculia, and Scuticociliatia. Evolutionary relationships among those Oligohymenophorea subclasses were discussed, also in the light of recent phylogenomic reconstructions based on a set of nuclear genes. Besides, as a little-known species, T. similis was also redescribed and neotypified based on data from two populations collected from wastewater treatment plants (WWTPs) in Brazil and Italy, by means of integrative methods (i.e., living observation, silver staining methods, scanning and transmission electron microscopy, and 18S rDNA phylogeny). After emended diagnosis, it is characterized by: (1) the sewage habitat; (2) the lorica with a single valve and small undulations; (3) the 7-22 µm-long inner stalk; and (4) the presence of only a single postciliary microtubule on the left side of the aciliferous row in the haplokinety. Among Vaginicolidae family, our 18S rRNA gene-based phylogenetic analysis revealed that Thuricola and Cothurnia are monophyletic genera, and Vaginicola could be a polyphyletic genus.

Morphology, ultrastructure, genomics, and phylogeny of Euplotes vanleeuwenhoeki sp. nov. and its ultra-reduced endosymbiont "Candidatus Pinguicoccus supinus" sp. nov

Taxonomy is the science of defining and naming groups of biological organisms based on shared characteristics and, more recently, on evolutionary relationships. With the birth of novel genomics/bioinformatics techniques and the increasing interest in microbiome studies, a further advance of taxonomic discipline appears not only possible but highly desirable. The present work proposes a new approach to modern taxonomy, consisting in the inclusion of novel descriptors in the organism characterization: (1) the presence of associated microorganisms (e.g.: symbionts, microbiome), (2) the mitochondrial genome of the host, (3) the symbiont genome. This approach aims to provide a deeper comprehension of the evolutionary/ecological dimensions of organisms since their very first description. Particularly interesting, are those complexes formed by the host plus associated microorganisms, that in the present study we refer to as "holobionts". We illustrate this approach through the description of the ciliate Euplotes vanleeuwenhoeki sp. nov. and its bacterial endosymbiont "Candidatus Pinguicoccus supinus" gen. nov., sp. nov. The endosymbiont possesses an extremely reduced genome (~ 163 kbp); intriguingly, this suggests a high integration between host and symbiont.

Ciliate mitoribosome illuminates evolutionary steps of mitochondrial translation

To understand the steps involved in the evolution of translation, we used Tetrahymena thermophila, a ciliate with high coding capacity of the mitochondrial genome, as the model organism and characterized its mitochondrial ribosome (mitoribosome) using cryo-EM. The structure of the mitoribosome reveals an assembly of 94-ribosomal proteins and four-rRNAs with an additional protein mass of ~700 kDa on the small subunit, while the large subunit lacks 5S rRNA. The structure also shows that the small subunit head is constrained, tRNA binding sites are formed by mitochondria-specific protein elements, conserved protein bS1 is excluded, and bacterial RNA polymerase binding site is blocked. We provide evidence for anintrinsic protein targeting system through visualization of mitochondria-specific mL105 by the exit tunnel that would facilitate the recruitment of a nascent polypeptide. Functional protein uS3m is encoded by three complementary genes from the nucleus and mitochondrion, establishing a link between genetic drift and mitochondrial translation. Finally, we reannotated nine open reading frames in the mitochondrial genome that code for mitoribosomal proteins.

Co-evolution in the Jungle: From Leafcutter Ant Colonies to Chromosomal Ends

Biological entities are multicomponent systems where each part is directly or indirectly dependent on the others. In effect, a change in a single component might have a consequence on the functioning of its partners, thus affecting the fitness of the entire system. In this article, we provide a few examples of such complex biological systems, ranging from ant colonies to a population of amino acids within a single-polypeptide chain. Based on these examples, we discuss one of the central and still challenging questions in biology: how do such multicomponent consortia co-evolve? More specifically, we ask how telomeres, nucleo-protein complexes protecting the integrity of linear DNA chromosomes, originated from the ancestral organisms having circular genomes and thus not dealing with end-replication and end-protection problems. Using the examples of rapidly evolving topologies of mitochondrial genomes in eukaryotic microorganisms, we show what means of co-evolution were employed to accommodate various types of telomere-maintenance mechanisms in mitochondria. We also describe an unprecedented runaway evolution of telomeric repeats in nuclei of ascomycetous yeasts accompanied by co-evolution of telomere-associated proteins. We propose several scenarios derived from research on telomeres and supported by other studies from various fields of biology, while emphasizing that the relevant answers are still not in sight. It is this uncertainty and a lack of a detailed roadmap that makes the journey through the jungle of biological systems still exciting and worth undertaking.

Comparative Studies on the Polymorphism and Copy Number Variation of mtSSU rDNA in Ciliates (Protista, Ciliophora): Implications for Phylogenetic, Environmental, and Ecological Research

While nuclear small subunit ribosomal DNA (nSSU rDNA) is the most commonly-used gene marker in studying phylogeny, ecology, abundance, and biodiversity of microbial eukaryotes, mitochondrial small subunit ribosomal DNA (mtSSU rDNA) provides an alternative. Recently, both copy number variation and sequence variation of nSSU rDNA have been demonstrated for diverse organisms, which can contribute to misinterpretation of microbiome data. Given this, we explore patterns for mtSSU rDNA among 13 selected ciliates (representing five classes), a major component of microbial eukaryotes, estimating copy number and sequence variation and comparing to that of nSSU rDNA. Our study reveals: (1) mtSSU rDNA copy number variation is substantially lower than that for nSSU rDNA; (2) mtSSU rDNA copy number ranges from 1.0 × 104 to 8.1 × 105; (3) a most common sequence of mtSSU rDNA is also found in each cell; (4) the sequence variation of mtSSU rDNA are mainly indels in poly A/T regions, and only half of species have sequence variation, which is fewer than that for nSSU rDNA; and (5) the polymorphisms between haplotypes of mtSSU rDNA would not influence the phylogenetic topology. Together, these data provide more insights into mtSSU rDNA as a powerful marker especially for microbial ecology studies.

The linear mitochondrial genome of the quarantine chytrid Synchytrium endobioticum; insights into the evolution and recent history of an obligate biotrophic plant pathogen

BACKGROUND

Chytridiomycota species (chytrids) belong to a basal lineage in the fungal kingdom. Inhabiting terrestrial and aquatic environments, most are free-living saprophytes but several species cause important diseases: e.g. Batrachochytrium dendrobatidis, responsible for worldwide amphibian decline; and Synchytrium endobioticum, causing potato wart disease. S. endobioticum has an obligate biotrophic lifestyle and isolates can be further characterized as pathotypes based on their virulence on a differential set of potato cultivars. Quarantine measures have been implemented globally to control the disease and prevent its spread. We used a comparative approach using chytrid mitogenomes to determine taxonomical relationships and to gain insights into the evolution and recent history of introductions of this plant pathogen.

RESULTS

We assembled and annotated the complete mitochondrial genome of 30 S. endobioticum isolates and generated mitochondrial genomes for five additional chytrid species. The mitochondrial genome of S. endobioticum is linear with terminal inverted repeats which was validated by tailing and PCR amplifying the telomeric ends. Surprisingly, no conservation in organisation and orientation of mitochondrial genes was observed among the Chytridiomycota except for S. endobioticum and its sister species Synchytrium microbalum. However, the mitochondrial genome of S. microbalum is circular and comprises only a third of the 72.9 Kbp found for S. endobioticum suggesting recent linearization and expansion. Four mitochondrial lineages were identified in the S. endobioticum mitochondrial genomes. Several pathotypes occur in different lineages, suggesting that these have emerged independently. In addition, variations for polymorphic sites in the mitochondrial genome of individual isolates were observed demonstrating that S. endobioticum isolates represent a community of different genotypes. Such communities were shown to be complex and stable over time, but we also demonstrate that the use of semi-resistant potato cultivars triggers a rapid shift in the mitochondrial haplotype associated with increased virulence.

CONCLUSIONS

Mitochondrial genomic variation shows that S. endobioticum has been introduced into Europe multiple times, that several pathotypes emerged multiple times, and that isolates represent communities of different genotypes. Our study represents the most comprehensive dataset of chytrid mitogenomes, which provides new insights into the extraordinary dynamics and evolution of mitochondrial genomes involving linearization, expansion and reshuffling.

If the cap fits, wear it: an overview of telomeric structures over evolution

Genome organization into linear chromosomes likely represents an important evolutionary innovation that has permitted the development of the sexual life cycle; this process has consequently advanced nuclear expansion and increased complexity of eukaryotic genomes. Chromosome linearity, however, poses a major challenge to the internal cellular machinery. The need to efficiently recognize and repair DNA double-strand breaks that occur as a consequence of DNA damage presents a constant threat to native chromosome ends known as telomeres. In this review, we present a comparative survey of various solutions to the end protection problem, maintaining an emphasis on DNA structure. This begins with telomeric structures derived from a subset of prokaryotes, mitochondria, and viruses, and will progress into the typical telomere structure exhibited by higher organisms containing TTAGG-like tandem sequences. We next examine non-canonical telomeres from Drosophila melanogaster, which comprise arrays of retrotransposons. Finally, we discuss telomeric structures in evolution and possible switches between canonical and non-canonical solutions to chromosome end protection.

Diversity of mitochondrial genome organization

In this review, we discuss types of mitochondrial genome structural organization (architecture), which includes the following characteristic features: size and the shape of DNA molecule, number of encoded genes, presence of cryptogenes, and editing of primary transcripts.

The Oxytricha trifallax mitochondrial genome

The Oxytricha trifallax mitochondrial genome contains the largest sequenced ciliate mitochondrial chromosome (~70 kb) plus a ~5-kb linear plasmid bearing mitochondrial telomeres. We identify two new ciliate split genes (rps3 and nad2) as well as four new mitochondrial genes (ribosomal small subunit protein genes: rps- 2, 7, 8, 10), previously undetected in ciliates due to their extreme divergence. The increased size of the Oxytricha mitochondrial genome relative to other ciliates is primarily a consequence of terminal expansions, rather than the retention of ancestral mitochondrial genes. Successive segmental duplications, visible in one of the two Oxytricha mitochondrial subterminal regions, appear to have contributed to the genome expansion. Consistent with pseudogene formation and decay, the subtermini possess shorter, more loosely packed open reading frames than the remainder of the genome. The mitochondrial plasmid shares a 251-bp region with 82% identity to the mitochondrial chromosome, suggesting that it most likely integrated into the chromosome at least once. This region on the chromosome is also close to the end of the most terminal member of a series of duplications, hinting at a possible association between the plasmid and the duplications. The presence of mitochondrial telomeres on the mitochondrial plasmid suggests that such plasmids may be a vehicle for lateral transfer of telomeric sequences between mitochondrial genomes. We conjecture that the extreme divergence observed in ciliate mitochondrial genomes may be due, in part, to repeated invasions by relatively error-prone DNA polymerase-bearing mobile elements.

Evolution of linear chromosomes and multipartite genomes in yeast mitochondria

Mitochondrial genome diversity in closely related species provides an excellent platform for investigation of chromosome architecture and its evolution by means of comparative genomics. In this study, we determined the complete mitochondrial DNA sequences of eight Candida species and analyzed their molecular architectures. Our survey revealed a puzzling variability of genome architecture, including circular- and linear-mapping and multipartite linear forms. We propose that the arrangement of large inverted repeats identified in these genomes plays a crucial role in alterations of their molecular architectures. In specific arrangements, the inverted repeats appear to function as resolution elements, allowing genome conversion among different topologies, eventually leading to genome fragmentation into multiple linear DNA molecules. We suggest that molecular transactions generating linear mitochondrial DNA molecules with defined telomeric structures may parallel the evolutionary emergence of linear chromosomes and multipartite genomes in general and may provide clues for the origin of telomeres and pathways implicated in their maintenance.

A linear molecule with two large inverted repeats: the mitochondrial genome of the stramenopile Proteromonas lacertae

Mitochondrial evolution has given rise to a complex array of organelles, ranging from classical aerobic mitochondria to mitochondrial remnants known as hydrogenosomes and mitosomes. The latter are found in anaerobic eukaryotes, and these highly derived organelles often retain only scant evidence of their mitochondrial origins. Intermediate evolutionary stages have also been reported as facultatively or even strictly anaerobic mitochondria, and hydrogenosomes that still retain some mitochondrial features. However, the diversity among these organelles with transitional features remains rather unclear and barely studied. Here, we report the sequence, structure, and gene content of the mitochondrial DNA of the anaerobic stramenopile Proteromonas lacertae. It has a linear genome with a unique central region flanked by two identical large inverted repeats containing numerous genes and "telomeres" with short inverted repeats. Comparison with the organelle genome of the strictly anaerobic human parasite Blastocystis reveals that, despite the close similarity of the sequences, features such as the genome structure display striking differences. It remains unclear whether the virtually identical gene repertoires are the result of convergence or descent.

Structure and Evolution of the Atypical Mitochondrial Genome of Armadillidium vulgare (Isopoda, Crustacea)

The crustacean isopod Armadillidium vulgare is characterized by an unusual approximately 42-kb-long mitochondrial genome consisting of two molecules co-occurring in mitochondria: a circular approximately 28-kb dimer formed by two approximately 14-kb monomers fused in opposite polarities and a linear approximately 14-kb monomer. Here we determined the nucleotide sequence of the fundamental monomeric unit of A. vulgare mitochondrial genome, to gain new insight into its structure and evolution. Our results suggest that the junction zone between monomers of the dimer structure is located in or near the control region. Direct sequencing indicated that the nucleotide sequences of the different monomer units are virtually identical. This suggests that gene conversion and/or replication processes play an important role in shaping nucleotide sequence variation in this mitochondrial genome. The only heteroplasmic site we identified predicts an alloacceptor tRNA change from tRNA(Ala) to tRNA(Val). Therefore, in A. vulgare, tRNA(Ala) and tRNA(Val) are found at the same locus in different monomers, ensuring that both tRNAs are present in mitochondria. The presence of this heteroplasmic site in all sequenced individuals suggests that the polymorphism is selectively maintained, probably because of the necessity of both tRNAs for maintaining proper mitochondrial functions. Thus, our results provide empirical evidence for the tRNA gene recruitment model of tRNA evolution. Moreover, interspecific comparisons showed that the A. vulgare mitochondrial gene order is highly derived compared to the putative ancestral arthropod type. By contrast, an overall high conservation of mitochondrial gene order is observed within crustacean isopods.

Mitochondrial genome diversity: evolution of the molecular architecture and replication strategy

Mitochondrial genomes in organisms from diverse phylogenetic groups vary in both size and molecular form. Although the types of mitochondrial genome appear very dissimilar, several lines of evidence argue that they do not differ radically. This would imply that interconversion between different types of mitochondrial genome might have occurred via relatively simple mechanisms. We exemplify this scenario on patterns accompanying evolution of mitochondrial telomeres. We propose that mitochondrial telomeres are derived from mobile elements (transposons or plasmids) that invaded mitochondria, integrated into circular or polydisperse linear mitochondrial DNAs (mtDNAs) and subsequently enabled precise resolution of the linear genophore. Simply, the selfish elements generated a problem - how to maintain the ends of a linear DNA - and, at the same time, made themselves essential by providing its solution. This scenario implies that insertion or deletion of such resolution elements may represent relatively simple routes for interconversion between different forms of the mitochondrial genome.

Complete sequence of the mitochondrial genome of Tetrahymena thermophila and comparative methods for identifying highly divergent genes

The complete sequence of the mitochondrial genome of Tetrahymena thermophila has been determined and compared with the mitochondrial genome of Tetrahymena pyriformis. The sequence similarity clearly indicates homology of the entire T.thermophila and T.pyriformis mitochondrial genomes. The T.thermophila genome is very compact, most of the intergenic regions are short (only three are longer than 63 bp) and comprise only 3.8% of the genome. The nad9 gene is tandemly duplicated in T.thermophila. Long terminal inverted repeats and the nad9 genes are undergoing concerted evolution. There are 55 putative genes: three ribosomal RNA genes, eight transfer RNA genes, 22 proteins with putatively assigned functions and 22 additional open reading frames of unknown function. In order to extend indications of homology beyond amino acid sequence similarity we have examined a number of physico-chemical properties of the mitochondrial proteins, including theoretical pI, molecular weight and particularly the predicted transmembrane spanning regions. This approach has allowed us to identify homologs to ymf58 (nad4L), ymf62 (nad6) and ymf60 (rpl6).

The curious history of yeast mitochondrial DNA

Forty years ago, soon after yeast mitochondrial DNA (mtDNA) was recognized, some animal versions of mtDNA were shown to comprise circular molecules. Supporting an idea that mitochondria had evolved from bacteria, this finding generated a dogmatic belief that yeast mtDNA was also circular, and the endless linear molecules actually observed in yeast were regarded as broken circles. This concept persisted for 30 years and has distorted our understanding of the true nature of the molecule.

Mitochondrial DNA from the liverwort Marchantia polymorpha: circularly permuted linear molecules, head-to-tail concatemers, and a 5' protein

Mapping predicts that the mitochondrial genome of the liverwort Marchantia polymorpha exists as a circular molecule, although nearly all the mitochondrial DNA (mtDNA) is found as genome-sized and multigenomic molecules in linear and branched form. We used restriction enzymes with one recognition site per genome, end-specific exonucleases and pulsed-field gel electrophoresis (PFGE) to analyze the arrangement of genomic units and the terminal structure of the molecules. We find a head-to-tail arrangement in the concatemers and circular permutation in both the monomeric and multigenomic molecules. The termini contain covalently bound protein at the 5' end and an open (unblocked) 3' end. We find that the standard in-gel procedure used to prepare large DNA molecules for PFGE may introduce extraction artifacts leading to erroneous conclusions about the termini. These artifacts can be reduced by omitting high salt (high EDTA) and protease during mitochondrial lysis. Our results suggest that the mtDNA may use a T4 phage-like mechanism of replication and that the linear molecules may be due to strand breaks mediated by type II topoisomerase.

Complete sequence of the mitochondrial genome of Tetrahymena pyriformis and comparison with Paramecium aurelia mitochondrial DNA

We report the complete nucleotide sequence of the Tetrahymena pyriformis mitochondrial genome and a comparison of its gene content and organization with that of Paramecium aurelia mtDNA. T. pyriformis mtDNA is a linear molecule of 47,172 bp (78.7 % A+T) excluding telomeric sequences (identical tandem repeats of 31 bp at each end of the genome). In addition to genes encoding the previously described bipartite small and large subunit rRNAs, the T. pyriformis mitochondrial genome contains 21 protein-coding genes that are clearly homologous to genes of defined function in other mtDNAs, including one (yejR) that specifies a component of a cytochrome c biogenesis pathway. As well, T. pyriformis mtDNA contains 22 open reading frames of unknown function larger than 60 codons, potentially specifying proteins ranging in size from 74 to 1386 amino acid residues. A total of 13 of these open reading frames ("ciliate-specific") are found in P. aurelia mtDNA, whereas the remaining nine appear to be unique to T. pyriformis; however, of the latter, five are positionally equivalent and of similar size in the two ciliate mitochondrial genomes, suggesting they may also be homologous, even though this is not evident from sequence comparisons. Only eight tRNA genes encoding seven distinct tRNAs are found in T. pyriformis mtDNA, formally confirming a long-standing proposal that most T. pyriformis mitochondrial tRNAs are nucleus-encoded species imported from the cytosol. Atypical features of mitochondrial gene organization and expression in T. pyriformis mtDNA include split and rearranged large subunit rRNA genes, as well as a split nad1 gene (encoding subunit 1 of NADH dehydrogenase of respiratory complex I) whose two segments are located on and transcribed from opposite strands, as is also the case in P. aurelia. Gene content and arrangement are very similar in T. pyriformis and P. aurelia mtDNAs, the two differing by a limited number of duplication, inversion and rearrangement events. Phylogenetic analyses using concatenated sequences of several mtDNA-encoded proteins provide high bootstrap support for the monophyly of alveolates (ciliates, dinoflagellates and apicomplexans) and slime molds.

Linear mitochondrial genomes: 30 years down the line

At variance with the earlier belief that mitochondrial genomes are represented by circular DNA molecules, a large number of organisms have been found to carry linear mitochondrial DNA. Studies of linear mitochondrial genomes might provide a novel view on the evolutionary history of organelle genomes and contribute to delineating mechanisms of maintenance and functioning of telomeres. Because linear mitochondrial DNA is present in a number of human pathogens, its replication mechanisms might become a target for drugs that would not interfere with replication of human circular mitochondrial DNA.

Mitochondrial import of only one of three nuclear-encoded glutamine tRNAs in Tetrahymena thermophila

The mitochondrial genome of Tetrahymena does not appear to encode enough tRNAs to perform mitochondrial protein synthesis. It has therefore been proposed that nuclear-encoded tRNAs are imported into the mitochondria. T.thermophila has three major glutamine tRNAs: tRNA(Gln)(UUG), tRNA(Gln)(UUA) and tRNA(Gln)(CUA). Each of these tRNAs functions in cytosolic translation. However, due to differences between the Tetrahymena nuclear and mitochondrial genetic codes, only tRNA(Gln)(UUG) has the capacity to function in mitochondrial translation as well. Here we show that approximately 10-20% of the cellular complement of tRNA(Gln)(UUG) is present in mitochondrial RNA fractions, compared with 1% or less for the other two glutamine tRNAs. Furthermore, this glutamine tRNA is encoded only by a family of nuclear genes, the sequences of several of which are presented. Finally, when marked versions of tRNA(Gln)(UUG) and tRNA(Gln)(UUA) flanked by identical sequences are expressed in the macronucleus, only the former undergoes mitochondrial import; thus sequences within tRNA(Gln)(UUG) direct import. Because tRNA(Gln)(UUG) is a constituent of mitochondrial RNA fractions and is encoded only by nuclear genes, and because ectopically expressed tRNA(Gln)(UUG) fractionates with mitochondria like its endogenous counterpart, we conclude that it is an imported tRNA in T.thermophila.

Linear mitochondrial genome organization in vivo in the genus Pythium

Pulsed-field gel electrophoresis (PFGE) of isolates of Pythium oligandrum with linear mitochondrial genomes revealed a distinct band in ethidium bromide-stained gels similar in size to values estimated by restriction mapping of mitochondrial DNA (mtDNA). Southern analysis confirmed that these bands were mtDNA and indicated that linear genomes were present in unit-length size as well as multimers. Isolates of this species with circular mtDNA restriction maps also had low levels of linear mono- and multimers visualized by Southern analysis of PFGE gels. Examination of 17 additional species revealed similar results; three species had distinct linear mtDNA bands in ethidium bromide-stained gels while the remainder had linear mono- and multi-mers in lower amounts detected only by Southern analysis. Sequence analysis of an isolate of P. oligandrum with a primarily circular mitochondrial genomic map and a low amount of linear molecules revealed that the small unique region of the circular map (which corresponded to the terminal region of linear genomes) was flanked by palindromic intrastrand complementary sequences separated by a unique 194-bp sequence. Sequences with similarity to ATPase9 coding regions from other organisms were located adjacent to this region. Sequences with similarity to mitochondrial origins of replication and autonomously replicating sequences were also located in this region: their potential involvement in the generation of linear molecules is discussed.

Mitochondrial DNA of Chlamydomonas reinhardtii: the structure of the ends of the linear 15.8-kb genome suggests mechanisms for DNA replication

The mitochondrial genome of Chlamydomonas reinhardtii is a linear double-stranded DNA of 15.8 kb. With the exception of the termini its DNA sequence has been published. Here we describe the unique structure of the two termini determined from cloned fragments or, for the very terminal sequences, by the Maxam and Gilbert method after 5' labeling of uncloned terminal fragments. The 15.8-kb DNA is characterized by terminal inverted repeats of 531 or 532 bp in length including long 3' extensions. The 3' single-stranded extensions of the left and right ends are non-complementary, identical in sequence, and comprise 39 to 41 nucleotides. Remarkably, the linear genome possesses in addition an internal 86-bp repeat of the two outermost sequences. The unusual structure of the 15.8-kb DNA termini is compared with those of other linear mitochondrial DNAs. Possible mechanisms of 15.8-kb DNA replication are discussed.

Linear mitochondrial DNAs of yeasts: frequency of occurrence and general features

In most yeast species, the mitochondrial DNA (mtDNA) has been reported to be a circular molecule. However, two cases of linear mtDNA with specific termini have previously been described. We examined the frequency of occurrence of linear forms of mtDNA among yeasts by pulsed-field gel electrophoresis. Among the 58 species from the genera Pichia and Williopsis that we examined, linear mtDNA was found with unexpectedly high frequency. Thirteen species contained a linear mtDNA, as confirmed by restriction mapping, and labeling, and electron microscopy. The mtDNAs from Pichia pijperi, Williopsis mrakii, and P. jadinii were studied in detail. In each case, the left and right terminal fragments shared homologous sequences. Between the terminal repeats, the order of mitochondrial genes was the same in all of the linear mtDNAs examined, despite a large variation of the genome size. This constancy of gene order is in contrast with the great variation of gene arrangement in circular mitochondrial genomes of yeasts. The coding sequences determined on several genes were highly homologous to those of the circular mtDNAs, suggesting that these two forms of mtDNA are not of distant origins.

Linear mitochondrial DNAs of yeasts: frequency of occurrence and general features

In most yeast species, the mitochondrial DNA (mtDNA) has been reported to be a circular molecule. However, two cases of linear mtDNA with specific termini have previously been described. We examined the frequency of occurrence of linear forms of mtDNA among yeasts by pulsed-field gel electrophoresis. Among the 58 species from the genera Pichia and Williopsis that we examined, linear mtDNA was found with unexpectedly high frequency. Thirteen species contained a linear mtDNA, as confirmed by restriction mapping, and labeling, and electron microscopy. The mtDNAs from Pichia pijperi, Williopsis mrakii, and P. jadinii were studied in detail. In each case, the left and right terminal fragments shared homologous sequences. Between the terminal repeats, the order of mitochondrial genes was the same in all of the linear mtDNAs examined, despite a large variation of the genome size. This constancy of gene order is in contrast with the great variation of gene arrangement in circular mitochondrial genomes of yeasts. The coding sequences determined on several genes were highly homologous to those of the circular mtDNAs, suggesting that these two forms of mtDNA are not of distant origins.

Amplification of tandemly repeated origin control sequences confers a replication advantage on rDNA replicons in Tetrahymena thermophila

The macronuclear rRNA genes (rDNA) in the ciliate Tetrahymena thermophila are normally palindromic linear replicons, containing two copies of the replication origin region in inverted orientation. A circular plasmid containing a single Tetrahymena rRNA gene (one half palindrome) joined to a tandem repeat of a 1.9-kilobase (kb) rDNA segment encompassing the rDNA replication origin and known replication control elements was used to transform Tetrahymena macronuclei by microinjection. This plasmid was shown previously to have a replication advantage over the rDNA allele of the recipient cell strain (G.-L. Yu and E. H. Blackburn, Proc. Natl. Acad. Sci. USA 86:8487-8491, 1990). During vegetative cell divisions, the circular and palindromic rDNAs were rapidly replaced by novel, successively longer linear rDNAs that eventually contained up to 30 tandem 1.9-kb repeats, resulting from homologous but unequal crossovers between the 1.9-kb repeats. We present evidence to show that increasing the number of copies of the replication control regions increases the replicative advantage of the rDNA, the first such situation for a cellular nuclear replicon in a eucaryote.

Amplification of tandemly repeated origin control sequences confers a replication advantage on rDNA replicons in Tetrahymena thermophila

The macronuclear rRNA genes (rDNA) in the ciliate Tetrahymena thermophila are normally palindromic linear replicons, containing two copies of the replication origin region in inverted orientation. A circular plasmid containing a single Tetrahymena rRNA gene (one half palindrome) joined to a tandem repeat of a 1.9-kilobase (kb) rDNA segment encompassing the rDNA replication origin and known replication control elements was used to transform Tetrahymena macronuclei by microinjection. This plasmid was shown previously to have a replication advantage over the rDNA allele of the recipient cell strain (G.-L. Yu and E. H. Blackburn, Proc. Natl. Acad. Sci. USA 86:8487-8491, 1990). During vegetative cell divisions, the circular and palindromic rDNAs were rapidly replaced by novel, successively longer linear rDNAs that eventually contained up to 30 tandem 1.9-kb repeats, resulting from homologous but unequal crossovers between the 1.9-kb repeats. We present evidence to show that increasing the number of copies of the replication control regions increases the replicative advantage of the rDNA, the first such situation for a cellular nuclear replicon in a eucaryote.

Mitochondrial DNA of Physarum polycephalum: physical mapping, cloning and transcription mapping

Mitochondrial DNA (mtDNA) has been isolated from four strains of Physarum polycephalum and a restriction site map has been determined using nine restriction enzymes. The restriction site maps of the four strains are similar but each strain is distinguished by insertions, deletions and restriction enzyme site polymorphisms. The sum of the restriction fragments gives mitochondrial genome sizes which vary from about 56 kb to 62 kb. In all four strains the composite map of the restriction enzyme sites for the mtDNA is circular. Knowledge of the restriction enzyme map has enabled cloning of mtDNA fragments representing the entire mtDNA of strain M3. The cloned fragments have been used to create a transcription map of the mtDNA.

Telomeric repeats of Tetrahymena malaccensis mitochondrial DNA: a multimodal distribution that fluctuates erratically during growth

The linear mitochondrial DNA (mtDNA) of Tetrahymena malaccensis has tandem 52-base-pair repeats at its telomeres. The mtDNA has a multimodal distribution of telomeres. Different groups in the distribution have different numbers of telomeric repeats. The standard deviation of the size of each end group is independent of the mean size of the end group. The two sides of the mtDNA have different multimodal distributions of repeats. Cloned cell lines have multimodal distributions of mtDNA telomeres distinct from that of the original cell line. The number of telomere end groups and the average size of the end groups change in an erratic fashion as the cells are passaged and do not reach a stable equilibrium distribution in 185 generations. We propose that the mean size of a telomere end group and the size distribution of an end group are independently regulated. The system controlling the average size of end groups may be defective in T. malaccensis, since a closely related species (T. thermophila) does not have a multimodal distribution of mtDNA telomeres. T. hyperangularis, which has different telomeric repeats on each side of its mtDNA, has a multimodal distribution of mtDNA telomeres on only one side, suggesting that the mechanism controlling the average number of repeats in an end group can be sequence specific. These mitochondrial telomeres provide a new example of the more general phenomenon of expansion and contraction of arrays of repeated sequences seen, for example, with simple-sequence "satellite" DNAs; however, the mitochondrial telomeres change on a very short time scale.

Telomeric repeats of Tetrahymena malaccensis mitochondrial DNA: a multimodal distribution that fluctuates erratically during growth

The linear mitochondrial DNA (mtDNA) of Tetrahymena malaccensis has tandem 52-base-pair repeats at its telomeres. The mtDNA has a multimodal distribution of telomeres. Different groups in the distribution have different numbers of telomeric repeats. The standard deviation of the size of each end group is independent of the mean size of the end group. The two sides of the mtDNA have different multimodal distributions of repeats. Cloned cell lines have multimodal distributions of mtDNA telomeres distinct from that of the original cell line. The number of telomere end groups and the average size of the end groups change in an erratic fashion as the cells are passaged and do not reach a stable equilibrium distribution in 185 generations. We propose that the mean size of a telomere end group and the size distribution of an end group are independently regulated. The system controlling the average size of end groups may be defective in T. malaccensis, since a closely related species (T. thermophila) does not have a multimodal distribution of mtDNA telomeres. T. hyperangularis, which has different telomeric repeats on each side of its mtDNA, has a multimodal distribution of mtDNA telomeres on only one side, suggesting that the mechanism controlling the average number of repeats in an end group can be sequence specific. These mitochondrial telomeres provide a new example of the more general phenomenon of expansion and contraction of arrays of repeated sequences seen, for example, with simple-sequence "satellite" DNAs; however, the mitochondrial telomeres change on a very short time scale.

Linear mitochondrial deoxyribonucleic acid from the yeast Hansenula mrakii

The mitochondrial deoxyribonucleic acid (mtDNA) from a petite-negative yeast, Hansenula mrakii, was studied. A linear restriction map was constructed with 11 restriction enzymes. The linearity of the genome was confirmed by direct end labeling of the molecule, followed by restriction analysis. The molecular weight of the DNA was found to be 55,000 base pairs. This is the first linear mtDNA found in yeast species. Using specific gene probes obtained from Saccharomyces cerevisiae mtDNA, we have constructed a gene map of H. mrakii mtDNA. The arrangement of genes in this linear genome was very different from the circular mtDNA of other known yeasts.

Phylogenetic relationships and altered genome structures among Tetrahymena mitochondrial DNAs

Tetrahymena thermophila mitochondrial DNA is a linear molecule with two tRNAs, large subunit beta (LSU beta) rRNA (21S rRNA) and LSU alpha rRNA (5.8S-like RNA) encoded near each terminus. The DNA sequence of approximately 550 bp of this region was determined in six species of Tetrahymena. In three species the LSU beta rRNA and tRNA(leu) genes were not present on one end of the DNA, demonstrating a mitochondrial genome organization different from that of T. thermophila. The DNA sequence of the LSU alpha rRNA was used to construct a mitochondrial phylogenetic tree, which was found to be topologically equivalent to a phylogenetic tree based on nuclear small subunit rRNA sequences (Sogin et al. (1986) EMBO J. 5, 3625-3630). The mitochondrial rRNA gene was found to accumulate base-pair substitutions considerably faster than the nuclear rRNA gene, the rate difference being similar to that observed for mammals.

The cytochrome oxidase subunit I gene of Tetrahymena: a 57 amino acid NH2-terminal extension and a 108 amino acid insert

The gene sequence for cytochrome oxidase subunit I (COI) in the ciliate Tetrahymena mitochondrial DNA has been determined and shown to be coded by the same strand as codes the genes (in order) for 14S rRNA, tRNA(trp), tRNA(glu), 21S rRNA, tRNA(leu) and tRNA(met). The predicted protein has 698 amino acids, including an NH2-terminal 57 amino acid extension and a 108 amino acid insert originally found in Paramecium COI. These extension and insert segments are not highly hydrophobic but are relatively rich in lysine, arginine and serine. In analogy with the presequence of nuclear-encoded mitochondrial proteins, they might function as a transmembrane signal. The remaining polypeptide segments show a hydrophobicity characteristic of membrane spanning proteins. TCOI shows a 64% amino acid identity with Paramecium COI but less than a 38% amino acid conservation with human COI. The Tetrahymena mitochondrial code is analogous with the mammalian mitochondrial code; but differs from the Tetrahymena nuclear genetic code; TGA is exclusively translated as tryptophan; ATA is used as an initiation codon probably for methionine, and TAA as a stop codon; the arginine codons (CGN) are not used. The use of the leucine codon TTA in TCOI is contradictory to the codon recognition pattern previously obtained from the isolated tRNA(leu) isoacceptors recognizing only the CUN codons, but consistent with the tRNA(leu) (anticodon UAA) gene encoded in the genome. The reason for this inconsistency has not been resolved.

A fine restriction map of the linear mitochondrial DNA of Tetrahymena pyriformis: genome size, map locations of rRNA and tRNA genes, terminal inversion repeat, and restriction site polymorphism

A fine restriction map of the linear mitochondrial DNA of Tetrahymena pyriformis strain ST is presented. 1. Based on agarose gel electrophoresis data together with limited nucleotide sequences available on some restriction fragments, we estimate the actual size of this genome to be about 55,000 base pairs. 2. Seven tRNA gene locations have been assigned, which are scattered along the genome length. Six of these locations encode the genes for tRNA(phe), tRNA(his), tRNA(trp), and tRNA(glu), and the duplicate tRNA(tyr) genes which are located at the inverted terminal repeat segments. The tRNA gene(s) encoded in one location has not been identified. We have not yet found the tRNA(leu) and tRNA(met) genes, which were previously shown to be encoded in the genome (Chiu et al. 1974; Suyama 1982). 3. We have mapped the 14S rRNA gene by sequencing the 170 bp segment of EcoRI fragment 8 and by aligning its sequence with E. coli 16S rRNA. From our recent complete sequence data the gene size was found to be about 1,650 bp, which is unexpectedly large for the 14S rRNA which has an estimated size of 1,300 bp. The 14S rRNA is probably a cleavage product of the larger primary transcript of which 200-300 bases of the 5' end are missing. 4. The duplicate copies of the 21S rRNA gene at the terminal duplication inversion segments were analyzed. ClaI fragment 7 (1,500 bp) corresponds in sequence from base position 850 to 2,390 of the 20S rRNA gene of Paramecium mitochondrial DNA (Seilhamer et al. 1984b). The 21S gene is approximately 2,500 bp long. The presence of some restriction site polymorphism is apparent in this segment. 5. Each of the 21S gene copies precedes the tRNA(tyr) gene, but the space flanking one tRNA(tyr) gene differs in size and restriction sites from the space flanking another tRNA(tyr) gene. Thus, this space corresponds to the segment of an imperfect match in the terminal duplication inversion of Goldbach et al. (1978a). 6. Saccharomyces cerevisiae mitochondrial probes including Cob, ATPase VI and IX, and cytochrome oxidase I gene sequences, 21S and 15S rRNAs, and mouse mitochondrial DNA showed no significant hybridization with any restriction fragments of Tetrahymena mitochondrial DNA. The results are in accordance with an extensive sequence divergence previously found in the Tetrahymena mitochondrial genome (Goldbach et al. 1977).

A review of the structure and organization of the mitochondrial genome of higher plants

The structure, organisation and functions of the mitochondrial genome of most groups of eukaryotic organisms are known to varying degrees with the notable exception of higher plants. With the recent completion of physical mapping studies on the Brassica campestris (Chinese cabbage, turnip) and Zea mays (maize) mitochondrial genomes, many of the apparent problems of plant mitochondrial genome structure can now be answered. In this manuscript I review the literature relating to the physical observations on plant mitochondrial DNA (mtDNA) and assess the data in relation to our current understanding of mitochondrial genome structure.

Organization and closing of mitochondrial deoxyribonucleic acid from Paramecium tetraaurelia and Paramecium primaurelia

Previously we showed that the mitochondrial deoxyribonucleic acid (DNA) from Paramecium aurelia consists of a linear genome and that replication of this genome is initiated at one terminus and proceeds unidirectionally to the other terminus. Analyses of mitochondria from four closely related species (1, 4, 5, and 7) indicated that the species 1, 5, and 7 DNAs are essentially completely homologous but that the species 4 mitochondrial DNA is only 40 to 50% homologous with that from species 1. The major regions of homology are those containing the genes for ribosomal ribonucleic acid (RNA). To understand the replication and organization of the linear mitochondrial genome better, we compared species 1 (Paramecium primaurelia) and 4 (Paramecium tetraaurelia) DNAs with regard to restriction fragment mapping and homology between initiation regions; we also identified the sites of the genes for ribosomal RNA. In general, the structures of the species 1 and 4 mitochondrial genomes were quite similar. Each ribosomal RNA gene was present in one copy per genome, with the large ribosomal RNA gene located near the terminal region of replication and the small ribosomal RNA gene located more centrally. These two genes were separated by about 10 kilobases in the species 1 genome and by about 12 kilobases in the species 4 genome. In contrast to our previous findings, by using nonstringent hybridization conditions we detected homology between the species 1 and 4 DNA fragments containing the initiation regions. We constructed recombinant DNA clones for many fragments, especially those containing the initiation region and the ribosomal RNA genes. We also constructed restriction enzyme maps for six enzymes for both P. primaurelia and P. tetraaurelia.

Organization and closing of mitochondrial deoxyribonucleic acid from Paramecium tetraaurelia and Paramecium primaurelia

Previously we showed that the mitochondrial deoxyribonucleic acid (DNA) from Paramecium aurelia consists of a linear genome and that replication of this genome is initiated at one terminus and proceeds unidirectionally to the other terminus. Analyses of mitochondria from four closely related species (1, 4, 5, and 7) indicated that the species 1, 5, and 7 DNAs are essentially completely homologous but that the species 4 mitochondrial DNA is only 40 to 50% homologous with that from species 1. The major regions of homology are those containing the genes for ribosomal ribonucleic acid (RNA). To understand the replication and organization of the linear mitochondrial genome better, we compared species 1 (Paramecium primaurelia) and 4 (Paramecium tetraaurelia) DNAs with regard to restriction fragment mapping and homology between initiation regions; we also identified the sites of the genes for ribosomal RNA. In general, the structures of the species 1 and 4 mitochondrial genomes were quite similar. Each ribosomal RNA gene was present in one copy per genome, with the large ribosomal RNA gene located near the terminal region of replication and the small ribosomal RNA gene located more centrally. These two genes were separated by about 10 kilobases in the species 1 genome and by about 12 kilobases in the species 4 genome. In contrast to our previous findings, by using nonstringent hybridization conditions we detected homology between the species 1 and 4 DNA fragments containing the initiation regions. We constructed recombinant DNA clones for many fragments, especially those containing the initiation region and the ribosomal RNA genes. We also constructed restriction enzyme maps for six enzymes for both P. primaurelia and P. tetraaurelia.