Computational DNA hole spectroscopy: A new tool to predict mutation hotspots, critical base pairs, and disease 'driver' mutations

Estimating the Vertical Ionization Potential of Single-Stranded DNA Molecules

The electronic properties of DNA molecules, defined by the sequence-dependent ionization potentials of nucleobases, enable long-range charge transport along the DNA stacks. This has been linked to a range of key physiological processes in the cells and to the triggering of nucleobase substitutions, some of which may cause diseases. To gain molecular-level understanding of the sequence dependence of these phenomena, we estimated the vertical ionization potential (vIP) of all possible nucleobase stacks in B-conformation, containing one to four Gua, Ade, Thy, Cyt, or methylated Cyt. To do this, we used quantum chemistry calculations and more precisely the second-order Møller-Plesset perturbation theory (MP2) and three double-hybrid density functional theory methods, combined with several basis sets for describing atomic orbitals. The calculated vIP of single nucleobases were compared to experimental data and those of nucleobase pairs, triplets, and quadruplets, to observed mutability frequencies in the human genome, reported to be correlated with vIP values. This comparison selected MP2 with the 6-31G* basis set as the best of the tested calculation levels. These results were exploited to set up a recursive model, called vIPer, which estimates the vIP of all possible single-stranded DNA sequences of any length based on the calculated vIPs of overlapping quadruplets. vIPer's vIP values correlate well with oxidation potentials measured by cyclic voltammetry and activities obtained through photoinduced DNA cleavage experiments, further validating our approach. vIPer is freely available on the github.com/3BioCompBio/vIPer repository.

The exposure levels and health risk assessment of antibiotics in urine and its association with platelet mitochondrial DNA methylation in adults from Tianjin, China: A preliminary study

There has been extensive research on antibiotics exposure in adults by biomonitoring, but the biological mechanisms and potential risks to human health remain limited. In this study, 102 adults aged 26-44 years in Tianjin were studied and 23 common antibiotics in urine were analyzed by Liquid chromatography-mass spectrometry (LC-MS). All antibiotics were detected in urine, with an overall detection frequency of 40.4% (the detection frequencies of phenothiazines, quinolones, sulfonamides, tetracyclines, and chloramphenicol were 77%, 54%, 24%, 28%, and 49%, respectively.). Ofloxacin and enrofloxacin had the highest detection frequencies (85% and 81%), with median concentrations of 0.26 (IQR: 0.05-1.36) and 0.09 (IQR: 0.03-0.14) ng/mL, respectively. Based on health risk assessment, the predicted estimated daily exposures (EDEs) ranged from 0 μg/kg/day to 13.98 μg/kg/day. The hazard quotient (HQ) values of all the antibiotics except ofloxacin and ciprofloxacin were bellow one, which are considered safe. For all blood samples, the mitochondrial DNA (mtDNA) methylation levels in the MT-ATP6 (ranging between 3.86% and 34.18%) were slightly higher than MT-ATP8 and MT-ND5 (ranging between 0.57% and 9.32%, 1.08% and 19.62%, respectively). Furthermore, mtDNA methylation from MT-ATP6, MT-ATP8 and MT-ND5 were measured by bisulfite-PCR pyrosequencing. The association (P < 0.05) was found between mtDNA methylation level (MT-ATP8 and MT-ND5) and individual antibiotics including chlorpromazine, ciprofloxacin, enrofloxacin, norfloxacin, pefloxacin, sulfaquinoxaline, sulfachloropyridazine, chloramphenicol, and thiamphenicol, indicating that persistent exposure to low-dose multiple antibiotics may affect the mtDNA methylation level and in turn pose health risks.

Low frequency mitochondrial DNA heteroplasmy SNPs in blood, retina, and [RPE+choroid] of age-related macular degeneration subjects

Purpose

Mitochondrial (mt) DNA damage is associated with age-related macular degeneration (AMD) and other human aging diseases. This study was designed to quantify and characterize mtDNA low-frequency heteroplasmy single nucleotide polymorphisms (SNPs) of three different tissues isolated from AMD subjects using Next Generation Sequencing (NGS) technology.

Methods

DNA was extracted from neural retina, [RPE+choroid] and blood from three deceased age-related macular degeneration (AMD) subjects. Entire mitochondrial genomes were analyzed for low-frequency heteroplasmy SNPs using NGS technology that independently sequenced both mtDNA strands. This deep sequencing method (average sequencing depth of 30,000; range 1,000–100,000) can accurately differentiate low-frequency heteroplasmy SNPs from DNA modification artifacts. Twenty-three ‘hot-spot’ heteroplasmy mtDNA SNPs were analyzed in 222 additional blood samples.

Results

Germline homoplasmy SNPs that defined mtDNA haplogroups were consistent in the three tissues of each subject. Analyses of SNPs with <40% heteroplasmy revealed the blood had significantly greater numbers of heteroplasmy SNPs than retina alone (p≤0.05) or retina+choroid combined (p = 0.008). Twenty-three ‘hot-spot’ mtDNA heteroplasmy SNPs were present, with three being non-synonymous (amino acid change). Four ‘hot-spot’ heteroplasmy SNPs (m.1120C>T, m.1284T>C, m.1556C>T, m.7256C>T) were found in additional samples (n = 222). Five heteroplasmy SNPs (m.4104A>G, m.5320C>T, m.5471G>A, m.5474A>G, m.5498A>G) declined with age. Two heteroplasmy SNPs (m.13095T>C, m.13105A>G) increased in AMD compared to Normal samples. In the heteroplasmy SNPs, very few transversion mutations (purine to pyrimidine or vice versa, associated with oxidative damage) were found and the majority were transition changes (purine to purine or pyrimidine to pyrimidine, associated with replication errors).

Conclusion

Within an individual, the blood, retina and [RPE+choroid] contained identical homoplasmy SNPs representing inherited germline mtDNA haplogroup. NGS methodology showed significantly more mtDNA heteroplasmy SNPs in blood compared to retina and [RPE+choroid], suggesting the latter tissues have substantial protection. Significantly higher heteroplasmy levels of m.13095T>C and m.13105A>G may represent potential AMD biomarkers. Finally, high levels of transition mutations suggest that accumulation of heteroplasmic SNPs may occur through replication errors rather than oxidative damage.

Relation between DNA ionization potentials, single base substitutions and pathogenic variants

Background

It is nowadays clear that single base substitutions that occur in the human genome, of which some lead to pathogenic conditions, are non-random and influenced by their flanking nucleobase sequences. However, despite recent progress, the understanding of these "non-local" effects is still far from being achieved.

Results

To advance this problem, we analyzed the relationship between the base mutability in specific gene regions and the electron hole transport along the DNA base stacks, as it is one of the mechanisms that have been suggested to contribute to these effects. More precisely, we studied the connection between the normalized frequency of single base substitutions and the vertical ionization potential of the base and its flanking sequence, estimated using MP2/6-31G* ab initio quantum chemistry calculations. We found a statistically significant overall anticorrelation between these two quantities: the lower the vIP value, the more probable the substitution. Moreover, the slope of the regression lines varies. It is larger for introns than for exons and untranslated regions, and for synonymous than for missense substitutions. Interestingly, the correlation appears to be more pronounced when considering the flanking sequence of the substituted base in the 3’ rather than in the 5’ direction, which corresponds to the preferred direction of charge migration. A weaker but still statistically significant correlation is found between the ionization potentials and the pathogenicity of the base substitutions. Moreover, pathogenicity is also preferentially associated with larger changes in ionization potentials upon base substitution.

Conclusions

With this analysis we gained new insights into the complex biophysical mechanisms that are at the basis of mutagenesis and pathogenicity, and supported the role of electron-hole transport in these matters.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5867-y) contains supplementary material, which is available to authorized users.

MtBrowse: An integrative genomics browser for human mitochondrial DNA

The human mitochondrion is a unique semi-autonomous organelle with a genome of its own and also requires nuclear encoded components to carry out its functions. In addition to being the powerhouse of the cell, mitochondria plays a central role in several metabolic pathways. It is therefore challenging to delineate the cause-effect relationship in context of mitochondrial dysfunction. Several studies implicate mutations in mitochondrial DNA (mtDNA) in various complex diseases. The human mitochondrial DNA (mtDNA) encodes a set of 37 genes, 13 protein coding, 22 tRNAs and two ribosomal RNAs, which are essential structural and functional components of the electron transport chain. As mentioned above, variations in these genes have been implicated in a broad spectrum of diseases and are extensively reported in literature and various databases. A large number of databases and prediction methods have been published to elucidate the role of human mitochondrial DNA in various disease phenotypes. However, there is no centralized resource to visualize this genotype-phenotype data. Towards this, we have developed MtBrowse: an integrative genomics browser for human mtDNA. As of now, MtBrowse has four categories - Gene, Disease, Reported variation and Variation prediction. These categories have 105 tracks and house data on mitochondrial reference genes, around 600 variants reported in literature with respect to various disease phenotypes and predictions for potential pathogenic variations in protein-coding genes. MtBrowse also hosts genomic variation data from over 5000 individuals on 22 disease phenotypes. MtBrowse may be accessed at http://ab-openlab.csir.res.in/cgi-bin/gb2/gbrowse.

Influence of Electron-Holes on DNA Sequence-Specific Mutation Rates

Biases in mutation rate can influence molecular evolution, yielding rates of evolution that vary widely in different parts of the genome and even among neighboring nucleotides. Here, we explore one possible mechanism of influence on sequence-specific mutation rates, the electron–hole, which can localize and potentially trigger a replication mismatch. A hole is a mobile site of positive charge created during one-electron oxidation by, for example, radiation, contact with a mutagenic agent, or oxidative stress. Its quantum wavelike properties cause it to localize at various sites with probabilities that vary widely, by orders of magnitude, and depend strongly on the local sequence. We find significant correlations between hole probabilities and mutation rates within base triplets, observed in published mutation accumulation experiments on four species of bacteria. We have also computed hole probability spectra for hypervariable segment I of the human mtDNA control region, which contains several mutational hotspots, and for heptanucleotides in noncoding regions of the human genome, whose polymorphism levels have recently been reported. We observe significant correlations between hole probabilities, and context-specific mutation and substitution rates. The correlation with hole probability cannot be explained entirely by CpG methylation in the heptanucleotide data. Peaks in hole probability tend to coincide with mutational hotspots, even in mtDNA where CpG methylation is rare. Our results suggest that hole-enhanced mutational mechanisms, such as oxidation-stabilized tautomerization and base deamination, contribute to molecular evolution.