CpG Suppression Research Articles

Alphavirus Evasion of Zinc Finger Antiviral Protein (ZAP) Correlates with CpG Suppression in a Specific Viral nsP2 Gene Sequence.

Certain re-emerging alphaviruses, such as chikungunya virus (CHIKV), cause serious disease and widespread epidemics. To develop virus-specific therapies, it is critical to understand the determinants of alphavirus pathogenesis and virulence. One major determinant is viral evasion of the host interferon response, which upregulates antiviral effectors, including zinc finger antiviral protein (ZAP). Here, we demonstrated that Old World alphaviruses show differential sensitivity to endogenous ZAP in 293T cells: Ross River virus (RRV) and Sindbis virus (SINV) are more sensitive to ZAP than o'nyong'nyong virus (ONNV) and CHIKV. We hypothesized that the more ZAP-resistant alphaviruses evade ZAP binding to their RNA. However, we did not find a correlation between ZAP sensitivity and binding to alphavirus genomic RNA. Using a chimeric virus, we found the ZAP sensitivity determinant lies mainly within the alphavirus non-structural protein (nsP) gene region. Surprisingly, we also did not find a correlation between alphavirus ZAP sensitivity and binding to nsP RNA, suggesting ZAP targeting of specific regions in the nsP RNA. Since ZAP can preferentially bind CpG dinucleotides in viral RNA, we identified three 500-bp sequences in the nsP region where CpG content correlates with ZAP sensitivity. Interestingly, ZAP binding to one of these sequences in the nsP2 gene correlated to sensitivity, and we confirmed that this binding is CpG-dependent. Our results demonstrate a potential strategy of alphavirus virulence by localized CpG suppression to evade ZAP recognition.

Less is more: Biased loss of CpG dinucleotides strengthens antiviral immunity.

Viruses may not only affect our daily lives but also shape our genome evolution. A recent study shows that the zinc-finger antiviral protein (ZAP) drives CpG suppression in a biased manner. Genes involved in the defense against viral invaders are particularly CpG suppressed to avoid self-targeting and to promote an effective immune response.

Minimal impact of ZAP on lentiviral vector production and transduction efficiency

The antiviral protein ZAP binds CpG dinucleotides in viral RNA to inhibit replication. This has likely led to the CpG suppression observed in many RNA viruses, including retroviruses. Sequences added to retroviral vector genomes, such as internal promoters, transgenes, or regulatory elements, substantially increase CpG abundance. Because these CpGs could allow retroviral vector RNA to be targeted by ZAP, we analyzed whether it restricts vector production, transduction efficiency, and transgene expression. Surprisingly, even though CpG-high HIV-1 was efficiently inhibited by ZAP in HEK293T cells, depleting ZAP did not substantially increase lentiviral vector titer using several packaging and genome plasmids. ZAP overexpression also did not inhibit lentiviral vector titer. In addition, decreasing CpG abundance in a lentiviral vector genome did not increase its titer, and a gammaretroviral vector derived from murine leukemia virus was not substantially restricted by ZAP. Overall, we show that the increased CpG abundance in retroviral vectors relative to the wild-type retroviruses they are derived from does not intrinsically sensitize them to ZAP. Further understanding of how ZAP specifically targets transcripts to inhibit their expression may allow the development of CpG sequence contexts that efficiently recruit or evade this antiviral system.

Origin and evolution of the zinc finger antiviral protein.

The human zinc finger antiviral protein (ZAP) recognizes RNA by binding to CpG dinucleotides. Mammalian transcriptomes are CpG-poor, and ZAP may have evolved to exploit this feature to specifically target non-self viral RNA. Phylogenetic analyses reveal that ZAP and its paralogue PARP12 share an ancestral gene that arose prior to extensive eukaryote divergence, and the ZAP lineage diverged from the PARP12 lineage in tetrapods. Notably, the CpG content of modern eukaryote genomes varies widely, and ZAP-like genes arose subsequent to the emergence of CpG-suppression in vertebrates. Human PARP12 exhibited no antiviral activity against wild type and CpG-enriched HIV-1, but ZAP proteins from several tetrapods had antiviral activity when expressed in human cells. In some cases, ZAP antiviral activity required a TRIM25 protein from the same or related species, suggesting functional co-evolution of these genes. Indeed, a hypervariable sequence in the N-terminal domain of ZAP contributed to species-specific TRIM25 dependence in antiviral activity assays. Crosslinking immunoprecipitation coupled with RNA sequencing revealed that ZAP proteins from human, mouse, bat and alligator exhibit a high degree of CpG-specificity, while some avian ZAP proteins appear more promiscuous. Together, these data suggest that the CpG- rich RNA directed antiviral activity of ZAP-related proteins arose in tetrapods, subsequent to the onset of CpG suppression in certain eukaryote lineages, with subsequent species-specific adaptation of cofactor requirements and RNA target specificity.

The dinucleotide composition of the Zika virus genome is shaped by conflicting evolutionary pressures in mammalian hosts and mosquito vectors.

Most vertebrate RNA viruses show pervasive suppression of CpG and UpA dinucleotides, closely resembling the dinucleotide composition of host cell transcriptomes. In contrast, CpG suppression is absent in both invertebrate mRNA and RNA viruses that exclusively infect arthropods. Arthropod-borne (arbo) viruses are transmitted between vertebrate hosts by invertebrate vectors and thus encounter potentially conflicting evolutionary pressures in the different cytoplasmic environments. Using a newly developed Zika virus (ZIKV) model, we have investigated how demands for CpG suppression in vertebrate cells can be reconciled with potentially quite different compositional requirements in invertebrates and how this affects ZIKV replication and transmission. Mutant viruses with synonymously elevated CpG or UpA dinucleotide frequencies showed attenuated replication in vertebrate cell lines, which was rescued by knockout of the zinc-finger antiviral protein (ZAP). Conversely, in mosquito cells, ZIKV mutants with elevated CpG dinucleotide frequencies showed substantially enhanced replication compared to wild type. Host-driven effects on virus replication attenuation and enhancement were even more apparent in mouse and mosquito models. Infections with CpG- or UpA-high ZIKV mutants in mice did not cause typical ZIKV-induced tissue damage and completely protected mice during subsequent challenge with wild-type virus, which demonstrates their potential as live-attenuated vaccines. In contrast, the CpG-high mutants displayed enhanced replication in Aedes aegypti mosquitoes and a larger proportion of mosquitoes carried infectious virus in their saliva. These findings show that mosquito cells are also capable of discriminating RNA based on dinucleotide composition. However, the evolutionary pressure on the CpG dinucleotides of viral genomes in arthropod vectors directly opposes the pressure present in vertebrate host cells, which provides evidence that an adaptive compromise is required for arbovirus transmission. This suggests that the genome composition of arbo flaviviruses is crucial to maintain the balance between high-level replication in the vertebrate host and persistent replication in the mosquito vector.

Use of a small DNA virus model to investigate mechanisms of CpG dinucleotide-induced attenuation of virus replication.

Suppression of the CpG dinucleotide is widespread in RNA viruses infecting vertebrates and plants, and in the genomes of retroviruses and small mammalian DNA viruses. The functional basis for CpG suppression in the latter was investigated through the construction of mutants of the parvovirus, minute virus of mice (MVM) with increased CpG or TpA dinucleotides in the VP gene. CpG-high mutants displayed extraordinary attenuation in A9 cells compared to wild-type MVM (>six logs), while TpA elevation showed no replication effect. Attenuation was independent of Toll-like receptor 9 and STING-mediated DNA recognition pathways and unrelated to effects on translation efficiency. While translation from codon-optimized VP RNA was enhanced in a cell-free assay, MVM containing this sequence was highly attenuated. Further mutational analysis indicated that this arose through its increased numbers of CpG dinucleotides (7→70) and separately from its increased G+C content (42.3→57.4 %), which independently attenuated replication. CpG-high viruses showed impaired NS mRNA expression by qPCR and reduced NS and particularly VP protein expression detected by immunofluorescence and replication in A549 cells, effects reversed in zinc antiviral protein (ZAP) knockout cells, even though nuclear relocalization of VP remained defective. The demonstrated functional basis for CpG suppression in MVM and potentially other small DNA viruses and the observed intolerance of CpGs in coding sequences, even after codon optimization, has implications for the use of small DNA virus vectors in gene therapy and immunization.

Intra-genome variability in the dinucleotide composition of SARS-CoV-2.

CpG dinucleotides are under-represented in the genomes of single-stranded RNA viruses, and SARS-CoV-2 is no exception to this. Artificial modification of CpG frequency is a valid approach for live attenuated vaccine development; if this is to be applied to SARS-CoV-2, we must first understand the role CpG motifs play in regulating SARS-CoV-2 replication. Accordingly, the CpG composition of the SARS-CoV-2 genome was characterised. CpG suppression among coronaviruses does not differ between virus genera but does vary with host species and primary replication site (a proxy for tissue tropism), supporting the hypothesis that viral CpG content may influence cross-species transmission. Although SARS-CoV-2 exhibits overall strong CpG suppression, this varies considerably across the genome, and the Envelope (E) open reading frame (ORF) and ORF10 demonstrate an absence of CpG suppression. Across the Coronaviridae, E genes display remarkably high variation in CpG composition, with those of SARS and SARS-CoV-2 having much higher CpG content than other coronaviruses isolated from humans. This is an ancestrally derived trait reflecting their bat origins. Conservation of CpG motifs in these regions suggests that they have a functionality which over-rides the need to suppress CpG; an observation relevant to future strategies towards a rationally attenuated SARS-CoV-2 vaccine.

Synonymous Dinucleotide Usage: A Codon-Aware Metric for Quantifying Dinucleotide Representation in Viruses.

Distinct patterns of dinucleotide representation, such as CpG and UpA suppression, are characteristic of certain viral genomes. Recent research has uncovered vertebrate immune mechanisms that select against specific dinucleotides in targeted viruses. This evidence highlights the importance of systematically examining the dinucleotide composition of viral genomes. We have developed a novel metric, called synonymous dinucleotide usage (SDU), for quantifying dinucleotide representation in coding sequences. Our method compares the abundance of a given dinucleotide to the null hypothesis of equal synonymous codon usage in the sequence. We present a Python3 package, DinuQ, for calculating SDU and other relevant metrics. We have applied this method on two sets of invertebrate- and vertebrate-specific flaviviruses and rhabdoviruses. The SDU shows that the vertebrate viruses exhibit consistently greater under-representation of CpG dinucleotides in all three codon positions in both datasets. In comparison to existing metrics for dinucleotide quantification, the SDU allows for a statistical interpretation of its values by comparing it to a null expectation based on the codon table. Here we apply the method to viruses, but coding sequences of other living organisms can be analysed in the same way.

CpG Dinucleotides Inhibit HIV-1 Replication through Zinc Finger Antiviral Protein (ZAP)-Dependent and -Independent Mechanisms.

CpG dinucleotides are suppressed in the genomes of many vertebrate RNA viruses, including HIV-1. The cellular antiviral protein ZAP (zinc finger antiviral protein) binds CpGs and inhibits HIV-1 replication when CpGs are introduced into the viral genome. However, it is not known if ZAP-mediated restriction is the only mechanism driving CpG suppression. To determine how CpG dinucleotides affect HIV-1 replication, we increased their abundance in multiple regions of the viral genome and analyzed the effect on RNA expression, protein abundance, and infectious-virus production. We found that the antiviral effect of CpGs was not correlated with their abundance. Interestingly, CpGs inserted into some regions of the genome sensitize the virus to ZAP antiviral activity more efficiently than insertions into other regions, and this sensitivity can be modulated by interferon treatment or ZAP overexpression. Furthermore, the sensitivity of the virus to endogenous ZAP was correlated with its sensitivity to the ZAP cofactor KHNYN. Finally, we show that CpGs in some contexts can also inhibit HIV-1 replication by ZAP-independent mechanisms, and one of these is the activation of a cryptic splice site at the expense of a canonical splice site. Overall, we show that the location and sequence context of the CpG in the viral genome determines its antiviral activity.IMPORTANCE Some RNA virus genomes are suppressed in the nucleotide combination of a cytosine followed by a guanosine (CpG), indicating that they are detrimental to the virus. The antiviral protein ZAP binds viral RNA containing CpGs and prevents the virus from multiplying. However, it remains unknown how the number and position of CpGs in viral genomes affect restriction by ZAP and whether CpGs have other antiviral mechanisms. Importantly, manipulating the CpG content in viral genomes could help create new vaccines. HIV-1 shows marked CpG suppression, and by introducing CpGs into its genome, we show that ZAP efficiently targets a specific region of the viral genome, that the number of CpGs does not predict the magnitude of antiviral activity, and that CpGs can inhibit HIV-1 gene expression through a ZAP-independent mechanism. Overall, the position of CpGs in the HIV-1 genome determines the magnitude and mechanism through which they inhibit the virus.

CpG Frequency in the 5' Third of the env Gene Determines Sensitivity of Primary HIV-1 Strains to the Zinc-Finger Antiviral Protein.

CpG dinucleotide suppression has been reported to allow HIV-1 to evade inhibition by the zinc-finger antiviral protein (ZAP). Here, we show that primate lentiviruses display marked differences in CpG frequencies across their genome, ranging from 0.44% in simian immunodeficiency virus SIVwrc from Western red colobus to 2.3% in SIVmon infecting mona monkeys. Moreover, functional analyses of a large panel of human and simian immunodeficiency viruses revealed that the magnitude of CpG suppression does not correlate with their susceptibility to ZAP. However, we found that the number of CpG dinucleotides within a region of ∼700 bases at the 5' end of the env gene determines ZAP sensitivity of primary HIV-1 strains but not of HIV-2. Increased numbers of CpGs in this region were associated with reduced env mRNA expression and viral protein production. ZAP sensitivity profiles of chimeric simian-human immunodeficiency viruses (SHIVs) expressing different HIV-1 env genes were highly similar to those of the corresponding HIV-1 strains. The frequency of CpGs in the identified env region correlated with differences in clinical progression rates. Thus, the CpG frequency in a specific part of env, rather than the overall genomic CpG content, governs the susceptibility of HIV-1 to ZAP and might affect viral pathogenicity in vivoIMPORTANCE Evasion of the zinc-finger antiviral protein (ZAP) may drive CpG dinucleotide suppression in HIV-1 and many other viral pathogens but the viral determinants of ZAP sensitivity are poorly defined. Here, we examined CpG suppression and ZAP sensitivity in a large number of primate lentiviruses and demonstrate that their genomic frequency of CpGs varies substantially and does not correlate with ZAP sensitivity. We further show that the number of CpG residues in a defined region at the 5' end of the env gene together with structural features plays a key role in HIV-1 susceptibility to ZAP and correlates with differences in clinical progression rates in HIV-1-infected individuals. Our identification of a specific part of env as a major determinant of HIV-1 susceptibility to ZAP restriction provides a basis for future studies of the underlying inhibitory mechanisms and their potential relevance in the pathogenesis of AIDS.

Knockdown of Plakophilin 2 Downregulates miR-184 Through CpG Hypermethylation and Suppression of the E2F1 Pathway and Leads to Enhanced Adipogenesis In Vitro.

PKP2, encoding plakophilin 2 (PKP2), is the most common causal gene for arrhythmogenic cardiomyopathy. To characterize miRNA expression profile in PKP2-deficient cells. Control and PKP2-knockdown HL-1 (HL-1(Pkp2-shRNA)) cells were screened for 750 miRNAs using low-density microfluidic panels. Fifty-nine miRNAs were differentially expressed. MiR-184 was the most downregulated miRNA. Expression of miR-184 in the heart and cardiac myocyte was developmentally downregulated and was low in mature myocytes. MicroRNA-184 was predominantly expressed in cardiac mesenchymal progenitor cells. Knockdown of Pkp2 in cardiac mesenchymal progenitor cells also reduced miR-184 levels. Expression of miR-184 was transcriptionally regulated by the E2F1 pathway, which was suppressed in PKP2-deficient cells. Activation of E2F1, on overexpression of its activator CCND1 (cyclin D1) or knockdown of its inhibitor retinoblastoma 1, partially rescued miR-184 levels. In addition, DNA methyltransferase-1 was recruited to the promoter region of miR-184, and the CpG sites at the upstream region of miR-184 were hypermethylated. Treatment with 5-aza-2'-deoxycytidine, a demethylation agent, and knockdown of DNA methyltransferase-1 partially rescued miR-184 level. Pathway analysis of paired miR-184:mRNA targets identified cell proliferation, differentiation, and death as the main affected biological processes. Knockdown of miR-184 in HL-1 cells and mesenchymal progenitor cells induced and, conversely, its overexpression attenuated adipogenesis. PKP2 deficiency leads to suppression of the E2F1 pathway and hypermethylation of the CpG sites at miR-184 promoter, resulting in downregulation of miR-184 levels. Suppression of miR-184 enhances and its activation attenuates adipogenesis in vitro. Thus, miR-184 contributes to the pathogenesis of adipogenesis in PKP2-deficient cells.

Discovery and Genomic Characterization of a Novel Bat Sapovirus with Unusual Genomic Features and Phylogenetic Position

Sapovirus is a genus of caliciviruses that are known to cause enteric disease in humans and animals. There is considerable genetic diversity among the sapoviruses, which are classified into different genogroups based on phylogenetic analysis of the full-length capsid protein sequence. While several mammalian species, including humans, pigs, minks, and dogs, have been identified as animal hosts for sapoviruses, there were no reports of sapoviruses in bats in spite of their biological diversity. In this report, we present the results of a targeted surveillance study in different bat species in Hong Kong. Five of the 321 specimens from the bat species, Hipposideros pomona, were found to be positive for sapoviruses by RT-PCR. Complete or nearly full-length genome sequences of approximately 7.7 kb in length were obtained for three strains, which showed similar organization of the genome compared to other sapoviruses. Interestingly, they possess many genomic features atypical of most sapoviruses, like high G+C content and minimal CpG suppression. Phylogenetic analysis of the viral proteins suggested that the bat sapovirus descended from an ancestral sapovirus lineage and is most closely related to the porcine sapoviruses. Codon usage analysis showed that the bat sapovirus genome has greater codon usage bias relative to other sapovirus genomes. In summary, we report the discovery and genomic characterization of the first bat calicivirus, which appears to have evolved under different conditions after early divergence from other sapovirus lineages.

Genomic DNA k-mer spectra: models and modalities

Tetrapods, unlike other organisms, have multimodal spectra of k-mers in their genomes

B7-02: Identification of differentially methylated CpG Islands in the early stage of human pulmonary adenocarcinoma

In the development of cancer, a series of tumor suppressor genes is inactivated by point mutation and chromosomal deletion. Aberrant methylation of CpG islands has recently been shown to serve as an alternative way of inactivating such genes in cancer. Methylation of cytosine residues in CpG dinucleotides is an extremely important epigenetic modification of the eukaryotic genome that affects various cell processes. Methylation is known to play an important role in neoplasia by inactivating tumor suppressor genes such as Rb, p16, and estrogen receptor. In this study we used CpG island amplification and suppression subtractive hybridization to identify the aberrant methylation of CpG islands in A/J mouse lung adenoma, in order to clarify whether a hypermethylated homologue gene exists in human pulmonary adenocarcinoma. NNK (4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone)-induced A/J mouse lung adenoma and normal tissue were used in this study. The histology of A/J mouse adenoma mimics the early stage of human pulmonary adenocarcinoma. The genomic DNA was extracted using standard procedures. The PCR-based methylated CpG island amplification (MCA) technique was used for detection of methylated CpG islands. Suppression subtractive hybridization (SSH) and differential screening (DS) were used to identify the differentially methylated sequence in A/J mouse lung adenoma tissue. For the genes selected after SSH and DS analysis, quantitative Real-time PCR was used to check their expression level in A/J mouse lung adenoma tissue and normal tissue. Using homoloGene software, human homologue genes were then detected. Real-time PCR was also performed to check the expression level of these genes in human lung adenocarcinoma and normal lung tissue. Bisulfite genomic sequencing was done to confirm the methylation status of these genes in the promotor area. After MCA-SSH and DS, 117 clones from a total of 1000 were selected for sequence analysis. Four genes were identified in the DNA fragments from 117 clones from A/J mouse lung tumor tissue. These were KIF21A (kinesin family member 21A), Samd14 (sterile alpha motif domain containing 14), EG436235 (similar to regulator of G-protein signal 3) and Vwa1 (von Willebrand factor A domain containing 1). Among these four genes, two gene fragments had a genomic region that met the criteria for a CpG island. RT-PCR showed differences in expression level between A/J mouse normal lung tissue and lung adenoma tissue for the KIF21A (p<0.005), Samd14 (p=0.005) and EG436235 (p<0.01) genes. Further study of the human homologues of KIF21A (p<0.01) and Samd14 (p<0.002) showed that they were expressed in normal lung tissue but markedly down-regulated in lung adenocarcinoma tissue. Bisulfite sequencing revealed that the promtor region of Samd14 in human lung adenocarcinoma was methylated more frequently than that of its normal counterpart. The expression level of the KIF21A (location: 12q12) and Samd14 (location: 17q21.33) genes is down regulated in human pulmonary adenocarcinoma tissue. The methylation frequency of the Samd14 gene in the promotor region is higher in human pulmonary adenocarcinoma than in normal lung tissue. These data suggest that promotor methylation of these genes is a specific event in pulmonary adenocarcinogenesis, and that their down-regulation may be functionally associated with malignant progression of lung adenocarcinoma

Intron length and accelerated 3′ gene evolution

Genetic evolution depends in part upon a balance between negative selection and environmentally driven mutation. To explore whether this balance is affected by gene structure, we have used phylogenetic data mining to compare gene compositions across a range of species. Here we show that genomes of higher species exhibit a greater frequency of 5′ CpG islands and of CpG→TpG/CpA transitions. This latter mutational pattern exhibits a 5′-to-3′ trend in higher species, consistent with a length-dependent effect on methylation-dependent CpG suppression. Associated strand asymmetry (TpG>CpA) declines with gene length, implying attenuation of transcription-coupled repair 3′ to introns. A sharp 3′ rise in coding region single-nucleotide polymorphism frequency further supports a mechanistic role for intron length in promoting genetic variation by reducing repair and/or weakening negative selection. Consistent with this, the Ka/Ks ratio of 3′ exons exceeds that of centrally located exons in intron-containing, but not in intronless, genes (p<0.0003). We conclude that the efficiency of transcription-coupled repair decreases with gene length, suggesting in turn that 3′ gene evolution is accelerated both by introns and by gene methylation.

CodonO: a new informatics method for measuring synonymous codon usage bias within and across genomes

Codon usage bias reveals the silent molecular evolution pattern. Previous research showed the codon usage bias was associated with many biological processes, such as protein expression level, genetic coding error minimization, mRNA stability, codon context, tRNA richness, CpG suppression, DNA methylation, and tissue or organ specificity. In this paper we reviewed major methods of codon usage bias measurement. Since most of existing methods are only suitable for the comparison of codon usage bias within a single genome, we introduced a new informatics index, referred to as synonymous codon usage order (SCUO), to measure synonymous codon usage bias within and across genomes. In this method, Shannon informational theory was applied to describe the SCUO of each gene using a value ranging from 0 to 1, with larger values associated with greater codon usage bias. We compared our method with the codon adaptation index (CAI) method for measuring codon bias in Escherichia coli and Sacharomyces cerevisiae. We also studied the correlation between SCUO and CAI, and the relation of SCUO with gene length, gene function, and GC composition. In addition, we explored the correlation between SCUO and mRNA abundance in S. cerevisiae using SAGE expression data. The software package codonO is freely available at http://digbio.missouri.edu/∼wanx/cu/codonO/.

The ACIN1 Gene is Hypermethylated in Early Stage Lung Adenocarcinoma

In recent years, many studies have performed genome-wide searching for differentially methylated genes in cancer. We hypothesized that characteristic aberrant hypermethylation of CpG islands of certain genes may exist in the early stages of lung adenocarcinoma and that such alterations may be useful in the detection and treatment of early lung adenocarcinoma. A pair of immortalized cell lines originating from atypical adenomatous hyperplasia (PL16T) and from the resected end of the bronchus of the same patient (PL16B) was searched for aberrantly and differentially hypermethylated DNA fragments by a combination of the methylated CpG island amplification and suppression subtractive hybridization methods. From 229 clones, we selected 15 fragments that had a genomic region meeting the criteria for a CpG island. We identified a gene, apoptotic chromatin condensation inducer 1 (ACIN1), that was hypermethylated in PL16T. A higher frequency of hypermethylation at a locus at the 5': end of the DNA fragment isolated from the ACIN1 gene was found in small-sized adenocarcinoma (2 cm or less) (30/37, 81%) compared with normal lung tissue (9/37, 24%, p < 0.05). Interestingly, hypermethylation of ACIN1 was detected relatively frequently in the normal counterpart of adenocarcinoma without bronchioloalveolar carcinoma (BAC) component (7/16, 44%), but was rare in the normal counterpart of adenocarcinoma with BAC component (2/21, 10%, P < 0.05). We found hypermethylation of the ACIN1 gene in early stage lung adenocarcinoma. The role of methylation status in the development and malignant transformation of lung adenocarcinoma requires clarification.

The molecular basis for the lack of immunostimulatory activity of vertebrate DNA.

Macrophages and B cells are activated by unmethylated CpG-containing sequences in bacterial DNA. The lack of activity of self DNA has generally been attributed to CpG suppression and methylation, although the role of methylation is in doubt. The frequency of CpG in the mouse genome is 12.5% of Escherichia coli, with unmethylated CpG occurring at approximately 3% the frequency of E. coli. This suppression of CpG alone is insufficient to explain the inactivity of self DNA; vertebrate DNA was inactive at 100 micro g/ml, 3000 times the concentration at which E. coli DNA activity was observed. We sought to resolve why self DNA does not activate macrophages. Known active CpG motifs occurred in the mouse genome at 18% of random occurrence, similar to general CpG suppression. To examine the contribution of methylation, genomic DNAs were PCR amplified. Removal of methylation from the mouse genome revealed activity that was 23-fold lower than E. coli DNA, although there is only a 7-fold lower frequency of known active CpG motifs in the mouse genome. This discrepancy may be explained by G-rich sequences such as GGAGGGG, which potently inhibited activation and are found in greater frequency in the mouse than the E. coli genome. In summary, general CpG suppression, CpG methylation, inhibitory motifs, and saturable DNA uptake combined to explain the inactivity of self DNA. The immunostimulatory activity of DNA is determined by the frequency of unmethylated stimulatory sequences within an individual DNA strand and the ratio of stimulatory to inhibitory sequences.

Genomics: more than the sum of the parts.

It has been known for some time that DNA composition varies across a given genome as well as between genomes (Filipski et al. 1973;Wagner and Capesius 1981). Genomic sequencing projects allow this observation to be confirmed at the sequence level (The Arabidopsis Genome Initiative 2000; Ashikawa 2001). However, the cause and function of these compositional differences are still obscure. Among the theories that may explain these phenomena (Eyre-Walker and Hurst 2001), mutation bias from C to T due to deamination of methylated C has been commonly used to account for them (Coulondre et al. 1978). As methylation is probably involved in a mechanism to silence transposable elements (Martienssen 1998), it makes sense that inactive methylated transposons can easily undergo C to T transition because they are under no selective pressure. However, this theory cannot explain other related compositional biases such as the CpG suppression observed in animal mitochondria (Cardon et al. 1994), where there is no DNA methylation.

Comprehensive analysis of CpG islands in human chromosomes 21 and 22.

CpG islands are useful markers for genes in organisms containing 5-methylcytosine in their genomes. In addition, CpG islands located in the promoter regions of genes can play important roles in gene silencing during processes such as X-chromosome inactivation, imprinting, and silencing of intragenomic parasites. The generally accepted definition of what constitutes a CpG island was proposed in 1987 by Gardiner-Garden and Frommer [Gardiner-Garden, M. & Frommer, M. (1987) J. Mol. Biol. 196, 261-282] as being a 200-bp stretch of DNA with a C+G content of 50% and an observed CpG/expected CpG in excess of 0.6. Any definition of a CpG island is somewhat arbitrary, and this one, which was derived before the sequencing of mammalian genomes, will include many sequences that are not necessarily associated with controlling regions of genes but rather are associated with intragenomic parasites. We have therefore used the complete genomic sequences of human chromosomes 21 and 22 to examine the properties of CpG islands in different sequence classes by using a search algorithm that we have developed. Regions of DNA of greater than 500 bp with a G+C equal to or greater than 55% and observed CpG/expected CpG of 0.65 were more likely to be associated with the 5' regions of genes and this definition excluded most Alu-repetitive elements. We also used genome sequences to show strong CpG suppression in the human genome and slight suppression in Drosophila melanogaster and Saccharomyces cerevisiae. This finding is compatible with the recent detection of 5-methylcytosine in Drosophila, and might suggest that S. cerevisiae has, or once had, CpG methylation.