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

Distribution of Blastocystis subtypes isolated from various animal hosts in Thailand

Blastocystis is a parasitic protist of a variety of hosts, including humans. Mapping the distribution of Blastocystis and its genetic variants across different host species can help us understand the epidemiology of this organism and its role in health and disease. This study aimed to identify subtypes of Blastocystis detected in different animal hosts in Thailand. A total of 825 fecal samples belonging to 18 vertebrate orders, 36 families, 68 genera, and 80 species were collected. Of these, 111 specimens were Blastocystis-positive by culture. Seventy-nine samples were subjected to small subunit (SSU) ribosomal DNA amplification by PCR, and reliable subtype data were obtained for 61 specimens. At least 14 subtypes (ST), namely ST1 to ST10, ST14/ST24/ST25 complex, ST23, ST26, and ST29 were detected. In addition, Blastocystis was found in tortoises. ST1 (3.2%) and ST5 (11.5%) were found in pigs, ST2 (1.6%) and ST3 (3.2%) in non-human primates, ST4 (14.7%) in rodents and ruminants, ST6 (4.9%), ST7 (30%), ST9 (1.6%), and ST29 (1.6%) in birds, ST8 (6.6%) in Green peafowl and East Asian Porcupine, and ST10 (4.9%), ST14/ST24/ST25 (9.8%), ST23 (1.6%) and ST26 (1.6%) in ruminants. The sequence recovered from the elongated tortoises (Indotestudo elongata) (3.2%) was phylogenetically placed within the reptilian cluster of Blastocystis, for which no subtype system is available yet. Of note, we did not obtain Blastocystis sequences from any of the many canids and felids sampled in the study, and our data are in support of host specificity of Blastocystis, according to both colonization and subtype distribution.

Evolutionary analysis of cellular reduction and anaerobicity in the hyper-prevalent gut microbe Blastocystis

Blastocystis is the most prevalent microbial eukaryote in the human and animal gut, yet its role as commensal or parasite is still under debate. Blastocystis has clearly undergone evolutionary adaptation to the gut environment and possesses minimal cellular compartmentalization, reduced anaerobic mitochondria, no flagella, and no reported peroxisomes. To address this poorly understood evolutionary transition, we have taken a multi-disciplinary approach to characterize Proteromonas lacertae, the closest canonical stramenopile relative of Blastocystis. Genomic data reveal an abundance of unique genes in P. lacertae but also reductive evolution of the genomic complement in Blastocystis. Comparative genomic analysis sheds light on flagellar evolution, including 37 new candidate components implicated with mastigonemes, the stramenopile morphological hallmark. The P. lacertae membrane-trafficking system (MTS) complement is only slightly more canonical than that of Blastocystis, but notably, we identified that both organisms encode the complete enigmatic endocytic TSET complex, a first for the entire stramenopile lineage. Investigation also details the modulation of mitochondrial composition and metabolism in both P. lacertae and Blastocystis. Unexpectedly, we identify in P. lacertae the most reduced peroxisome-derived organelle reported to date, which leads us to speculate on a mechanism of constraint guiding the dynamics of peroxisome-mitochondrion reductive evolution on the path to anaerobiosis. Overall, these analyses provide a launching point to investigate organellar evolution and reveal in detail the evolutionary path that Blastocystis has taken from a canonical flagellated protist to the hyper-divergent and hyper-prevalent animal and human gut microbe.

Microheliella maris possesses the most gene-rich mitochondrial genome in Diaphoretickes

The mitochondrial genomes are very diverse, but their evolutionary history is unclear due to the lack of efforts to sequence those of protists (unicellular eukaryotes), which cover a major part of the eukaryotic tree. Cryptista comprises cryptophytes, goniomonads, kathablepharids, and Palpitomonas bilix, and their mitochondrial genomes (mt-genomes) are characterized by various gene contents, particularly the presence/absence of an ancestral (bacterial) system for the cytochrome c maturation system. To shed light on mt-genome evolution in Cryptista, we report the complete mt-genome of Microheliella maris, which was recently revealed to branch at the root of Cryptista. The M. maris mt-genome was reconstructed as a circular mapping chromosome of 61.2 kbp with a pair of inverted repeats (12.9 kbp) and appeared to be the most gene-rich among the mt-genomes of the members of Diaphoretickes (a mega-scale eukaryotic assembly including Archaeplastida, Cryptista, Haptista, and SAR) studied so far, carrying 53 protein-coding genes. With this newly sequenced mt-genome, we inferred and discussed the evolution of the mt-genome in Cryptista and Diaphoretickes.

The Genome of the Mitochondrion-Related Organelle in Cepedea longa, a Large Endosymbiotic Opalinid Inhabiting the Recta of Frogs

Mitochondrion-related organelles (MROs) are loosely defined as degenerated mitochondria in anaerobic and microaerophilic lineages. Opalinids are commonly regarded as commensals in the guts of cold-blooded amphibians. It may represent an intermediate adaptation stage between the conventional aerobic mitochondria and derived anaerobic MROs. In the present study, we sequenced and analyzed the MRO genome of Cepedea longa. It has a linear MRO genome with large inverted repeat gene regions at both ends. Compared to Blastocystis and Proteromonas lacertae, the MRO genome of C. longa has a higher G + C content and repeat sequences near the central region. Although three Opalinata species have different morphological characteristics, phylogenetic analyses based on eight concatenated nad genes indicate that they are close relatives. The phylogenetic analysis showed that C. longa clustered with P. lacertae with strong support. The 18S rRNA gene-based phylogeny resolved the Opalinea clade as a sister clade to Karotomorpha, which then further grouped with Proteromonas. The paraphyly of Proteromonadea needs to be verified due to the lack of MRO genomes for key species, such as Karotomorpha, Opalina and Protoopalina. Besides, our dataset and analyses offered slight support for the paraphyly of Bigyra.

Dietary effects on gut microbiota of the mesquite lizard Sceloporus grammicus (Wiegmann, 1828) across different altitudes

BACKGROUND

High-altitude ecosystems are extreme environments that generate specific physiological, morphological, and behavioral adaptations in ectotherms. The shifts in gut microbiota of the ectothermic hosts as an adaptation to environmental changes are still largely unknown. We investigated the food ingested and the bacterial, fungal, and protistan communities in feces of the lizard Sceloporus grammicus inhabiting an altitudinal range using metabarcoding approaches.

RESULTS

The bacterial phyla Bacteroidetes and Firmicutes, and the genera Bacteroides and Parabacteroides dominated the core fecal bacteriome, while Zygomycota and Ascomycota, and the species Basidiobolus ranarum and Basidiobolus magnus dominated the core fecal mycobiome. The diet of S. grammicus included 29 invertebrate families belonging to Arachnida, Chilopoda, and Insecta. The diversity and abundance of its diet decreased sharply at high altitudes, while the abundance of plant material and Agaricomycetes was significantly higher at the highest site. The composition of the fecal microbiota of S. grammicus was different at the three altitudes, but not between females and males. Dietary restriction in S. grammicus at 4150 m might explain the high fecal abundance of Akkermansia and Oscillopira, bacteria characteristic of long fasting periods, while low temperature favored B. magnus. A high proportion of bacterial functions were digestive in S. grammicus at 2600 and 3100, while metabolism of aminoacids, vitamins, and key intermediates of metabolic pathways were higher at 4150 m. Different assemblages of fungal species in the lizard reflect differences in the environments at different elevations. Pathogens were more prevalent at high elevations than at the low ones.

CONCLUSIONS

Limiting food resources at high elevations might oblige S. grammicus to exploit other food resources and its intestinal microbiota have degradative and detoxifying capacities. Sceloporus grammicus might have acquired B. ranarum from the insects infected by the fungus, but its commensal relationship might be established by the quitinolytic capacities of B. ranarum. The mycobiome participate mainly in digestive and degradative functions while the bacteriome in digestive and metabolic functions.

Mitochondrial Genome Evolution and a Novel RNA Editing System in Deep-Branching Heteroloboseids

Discoba (Excavata) is an evolutionarily important group of eukaryotes that includes Jakobida, with the most bacterial-like mitochondrial genomes known, and Euglenozoa, many of which have extensively fragmented mitochondrial genomes. However, little is known about the mitochondrial genomes of Heterolobosea, the third main group of Discoba. Here, we studied two heteroloboseids—an undescribed amoeba “BB2” and Pharyngomonas kirbyi. Phylogenomic analysis revealed that they form a clade that is a sister group to all other Heterolobosea. We characterized the mitochondrial genomes of BB2 and P. kirbyi, which encoded 44 and 48 putative protein-coding genes respectively. Their gene contents were similar to that of Naegleria. In BB2, mitochondrially encoded RNAs were heavily edited, with ∼500 mononucleotide insertion events, mostly guanosines. These insertions always have the same identity as an adjacent nucleotide. Editing occurs in all ribosomal RNAs and protein-coding transcripts except one, and half of the transfer RNAs. Analysis of Illumina deep-sequencing data suggested that this RNA editing is very accurate and efficient, and most likely co-transcriptional. The dissimilarity of this editing process to other RNA editing phenomena in discobids, as well as its apparent absence in P. kirbyi, suggest that this remarkably extensive system of insertional editing evolved independently in the BB2 lineage, after its divergence from the P. kirbyi lineage.

Blastocystis Mitochondrial Genomes Appear to Show Multiple Independent Gains and Losses of Start and Stop Codons

Complete mitochondrion-related organelle (MRO) genomes of several subtypes (STs) of the unicellular stramenopile Blastocystis are presented. Complete conservation of gene content and synteny in gene order is observed across all MRO genomes, comprising 27 protein coding genes, 2 ribosomal RNA genes, and 16 transfer RNA (tRNA) genes. Despite the synteny, differences in the degree of overlap between genes were observed between subtypes and also between isolates within the same subtype. Other notable features include unusual base-pairing mismatches in the predicted secondary structures of some tRNAs. Intriguingly, the rps4 gene in some MRO genomes is missing a start codon and, based on phylogenetic relationships among STs, this loss has happened twice independently. One unidentified open reading frame (orf160) is present in all MRO genomes. However, with the exception of ST4 where the feature has been lost secondarily, orf160 contains variously one or two in-frame stop codons. The overall evidence suggests that both the orf160 and rps4 genes are functional in all STs, but how they are expressed remains unclear.

Mitochondrial Genome of Palpitomonas bilix: Derived Genome Structure and Ancestral System for Cytochrome c Maturation

We here reported the mitochondrial (mt) genome of one of the heterotrophic microeukaryotes related to cryptophytes, Palpitomonas bilix. The P. bilix mt genome was found to be a linear molecule composed of “single copy region” (∼16 kb) and repeat regions (∼30 kb) arranged in an inverse manner at both ends of the genome. Linear mt genomes with large inverted repeats are known for three distantly related eukaryotes (including P. bilix), suggesting that this particular mt genome structure has emerged at least three times in the eukaryotic tree of life. The P. bilix mt genome contains 47 protein-coding genes including ccmA, ccmB, ccmC, and ccmF, which encode protein subunits involved in the system for cytochrome c maturation inherited from a bacterium (System I). We present data indicating that the phylogenetic relatives of P. bilix, namely, cryptophytes, goniomonads, and kathablepharids, utilize an alternative system for cytochrome c maturation, which has most likely emerged during the evolution of eukaryotes (System III). To explain the distribution of Systems I and III in P. bilix and its phylogenetic relatives, two scenarios are possible: (i) System I was replaced by System III on the branch leading to the common ancestor of cryptophytes, goniomonads, and kathablepharids, and (ii) the two systems co-existed in their common ancestor, and lost differentially among the four descendants.

The Earliest Stages of Mitochondrial Adaptation to Low Oxygen Revealed in a Novel Rhizarian

Mitochondria exist on a functional and evolutionary continuum that includes anaerobic mitochondrion-related organelles (MROs), such as hydrogenosomes. Hydrogenosomes lack many classical mitochondrial features, including conspicuous cristae, mtDNA, the tricarboxylic acid (TCA) cycle, and ATP synthesis powered by an electron transport chain (ETC); instead, they produce ATP anaerobically, liberating H₂ and CO₂ gas in the process. However, our understanding of the evolutionary transformation from aerobic mitochondria to various MRO types remains incomplete. Here we describe a novel MRO from a cercomonad (Brevimastigomonas motovehiculus n. sp.; Rhizaria). We have sequenced its 30,608-bp mtDNA and characterized organelle function through a combination of transcriptomic, genomic, and cell biological approaches. B. motovehiculus MROs are metabolically versatile, retaining mitochondrial metabolic pathways, such as a TCA cycle and ETC-driven ATP synthesis, but also possessing hydrogenosomal-type pyruvate metabolism and substrate-level phosphorylation. Notably, the B. motovehiculus ETC is degenerate and appears to be losing cytochrome-based electron transport (complexes III and IV). Furthermore, the F₁F_o ATP synthase (complex V) is unique, with the highly conserved Atpα subunit fragmented into four separate pieces. The B. motovehiculus MRO appears to be in the process of losing aerobic metabolic capacities. Our findings shed light on the transition between organelle types, specifically the early stages of mitochondrial adaptation to anaerobiosis.

Inverted repeats and genome architecture conversions of terrestrial isopods mitochondrial DNA

Mitochondrial DNA (mtDNA) is usually depicted as a circular molecule, however, there is increasing evidence that linearization of mtDNA evolved independently many times in organisms such as fungi, unicellular eukaryotes, and animals. Recent observations in various models with linear mtDNA revealed the presence of conserved inverted repeats (IR) at both ends that, when they become single-stranded, may be able to fold on themselves to create telomeric-hairpins involved in genome architecture conversions. The atypical mtDNA of terrestrial isopods (Crustacea: Oniscidea) composed of linear monomers and circular dimers is an interesting model to study genome architecture conversions. Here, we present the mtDNA control region sequences of two species of the genus Armadillidium: A. vulgare and A. pelagicum. All features of arthropods mtDNA control regions are present (origin of replication, poly-T stretch, GA and TA-rich blocks and one variable domain), plus a conserved IR. This IR can potentially fold into a hairpin structure and is present in two different orientations among the A. vulgare populations: either in one sense or in its reverse complement. This polymorphism, also observed in a single individual (heteroplasmy), might be a signature of genome architecture conversions from linear to circular monomeric mtDNA via successive opening and closing of the molecules.

Sequencing and characterization of the complete mitochondrial genomes of three Pneumocystis species provide new insights into divergence between human and rodent Pneumocystis

Pneumocystis jirovecii is an important opportunistic pathogen associated with AIDS and other immunodeficient conditions. Currently, very little is known about its nuclear and mitochondrial genomes. In this study, we sequenced the complete mitochondrial genome (mtDNA) of this organism and its closely related species Pneumocystis carinii and Pneumocystis murina by a combination of sequencing technologies. Our study shows that P. carinii and P. murina mtDNA share a nearly identical number and order of genes in a linear configuration, whereas P. jirovecii has a circular mtDNA containing nearly the same set of genes but in a different order. Detailed studies of the mtDNA terminal structures of P. murina and P. carinii suggest a unique replication mechanism for linear mtDNA. Phylogenetic analysis supports a close association of Pneumocystis species with Taphrina, Saitoella, and Schizosaccharomyces, and divergence within Pneumocystis species, with P. murina and P. carinii being more closely related to each other than either is to P. jirovecii. Comparative analysis of four complete P. jirovecii mtDNA sequences in this study and previously reported mtDNA sequences for diagnosing and genotyping suggests that the current diagnostic and typing methods can be improved using the complete mtDNA data. The availability of the complete P. jirovecii mtDNA also opens the possibility of identifying new therapeutic targets.

Mitochondrial evolution

Viewed through the lens of the genome it contains, the mitochondrion is of unquestioned bacterial ancestry, originating from within the bacterial phylum α-Proteobacteria (Alphaproteobacteria). Accordingly, the endosymbiont hypothesis--the idea that the mitochondrion evolved from a bacterial progenitor via symbiosis within an essentially eukaryotic host cell--has assumed the status of a theory. Yet mitochondrial genome evolution has taken radically different pathways in diverse eukaryotic lineages, and the organelle itself is increasingly viewed as a genetic and functional mosaic, with the bulk of the mitochondrial proteome having an evolutionary origin outside Alphaproteobacteria. New data continue to reshape our views regarding mitochondrial evolution, particularly raising the question of whether the mitochondrion originated after the eukaryotic cell arose, as assumed in the classical endosymbiont hypothesis, or whether this organelle had its beginning at the same time as the cell containing it.

First complete mitochondrial genome sequence from a box jellyfish reveals a highly fragmented linear architecture and insights into telomere evolution

Animal mitochondrial DNAs (mtDNAs) are typically single circular chromosomes, with the exception of those from medusozoan cnidarians (jellyfish and hydroids), which are linear and sometimes fragmented. Most medusozoans have linear monomeric or linear bipartite mitochondrial genomes, but preliminary data have suggested that box jellyfish (cubozoans) have mtDNAs that consist of many linear chromosomes. Here, we present the complete mtDNA sequence from the winged box jellyfish Alatina moseri (the first from a cubozoan). This genome contains unprecedented levels of fragmentation: 18 unique genes distributed over eight 2.9- to 4.6-kb linear chromosomes. The telomeres are identical within and between chromosomes, and recombination between subtelomeric sequences has led to many genes initiating or terminating with sequences from other genes (the most extreme case being 150 nt of a ribosomal RNA containing the 5′ end of nad2), providing evidence for a gene conversion–based model of telomere evolution. The silent-site nucleotide variation within the A. moseri mtDNA is among the highest observed from a eukaryotic genome and may be associated with elevated rates of recombination.

The organellar genome and metabolic potential of the hydrogen-producing mitochondrion of Nyctotherus ovalis

It is generally accepted that hydrogenosomes (hydrogen-producing organelles) evolved from a mitochondrial ancestor. However, until recently, only indirect evidence for this hypothesis was available. Here, we present the almost complete genome of the hydrogen-producing mitochondrion of the anaerobic ciliate Nyctotherus ovalis and show that, except for the notable absence of genes encoding electron transport chain components of Complexes III, IV, and V, it has a gene content similar to the mitochondrial genomes of aerobic ciliates. Analysis of the genome of the hydrogen-producing mitochondrion, in combination with that of more than 9,000 genomic DNA and cDNA sequences, allows a preliminary reconstruction of the organellar metabolism. The sequence data indicate that N. ovalis possesses hydrogen-producing mitochondria that have a truncated, two step (Complex I and II) electron transport chain that uses fumarate as electron acceptor. In addition, components of an extensive protein network for the metabolism of amino acids, defense against oxidative stress, mitochondrial protein synthesis, mitochondrial protein import and processing, and transport of metabolites across the mitochondrial membrane were identified. Genes for MPV17 and ACN9, two hypothetical proteins linked to mitochondrial disease in humans, were also found. The inferred metabolism is remarkably similar to the organellar metabolism of the phylogenetically distant anaerobic Stramenopile Blastocystis. Notably, the Blastocystis organelle and that of the related flagellate Proteromonas lacertae also lack genes encoding components of Complexes III, IV, and V. Thus, our data show that the hydrogenosomes of N. ovalis are highly specialized hydrogen-producing mitochondria.

Complex repeat structures and novel features in the mitochondrial genomes of the diatoms Phaeodactylum tricornutum and Thalassiosira pseudonana

The mitochondrial genome of the raphid pennate diatom Phaeodactylum tricornutum has several novel features compared with the mitochondrial genomes of the centric diatom Thalassiosira pseudonana and the araphid pennate diatom Synedra acus. It is almost double the size (77,356 bp) due to a 35,454 bp sequence block consisting of an elaborate combination of direct repeats, making it the largest stramenopile (heterokont) mitochondrial genome known. In addition, the cox1 gene has a +1 translational frameshift involving Pro codons CCC and CCT, the first translational frameshift to be detected in an algal mitochondrial genome. The nad9 and rps14 genes are fused by the insertion of an in-frame sequence and cotranscribed. The nad11 gene is split into two parts corresponding to the FeS and molybdate-binding domains, but both parts are still on the mitochondrial genome, in contrast to the brown algae where the second domain appears to have been transferred to the nucleus. In contrast to P. tricornutum, the repeat region of T. pseudonana consists of a much smaller 4790 bp string of almost identical double-hairpin elements, evidence of slipped-strand mispairing and active gene conversion. The diatom mitochondrial genomes have undergone considerable gene rearrangement since the three lineages of diatoms diverged, but all three have kept their repeat regions segregated from their relatively compact coding regions.

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.