The tendency to recreate ancestral CG dinucleotides in the human genome

Background  

The CG dinucleotides are known to be deficient in the human genome, due to a high mutation rate from 5-methylated CG to TG and its complementary pair CA. Meanwhile, many cellular functions rely on these CG dinucleotides, such as gene expression controlled by cytosine methylation status. Thus, CG dinucleotides that provide essential functional substrates should be retained in genomes. How these two conflicting processes regarding the fate of CG dinucleotides - i.e., high mutation rate destroying CG dinucleotides, vs. functional processes that require their preservation remains an unsolved question.     

Structure–activity relationships of dinucleotides: Potent and selective agonists of P2Y receptors

Dinucleoside polyphosphates act as agonists on purinergic P2Y receptors to mediate a variety of cellular processes. Symmetrical, naturally occurring purine dinucleotides are found in most living cells and their actions are generally known. Unsymmetrical purine dinucleotides and all pyrimidine containing dinucleotides, however, are not as common and therefore their actions are not well understood. To carry out a thorough examination of the activities and specificities of these dinucleotides, a robust method of synthesis was developed to allow manipulation of either nucleoside of the dinucleotide as well as the phosphate chain lengths. Adenosine containing dinucleotides exhibit some level of activity on P2Y₁ while uridine containing dinucleotides have some level of agonist response on P2Y₂ and P2Y₆. The length of the linking phosphate chain determines a different specificity; diphosphates are most accurately mimicked by dinucleoside triphosphates and triphosphates most resemble dinucleoside tetraphosphates. The pharmacological activities and relative metabolic stabilities of these dinucleotides are reported with their potential therapeutic applications being discussed.     

CpG Dinucleotide Frequencies Reveal the Role of Host Methylation Capabilities in Parvovirus Evolution

Parvoviruses are rapidly evolving viruses that infect a wide range of hosts, including vertebrates and invertebrates. Extensive methylation of the parvovirus genome has been recently demonstrated. A global pattern of methylation of CpG dinucleotides is seen in vertebrate genomes, compared to “fractional” methylation patterns in invertebrate genomes. It remains unknown if the loss of CpG dinucleotides occurs in all viruses of a given DNA virus family that infect host species spanning across vertebrates and invertebrates. We investigated the link between the extent of CpG dinucleotide depletion among autonomous parvoviruses and the evolutionary lineage of the infected host. We demonstrate major differences in the relative abundance of CpG dinucleotides among autonomous parvoviruses which share similar genome organization and common ancestry, depending on the infected host species. Parvoviruses infecting vertebrate hosts had significantly lower relative abundance of CpG dinucleotides than parvoviruses infecting invertebrate hosts. The strong correlation of CpG dinucleotide depletion with the gain in TpG/CpA dinucleotides and the loss of TpA dinucleotides among parvoviruses suggests a major role for CpG methylation in the evolution of parvoviruses. Our data present evidence that links the relative abundance of CpG dinucleotides in parvoviruses to the methylation capabilities of the infected host. In sum, our findings support a novel perspective of host-driven evolution among autonomous parvoviruses.     

A comprehensive survey of human polymorphisms at conserved splice dinucleotides and its evolutionary relationship with alternative splicing

Background  

Alternative splicing (AS) is a key molecular process that endows biological functions with diversity and complexity. Generally, functional redundancy leads to the generation of new functions through relaxation of selective pressure in evolution, as exemplified by duplicated genes. It is also known that alternatively spliced exons (ASEs) are subject to relaxed selective pressure. Within consensus sequences at the splice junctions, the most conserved sites are dinucleotides at both ends of introns (splice dinucleotides). However, a small number of single nucleotide polymorphisms (SNPs) occur at splice dinucleotides. An intriguing question relating to the evolution of AS diversity is whether mutations at splice dinucleotides are maintained as polymorphisms and produce diversity in splice patterns within the human population. We therefore surveyed validated SNPs in the database dbSNP located at splice dinucleotides of all human genes that are defined by the H-Invitational Database.     

Sequence structure of human nucleosome DNA

Positional distributions of various dinucleotides in experimentally derived human nucleosome DNA sequences are analyzed. Nucleosome positioning in this species is found to depend largely on GG and CC dinucleotides periodically distributed along the nucleosome DNA sequence, with the period of 10.4 bases. The GG and CC dinucleotides oscillate counterphase, i.e., their respective preferred positions are shifted about a half-period from one another, as it was observed earlier for AA and TT dinucleotides. Other purine-purine and pyrimidine-pyrimidine dinucleotides (RR and YY) display the same periodical and counterphase pattern. The dominance of oscillating GG and CC dinucleotides in human nucleosomes and the contribution of AG(CT), GA(TC), and AA(TT) suggest a general nucleosome DNA sequence pattern - counterphase oscillation of RR and YY dinucleotides. AA and TT dinucleotides, commonly accepted as major players, are only weak contributors in the case of human nucleosomes.     

Sequence Structure of Hidden 10.4-base Repeat in the Nucleosomes of C. elegans

Abstract

By measuring prevailing distances between YY, YR, RR, and RY dinucleotides in the large database of the nucleosome DNA fragments from C. elegans, the consensus sequence structure of the nucleosome DNA repeat of C. elegans was reconstructed: (YYYYYRRRRR)_n. An actual period was estimated to be 10.4 bases. The pattern is fully consistent with the nucleosome DNA patterns of other eukaryotes, as established earlier, and, thus, the YYYYYRRRRR repeat can be considered as consensus nucleosome DNA sequence repeat across eukaryotic species. Similar distance analysis for [A, T] dinucleotides suggested the related pattern (TTTYTARAAA)_n where the TT and AA dinucleotides display rather out of phase behavior, contrary to the “AA or TT” in-phase periodicity, considered in some publications. A weak 5-base periodicity in the distribution of TA dinucleotides was detected.     

Sequence-specific ultrasonic cleavage of DNA

We investigated the phenomenon of ultrasonic cleavage of DNA by analyzing a large set of cleavage patterns of DNA restriction fragments using polyacrylamide gel electrophoresis. The cleavage intensity of individual phosphodiester bonds was found to depend on the nucleotide sequence and the position of the bond with respect to the ends of the fragment. The relative intensities of cleavage of the central phosphodiester bond in 16 dinucleotides and 256 tetranucleotides were determined by multivariate statistical analysis. We observed a remarkable enhancement of the mean values of the relative intensities of cleavage (cleavage rates) in phosphodiester bonds following deoxycytidine, which diminished in the row of dinucleotides: d(CpG) > d(CpA) > d(CpT) >> d(CpC). The cleavage rates for all pairs of complementary dinucleotides were significantly different from each other. The effect of flanking nucleotides in tetranucleotides on cleavage rates of all 16 types of central dinucleotides was also statistically significant. The sequence-dependent ultrasonic cleavage rates of dinucleotides are consistent with reported data on the intensity of the conformational motion of their 5′-deoxyribose. As a measure of local conformational dynamics, cleavage rates may be useful for characterizing functional regions of the genome.     

REMARKABLE SELECTIVE CONSTRAINTS ON EXONIC DINUCLEOTIDE REPEATS

Long dinucleotide repeats found in exons present a substantial mutational hazard: mutations at these loci occur often and generate frameshifts. Here, we provide clear and compelling evidence that exonic dinucleotides experience strong selective constraint. In humans, only 18 exonic dinucleotides have repeat lengths greater than six, which contrasts sharply with the genome‐wide distribution of dinucleotides. We genotyped each of these dinucleotides in 200 humans from eight 1000 Genomes Project populations and found a near‐absence of polymorphism. More remarkably, divergence data demonstrate that repeat lengths have been conserved across the primate phylogeny in spite of what is likely considerable mutational pressure. Coalescent simulations show that even a very low mutation rate at these loci fails to explain the anomalous patterns of polymorphism and divergence. Our data support two related selective constraints on the evolution of exonic dinucleotides: a short‐term intolerance for any change to repeat length and a long‐term prevention of increases to repeat length. In general, our results implicate purifying selection as the force that eliminates new, deleterious mutants at exonic dinucleotides. We briefly discuss the evolution of the longest exonic dinucleotide in the human genome—a 10 x CA repeat in fibroblast growth factor receptor‐like 1 (FGFRL1)—that should possess a considerably greater mutation rate than any other exonic dinucleotide and therefore generate a large number of deleterious variants.     

Systematic CpT (ApG) Depletion and CpG Excess Are Unique Genomic Signatures of Large DNA Viruses Infecting Invertebrates

Differences in the relative abundance of dinucleotides, if any may provide important clues on host-driven evolution of viruses. We studied dinucleotide frequencies of large DNA viruses infecting vertebrates (n = 105; viruses infecting mammals = 99; viruses infecting aves = 6; viruses infecting reptiles = 1) and invertebrates (n = 88; viruses infecting insects = 84; viruses infecting crustaceans = 4). We have identified systematic depletion of CpT(ApG) dinucleotides and over-representation of CpG dinucleotides as the unique genomic signature of large DNA viruses infecting invertebrates. Detailed investigation of this unique genomic signature suggests the existence of invertebrate host-induced pressures specifically targeting CpT(ApG) and CpG dinucleotides. The depletion of CpT dinucleotides among large DNA viruses infecting invertebrates is at least in part, explained by non-canonical DNA methylation by the infected host. Our findings highlight the role of invertebrate host-related factors in shaping virus evolution and they also provide the necessary framework for future studies on evolution, epigenetics and molecular biology of viruses infecting this group of hosts.     

Repertoires of the Nucleosome-Positioning Dinucleotides

It is generally accepted that the organization of eukaryotic DNA into chromatin is strongly governed by a code inherent in the genomic DNA sequence. This code, as well as other codes, is superposed on the triplets coding for amino acids. The history of the chromatin code started three decades ago with the discovery of the periodic appearance of certain dinucleotides, with AA/TT and RR/YY giving the strongest signals, all with a period of 10.4 bases. Every base-pair stack in the DNA duplex has specific deformation properties, thus favoring DNA bending in a specific direction. The appearance of the corresponding dinucleotide at the distance 10.4 xn bases will facilitate DNA bending in that direction, which corresponds to the minimum energy of DNA folding in the nucleosome. We have analyzed the periodic appearances of all 16 dinucleotides in the genomes of thirteen different eukaryotic organisms. Our data show that a large variety of dinucleotides (if not all) are, apparently, contributing to the nucleosome positioning code. The choice of the periodical dinucleotides differs considerably from one organism to another. Among other 10.4 base periodicities, a strong and very regular 10.4 base signal was observed for CG dinucleotides in the genome of the honey bee A. mellifera. Also, the dinucleotide CG appears as the only periodical component in the human genome. This observation seems especially relevant since CpG methylation is well known to modulate chromatin packing and regularity. Thus, the selection of the dinucleotides contributing to the chromatin code is species specific, and may differ from region to region, depending on the sequence context.     

Human nucleosomes: special role of CG dinucleotides and Alu-nucleosomes

Background

The periodical occurrence of dinucleotides with a period of 10.4 bases now is undeniably a hallmark of nucleosome positioning. Whereas many eukaryotic genomes contain visible and even strong signals for periodic distribution of dinucleotides, the human genome is rather featureless in this respect. The exact sequence features in the human genome that govern the nucleosome positioning remain largely unknown.

Results

When analyzing the human genome sequence with the positional autocorrelation method, we found that only the dinucleotide CG shows the 10.4 base periodicity, which is indicative of the presence of nucleosomes. There is a high occurrence of CG dinucleotides that are either 31 (10.4 × 3) or 62 (10.4 × 6) base pairs apart from one another - a sequence bias known to be characteristic of Alu-sequences. In a similar analysis with repetitive sequences removed, peaks of repeating CG motifs can be seen at positions 10, 21 and 31, the nearest integers of multiples of 10.4.

Conclusions

Although the CG dinucleotides are dominant, other elements of the standard nucleosome positioning pattern are present in the human genome as well. The positional autocorrelation analysis of the human genome demonstrates that the CG dinucleotide is, indeed, one visible element of the human nucleosome positioning pattern, which appears both in Alu sequences and in sequences without repeats. The dominant role that CG dinucleotides play in organizing human chromatin is to indicate the involvement of human nucleosomes in tuning the regulation of gene expression and chromatin structure, which is very likely due to cytosine-methylation/-demethylation in CG dinucleotides contained in the human nucleosomes. This is further confirmed by the positions of CG-periodical nucleosomes on Alu sequences. Alu repeats appear as monomers, dimers and trimers, harboring two to six nucleosomes in a run. Considering the exceptional role CG dinucleotides play in the nucleosome positioning, we hypothesize that Alu-nucleosomes, especially, those that form tightly positioned runs, could serve as "anchors" in organizing the chromatin in human cells.     

Frequency of abnormal human haemoglobins caused by C----T transitions in CpG dinucleotides

A large part of human genetic disease apparently arises from deamination of cytosine residues in methylated CpG dinucleotides. Their mutation rate is known to be high when C is present as 5-methyl-cytosine, but is believed to be normal when it is unmethylated. The beta-globin gene contains five, the gamma-globin gene two, and each of the alpha-globin genes contains 35 CpG dinucleotides. The CpG dinucleotides in the beta and gamma-globin genes are methylated, while those in the alpha-globin genes are under-methylated. One would therefore have expected the CpG dinucleotides to be a frequent source of mutations in the beta and gamma-globin genes, but not in the alpha-globin genes. In fact, the evidence points to CpG dinucleotides being a frequent source of mutations in both the alpha and beta-globin genes. This suggests either that the mutation rates of both methylated and unmethylated CpG dinucleotides are abnormally high, which conflicts with published evidence, or that there is a finite chance of some of these in the alpha-globin genes of certain individuals being methylated and therefore subject to mutation.     

The Universal Trend of Amino Acid Gain–Loss is Caused by CpG Hypermutability

Understanding the cause of the changes in the amino acid composition of proteins is essential for understanding the evolution of protein functions. Since the early 1970s, it has been known that the frequency of some amino acids in protein sequences is increasing and that of others is decreasing. Recently, it was found that the trends of amino acid changes were similar in 15 taxa representing Bacteria, Archaea, and Eukaryota. However, the cause of this similarity in the trend of the gains and losses of amino acids continued to be debated. Here, we show that this trend of the gain and loss of amino acids can be simply explained by CpG hypermutability. We found that the frequency of amino acids coded by codons with TpG dinucleotides and those with CpA dinucleotides is increasing, while that of amino acids coded by codons with CpG dinucleotides is decreasing. We also found that organisms that lack DNA methyltransferase show different trends of the gain and loss of amino acids. DNA methyltransferase methylates CpG dinucleotides and induces CpG hypermutability. The incorporation of CpG hypermutability into models of protein evolution will improve studies on protein evolution in different organisms. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

PCR-based evaluation of sequence specificity of DNA fragmentation by ultrasound

Ultrasonic fragmentation, which is a simple and convenient method for the mechanical degradation of DNA, is widely used in modern genome studies as one of the sample preparation steps. It has been recently found that the DNA breaks occur more often in the regions containing 5'-CG-3' dinucleotides. We studied the influence of the 5'-CG-3' dinucleotides on the efficiency of the 28S rRNA gene amplification during PCR with sonicated DNA of Mantis religiosa. It was shown that the amplification rate depends on the template length and the number of 5'-CG-3' dinucleotides. Amplification of the DNA regions with a higher 5'-CG-3' density is less efficient because of their higher sensitivity to ultrasound. The amount of the amplified DNA templates is inversely proportional to the 5'-CG-3'number.     

Human Monocyte Recognition of Adenosine-Based Cyclic Dinucleotides Unveils the A2a Gαs Protein-Coupled Receptor Tonic Inhibition of Mitochondrially Induced Cell Death

Cyclic dinucleotides are important messengers for bacteria and protozoa and are well-characterized immunity alarmins for infected mammalian cells through intracellular binding to STING receptors. We sought to investigate their unknown extracellular effects by adding cyclic dinucleotides to the culture medium of freshly isolated human blood cells in vitro. Here we report that adenosine-containing cyclic dinucleotides induce the selective apoptosis of monocytes through a novel apoptotic pathway. We demonstrate that these compounds are inverse agonist ligands of A2a, a G_αs-coupled adenosine receptor selectively expressed by monocytes. Inhibition of monocyte A2a by these ligands induces apoptosis through a mechanism independent of that of the STING receptors. The blockade of basal (adenosine-free) signaling from A2a inhibits protein kinase A (PKA) activity, thereby recruiting cytosolic p53, which opens the mitochondrial permeability transition pore and impairs mitochondrial respiration, resulting in apoptosis. A2a antagonists and inverse agonist ligands induce apoptosis of human monocytes, while A2a agonists are antiapoptotic. In vivo, we used a mock developing human hematopoietic system through NSG mice transplanted with human CD34⁺ cells. Treatment with cyclic di-AMP selectively depleted A2a-expressing monocytes and their precursors via apoptosis. Thus, monocyte recognition of cyclic dinucleotides unravels a novel proapoptotic pathway: the A2a G_αs protein-coupled receptor (GPCR)-driven tonic inhibitory signaling of mitochondrion-induced cell death.     

Cyclic‐di‐GMP and cyclic‐di‐AMP activate the NLRP3 inflammasome

The cyclic dinucleotides 3'‐5'diadenylate (c‐diAMP) and 3'‐5' diguanylate (c‐diGMP) are important bacterial second messengers that have recently been shown to stimulate the secretion of type I Interferons (IFN‐Is) through the c‐diGMP‐binding protein MPYS/STING. Here, we show that physiologically relevant levels of cyclic dinucleotides also stimulate a robust secretion of IL‐1β through the NLRP3 inflammasome. Intriguingly, this response is independent of MPYS/STING. Consistent with most NLRP3 inflammasome activators, the response to c‐diGMP is dependent on the mobilization of potassium and calcium ions. However, in contrast to other NLRP3 inflammasome activators, this response is not associated with significant changes in mitochondrial potential or the generation of mitochondrial reactive oxygen species. Thus, cyclic dinucleotides activate the NLRP3 inflammasome through a unique pathway that could have evolved to detect pervasive bacterial pathogen‐associated molecular patterns associated with intracellular infections.