Complete mtDNA Research Articles

Comprehensive Molecular Characterization of the Mitochondrial Genome of the Takin Lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae).

The takin lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae) causes lethal bronchopneumonia and represents severe threats to captive and wild populations. However, until now there has been very limited information available concerning the molecular epidemiology and evolutionary biology of V. eleguneniensis. Mitochondrial genomes (mtDNAs) can provide resources for investigations in these areas and, therefore, can assist with the surveillance and control of this lungworm. Herein, the complete mtDNA of V. eleguneniensis was sequenced and characterized with Illumina pipeline analyses. This circular genome (13,625 bp) encoded twelve protein-coding genes (PCGs), two rRNAs, and twenty-two tRNAs, with notable levels of AT and GC skews. Comparative genomics revealed a purifying selection among PCGs, with cox1 and nad6 having the lowest and the highest evolutionary rate, respectively. Genome-wide phylogenies showed a close relationship between V. eleguneniensis and Protostrongylus rufescens in Strongylida. Single gene (PCGs or rRNAs)-based phylogenies indicated that cox1 and nad5 genes shared the same family-level topology with that inferred from genomic datasets, suggesting that both genes could be suitable genetic markers for evolutionary and phylogenetic studies of Strongylida species. This was the first mtDNA of any member of the genus Varestrongylus, and its comprehensive molecular characterization represents a new resource for systematic, population genetic and evolutionary biological studies of Varestrongylus lungworms in wildlife.

The complete mitochondrial genome of Numenius minutus (Charadriiformes: Scolopacidae): comparative and phylogenetic analysis

In this study, we used the next-generation sequencing method to obtain mitochondrial DNA (mtDNA) of Numenius minutus Gould 1841 in Scolopacidae, after which we analyzed the structure and phylogeny of the Charadriiformes. The complete mtDNA was 17,047 bp in length, and contained 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes, and one control region (CR). The gene structure and arrangement of the mitochondrial genome of 64 Charadriiformes species were basically the same as most birds. The reconstructed phylogenetic tree demonstrated that Numenius species were sister groups and monophyletic in Scolopacidae.

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

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

Phylogeography of the striped field mouse (Apodemus agrarius Pallas, 1771) in light of new data from central part of Northern Eurasia.

A phylogeographic analysis of A. agrarius based on the complete mtDNA cytochrome b and control region sequences has been performed using data obtained for the first time for the species from large regions of the central part of Northern Eurasia (23 localities of Altai, Western Siberia, and the Urals). The obtained results have demonstrated a complex intraspecific differentiation of A. agrarius, which has manifested not only in the isolation of the isles populations in Southeast Asia (Jeju and Taiwan), but also in the genetic heterogeneity of mainland populations, which has reflected the history of the modern intraspecific genetic diversity formation against the background of changing physiographic conditions of Eurasia in the Quaternary. The divergence of genetic lineages has taken place apparently simultaneously (in mid-Pleistocene) on the territory of the Eastern part of the modern disjunctive range, where all the identified lineages are present today. The demographic history and possible evolutionary scenarios for A. agrarius in the Western part of the range have been considered. TMRC reconstructions have shown that the lifetime of the common ancestor of the lineage that expanded in the Western Palearctic is about 17.7 [95% HPD 13.2–22.5] kyr. This suggests that the transcontinental expansion of A. agrarius is a relatively recent event that has occurred after the LGM.

The complete paternally inherited mitochondrial genomes of three clam species in genus Macridiscus (Bivalvia: Veneridae): A TDRL model of dimer-mitogenome rearrangement of doubly uniparental inheritance

One of the most striking exceptions to strict maternal inheritance of mitochondrial DNA (mtDNA) in the animal kingdom is a system called doubly uniparental inheritance (DUI), which exists in several bivalve species. DUI is characterized by the presence of two distinct sex-associated mitochondrial lineages: one transmitted through eggs (F-type mtDNA) and the other through sperm (M-type mtDNA). Presently, most known species exhibiting DUI belong to the freshwater bivalve order Unionoida. Other groups with species exhibiting DUI include the orders Mytiloida, Veneroida, and Nuculanoida. In Veneroida, the complete M-type mtDNA is available for two species. We report the presence of DUI in three species belonging to genus Macridiscus (Macridiscus melanaegis, Macridiscus multifarious, and Macridiscus semicancellata), in the order Veneroida, further obtaining their complete M-type mitogenomes. The M-type mitogenome sizes for M. melanaegis, M. multifarious, and M. semicancellata were 19,019 bp, 18,694 bp, and 18,726 bp, respectively, and the mean nucleotide difference between M-type and F-type mitogenomes was 21–23%. We compared the M-type and F-type mitogenomes and found that they show roughly the same genome features, except for gene order. In phylogenetic analyses of Veneroida, a “gender-joining” pattern was revealed within Macridiscus, similar to the pattern of “partial” Mytilus complex (except Unionida). This new insight provides novel evidence supporting the theory that Veneroida and Mytiloida have a more similar DUI pattern than Unionida. A large-scale rearrangement between the sex-linked mitogenomes of the three Macridiscus species was reported. From the observed rearrangement patterns, gene rearrangement between the two sex-linked mitogenomes could be explained by the tandem duplication and random loss (TDRL) model of dimer-mitogenome. This is the first report of heterogeneous genomes with two types of large-scale arrangements in the same organism, and may be contribute significantly to the study of mitochondrial recombination mechanisms.

Origins of Polynesian Pigs Revealed by Mitochondrial Whole Genome Ancient DNA

Simple SummaryRetracing the ancient human migration routes in the remote islands of the Pacific relies on robust models of the origins and spread of animals that were commensal to long-distance ocean voyages. Domestic pigs (Sus scrofa) in Polynesia belong to a rare mitochondrial DNA group whose geographic origins are disputed. We report new complete genome ancient DNA that suggests all founding populations of pigs in Polynesia, first settled by people about 2800–700 years ago, can be traced back to northern peninsular Southeast Asia.Domestic pigs (Sus scrofa) were first transported to Polynesia through a series of long-distance voyages ultimately linked to the Neolithic expansion of Austronesian-speaking people out of Asia. The descendants of the founding pigs belong to a rare mtDNA group referred to as the “Pacific Clade” that may have originated in peninsular or island Southeast Asia. We report the first whole genome mtDNA from domestic pigs from any of the remote islands of the Pacific. In this brief report, we describe the close link we discovered between ancient mtDNA from archaeological specimens from across Polynesia and from that of modern pigs in northern peninsular Southeast Asia, specifically southern China’s Yunnan Province. More complete mtDNA coverage in commensal animals is necessary to improve our picture of the settlement of Polynesia (ca. 2800–700 years before the present) and specify the route, or routes, that pigs took from northern peninsular Southeast Asia.

Intraspecies variation of the mitochondrial genome: An evaluation for phylogenetic approaches based on the conventional choices of genes and segments on mitogenome.

Intraspecies nucleotide sequence variation is a key to understanding the evolutionary history of a species, such as the geographic distribution and population structure. To date, numerous phylogenetic and population genetics studies have been conducted based on the sequences of a gene or an intergenic region on the mitochondrial genome (mtDNA), such as cytochrome c oxidase subunits or the D-loop. To evaluate the credibility of the usage of such ‘classic’ markers, we compared the phylogenetic inferences based on the analyses of the partial and entire mtDNA sequences. Importantly, the phylogenetic reconstruction based on the short marker sequences did not necessarily reproduce the tree topologies based on the analyses of the entire mtDNA. In addition, analyses on the datasets of various organisms revealed that the analyses based on the classic markers yielded phylogenetic trees with poor confidence in all tested cases compared to the results based on full-length mtDNA. These results demonstrated that phylogenetic analyses based on complete mtDNA sequences yield more insightful results compared to those based on mitochondrial genes and segments. To ameliorate the shortcomings of the classic markers, we identified a segment of mtDNA that may be used as an ‘approximate marker’ to closely reproduce the phylogenetic inference obtained from the entire mtDNA in the case of mammalian species, which can be utilized to design amplicon-seq-based studies. Our study demonstrates the importance of the choice of mitochondrial markers for phylogenetic analyses and proposes a novel approach to choosing appropriate markers for mammalian mtDNA that reproduces the phylogenetic inferences obtained from full-length mtDNA.

Analysis of the earliest complete mtDNA genome of a Caribbean colonial horse (Equus caballus) from 16th-century Haiti.

Unlike other European domesticates introduced in the Americas after the European invasion, equids (Equidae) were previously in the Western Hemisphere but were extinct by the late Holocene era. The return of equids to the Americas through the introduction of the domestic horse (Equus caballus) is documented in the historical literature but is not explored fully either archaeologically or genetically. Historical documents suggest that the first domestic horses were brought from the Iberian Peninsula to the Caribbean in the late 15th century CE, but archaeological remains of these early introductions are rare. This paper presents the mitochondrial genome of a late 16th century horse from the Spanish colonial site of Puerto Real (northern Haiti). It represents the earliest complete mitogenome of a post-Columbian domestic horse in the Western Hemisphere offering a unique opportunity to clarify the phylogeographic history of this species in the Americas. Our data supports the hypothesis of an Iberian origin for this early translocated individual and clarifies its phylogenetic relationship with modern breeds in the Americas.

Complete human mtDNA genome sequences revealed haplotype features of the Hmong-Mien language family in Vietnam

Vietnam is the homeland of 54 different ethnicities that belong to 5 major language families of the world, including Austroasiatic, Tai-Kadai, Hmong-Mien, Sino-Tibetan and Austronesian. Hmong-Mien, an ethnolinguistic family, presumably stemmed from Southern China and later spread to the Southeast Asia region. In this study, we analyzed the mitochondrial DNA sequences taken from 120 males belonging Hmong-Mien (HM) language family in Vietnam: Dao, Hmong, and Pathen, revealing 352 unique variants. Dao has the most number of polymorphisms (230 unique SNPs occurring 1469 times), followed by Hmong (181 unique SNPs occurring 1367 times) and Pathen (159 SNPs occurring 1243 times). Genetic variations within each population and among Hmong-Mien speakers were further measured by computations of haplotype diversity (H), nucleotide diversity (π) and fixation index (FST). There are nine major haplogroups (A, B, C, D, F, G, M, N9, and R) detected, with F and B making up over half of each population (Hmong: 56.09% (23/41), Pathen: 58.33% (21/36), Dao: 62.79% (27/43)). Haplotype classification was further divided into 30 haplogroups, of which 80% of them were specific to a single minority. Dao remained the most genetically diverse (H=0.957), while Pathen was the most homogeneous population (H=0.900). In terms of genetic distance, Dao and Hmong were more distinguished from each other, while Hmong and Pathen were more related. Complete mtDNA sequences of Viet HM speakers increased the mtDNA depository, improving the understanding of the genetic structure underlying this language family.

Complete mitochondrial genome and phylogenetic analysis of Mystacoleucus lepturus (Teleostei, Cypriniformes, Cyprinidae)

In this study, we aimed to sequence and annotate the complete mtDNA genome sequence of Mystacoleucus lepturus, which were collected from Luosuojiang River in Menglun area, Yunnan Province, China. The mitochondrial genome was 16,592 bp in length, comprising 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and two non-coding regions (origin of L-strand replication and control region). The whole genome contained C (26.5%), A (32.5%), T (25.3%), and G (15.7%), with an obvious A + T bias (57.8%). Based on the concatenated amino acids sequences of 13 PCGs of M. lepturus and other 22 fishes, a phylogenetic tree was reconstructed using the maximum-likelihood method. The result of phylogenetic analysis supported a close relationship between M. lepturus and M. marginatus. The fundamental genetic data of M. lepturus would be useful for conservation and phylogeny.

Whole exome sequencing is an alternative method in the diagnosis of mitochondrial DNA diseases.

BackgroundMitochondrial disease (MD) is genetically a heterogeneous group of disorders with impairment in respiratory chain complexes or pathways associated with the mitochondrial function. Nowadays, it is still a challenge for the genetic screening of MD due to heteroplasmy of mitochondrial genome and the complex model of inheritance. This study was designed to investigate the feasibility of whole exome sequencing (WES)‐based testing as an alternative option for the diagnosis of MD.MethodsA Chinese Han cohort of 48 patients with suspect MD features was tested using nanoWES, which was a self‐designed WES technique that covered the complete mtDNA genome and 21,019 nuclear genes. Fourteen patients were identified with a single genetic variant and three with single deletion in mtDNA.ResultsThe heteroplasmy levels of variants in mitochondrial genome range from 11% to 100%. NanoWES failed to identify multiple deletions in mtDNA compared with long range PCR and massively parallel sequencing (LR‐PCR/MPS). However, our testing showed obvious advantages in identifying variations in nuclear DNA. Based on nanoWES, we identified two patients with nuclear DNA variation. One of them showed Xp22.33‐q28 duplication, which indicated a possibility of Klinefelter syndrome.ConclusionNanoWES yielded a diagnostic rate of 35.4% for MD. With the rapid advances of next generation sequencing technique and decrease in cost, we recommend the usage of nanoWES as a first‐line method in clinical diagnosis.

Genetic characterization of populations in the Marquesas Archipelago in the context of the Austronesian expansion

Our exploration of the genetic constitution of Nuku Hiva (n = 51), Hiva Oa (n = 28) and Tahuata (n = 8) of the Marquesas Archipelago based on the analyses of genome-wide autosomal markers as well as high-resolution genotyping of paternal and maternal lineages provides us with information on the origins and settlement of these islands at the fringe of the Austronesian expansion. One widespread theme that emerges from this study is the genetic uniformity and relative isolation exhibited by the Marquesas and Society populations. This genetic homogeneity within East Polynesia groups is reflected in their limited average heterozygosity, uniformity of constituents in the Structure analyses, reiteration of complete mtDNA sequences, marked separation from Asian and other Oceanic populations in the PC analyses, limited differentiation in the PCAs and large number of IBD segments in common. Both the f3 and the Outgroup f3 results provide indications of intra-East Polynesian gene flow that may have promoted the observed intra-East Polynesia genetic homogeneity while ALDER analyses indicate that East Polynesia experienced two gene flow episodes, one relatively recent from Europe that coincides roughly with the European incursion into the region and an early one that may represent the original settlement of the islands by Austronesians. Median Network analysis based on high-resolution Y-STR loci under C2a-M208 generates a star-like topology with East Polynesian groups (especially from the Society Archipelago) in central stem positions and individuals from the different populations radiating out one mutational step away while several Samoan and outlier individuals occupy peripheral positions. This arrangement of populations is congruent with dispersals of C2a-M208 Y chromosomes from East Polynesia as a migration hub signaling dispersals in various directions. The equivalent ages of the C2a-M208 lineage of the populations in the Network corroborate an east to west flow of the most abundant Polynesian Y chromosome.

Diatom Biodiversity and Speciation Revealed by Comparative Analysis of Mitochondrial Genomes

Diatoms (Bacillariophyta) constitute one of the most diverse and ecologically significant groups of phytoplankton, comprising 100,000–200,000 species in three classes Bacillariophyceae, Mediophyceae, and Coscinodiscophyceae. However, due to the limited resolution of common molecular markers including 18S rDNA, 28S rDNA, ITS, rbcL, and cox1, diatom biodiversity has not been adequately ascertained. Organelle genomes including mitochondrial genomes (mtDNAs) have been proposed to be “super barcodes” for distinguishing diatom species because of their rich genomic content, and the rapid progress of DNA sequencing technologies that has made it possible to construct mtDNAs with increasing throughout and decreasing cost. Here, we constructed complete mtDNAs of 15 diatom species including five Coscinodiscophyceae species (Guinardia delicatula, Guinardia striata, Stephanopyxis turris, Paralia sulcata, and Actinocyclus sp.), four Mediophyceae species (Hemiaulus sinensis, Odontella aurita var. minima, Lithodesmioides sp., and Helicotheca tamesis), and six Bacillariophyceae species (Nitzschia ovalis, Nitzschia sp., Nitzschia traheaformis, Cylindrotheca closterium, Haslea tsukamotoi, and Pleurosigma sp.) to test the practicality of using mtDNAs as super barcodes. We found that mtDNAs have much higher resolution compared to common molecular markers as expected. Comparative analysis of mtDNAs also suggested that mtDNAs are valuable in evolutionary studies by revealing extensive genome rearrangement events with gene duplications, gene losses, and gains and losses of introns. Synteny analyses of mtDNAs uncovered high conservation among species within an order, but extensive rearrangements including translocations and/or inversions between species of different orders within a class. Duplication of cox1 was discovered for the first time in diatoms in Nitzschia traheaformis and Haslea tsukamotoi. Molecular dating analysis revealed that the three diatom classes split 100 Mya and many diatom species appeared since 50 Mya. In conclusion, more diatom mtDNAs representing different orders will play great dividends to explore biodiversity and speciation of diatoms in different ecological regions.

Genetic variation of the mitochondrial DNA control region across plains bison herds in USA and Canada

Once numbering in the tens of millions, bison (Bison bison) have faced a broad range of challenges over the past century, including genetic impacts from the population bottleneck in the 1800s, and significant loss and fragmentation of habitat resulting in little opportunity for population growth on remaining small, geographically isolated reserves. To identify best practices for bison stewardship against this backdrop, managers must understand the genetic composition of existing conservation herds. This study characterized 14 plains bison (Bison bison bison) herds managed by the U.S. Department of Interior and Parks Canada Agency based on complete mtDNA control region sequences. Among 209 bison, we detected 11 major mtDNA control region haplotypes based on nucleotide substitutions and 23 sub-haplotypes where indels are considered. We determined matrilineal relationships between the herds and compared our genetic findings to historic records. The recent common ancestry of modern bison deriving from small, scattered groups combined with gene flow through foundation and translocation events between herds during the last 100 years, is reflected in Fst value (0.21), haplotype (0.48 ± 0.04) and nucleotide (0.004 ± 0.002) diversities, and mean number of pairwise differences (3.38 ± 1.74). Genetic diversity was distributed unevenly among herds, with 21.41% of genetic variation observed between herds. Median joining network, together with trends in the Tajima's D and Fs tests, revealed two patterns in the recent evolution of mtDNA sequences in bison: mutational process has generated diversity with a Hap 1 haplotype epicenter, and missed mtDNA haplotypes exist in the network due to bottleneck, loss through management practices, or incomplete sampling of specimens across conservation herds. This work significantly expands characterization of the genetic diversity among bison conservation herds, providing additional decision support for managers considering restoring gene flow to achieve long-term species viability.

A mitochondrial genome assembly of the opal chimaera, Chimaera opalescens Luchetti, Iglésias et Sellos 2011, using PacBio HiFi long reads

Chondrichthyans (sharks, rays and chimeras) are a fascinating and highly vulnerable group of early branching gnathostomes. However, they remain comparatively poorly sampled from the point of view of molecular resources, with deep water taxa being particularly data deficient. The development of long-read sequencing technologies enables the analysis of phylogenetic relationships through a precise and reliable assembly of complete mtDNA genomes. The sequencing and characterization of the complete mitogenome of the opal chimera Chimera opalescens Luchetti, Iglésias et Sellos 2011, using the long-read technique PacBio HiFi is presented. The entire mitogenome was 23,411 bp long and shows the same overall content, i.e. 13 protein-coding genes, 22 transfer RNA and 2 ribosomal RNA genes, as all other examined Chondrichthyan mitogenomes. Phylogenetic reconstructions using all available Chondrichthyan mitogenomes, including 11 Holocephali (chimeras and ratfishes), places C. opalescens within the Chimaeridae family. Furthermore, the results reinforce previous findings, showing the genus Chimera as paraphyletic and thus highlighting the need to expand molecular approaches in this group of cartilaginous fishes.

TFB2M and POLRMT are essential for mammalian mitochondrial DNA replication

Two classes of replication intermediates have been observed from mitochondrial DNA (mtDNA) in many mammalian tissue and cells with two-dimensional agarose gel electrophoresis. One is assigned to leading-strand synthesis in the absence of synchronous lagging-strand synthesis (strand-asynchronous replication), and the other has properties of coupled leading- and lagging-strand synthesis (strand-coupled replication). While strand-asynchronous replication is primed by long noncoding RNA synthesized from a defined transcription initiation site, little is known about the commencement of strand-coupled replication. To investigate it, we attempted to abolish strand-asynchronous replication in cultured human cybrid cells by knocking out the components of the transcription initiation complexes, mitochondrial transcription factor B2 (TFB2M/mtTFB2) and mitochondrial RNA polymerase (POLRMT/mtRNAP). Unexpectedly, removal of either protein resulted in complete mtDNA loss, demonstrating for the first time that TFB2M and POLRMT are indispensable for the maintenance of human mtDNA. Moreover, a lack of TFB2M could not be compensated for by mitochondrial transcription factor B1 (TFB1M/mtTFB1). These findings indicate that TFB2M and POLRMT are crucial for the priming of not only strand-asynchronous but also strand-coupled replication, providing deeper insights into the molecular basis of mtDNA replication initiation.

The complete mitochondrial genome for Stanley’s slug snake Pareas stanleyi (Serpentes: Pareidae) by next-generation sequencing

The complete mitochondrial DNA (mtDNA) for Pareas stanleyi was determined in this study. The length of complete mtDNA was 18,932 bp, including 13 protein-coding genes (PCGs) (COI-III, ND1-6, ND4L, ATP6, ATP8 and CYTB), 23 tRNA genes, 2 rRNA genes, a L-chain replication-initiating non-coding region (NCR) and 2 control regions. The overall base composition of the sequence is 24.76% of T, 29.20% of C, 30.87% of A, and 15.16% of G, with a total A + T content of 55.63%. The phylogenetic tree showed that P. stanleyi had a close relationship with the other two species (P. boulengeri and P. formosensis) from the genus Pareas.

A globally diverse reference alignment and panel for imputation of mitochondrial DNA variants

BackgroundVariation in mitochondrial DNA (mtDNA) identified by genotyping microarrays or by sequencing only the hypervariable regions of the genome may be insufficient to reliably assign mitochondrial genomes to phylogenetic lineages or haplogroups. This lack of resolution can limit functional and clinical interpretation of a substantial body of existing mtDNA data. To address this limitation, we developed and evaluated a large, curated reference alignment of complete mtDNA sequences as part of a pipeline for imputing missing mtDNA single nucleotide variants (mtSNVs). We call our reference alignment and pipeline MitoImpute.ResultsWe aligned the sequences of 36,960 complete human mitochondrial genomes downloaded from GenBank, filtered and controlled for quality. These sequences were reformatted for use in imputation software, IMPUTE2. We assessed the imputation accuracy of MitoImpute by measuring haplogroup and genotype concordance in data from the 1000 Genomes Project and the Alzheimer’s Disease Neuroimaging Initiative (ADNI). The mean improvement of haplogroup assignment in the 1000 Genomes samples was 42.7% (Matthew’s correlation coefficient = 0.64). In the ADNI cohort, we imputed missing single nucleotide variants.ConclusionThese results show that our reference alignment and panel can be used to impute missing mtSNVs in existing data obtained from using microarrays, thereby broadening the scope of functional and clinical investigation of mtDNA. This improvement may be particularly useful in studies where participants have been recruited over time and mtDNA data obtained using different methods, enabling better integration of early data collected using less accurate methods with more recent sequence data.

Mitochondrial Genetic Heterogeneity in Leber's Hereditary Optic Neuropathy: Original Study with Meta-Analysis.

Leber’s hereditary optic neuropathy (LHON) is a mitochondrial disorder that causes loss of central vision. Three primary variants (m.3460G>A, m.11778G>A, and m.14484T>C) and about 16 secondary variants are responsible for LHON in the majority of the cases. We investigated the complete mitochondrial DNA (mtDNA) sequences of 189 LHON patients and found a total of 54 disease-linked pathogenic variants. The primary variants m.11778G>A and m.14484T>C were accountable for only 14.81% and 2.64% cases, respectively. Patients with these two variants also possessed additional disease-associated variants. Among 156 patients who lacked the three primary variants, 16.02% harboured other LHON-associated variants either alone or in combination with other disease-associated variants. Furthermore, we observed that none of the haplogroups were explicitly associated with LHON. We performed a meta-analysis of m.4216T>C and m.13708G>A and found a significant association of these two variants with the LHON phenotype. Based on this study, we recommend the use of complete mtDNA sequencing to diagnose LHON, as we found disease-associated variants throughout the mitochondrial genome.

Revisiting the Neurospora crassa mitochondrial genome.

The mitochondrial genome of Neurospora crassa has been less studied than its nuclear counterpart, yet it holds great potential for understanding the diversity and evolution of this important fungus. Here we describe a new mitochondrial DNA (mtDNA) complete sequence of a N. crassa wild type strain. The genome with 64839bp revealed 21 protein-coding genes and several hypothetical open reading frames with no significant homology to any described gene. Five large repetitive regions were identified across the genome, including partial or complete genes. The largest repeated region holds a partial nd2 section that was also detected in Neurospora intermedia, suggesting a rearrangement that occurred before the N. crassa speciation. Interestingly, N. crassa has a palindrome adjacent to the partial nd2 repeated region possibly related to the genomic rearrangement, which is absent in N. intermedia. Finally, we compared the sequences of the three available N. crassa complete mtDNAs and found low levels of intraspecific variability. Most differences among strains were due to small indels in noncoding regions. The revisiting of the N. crassa mtDNA forms the basis for future studies on mitochondrial genome organization and variability.