Rapid mitochondrial genome sequencing based on Oxford Nanopore Sequencing and a proxy for vertebrate species identification

Molecular information is crucial for species identification when facing challenging morphology‐based specimen identifications. The use of DNA barcodes partially solves this problem, but in some cases when PCR is not an option (i.e., primers are not available, problems in reaction standardization), amplification‐free approaches could be an optimal alternative. Recent advances in DNA sequencing, like the MinION device from Oxford Nanopore Technologies (ONT), allow to obtain genomic data with low laboratory and technical requirements, and at a relatively low cost. In this study, we explore ONT sequencing for molecular species identification from a total DNA sample obtained from a neotropical rodent and we also test the technology for complete mitochondrial genome reconstruction via genome skimming. We were able to obtain “de novo” the complete mitogenome of a specimen from the genus Melanomys (Cricetidae: Sigmodontinae) with average depth coverage of 78X using ONT‐only data and by combining multiple assembly routines. Our pipeline for an automated species identification was able to identify the sample using unassembled sequence data (raw) in a reasonable computing time, which was substantially reduced when a priori information related to the organism identity was known. Our findings suggest ONT sequencing as a suitable candidate to solve species identification problems in metazoan nonmodel organisms and generate complete mtDNA datasets.     

Affordable de novo generation of fish mitogenomes using amplification-free enrichment of mitochondrial DNA and deep sequencing of long fragments

Biomonitoring surveys make use of metabarcoding tools to describe the community composition. These studies match their sequencing results against public genomic databases to identify the species. However, mitochondrial genomic reference data are yet incomplete, only a few genes may be available, or the suitability of existing sequence data is suboptimal for species level resolution. Here, we present a dedicated and cost-effective workflow with no DNA amplification for generating complete fish mitogenomes for the purpose of strengthening fish mitochondrial databases. Two different strategies using long fragment sequencing with Oxford Nanopore technology coupled with mitochondrial DNA enrichment were used. One where the enrichment is achieved by preferential isolation of mitochondria followed by DNA extraction and nuclear DNA depletion (“mitoenrichment”). A second enrichment approach takes advantage of the CRISPR Cas9 targeted scission on previously dephosphorylated DNA (“targeted mitosequencing”). The sequencing results varied between tissue, species, and integrity of the DNA. The mitoenrichment method yielded 0.17%–12.33% of sequences on target and a mean coverage ranging from 74.9 to 805-fold. The targeted mitosequencing experiment from native genomic DNA yielded 1.83%–55% of sequences on target and a 38 to 2123-fold mean coverage. These produced complete mitogenomes of species with homopolymeric regions, tandem repeats, and gene rearrangements. We demonstrate that deep sequencing of long fragments of native fish DNA can be achieved with low computational resources in a cost-effective manner, opening the discovery of mitogenomes of nonmodel or understudied fish taxa to a broad range of laboratories worldwide.     

Sequencing of whole plastid genomes and nuclear ribosomal DNA of Diospyros species (Ebenaceae) endemic to New Caledonia: many species,little divergence

Background and Aims Some plant groups, especially on islands, have been shaped by strong ancestral bottlenecks and rapid, recent radiation of phenotypic characters. Single molecular markers are often not informative enough for phylogenetic reconstruction in such plant groups. Whole plastid genomes and nuclear ribosomal DNA (nrDNA) are viewed by many researchers as sources of information for phylogenetic reconstruction of groups in which expected levels of divergence in standard markers are low. Here we evaluate the usefulness of these data types to resolve phylogenetic relationships among closely related Diospyros species.Methods Twenty-two closely related Diospyros species from New Caledonia were investigated using whole plastid genomes and nrDNA data from low-coverage next-generation sequencing (NGS). Phylogenetic trees were inferred using maximum parsimony, maximum likelihood and Bayesian inference on separate plastid and nrDNA and combined matrices.Key Results The plastid and nrDNA sequences were, singly and together, unable to provide well supported phylogenetic relationships among the closely related New Caledonian Diospyros species. In the nrDNA, a 6-fold greater percentage of parsimony-informative characters compared with plastid DNA was found, but the total number of informative sites was greater for the much larger plastid DNA genomes. Combining the plastid and nuclear data improved resolution. Plastid results showed a trend towards geographical clustering of accessions rather than following taxonomic species.Conclusions In plant groups in which multiple plastid markers are not sufficiently informative, an investigation at the level of the entire plastid genome may also not be sufficient for detailed phylogenetic reconstruction. Sequencing of complete plastid genomes and nrDNA repeats seems to clarify some relationships among the New Caledonian Diospyros species, but the higher percentage of parsimony-informative characters in nrDNA compared with plastid DNA did not help to resolve the phylogenetic tree because the total number of variable sites was much lower than in the entire plastid genome. The geographical clustering of the individuals against a background of overall low sequence divergence could indicate transfer of plastid genomes due to hybridization and introgression following secondary contact.     

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30–35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.     

Decrease of MtDNA copy number affects mitochondrial function and involves in the pathological consequences of ischaemic stroke

The mtDNA copy number can affect the function of mitochondria and play an important role in the development of diseases. However, there are few studies on the mechanism of mtDNA copy number variation and its effects in IS. The specific mechanism of mtDNA copy number variation is still unclear. In this study, mtDNA copy number of 101 IS patients and 101 normal controls were detected by qRT‐PCR, the effect of D‐loop variation on mtDNA copy number of IS patients was explored. Then, a TFAM gene KD‐OE PC12 cell model was constructed to explore the effect of mtDNA copy number variation on mitochondrial function. The results showed that the mtDNA copy number level of the IS group was significantly lower than that of the normal control group (p < 0.05). The relative expression of TFAM gene mRNA in the cells of the OGD/R treatment group was significantly lower than that of the control group (p < 0.05). In addition, after TFAM gene knockdown and over‐expression plasmids were transfected into HEK 293T cells, mtDNA copy number and ATP production level of Sh‐TFAM transfection group was significantly decreased (p < 0.05), while mtDNA copy number and ATP production level of OE‐TFAM transfected group were significantly higher than that of blank control group and OE‐ctrl negative control group (p < 0.01). Our study demonstrated that mitochondrial D‐loop mutation and TFAM gene dysfunction can cause the decrease of mtDNA copy number, thus affecting the mitochondrial metabolism and function of nerve cells, participating in the pathological damage mechanism of IS.     

Influence of cell distribution and diabetes status on the association between mitochondrial DNA copy number and aging phenotypes in the InCHIANTI study

Mitochondrial DNA copy number (mtDNA‐CN) estimated in whole blood is a novel marker of mitochondrial mass and function that can be used in large population‐based studies. Analyses that attempt to relate mtDNA‐CN to specific aging phenotypes may be confounded by differences in the distribution of blood cell types across samples. Also, low or high mtDNA‐CN may have a different meaning given the presence of diseases associated with mitochondrial damage. We evaluated the impact of blood cell type distribution and diabetes status on the association between mtDNA‐CN and aging phenotypes, namely chronologic age, interleukin‐6, hemoglobin, and all‐cause mortality, among 672 participants of the InCHIANTI study. After accounting for white blood cell count, platelet count, and white blood cell proportions in multivariate models, associations of mtDNA‐CN with age and interleukin‐6 were no longer statistically significant. Evaluation of a statistical interaction by diabetes status suggested heterogeneity of effects in the analysis of mortality (P < 0.01). The magnitude and direction of associations between mtDNA‐CN estimated from blood samples and aging phenotypes are influenced by the sample cell type distribution and disease status. Therefore, accounting for these factors may aid understanding of the relevance of mitochondrial DNA copy number to health and aging.     

Decreased mitochondrial DNA copy number in the hippocampus and peripheral blood during opiate addiction is mediated by autophagy and can be salvaged by melatonin

Drug addiction is a chronic brain disease that is a serious social problem and causes enormous financial burden. Because mitochondrial abnormalities have been associated with opiate addiction, we examined the effect of morphine on mtDNA levels in rat and mouse models of addiction and in cultured cells. We found that mtDNA copy number was significantly reduced in the hippocampus and peripheral blood of morphine-addicted rats and mice compared with control animals. Concordantly, decreased mtDNA copy number and elevated mtDNA damage were observed in the peripheral blood from opiate-addicted patients, indicating detrimental effects of drug abuse and stress. In cultured rat pheochromocytoma (PC12) cells and mouse neurons, morphine treatment caused many mitochondrial defects, including a reduction in mtDNA copy number that was mediated by autophagy. Knockdown of the Atg7 gene was able to counteract the loss of mtDNA copy number induced by morphine. The mitochondria-targeted antioxidant melatonin restored mtDNA content and neuronal outgrowth and prevented the increase in autophagy upon morphine treatment. In mice, coadministration of melatonin with morphine ameliorated morphine-induced behavioral sensitization, analgesic tolerance and mtDNA content reduction. During drug withdrawal in opiate-addicted patients and improvement of protracted abstinence syndrome, we observed an increase of serum melatonin level. Taken together, our study indicates that opioid addiction is associated with mtDNA copy number reduction and neurostructural remodeling. These effects appear to be mediated by autophagy and can be salvaged by melatonin.     

Mitochondrial DNA copy number is maintained during spermatogenesis and in the development of male larvae to sustain the doubly uniparental inheritance of mitochondrial DNA system in the blue mussel Mytilus galloprovincialis

Doubly uniparental inheritance (DUI) of mitochondrial (mt) DNA has been reported in the blue mussel Mytilus galloprovincialis. In DUI, males inherit both paternal (M type) and maternal (F type) mtDNA. Here we investigated changes in M type mtDNA copy numbers and mitochondrial mass in testicular cells by real‐time polymerase chain reaction and flow cytometry. The ratios of M type mtDNA copy numbers to nuclear DNA content were not different between haploid (1n), diploid (2n) and tetraploid (4n) spermatogenic cells. The mitochondrial mass decreased gradually during spermatogenesis. These results suggest that mtDNA and mitochondrial mass are maintained during spermatogenesis. We then traced M type mtDNA in larvae after fertilization. M type mtDNA was maintained up to 24 h after fertilization in the male‐biased crosses, but decreased significantly in female‐biased crosses (predicted by Mito Tracker staining pattern). These results are strikingly different from those reported for mammals and fish, where it is well known that the mitochondria and mtDNA are reduced during spermatogenesis and that sperm mitochondria and mtDNA are eliminated soon after fertilization. Thus, the M type mtDNA copy number is maintained during spermatogenesis and in the development of male larvae to sustain the DUI system in the blue mussel.     

Nanopore sequencing identifies a higher frequency and expanded spectrum of mitochondrial DNA deletion mutations in human aging

Mitochondrial DNA (mtDNA) deletion mutations cause many human diseases and are linked to age-induced mitochondrial dysfunction. Mapping the mutation spectrum and quantifying mtDNA deletion mutation frequency is challenging with next-generation sequencing methods. We hypothesized that long-read sequencing of human mtDNA across the lifespan would detect a broader spectrum of mtDNA rearrangements and provide a more accurate measurement of their frequency. We employed nanopore Cas9-targeted sequencing (nCATS) to map and quantitate mtDNA deletion mutations and develop analyses that are fit-for-purpose. We analyzed total DNA from vastus lateralis muscle in 15 males ranging from 20 to 81 years of age and substantia nigra from three 20-year-old and three 79-year-old men. We found that mtDNA deletion mutations detected by nCATS increased exponentially with age and mapped to a wider region of the mitochondrial genome than previously reported. Using simulated data, we observed that large deletions are often reported as chimeric alignments. To address this, we developed two algorithms for deletion identification which yield consistent deletion mapping and identify both previously reported and novel mtDNA deletion breakpoints. The identified mtDNA deletion frequency measured by nCATS correlates strongly with chronological age and predicts the deletion frequency as measured by digital PCR approaches. In substantia nigra, we observed a similar frequency of age-related mtDNA deletions to those observed in muscle samples, but noted a distinct spectrum of deletion breakpoints. NCATS-mtDNA sequencing allows the identification of mtDNA deletions on a single-molecule level, characterizing the strong relationship between mtDNA deletion frequency and chronological aging.     

Establishment of a mitochondrial DNA sequence database for the identification of fish species commercially available in South Africa

The limitations intrinsic to morphology-based identification systems have created an urgent need for reliable genetic methods that enable the unequivocal recognition of fish species, particularly those that are prone to overexploitation and/or market substitution. The aim of this study was to develop a comprehensive reference library of DNA sequence data to allow the explicit identification of 53 commercially available fish species in South Africa, most of which were locally caught marine species. Sequences of approximately 655 base pairs were generated for all species from the cytochrome c oxidase I (COI) gene, the region widely adopted for DNA barcoding. Specimens of the genus Thunnus were examined in further detail, employing additional mitochondrial DNA control region sequencing. Cumulative analysis of the sequences from the COI region revealed mean conspecific, congeneric and confamilial Kimura 2-parameter distances of 0.10%, 4.58% and 15.43%, respectively. The results showed that the vast majority (98%) of fish species examined could be readily differentiated by their COI barcodes, but that supplementary control region sequencing was more useful for the discrimination of three Thunnus species. Additionally, the analysis of COI data raised the prospect that Thyrsites atun (snoek) could constitute a species pair. The present study has established the necessary genetic information to permit the unambiguous identification of 53 commonly marketed fish species in South Africa, the applications of which hold a plethora of benefits relating to ecology research, fisheries management and control of commercial practices.     

Nucleotide sequence of the soybean mitochondrial 18S rRNA gene: evidence for a slow rate of divergence in the plant mitochondrial genome

Summary The nucleotide sequence of the 18S rRNA gene from soybean mitochondria was determined and is presented here in comparison to the 18S rRNA genes from wheat and maize mitochondria. All three genes exhibit remarkable sequence similarity supporting the proposal that there is a slower rate of nucleotide divergence in plant mitochondrial DNA (mtDNA) as compared to the mtDNA of animals. A lower degree of sequence similarity is observed between the dicotyledenous plant soybean and either wheat (84%) or maize (85%) than between the two monocots (96%). A possible secondary structure for the soybean 18S rRNA is presented that is analogous to the proposed structure for the E. coli 16S rRNA.     

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.     

Phylogeography and population structure of the Atlantic and Mediterranean green turtle Chelonia mydas: a mitochondrial DNA control region sequence assessment

Mitochondrial (mt) DNA sequences were analysed to resolve the phylogeography and population genetic structure of Atlantic and Mediterranean populations of green turtles ( Chelonia mydas ). Analysis of sequence variation over 487 base pairs of the control (D-loop) region identified 18 haplotypes among 147 individuals from nine nesting populations. Pairwise comparisons of haplotype frequencies distinguished most nesting colonies, indicating significant genetic differentiation among rookeries and a strong propensity for natal homing behaviour by nesting females. Comparison of control region sequence data to earlier restriction fragment length polymorphism (RFLP) data for the same individuals demonstrates approximately a sixfold higher substitution rate in the 5' end of the control region. The sequence data provide higher resolution both in terms of the number of mtDNA genotype variants and the phylogeographic relationships detected within the Atlantic region, and reveal a gene genealogy that distinguishes two groups of haplotypes corresponding to (i) the western Caribbean and Mediterranean, and (ii) eastern Caribbean, South Atlantic and West Africa. The data suggest that phylogeographic patterns in the Atlantic Ocean may be interpreted in terms of female nest site fidelity and episodic dispersal events. The distribution of mtDNA haplotypes within the region is thus explained by the geological and climatic alternations (glacial and interglacial) over the last million years.     

Doubly uniparental inheritance of mitochondrial DNA in freshwater mussels: History and status of the European species

Doubly uniparental inheritance (DUI) is a mode of inheriting mitochondrial DNA that is distinct from strictly maternal inheritance. It has been described in nine and three families of marine and freshwater mussels, respectively, including the European margaritiferids and unionids. Among the 16 freshwater species of Unionida inhabiting Europe, DUI has been described in 9 species of dioecious mussels and was absent from a single hermaphroditic species and from secondary hermaphroditic specimens. The DUI freshwater mussels include two vastly genetically different mitochondrial genomes: maternal (F genome) and paternal (M genome), which coexist within the same specimen but in different tissues. The F genome is present in all female tissues and somatic male tissues. It is inherited in the typical, maternal, manner. Conversely, the M genome is located primarily in the male gonads and generative cells, and is inherited paternally. Dioecious Unionidae display unique characteristics that have been interrelated for over 200 million years: a high fidelity of the transmission of the F and M genomes in DUI and two paths of spermatogenesis–the typical path that produces sperm cells containing mitochondria with the F genome and the atypical path that produces sperm cells with the M genome. The mitogenomes of freshwater mussels display unique features that are not present in any other animal, that is, an additional, gender-specific gene and an elongated cox2 gene occurring exclusively in the M genome. These features mean that the mitochondria, in addition to their basic function of producing energy, also may take part in determining sex in these dioecious organisms.     

Shotgun mitogenomics across body size classes in a local assemblage of tropical Diptera: Phylogeny,species diversity and mitochondrial abundance spectrum

Mitochondrial genomes can be assembled readily from shotgun‐sequenced DNA mixtures of mass‐trapped arthropods (“mitochondrial metagenomics”), speeding up the taxonomic characterization. Bulk sequencing was conducted on some 800 individuals of Diptera obtained by canopy fogging of a single tree in Borneo dominated by small (<1.5 mm) individuals. Specimens were split into five body size classes for DNA extraction, to equalize read numbers across specimens and to study how body size, a key ecological trait, interacts with species and phylogenetic diversity. Genome assembly produced 304 orthologous mitochondrial contigs presumed to each represent a different species. The small‐bodied fraction was the by far most species‐rich (187 contigs). Identification of contigs was through phylogenetic analysis together with 56 reference mitogenomes, which placed most of the Bornean community into seven clades of small‐bodied species, indicating phylogenetic conservation of body size. Mapping of shotgun reads against the mitogenomes showed wide ranges of read abundances within each size class. Ranked read abundance plots were largely log‐linear, indicating a uniformly filled abundance spectrum, especially for small‐bodied species. Small‐bodied species differed greatly from other size classes in neutral metacommunity parameters, exhibiting greater levels of immigration, besides greater total community size. We suggest that the established uses of mitochondrial metagenomics for analysis of species and phylogenetic diversity can be extended to parameterize recent theories of community ecology and biodiversity, and by focusing on the number mitochondria, rather than individuals, a new theoretical framework for analysis of mitochondrial abundance spectra can be developed that incorporates metabolic activity approximated by the count of mitochondria.     

Somatic mosaicism for copy-neutral loss of heterozygosity and DNA copy number variations in the human genome

Background

Somatic mosaicism denotes the presence of genetically distinct populations of somatic cells in one individual who has developed from a single fertilised oocyte. Mosaicism may result from a mutation that occurs during postzygotic development and is propagated to only a subset of the adult cells. Our aim was to investigate both somatic mosaicism for copy-neutral loss of heterozygosity (cn-LOH) events and DNA copy number variations (CNVs) in fully differentiated tissues.

Results

We studied panels of tissue samples (11–12 tissues per individual) from four autopsy subjects using high-resolution Illumina HumanOmniExpress-12 BeadChips to reveal the presence of possible intra-individual tissue-specific cn-LOH and CNV patterns.We detected five mosaic cn-LOH regions >5 Mb in some tissue samples in three out of four individuals. We also detected three CNVs that affected only a portion of the tissues studied in one out of four individuals. These three somatic CNVs range from 123 to 796 kb and are also found in the general population. An attempt was made to explain the succession of genomic events that led to the observed somatic genetic mosaicism under the assumption that the specific mosaic patterns of CNV and cn-LOH changes reflect their formation during the postzygotic embryonic development of germinal layers and organ systems.

Conclusions

Our results give further support to the idea that somatic mosaicism for CNVs, and also cn-LOHs, is a common phenomenon in phenotypically normal humans. Thus, the examination of only a single tissue might not provide enough information to diagnose potentially deleterious CNVs within an individual. During routine CNV and cn-LOH analysis, DNA derived from a buccal swab can be used in addition to blood DNA to get information about the CNV/cn-LOH content in tissues of both mesodermal and ectodermal origin. Currently, the real frequency and possible phenotypic consequences of both CNVs and cn-LOHs that display somatic mosaicism remain largely unknown. To answer these questions, future studies should involve larger cohorts of individuals and a range of tissues.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1916-3) contains supplementary material, which is available to authorized users.