Mitochondrial Evolution Research Articles

The Complete Mitogenome of Apostasia fujianica Y.Li & S.Lan and Comparative Analysis of Mitogenomes across Orchidaceae.

Apostasia fujianica belongs to the genus Apostasia and is part of the basal lineage in the phylogenetic tree of the Orchidaceae. Currently, there are only ten reported complete mitochondrial genomes in orchids, which greatly hinders the understanding of mitochondrial evolution in Orchidaceae. Therefore, we assembled and annotated the mitochondrial genome of A. fujianica, which has a length of 573,612 bp and a GC content of 44.5%. We annotated a total of 44 genes, including 30 protein-coding genes, 12 tRNA genes, and two rRNA genes. We also performed relative synonymous codon usage (RSCU) analysis, repeat sequence analysis, intergenomic transfer (IGT) analysis, and Ka/Ks analysis for A. fujianica and conducted RNA editing site analysis on the mitochondrial genomes of eight orchid species. We found that most protein-coding genes are under purifying selection, but nad6 is under positive selection, with a Ka/Ks value of 1.35. During the IGT event in A. fujianica's mitogenome, the trnN-GUU, trnD-GUC, trnW-CCA, trnP-UGG, and psaJ genes were identified as having transferred from the plastid to the mitochondrion. Compared to other monocots, the family Orchidaceae appears to have lost the rpl10, rpl14, sdh3, and sdh4 genes. Additionally, to further elucidate the evolutionary relationships among monocots, we constructed a phylogenetic tree based on the complete mitogenomes of monocots. Our study results provide valuable data on the mitogenome of A. fujianica and lay the groundwork for future research on genetic variation, evolutionary relationships, and breeding of Orchidaceae.

De novo assembly and comprehensive analysis of the mitochondrial genome of Taxus wallichiana reveals different repeats mediate recombination to generate multiple conformations

Taxus plants are the exclusive source of paclitaxel, an anticancer drug with significant medicinal and economic value. Interspecies hybridization and gene introgression during evolution have obscured distinctions among Taxus species, complicating their phylogenetic classification. While the chloroplast genome of Taxus wallichiana, a widely distributed species in China, has been sequenced, its mitochondrial genome (mitogenome) remains uncharacterized.We sequenced and assembled the T. wallichiana mitogenome using BGI short reads and Nanopore long reads, facilitating comparisons with other gymnosperm mitogenomes. The T. wallichiana mitogenome spanning 469,949 bp, predominantly forms a circular configuration with a GC content of 50.51%, supplemented by 3 minor configurations mediated by one pair of LRs and two pairs of IntRs. It includes 32 protein-coding genes, 7 tRNA genes, and 3 rRNA genes, several of which exist in multiple copies.We detailed the mitogenome's structure, codon usage, RNA editing, and sequence migration between organelles, constructing a phylogenetic tree to elucidate evolutionary relationships. Unlike typical gymnosperm mitochondria, T. wallichiana shows no evidence of mitochondrial-plastid DNA transfer (MTPT), highlighting its unique genomic architecture. Synteny analysis indicated extensive genomic rearrangements in T. wallichiana, likely driven by recombination among abundant repetitive sequences. This study offers a high-quality T. wallichiana mitogenome, enhancing our understanding of gymnosperm mitochondrial evolution and supporting further cultivation and utilization of Taxus species.

Accelerated mitochondrial evolution and asymmetric fitness of hybrids contribute to the persistence of Helix thessalica in the Helix pomatia range.

Interbreeding and introgression between recently diverged species is common. However, the processes that prevent these species from merging where they co-occur are not well understood. We studied the mechanisms that allowed an isolated group of populations of the snail Helix thessalica to persist within the range of the related Helix pomatia despite high gene flow. Using genomic cline analysis, we found that the nuclear gene flow between the two taxa across the mosaic hybrid zone was not different from that expected under neutral admixture, but that the exchange of mtDNA was asymmetric. Tests showed that there is relaxed selection in the mitochondrial genome of H. thessalica and that the substitution rate is elevated compared to that of H. pomatia. A lack of hybrids that combine the mtDNA of H. thessalica with a mainly (>46%) H. pomatia genomic background indicates that the nuclear-encoded mitochondrial proteins of H. pomatia are not well adapted to the more rapidly evolving proteins and RNAs encoded by the mitochondrion of H. thessalica. The presumed reduction of fitness of hybrids with the fast-evolving mtDNA of H. thessalica and a high H. pomatia ancestry, similar to 'Darwin's Corollary to Haldane's rule', resulted in a relative loss of H. pomatia nuclear ancestry compared to H. thessalica ancestry in the hybrid zone. This probably prevents the H. thessalica populations from merging quickly with the surrounding H. pomatia populations and supports the hypothesis that incompatibilities between rapidly evolving mitochondrial genes and nuclear genes contribute to speciation.

Joint inference of cell lineage and mitochondrial evolution from single-cell sequencing data.

Eukaryotic cells contain organelles called mitochondria that have their own genome. Most cells contain thousands of mitochondria which replicate, even in nondividing cells, by means of a relatively error-prone process resulting in somatic mutations in their genome. Because of the higher mutation rate compared to the nuclear genome, mitochondrial mutations have been used to track cellular lineage, particularly using single-cell sequencing that measures mitochondrial mutations in individual cells. However, existing methods to infer the cell lineage tree from mitochondrial mutations do not model "heteroplasmy," which is the presence of multiple mitochondrial clones with distinct sets of mutations in an individual cell. Single-cell sequencing data thus provide a mixture of the mitochondrial clones in individual cells, with the ancestral relationships between these clones described by a mitochondrial clone tree. While deconvolution of somatic mutations from a mixture of evolutionarily related genomes has been extensively studied in the context of bulk sequencing of cancer tumor samples, the problem of mitochondrial deconvolution has the additional constraint that the mitochondrial clone tree must be concordant with the cell lineage tree. We formalize the problem of inferring a concordant pair of a mitochondrial clone tree and a cell lineage tree from single-cell sequencing data as the Nested Perfect Phylogeny Mixture (NPPM) problem. We derive a combinatorial characterization of the solutions to the NPPM problem, and formulate an algorithm, MERLIN, to solve this problem exactly using a mixed integer linear program. We show on simulated data that MERLIN outperforms existing methods that do not model mitochondrial heteroplasmy nor the concordance between the mitochondrial clone tree and the cell lineage tree. We use MERLIN to analyze single-cell whole-genome sequencing data of 5220 cells of a gastric cancer cell line and show that MERLIN infers a more biologically plausible cell lineage tree and mitochondrial clone tree compared to existing methods. https://github.com/raphael-group/MERLIN.

Complete mitochondrial genome of an oleaginous microalga Vischeria punctata (Eustigmatophyceae: Chlorobotryaceae) and phylogenetic analysis

Vischeria punctata, as first described by Vischer in 1945, is a member of the family Chlorobotryaceae, within the order Eustigmatales. This species is recognized for its potential as a source of biofuels and other high-value products. In the present investigation, the whole genome of V. punctata was sequenced utilizing the Illumina HiSeq 4000 platform, enabling the assembly and annotation of its complete mitochondrial genome. The resulting circular genome spans 41,528 base pairs (bp) with a guanine-cytosine (GC) content of 27.3%. This genome encompasses 36 protein-coding genes, alongside 28 transfer RNA (tRNA), and three ribosomal RNA (rRNA) genes. The evolutionary trajectory of V. punctata was further explored by constructing a phylogenetic tree derived from the mitochondrial 33 gene dataset of 16 Ochrophyta species. Comparative analysis reveals that V. punctata bears closer ties to Vischeria sp. CAUP Q202 than to Vischeria stellata strain SAG 33.83, suggesting shared evolutionary pathways and phenotypic traits. This investigation constitutes the inaugural study into the mitochondrial evolution and phylogenetic patterning of the mitogenome in V. punctata. The outcomes from this research bolster our understanding of the genetic diversity and evolutionary processes within the class Eustigmatophyceae. In particular, the mitochondrial genome of V. punctata serves as a valuable resource in elucidating these aspects.

First complete mitochondrial genome of the Alashan ground squirrel (Spermophilus alashanicus) (Rodentia: Sciuridae) from Ningxia, China

The Alashan ground squirrel (Spermophilus alashanicus) is primarily distributed in the regions of Inner Mongolia and Ningxia, China. In this study, we present the first complete mitochondrial genome of S. alashanicus. The genome spans 16,464 base pairs and comprises 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and a single control region with a marked AT bias. The overall GC content is 35.4%. Phylogenetic analyses indicate that S. alashanicus clusters are closely associated with S. dauricus. This comprehensive characterization of the S. alashanicus mitochondrial genome serves as a foundational resource for future studies on mitochondrial evolution, species identification, population genomics, and phylogenetics.

The pace of mitochondrial molecular evolution varies with seasonal migration distance.

Animals that engage in long-distance seasonal migration experience strong selective pressures on their metabolic performance and life history, with potential consequences for molecular evolution. Species with slow life histories typically show lower rates of synonymous substitution (dS) than "fast" species. Previous research suggests long-distance seasonal migrants have a slower life history strategy than short-distance migrants, raising the possibility that rates of molecular evolution may covary with migration distance. Additionally, long-distance migrants may face strong selection on metabolically-important mitochondrial genes due to their long-distance flights. Using over 1,000 mitochondrial genomes, we assessed the relationship between migration distance and mitochondrial molecular evolution in 39 boreal-breeding migratory bird species. We show that migration distance correlates negatively with dS, suggesting that the slow life history associated with long-distance migration is reflected in rates of molecular evolution. Mitochondrial genes in every study species exhibited evidence of purifying selection, but the strength of selection was greater in short-distance migrants, contrary to our predictions. This result may indicate effects of selection for cold tolerance on mitochondrial evolution among species overwintering at high latitudes. Our study demonstrates that the pervasive correlation between life history and molecular evolutionary rates exists in the context of differential adaptations to seasonality.

ATP synthase-associated coiled-coil-helix-coiled-coil-helix (CHCH) domain-containing proteins are critical for mitochondrial function in Toxoplasma gondii.

Members of the coiled-coil-helix-coiled-coil-helix (CHCH) domain protein family are transported into the mitochondrial intermembrane space, where they play important roles in the biogenesis and function of the organelle. Unexpectedly, the ATP synthase of the apicomplexan Toxoplasma gondii harbors CHCH domain-containing subunits of unknown function. As no other ATP synthase studied to date contains this class of proteins, characterizing their function will be of broad interest to the fields of molecular parasitology and mitochondrial evolution. Here, we demonstrate that that two T. gondii ATP synthase subunits containing CHCH domains are required for parasite survival and for stability and function of the ATP synthase. We also show that knockdown disrupts multiple aspects of the mitochondrial morphology of T. gondii and that mutation of key residues in the CHCH domains caused mis-localization of the proteins. This work provides insight into the unique features of the apicomplexan ATP synthase, which could help to develop therapeutic interventions against this parasite and other apicomplexans, such as the malaria-causing parasite Plasmodium falciparum.

Assembly and comparative analysis of the complete mitochondrial genome of Isopyrum anemonoides (Ranunculaceae).

Ranunculaceae is a large family of angiosperms comprising 2500 known species-a few with medicinal and ornamental values. Despite this, only two mitochondrial genomes (mitogenomes) of the family have been released in GenBank. Isopyrum anemonoides is a medicinal plant belonging to the family Ranunculaceae, and its chloroplast genome has recently been reported; however, its mitogenome remains unexplored. In this study, we assembled and analyzed the complete mitochondrial genome of I. anemonoides and performed a comparative analysis against different Ranunculaceae species, reconstructing the phylogenetic framework of Isopyrum. The circular mitogenome of I. anemonoides has a length of 206,722 bp, with a nucleotide composition of A (26.4%), T (26.4%), C (23.6%), and G (23.6%), and contains 62 genes, comprising 37 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and three ribosomal RNA (rRNA) genes. Abundantly interspersed repetitive and simple sequence repeat (SSR) loci were detected in the I. anemonoides mitogenome, with tetranucleotide repeats accounting for the highest proportion of SSRs. By detecting gene migration, we observed gene exchange between the chloroplast and mitogenome in I. anemonoides, including six intact tRNA genes, six PCG fragments, and fragments from two rRNA genes. Comparative mitogenome analysis of three Ranunculaceae species indicated that the PCG contents were conserved and the GC contents were similar. Selective pressure analysis revealed that only two genes (nad1 and rpl5) were under positive selection during their evolution in Ranunculales, and two specific RNA editing sites (atp6 and mttB) were detected in the I. anemonoides mitogenome. Moreover, a phylogenetic analysis based on the mitogenomes of I. anemonoides and the other 15 taxa accurately reflected the evolutionary and taxonomic status of I. anemonoides. Overall, this study provides new insights into the genetics, systematics, and evolution of mitochondrial evolution in Ranunculaceae, particularly I. anemonoides.

Worldwide Population Genomics Reveal Long-Term Stability of the Mitochondrial Genome Architecture in a Keystone Marine Plant.

Mitochondrial genomes (mitogenomes) of flowering plants are composed of multiple chromosomes. Recombination within and between the mitochondrial chromosomes may generate diverse DNA molecules termed isoforms. The isoform copy number and composition can be dynamic within and among individual plants due to uneven replication and homologous recombination. Nonetheless, despite their functional importance, the level of mitogenome conservation within species remains understudied. Whether the ontogenetic variation translates to evolution of mitogenome composition over generations is currently unknown. Here we show that the mitogenome composition of the seagrass Zostera marina is conserved among worldwide populations that diverged ca. 350,000 years ago. Using long-read sequencing, we characterized the Z. marina mitochondrial genome and inferred the repertoire of recombination-induced configurations. To characterize the mitochondrial genome architecture worldwide and study its evolution, we examined the mitogenome in Z. marina meristematic region sampled in 16 populations from the Pacific and Atlantic oceans. Our results reveal a striking similarity in the isoform relative copy number, indicating a high conservation of the mitogenome composition among distantly related populations and within the plant germline, despite a notable variability during individual ontogenesis. Our study supplies a link between observations of dynamic mitogenomes at the level of plant individuals and long-term mitochondrial evolution.

Mitochondrial evolution: Gene shuffling, endosymbiosis, and signaling.

Genes for cardiolipin and ceramide synthesis occur in some alphaproteobacterial genomes. They shed light on mitochondrial origin and signaling in the first eukaryotic cells.

Maturases and Group II Introns in the Mitochondrial Genomes of the Deepest Jakobid Branch.

Ophirinina is a recently described suborder of jakobid protists (Excavata) with only one described species to date, Ophirina amphinema. Despite the acquisition and analysis of massive transcriptomic and mitogenomic sequence data from O. amphinema, its phylogenetic position among excavates remained inconclusive, branching as sister group either to all Jakobida or to all Discoba. From a morphological perspective, it has not only several typical jakobid features but also unusual traits for this group, including the morphology of mitochondrial cristae (sac-shaped to flattened-curved cristae) and the presence of two flagellar vanes. In this study, we have isolated, morphologically characterized, and sequenced genome and transcriptome data of two new Ophirinina species: Ophirina chinija sp. nov. and Agogonia voluta gen. et sp. nov. Ophirina chinija differs from O. amphinema in having rounded cell ends, subapically emerging flagella and a posterior cell protrusion. The much more distantly related A. voluta has several unique ultrastructural characteristics, including sac-shaped mitochondrial cristae and a complex "B" fiber. Phylogenomic analyses with a large conserved-marker dataset supported the monophyly of Ophirina and Agogonia within the Ophirinina and, more importantly, resolved the conflicting position of ophirinids as the sister clade to all other jakobids. The characterization of the mitochondrial genomes showed that Agogonia differs from all known gene-rich jakobid mitogenomes by the presence of two group II introns and their corresponding maturase protein genes. A phylogenetic analysis of the diversity of known maturases confirmed that the Agogonia proteins are highly divergent from each other and define distant families among the prokaryotic and eukaryotic maturases. This opens the intriguing possibility that, compared to other jakobids, Ophirinina may have retained additional mitochondrial elements that may help to understand the early diversification of eukaryotes and the evolution of mitochondria.

The first two mitochondrial genomes from Apiotrichum reveal mitochondrial evolution and different taxonomic assignment of Trichosporonales

Apiotrichum is a diverse anamorphic basidiomycetous yeast genus, and its mitogenome characterization has not been revealed. In this study, we assembled two Apiotrichum mitogenomes and compared them with mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. The mitogenomes of Apiotrichum gracile and A. gamsii comprised circular DNA molecules, with sizes of 34,648 bp and 38,096 bp, respectively. Intronic regions were found contributed the most to the size expansion of A. gamsii mitogenome. Comparative mitogenomic analysis revealed that 6.85–38.89% of nucleotides varied between tRNAs shared by the two Apiotrichum mitogenomes. The GC content of all core PCGs in A. gamsii was lower than that of A. gracile, with an average low value of 4.97%. The rps3 gene differentiated the most among Agaricomycotina, Pucciniomycotina and Ustilaginomycotina species, while nad4L gene was the most conserved in evolution. The Ka/Ks values for cob and rps3 genes were > 1, indicating the two genes may be subjected to positive selection in Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. Frequent intron loss/gain events and potential intron transfer events have been detected in evolution of Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. We further detected large-scale gene rearrangements between the 19 mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina, and fifteen of the 17 mitochondrial genes shared by Apiotrichum varied in gene arrangements. Phylogenetic analyses based on maximum likelihood and Bayesian inference methods using a combined mitochondrial gene dataset revealed different taxonomic assignment of two Apiotrichum species, wherein A. gamsii had a more closely relationship with Trichosporon asahii. This study served as the first report on mitogenomes from the genus Apiotrichum, which promotes the understanding of evolution, genomics, and phylogeny of Apiotrichum.

High-throughput sequencing yields a complete mitochondrial genome of the rice thrips, Stenchaetothrips biformis (Thysanoptera: Thripidae)

Rice thrips, Stenchaetothrips biformis (Bagnall, 1913), are one of the destructive pests of rice. Here, the complete mitochondrial genome of S. biformis was sequenced using high-throughput sequencing. The mitogenome is 15,359 bp long with an A + T content of 76.94%, which contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNAs), 2 ribosomal RNA genes (rRNAs) and 2 putative control regions (CRs). The phylogenetic analysis showed that S. biformis is closely related to Thrips imaginis and Thrips palmi. This new mitochondrial genome data can be better used to provide a basis for studies of the mitochondrial evolution of Thysanoptera.

Comparative structural analysis on the mitochondrial DNAs from various strains of Lentinula edodes.

The evolution of mitochondria through variations in mitochondrial DNA (mtDNA) is one of the intriguing questions in eukaryotic cells. In order to assess the causes of the variations in mitochondria, the mtDNAs of the 21 strains of Lentinula edodes were assembled for this study, and analyzed together with four published mtDNA sequences. The mtDNAs were within the sizes of 117 kb ~ 122 kb. The gene number was observed consistent except for two mtDNAs, which carry a duplicated trnG1-trnG2 unit or a putative gene deletion. The size variation was largely attributed to the number of introns, repeated sequences, transposable elements (TEs), and plasmid-related sequences. Intron loss and gain were found from cox1, rnl, and rns of three mtDNAs. Loss of two introns in cox1 of KY217797.1 reduced its size by 2.7 kb, making it the smallest cox1 gene (8.4 kb) among the cox1s of the 25 mtDNAs, whereas gain of a Group II intron (2.65 kb) and loss of a Group I intron (1.7 kb) in cox1 of MF774813.1 resulted in the longest cox1 (12 kb). In rnl of L. edodes, we discovered four intron insertion consensus sequences which were unique to basidiomycetes but not ascomycetes. Differential incorporation of introns was the primary cause of the rnl size polymorphism. Homing endonucleases (HEGs) were suggestively involved in the mobilization of the introns because all of the introns have HEG genes of the LAGRIDADG or GIY-YIG families with the conserved HEG cleavage sites. TEs contributed to 11.04% of the mtDNA size in average, of which 7.08% was LTR-retrotransposon and 3.96% was DNA transposon, whereas the repeated sequences covered 4.6% of the mtDNA. The repeat numbers were variable in a strain-dependent manner. Both the TEs and repeated sequences were mostly found in the intronic and intergenic regions. Lastly, two major deletions were found in the plasmid-related sequence regions (pol2-pol3 and pol1-atp8) in the five mtDNAs. Particularly, the 6.8 kb-long deletion at pol2-pol3 region made MF774813.1 the shortest mtDNA of all. Our results demonstrate that mtDNA is a dynamic molecule that persistently evolves over a short period of time by insertion/deletion and repetition of DNA segments at the strain level.

Mitochondrial effects on fertility and longevity in Tigriopus californicus contradict predictions of the mother's curse hypothesis.

Strict maternal inheritance of mitochondria favours the evolutionary accumulation of sex-biased fitness effects, as mitochondrial evolution occurs exclusively in female lineages. The 'mother's curse' hypothesis proposes that male-harming mutations should accumulate in mitochondrial genomes when they have neutral or beneficial effects on female fitness. Rigorous empirical tests have largely focused on Drosophila, where support for the predictions of mother's curse has been mixed. We investigated the impact of mother's curse mutations in Tigriopus californicus, a minute crustacean. Using non-recombinant backcrosses, we introgressed four divergent mitochondrial haplotypes into two nuclear backgrounds and recorded measures of fertility and longevity. We found that the phenotypic effects of mitochondrial mutations were context dependent, being influenced by the nuclear background in which they were expressed, as well as the sex of the individual and rearing temperature. Mitochondrial haplotype effects were greater for fertility than longevity, and temperature effects were greater for longevity. However, in opposition to mother's curse expectations, females had higher mitochondrial genetic variance than males for fertility and longevity, little evidence of sexual antagonism favouring females was found, and the impacts of mitonuclear mismatch harmed females but not males. Together, this indicates that selection on mitochondrial variation has not resulted in the accumulation of male mutation load in Tigriopus californicus.

The Complete Mitochondrial Genome of Hyotissa hyotis (Bivalvia: Gryphaeidae) Reveals a Unique Gene Order within Ostreoidea

The mitochondrial (mt) genome is an important tool when studying the evolution of metazoan animals. The oyster family Gryphaeidae, together with Ostreidae, is one of the two extant taxa of superfamily Ostreoidea. Up until now, the available mitochondrial genomes of oysters were all limited to family Ostreidae. In the present study, the first complete mtDNA of family Gryphaeidae represented by Hyotissa hyotis was sequenced and compared with other available ostreoid mtDNA. The mtDNA of H. hyotis is 22,185 bp in length, encoding 13 protein-coding-genes (PCGs), two ribosomal RNA (rRNA) and 23 transfer RNA (tRNA) genes. Within all the intergenic regions that range from 2 to 1528 bp, two large non-coding regions were identified. The first large non-coding region, located between Cox1 and trnA, contains 1528 nucleotides, while the second one is 1191 bp in length and positioned between Cytb and Nad2. The nucleotide composition of the whole mtDNA is A + T biased, accounting for 59.2%, with a negative AT skew value of −0.20 and a positive GC skew value of 0.33. In contrast to the mtDNA of Ostreidae, neither the split of rrnL nor rrnS was detected in that of H. hyotis. The duplication of trnW of H. hyotis was also discovered for the first time within Ostreoidea. The gene order of H. hyotis is quite different from those of ostreids, indicating extensive rearrangements within superfamily Ostreoidea. The reconstructed phylogeny supported H. hyotis as sister to Ostreidae, with the latter clade formed by Ostrea + (Saccostrea + Crassostrea). This study could provide important information for further understanding the mitochondrial evolution of oysters.

Complete mitochondrial genome of Phoneutria depilata (Araneae, Ctenidae): New insights into the phylogeny and evolution of spiders

Spiders (Araneae) are the most abundant terrestrial predators and megadiverse on earth. In recent years, the mitochondrial genome of a great diversity of species has been sequenced, mainly for ecological and commercial purposes. These studies have uncovered the existence of a variety of mitochondrial genome rearrangements. However, there is poor genetic information in several taxonomic families of spiders. We have sequenced the complete genome of Phoneutria depilata (Ctenidae) and, based on this, extract the mitogenomes of other ctenid species from published transcriptomes to perform a comparative study among spider species to determine the relationship between the level of mitochondrial rearrangements and its possible relationship with molecular variability in spiders. Complete mitochondrial genomes of eighteen spiders (including eight Ctenidae species) were obtained by two different methodologies (sequencing and transcriptome extraction). Fifty-eight spider mitochondrial genomes were downloaded from the NCBI database for gene order analysis. After verifying the annotation of each mitochondrial gene, a phylogenetic and a gene order analysis from 76 spider mitochondrial genomes were carried out. Our results show a high rate of annotation error in the published spider mitochondrial genomes, which could lead to errors in phylogenetic inference. Moreover, to provide new mitochondrial genomes in spiders by two different methodologies to obtain them, our analysis identifies six different mitochondrial architectures among all spiders. Translocation or tandem duplication random loss (TDRL) events in tRNA genes were identified to explain the evolution of the spider mitochondrial genome. In addition, our findings provide new insights into spider mitochondrial evolution.

A Draft Mitogenome of Plectus Murrayi.

Plectus murrayi is a free-living microbivorous nematode endemic to Antarctic soils. Our draft assembly of its mitogenome was 15,656 bp long, containing 12 protein-coding, eight transfer RNA (tRNA), and two ribosomal RNA (rRNA) genes. Mitophylogenomic analyses extend our understanding of mitochondrial evolution in Nematoda.

Independent evolution of highly variable, fragmented mitogenomes of parasitic lice

The mitochondrial genomes (mitogenomes) of bilaterian animals are highly conserved structures that usually consist of a single circular chromosome. However, several species of parasitic lice (Insecta: Phthiraptera) possess fragmented mitogenomes, where the mitochondrial genes are present on separate, circular chromosomes. Nevertheless, the extent, causes, and consequences of this structural variation remain poorly understood. Here, we combined new and existing data to better understand the evolution of mitogenome fragmentation in major groups of parasitic lice. We found strong evidence that fragmented mitogenomes evolved many times within parasitic lice and that the level of fragmentation is highly variable, including examples of heteroplasmic arrangements. We also found a significant association between mitochondrial fragmentation and signatures of relaxed selection. Mitochondrial fragmentation was also associated with changes to a lower AT%, possibly due to differences in mutation biases. Together, our results provide a significant advance in understanding the process of mitogenome fragmentation and provide an important perspective on mitochondrial evolution in eukaryotes.