The use of next generation sequencing technology to study the effect of radiation therapy on mitochondrial DNA mutation. Mutat Res. 2012;744:154-160. [PMID: 22387842 DOI: 10.1016/j.mrgentox.2012.02.006]

The quality and detection limits of mitochondrial heteroplasmy by long read nanopore sequencing

This study evaluates long-read and short-read sequencing for mitochondrial DNA (mtDNA) heteroplasmy detection. 592,315 bootstrapped datasets generated from two single-nucleotide mismatched ultra-deep sequenced mtDNA samples were used to assess basecalling error and accuracy, limit of heteroplasmy detection, and heteroplasmy detection across various coverage depths. Results showed high Phred scores of data with GC-rich sequence bias for long reads. Limit of detection of 12% heteroplasmy was identified, showing strong correlation (R2 ≥ 0.955) with expected heteroplasmy but underreporting tendency of high-level variants. Nanopore sequencing shows potential for direct applicability in mitochondrial diseases diagnostics, but stringent validation processes to ensure diagnostic result quality are required.

Genetic Diversity, Heteroplasmy, and Recombination in Mitochondrial Genomes of Daphnia pulex, Daphnia pulicaria, and Daphnia obtusa

Genetic variants of mitochondrial DNA at the individual (heteroplasmy) and population (polymorphism) levels provide insight into their roles in multiple cellular and evolutionary processes. However, owing to the paucity of genome-wide data at the within-individual and population levels, the broad patterns of these two forms of variation remain poorly understood. Here, we analyze 1,804 complete mitochondrial genome sequences from Daphnia pulex, Daphnia pulicaria, and Daphnia obtusa. Extensive heteroplasmy is observed in D. obtusa, where the high level of intraclonal divergence must have resulted from a biparental-inheritance event, and recombination in the mitochondrial genome is apparent, although perhaps not widespread. Global samples of D. pulex reveal remarkably low mitochondrial effective population sizes, <3% of those for the nuclear genome. In addition, levels of population diversity in mitochondrial and nuclear genomes are uncorrelated across populations, suggesting an idiosyncratic evolutionary history of mitochondria in D. pulex. These population-genetic features appear to be a consequence of background selection associated with highly deleterious mutations arising in the strongly linked mitochondrial genome, which is consistent with polymorphism and divergence data suggesting a predominance of strong purifying selection. Nonetheless, the fixation of mildly deleterious mutations in the mitochondrial genome also appears to be driving positive selection on genes encoded in the nuclear genome whose products are deployed in the mitochondrion.

UMI-Varcal: A Low-Frequency Variant Caller for UMI-Tagged Paired-End Sequencing Data

The rapid transition from traditional sequencing methods to Next-Generation Sequencing (NGS) has allowed for a faster and more accurate detection of somatic variants (Single-Nucleotide Variant (SNV) and Copy Number Variation (CNV)) in tumor cells. NGS technologies require a succession of steps during which false variants can be silently added at low frequencies. Filtering these artifacts can be a rather difficult task especially when the experiments are designed to look for very low frequency variants. Recently, adding unique molecular barcodes called UMI (Unique Molecular Identifier) to the DNA fragments appears to be a very effective strategy to specifically filter out false variants from the variant calling results (Kukita et al. DNA Res 22(4):269-277, 2015; Newman et al. Nat Biotechnol 34(5):547-555, 2016; Schmitt et al. Proc Natl Acad Sci U S A 109(36):14508-14513). Here, we describe UMI-VarCal (Sater et al. Bioinformatics 36:2718-2724, 2020), which can use the UMI information from UMI-tagged reads to offer a faster and more accurate variant calling analysis.

Characterization of the Different Levels of Variation in 45S rRNA Genes

Although commonly regarded as nearly identical in sequence, 45S ribosomal RNA genes represent a massive source of genetic variation at different levels. Unfortunately, due to their repetitive nature and the difficulty to assemble their multiple copies in tandem, these important genomic elements remain largely unexplored in sequencing projects. Here, I describe how to exploit next generation sequencing data to estimate their copy number in an organism and detect true polymorphic sites within and among individuals. Furthermore, for species that carry multiple 45S ribosomal RNA gene clusters, I show how to make use of experimental populations to assign some of these variants to their cluster of origin.

Genetic landscape of human mitochondrial genome using whole genome sequencing

Increasing evidences suggest that mitochondrial dysfunction is implicated in diseases and aging, and whole-genome sequencing (WGS) is the most unbiased method in analyzing the mitochondrial genome (mtDNA). However, the genetic landscape of mtDNA in the Chinese population has not been fully examined. Here, we described the genetic landscape of mtDNA using WGS data from Chinese individuals (n = 3241). We identified 3892 mtDNA variants, of which 3349 (86%) were rare variants. Interestingly, we observed a trend toward extreme heterogeneity of mtDNA variants. Our study observed a distinct purifying selection on mtDNA, which inhibits the accumulation of harmful heteroplasmies at the individual level: (1) mitochondrial dN/dS ratios were much less than 1; (2) the dN/dS ratio of heteroplasmies was higher than homoplasmies; (3) heteroplasmies had more indels and predicted deleterious variants than homoplasmies. Furthermore, we found that haplogroup M (20.27%) and D (20.15%) had the highest frequencies in the Chinese population, followed by B (18.51%) and F (16.45%). The number of variants per individual differed across haplogroup groups, with a higher number of homoplasmies for the M lineage. Meanwhile, mtDNA copy number was negatively correlated with age but positively correlated with the female sex. Finally, we developed an mtDNA variation database of Chinese populations called MTCards (http://genemed.tech/mtcards/) to facilitate the query of mtDNA variants in this study. In summary, these findings contribute to different aspects of understanding mtDNA, providing a better understanding of the genetic basis of mitochondrial-related diseases.

Ideafix: a decision tree-based method for the refinement of variants in FFPE DNA sequencing data

Increasingly, treatment decisions for cancer patients are being made from next-generation sequencing results generated from formalin-fixed and paraffin-embedded (FFPE) biopsies. However, this material is prone to sequence artefacts that cannot be easily identified. In order to address this issue, we designed a machine learning-based algorithm to identify these artefacts using data from >1 600 000 variants from 27 paired FFPE and fresh-frozen breast cancer samples. Using these data, we assembled a series of variant features and evaluated the classification performance of five machine learning algorithms. Using leave-one-sample-out cross-validation, we found that XGBoost (extreme gradient boosting) and random forest obtained AUC (area under the receiver operating characteristic curve) values >0.86. Performance was further tested using two independent datasets that resulted in AUC values of 0.96, whereas a comparison with previously published tools resulted in a maximum AUC value of 0.92. The most discriminating features were read pair orientation bias, genomic context and variant allele frequency. In summary, our results show a promising future for the use of these samples in molecular testing. We built the algorithm into an R package called Ideafix (DEAmination FIXing) that is freely available at https://github.com/mmaitenat/ideafix.

UMI-VarCal: a new UMI-based variant caller that efficiently improves low-frequency variant detection in paired-end sequencing NGS libraries

MOTIVATION

Next-generation sequencing has become the go-to standard method for the detection of single-nucleotide variants in tumor cells. The use of such technologies requires a PCR amplification step and a sequencing step, steps in which artifacts are introduced at very low frequencies. These artifacts are often confused with true low-frequency variants that can be found in tumor cells and cell-free DNA. The recent use of unique molecular identifiers (UMI) in targeted sequencing protocols has offered a trustworthy approach to filter out artefactual variants and accurately call low-frequency variants. However, the integration of UMI analysis in the variant calling process led to developing tools that are significantly slower and more memory consuming than raw-reads-based variant callers.

RESULTS

We present UMI-VarCal, a UMI-based variant caller for targeted sequencing data with better sensitivity compared to other variant callers. Being developed with performance in mind, UMI-VarCal stands out from the crowd by being one of the few variant callers that do not rely on SAMtools to do their pileup. Instead, at its core runs an innovative homemade pileup algorithm specifically designed to treat the UMI tags in the reads. After the pileup, a Poisson statistical test is applied at every position to determine if the frequency of the variant is significantly higher than the background error noise. Finally, an analysis of UMI tags is performed, a strand bias and a homopolymer length filter are applied to achieve better accuracy. We illustrate the results obtained using UMI-VarCal through the sequencing of tumor samples and we show how UMI-VarCal is both faster and more sensitive than other publicly available solutions.

AVAILABILITY AND IMPLEMENTATION

The entire pipeline is available at https://gitlab.com/vincent-sater/umi-varcal-master under MIT license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Likelihood of Adverse Pregnancy Outcomes and Genetic Disease (Transgenerational Effects) from Exposure to Radioactive Fallout from the 1945 Trinity Atomic Bomb Test

The potential health consequences of the Trinity nuclear weapon test of 16 July 1945 at Alamogordo, New Mexico, are challenging to assess. Population data are available for mortality but not for cancer incidence for New Mexico residents for the first 25 y after the test, and the estimates of radiation dose to the nearby population are lower than the cumulative dose received from ubiquitous natural background radiation. Despite the estimates of low population exposures, it is believed by some that cancer rates in counties near the Trinity test site (located in Socorro County) are elevated compared with other locations across the state. Further, there is a concern about adverse pregnancy outcomes and genetic diseases (transgenerational or heritable effects) related to population exposure to fallout radiation. The possibility of an intergenerational effect has long been a concern of exposed populations, e.g., Japanese atomic bomb survivors, survivors of childhood and adolescent cancer, radiation workers, and environmentally exposed groups. In this paper, the likelihood of discernible transgenerational effects is discounted because (1) in all large-scale comprehensive studies of exposed populations, no heritable genetic effects have been demonstrated in children of exposed parents; (2) the distribution of estimated doses from Trinity is much lower than in other studied populations where no transgenerational effects have been observed; and (3) there is no evidence of increased cancer rates among the scientific, military, and professional participants at the Trinity test and at other nuclear weapons tests who received much higher doses than New Mexico residents living downwind of the Trinity site.

Role of Mitochondrial DNA (mtDNA) Variations in Cancer Development: A Systematic Review

mtDNA is the closed circular, ds-DNA present in mitochondria of eukaryotic cells and are inherited maternally. Besides being the power house of the cell, mitochondria are also responsible for the regulation of redox homeostasis, signaling, metabolism, immunity, survival and apoptosis. Lack of a 'Systematic Review' on mtDNA variations and cancers encouraged us to perform the present study. Pubmed', 'Embase' and 'Cochrane Library' databases were searched using keywords 'Mitochondrial DNA' OR 'mtDNA' OR 'mDNA' AND 'polymorphism' AND 'cancer' AND 'risk' to retrieve literature. Polymorphisms occupy first rank among mtDNA variations followed by CNV, MSI, mutations and hold a great potential to emerge as key predictors for human cancers.

The aging heart

As the elderly segment of the world population increases, it is critical to understand the changes in cardiac structure and function during the normal aging process. In this review, we outline the key molecular pathways and cellular processes that underlie the phenotypic changes in the heart and vasculature that accompany aging. Reduced autophagy, increased mitochondrial oxidative stress, telomere attrition, altered signaling in insulin-like growth factor, growth differentiation factor 11, and 5'- AMP-activated protein kinase pathways are among the key molecular mechanisms underlying cardiac aging. Aging promotes structural and functional changes in the atria, ventricles, valves, myocardium, pericardium, the cardiac conduction system, and the vasculature. We highlight the factors known to accelerate and attenuate the intrinsic aging of the heart and vessels in addition to potential preventive and therapeutic avenues. A greater understanding of the processes involved in cardiac aging may facilitate our ability to mitigate the escalating burden of CVD in older individuals and promote healthy cardiac aging.

The first in vivo multiparametric comparison of different radiation exposure biomarkers in human blood

The increasing risk of acute large-scale radiological/nuclear exposures of population underlines the necessity of developing new, rapid and high throughput biodosimetric tools for estimation of received dose and initial triage. We aimed to compare the induction and persistence of different radiation exposure biomarkers in human peripheral blood in vivo. Blood samples of patients with indicated radiotherapy (RT) undergoing partial body irradiation (PBI) were obtained soon before the first treatment and then after 24 h, 48 h, and 5 weeks; i.e. after 1, 2, and 25 fractionated RT procedures. We collected circulating peripheral blood from ten patients with tumor of endometrium (1.8 Gy per fraction) and eight patients with tumor of head and neck (2.0–2.121 Gy per fraction). Incidence of dicentrics and micronuclei was monitored as well as determination of apoptosis and the transcription level of selected radiation-responsive genes. Since mitochondrial DNA (mtDNA) has been reported to be a potential indicator of radiation damage in vitro, we also assessed mtDNA content and deletions by novel multiplex quantitative PCR. Cytogenetic data confirmed linear dose-dependent increase in dicentrics (p < 0.01) and micronuclei (p < 0.001) in peripheral blood mononuclear cells after PBI. Significant up-regulations of five previously identified transcriptional biomarkers of radiation exposure (PHPT1, CCNG1, CDKN1A, GADD45, and SESN1) were also found (p < 0.01). No statistical change in mtDNA deletion levels was detected; however, our data indicate that the total mtDNA content decreased with increasing number of RT fractions. Interestingly, the number of micronuclei appears to correlate with late radiation toxicity (r² = 0.9025) in endometrial patients suggesting the possibility of predicting the severity of RT-related toxicity by monitoring this parameter. Overall, these data represent, to our best knowledge, the first study providing a multiparametric comparison of radiation biomarkers in human blood in vivo, which have potential for improving biological dosimetry.

Recent Advances in Detecting Mitochondrial DNA Heteroplasmic Variations

The co-existence of wild-type and mutated mitochondrial DNA (mtDNA) molecules termed heteroplasmy becomes a research hot point of mitochondria. In this review, we listed several methods of mtDNA heteroplasmy research, including the enrichment of mtDNA and the way of calling heteroplasmic variations. At the present, while calling the novel ultra-low level heteroplasmy, high-throughput sequencing method is dominant while the detection limit of recorded mutations is accurate to 0.01% using the other quantitative approaches. In the future, the studies of mtDNA heteroplasmy may pay more attention to the single-cell level and focus on the linkage of mutations.

Mitochondrial mutations drive prostate cancer aggression

Nuclear mutations are well known to drive tumor incidence, aggression and response to therapy. By contrast, the frequency and roles of mutations in the maternally inherited mitochondrial genome are poorly understood. Here we sequence the mitochondrial genomes of 384 localized prostate cancer patients, and identify a median of one mitochondrial single-nucleotide variant (mtSNV) per patient. Some of these mtSNVs occur in recurrent mutational hotspots and associate with aggressive disease. Younger patients have fewer mtSNVs than those who diagnosed at an older age. We demonstrate strong links between mitochondrial and nuclear mutational profiles, with co-occurrence between specific mutations. For example, certain control region mtSNVs co-occur with gain of the MYC oncogene, and these mutations are jointly associated with patient survival. These data demonstrate frequent mitochondrial mutation in prostate cancer, and suggest interplay between nuclear and mitochondrial mutational profiles in prostate cancer.

In prostate cancer, the role of mutations in the maternally-inherited mitochondrial genome are not well known. Here, the authors demonstrate frequent, age-dependent mitochondrial mutation in prostate cancer. Strong links between mitochondrial and nuclear mutational profiles are associated with clinical aggressivity.

A comprehensive analysis of mitochondrial genes variants and their association with antipsychotic-induced weight gain

Antipsychotic Induced Weight Gain (AIWG) is a common and severe side effect of many antipsychotic medications. Mitochondria play a vital role for whole-body energy homeostasis and there is increasing evidence that antipsychotics modulate mitochondrial function. This study aimed to examine the role of variants in nuclear-encoded mitochondrial genes and the mitochondrial DNA (mtDNA) in conferring risk for AIWG. We selected 168 European-Caucasian individuals from the CATIE sample based upon meeting criteria of multiple weight measures while taking selected antipsychotics (risperidone, quetiapine or olanzapine). We tested the association of 670 nuclear-encoded mitochondrial genes with weight change (%) using MAGMA software. Thirty of these genes showed nominally significant P-values (<0.05). We were able to replicate the association of three genes, CLPB, PARL, and ACAD10, with weight change (%) in an independent prospectively assessed AIWG sample. We analyzed mtDNA variants in a subset of 74 of these individuals using next-generation sequencing. No common or rare mtDNA variants were found to be significantly associated with weight change (%) in our sample. Additionally, analysis of mitochondrial haplogroups showed no association with weight change (%). In conclusion, our findings suggest nuclear-encoded mitochondrial genes play a role in AIWG. Replication in larger sample is required to validate our initial report of mtDNA variants in AIWG.

Bivariate Poisson models with varying offsets: an application to the paired mitochondrial DNA dataset

To assess the effect of chemotherapy on mitochondrial genome mutations in cancer survivors and their offspring, a study sequenced the full mitochondrial genome and determined the mitochondrial DNA heteroplasmic (mtDNA) mutation rate. To build a model for counts of heteroplasmic mutations in mothers and their offspring, bivariate Poisson regression was used to examine the relationship between mutation count and clinical information while accounting for the paired correlation. However, if the sequencing depth is not adequate, a limited fraction of the mtDNA will be available for variant calling. The classical bivariate Poisson regression model treats the offset term as equal within pairs; thus, it cannot be applied directly. In this research, we propose an extended bivariate Poisson regression model that has a more general offset term to adjust the length of the accessible genome for each observation. We evaluate the performance of the proposed method with comprehensive simulations, and the results show that the regression model provides unbiased parameter estimations. The use of the model is also demonstrated using the paired mtDNA dataset.

Improvements and impacts of GRCh38 human reference on high throughput sequencing data analysis

Analyses of high throughput sequencing data starts with alignment against a reference genome, which is the foundation for all re-sequencing data analyses. Each new release of the human reference genome has been augmented with improved accuracy and completeness. It is presumed that the latest release of human reference genome, GRCh38 will contribute more to high throughput sequencing data analysis by providing more accuracy. But the amount of improvement has not yet been quantified. We conducted a study to compare the genomic analysis results between the GRCh38 reference and its predecessor GRCh37. Through analyses of alignment, single nucleotide polymorphisms, small insertion/deletions, copy number and structural variants, we show that GRCh38 offers overall more accurate analysis of human sequencing data. More importantly, GRCh38 produced fewer false positive structural variants. In conclusion, GRCh38 is an improvement over GRCh37 not only from the genome assembly aspect, but also yields more reliable genomic analysis results.

Supply and Demand of Energy in the Oocyte and the Role of Mitochondria

The sole purpose of any mammalian oocyte is to combine with a spermatozoon and form a viable embryo that implants into the uterus and forms a viable foetus. Most of the structures and mechanisms for this reside within the oocyte itself. The sperm limits itself to fertilisation of the oocyte; apart from this, its only contribution is the male genome and the centrosome, required for cell division. Both intrinsic and extrinsic factors determine the formation of a viable embryo. However, the fundamental necessity for successful reproduction resides within the capacity for the developing embryo to generate sufficient levels of energy for optimal development to occur. Energy is generated principally within mitochondria. In this chapter, we discuss some of the fundamental processes of preimplantation embryo development and the role of mitochondria in providing sufficient energy for the successful completion of these processes. We discuss mitochondrial genetics, replication and energy production. Ageing appears to affect the capacity of the mitochondrion to produce sufficient energy to balance the requirements of the embryo. We discuss some of the theories of the effect of maternal age on mitochondrial physiology and the role this plays in reproduction. We propose that maternal age has longer-term effects on individuals than simply on the efficiency of reproduction. We also discuss some of the procedures assisted reproduction has proposed to alleviate the effect of maternal age on reproduction.

Practicability of mitochondrial heteroplasmy detection through an Illumina genotyping array

MOTIVATION

High-throughput genomic data often contain unexpected information that can be mined for alternative applications. Despite the rise of high-throughput sequencing, Illumina genotyping arrays remain a driving force in large scale genetic and epidemiology studies. By processing and analyzing genotyping data of over 100,000 samples genotyped on Illumina genotyping arrays, we discovered evidence that indicates that mitochondrial heteroplasmy can be estimated from the fluorescence intensity data of the array. To verify our hypothesis, we conducted a sequencing validation study.

RESULT

Mitochondrial DNA targeted sequencing was performed on three samples that had been genotyped using the Illumina exome genotyping array. In each sample chosen, one heteroplasmy target was identified from the genotyping array, and sequencing data verified all three putative heteroplasmic sites. The estimated heteroplasmy level difference between that estimated from the genotyping fluorescence intensity and that directly measured from sequencing was 3.2% on average. Our analysis showed that an Illumina genotyping array can accurately and reliably estimate high-level heteroplasmy (>40%); however, intensity data from a genotyping array is not suitable for estimating low level heteroplasmy (<25%).

Obstetric and offspring risks of women's morbid conditions linked to prior anticancer treatments

BACKGROUND

Literature shows the effects of type of cancer and/or anticancer treatment on live birth percentages and/or pregnancy and neonatal complications in female cancer survivors. However, studies analyzing the obstetric and offspring risks of the morbid conditions associated with previous anti-cancer treatments are missing. The present review aims to uncover these risks.

METHODS

A literature search based on publications up to March 2016 identified by PubMed and references cited in relevant articles.

RESULTS

The morbid conditions associated with prior anticancer treatments including chemotherapy, radiotherapy, surgery, and/or hematopoietic stem-cell transplant may induce not only obstetric and neonatal complications but also long-term effects on offspring. Whereas some risks are predominantly evidenced in untreated women others are observed in both treated and untreated women. These risks may be superimposed on those induced by the current women's trend in Western societies to postpone maternity.

CONCLUSIONS

Medical professionals should be aware and inform female cancer survivors wishing to have a child not only of the short- and long-term risks to themselves and their prospective offspring of previous anticancer treatments, fertility-preservation technologies, and pregnancy itself, but also of those risks linked to the morbid conditions induced by prior anticancer treatments. Once female cancer survivors wishing to have a child have been properly informed about the risks of reproduction, they will be best placed to make decisions of whether or not to have a biological or donor-conceived child. In addition, when medical professionals be aware of these risks, they will be also best placed to provide appropriate treatments before/during pregnancy in order to prevent or alleviate the impact of these morbid conditions on maternal and offspring health.

mtDNA-Server: next-generation sequencing data analysis of human mitochondrial DNA in the cloud

Next generation sequencing (NGS) allows investigating mitochondrial DNA (mtDNA) characteristics such as heteroplasmy (i.e. intra-individual sequence variation) to a higher level of detail. While several pipelines for analyzing heteroplasmies exist, issues in usability, accuracy of results and interpreting final data limit their usage. Here we present mtDNA-Server, a scalable web server for the analysis of mtDNA studies of any size with a special focus on usability as well as reliable identification and quantification of heteroplasmic variants. The mtDNA-Server workflow includes parallel read alignment, heteroplasmy detection, artefact or contamination identification, variant annotation as well as several quality control metrics, often neglected in current mtDNA NGS studies. All computational steps are parallelized with Hadoop MapReduce and executed graphically with Cloudgene. We validated the underlying heteroplasmy and contamination detection model by generating four artificial sample mix-ups on two different NGS devices. Our evaluation data shows that mtDNA-Server detects heteroplasmies and artificial recombinations down to the 1% level with perfect specificity and outperforms existing approaches regarding sensitivity. mtDNA-Server is currently able to analyze the 1000G Phase 3 data (n = 2,504) in less than 5 h and is freely accessible at https://mtdna-server.uibk.ac.at.

Mitochondria sequence mapping strategies and practicability of mitochondria variant detection from exome and RNA sequencing data

The rapid progress in high-throughput sequencing has significantly enriched our capacity for studying the mitochondrial DNA (mtDNA). In addition to performing specific mitochondrial targeted sequencing, an increasingly popular alternative approach is using the off-target reads from exome sequencing to infer mtDNA variants, including single nucleotide polymorphisms (SNPs) and heteroplasmy. However, the effectiveness and practicality of this approach have not been tested. Recently, RNAseq data have also been suggested as a good source for alternative data mining, but whether mitochondrial variants can be detected from RNAseq data has not been validated. We designed a study to evaluate the practicability of mtDNA variant detection using exome and RNA sequencing data. Five breast cancer cell lines were sequenced through mitochondrial targeted, exome, and RNA sequencing. Mitochondrial targeted sequencing was used as the gold standard to compute the validation and false discovery rates of SNP and heteroplasmy detection in exome and RNAseq data. We found that exome and RNA sequencing can accurately detect mitochondrial SNPs. However, the lower false discovery rate makes exome sequencing a better choice for heteroplasmy detection than RNAseq. Furthermore, we examined three alignment strategies and found that aligning reads directly to the mitochondrial reference genome or aligning reads to the nuclear and mitochondrial references genomes simultaneously produced the best results, and that aligning to the nuclear genome first and afterwards to the mitochondrial genome performed poorly. In conclusion, our study provides important guidelines for future studies that intend to use either exome sequencing or RNAseq data to infer mitochondrial SNPs and heteroplasmy.

Mitochondrial mutations in subjects with psychiatric disorders

A considerable body of evidence supports the role of mitochondrial dysfunction in psychiatric disorders and mitochondrial DNA (mtDNA) mutations are known to alter brain energy metabolism, neurotransmission, and cause neurodegenerative disorders. Genetic studies focusing on common nuclear genome variants associated with these disorders have produced genome wide significant results but those studies have not directly studied mtDNA variants. The purpose of this study is to investigate, using next generation sequencing, the involvement of mtDNA variation in bipolar disorder, schizophrenia, major depressive disorder, and methamphetamine use. MtDNA extracted from multiple brain regions and blood were sequenced (121 mtDNA samples with an average of 8,800x coverage) and compared to an electronic database containing 26,850 mtDNA genomes. We confirmed novel and rare variants, and confirmed next generation sequencing error hotspots by traditional sequencing and genotyping methods. We observed a significant increase of non-synonymous mutations found in individuals with schizophrenia. Novel and rare non-synonymous mutations were found in psychiatric cases in mtDNA genes: ND6, ATP6, CYTB, and ND2. We also observed mtDNA heteroplasmy in brain at a locus previously associated with schizophrenia (T16519C). Large differences in heteroplasmy levels across brain regions within subjects suggest that somatic mutations accumulate differentially in brain regions. Finally, multiplasmy, a heteroplasmic measure of repeat length, was observed in brain from selective cases at a higher frequency than controls. These results offer support for increased rates of mtDNA substitutions in schizophrenia shown in our prior results. The variable levels of heteroplasmic/multiplasmic somatic mutations that occur in brain may be indicators of genetic instability in mtDNA.

Protection of the gametes embryo/fetus from prenatal radiation exposure

There is no convincing evidence of germline mutation manifest as heritable disease in the offspring of humans attributable to ionizing radiation, yet radiation clearly induces mutations in microbes and somatic cells of rodents and humans. Doses to the embryo estimated to be in the range of 0.15-0.2 Gy during the pre-implantation and pre-somite stages may increase the risk of embryonic loss. However, an increased risk of congenital malformations or growth retardation has not been observed in the surviving embryos. These results are primarily derived from mammalian animal studies and are referred to as the "all-or-none phenomenon." The tissue reaction effects of ionizing radiation (previously referred to as deterministic effects) are congenital malformations, mental retardation, decreased intelligence quotient, microcephaly, neurobehavioral effects, convulsive disorders, growth retardation (height and weight), and embryonic and fetal death (miscarriage, stillbirth). All these effects are consistent with having a threshold dose below which there is no increased risk. The risk of cancer in offspring that have been exposed to diagnostic x-ray procedures while in utero has been debated for 55 y. High doses to the embryo or fetus (e.g., >0.5 Gy) increase the risk of cancer. Most pregnant women exposed to x-ray procedures and other forms of ionizing radiation today received doses to the embryo or fetus <0.1 Gy. The risk of cancer in offspring exposed in utero at exposures <0.1 Gy is controversial and has not been fully resolved. Diagnostic imaging procedures using ionizing radiation that are clinically indicated for the pregnant patient and her fetus should be performed because the clinical benefits outweigh the potential oncogenic risks.

Aging: A mitochondrial DNA perspective, critical analysis and an update

The mitochondrial theory of aging, a mainstream theory of aging which once included accumulation of mitochondrial DNA (mtDNA) damage by reactive oxygen species (ROS) as its cornerstone, has been increasingly losing ground and is undergoing extensive revision due to its inability to explain a growing body of emerging data. Concurrently, the notion of the central role for mtDNA in the aging process is being met with increased skepticism. Our progress in understanding the processes of mtDNA maintenance, repair, damage, and degradation in response to damage has largely refuted the view of mtDNA as being particularly susceptible to ROS-mediated mutagenesis due to its lack of “protective” histones and reduced complement of available DNA repair pathways. Recent research on mitochondrial ROS production has led to the appreciation that mitochondria, even in vitro, produce much less ROS than previously thought, automatically leading to a decreased expectation of physiologically achievable levels of mtDNA damage. New evidence suggests that both experimentally induced oxidative stress and radiation therapy result in very low levels of mtDNA mutagenesis. Recent advances provide evidence against the existence of the “vicious” cycle of mtDNA damage and ROS production. Meta-studies reveal no longevity benefit of increased antioxidant defenses. Simultaneously, exciting new observations from both comparative biology and experimental systems indicate that increased ROS production and oxidative damage to cellular macromolecules, including mtDNA, can be associated with extended longevity. A novel paradigm suggests that increased ROS production in aging may be the result of adaptive signaling rather than a detrimental byproduct of normal respiration that drives aging. Here, we review issues pertaining to the role of mtDNA in aging.

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.

High-throughput sequencing in mitochondrial DNA research

Next-generation sequencing, also known as high-throughput sequencing, has greatly enhanced researchers' ability to conduct biomedical research on all levels. Mitochondrial research has also benefitted greatly from high-throughput sequencing; sequencing technology now allows for screening of all 16,569 base pairs of the mitochondrial genome simultaneously for SNPs and low level heteroplasmy and, in some cases, the estimation of mitochondrial DNA copy number. It is important to realize the full potential of high-throughput sequencing for the advancement of mitochondrial research. To this end, we review how high-throughput sequencing has impacted mitochondrial research in the categories of SNPs, low level heteroplasmy, copy number, and structural variants. We also discuss the different types of mitochondrial DNA sequencing and their pros and cons. Based on previous studies conducted by various groups, we provide strategies for processing mitochondrial DNA sequencing data, including assembly, variant calling, and quality control.

Very low-level heteroplasmy mtDNA variations are inherited in humans

Little is known about the inheritance of very low heteroplasmy mitochondria DNA (mtDNA) variations. Even with the development of new next-generation sequencing methods, the practical lower limit of measured heteroplasmy is still about 1% due to the inherent noise level of the sequencing. In this study, we sequenced the mitochondrial genome of 44 individuals using Illumina high-throughput sequencing technology and obtained high-coverage mitochondria sequencing data. Our study population contains many mother-offspring pairs. This unique study design allows us to bypass the usual heteroplasmy limitation by analyzing the correlation of mutation levels at each position in the mtDNA sequence between maternally related pairs and non-related pairs. The study showed that very low heteroplasmy variants, down to almost 0.1%, are inherited maternally and that this inheritance begins to decrease at about 0.5%, corresponding to a bottleneck of about 200 mtDNA.

Effects of ionizing radiation on mitochondria

The current concept of radiobiology posits that damage to the DNA in the cell nucleus is the primary cause for the detrimental effects of radiation. However, emerging experimental evidence suggests that this theoretical framework is insufficient for describing extranuclear radiation effects, particularly the response of the mitochondria, an important site of extranuclear, coding DNA. Here, we discuss experimental observations of the effects of ionizing radiation on the mitochondria at (1) the DNA and (2) functional levels. The roles of mitochondria in (3) oxidative stress and (4) late radiation effects are discussed. In this review, we summarize the current understanding of targets for ionizing radiation outside the cell nucleus. Available experimental data suggest that an increase in the tumoricidal efficacy of radiation therapy might be achievable by targeting mitochondria. Likewise, more specific protection of mitochondria and its coding DNA should reduce damage to healthy cells exposed to ionizing radiation.

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.

Finding the lost treasures in exome sequencing data

Exome sequencing is one of the most cost-efficient sequencing approaches for conducting genome research on coding regions. However, significant portions of the reads obtained in exome sequencing come from outside of the designed target regions. These additional reads are generally ignored, potentially wasting an important source of genomic data. There are three major types of unintentionally sequenced read that can be found in exome sequencing data: reads in introns and intergenic regions, reads in the mitochondrial genome, and reads originating in viral genomes. All of these can be used for reliable data mining, extending the utility of exome sequencing. Large-scale exome sequencing data repositories, such as The Cancer Genome Atlas (TCGA), the 1000 Genomes Project, National Heart, Lung, and Blood Institute (NHLBI) Exome Sequencing Project, and The Sequence Reads Archive, provide researchers with excellent secondary data-mining opportunities to study genomic data beyond the intended target regions.

MitoSeek: extracting mitochondria information and performing high-throughput mitochondria sequencing analysis

MOTIVATION

Exome capture kits have capture efficiencies that range from 40 to 60%. A significant amount of off-target reads are from the mitochondrial genome. These unintentionally sequenced mitochondrial reads provide unique opportunities to study the mitochondria genome.

RESULTS

MitoSeek is an open-source software tool that can reliably and easily extract mitochondrial genome information from exome and whole genome sequencing data. MitoSeek evaluates mitochondrial genome alignment quality, estimates relative mitochondrial copy numbers and detects heteroplasmy, somatic mutation and structural variants of the mitochondrial genome. MitoSeek can be set up to run in parallel or serial on large exome sequencing datasets.

AVAILABILITY

https://github.com/riverlee/MitoSeek

Comprehensive next-generation sequence analyses of the entire mitochondrial genome reveal new insights into the molecular diagnosis of mitochondrial DNA disorders

PURPOSE

The application of massively parallel sequencing technology to the analysis of the mitochondrial genome has demonstrated great improvement in the molecular diagnosis of mitochondrial DNA-related disorders. The objective of this study was to investigate the performance characteristics and to gain new insights into the analysis of the mitochondrial genome.

METHODS

The entire mitochondrial genome was analyzed as a single amplicon using a long-range PCR-based enrichment approach coupled with massively parallel sequencing. The interference of the nuclear mitochondrial DNA homologs was distinguished from the actual mitochondrial DNA sequences by comparison with the results obtained from conventional PCR-based Sanger sequencing using multiple pairs of primers.

RESULTS

Our results demonstrated the uniform coverage of the entire mitochondrial genome. Massively parallel sequencing of the single amplicon revealed the presence of single-nucleotide polymorphisms and nuclear homologs of mtDNA sequences that cause the erroneous and inaccurate variant calls when PCR/Sanger sequencing approach was used. This single amplicon massively parallel sequencing strategy provides an accurate quantification of mutation heteroplasmy as well as the detection and mapping of mitochondrial DNA deletions.

CONCLUSION

The ability to quantitatively and qualitatively evaluate every single base of the entire mitochondrial genome is indispensible to the accurate molecular diagnosis and genetic counseling of mitochondrial DNA-related disorders. This new approach may be considered as first-line testing for comprehensive analysis of the mitochondrial genome.Genet Med 2013:15(5):388-394.

The effect of strand bias in Illumina short-read sequencing data

BACKGROUND

When using Illumina high throughput short read data, sometimes the genotype inferred from the positive strand and negative strand are significantly different, with one homozygous and the other heterozygous. This phenomenon is known as strand bias. In this study, we used Illumina short-read sequencing data to evaluate the effect of strand bias on genotyping quality, and to explore the possible causes of strand bias.

RESULT

We collected 22 breast cancer samples from 22 patients and sequenced their exome using the Illumina GAIIx machine. By comparing the consistency between the genotypes inferred from this sequencing data with the genotypes inferred from SNP chip data, we found that, when using sequencing data, SNPs with extreme strand bias did not have significantly lower consistency rates compared to SNPs with low or no strand bias. However, this result may be limited by the small subset of SNPs present in both the exome sequencing and the SNP chip data. We further compared the transition and transversion ratio and the number of novel non-synonymous SNPs between the SNPs with low or no strand bias and those with extreme strand bias, and found that SNPs with low or no strand bias have better overall quality. We also discovered that the strand bias occurs randomly at genomic positions across these samples, and observed no consistent pattern of strand bias location across samples. By comparing results from two different aligners, BWA and Bowtie, we found very consistent strand bias patterns. Thus strand bias is unlikely to be caused by alignment artifacts. We successfully replicated our results using two additional independent datasets with different capturing methods and Illumina sequencers.

CONCLUSION

Extreme strand bias indicates a potential high false-positive rate for SNPs.