Transition and transversion rate in the evolution of animal mitochondrial DNA

Single-mitochondrion sequencing uncovers distinct mutational patterns and heteroplasmy landscape in mouse astrocytes and neurons

BACKGROUND

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting "normal" mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes, we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes. In this study, we amplified mt-genomes from 1645 single mitochondria isolated from mouse single astrocytes and neurons to (1) determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, (2) assess differences in mtDNA SNVs between neurons and astrocytes, and (3) study co-segregation of variants in the mouse mtDNA.

RESULTS

(1) The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. (2) Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G > A and 9419:C > T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. (3) Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

CONCLUSIONS

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

Single-Mitochondrion Sequencing Uncovers Distinct Mutational Patterns and Heteroplasmy Landscape in Mouse Astrocytes and Neurons

Background

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting 'normal' mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes In this study we amplified mt-genomes from 1,645 single mitochondria (mts) isolated from mouse single astrocytes and neurons to 1. determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, 2. assess differences in mtDNA SNVs between neurons and astrocytes, and 3. Study cosegregation of variants in the mouse mtDNA.

Results

1. The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. 2. Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G>A and 9419:C>T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. 3. Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

Conclusion

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

Mitochondrial DNA sequencing of Kehilan and Hamdani horses from Saudi Arabia

The Arabian horse breed is well known for its purity and played a key role in the genetic improvement of other horses worldwide. The mitochondrial genome plays a vital role in maternal inheritance and it's helpful to evaluate its genetic diversity and conservation. It has higher mutation rates than nuclear DNA in vertebrates and therefore reveals phylogenetic relationships and haplotypes. In this study, the mitochondrial genome mutations in two Saudi horse strains, Kehilan and Hamdani demonstrated various changes in the gene and amino acid levels and included two other Saudi horses (Hadban and Seglawi) from the previous study for phylogenetic comparison. The whole mitochondrial genome sequencing resulted in intra and inter mtDNA variations between the studied horses. Interestingly, the Hamdani horse has nucleotide substitutions similar to those of the Hadban horse, which is reflected in the phylogenetic tree as a significantly close relationship. This type of study provides a better understanding of mitogenome structure and conservation of livestock species genetic data.

The frequency of the known mitochondrial variants associated with drug-induced toxicity in a Korean population

BACKGROUND

Few studies have annotated the whole mitochondrial DNA (mtDNA) genome associated with drug responses in Asian populations. This study aimed to characterize mtDNA genetic profiles, especially the distribution and frequency of well-known genetic biomarkers associated with diseases and drug-induced toxicity in a Korean population.

METHOD

Whole mitochondrial genome was sequenced for 118 Korean subjects by using a next-generation sequencing approach. The bioinformatic pipeline was constructed for variant calling, haplogroup classification and annotation of mitochondrial mutation.

RESULTS

A total of 681 variants was identified among all subjects. The MT-TRNP gene and displacement loop showed the highest numbers of variants (113 and 74 variants, respectively). The m.16189T > C allele, which is known to reduce the mtDNA copy number in human cells was detected in 25.4% of subjects. The variants (m.2706A > G, m.3010A > G, and m.1095T > C), which are associated with drug-induced toxicity, were observed with the frequency of 99.15%, 30.51%, and 0.08%, respectively. The m.2150T > A, a genotype associated with highly disruptive effects on mitochondrial ribosomes, was identified in five subjects. The D and M groups were the most dominant groups with the frequency of 34.74% and 16.1%, respectively.

CONCLUSIONS

Our finding was consistent with Korean Genome Project and well reflected the unique profile of mitochondrial haplogroup distribution. It was the first study to annotate the whole mitochondrial genome with drug-induced toxicity to predict the ADRs event in clinical implementation for Korean subjects. This approach could be extended for further study for validation of the potential ethnic-specific mitochondrial genetic biomarkers in the Korean population.

A new method for long-read sequencing of animal mitochondrial genomes: application to the identification of equine mitochondrial DNA variants

BACKGROUND

Mitochondrial DNA is remarkably polymorphic. This is why animal geneticists survey mitochondrial genomes variations for fundamental and applied purposes. We present here an approach to sequence whole mitochondrial genomes using nanopore long-read sequencing. Our method relies on the selective elimination of nuclear DNA using an exonuclease treatment and on the amplification of circular mitochondrial DNA using a multiple displacement amplification step.

RESULTS

We optimized each preparative step to obtain a 100 million-fold enrichment of horse mitochondrial DNA relative to nuclear DNA. We sequenced these amplified mitochondrial DNA using nanopore sequencing technology and obtained mitochondrial DNA reads that represented up to half of the sequencing output. The sequence reads were 2.3 kb of mean length and provided an even coverage of the mitochondrial genome. Long-reads spanning half or more of the whole mtDNA provided a coverage that varied between 118X and 488X. We evaluated SNPs identified using these long-reads by Sanger sequencing as ground truth and found a precision of 100.0%; a recall of 93.1% and a F1-score of 0.964 using the Twilight horse mtDNA reference. The choice of the mtDNA reference impacted variant calling efficiency with F1-scores varying between 0.947 and 0.964.

CONCLUSIONS

Our method to amplify mtDNA and to sequence it using the nanopore technology is usable for mitochondrial DNA variant analysis. With minor modifications, this approach could easily be applied to other large circular DNA molecules.

Mitochondrial apurinic/apyrimidinic endonuclease 1 enhances mtDNA repair contributing to cell proliferation and mitochondrial integrity in early stages of hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is the leading cause of primary liver cancers. Surveillance of individuals at specific risk of developing HCC, early diagnostic markers, and new therapeutic approaches are essential to obtain a reduction in disease-related mortality. Apurinic/apyrimidinic endonuclease 1 (APE1) expression levels and its cytoplasmic localization have been reported to correlate with a lower degree of differentiation and shorter survival rate. The aim of this study is to fully investigate, for the first time, the role of the mitochondrial form of APE1 in HCC.

METHODS

As a study model, we analyzed samples from a cohort of patients diagnosed with HCC who underwent surgical resection. Mitochondrial APE1 content, expression levels of the mitochondrial import protein Mia40, and mtDNA damage of tumor tissue and distal non-tumor liver of each patient were analyzed. In parallel, we generated a stable HeLa clone for inducible silencing of endogenous APE1 and re-expression of the recombinant shRNA resistant mitochondrially targeted APE1 form (MTS-APE1). We evaluated mtDNA damage, cell growth, and mitochondrial respiration.

RESULTS

APE1's cytoplasmic positivity in Grades 1 and 2 HCC patients showed a significantly higher expression of mitochondrial APE1, which accounted for lower levels of mtDNA damage observed in the tumor tissue with respect to the distal area. In the contrast, the cytoplasmic positivity in Grade 3 was not associated with APE1's mitochondrial accumulation even when accounting for the higher number of mtDNA lesions measured. Loss of APE1 expression negatively affected mitochondrial respiration, cell viability, and proliferation as well as levels of mtDNA damage. Remarkably, the phenotype was efficiently rescued in MTS-APE1 clone, where APE1 is present only within the mitochondrial matrix.

CONCLUSIONS

Our study confirms the prominent role of the mitochondrial form of APE1 in the early stages of HCC development and the relevance of the non-nuclear fraction of APE1 in the disease progression. We have also confirmed overexpression of Mia40 and the role of the MIA pathway in the APE1 import process. Based on our data, inhibition of the APE1 transport by blocking the MIA pathway could represent a new therapeutic approach for reducing mitochondrial metabolism by preventing the efficient repair of mtDNA.

ForestQC: Quality control on genetic variants from next-generation sequencing data using random forest

Next-generation sequencing technology (NGS) enables the discovery of nearly all genetic variants present in a genome. A subset of these variants, however, may have poor sequencing quality due to limitations in NGS or variant callers. In genetic studies that analyze a large number of sequenced individuals, it is critical to detect and remove those variants with poor quality as they may cause spurious findings. In this paper, we present ForestQC, a statistical tool for performing quality control on variants identified from NGS data by combining a traditional filtering approach and a machine learning approach. Our software uses the information on sequencing quality, such as sequencing depth, genotyping quality, and GC contents, to predict whether a particular variant is likely to be false-positive. To evaluate ForestQC, we applied it to two whole-genome sequencing datasets where one dataset consists of related individuals from families while the other consists of unrelated individuals. Results indicate that ForestQC outperforms widely used methods for performing quality control on variants such as VQSR of GATK by considerably improving the quality of variants to be included in the analysis. ForestQC is also very efficient, and hence can be applied to large sequencing datasets. We conclude that combining a machine learning algorithm trained with sequencing quality information and the filtering approach is a practical approach to perform quality control on genetic variants from sequencing data.

DNA DIVERGENCE AMONG HOMINOIDS

We have determined the degree of single-copy DNA divergence among the extant members of the Hominoidea employing the technique of DNA-DNA hybridization. The species studied include humans, two species of chimpanzees, gorillas, two subspecies of orangutans, and two species of gibbons; as an outgroup we have used a member of the Old World monkeys (Cercopithecidae), the baboon. Our methods are different from those previously used and allow us to control for two factors other than base-pair mismatch that can affect the thermal stability of DNA duplexes: the base composition and duplex length. In addition, we have studied more than one individual for most species and thus are able to assess the effect of intraspecific variation on phylogenetic conclusions. The results indicate that the closest extant relatives of humans are the chimpanzees. Gorillas are the next closest, followed by orangutans and gibbons. This result is strongly supported statistically, as there is virtually no overlap in measurements between different taxa. Our conclusions are in agreement with a growing amount of molecular evidence supporting this pattern of relatedness. The data behave as a reasonably good molecular clock, and we do not see an indication of slowdown in molecular evolution in the clade containing humans and African apes, contrary to what has been documented for protein-coding regions. Because of the clocklike nature of the results, we have estimated that the divergence of humans and chimpanzees occurred about 6-8 million years ago. Results from orangutans indicate that the Borneo and Sumatra populations are genetically distinct, about as different as the named species of chimpanzees.

Genome measures used for quality control are dependent on gene function and ancestry

MOTIVATION

The transition/transversion (Ti/Tv) ratio and heterozygous/nonreference-homozygous (het/nonref-hom) ratio have been commonly computed in genetic studies as a quality control (QC) measurement. Additionally, these two ratios are helpful in our understanding of the patterns of DNA sequence evolution.

RESULTS

To thoroughly understand these two genomic measures, we performed a study using 1000 Genomes Project (1000G) released genotype data (N=1092). An additional two datasets (N=581 and N=6) were used to validate our findings from the 1000G dataset. We compared the two ratios among continental ancestry, genome regions and gene functionality. We found that the Ti/Tv ratio can be used as a quality indicator for single nucleotide polymorphisms inferred from high-throughput sequencing data. The Ti/Tv ratio varies greatly by genome region and functionality, but not by ancestry. The het/nonref-hom ratio varies greatly by ancestry, but not by genome regions and functionality. Furthermore, extreme guanine + cytosine content (either high or low) is negatively associated with the Ti/Tv ratio magnitude. Thus, when performing QC assessment using these two measures, care must be taken to apply the correct thresholds based on ancestry and genome region. Failure to take these considerations into account at the QC stage will bias any following analysis.

CONTACT

yan.guo@vanderbilt.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Three-stage quality control strategies for DNA re-sequencing data

Advances in next-generation sequencing (NGS) technologies have greatly improved our ability to detect genomic variants for biomedical research. In particular, NGS technologies have been recently applied with great success to the discovery of mutations associated with the growth of various tumours and in rare Mendelian diseases. The advance in NGS technologies has also created significant challenges in bioinformatics. One of the major challenges is quality control of the sequencing data. In this review, we discuss the proper quality control procedures and parameters for Illumina technology-based human DNA re-sequencing at three different stages of sequencing: raw data, alignment and variant calling. Monitoring quality control metrics at each of the three stages of NGS data provides unique and independent evaluations of data quality from differing perspectives. Properly conducting quality control protocols at all three stages and correctly interpreting the quality control results are crucial to ensure a successful and meaningful study.

The use of next generation sequencing technology to study the effect of radiation therapy on mitochondrial DNA mutation

The human mitochondrial genome has an exclusively maternal mode of inheritance. Mitochondrial DNA (mtDNA) is particularly vulnerable to environmental insults due in part to an underdeveloped DNA repair system, limited to base excision and homologous recombination repair. Radiation exposure to the ovaries may cause mtDNA mutations in oocytes, which may in turn be transmitted to offspring. We hypothesized that the children of female cancer survivors who received radiation therapy may have an increased rate of mtDNA heteroplasmy mutations, which conceivably could increase their risk of developing cancer and other diseases. We evaluated 44 DNA blood samples from 17 Danish and 1 Finnish families (18 mothers and 26 children). All mothers had been treated for cancer as children and radiation doses to their ovaries were determined based on medical records and computational models. DNA samples were sequenced for the entire mitochondrial genome using the Illumina GAII system. Mother's age at sample collection was positively correlated with mtDNA heteroplasmy mutations. There was evidence of heteroplasmy inheritance in that 9 of the 18 families had at least one child who inherited at least one heteroplasmy site from his or her mother. No significant difference in single nucleotide polymorphisms between mother and offspring, however, was observed. Radiation therapy dose to ovaries also was not significantly associated with the heteroplasmy mutation rate among mothers and children. No evidence was found that radiotherapy for pediatric cancer is associated with the mitochondrial genome mutation rate in female cancer survivors and their children.

SeqVis: a tool for detecting compositional heterogeneity among aligned nucleotide sequences

Compositional heterogeneity is a poorly appreciated attribute of aligned nucleotide and amino acid sequences. It is a common property of molecular phylogenetic data, and it has been found to occur across sequences and/or across sites. Most molecular phylogenetic methods assume that the sequences have evolved under globally stationary, reversible, and homogeneous conditions, implying that the sequences should be compositionally homogeneous. The presence of the above-mentioned compositional heterogeneity implies that the sequences must have evolved under more general conditions than is commonly assumed. Consequently, there is a need for reliable methods to detect under what conditions alignments of nucleotides or amino acids may have evolved. In this chapter, we describe one such program. SeqVis is designed to survey aligned nucleotide sequences. We discuss pros-et-cons of this program in the context of other methods to detect compositional heterogeneity and violated phylogenetic assumptions. The benefits provided by SeqVis are demonstrated in two studies of alignments of nucleotides, one of which contained 7542 nucleotides from 53 species.

The diversity present in 5140 human mitochondrial genomes

We analyzed the current status (as of the end of August 2008) of human mitochondrial genomes deposited in GenBank, amounting to 5140 complete or coding-region sequences, in order to present an overall picture of the diversity present in the mitochondrial DNA of the global human population. To perform this task, we developed mtDNA-GeneSyn, a computer tool that identifies and exhaustedly classifies the diversity present in large genetic data sets. The diversity observed in the 5140 human mitochondrial genomes was compared with all possible transitions and transversions from the standard human mitochondrial reference genome. This comparison showed that tRNA and rRNA secondary structures have a large effect in limiting the diversity of the human mitochondrial sequences, whereas for the protein-coding genes there is a bias toward less variation at the second codon positions. The analysis of the observed amino acid variations showed a tolerance of variations that convert between the amino acids V, I, A, M, and T. This defines a group of amino acids with similar chemical properties that can interconvert by a single transition.

Ratios of radical to conservative amino acid replacement are affected by mutational and compositional factors and may not be indicative of positive Darwinian selection

The ratio of radical to conservative amino acid replacements is frequently used to infer positive Darwinian selection. This method is based on the assumption that radical replacements are more likely than conservative replacements to improve the function of a protein. Therefore, if positive selection plays a major role in the evolution of a protein, one would expect the radical-conservative ratio to exceed the expectation under neutrality. Here, we investigate the possibility that factors unrelated to selection, i.e., transition-transversion ratio, codon usage, genetic code, and amino acid composition, influence the radical-conservative replacement ratio. All factors that have been studied were found to affect the radical-conservative replacement ratio. In particular, amino acid composition and transition-transversion ratio are shown to have the most profound effects. Because none of the studied factors had anything to do with selection (positive or otherwise) and also because all of them (singly or in combination) affected a measure that was supposed to be indicative of positive selection, we conclude that selectional inferences based on radical-conservative replacement ratios should be treated with suspicion.

Expanding the functional human mitochondrial DNA database by the establishment of primate xenomitochondrial cybrids

The nuclear and mitochondrial genomes coevolve to optimize approximately 100 different interactions necessary for an efficient ATP-generating system. This coevolution led to a species-specific compatibility between these genomes. We introduced mitochondrial DNA (mtDNA) from different primates into mtDNA-less human cells and selected for growth of cells with a functional oxidative phosphorylation system. mtDNA from common chimpanzee, pigmy chimpanzee, and gorilla were able to restore oxidative phosphorylation in the context of a human nuclear background, whereas mtDNA from orangutan, and species representative of Old-World monkeys, New-World monkeys, and lemurs were not. Oxygen consumption, a sensitive index of respiratory function, showed that mtDNA from chimpanzee, pigmy chimpanzee, and gorilla replaced the human mtDNA and restored respiration to essentially normal levels. Mitochondrial protein synthesis was also unaltered in successful "xenomitochondrial cybrids." The abrupt failure of mtDNA from primate species that diverged from humans as recently as 8-18 million years ago to functionally replace human mtDNA suggests the presence of one or a few mutations affecting critical nuclear-mitochondrial genome interactions between these species. These cellular systems provide a demonstration of intergenus mtDNA transfer, expand more than 20-fold the number of mtDNA polymorphisms that can be analyzed in a human nuclear background, and provide a novel model for the study of nuclear-mitochondrial interactions.

Glutamine synthetase gene evolution: a good molecular clock

Glutamine synthetase (EC 6.3.1.2) gene evolution in various animals, plants, and bacteria was evaluated by a general stationary Markov model. The evolutionary process proved to be unexpectedly regular even for a time span as long as that between the divergence of prokaryotes from eukaryotes. This enabled us to draw phylogenetic trees for species whose phylogeny cannot be easily reconstructed from the fossil record. Our calculation of the times of divergence of the various organelle-specific enzymes led us to hypothesize that the pea and bean chloroplast genes for these enzymes originated from the duplication of nuclear genes as a result of the different metabolic needs of the various species. Our data indicate that the duplication of plastid glutamine synthetase genes occurred long after the endosymbiotic events that produced the organelles themselves.

Mitochondrial genome in animal cells. Structure, organization, and evolution

In the past decade, the development of new DNA, RNA, and protein technologies has greatly incremented the knowledge about the organization and expression of mitochondrial DNA. The complete base sequence of mitochondrial DNA of several animals is known and many data are rapidly accumulating on the mitochondrial genomes of other systems. Here we discuss the results so far obtained that disclosed unexpected features of mitochondrial genetics. Furthermore, mitochondrial DNA has become established as a powerful tool for evolutionary studies in animals. Evidences are presented demonstrating that the evolution of mitochondrial DNA has proceeded in different ways in the various taxonomic groups. Data on heteroplasmic animals, which demonstrate the rapid evolution of mitochondrial DNA, are also presented.

A simple quantitative model of the molecular clock

We present the ideas, and their motivation, at the basis of a simple model of nucleic acid evolution: the stationary Markov process, or Markov clock. After a brief review of its relevant mathematical properties, the Markov clock is applied to nucleotide sequences from mitochondrial and nuclear genes of different species. Particular emphasis is given to the necessity of carrying out a correct statistical analysis, which allows us to check quantitatively the applicability of our model. We find evidence that the Markov clock ticks in many different processes, and that its limitations can be understood in terms of a simple idea that we call the "base-drift" hypothesis. This hypothesis correlates the deviations from the stationarity of the Markov process to the evolutionary distance dAB(p) of two species A and B, relative to the process P. We conclude by discussing the implications of our findings for future work.

DNA hybridization evidence of hominoid phylogeny: results from an expanded data set

The living hominoids are human, the two species of chimpanzees, gorilla, orangutan, and nine species of gibbons. The cercopithecoids (Old World monkeys) are the sister group of the hominoids. A consensus about the phylogeny of the hominoids has been reached for the branching order of the gibbons (earliest) and the orangutan (next earliest), but the branching order among gorilla, chimpanzees, and human remains in contention. In 1984 we presented DNA-DNA hybridization data, based on 183 DNA hybrids, that we interpreted as evidence that the branching order, from oldest to most recent, was gibbons, orangutan, gorilla, chimpanzees, and human. In the present paper we report on an expanded data set totaling 514 DNA hybrids, which supports the branching order given above. The ranges for the datings of divergence nodes are Old World monkeys, 25-34 million years (Myr) ago; gibbons, 16.4-23 Myr ago; orangutan, 12.2-17 Myr ago; gorilla, 7.7-11 Myr ago; chimpanzees-human, 5.5-7.7 Myr ago. The possible effects of differences in age at first breeding are discussed, and some speculations about average genomic rates of evolution are presented.