Mitochondrial DNA Research Articles

Organellar genome dynamics of exogenous stages of Eimeria tenella

BackgroundCoccidia are a group of intracellular protozoal parasites within the phylum Apicomplexa. Eimeria tenella, one of the species that cause intestinal coccidiosis in poultry, can cause significant mortality and morbidity. Diploid oocysts of Eimeria species are shed in the feces of an infected host and must sporulate to achieve infectivity. This process results in eight haploid infectious units, called sporozoites, held within a single oocyst. Each Eimeria spp. parasite possesses a single apicoplast and a single mitochondrion, both of which carry multiple copies of their respective organellar genomes. Reports of copy numbers of organellar genomes have varied widely.MethodsWe report the application of quantitative polymerase chain reaction (qPCR), supported by next-generation sequencing, for the quantification of the extranuclear genomes relative to the nuclear genome over the course of sporulation and following its completion.ResultsAt 64 elapsed hours, 93.0% of oocysts were fully sporulated; no increase in percent sporulation was observed after this time. Apicoplast relative genome copy number showed several significant shifts up to 72 elapsed hours, after which no significant shifts were observed. Oocysts were shed with approximately 60% the amount of apicoplast DNA present at 72 h, after which point no significant shifts in apicoplast genome relative abundance occurred. Mitogenome relative copy number showed only two significant shifts, from 16 to 24 elapsed hours and from 24 to 32 elapsed hours. Oocysts were shed with approximately 28% the amount of mitochondrial DNA that was present at the time sporulation was deemed morphologically complete, at 64 elapsed hours.ConclusionsThe characterization of the dynamics of genome abundance in exogenous stages sheds new light on the basic biology of Eimeria spp. and supports the use of extranuclear targets for molecular modes of parasite quantification and identification with improved sensitivity and accuracy.Graphical

Seeing the unseen in characterizing RNA editome during rice endosperm development

Rice (Oryza sativa L.) endosperm is essential to provide nutrients for seed germination and determine grain yield. RNA editing, a post-transcriptional modification essential for plant development, unfortunately, is not fully characterized during rice endosperm development. Here, we perform systematic analyses to characterize RNA editome during rice endosperm development. We find that most editing sites are C-to-U CDS-recoding in mitochondria, leading to increased hydrophobic amino acids and changed structures of mitochondrial proteins. Comparative analysis of RNA editome reveals that CDS-recoding sites present higher editing frequencies with lower variabilities and their resultant recoded amino acids tend to exhibit stronger evolutionary conservation across many land plants. Furthermore, we classify mitochondrial genes into three groups, presenting distinct patterns in terms of CDS-recoding events. Besides, we conduct genome-wide screening to detect pentatricopeptide repeat (PPR) proteins and construct PPR-RNA binding profiles, yielding candidate PPR editing factors related to rice endosperm development. Taken together, our findings provide valuable insights for deciphering fundamental mechanisms of rice endosperm development underlying RNA editing machinery.

Comparative analysis of the whole mitochondrial genomes of four species in sect. Chrysantha (Camellia L.), endemic taxa in China

BackgroundThe sect. Chrysantha Chang of plants with yellow flowers of Camellia species as the “Queen of the Tea Family”, most of these species are narrowly distributed endemics of China and are currently listed Grde-II in National Key Protected Wild Plant of China. They are commercially important plants with horticultural medicinal and scientific research value. However, the study of the sect. Chrysantha species genetics are still in its infancy, to date, the mitochondrial genome in sect. Chrysantha has been still unexplored.ResultsIn this study, we provide a comprehensive assembly and annotation of the mitochondrial genomes for four species within the sect. Chrysantha. The results showed that the mitochondrial genomes were composed of closed-loop DNA molecules with sizes ranging from 850,836 bp (C. nitidissima) to 1,098,121 bp (C. tianeensis) with GC content of 45.71–45.78% and contained 48–58 genes, including 28–37 protein-coding genes, 17–20 tRNA genes and 2 rRNA genes. We also examined codon usage, sequence repeats, RNA editing and selective pressure in the four species. Then, we performed a comprehensive comparison of their basic structures, GC contents, codon preferences, repetitive sequences, RNA editing sites, Ka/Ks ratios, haplotypes, and RNA editing sites. The results showed that these plants differ little in gene type and number. C. nitidissima has the greatest variety of genes, while C. tianeensis has the greatest loss of genes. The Ka/Ks values of the atp6 gene in all four plants were greater than 1, indicating positive selection. And the codons ending in A and T were highly used. In addition, the RNA editing sites differed greatly in number, type, location, and efficiency. Twelve, six, five, and twelve horizontal gene transfer (HGT) fragments were found in C. tianeensis, Camellia huana, Camellia liberofilamenta, and C. nitidissima, respectively. The phylogenetic tree clusters the four species of sect. Chrysantha plants into one group, and C. huana and C. liberofilamenta have closer affinities.ConclusionsIn this study, the mitochondrial genomes of four sect. Chrysantha plants were assembled and annotated, and these results contribute to the development of new genetic markers, DNA barcode databases, genetic improvement and breeding, and provide important references for scientific research, population genetics, and kinship identification of sect. Chrysantha plants.

The complete mitochondrial genome of the endangered Atlantic Pygmy Devil Ray, Mobula hypostoma (Bancroft, 1831), from Brazil.

Mobulidae is a monophyletic family within the Myliobatiformes that comprises pelagic species represented by manta and devil rays. Among the genus Mobula, the Atlantic Pygmy Devil Ray - Mobula hypostoma - is reported in coastal regions exclusively in tropical and subtropical Atlantic Ocean from 1 to 100m deep. In Brazil, M. hypostoma is one of the least studied Mobula species. It is regularly misidentified, especially as Mobula thurstoni, and is commonly listed as bycatch, in fishery inventories, or related to opportunistic sightings in the national territory. Here, we describe the complete nucleotide sequence of the mitochondrial genome (mitogenome) from Mobula hypostoma, which is 18,141bp in length and comprises 13 protein-coding, two ribosomal RNA, and 22 transfer RNA genes. The M. hypostoma mitochondrial genes organisation and mitochondrial genome length are similar to other Mobula species, and the phylogenetic reconstruction indicates M. hypostoma as closely related to Mobula munkiana. The Brazilian mitogenome of M. hypostoma is expected to be a valuable resource for molecular-based species identification, and evolutionary and phylogeography studies.

The genome sequence of the Sprawler moth, Asteroscopus sphinx Hufnagel, 1766

We present a genome assembly from an individual male Asteroscopus sphinx (the Sprawler moth; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence has a total length of 857.30 megabases. Most of the assembly is scaffolded into 32 chromosomal pseudomolecules, including the Z sex chromosome and a putative B chromosome. The mitochondrial genome has also been assembled and is 15.35 kilobases in length.

The genome sequence of the Currant Clearwing moth, Synanthedon tipuliformis (Clerck, 1759)

We present a genome assembly from an individual male Synanthedon tipuliformis (the Currant Clearwing; Arthropoda; Insecta; Lepidoptera; Sesiidae). The genome sequence is 295.8 megabases in span. Most of the assembly (99.98%) is scaffolded into 31 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 27.05 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,878 protein-coding genes.

The genome sequence of a tachinid fly, Germaria angustata (Zetterstedt, 1844)

We present a genome assembly from an individual female tachinid fly Germaria angustata (Arthropoda; Insecta; Diptera; Tachinidae). The genome sequence has a total length of 586.50 megabases. Most of the assembly is scaffolded into 6 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 16.77 kilobases in length.

The genome sequence of the Orange-tipped sea squirt, Corella eumyota Traustedt, 1882

We present a genome assembly from an individual specimen of Corella eumyota (the Orange-tipped sea squirt; Chordata; Ascidiacea; Phlebobranchia; Corellidae). The genome sequence is 129.3 megabases in span. Most of the assembly is scaffolded into 7 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 14.53 kilobases in length.

The genome sequence of the Satin Lutestring moth, Tetheella fluctuosa (Hübner, 1803)

We present a genome assembly from an individual female Tetheella fluctuosa (the Satin Lutestring moth; Arthropoda; Insecta; Lepidoptera; Drepanidae). The genome sequence has a total length of 369.10 megabases. Most of the assembly is scaffolded into 32 chromosomal pseudomolecules, including the Z and W sex chromosomes. The mitochondrial genome has also been assembled and is 15.41 kilobases in length. Gene annotation of this assembly on Ensembl identified 18,318 protein-coding genes

The genome sequence of the common green Tenthredo, Tenthredo mesomela Linnaeus, 1758

We present a genome assembly from an individual female Tenthredo mesomela (the common green Tenthredo; Arthropoda; Insecta; Hymenoptera; Tenthredinidae). The genome sequence is 392.8 megabases in span. Most of the assembly is scaffolded into 10 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 15.6 kilobases in length. Gene annotation of this assembly on Ensembl has identified 11,086 protein coding genes.

Communicator whistles: A Trek through the taxonomy of the Boophis marojezensis complex reveals seven new, morphologically cryptic treefrogs from Madagascar (Amphibia: Anura: Mantellidae)

Abstract The Malagasy stream-breeding treefrog species Boophis marojezensis contains bioacoustically and genetically highly divergent populations. Some of these populations have been defined as candidate species and emit somewhat bizarre advertisement calls consisting of multiple whistle-notes. We here enable a long-overdue taxonomic revision of this species complex by applying a museomics approach to sequence DNA from the holotype of B. marojezensis. Based on an integrative approach that combines divergence levels in mitochondrial DNA and in three nuclear-encoded genes, morphological data, and bioacoustic comparisons, we conclude that eight different species exist in this complex, seven of which are formally described herein as new. Although morphological differences between species are small and mainly separate small-sized from larger-sized species, conclusive evidence for the new species comes from their sympatric and sometimes syntopic occurrence without haplotype sharing in three nuclear genes and under maintenance of bioacoustic differences. Uncorrected genetic divergences in the mitochondrial 16S rRNA gene are >3% in almost all cases, and in some cases up to 8%. In reference to the otherworldly sounds by which these frogs fill Malagasy rainforests, some of them reminiscent of sounds of technical equipment in the fictional “Star Trek” universe, we here name and describe the seven new species in honor of fictional captains of starships, namely B. kirkisp. nov., B. picardisp. nov., B. siskoisp. nov., B. janewayaesp. nov., B. archerisp. nov., B. pikeisp. nov., and B. burnhamaesp. nov. The majority of these species occur in northern Madagascar, where up to three species can occur in immediate geographical proximity, e.g., B. marojezensis, B. burnhamaesp. nov. and B. pikeisp. nov. at different elevations in the Marojejy Massif. South of 16°S latitude, only B. janewayaesp. nov., B. picardisp. nov., and B. kirkisp. nov. are found, with the latter extending southwards to Ranomafana National Park. Our study confirms the existence of numerous morphologically cryptic and microendemic species among Madagascar’s amphibians, some of which are known only from unprotected sites and require adequate conservation management.

High mitochondrial DNA content is a key determinant of stemness, proliferation, cell migration, and cancer metastasis in vivo

Here, we examined the potential role of mitochondrial DNA (mtDNA) levels in conveying aggressive phenotypes in cancer cells, using two widely-used breast cell lines as model systems (MCF7[ER+] and MDA-MB-231[ER-]). These human breast cancer cell lines were fractionated into mtDNA-high and mtDNA-low cell sub-populations by flow cytometry, using SYBR Gold as a vital probe to stain mitochondrial nucleoids in living cells. Enrichment of mtDNA-high and mtDNA-low cell sub-populations was independently validated, using a specific DNA-binding mAb probe (AC-30-10), and mitochondrial-based functional assays. As predicted, mtDNA-high MCF7 cells showed significant increases in mitochondrial mass, membrane potential, and superoxide production, as well as increased mitochondrial respiration and ATP production. Moreover, mtDNA-high MCF7 cells demonstrated increases in stemness features, such as anchorage-independent growth and CD44 levels, as well as drug-resistance to Gemcitabine and Tamoxifen. Proliferation rates were also significantly increased, with a dramatic shift towards the S- and G2/M-phases of the cell cycle; this was indeed confirmed by RNA-Seq analysis. Complementary results were obtained with MDA-MB-231 cells. More specifically, mtDNA-high MDA-MB-231 cells showed increases in stemness features and ATP production, as well as rapid cell cycle progression. Moreover, mtDNA-high MDA-MB-231 cells also exhibited increases in both cell migration and invasion, suggesting a role for mtDNA in distant metastasis. To test this hypothesis more directly, a preclinical in vivo model was utilized. For this purpose, MDA-MB-231 tumour cell grafts were treated with an established mtDNA synthesis inhibitor, namely Alovudine (3’-deoxy-3’-fluorothymidine). As expected, drug-induced depletion of mtDNA led to a shift from mitochondrial to glycolytic metabolism. Interestingly, Alovudine very effectively reduced the formation of spontaneous metastases by nearly 70%, but minimally inhibited tumour growth by approximately 20%. Taken together, these data suggest that high mtDNA content is a key driver of stemness, proliferation, and migration, as well as cancer cell metastasis.

Diversity, divergence time, and biogeography of the genus Albatrellus (Agaricomycetes, Russulales)

ABSTRACT The genus Albatrellus is an important group of stipitate terrestrial fungi in the order Russulales. Some species in the genus form ectomycorrhizae, mostly with trees of Pinaceae; some are well-known edible mushrooms. However, its diversity and biogeography are unclear. Taxonomic and phylogenetic studies on Albatrellus were carried out by morphological examination, which included detailed observations of the fruiting body, spore shape and size, and other key features, together with potential hosts. These observations were then compared and analysed using multi-locus molecular phylogenetic analyses, including the internal transcribed spacer regions (ITS), the large subunit nuclear ribosomal RNA gene (nLSU), the translation elongation factor 1-α gene (tef1), the largest subunit of RNA polymerase II (rpb1), the second largest subunit of RNA polymerase II (rpb2), the small subunit mitochondrial rRNA gene (mtSSU), and the small subunit of the nuclear ribosomal RNA gene (nucSSU). The results demonstrated that the species of Albatrellus formed eight clades. Nine new species are described and illustrated, and two new combinations are proposed. A total of 38 species are accepted in Albatrellus worldwide. Of those species, 26, 7, and 8 species are distributed in Asia, Europe, and North America, respectively. The divergence time indicated that the maximum crown age of Albatrellus was approximately 70.5 million years ago, and East Asia and North America are the likely ancestral areas. Dispersal and differentiation to other continents occurred during the late Paleocene and Miocene. Three kinds of dispersal routes are proposed: East Asia and Europe, East Asia and North America, and Europe and North America.

Nuclear compensatory evolution driven by mito-nuclear incompatibilities

Mitochondrial function relies on the coordinated expression of mitochondrial and nuclear genes, exhibiting remarkable resilience despite high mitochondrial mutation rates. The nuclear compensation mechanism suggests deleterious mitochondrial alleles drive compensatory nuclear mutations to preserve mito-nuclear compatibility. However, prevalence and factors conditioning this phenomenon remain debated due to its conflicting evidence. Here, we investigate how mito-nuclear incompatibilities impact substitutions in a model for species radiation. Mating success depends on genetic compatibility (nuclear DNA) and spatial proximity. Populations evolve from partially compatible mito-nuclear states, simulating mitochondrial DNA (mtDNA) introgression. Mutations do not confer advantages nor disadvantages, but individual fecundity declines with increasing incompatibilities, selecting for mito-nuclear coordination. We find that selection for mito-nuclear compatibility affects each genome differently based on their initial state. In compatible gene pairs, selection reduces substitutions in both genomes, while in incompatible nuclear genes, it consistently promotes compensation, facilitated by more mismatches. Interestingly, high mitochondrial mutation rates can reduce nuclear compensation by increasing mtDNA rectification, while substitutions in initially compatible nuclear gene are boosted. Finally, the presence of incompatibilities accelerates species radiation, but equilibrium richness is not directly correlated to substitution rates, revealing the complex dynamics triggered by mitochondrial introgression and mito-nuclear coevolution. Our study provides a perspective on nuclear compensation and the role of mito-nuclear incompatibilities in speciation by exploring extreme scenarios and identifying trends that empirical data alone cannot reveal. We emphasize the challenges in detecting these dynamics and propose analyzing specific genomic signatures could shed light on this evolutionary process.

Complex Models of Sequence Evolution Improve Fit, but not Gene Tree Discordance, for Tetrapod Mitogenomes.

Variation in gene tree estimates is widely observed in empirical phylogenomic data and is often assumed to be the result of biological processes. However, a recent study using tetrapod mitochondrial genomes to control for biological sources of variation due to their haploid, uniparentally inherited, and non-recombining nature found that levels of discordance among mitochondrial gene trees were comparable to those found in studies that assume only biological sources of variation. Additionally, they found that several of the models of sequence evolution chosen to infer gene trees were doing an inadequate job fitting the sequence data. These results indicated that significant amounts of gene tree discordance in empirical data may be due to poor fit of sequence evolution models, and that more complex and biologically realistic models may be needed. To test how the fit of sequence evolution models relates to gene tree discordance, we analyzed the same mitochondrial datasets as the previous study using two additional, more complex models of sequence evolution that each includes a different biologically realistic aspect of the evolutionary process: a covarion model to incorporate site-specific rate variation across lineages (heterotachy), and a partitioned model to incorporate variable evolutionary patterns by codon position. Our results show that both additional models fit the data better than the models used in the previous study, with the covarion being consistently and strongly preferred as tree size increases. However, even these more preferred models still inferred highly discordant mitochondrial gene trees, thus deepening the mystery around what we label the "Mito-Phylo Paradox" and leading us to ask whether the observed variation could, in fact, be biological in nature after all.

A new species of Metaurus Stål, 1866 (Hemiptera, Fulgoromorpha, Dictyopharidae), supplemented with mitogenome data from China

A new planthopper species, Metaurus mohanensis Zheng & Chen, sp. nov., is described and illustrated from Yunnan, China. A key to differentiate species within the genus Metaurus is provided. The geographical distribution of Metaurus species and mitochondrial genome data of the newly described species are also included.

Mitochondrial haplogroup R offers protection against obesity in Kuwaiti and Qatari populations

BackgroundThe Kuwaiti and Qatari populations have a high prevalence of obesity, a major risk factor for various metabolic disorders. Previous studies have independently explored mitochondrial DNA (mtDNA) variations and their association with obesity in these populations. This study aims to investigate the role of mtDNA haplogroups and variants in obesity risk among these Gulf populations.MethodsWhole exome sequencing data from 1,112 participants (348 Kuwaitis and 764 Qataris) were analyzed for mtDNA variants. Participants were classified as obese or non-obese based on body mass index (BMI). Association analyses were performed to examine the relationship between mtDNA haplogroups and obesity, adjusting for covariates such as age and sex.ResultsHaplogroup R was found to be protective against obesity, with an odds ratio (OR) of 0.69 (p = 0.045). This association remained significant after adjusting for age and sex (OR = 0.694; 95% CI: 0.482-0.997; p = 0.048). Several mtDNA variants, particularly those involved in mitochondrial energy metabolism, showed nominal associations with obesity, but these did not remain significant after correcting for multiple testing.ConclusionHaplogroup R consistently demonstrates a protective association against obesity in both Kuwaiti and Qatari populations, highlighting its potential as a biomarker for obesity risk in the Gulf region. However, further research with larger sample sizes is needed to validate these findings and clarify the role of mtDNA variants in obesity.

Rapid evolution of mitochondrion-related genes in haplodiploid arthropods

BackgroundMitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes.ResultsWe show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants.ConclusionsOur results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans.

Assembly and comparative analysis of the first complete mitochondrial genome of zicaitai (Brassica rapa var. Purpuraria): insights into its genetic architecture and evolutionary relationships

Assembly and comparative analysis of the first complete mitochondrial genome of zicaitai (Brassica rapa var. Purpuraria): insights into its genetic architecture and evolutionary relationships

Clinical features and genetic analysis of child with Progressive external ophthalmoplegia with mitochondrial DNA deletions, autosomal dominant 6 due to variant of DNA2 gene

To explore the genetic etiology for a child with Progressive external ophthalmoplegia with mitochondrial DNA deletions, autosomal dominant 6 (PEOA6). A child who had attended the Women and Children's Hospital Affiliated to Ningbo University on 7 August, 2023 was selected as the study subject. Clinical data of the child were analyzed retrospectively. The child and her parents were subjected to whole exome sequencing (WES), and candidate variant was verified by Sanger sequencing and bioinformatic analysis. This study was approved by the Women and Children's Hospital Affiliated to Ningbo University (Ethics No. EC2020-048). The child, a 7-year-old female, had presented with limb muscle pain, amyosthenia, significantly increased creatine kinase, congenital diaphragmatic hernia and recurrent respiratory tract infections. WES revealed that the she has harbored a heterozygous c.1590G>C (p.L530F) variant of the DNA2 gene, which was verified to have a de novo origin by Sanger sequencing. Based on the guidelines from the American College of Medical Genetics and Genomics (ACMG), the c.1590G>C was rated as a likely pathogenic variant (PS2+PM2_Supporting+PP3). The c.1590G>C (p.L530F) variant of the DNA2 gene probably underlay the PEOA6 in this child.