Complete mtGenome Research Articles

Comparative mitogenome research revealed the phylogenetics and evolution of the superfamily Tenebrionoidea (Coleoptera: Polyphage).

Despite the worldwide distribution and rich diversity of the superfamily Tenebrionoidea, the knowledge of the mitochondrial genomes (mtgenome) characteristics of the superfamily is still very limited, and its phylogenetics and evolution remain unresolved. In the present study, we newly sequenced mtgenomes from 19 species belonging to Tenebrionoidea, and a total of 90 mitochondrial genomes from 16 families of Tenebrionoidea were used for phylogenetic analysis. There exist 37 genes for all 82 species of complete mtgenomes of 16 families investigated, and their characteristics are identical as reported mtgenomes of other Tenebrionoids. The Ka/Ks analysis suggests that all 13 PCGs have undergone a strong purifying selection. The phylogenetic analysis suggests the monophyly of Mordellidae, Meloidae, Oedemeridae, Pyrochroidae, Salpingidae, Scraptiidae, Lagriidae, and Tenebrionidae, and the Mordellidae is close to the Ripiphoridae. The "Tenebrionidae clade" and "Meloidae clade" are monophyletic, and both of them are sister groups. In the "Meloidae clade," Meloidae is close to Anthicidae. In the "Tenebrionidae clade," the family Lagriidae and Tenebrionidae are sister groups. The divergence time analysis suggests that Tenebrionoidea originated in the late Jurassic, Meloidae Mordellidae, Lagriidae, and Tenebrionidae in the Cretaceous, Oedemeridae in Paleogene. The work lays a base for the study of mtgenome, phylogenetics, and evolution of the superfamily Tenebrionoidea.

Rapid and collective determination of the complete "hot-spring frog" mitochondrial genome containing long repeat regions using Nanopore sequencing.

The mitochondrial genome (mt-genome) is one of the promising molecular markers for phylogenetics and population genetics. Recently, various mt-genomes have been determined rapidly by using massively parallel sequencers. However, the control region (CR, also called D-loop) in mt-genomes remain difficult to precisely determine due to the presence of repeat regions. Here, using Nanopore sequencing, we succeeded in rapid and collective determination of complete mt-genome of the hot-spring frog, Buergeria japonica, and found that its mt-genome size was 22,274 bp including CR (6,929 bp) with two types of tandem repeat motifs forming repeat regions. Comparison of assembly strategies revealed that the long- and short-read data combined together enabled efficient determination of the CR, but the short-read data alone did not. The B. japonica CR was longer than that of a congenic species inhabiting cooler climate areas, Buergeria buergeri, because of the long repeat regions in the former. During the thermal adaptation of B. japonica, the longer repeat regions in its CR may have accumulated within a period after divergence from B. buergeri.

Whole Mitochondrial Genome Detection and Analysis of Two- to Four-Generation Maternal Pedigrees Using a New Massively Parallel Sequencing Panel.

Mitochondrial DNA (mtDNA) is an effective genetic marker in forensic practice, especially for aged bones and hair shafts. Detection of the whole mitochondrial genome (mtGenome) using traditional Sanger-type sequencing is laborious and time-consuming. Additionally, its ability to distinguish point heteroplasmy (PHP) and length heteroplasmy (LHP) is limited. The application of massively parallel sequencing in mtDNA detection helps researchers to study the mtGenome in-depth. The ForenSeq mtDNA Whole Genome Kit, which contains a total of 245 short amplicons, is one of the multiplex library preparation kits for the mtGenome. We used this system to detect the mtGenome in the blood samples and hair shafts of thirty-three individuals from eight two-generation pedigrees, one three-generation pedigree, and one four-generation pedigree. High-quality sequencing results were obtained. Ten unique mtGenome haplotypes were observed in the mothers from the ten pedigrees. A total of 26 PHPs were observed using the interpretation threshold of 6%. Eleven types of LHPs in six regions were evaluated in detail. When considering homoplasmic variants only, consistent mtGenome haplotypes were observed between the twice-sequenced libraries and between the blood and hair shafts from the same individual and among maternal relatives in the pedigrees. Four inherited PHPs were observed, and the remainder were de novo/disappearing PHPs in the pedigrees. Our results demonstrate the effective capability of the ForenSeq mtDNA Whole Genome Kit to generate the complete mtGenome in blood and hair shafts, as well as the complexity of mtDNA haplotype comparisons between different types of maternal relatives when heteroplasmy is considered.

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

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

Characterization of the complete mitochondrial genome of Miamiensis avidus causing flatfish scuticociliatosis.

Miamiensis avidus is a parasitic pathogen that causes the disease scuticociliatosis in teleost fish species. It is a ciliate and a free-living marine protozoan belonging to the order Philasterida, subclass Scuticociliatida, class Oligohymenophorea, and phylum Ciliophora. The complete mt-genome of M. avidus was linear and 38,695bp in length with 47 genes, including 40 protein-coding genes, two ribosomal RNA (rRNA) genes, and five transfer RNA (tRNA) genes. Of these, 20 genes typically belong to the clusters of orthologous groups, playing roles in energy production and conversion, translation, ribosomal structure and biogenesis, and defense mechanisms. This is the first report of sequencing and characterization of the mt-genome of M. avidus, which was observed to be linear and possessing the typical ciliate mitochondrial genome organization and phylogenetic relationships. Remarkable differences were observed between M. avidus and other ciliates in the mitochondrially encoded rRNAs, extensive gene loss in ribosomal genes and tRNAs, terminal repeat sequences, and stop codon usage. A comparative and phylogenetic analysis of M. avidus and Uronema marinum of the order Hymenostomatida, which is most closely related to the order Philasterida, signified the promise of the mitogenome data of M. avidus as a valuable genetic marker in species detection and taxonomic research. The present study has potential applications in epidemiological studies and host-parasite interaction investigations facilitating disease control.

Large dataset of octocoral mitochondrial genomes provides new insights into mt-mutS evolution and function

All studied octocoral mitochondrial genomes (mt-genomes) contain a homologue of the Escherichia coli mutS gene, a member of a gene family encoding proteins involved in DNA mismatch repair, other types of DNA repair, meiotic recombination, and other functions. Although mutS homologues are found in all domains of life, as well as viruses, octocoral mt-mutS is the only such gene found in an organellar genome. While the function of mtMutS is not known, its domain architecture, conserved sequence, and presence of several characteristic residues suggest its involvement in mitochondrial DNA repair. This inference is supported by exceptionally low rates of mt-sequence evolution observed in octocorals. Previous studies of mt-mutS have been limited by the small number of octocoral mt-genomes available. We utilized sequence-capture data from the recent Quattrini et al. 2020 study [Nature Ecology & Evolution 4:1531–1538] to assemble complete mt-genomes for 94 species of octocorals. Combined with sequences publicly available in GenBank, this resulted in a dataset of 184 complete mt-genomes, which we used to re-analyze the conservation and evolution of mt-mutS. In our analysis, we discovered the first case of mt-mutS loss among octocorals in one of the two Pseudoanthomastus spp. assembled from Quattrini et al. data. This species displayed accelerated rate and changed patterns of nucleotide substitutions in mt-genome, which we argue provide additional evidence for the role of mtMutS in DNA repair. In addition, we found accelerated mt-sequence evolution in the presence of mt-mutS in several octocoral lineages. This accelerated evolution did not appear to be the result of relaxed selection pressure and did not entail changes in patterns of nucleotide substitutions. Overall, our results support previously reported patterns of conservation in mt-mutS and suggest that mtMutS is involved in DNA repair in octocoral mitochondria. They also indicate that the presence of mt-mutS contributes to, but does not fully explain, the low rates of sequence evolution in octocorals

First report of complete mitochondrial genome in the subfamily Alleculinae and mitochondrial genome-based phylogenetics in Tenebrionidae (Coleoptera: Tenebrionoidea).

Despite worldwide distribution and rich diversity, the knowledge of mitochondrial genome (mtgenome) characteristics within the family Tenebrionidae (Coleoptera) is still very limited, and phylogenetics remain unresolved for the family. In this study, the mtgenomes for 12 species are newly sequenced and annotated. Three of the species represent the first report of complete mtgenomes from the subfamily Alleculinae. Comparative analysis of 36 tenebrionid mtgenomes finds that gene composition and order are the same as a theoretical ancestral insect mtgenome, and AT bias, length variation, and codon usage are consistent with other reported beetle mtgenomes. Some intergenic overlap and gap sequences may contain phylogenetically informative information, whereas neither the conserved nor repeat sequences in the control region do. The subfamilies Lagriinae, Pimeliinae, Stenochiinae, and Alleculinae are found to be monophyletic, and the subfamilies Diaperinae and Tenebrioninae polyphyletic in our study. Furthermore, Lagriinae is sister to the rest of the subfamilies. At the tribal level, the tribes Lagriini in the subfamily Lagrrinae; Strongyliini in Stenochiinae; Cteniopodini in Alleculinae; and Triboliini, Opatrini, and Amarygmini in Tenebrioninae are monophyletic; while Diaperini in Diaperinae, and Tenebrionini in Tenebrioninae are polyphyletic.

Raising names from the dead: A time-calibrated phylogeny of frog shells (Bursidae, Tonnoidea, Gastropoda) using mitogenomic data

With 59 Recent species, Bursidae, known as «frog shells», are a small but widely distributed group of tropical and subtropical gastropods that are most diverse in the Indo-West Pacific. The present study is aimed at reconstructing phylogenetic relationships of bursid gastropods based on extensive and representative taxon sampling. Five genetic markers (cytochrome c oxidase subunit I (cox1), 16 s and 12 s rRNA mitochondrial genes, 28 s rRNA and Histone H3 nuclear gene) were sequenced for over 30 species in every known genus but Crossata. Furthermore, we sequenced the complete mt-genome of 9 species (10 specimens) (Aspa marginata, Marsupina bufo, Korrigania quirihorai, Korrigania fijiensis, Tutufa rubeta, Bursa lamarckii, Lampasopsis rhodostoma (twice), Bufonaria perelegans and Bursa aff. tuberosissima). Our analysis recovered Bursidae as a monophyletic group, whereas the genus Bursa was found to be polyphyletic. The genera Talisman and Dulcerana are resurrected and the genera Alanbeuella gen. nov. and Korrigania gen. nov. are described. Dating analysis using 21 extinct taxa for node and simplified tip calibrations was performed, showing a diversification of the group in two phases. Diversification may be linked to tectonic events leading to biodiversity relocation from the western Tethys toward the Indo-Pacific.

The Complete Mitochondrial Genome of Four Hylicinae (Hemiptera: Cicadellidae): Structural Features and Phylogenetic Implications.

Simple SummaryHylicinae, containing 43 described species in 13 genera of two tribes, is one of the most morphologically unique subfamilies of Cicadellidae. Phylogenetic studies on this subfamily were mainly based on morphological characters or several gene fragments and just involved single or two taxa. No mitochondrial genome was reported in Hylicinae before. Therefore, we sequenced and analyzed four complete mtgenomes of Hylicinae (Nacolus tuberculatus, Hylica paradoxa, Balala fujiana, and Kalasha nativa) for the first time to reveal mtgenome characterizations and reconstruct phylogenetic relationships of this group. The comparative analyses showed the mtgenome characterizations of Hylicinae are similar to members of Membracoidea. In phylogenetic results, Hylicinae was recovered as a monophyletic group in Cicadellidae and formed to the sister group of Coelidiinae + Iassinae. These results provide the comprehensive framework and worthy information toward the future researches of this subfamily.To reveal mtgenome characterizations and reconstruct phylogenetic relationships of Hylicinae, the complete mtgenomes of four hylicine species, including Nacolus tuberculatus, Hylica paradoxa, Balala fujiana, and Kalasha nativa, were sequenced and comparatively analyzed for the first time. We also carried out the richest (11) subfamily sampling of Cicadellidae to date, and reconstructed phylogenetic relationships of Membracoidea among 61 species based on three datasets using maximum likelihood and Bayesian inference analyses. All new sequenced mtgenomes are molecules ranging from 14,918 to 16,221 bp in length and are double stranded, circular in shape. The gene composition and arrangement of these mtgenomes are consistent with members of Membracoidea. Among 13 protein-coding genes, most show typical ATN start codons and TAR (TAA/TAG) or an incomplete stop codon T–, and several genes start by TTG/GTG. Results of the analysis for sliding window, nucleotide diversity, and nonsynonymous substitution/synonymous substitution indicate cox1 is a comparatively slower-evolving gene while atp8 is the fastest gene. In line with previous researches, phylogenetic results indicate that treehopper families are paraphyletic with respect to family Cicadellidae and also support the monophyly of all involved subfamilies including Hylicinae. Relationships among the four hylicine genera were recovered as (Hylica + (Nacolus + (Balala + Kalasha))).

Searching the undetected mtDNA variants in forensic MPS data

The efficiency of MPS in forensic mtDNA analysis has been thoroughly proven, although a reliable and well established data evaluation still remains a critical point. Numerous bioinformatics tools have been developed, but most of them require specific operating systems and high costs, while free open-source programs with user-friendly interfaces are few. In this study, 43 full mtGenomes were sequenced using the Ion Personal Genome Machine™ (PGM™) System and analyzed utilizing the plug-in Variant Caller (TVC) of the Ion Torrent Software Suite and the mtDNA-Server (mDS), a free web-based mitochondrial analysis tool for MPS data. The outcomes of these two different analysis tools were compared to variants noted after manual inspection of the aligned reads performed using Integrative Genomics Viewer (IGV). The comparison highlighted the presence of thirty-nine discordant variant calls, which were resolved by Sanger sequencing that confirmed the presence of all variants, except for 7 deletions. The combined adoption of IGV and Sanger type sequencing confirmatory steps, in addition of TVC and mDS analysis, resulted in a more accurate variants assignment with the detection of 32 additional true polymorphisms, which were noted in the final dataset. Regarding the heteroplasmy issue, out of a total of thirty heteroplasmic variants, twenty-eight were detected by the TVC, while the mDS detected twenty-two.Overall, none of the used bioinformatics tools were the perfect choice and a secondary analysis with an expert's opinion in complete mtGenome MPS data evaluation is still required in forensic genetic analysis.

Characterization of Two Complete Mitochondrial Genomes of Ledrinae (Hemiptera: Cicadellidae) and Phylogenetic Analysis.

Simple SummaryLedrinae is a small subfamily with many unique characteristics and comprises 5 tribes with 39 genera including approximately 300 species. The monophyly of Ledrinae and the phylogenetic relationships among cicadellid subfamilies remain controversial. To provide further insight into the taxonomic status and phylogenetic status of Ledrinae, two additional complete mitochondrial genomes of Ledrinae species (Tituria sagittata and Petalocephala chlorophana) are newly sequenced and comparatively analyzed. The results showed the sequenced genes of Ledrinae retain the putative ancestral order for insects. In this study, phylogenetic analyses based on expanded sampling and gene data from GenBank indicated that Ledrinae appeared as monophyletic with maximum bootstrap support values and maximum Bayesian posterior probabilities. Bayesian inference and maximum likelihood analysis of concatenated alignments of three datasets produced a well-resolved framework of Cicadellidae and valuable data toward future study in this subfamily.Mitochondrial genomes are widely used for investigations into phylogeny, phylogeography, and population genetics. More than 70 mitogenomes have been sequenced for the diverse hemipteran superfamily Membracoidea, but only one partial and two complete mtgenomes mitochondrial genomes have been sequenced for the included subfamily Ledrinae. Here, the complete mitochondrial genomes (mitogenomes) of two additional Ledrinae species are newly sequenced and comparatively analyzed. Results show both mitogenomes are circular, double-stranded molecules, with lengths of 14,927 bp (Tituria sagittata) and 14,918 bp (Petalocephala chlorophana). The gene order of these two newly sequenced Ledrinae is highly conserved and typical of members of Membracoidea. Similar tandem repeats in the control region were discovered in Ledrinae. Among 13 protein-coding genes (PCGs) of reported Ledrinae mitogenomes, analyses of the sliding window, nucleotide diversity, and nonsynonymous substitution (Ka)/synonymous substitution (Ks) indicate atp8 is a comparatively fast-evolving gene, while cox1 is the slowest. Phylogenetic relationships were also reconstructed for the superfamily Membracoidea based on expanded sampling and gene data from GenBank. This study shows that all subfamilies (sensu lato) are recovered as monophyletic. In agreement with previous studies, these results indicate that leafhoppers (Cicadellidae) are paraphyletic with respect to the two recognized families of treehoppers (Aetalionidae and Membracidae). Relationships within Ledrinae were recovered as (Ledra + (Petalocephala + Tituria)).

Evolutionary implications of analyses of complete mitochondrial genomes across order Zoantharia (Cnidaria: Hexacorallia)

Cnidarians are early-diverging metazoans, but evolutionary aspects of some taxa are still poorly understood, as in the order Zoantharia (Anthozoa: Hexacorallia). Zoantharians have been divided into two suborders based on the arrangement of the fifth septae as complete (Macrocnemina) or incomplete (Brachycnemina). Previous molecular phylogenetic analyses have indicated the need for re-evaluation as Macrocnemina has been found to be paraphyletic. Despite many phylogenetic studies, the recovery of complete mitochondrial genomes (mt-genomes) for systematic and evolutionary studies of zoantharians has been limited. The present study represents the first to sequence the complete mt-genomes of members of eight of nine zoantharian families. Although all examined mt-genomes had the same gene order arrangement, there were variations among mt-genomes' sizes, nucleotide substitution rates, and introns. Only two species did not have the cox1 intron, which harbors a gene coding a homing endonuclease of the LAGLIDADG type. Our mitogenomic analyses also showed relatively high nucleotide diversity in mt-DNA regions other than the standard regions traditionally considered for DNA barcoding of this group. Phylogenetic analyses using 13 mt-genome protein-coding genes recovered a fully resolved tree with clear separation between macrocnemic representatives. Ancestral state reconstruction analyses revealed three main transitions in arrangement of the marginal musculature through the evolutionary history of the order. An “early” transition from reticulate mesogleal to a cteniform endodermal arrangement was followed by transitions that occurred in the common ancestor of the Brachycnemina and family Hydrozoanthidae. Our results indicate the need for clarification of higher-level phylogeny and taxonomy of Zoantharia.

The complete mitochondrial genome of a parasite at the animal-fungal boundary

BackgroundSphaerothecum destruens is an obligate intracellular fish parasite which has been identified as a serious threat to freshwater fishes. Taxonomically, S. destruens belongs to the order Dermocystida within the class Ichthyosporea (formerly referred to as Mesomycetozoea), which sits at the animal-fungal boundary. Mitochondrial DNA (mtDNA) sequences can be valuable genetic markers for species detection and are increasingly used in environmental DNA (eDNA) based species detection. Furthermore, mtDNA sequences can be used in epidemiological studies by informing detection, strain identification and geographical spread.MethodsWe amplified the entire mitochondrial (mt) genome of S. destruens in two overlapping long fragments using primers designed based on the cox1, cob and nad5 partial sequences. The mt-genome architecture of S. destruens was then compared to close relatives to gain insights into its evolution.ResultsThe complete mt-genome of Sphaerothecum destruens is 23,939 bp in length and consists of 47 genes including 21 protein-coding genes, 2 rRNA, 22 tRNA and two unidentified open reading frames. The mitochondrial genome of S. destruens is intronless and compact with a few intergenic regions and includes genes that are often missing from animal and fungal mt-genomes, such as, the four ribosomal proteins (small subunit rps13 and 14; large subunit rpl2 and 16), tatC (twin-arginine translocase component C), and ccmC and ccmF (cytochrome c maturation protein ccmC and heme lyase).ConclusionsWe present the first mt-genome of S. destruens which also represents the first mt-genome for the order Dermocystida. The availability of the mt-genome can assist the detection of S. destruens and closely related parasites in eukaryotic diversity surveys using eDNA and assist epidemiological studies by improving molecular detection and tracking the parasite’s spread. Furthermore, as the only representative of the order Dermocystida, its mt-genome can be used in the study of mitochondrial evolution of the unicellular relatives of animals.

Nanopore sequencing: An enrichment-free alternative to mitochondrial DNA sequencing.

Mitochondrial DNA sequence data are often utilized in disease studies, conservation genetics and forensic identification. The current approaches for sequencing the full mtGenome typically require several rounds of PCR enrichment during Sanger or MPS protocols followed by fairly tedious assembly and analysis. Here we describe an efficient approach to sequencing directly from genomic DNA samples without prior enrichment or extensive library preparation steps. A comparison is made between libraries sequenced directly from native DNA and the same samples sequenced from libraries generated with nine overlapping mtDNA amplicons on the Oxford Nanopore MinION™ device. The native and amplicon library preparation methods and alternative base calling strategies were assessed to establish error rates and identify trends of discordance between the two library preparation approaches. For the complete mtGenome, 16 569 nucleotides, an overall error rate of approximately 1.00% was observed. As expected with mtDNA, the majority of error was detected in homopolymeric regions. The use of a modified basecaller that corrects for ambiguous signal in homopolymeric stretches reduced the error rate for both library preparation methods to approximately 0.30%. Our study indicates that direct mtDNA sequencing from native DNA on the MinION™ device provides comparable results to those obtained from common mtDNA sequencing methods and is a reliable alternative to approaches using PCR‐enriched libraries.

Concordance of mitochondrial DNA sequencing methods on bloodstains using Ion PGM™

In this study, the complete mitochondrial genome (mtGenome) of six samples from three forensic cases was sequenced using the Ion Torrent Personal Genome Machine (PGM). The analyzed samples from forensic cases included bloodstains from several materials, such as gauze, Flinder's Technology Associates (FTA) cards and swabs. The age of the samples ranged from two months to twelve years. The complete mtGenomes were amplified using the tiling sequencing strategy which divided the whole mtGenome into 162 amplicons. All amplicons were successfully recovered. A phylogenetic analysis was performed to determine the accuracy of the PGM data, and which were compared to partial Sanger-based sequencing data. The average coverage of the PGM data were above 4000× in all case samples, and 99.86% concordance was observed using both sequencing methods. In conclusion, we demonstrate the ability to recover the complete mtGenome from bloodstains with relatively poor DNA quality by PGM. Moreover, the results are concordant with Sanger sequencing data. This new method has potential use in forensic practice.

Complete Mitochondrial Genome Sequencing of a Burial from a Romano-Christian Cemetery in the Dakhleh Oasis, Egypt: Preliminary Indications.

The curse of ancient Egyptian DNA was lifted by a recent study which sequenced the mitochondrial genomes (mtGenome) of 90 ancient Egyptians from the archaeological site of Abusir el-Meleq. Surprisingly, these ancient inhabitants were more closely related to those from the Near East than to contemporary Egyptians. It has been accepted that the timeless highway of the Nile River seeded Egypt with African genetic influence, well before pre-Dynastic times. Here we report on the successful recovery and analysis of the complete mtGenome from a burial recovered from a remote Romano–Christian cemetery, Kellis 2 (K2). K2 serviced the ancient municipality of Kellis, a village located in the Dakhleh Oasis in the southwest desert in Egypt. The data were obtained by high throughput sequencing (HTS) performed independently at two ancient DNA facilities (Armed Forces DNA Identification Laboratory, Dover, DE, USA and Carl R. Woese Institute for Genomic Biology, University of Illinois Urbana-Champaign, Urbana, IL, USA). These efforts produced concordant haplotypes representing a U1a1a haplogroup lineage. This result indicates that Near Eastern maternal influence previously identified at Abusir el-Meleq was also present further south, in ancient Kellis during the Romano–Christian period.

Entire Mitochondrial DNA Sequencing on Massively Parallel Sequencing for the Korean Population.

Mitochondrial DNA (mtDNA) genome analysis has been a potent tool in forensic practice as well as in the understanding of human phylogeny in the maternal lineage. The traditional mtDNA analysis is focused on the control region, but the introduction of massive parallel sequencing (MPS) has made the typing of the entire mtDNA genome (mtGenome) more accessible for routine analysis. The complete mtDNA information can provide large amounts of novel genetic data for diverse populations as well as improved discrimination power for identification. The genetic diversity of the mtDNA sequence in different ethnic populations has been revealed through MPS analysis, but the Korean population not only has limited MPS data for the entire mtGenome, the existing data is mainly focused on the control region. In this study, the complete mtGenome data for 186 Koreans, obtained using Ion Torrent Personal Genome Machine (PGM) technology and retrieved from rather common mtDNA haplogroups based on the control region sequence, are described. The results showed that 24 haplogroups, determined with hypervariable regions only, branched into 47 subhaplogroups, and point heteroplasmy was more frequent in the coding regions. In addition, sequence variations in the coding regions observed in this study were compared with those presented in other reports on different populations, and there were similar features observed in the sequence variants for the predominant haplogroups among East Asian populations, such as Haplogroup D and macrohaplogroups M9, G, and D. This study is expected to be the trigger for the development of Korean specific mtGenome data followed by numerous future studies.

Concordance and reproducibility of a next generation mtGenome sequencing method for high-quality samples using the Illumina MiSeq

Sanger-type sequencing (STS) of mitochondrial DNA (mtDNA), specifically the control region (CR), is routinely employed in forensics in human identification and missing persons scenarios. Yet next-generation sequencing (NGS) has the potential to overcome some of the major limitations of STS processing, permitting reasonable paths forward for full mitochondrial genome (mtGenome) sequencing, while also offering higher-throughput and higher sensitivity capabilities. To establish the accuracy and reproducibility of NGS for the development of mtDNA data, 90 DNA extracts that were previously used to generate forensic quality full mtGenomes using STS were sequenced using Nextera XT library preparation and the Illumina MiSeq. Using the same amplicon product, replicate library sets were generated and sequenced at different laboratories, and analysis was performed in replicate using the CLC Genomics Workbench. Both sequencing sets resulted in 99.998% of positions with greater than 10X coverage when 96 samples (including controls) were multiplexed. Overall, 99.9996% concordance was observed between the NGS data and the STS data for the full mtGenome. The only “discordant” calls involved low level point heteroplasmies, with the differences resulting from stochastic variation and/or the increased sensitivity of NGS. Higher sensitivity also allowed for the detection of a mixed sample previously not detected with STS. Additionally, variant calls were reproducible between sequencing sets and between software analysis versions with the variant frequency only differing by 0.23% and 0.01%, respectively. Further validation studies and specialized software functionality tailored to forensic practice should facilitate the incorporation of NGS processing into standard casework applications. The data herein comprise the largest, and likely most thoroughly examined, complete mtGenome STS-NGS concordance dataset available.

The complete mitochondrial genome of Anopheles minimus (Diptera: Culicidae) and the phylogenetics of known Anopheles mitogenomes.

Anopheles minimus is an important vector of human malaria in southern China and Southeast Asia. The phylogenetics of mosquitoes has not been well resolved, and the mitochondrial genome (mtgenome) has proven to be an important marker in the study of evolutionary biology. In this study, the complete mtgenome of An. minimus was sequenced for the first time. It is 15395 bp long and encodes 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs) and a non-coding region. The gene organization is consistent with those of known Anopheles mtgenomes. The mtgenome performs a clear bias in nucleotide composition with a positive AT-skew and a negative GC-skew. All 13 PCGs prefer to use the codon UUA (Leu), ATN as initiation codon but cytochrome-oxidase subunit 1 (COI) and ND5, with TCG and GTG, and TAA as termination codon, but COI, COII, COIII and ND4, all with the incomplete T. tRNAs have the typical clover-leaf structure, but tRNA(Ser(AGN)) is consistent with known Anopheles mtgenomes. The control region includes a conserved T-stretch and a (TA)n stretch, and has the highest A+T content at 93.1%. The phylogenetics of An. minimus with 18 other Anopheles species was constructed by maximum likelihood and Bayesian inference, based on concatenated PCG sequences. The subgenera, Cellia and Anopheles, and Nyssorhynchus and Kerteszia have mutually close relationships, respectively. The Punctulatus group and Leucosphyrus group of Neomyzomyia Series, and the Albitarsis group of Albitarsis Series were suggested to be monophyletic. The monophyletic status of the subgenera, Cellia, Anopheles, Nyssorhynchus and Kerteszia need to be further elucidated.

Strategies for complete mitochondrial genome sequencing on Ion Torrent PGM™ platform in forensic sciences

Next generation sequencing (NGS) is a time saving and cost-efficient method to detect the complete mitochondrial genome (mtGenome) compared to Sanger sequencing. In this study we focused on developing strategies for mtGenome sequencing on the Ion Torrent PGM™ platform and NGS data analysis. With our experience, 4, 15 and 30 samples could be loaded onto Ion 314™, Ion 316™ and Ion 318™ chips respectively at a pooling concentration of 26pM, achieving to sufficient average coverage of ≥1500×and well strand balance of 1.05. Data processing software is essential to NGS mega data analysis. The in-house Perl scripts were developed for primary data analysis to screen out uncertain positions and samples from variant call format (VCF) reports and for pedigree study to perform pairwise comparisons. The Integrative Genomic Viewer (IGV) and the NextGENe software were introduced to secondary data analysis. The mthap and EMMA were employed for haplogroup assignment. The dataset was reviewed and approved by the EMPOP as the final version, which showed 2.66% error rate generated from the Torrent Variant Caller (TVC). Across the mtGenome, 4022 variants were found at 725 nucleotide positions, where ratio of transitions to transversions was estimated at 20.89:1 and 22.18% of variants was concentrated at hypervariable segments I and II (HVS-I and HVS-II). Totally, 107 complete mtGenome haplotypes were observed from 107 Northern Chinese Han and assigned to 88 haplogroups. The random match probability (RMP) of complete mtGenome was calculated as 0.009345794, decreasing 26.19% by comparison to that of HVS-I only, and the haplotype diversity (HD) was evaluated as 1, increasing 0.33% by comparison to that of HVS-I only. Principal component analysis (PCA) showed that our population was clustered to East and Southeast Asians. The strategies in this study are suitable for complete mtGenome sequencing on Ion Torrent PGM™ platform and Northern Chinese Han (EMP00670) is the first complete mtGenome dataset contributed to the EMPOP from East Asian.