Gene-associated CpG islands in plants as revealed by analyses of genomic sequences

We screened plant genome sequences, primarily from rice and Arabidopsis thaliana, for CpG islands, and identified DNA segments rich in CpG dinucleotides within these sequences. These CpG-rich clusters appeared in the analysed sequences as discrete peaks and occurred at the frequencies of one per 4.7 kb in rice and one per 4.0 kb in A. thaliana. In rice and A. thaliana, most of the CpG-rich clusters were associated with genes, which suggests that these clusters are useful landmarks in genome sequences for identifying genes in plants with small genomes. In contrast, in plants with larger genomes, only a few of the clusters were associated with genes. These plant CpG-rich clusters satisfied the criteria used for identifying human CpG islands, which suggests that these CpG clusters may be regarded as plant CpG islands. The position of each island relative to the 5'-end of its associated gene varied considerably. Genes in the analysed sequences were grouped into five classes according to the position of the CpG islands within their associated genes. A large proportion of the genes belonged to one of two classes, in which a CpG island occurred near the 5'-end of the gene or covered the whole gene region. The position of a plant CpG island within its associated gene appeared to be related to the extent of tissue-specific expression of the gene; the CpG islands of most of the widely expressed rice genes occurred near the 5'-end of the genes.     

CpGIF: an algorithm for the identification of CpG islands

CpG islands (CGIs) play a fundamental role in genome analysis and annotation, and contribute to improving the accuracy of promoter prediction. Besides, CGIs in promoter regions are abnormally methylated in cancer cells and thus can be used as tumor markers. However, current methods for identifying CGIs suffer from various drawbacks. We present a new algorithm for detecting CGIs, called CpG Island Finder (CpGIF), which combines the best features in the most commonly used algorithms and avoids their disadvantages as much as possible. Five public tools for CpG island searching are used to compare with CpGIF for the assessment of accuracy and computational efficiency. The results reveal that CpGIF has higher performance coefficient and correlation coefficient than these previous methods, which indicates that CpGIF is able to provide high sensitivity and specificity at the same time. CpGIF is also faster than those methods with comparable prediction accuracy.     

CpG islands are the second main factor shaping codon usage in human genes

A correspondence analysis of codon usage in human genes revealed, as expected, that the first axis is strongly correlated with the base composition at synonymous third codon positions. At one extreme of the second axis were localized genes with a high frequency of NCG and CGN codons. The great majority of these sequences were embedded in CpG islands, while the opposite is true for the genes placed at the other extreme. The two main conclusions of this paper are: (1) the influence of CpG islands on codon usage, and (2) since the second axis is orthogonal (and therefore independent) of the first, GC3-rich genes are not necessarily associated with CpG islands.     

Segmental distribution of genes harboring a CpG island-like region on rice chromosomes

In plant genomes, there exist discrete regions rich in CpG dinucleotides, namely CpG clusters. In rice, most of these CpG clusters are associated with genes. Rice genes are grouped into one of the five classes according to the position of an associated CpG cluster. Among them, class 1 genes, which harbor a CpG cluster at the 5′-terminus, share similarities with human genes having CpG islands. In the present study, by analyzing plant genome sequence data, primarily from rice, we investigated the chromosomal distribution of genes of each class, mainly class 1 genes. Class 1 genes were not uniformly distributed across the rice genome, but were clustered into discrete chromosomal segments. EST-based analysis of the distribution of expressed genes indicates that this segmental distribution of class 1 genes caused a preferential distribution of expressed genes within class 1 gene-rich segments. We then compared the methylation status of genes of each class to examine the possibility that differential DNA methylation, if any, is relevant to the observed differential expression level of genes inside and outside the class 1 segments. The difference in the methylation level between these genes was revealed to be fairly small, which does not support the above-mentioned possibility. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Mutations and CpG islands among hepatitis B virus genotypes in Europe

Background

Hepatitis B virus (HBV) genotypes have a distinct geographical distribution and influence disease progression and treatment outcomes. The purpose of this study was to investigate the distribution of HBV genotypes in Europe, the impact of mutation of different genotypes on HBV gene abnormalities, the features of CpG islands in each genotype and their potential role in epigenetic regulation.

Results

Of 383 HBV isolates from European patients, HBV genotypes A-G were identified, with the most frequent being genotype D (51.96%) in 12 countries, followed by A (39.16%) in 7 countries, and then E (3.66%), G (2.87%), B (1.57%), F (0.52%) and C (0.26%). A higher rate of mutant isolates were identified in those with genotype D (46.7%) followed by G (45.5%), and mutations were associated with structural and functional abnormalities of HBV genes. Conventional CpG island I was observed in genotypes A, B, C, D and E. Conventional islands II and III were detected in all A-G genotypes. A novel CpG island IV was found in genotypes A, D and E, and island V was only observed in genotype F. The A-G genotypes lacked the novel CpG island VI. “Split” CpG island I in genotypes D and E and “split” island II in genotypes A, D, E, F and G were observed. Two mutant isolates from genotype D and one from E were found to lack both CpG islands I and III.

Conclusions

HBV genotypes A-G were identified in European patients. Structural and functional abnormalities of HBV genes were caused by mutations leading to the association of genotypes D and G with increased severity of liver disease. The distribution, length and genetic traits of CpG islands were different between genotypes and their biological and clinical significances warrant further study, which will help us better understand the potential role of CpG islands in epigenetic regulation of the HBV genome.     

Sequence context at human single nucleotide polymorphisms: overrepresentation of CpG dinucleotide at polymorphic sites and suppression of variation in CpG islands

Human polymorphisms originate as mutations, and the influence of context on mutagenesis should be reflected in the distribution of sequences surrounding single nucleotide polymorphisms (SNPs). We have performed a computational survey of nearly two million human SNPs to determine if sequence-dependent hotspots for polymorphism exist in the human genome. Here we show that sequences containing CpG dinucleotides, which occur at low frequencies in the human genome, are 6.7-fold more abundant at polymorphic sites than expected. In contrast, polymorphisms in CpG sequences located within CpG islands, important regulatory regions that modulate gene expression, are 6.8-fold less prevalent than expected. The distribution of polymorphic alleles at CpGs in CpG islands is also significantly different from that in non-island regions. These data strongly support a role for 5-methylcytosine deamination in the generation of human variation, and suggest that variation at CpGs in islands is suppressed.     

Isolation and enrichment of human genomic CpG islands by methylation-sensitive mirror orientation selection

CpG islands (CGIs) in human genomic DNA are GC-rich fragments whose aberrant methylation is associated with human disease development. In the current study, methylation-sensitive mirror orientation selection (MS-MOS) was developed to efficiently isolate and enrich unmethylated CGIs from human genomic DNA. The unmethylated CGIs prepared by the MS-MOS procedure subsequently were used to construct a CGI library. Then the sequence characteristics of cloned inserts of the library were analyzed by bioinformatics tools, and the methylation status of CGI clones was analyzed by HpaII PCR. The results showed that the MS-MOS method could be used to isolate up to 0.001% of differentially existed unmethylated DNA fragments in two complex genomic DNA. In the CGI library, 34.1% of clones had insert sequences satisfying the minimal criteria for CGIs. Excluding duplicates, 22.0% of the 80,000 clones were unique CGI clones, representing 60% of all the predicted CGIs (about 29,000) in human genomic DNA, and most or all of the CGI clones were unmethylated in human normal cell DNA based on the HpaII PCR analysis results of randomly selected CGI clones. In conclusion, MS-MOS was an efficient way to isolate and enrich human genomic CGIs. The method has powerful potential application in the comprehensive identification of aberrantly methylated CGIs associated with human tumorigenesis to improve understanding of the epigenetic mechanisms involved.     

The clustering of CpG islands may constitute an important determinant of the 3D organization of interphase chromosomes

We used the 4C-Seq technique to characterize the genome-wide patterns of spatial contacts of several CpG islands located on chromosome 14 in cultured chicken lymphoid and erythroid cells. We observed a clear tendency for the spatial clustering of CpG islands present on the same and different chromosomes, regardless of the presence or absence of promoters within these CpG islands. Accordingly, we observed preferential spatial contacts between Sp1 binding motifs and other GC-rich genomic elements, including the DNA sequence motifs capable of forming G-quadruplexes. However, an anchor placed in a gene/CpG island-poor area formed spatial contacts with other gene/CpG island-poor areas on chromosome 14 and other chromosomes. These results corroborate the two-compartment model of the spatial organization of interphase chromosomes and suggest that the clustering of CpG islands constitutes an important determinant of the 3D organization of the eukaryotic genome in the cell nucleus. Using the ChIP-Seq technique, we mapped the genome-wide CTCF deposition sites in the chicken lymphoid and erythroid cells that were used for the 4C analysis. We observed a good correlation between the density of CTCF deposition sites and the level of 4C signals for the anchors located in CpG islands but not for an anchor located in a gene desert. It is thus possible that CTCF contributes to the clustering of CpG islands observed in our experiments.     

Construction of a bovine yeast artificial chromosome (YAC) library

We have constructed a bovine yeast artificial chromosome (YAC) library to provide a common resource for bovine genome research. We used leukocytes of a Japanese black bull ( Bos taurus ) as the DNA source, AB1380 for the yeast host, and pYAC4 for the vector. The library consists of 24 576 clones arranged in 256 96-well microtiter plates. An average insert size estimated from the analysis of 251 randomly selected clones was 480 kb. The rate of chimeric YACs evaluated by fluorescence in situ hybridization (FISH) analysis of 44 randomly selected clones was 36·4%. To estimate the number of genome equivalents, PCR-based screening was performed with 48 primer pairs and isolated 3·2 clones on average. In order to provide broad access for the scientific community, this library has been incorporated into the Reference Library system which provides high density filters for colony hybridization screening and a common database of the library.     

Construction and characterization of a bovine BAC library with four genome-equivalent coverage

A bovine artificial chromosome (BAC) library of 105 984 clones has been constructed in the vector pBeloBAC11 and organized in 3-dimension pools and high density membranes for screening by PCR and hybridization. The average insert size, determined after analysis of 388 clones, was estimated at 120 kb corresponding to a four genome coverage. Given the fact that a male was used to construct the library, the probability of finding any given autosomal and X or Y locus is respectively 0.98 and 0.86. The library was screened for 164 microsatellite markers and an average of 3.9 superpools was positive for each PCR system. None of the 50 or so BAC clones analysed by FISH was chimeric. This BAC library increases the international genome coverage for cattle to around 28 genome equivalents and extends the coverage of the ruminant genomes available at the Inra resource center to 15 genome equivalents.     

GIST: Genomic island suite of tools for predicting genomic islands in genomic sequences

Genomic Islands (GIs) are genomic regions that are originally from other organisms, through a process known as Horizontal Gene Transfer (HGT). Detection of GIs plays a significant role in biomedical research since such align genomic regions usually contain important features, such as pathogenic genes. We have developed a use friendly graphic user interface, Genomic Island Suite of Tools (GIST), which is a platform for scientific users to predict GIs. This software package includes five commonly used tools, AlienHunter, IslandPath, Colombo SIGI-HMM, INDeGenIUS and Pai-Ida. It also includes an optimization program EGID that ensembles the result of existing tools for more accurate prediction. The tools in GIST can be used either separately or sequentially. GIST also includes a downloadable feature that facilitates collecting the input genomes automatically from the FTP server of the National Center for Biotechnology Information (NCBI). GIST was implemented in Java, and was compiled and executed on Linux/Unix operating systems. AVAILABILITY: The database is available for free at http://www5.esu.edu/cpsc/bioinfo/software/GIST.     

Lead exposure induces changes in 5-hydroxymethylcytosine clusters in CpG islands in human embryonic stem cells and umbilical cord blood

Prenatal exposure to neurotoxicants such as lead (Pb) may cause stable changes in the DNA methylation (5mC) profile of the fetal genome. However, few studies have examined its effect on the DNA de-methylation pathway, specifically the dynamic changes of the 5-hydroxymethylcytosine (5hmC) profile. Therefore, in this study, we investigate the relationship between Pb exposure and 5mC and 5hmC modifications during early development. To study the changes in the 5hmC profile, we use a novel modification of the Infinium™ HumanMethylation450 assay (Illumina, Inc.), which we named HMeDIP-450K assay, in an in vitro human embryonic stem cell model of Pb exposure. We model Pb exposure-associated 5hmC changes as clusters of correlated, adjacent CpG sites, which are co-responding to Pb. We further extend our study to look at Pb-dependent changes in high density 5hmC regions in umbilical cord blood DNA from 48 mother-infant pairs from the Early Life Exposure in Mexico to Environmental Toxicants (ELEMENT) cohort. For our study, we randomly selected umbilical cord blood from 24 male and 24 female children from the 1st and 4th quartiles of Pb levels. Our data show that Pb-associated changes in the 5hmC and 5mC profiles can be divided into sex-dependent and sex-independent categories. Interestingly, differential 5mC sites are better markers of Pb-associated sex-dependent changes compared to differential 5hmC sites. In this study we identified several 5hmC and 5mC genomic loci, which we believe might have some potential as early biomarkers of prenatal Pb exposure.     

Construction and characterisation of a yeast artificial chromosome library containing three haploid maize genome equivalents

We have constructed a yeast artificial chromosome (YAC) library using high-molecular-weight DNA prepared from agarose-embedded leaf protoplasts of the maize inbred line UE95. This library contains 79 000 clones with an average insert size of 145 kb and should therefore represent approximately three haploid genome equivalents. The library is organised as an ordered array in duplicate microtitre plates. Forty-one pools of DNA from 1920 individual clones have been prepared for rapid screening of the library by the polymerase chain reaction (PCR). Using this approach, together with conventional colony hybridisation, we have been able to identify between one and eight positive clones for every probe used.