Unraveling heteroplasmy patterns with NOVOPlasty

Characterization of Complete Mitochondrial Genome of Badri Breed of Bos indicus (Bovidae: Bovinae): Selection Pressure and Comparative Analysis

High-altitude mammals are often subject to specific environmental obstacles, which exert selective pressure on their physiological and morphological traits, hence driving their evolutionary processes. It is anticipated that these circumstances will lead to the adaptive evolution of protein-coding genes (PCGs) in the mitochondrial genome, which play a crucial role in the oxidative phosphorylation system. In this study, we have generated the complete mitochondrial genome of the Badri breed of Bos indicus inhabiting a high-altitude environment to test the signatures of adaptive evolution on PCGs and their phylogenetic relationships. The complete mitogenome of the Badri breed is 16,339 bp and most tRNAs showed typical clover-leaf secondary structure with a few exceptions, like trnS1 and trnS2 without DHU arm and trnK without DHU loop. Comparative analysis of PCGs indicated that cox1 is the most conserved, while atp6 is the most variable gene. Moreover, the ratios of non-synonymous to synonymous substitution rates indicated the purifying selection (Ka/Ks < 1) in the protein-coding genes that shape the diversity in mitogenome of Bos indicus. Furthermore, Branch-site model (BSM) suggested that cox1, cox2, nad3, nad4L, and nad6 underwent stronger purifying selection (ω < 1) than other PCGs in 15 breeds of 4 species, including Badri. BSM also detected 10 positive sites in PCGs and one in 13 PCGs concatenated dataset. Additional analyses in Datamonkey indicated 11 positive sites and 23 purifying sites in the concatenated dataset, a relaxation of selection strength in nad3, and no evidence of episodic diversifying selection in any PCGs. Phylogeny revealed the sister relationship of the Badri with other breeds of Bos indicus as well as Bos frontalis (Gayal-2). The mitogenome of the Badri breed is an important genomic resource for conservation genetics of this species and also contributes to the understanding of the adaptive evolution of mitochondrial protein coding genes.

Mitogenomics providing new insights into the phylogenetic structure of subfamily Panchaetothripinae (Thripidae: Terebrantia)

Complete mitochondrial genome of two species of subfamily Panchaetothripinae, Astrothrips tumiceps (16,467 bp) and Monilothrips kempi (14,773 bp) are generated by Next-Generation Sequencing Method. In this study, the detailed annotation of these mitogenomes as well as comparative analyses are carried out to explore the codon usage, gene composition, and phylogenetic relationship of subfamilies of family Thripidae. Moreover, the gene rearrangement of subfamily Panchaetothripinae of family Thripidae is also studied. Both the mitogenomes featured by 37 genes including 13 PCGs, 22 tRNAs, 2 rRNAs and with single putative control region with a positive AT-skew and negative GC-skew. trnS1 without DHU arm in both species, trnV without DHU arm in M. kempi, and trnE without TΨC loop in As. tumiceps. Further, codon based comparative analysis depicted the existence of natural selection pressure on all the PCGs in all the subfamilies of family Thripidae. The phylogenetic analyses, using the Bayesian inference (BI) and Maximum likelihood (ML) supported the monophyly of two suborders and family Phlaeothripidae. The family Thripidae is recovered as paraphyletic and subfamily Panchaetothripinae is in sister relationship with family Aeolothripidae and Stenurothripidae rather than the other subfamilies of family Thripidae. The gene order of the order Thysanoptera is highly rearranged, while few members of the subfamily Panchaetothripinae showed similar gene order to family Stenurothripidae. Therefore, this study suggests that the phylogenetic relationship between the subfamily Panchaetothripinae and other families is uncertain, necessitating a whole genome-based study to clarify the position of Panchaetothripinae within the suborder Terebrantia.

Characterization and trans-generation dynamics of mitogene pool in the silver carp (Hypophthalmichthys molitrix)

Multicopied mitogenome are prone to mutation during replication often resulting in heteroplasmy. The derived variants in a cell, organ, or an individual animal constitute a mitogene pool. The individual mitogene pool is initiated by a small fraction of the egg mitogene pool. However, the characteristics and relationship between them has not yet been investigated. This study quantitatively analyzed the heteroplasmy landscape, genetic loads, and selection strength of the mitogene pool of egg and hatchling in the silver carp (Hypophthalmichthys molitrix) using high-throughput resequencing. The results showed heteroplasmic sites distribute across the whole mitogenome in both eggs and hatchlings. The dominant substitution was Transversion in eggs and Transition in hatching accounting for 95.23%±2.07% and 85.38%±6.94% of total HP sites, respectively. The total genetic loads were 0.293±0.044 in eggs and 0.228±0.022 in hatchlings (P=0.048). The dN/dS ratio was 58.03±38.98 for eggs and 9.44±3.93 for hatchlings (P=0.037). These results suggest that the mitogenomes were under strong positive selection in eggs with tolerance to variants with deleterious effects, while the selection was positive but much weaker in hatchlings showing marked quality control. Based on these findings, we proposed a trans-generation dynamics model to explain differential development mode of the two mitogene pool between oocyte maturation and ontogenesis of offspring. This study sheds light on significance of mitogene pool for persistence of populations and subsequent integration in ecological studies and conservation practices.

A Mysterious Asian Firefly Genus, Oculogryphus Jeng, Engel & Yang (Coleoptera, Lampyridae): The First Complete Mitochondrial Genome and Its Phylogenetic Implications

The firefly genus Oculogryphus Jeng, Engel & Yang, 2007 is a rare-species group endemic to Asia. Since its establishment, its position has been controversial but never rigorously tested. To address this perplexing issue, we are the first to present the complete mitochondrial sequence of Oculogryphus, using the material of O. chenghoiyanae Yiu & Jeng, 2018 determined through a comprehensive morphological identification. Our analyses demonstrate that its mitogenome exhibits similar characteristics to that of Stenocladius, including a rearranged gene order between trnC and trnW, and a long intergenic spacer (702 bp) between the two rearranged genes, within which six remnants (29 bp) of trnW were identified. Further, we incorporated this sequence into phylogenetic analyses of Lampyridae based on different molecular markers and datasets using ML and BI analyses. The results consistently place Oculogryphus within the same clade as Stenocladius in all topologies, and the gene rearrangement is a synapomorphy for this clade. It suggests that Oculogryphus should be classified together with Stenocladius in the subfamily Ototretinae at the moment. This study provides molecular evidence confirming the close relationship between Oculogryphus and Stenocladius and discovers a new phylogenetic marker helpful in clarifying the monophyly of Ototretinae, which also sheds a new light on firefly evolution.

Complete mitochondrial genome of the worm snail Thylacodes adamsii (Littorinimorpha: Vermetidae) from South Korea

The worm snail Thylacodes adamsii (Mörch, 1859) (Littorinimorpha: Vermetidae) is a sessile gastropod that mainly inhabits rocky shores along the warm temperate to tropical ocean. Herein, the complete mitochondrial genome (mitogenome) of T. adamsii from South Korea was characterized. The genome is 14,913 bp in length and contains 13 protein-coding genes (PCGs), 22 tRNA genes, and 2 rRNA genes. The genome organization and base composition of T. adamsii are similar to those of other vermetids. A phylogenetic tree was reconstructed using maximum likelihood based on the nucleotide sequences of the 13 PCGs; this tree supported the monophyly of Vermetidae. The complete mitogenome of T. adamsii can assist with molecular species identification and vermetid phylogenetic research in the future.

Twelve newly assembled jasmine chloroplast genomes: unveiling genomic diversity, phylogenetic relationships and evolutionary patterns among Oleaceae and Jasminum species

BACKGROUND

Jasmine (Jasminum), renowned for its ornamental value and captivating fragrance, has given rise to numerous species and accessions. However, limited knowledge exists regarding the evolutionary relationships among various Jasminum species.

RESULTS

In the present study, we sequenced seven distinct Jasminum species, resulting in the assembly of twelve high-quality complete chloroplast (cp) genomes. Our findings revealed that the size of the 12 cp genomes ranged from 159 to 165 kb and encoded 134-135 genes, including 86-88 protein-coding genes, 38-40 tRNA genes, and 8 rRNA genes. J. nudiflorum exhibited a larger genome size compared to other species, mainly attributed to the elevated number of forward repeats (FRs). Despite the typically conservative nature of chloroplasts, variations in the presence or absence of accD have been observed within J. sambac. The calculation of nucleotide diversity (Pi) values for 19 cp genomes indicated that potential mutation hotspots were more likely to be located in LSC regions than in other regions, particularly in genes ycf2, rbcL, atpE, ndhK, and ndhC (Pi > 0.2). Ka/Ks values revealed strong selection pressure on the genes rps2, atpA, rpoA, rpoC1, and rpl33 when comparing J. sambac with the three most closely related species (J. auriculatum, J. multiflorum, and J. dichotomum). Additionally, SNP identification, along with the results of Structure, PCA, and phylogenetic tree analyses, divided the Jasminum cp genomes into six groups. Notably, J. polyanthum showed gene flow signals from both the G5 group (J. nudiflorum) and the G3 group (J. tortuosum and J. fluminense). Phylogenetic tree analysis reflected that most species from the same genus clustered together with robust support in Oleaceae, strongly supporting the monophyletic nature of cp genomes within the genus Jasminum.

CONCLUSION

Overall, this study provides comprehensive insights into the genomic composition, variation, and phylogenetic relationships among various Jasminum species. These findings enhance our understanding of the genetic diversity and evolutionary history of Jasminum.

The complete chloroplast genome sequence of Rhamnella wilsonii Schneid (Rhamnaceae)

Rhamnella wilsonii Schneid 1914 is a member of the Rhamnaceae and endemic to China. In this study, the complete chloroplast genome of R. wilsonii was sequenced and assembled. The complete chloroplast genome was 160,049 bp in length, including a pair of inverted repeats (IRs) of 26,502 bp, one large single copy (LSC) region of 88,274 bp and one small single copy (SSC) region of 18,771 bp. The genome contained 129 genes, including 8 rRNA genes, 37 tRNA genes and 84 protein-coding genes. The overall GC content of the complete chloroplast genome was 37.15%. The phylogenetic analysis demonstrated that R. wilsonii is closely related to R. martinni. This study provides basic information for further studies on the identification and evolution of R. wilsonii and Rhamnella from genomic perspective.

Organellar phylogenomics of Ophioglossaceae fern genera

Previous phylogenies showed conflicting relationships among the subfamilies and genera within the fern family Ophioglossaceae. However, their classification remains unsettled where contrasting classifications recognize four to 15 genera. Since these treatments are mostly based on phylogenetic evidence using limited, plastid-only loci, a phylogenomic understanding is actually necessary to provide conclusive insight into the systematics of the genera. In this study, we have therefore compiled datasets with the broadest sampling of Ophioglossaceae genera to date, including all fifteen currently recognized genera, especially for the first time the South African endemic genus Rhizoglossum. Notably, our comprehensive phylogenomic matrix is based on both plastome and mitogenome genes. Inferred from the coding sequences of 83 plastid and 37 mitochondrial genes, a strongly supported topology for these subfamilies is presented, and is established by analyses using different partitioning approaches and substitution models. At the generic level, most relationships are well resolved except for few within the subfamily Ophioglossoideae. With this new phylogenomic scheme, key morphological and genomic changes were further identified along this backbone. In addition, we confirmed numerous horizontally transferred (HGT) genes in the genera Botrypus, Helminthostachys, Mankyua, Sahashia, and Sceptridium. These HGT genes are most likely located in mitogenomes and are predominately donated from angiosperm Santalales or non-Ophioglossaceae ferns. By our in-depth searches of the organellar genomes, we also provided phylogenetic overviews for the plastid and mitochondrial MORFFO genes found in these Ophioglossaceae ferns.

Characterization of complete mitochondrial genome of three Horse flies of the genus Tabanus (Diptera: Tabanidae): comparative analysis

BACKGROUND

Tabanidae (Horse-Flies or Deer-Flies) are one of the most economically important as well as medically significant haematophagous insect family within the order Diptera. Members of this group are also responsible for the mortality of substantial number of live-stock every year. Due to their pathogen transmission potential and vector competencies makes them an important insect group to study. Till now, mitochondrial genome of 18 species of tabanids were available.

METHODS AND RESULTS

The complete mitogenome of three species T. diversifrons (15,809 bp), T. rubidus (15,878 bp) and T. tenens (15,872 bp) were generated by Next generation sequencing method. They consist 37 genes, with a positive AT skew and a negative GC skew. The gene order of these three species is similar to the typical gene arrangement of infra-order Tabanomorpha. Most of the tRNAs showed typical clover-leaf secondary structure except trnS1, which lacks the DHU arm. The sliding window analysis showed that the nad4L is the most conserved while atp8, and nad6 are the most variable genes. Moreover, the ratios of non-synonymous to synonymous substitution rates indicated that all PCGs under the purifying selection. Phylogeny revealed Chrysops and Haematopota are monophyletic while species of Hybomitra are nested within the polyphyletic clade of Tabanus. T. diversifrons exhibits sister relationship with Atylotus miser. Two morphologically divergent species T. rubidus and T. tenens are found to be genetically similar and indistinguishable by mitochondrial genome.

CONCLUSIONS

The hypervariable genes like atp8 and nad6 can be used as molecular markers for the identification of recently diverged lineages of family Tabanidae. Further, to address uncertainties arising from the two morphological divergent species, it is imperative to obtain data from nuclear gene markers.

Extensive mitogenomic heteroplasmy and its implications in the phylogeny of the fish genus Megalobrama

Megalobrama is China's most economically valuable fish genera. Four species make up this genus: M. amblycephala (MA), M. terminalis (MT), M. pellegrini (MP), and M. hoffmanni (MH). Many researchers have investigated the genetic relationship of Megalobrama based on mitochondrial DNA (mtDNA) and discovered that the branches of the phylogenetic tree for MT and MP are intertwined. We hypothesized that this occurs because mitogenomic heteroplasmy is overlooked when working with mtDNA, which causes MP and MT positions to intersect in phylogenetic trees. To eliminate the influence of nuclear mitochondrial DNA fragments (NUMTs) before analyzing mitogenomic heteroplasmy, we used PLastZ to identify NUMTs, which were then removed from the samples for the subsequent heteroplasmy analysis. Using the heteroplasmy caller icHET, we discovered 126, 339, 135, and 203 heteroplasmic variants in six MA, MT, MP, and MH samples. We reconstructed the Megalobrama fish genus's phylogenetic tree using the RY coding method and rejecting the third position on codons, which improved the performance of the phylogenetic tree by increasing the ratio of treeness to relative component variability from 100.02 ± 1.76 to 688.59 ± 190.56. Despite this, the RY coding method cannot alter the intersection of MP and MT positions in phylogenetic trees. We hypothesize that gene flow between MT and MP leads to intertwining mtDNA-based phylogenetic trees. In conclusion, our findings on the mitogenomic heteroplasmy of Megalobrama provide new insights into mtDNA-based phylogenetic studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03523-0.

A proposed metric set for evaluation of genome assembly quality

Quality control is essential for genome assemblies; however, a consensus has yet to be reached on what metrics should be adopted for the evaluation of assembly quality. N50 is widely used for contiguity measurement, but its effectiveness is constantly in question. Prevailing metrics for the completeness evaluation focus on gene space, yet challenging areas such as tandem repeats are commonly overlooked. Achieving correctness has become an indispensable dimension for quality control, while prevailing assembly releases lack scores reflecting this aspect. We propose a metric set with a set of statistic indexes for effective, comprehensive evaluation of assemblies and provide a score of a finished assembly for each metric, which can be utilized as a benchmark for achieving high-quality genome assemblies.

Multichromosomal Mitochondrial Genome of Paphiopedilum micranthum: Compact and Fragmented Genome, and Rampant Intracellular Gene Transfer

Orchidaceae is one of the largest families of angiosperms. Considering the large number of species in this family and its symbiotic relationship with fungi, Orchidaceae provide an ideal model to study the evolution of plant mitogenomes. However, to date, there is only one draft mitochondrial genome of this family available. Here, we present a fully assembled and annotated sequence of the mitochondrial genome (mitogenome) of Paphiopedilum micranthum, a species with high economic and ornamental value. The mitogenome of P. micranthum was 447,368 bp in length and comprised 26 circular subgenomes ranging in size from 5973 bp to 32,281 bp. The genome encoded for 39 mitochondrial-origin, protein-coding genes; 16 tRNAs (three of plastome origin); three rRNAs; and 16 ORFs, while rpl10 and sdh3 were lost from the mitogenome. Moreover, interorganellar DNA transfer was identified in 14 of the 26 chromosomes. These plastid-derived DNA fragments represented 28.32% (46,273 bp) of the P. micranthum plastome, including 12 intact plastome origin genes. Remarkably, the mitogenome of P. micranthum and Gastrodia elata shared 18% (about 81 kb) of their mitochondrial DNA sequences. Additionally, we found a positive correlation between repeat length and recombination frequency. The mitogenome of P. micranthum had more compact and fragmented chromosomes compared to other species with multichromosomal structures. We suggest that repeat-mediated homologous recombination enables the dynamic structure of mitochondrial genomes in Orchidaceae.

A comparison of 25 complete chloroplast genomes between sister mangrove species Kandelia obovata and Kandelia candel geographically separated by the South China Sea

In 2003, Kandelia obovata was identified as a new mangrove species differentiated from Kandelia candel. However, little is known about their chloroplast (cp) genome differences and their possible ecological significance. In this study, 25 whole cp genomes, with seven samples of K. candel from Malaysia, Thailand, and Bangladesh and 18 samples of K. obovata from China, were sequenced for comparison. The cp genomes of both species encoded 128 genes, namely 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes, but the cp genome size of K. obovata was ~2 kb larger than that of K. candle due to the presence of more and longer repeat sequences. Of these, tandem repeats and simple sequence repeats exhibited great differences. Principal component analysis based on indels, and phylogenetic tree analyses constructed with homologous protein genes from the single-copy genes, as well as 38 homologous pair genes among 13 mangrove species, gave strong support to the separation of the two species within the Kandelia genus. Homologous genes ndhD and atpA showed intraspecific consistency and interspecific differences. Molecular dynamics simulations of their corresponding proteins, NAD(P)H dehydrogenase chain 4 (NDH-D) and ATP synthase subunit alpha (ATP-A), predicted them to be significantly different in the functions of photosynthetic electron transport and ATP generation in the two species. These results suggest that the energy requirement was a pivotal factor in their adaptation to differential environments geographically separated by the South China Sea. Our results also provide clues for future research on their physiological and molecular adaptation mechanisms to light and temperature.

A chromosome-length reference genome for the endangered Pacific pocket mouse reveals recent inbreeding in a historically large population

High-quality reference genomes are fundamental tools for understanding population history, and can provide estimates of genetic and demographic parameters relevant to the conservation of biodiversity. The federally endangered Pacific pocket mouse (PPM), which persists in three small, isolated populations in southern California, is a promising model for studying how demographic history shapes genetic diversity, and how diversity in turn may influence extinction risk. To facilitate these studies in PPM, we combined PacBio HiFi long reads with Omni-C and Hi-C data to generate a de novo genome assembly, and annotated the genome using RNAseq. The assembly comprised 28 chromosome-length scaffolds (N50 = 72.6 MB) and the complete mitochondrial genome, and included a long heterochromatic region on chromosome 18 not represented in the previously available short-read assembly. Heterozygosity was highly variable across the genome of the reference individual, with 18% of windows falling in runs of homozygosity (ROH) >1 MB, and nearly 9% in tracts spanning >5 MB. Yet outside of ROH, heterozygosity was relatively high (0.0027), and historical N_e estimates were large. These patterns of genetic variation suggest recent inbreeding in a formerly large population. Currently the most contiguous assembly for a heteromyid rodent, this reference genome provides insight into the past and recent demographic history of the population, and will be a critical tool for management and future studies of outbreeding depression, inbreeding depression, and genetic load.

The complete chloroplast genome sequence of Nymphoides peltata (S.G.Gmel.) Kuntze

Nymphoides peltata is an aquatic floating weed widely distributed in the northern hemisphere of Eurasia. To better determine its phylogenetic relationships with other Nymphoides species and other aquatic plant species, the chloroplast genome of N. peltata was sequenced. The chloroplast genome size is 152,198 bp, consisting of a large single-copy region (84,223 bp) and a small single-copy region (17,817 bp) separated by a pair of inverted repeats with a length of 25,079 bp. The chloroplast genome contains 127 genes, including 85 protein-coding genes, 34 tRNA genes, and eight rRNAs. The maximum-likelihood phylogenetic tree showed that N. peltata is more closely related to other Nymphoides species, providing new insight into the evolution and genetic diversity of aquatic weeds.

Exploring Mitogenomes Diversity of Fusarium musae from Banana Fruits and Human Patients

Fusarium musae has recently been described as a cross-kingdom pathogen causing post-harvest disease in bananas and systemic and superficial infection in humans. The taxonomic identity of fungal cross-kingdom pathogens is essential for confirming the identification of the species on distant infected hosts. Understanding the level of variability within the species is essential to decipher the population homogeneity infecting human and plant hosts. In order to verify that F. musae strains isolated from fruits and patients are part of a common population and to estimate their overall diversity, we assembled, annotated and explored the diversity of the mitogenomes of 18 F. musae strains obtained from banana fruits and human patients. The mitogenomes showed a high level of similarity among strains with different hosts' origins, with sizes ranging from 56,493 to 59,256 bp. All contained 27 tRNA genes and 14 protein-coding genes, rps3 protein, and small and large ribosomal subunits (rns and rnl). Variations in the number of endonucleases were detected. A comparison of mitochondrial endonucleases distribution with a diverse set of Fusarium mitogenomes allowed us to specifically discriminate F. musae from its sister species F. verticillioides and the other Fusarium species. Despite the diversity in F. musae mitochondria, strains from bananas and strains from human patients group together, indirectly confirming F. musae as a cross-kingdom pathogen.

Benchmarking Low-Frequency Variant Calling With Long-Read Data on Mitochondrial DNA

Background: Sequencing quality has improved over the last decade for long-reads, allowing for more accurate detection of somatic low-frequency variants. In this study, we used mixtures of mitochondrial samples with different haplogroups (i.e., a specific set of mitochondrial variants) to investigate the applicability of nanopore sequencing for low-frequency single nucleotide variant detection.

Methods: We investigated the impact of base-calling, alignment/mapping, quality control steps, and variant calling by comparing the results to a previously derived short-read gold standard generated on the Illumina NextSeq. For nanopore sequencing, six mixtures of four different haplotypes were prepared, allowing us to reliably check for expected variants at the predefined 5%, 2%, and 1% mixture levels. We used two different versions of Guppy for base-calling, two aligners (i.e., Minimap2 and Ngmlr), and three variant callers (i.e., Mutserve2, Freebayes, and Nanopanel2) to compare low-frequency variants. We used F₁ score measurements to assess the performance of variant calling.

Results: We observed a mean read length of 11 kb and a mean overall read quality of 15. Ngmlr showed not only higher F₁ scores but also higher allele frequencies (AF) of false-positive calls across the mixtures (mean F₁ score = 0.83; false-positive allele frequencies < 0.17) compared to Minimap2 (mean F₁ score = 0.82; false-positive AF < 0.06). Mutserve2 had the highest F₁ scores (5% level: F₁ score >0.99, 2% level: F₁ score >0.54, and 1% level: F₁ score >0.70) across all callers and mixture levels.

Conclusion: We here present the benchmarking for low-frequency variant calling with nanopore sequencing by identifying current limitations.

The first complete mitochondrial genome of the micro-whip-scorpion Schizomus zhensis (Arachnida: Schizomida) and phylogenetic analysis

The first complete mitochondrial genome of the micro-whip-scorpion (Schizomus zhensis) was assembled from Illumina-based whole genome sequencing. The circular genome is 14,309 bp in length with 13 protein coding genes, 22 tRNA, 2 rRNAs, and a control region. The A + T content of the mitogenome is 76.6%. Maximum likelihood phylogenetic analyses placed Schizomus zhensis as a sister group of (Thelyphonida + Amblypygi), but more mitogenomes of Tetrapulmonata need to be included in further analyses for a conclusive decision about the group systematics.

Benchmarking the Effectiveness and Accuracy of Multiple Mitochondrial DNA Variant Callers: Practical Implications for Clinical Application

Mitochondrial DNA (mtDNA) mutations contribute to human disease across a range of severity, from rare, highly penetrant mutations causal for monogenic disorders to mutations with milder contributions to phenotypes. mtDNA variation can exist in all copies of mtDNA or in a percentage of mtDNA copies and can be detected with levels as low as 1%. The large number of copies of mtDNA and the possibility of multiple alternative alleles at the same DNA nucleotide position make the task of identifying allelic variation in mtDNA very challenging. In recent years, specialized variant calling algorithms have been developed that are tailored to identify mtDNA variation from whole-genome sequencing (WGS) data. However, very few studies have systematically evaluated and compared these methods for the detection of both homoplasmy and heteroplasmy. A publicly available synthetic gold standard dataset was used to assess four mtDNA variant callers (Mutserve, mitoCaller, MitoSeek, and MToolBox), and the commonly used Genome Analysis Toolkit “best practices” pipeline, which is included in most current WGS pipelines. We also used WGS data from 126 trios and calculated the percentage of maternally inherited variants as a metric of calling accuracy, especially for homoplasmic variants. We additionally compared multiple pathogenicity prediction resources for mtDNA variants. Although the accuracy of homoplasmic variant detection was high for the majority of the callers with high concordance across callers, we found a very low concordance rate between mtDNA variant callers for heteroplasmic variants ranging from 2.8% to 3.6%, for heteroplasmy thresholds of 5% and 1%. Overall, Mutserve showed the best performance using the synthetic benchmark dataset. The analysis of mtDNA pathogenicity resources also showed low concordance in prediction results. We have shown that while homoplasmic variant calling is consistent between callers, there remains a significant discrepancy in heteroplasmic variant calling. We found that resources like population frequency databases and pathogenicity predictors are now available for variant annotation but still need refinement and improvement. With its peculiarities, the mitochondria require special considerations, and we advocate that caution needs to be taken when analyzing mtDNA data from WGS data.

Complex mitogenomic rearrangements within the Pectinidae (Mollusca: Bivalvia)

BACKGROUND

Scallops (Bivalvia: Pectinidae) present extraordinary variance in both mitochondrial genome size, structure and content, even when compared to the extreme diversity documented within Mollusca and Bivalvia. In pectinids, mitogenome rearrangements involve protein coding and rRNA genes along with tRNAs, and different genome organization patterns can be observed even at the level of Tribes. Existing pectinid phylogenies fail to resolve some relationships in the family, Chlamydinae being an especially problematic group.

RESULTS

In our study, we sequenced, annotated and characterized the mitochondrial genome of a member of Chlamydinae, Mimachlamys varia-a species of commercial interest and an effective bioindicator-revealing yet another novel gene arrangement in the Pectinidae. The phylogeny based on all mitochondrial protein coding and rRNA genes suggests the paraphyly of the Mimachlamys genus, further commending the taxonomic revision of the classification within the Chlamydinae subfamily. At the scale of the Pectinidae, we found that 15 sequence blocks are involved in mitogenome rearrangements, which behave as separate units.

CONCLUSIONS

Our study reveals incongruities between phylogenies based on mitochondrial protein-coding versus rRNA genes within the Pectinidae, suggesting that locus sampling affects phylogenetic inference at the scale of the family. We also conclude that the available taxon sampling does not allow for understanding of the mechanisms responsible for the high variability of mitogenome architecture observed in the Pectinidae, and that unraveling these processes will require denser taxon sampling.

Telomere to Telomere Genome Assembly of Fusarium musae F31, Causal Agent of Crown Rot Disease of Banana

Fusarium musae causes crown rot of banana and it is also associated to clinical fusariosis. A chromosome-level genome assembly of F. musae F31 obtained combining Nanopore long reads and Illumina paired-end reads resulted in 12 chromosomes plus one contig with overall N₅₀ of 4.36 Mb, and is presented together with its mitochondrial genome (58,072 bp). The F31 genome includes telomeric regions for 11 of the 12 chromosomes representing one of the most complete genomes available in the Fusarium fujikuroi species complex. The high-quality assembly of the F31 genome will be a valuable resource for studying the pathogenic interactions occurring between F. musae and banana. Moreover, it represents an important resource for understanding the genome evolution in the F. fujikuroi species complex.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Next-generation sequencing reveals the mitogenomic heteroplasmy in the topmouth culter (Culter alburnus Basilewsky, 1855)

BACKGROUND

The mitogenomic heteroplasmy is the presence of multiple haplotypes in the mitochondria, which could cause genetic diseases and is also associated with many critical biological functions. The topmouth culter (Culter alburnus Basilewsky, 1855) is one of the most important freshwater fish in the family of Cyprinidae in China. At present, there are no reports on the topmouth culter's mtDNA heteroplasmy and the existence of which is not known.

METHODS AND RESULTS

This study aimed to analyze the mitogenomic heteroplasmy in the topmouth culter by the next-generation sequencing of the fins' total DNA. The results confirmed the existence of the heteroplasmy and indicated the presence of the extensive heteroplasmy in the topmouth culter's mitogenome. There were 38 heteroplasmic variations in the protein-coding genes from the three specimens, with 33 non-synonymous substitutions accounting for 86.84% and five synonymous substitutions accounting for 13.16%. Among them, the ND6 had the most heteroplasmic variations but only one synonymous substitution. After removing the putative nuclear mitochondrial DNA fragments, the ratio of primary haplotype in the three specimens was 43.89%, 74.72%, and 32.76%, respectively. The three specimens contained 21, 7, and 21 haplotypes of the mitogenomes, respectively. Due to the extensive heteroplasmy, we reconstructed the phylogenetic tree of the topmouth culter using the RY-coding method, which improved the performance of the phylogenetic tree to some extent.

CONCLUSIONS

This study reported the mitogenomic heteroplasmy in the topmouth culter and enhanced the knowledge regarding the mitogenomic heteroplasmy in phylogenetic studies. As the topmouth culter is a commercial species, the mitogenomic heteroplasmy is crucial for the fisheries management of the topmouth culter.

The complete chloroplast genome sequence of Prunus glandulosa Thunb. (Rosaceae)

Prunus glandulosa Thunb. is an ornamental species in the genus Prunus (Rosaceae) mainly distributed in eastern China. It is often cultivated in gardens, roadsides, or shrub clusters. It looks like a cherry but resembles to prunes as well. We obtained the complete chloroplast genome of P. glandulosa using next-generation sequencing technology. The chloroplast genome is 158.078 bp in length with typical tetrad structure. It includes two copies of inverted repeats (IRs, 26.385 bp), a large single copy (LSC, 86.269 bp) and a small single copy (SSC, 19.039 bp). The total GC content is 36.7%, including 85 protein-coding genes (PCGs), 36 transfer RNA genes (tRNA), and eight ribosomal RNA genes (rRNA). The maximum-likelihood phylogeny using the full length of chloroplast genomes indicates that P. glandulosa is closer to prunes than to cherries.

Comparative mitochondrial genome analysis of the Mongolian redfin, Chanodichthys mongolicus (Xenocyprididae) from China reveals heteroplasmy

This study determined the mitochondrial genome (mitogenome) of Chanodichthys mongolicus from China's Qiantang River and analyzed its phylogenetic history in the Subfamily Cultrinae. Next-generation sequencing was used to obtain the mitogenome of C. mongolicus, GenBank Accession Number MZ032228. The mitochondrial genome length of C. mongolicus from China's Qiantang River is 16,622 bp. The genome contains 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and two central noncoding regions (the control region and the origin of light strand replication). Based on BLAST comparisons, the sequence identity of C. mongolicus MZ032228 from China's Qiantang River was 99.84% to that of Ancherythroculter wangi MG783573 from China's Nei River, 99.75% to C. mongolicus AP009060 from Russia's Black River. The phylogenetic analysis is consistent with BLAST comparisons in confirming that A. wangi MG783573 and C. mongolicus MZ032228 show a high genetic similarity. This study also confirms mitochondrial DNA heteroplasmy in C. mongolicus for the first time and documents 35 heterogeneous loci that were detected.

Complete Chloroplast Genomes of Fagus sylvatica L. Reveal Sequence Conservation in the Inverted Repeat and the Presence of Allelic Variation in NUPTs

Growing amounts of genomic data and more efficient assembly tools advance organelle genomics at an unprecedented scale. Genomic resources are increasingly used for phylogenetic analyses of many plant species, but are less frequently used to investigate within-species variability and phylogeography. In this study, we investigated genetic diversity of Fagus sylvatica, an important broadleaved tree species of European forests, based on complete chloroplast genomes of 18 individuals sampled widely across the species distribution. Our results confirm the hypothesis of a low cpDNA diversity in European beech. The chloroplast genome size was remarkably stable (158,428 ± 37 bp). The polymorphic markers, 12 microsatellites (SSR), four SNPs and one indel, were found only in the single copy regions, while inverted repeat regions were monomorphic both in terms of length and sequence, suggesting highly efficient suppression of mutation. The within-individual analysis of polymorphisms showed >9k of markers which were proportionally present in gene and non-gene areas. However, an investigation of the frequency of alternate alleles revealed that the source of this diversity originated likely from nuclear-encoded plastome remnants (NUPTs). Phylogeographic and Mantel correlation analysis based on the complete chloroplast genomes exhibited clustering of individuals according to geographic distance in the first distance class, suggesting that the novel markers and in particular the cpSSRs could provide a more detailed picture of beech population structure in Central Europe.

Characterization and Comparative Analysis of Complete Chloroplast Genomes of Three Species From the Genus Astragalus (Leguminosae)

Astragalus is the largest genus in Leguminosae. Several molecular studies have investigated the potential adulterants of the species within this genus; nonetheless, the evolutionary relationships among these species remain unclear. Herein, we sequenced and annotated the complete chloroplast genomes of three Astragalus species—Astragalus adsurgens, Astragalus mongholicus var. dahuricus, and Astragalus melilotoides using next-generation sequencing technology and plastid genome annotator (PGA) tool. All species belonged to the inverted repeat lacking clade (IRLC) and had similar sequences concerning gene contents and characteristics. Abundant simple sequence repeat (SSR) loci were detected, with single-nucleotide repeats accounting for the highest proportion of SSRs, most of which were A/T homopolymers. Using Astragalus membranaceus var. membranaceus as reference, the divergence was evident in most non-coding regions of the complete chloroplast genomes of these species. Seven genes (atpB, psbD, rpoB, rpoC1, trnV, rrn16, and rrn23) showed high nucleotide variability (Pi), and could be used as DNA barcodes for Astragalus sp. cemA and rpl33 were found undergoing positive selection by the section patterns in the coded protein. Phylogenetic analysis showed that Astragalus is a monophyletic group closely related to the genus Oxytropis within the tribe Galegeae. The newly sequenced chloroplast genomes provide insight into the unresolved evolutionary relationships within Astragalus spp. and are expected to contribute to species identification.

mtDNA Heteroplasmy: Origin, Detection, Significance, and Evolutionary Consequences

Mitochondrial DNA (mtDNA) is predominately uniparentally transmitted. This results in organisms with a single type of mtDNA (homoplasmy), but two or more mtDNA haplotypes have been observed in low frequency in several species (heteroplasmy). In this review, we aim to highlight several aspects of heteroplasmy regarding its origin and its significance on mtDNA function and evolution, which has been progressively recognized in the last several years. Heteroplasmic organisms commonly occur through somatic mutations during an individual’s lifetime. They also occur due to leakage of paternal mtDNA, which rarely happens during fertilization. Alternatively, heteroplasmy can be potentially inherited maternally if an egg is already heteroplasmic. Recent advances in sequencing techniques have increased the ability to detect and quantify heteroplasmy and have revealed that mitochondrial DNA copies in the nucleus (NUMTs) can imitate true heteroplasmy. Heteroplasmy can have significant evolutionary consequences on the survival of mtDNA from the accumulation of deleterious mutations and for its coevolution with the nuclear genome. Particularly in humans, heteroplasmy plays an important role in the emergence of mitochondrial diseases and determines the success of the mitochondrial replacement therapy, a recent method that has been developed to cure mitochondrial diseases.

Comparison of whole genome sequencing and targeted sequencing for mitochondrial DNA

Mitochondrial dysfunction has emerged to be associated with a broad spectrum of diseases, and there is an increasing demand for accurate detection of mitochondrial DNA (mtDNA) variants. Whole genome sequencing (WGS) has been the dominant sequencing approach to identify genetic variants in recent decades, but most studies focus on variants on the nuclear genome. Whole genome sequencing is also costly and time consuming. Sequencing specifically targeted for mtDNA is commonly used in the diagnostic settings and has lower costs. However, there is a lack of pairwise comparisons between these two sequencing approaches for calling mtDNA variants on a population basis. In this study, we compared WGS and mtDNA-targeted sequencing (targeted-seq) in analyzing mitochondrial DNA from 1499 participants recruited into the Severe Asthma Research Program (SARP). Our study reveals that targeted-sequencing and WGS have comparable capacity to determine genotypes and to call haplogroups and homoplasmies on mtDNA. However, there exists a large variability in calling heteroplasmies, especially for low-frequency heteroplasmies, which indicates that investigators should be cautious about heteroplasmies acquired from different sequencing methods. Further research is highly desired to improve variant detection methods for mitochondrial DNA.

Mitochondrial DNA intra-individual variation in a bumblebee species: A challenge for evolutionary studies and molecular identification

Mitochondrial DNA (mtDNA) regions have been widely used as molecular markers in evolutionary studies and species identification. However, the presence of heteroplasmy and NUMTs may represent obstacles. Heteroplasmy is a state where an organism has different mitochondrial haplotypes. NUMTs are nuclear pseudogenes originating from mtDNA sequences transferred to nuclear DNA. Evidences of heteroplasmy were already verified in the bumblebee Bombus morio in an earlier study. The present work investigated in more detail the presence of intra-individual haplotypes variation in this species. Heteroplasmy was detected in individuals from all the ten sampled locations, with an average of six heteroplasmic haplotypes per individual. In addition, some of these heteroplasmic haplotypes were shared among individuals from different locations, suggesting the existence of stable heteroplasmy in B. morio. These results demonstrated that heteroplasmy is likely to affect inferences based on mtDNA analysis, especially in phylogenetic, phylogeographic and population genetics studies. In addition, NUMTs were also detected. These sequences showed divergence of 2.7% to 12% in relation to the mitochondrial haplotypes. These levels of divergence could mislead conclusions in evolutionary studies and affect species identification through DNA barcoding.