Genomic compliance with Chargaff's second parity rule may have originated non-adaptively, but stem-loops now function adaptively

Of Chargaff's four rules on DNA base quantity, his second parity rule (PR-2) is the most contentious. Various biometricians (e.g., Sueoka, Lobry) regarded PR-2 compliance as a non-adaptive feature of modern genomes that could be modeled through interrelations among mutation rates. However, PR-2 compliance with stem-loop potential was considered adaptively relevant by biochemists familiar with analyses of nucleic acid structure (e.g., of Crick) and of meiotic recombination (e.g., of Kleckner). Meanwhile, other biometricians had shown that PR-2 complementarity extended beyond individual bases (1-mers) to oligonucleotides (k-mers), possibly reflecting "advantageous DNA structure" (Nussinov). An "introns early" hypothesis (Reanney, Forsdyke) had suggested a primordial nucleic acid world with recombination-mediated error-correction requiring genome-wide stem-loop potential to have evolved prior to localized intrusions of protein-encoding potential (exons). Thus, a primordial genome was equivalent to one long intron. Indeed, when assessed as the base order-dependent component (correcting for local influences of GC%), modern genes, especially when evolving rapidly under positive Darwinian selection, display high intronic stem-loop potential. This suggests forced migration from neighboring exons by competing protein-encoding potential. PR-2 compliance may have first arisen non-adaptively. Primary prototypic structures were later strengthened by their adaptive contribution to recombination. Thus, contentious views may actually be in harmony.

Potential Achilles heels of SARS-CoV-2 are best displayed by the base order-dependent component of RNA folding energy

The base order-dependent component of folding energy has revealed a highly conserved region in HIV-1 genomes that associates with RNA structure. This corresponds to a packaging signal that is recognized by the nucleocapsid domain of the Gag polyprotein. Long viewed as a potential HIV-1 "Achilles heel," the signal can be targeted by a new antiviral compound. Although SARS-CoV-2 differs in many respects from HIV-1, the same technology displays regions with a high base order-dependent folding energy component, which are also highly conserved. This indicates structural invariance (SI) sustained by natural selection. While the regions are often also protein-encoding (e. g. NSP3, ORF3a), we suggest that their nucleic acid level functions can be considered potential "Achilles heels" for SARS-CoV-2, perhaps susceptible to therapies like those envisaged for AIDS. The ribosomal frameshifting element scored well, but higher SI scores were obtained in other regions, including those encoding NSP13 and the nucleocapsid (N) protein.

Neutralism versus selectionism: Chargaff's second parity rule, revisited

Of Chargaff's four "rules" on DNA base frequencies, the functional interpretation of his second parity rule (PR2) is the most contentious. Thermophile base compositions (GC%) were taken by Galtier and Lobry (1997) as favoring Sueoka's neutral PR2 hypothesis over Forsdyke's selective PR2 hypothesis, namely that mutations improving local within-species recombination efficiency had generated a genome-wide potential for the strands of duplex DNA to separate and initiate recombination through the "kissing" of the tips of stem-loops. However, following Chargaff's GC rule, base composition mainly reflects a species-specific, genome-wide, evolutionary pressure. GC% could not have consistently followed the dictates of temperature, since it plays fundamental roles in both sustaining species integrity and, through primarily neutral genome-wide mutation, fostering speciation. Evidence for a local within-species recombination-initiating role of base order was obtained with a novel technology that masked the contribution of base composition to nucleic acid folding energy. Forsdyke's results were consistent with his PR2 hypothesis, appeared to resolve some root problems in biology and provided a theoretical underpinning for alignment-free taxonomic analyses using relative oligonucleotide frequencies (k-mer analysis). Moreover, consistent with Chargaff's cluster rule, discovery of the thermoadaptive role of the "purine-loading" of open reading frames made less tenable the Galtier-Lobry anti-selectionist arguments.

Mutations in disordered proteins as early indicators of nucleic acid changes triggering speciation

In this study, we analyze the role of different structural variants of proteins in the speciation processes. We separate human and mouse proteomes (taken as a reference) into three previously defined variants of disorder: ordered proteins (ORDPs), structured proteins with intrinsically disordered protein regions (IDPRs), and intrinsically disordered proteins (IDPs). Then, using the representation we call here Forsdyke plot, we study the correlation of DNA divergence with the corresponding protein (phenotypic) divergence in the three variants, comparing human and mouse coding sequences with their homologs from 26 eukaryotes. The parameters of the correlation are related to the speciation process. We find that the three variants of disordered proteins are differently related to the speciation process. Specifically, IDPs phenotypically diverge earlier than ORDPs and IDPRs. ORDPs diverge later but are phenotypically more reactive to nucleotide mutations than IDPRs and IDPs. Finally, IDPRs appear to diverge phenotypically later than IDPs, like ORDPs, but they are prone to accept mutations with rates that are similar to those of IDPs. We conclude that IDPs are involved in the early stages of the speciation process, whereas mutations in ORDPs, once speciation is initiated, accelerate phenotypic divergence.

Molecular sex: The importance of base composition rather than homology when nucleic acids hybridize

On learning that nucleic acid hybridization had been achieved in a test tube, Huxley hailed the discovery of "molecular sex." The description was apt, since sex involves recombination, which requires hybridization that, in turn, depends on a successful homology search. Conversely, when the homology search fails, recombination fails. In yeast, this failure has been attributed to "simple sequence divergence." But sequence divergence does not impair nucleic acid hybridization simply. Most natural single-stranded nucleic acids are predisposed to adopt higher-order structures containing stem-loops. Tomizawa showed that the rate-limiting step in the hybridization of single-stranded sequences is an initial "kissing" exploration between complementary loops, which must first be appropriately extruded and aligned. Successful duplex formation requires successful synchronization of matching higher-ordered structures, which depends, not so much on the degree of similarity between their base sequences as on the closeness of their base compositions (GC%). In these terms, we can understand how the anti-recombinational effect of GC% differences supports the duplication both of genes within a genome and of genomes within a genus (speciation).

Evaluation of FORS-D analysis: a comparison with the statistically significant stem-loop potential

The stem-loop potential of a nucleic acid segment (expressed as a FONS value), decomposes into base composition-dependent and base order-dependent components. The latter, expressed as a FORS-D value, is derived by subtracting the value of the base composition-dependent component (FORS-M) from the FONS value. FORS-D analysis is the use of FORS-D values to estimate the potential of local base order to contribute to a stem-loop structure, and it has been used to investigate the relationship between stem-loop structure and other selective pressures on genomes. In the present study, we evaluated the reliability of FORS-D analysis by comparing it with statistically significant stem-loop potential, another robust method developed by Le and Maizel for examining stem-loop structure. We found that FORS-M values calculated using 10 randomized sequences are as reliable as those calculated using 100 randomized sequences. The resulting FORS-D values have a similar trend and distribution as statistically significant stem-loop potential, implying that FORS-D analysis is as reliable as the latter in measuring the distribution of base order-dependent stem-loop potential. Since the calculation of the FORS-M values is time consuming, the integrated program Bodslp developed by us will become a convenient tool for large-scale FORS-D analysis. The results also suggest that for some purposes the online program SigStb developed by Le and Maizel may be used as an alternative tool for FORS-D analysis.

Calculation of folding energies of single-stranded nucleic acid sequences: conceptual issues

The stability of a folded single-stranded nucleic acid depends on the composition and order of its constituent bases and may be assessed by taking into account the pairing energies of its constituent dinucleotides. To assess the possible biological significance of a computed structure, Maizel and coworkers in the 1980s compared the energy of folding of a natural single-stranded RNA sequence with the energies of several versions of the same sequence produced by shuffling base order. However, in the 2000s many took as self-evident the view that shuffling at the mononucleotide level (single bases) was conceptual wrong and should be replaced by shuffling at the level of dinucleotides (retaining pairs of adjacent bases). Folding energies then became indistinguishable from those of corresponding shuffled sequences and doubt was cast on the importance of secondary structures. Nevertheless, some continued productively to employ the single base shuffling approach, the justification for which is the topic of this paper. Because dinucleotide pairing energies are needed to calculate structure, it does not follow that shuffling should not disrupt dinucleotides. Base shuffling allows determination of the relative contributions of base composition and base order to total folding energy. The potential for secondary structure arises from pressures acting at both DNA and RNA levels, and is abundant throughout genomes-with a probable primary role in recombination. Within a gene the potential can often be accommodated, and base order and composition work together (values have the same negative sign) in contributing to total folding energy. But sometimes protein-coding pressure on base order conflicts with the pressure for secondary structure and the values have opposite signs. Total folding energy can be deemed of potential biological significance when the average of several readings is significantly less than zero.