Search

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.

Our genome contains about half a million sites capable of forming G-quadruplex (G4) structures. Such structural formations, often localised at important regulatory loci, have high capability of altering the predisposition of corresponding genomic spans to endogenous and exogenous DNA damage. In this work, we devised an approach to systematically enrich and zoom onto structure-driven effects on the propensity to undergo 9 types of DNA damage: ultraviolet radiation-induced pyrimidine-pyrimidone (6-4) photoproduct PP and cyclobutane pyrimidine dimer CPD couplings (two dyad-based subtypes in each), cisplatin-mediated G-G crosslinks, reactive oxygen species induced 8-oxoguanine damage, DNA fragmentation upon natural decay and fossilisation, breakages from artificial enzymatic cleavage and ultrasound sonication. Our results indicate that the structural effects on DNA damageability at G4 sites are not a simple combination of shielding (G4 strand) and de-shielding (opposite strand) against damaging factors, and the outcomes have different patterns and variation from one damage type to another, highly dependent on the G4 strength and relative strand localisation. The results are accompanied by electronic structure calculations, detailed structural parallels and considerations.