Single genome retrieval of context-dependent variability in mutation rates for human germline

Generalised interrelations among mutation rates drive the genomic compliance of Chargaff's second parity rule

Chargaff's second parity rule (PR-2), where the complementary base and k-mer contents are matching within the same strand of a double stranded DNA (dsDNA), is a phenomenon that invited many explanations. The strict compliance of nearly all nuclear dsDNA to PR-2 implies that the explanation should also be similarly adamant. In this work, we revisited the possibility of mutation rates driving PR-2 compliance. Starting from the assumption-free approach, we constructed kinetic equations for unconstrained simulations. The results were analysed for their PR-2 compliance by employing symbolic regression and machine learning techniques. We arrived to a generalised set of mutation rate interrelations in place in most species that allow for their full PR-2 compliance. Importantly, our constraints explain PR-2 in genomes out of the scope of the prior explanations based on the equilibration under mutation rates with simpler no-strand-bias constraints. We thus reinstate the role of mutation rates in PR-2 through its molecular core, now shown, under our formulation, to be tolerant to previously noted strand biases and incomplete compositional equilibration. We further investigate the time for any genome to reach PR-2, showing that it is generally earlier than the compositional equilibrium, and well within the age of life on Earth.

8-Oxo-7,8-dihydroguanine, friend and foe: Epigenetic-like regulator versus initiator of mutagenesis

A high flux of reactive oxygen species during oxidative stress results in oxidative modification of cellular components including DNA. Oxidative DNA "damage" to the heterocyclic bases is considered deleterious because polymerases may incorrectly read the modifications causing mutations. A prominent member in this class is the oxidized guanine base 8-oxo-7,8-dihydroguanine (OG) that is moderately mutagenic effecting G→T transversion mutations. Recent reports have identified that formation of OG in G-rich regulatory elements in the promoters of the VEGF, TNFα, and SIRT1 genes can increase transcription via activation of the base excision repair (BER) pathway. Work in our laboratory with the G-rich sequence in the promoter of VEGF concluded that BER drives a shift in structure to a G-quadruplex conformation leading to gene activation in mammalian cells. More specifically, removal of OG from the duplex context by 8-oxoguanine glycosylase 1 (OGG1) produces an abasic site (AP) that destabilizes the duplex, shifting the equilibrium toward the G-quadruplex fold because of preferential extrusion of the AP into a loop. The AP is bound but inefficiently cleaved by apurinic/apyrimidinic endoDNase I (APE1) that likely allows recruitment of activating transcription factors for gene induction. The ability of OG to induce transcription ascribes a regulatory or epigenetic-like role for this oxidatively modified base. We compare OG to the 5-methylcytosine (5mC) epigenetic pathway including its oxidized derivatives, some of which poise genes for transcription while also being substrates for BER. The mutagenic potential of OG to induce only ∼one-third the number of mutations (G→T) compared to deamination of 5mC producing C→T mutations is described. These comparisons blur the line between friendly epigenetic base modifications and those that are foes, i.e. DNA "damage," causing genetic mutations.