Separation of microbial deoxyribonucleic acids into complementary strands

Chargaff's second parity rule lies at the origin of additive genetic interactions in quantitative traits to make omnigenic selection possible

Background

Francis Crick's central dogma provides a residue-by-residue mechanistic explanation of the flow of genetic information in living systems. However, this principle may not be sufficient for explaining how random mutations cause continuous variation of quantitative highly polygenic complex traits. Chargaff's second parity rule (CSPR), also referred to as intrastrand DNA symmetry, defined as near-exact equalities G ≈ C and A ≈ T within a single DNA strand, is a statistical property of cellular genomes. The phenomenon of intrastrand DNA symmetry was discovered more than 50 years ago; at present, it remains unclear what its biological role is, what the mechanisms are that force cellular genomes to comply strictly with CSPR, and why genomes of certain noncellular organisms have broken intrastrand DNA symmetry. The present work is aimed at studying a possible link between intrastrand DNA symmetry and the origin of genetic interactions in quantitative traits.

Methods

Computational analysis of single-nucleotide polymorphisms in human and mouse populations and of nucleotide composition biases at different codon positions in bacterial and human proteomes.

Results

The analysis of mutation spectra inferred from single-nucleotide polymorphisms observed in murine and human populations revealed near-exact equalities of numbers of reverse complementary mutations, indicating that random genetic variations obey CSPR. Furthermore, nucleotide compositions of coding sequences proved to be statistically interwoven via CSPR because pyrimidine bias at the 3rd codon position compensates purine bias at the 1st and 2nd positions.

Conclusions

According to Fisher's infinitesimal model, we propose that accumulation of reverse complementary mutations results in a continuous phenotypic variation due to small additive effects of statistically interwoven genetic variations. Therefore, additive genetic interactions can be inferred as a statistical entanglement of nucleotide compositions of separate genetic loci. CSPR challenges the neutral theory of molecular evolution-because all random mutations participate in variation of a trait-and provides an alternative solution to Haldane's dilemma by making a gene function diffuse. We propose that CSPR is symmetry of Fisher's infinitesimal model and that genetic information can be transferred in an implicit contactless manner.

The FtsK family of DNA translocases finds the ends of circles

A global view of bacterial chromosome choreography during the cell cycle is emerging, highlighting as a next challenge the description of the molecular mechanisms and factors involved. Here, we review one such factor, the FtsK family of DNA translocases. FtsK is a powerful and fast translocase that reads chromosome polarity. It couples segregation of the chromosome with cell division and controls the last steps of segregation in time and space. The second model protein of the family SpoIIIE acts in the transfer of the Bacillus subtilis chromosome during sporulation. This review focuses on the molecular mechanisms used by FtsK and SpoIIIE to segregate chromosomes with emphasis on the latest advances and open questions.

Inverse symmetry in complete genomes and whole-genome inverse duplication

The cause of symmetry is usually subtle, and its study often leads to a deeper understanding of the bearer of the symmetry. To gain insight into the dynamics driving the growth and evolution of genomes, we conducted a comprehensive study of textual symmetries in 786 complete chromosomes. We focused on symmetry based on our belief that, in spite of their extreme diversity, genomes must share common dynamical principles and mechanisms that drive their growth and evolution, and that the most robust footprints of such dynamics are symmetry related. We found that while complement and reverse symmetries are essentially absent in genomic sequences, inverse-complement plus reverse-symmetry is prevalent in complex patterns in most chromosomes, a vast majority of which have near maximum global inverse symmetry. We also discovered relations that can quantitatively account for the long observed but unexplained phenomenon of -mer skews in genomes. Our results suggest segmental and whole-genome inverse duplications are important mechanisms in genome growth and evolution, probably because they are efficient means by which the genome can exploit its double-stranded structure to enrich its code-inventory.

A study in entire chromosomes of violations of the intra-strand parity of complementary nucleotides (Chargaff's second parity rule)

Chargaff's rule of intra-strand parity (ISP) between complementary mono/oligonucleotides in chromosomes is well established in the scientific literature. Although a large numbers of papers have been published citing works and discussions on ISP in the genomic era, scientists are yet to find all the factors responsible for such a universal phenomenon in the chromosomes. In the present work, we have tried to address the issue from a new perspective, which is a parallel feature to ISP. The compositional abundance values of mono/oligonucleotides were determined in all non-overlapping sub-chromosomal regions of specific size. Also the frequency distributions of the mono/oligonucleotides among the regions were compared using the Kolmogorov–Smirnov test. Interestingly, the frequency distributions between the complementary mono/oligonucleotides revealed statistical similarity, which we named as intra-strand frequency distribution parity (ISFDP). ISFDP was observed as a general feature in chromosomes of bacteria, archaea and eukaryotes. Violation of ISFDP was also observed in several chromosomes. Chromosomes of different strains belonging a species in bacteria/archaea (Haemophilus influenza, Xylella fastidiosa etc.) and chromosomes of a eukaryote are found to be different among each other with respect to ISFDP violation. ISFDP correlates weakly with ISP in chromosomes suggesting that the latter one is not entirely responsible for the former. Asymmetry of replication topography and composition of forward-encoded sequences between the strands in chromosomes are found to be insufficient to explain the ISFDP feature in all chromosomes. This suggests that multiple factors in chromosomes are responsible for establishing ISFDP.

CodonExplorer: an online tool for analyzing codon usage and sequence composition, scaling from genes to genomes

DNA composition in general, and codon usage in particular, is crucial for understanding gene function and evolution. CodonExplorer, available online at http://bmf.colorado.edu/codonexplorer/, is an online tool and interactive database that contains millions of genes, allowing rapid exploration of the factors governing gene and genome compositional evolution and exploiting GC content and codon usage frequency to identify genes with composition suggesting high levels of expression or horizontal transfer.

Polarisation of prokaryotic chromosomes

In many prokaryotes, asymmetrical mutational or selective pressures have caused compositional skews between complementary strands of replication arms, especially sensitive in the distribution of guanine and cytosine. In Escherichia coli, most of the guanine/cytosine skew is caused by mutation rates differing on leading and lagging strands, but contribution of skewed functionally important guanine-rich motifs (Chi and Rag sites), which control chromosome repair or positioning, is noticeable. Interference between replication and gene expression plays a minor role. The situation may be different in other bacteria. Studies of chromosome processing and bacterial taxonomy might profit from consideration of chromosome polarisation.

DNA G+C content of the third codon position and codon usage biases of human genes

The human genome, as in other eukaryotes, has a wide heterogeneity in the DNA base composition. The evolutionary basis for this heterogeneity has been unknown. A previous study of the human genome (846 genes analyzed) has shown that, in the major range of the G+C content in the third codon position (0.25-0.75), biases from the Parity Rule 2 (PR2) among the synonymous codons of the four-codon amino acids are similar except in the highest G+C range (Sueoka, N., 1999. Translation-coupled violation of Parity Rule 2 in human genes is not the cause of heterogeneity of the DNA G+C content of third codon position. Gene 238, 53-58.). PR2 is an intra-strand rule where A=T and G=C are expected when there are no biases between the two complementary strands of DNA in mutation and selection rates (substitution rates). In this study, 14,026 human genes were analyzed. In addition, the third codon positions of two-codon amino acids were analyzed. New results show the following: (a) The G+C contents of the third codon position of human genes are scattered in the G+C range of 0.22-0.96 in the third codon position. (b) The PR2 biases are similar in the range of 0.25-0.75, whereas, in the high G+C range (0.75-0.96; 13% of the genes), the PR2-bias fingerprints are different from those of the major range. (c) Unlike the PR2 biases, the G+C contents of the third codon position for both four-codon and two-codon amino acids are all correlated almost perfectly with the G+C content of the third codon position over the total G+C ranges. These results support the notion that the directional mutation pressure, rather than the directional selection pressure, is mainly responsible for the heterogeneity of the G+C content of the third codon position.

Chargaff difference analysis of the bithorax complex of Drosophila melanogaster

Much of the fruit fly genome is compact ("Escherichia coli mode"), indicating a genome-wide selection pressure against DNA with little adaptive function. However, in the bithorax complex (BX-C) homeodomain genes are widely dispersed with large introns ("mammalian mode"). Chargaff difference analysis of compact bacterial and viral genomes has shown that most mRNAs have the potential to form stem-loop structures with purine-rich loops. Thus, for many taxa if transcription is to the right, the top (mRNA synonymous) DNA strand has purine-rich loop potential, and if transcription is to the left, the top (template) strand has pyrimidine-rich loop potential. The best indicator bases for transcription direction are A and T for AT-rich genomes, and C and G for CG-rich genomes. Consistent with this, Chargaff difference analysis of BX-C genes and several non-BX-C genes shows that, whatever the mode, mRNAs have the potential to form stem-loop structures with A-rich loops. We confirm that many potential open reading frames in the BX-C are unlikely to be functional. Conversely, we suggest that a few unassigned open reading frames may actually be functional. Since apparent organization in the mammalian mode cannot be explained in terms of unacknowledged open reading frames, yet the fruit fly genome is under pressure to be compact, it is likely that many BX-C functions do not involve the encoding of proteins.

Chromatographically fractionated complementary strands of Bacillus subtilis deoxyribonucleic acid: biological properties

Biological, physical, and chromatographic properties of methylated albuminkieselguhr (MAK)-fractionated complementary strands, designated as light (L) and heavy (H), of Bacillus subtilis deoxyribonucleic acid (DNA) are presented. The pattern of transforming activity along the MAK elution profile of alkilidenatured DNA shows that the residually active molecules selectively fractionated ahead of the L strand fraction, whereas the most active self-annealed molecules fractionated preferentially at the trailing end of the H strand fraction. The restoration rate of transforming activity in the late-eluting H molecules was rapid and independent of concentration during the annealing reaction. The data suggest that the self-annealing activity in the H strand is due in part to the formation of intrastrand secondary structures. Hydroxyapatite chromatography of self-annealed L and H strands yielded a major fraction (I) of highly purified strand preparations devoid of transforming activity and hypochromicity, and a minor "nativelike" fraction (II). Sedimentation velocity measurements show that, in addition to the mutual complementary nature of the L and H fractions, they differ in molecular size and possibly configuration.

Distribution of pyrimidine oligonucleotides in strands L and H of Bacillus subtilis DNA

The distribution of pyrimidine oligodeoxynucleotide clusters in L and H strands of Bacillus subtilis DNA separated by methylated albumin-Kieselguhr has been determined. Preparations of native and single-stranded DNA were degraded with diphenylamine in formic acid, and the released isostichs with the general formula of Py(n)P(n+1) were separated on DEAE-cellulose by chain length. Eleven isostichs were obtained for strands L and H in unequal proportions. Each isostich fraction was subfractionated by base composition on DEAE-cellulose at pH 3.0. 61 Pyrimidine oligonucleotide clusters were separated from the H strand and only 46 from the L strand. The findings show a higher degree of asymmetry between the strands in the distribution of cytosine-rich clusters as compared with thymine-rich clusters. The longest cytosine oligodeoxynucleotide present in both strands is of chain length 5. There is an unusually high distribution of thymine oligodeoxynucleotides of length 5-11. Up to chain length 6, the distribution of thymine oligodeoxynucleotides between the strands is about equal; from chain length 7 to 11 they occur predominantly in the H strand.

Action of DNA polymerase I of Escherichia coli with DNA-RNA hybrids as templates

Experiments indicating the ability of the ribo strand of a DNA-RNA template to guide polydeoxynucleotide synthesis by highly purified DNA polymerase I of E. coli (EC 2.7.7.7) are presented. With poly(rA).poly(dT) as template, poly(dT) is formed with a high efficiency, but almost no poly(dA). The specific activity of the enzyme, when tested with this template under suitable conditions, is eight times greater than that found for the poly(dA-dT) template. Single-stranded DNA fractions, with no template activity for DNA polymerase, are converted to efficient templates after their transcription by RNA polymerase. A concerted polymerization reaction, in which the action of DNA polymerase is dependent on that of RNA polymerase, can also be demonstrated with synthetic polydeoxynucleotides and single-stranded fractions of denatured DNA as templates.

Asymmetric template function of microbial deoxyribonucleic acids: transcription of messenger ribonucleic acid

In Bacillus subtilis and Escherichia coli, pulse-labeled ribonucleic acid (RNA) synthesized during step-down growth hybridized preferentially with the heavy (H) strand of methylated albumin-Kieselguhr-fractionated deoxyribonucleic acid (DNA). At high RNA inputs, the ratio of RNA hybridized with the H strand to that hybridized with the light (L) strand was 8.7 for B. subtilis and 2.0 for E. coli. At high DNA inputs, the H/L hybridization ratio increased by a factor of two. This change in the hybridization ratio was attributable to the fraction of the pulse-labeled RNA which is in stable RNA components. The hybridization peak of pulse-labeled RNA was specifically located in the late-eluting region of the absorbance profile of the H strand. This region was considered to represent the most actively transcribing H strand templates.

Asymmetric distribution of guanine plus thymine between complementary strands of deoxyribonucleic acid of members of the Enterobacteriaceae

Analysis of the deoxyribonucleic acid prepared from Proteus mirabilis, Escherichia coli, and Serratia marcescens in an alkaline CsCl gradient has shown that there is an asymmetric distribution of guanine plus thymine residues between the complementary strands of the deoxyribonucleic acid.

Asymmetric template function of microbial deoxyribonucleic acids: transcription of ribosomal and soluble ribonucleic acids

In Bacillus subtilis and Escherichia coli, 16 and 23S ribosomal ribonucleic acid (rRNA) hybridize exclusively with the heavy (H) strand of methylated albuminkieselguhr (MAK)-fractionated complementary deoxyribonucleic acid (DNA) strands. All the soluble RNA (4S RNA) in B. subtilis and 66 to 75% of the 4S RNA in E. coli also hybridize with the H strand. Interspecific hybridization shows that E. coli 23S rRNA also binds selectively to the DNA H strand of Salmonella typhimurium. The hybridization peak for all three cellular RNA components is specifically located in the late-eluting region of the absorbance profile of the DNA H strand. The early-eluting region of the light (L) strand preferentially inhibits the hybridization between the peak region of the H strand and 23S rRNA. These regions are considered to represent the transcribing sequences and their complements for 23S rRNA in the separated H and L strands of DNA, respectively.

Translational complementarity of RNA transcripts of native and denatured B. subtilis DNA

Previous evidence has suggested that the L and H fractions into which denatured DNA of B. subtilis can be separated by chromatography are complementary to each other with regard to base compositon, nucleotide sequence, and template properties, and that they may be regarded as families of fragments of the respective DNA strands. The present study tests the proposition that these fractions should also exhibit features of translational complementarity with respect to the DNA-dependent or the RNA-dependent incorporation of amino acids into ribosomes. This has been shown to be the case with the use of two pairs of amino acids: (1) lysine and phenylalanine; (2) proline and glycine.

Complementarity of RNA produced by enzymic transcription of native and denatured B. subtilis DNA

This paper describes the preparation and some of the properties of the RNA specimens synthesized with the aid of the RNA polymerase of E. coli by transcription of the following DNA templates: (a) undenatured B. subtilis DNA (yielding N-RNA); (b) separated strand fractions L and H isolated by chromatography of the denatured DNA on methylated albumin (yielding L-RNA and H-RNA, respectively). The study of the hybridization behavior of the various RNA products showed that N-RNA, though able to form hybrids with either strand, hybridized with H-DNA to twice as great an extent as with L-DNA. The transcripts of the separated L and H fractions exhibited complete specificity with respect to complexing: L-RNA hybridized only with L-DNA, H-RNA only with H-DNA.

Template properties of complementary fractions of denatured microbial deoxyribonucleic acids

DNA preparations from seven bacterial species and from E. coli phage T4, and also the complementary L and H fractions into which these DNA specimens, after denaturation, were separated by chromatography on methylated albumin-kieselguhr columns, were studied as templates in the RNA polymerase system, and the nucleotide composition of the RNA products was determined. The RNA transcripts of the separated L and H fractions were found to be faithful copies of the respective DNA fractions. This suggests "transcription analysis" as a sensitive and reliable analytical technique for the determination of the base composition of denatured DNA. The L and H fractions of T4 DNA were shown, both by temperature-absorbance profiles and by transcription analysis, to be mutually complementary. The RNA products formed with intact DNA as the template were not exact copies of the latter; their composition indicated that under our experimental conditions the "heavy" DNA strand is transcribed preferentially.

Complementary fractions of denatured DNA of coliphage T3 as templates for transcription

This paper discusses the evidence that two fractions obtained by the chromatography of denatured DNA of phage T3 on columns of methylated albumin-kieselguhr represent the complementary strands. This evidence derives from a variety of temperature-absorbance measurements and from the base composition of the RNA products synthesized by RNA polymerase under the direction of the two single-stranded DNA templates.