CpG island density and its correlations with genomic features in mammalian genomes

Germline methylation patterns determine the distribution of recombination events in the dog genome

The positive-regulatory domain containing nine gene, PRDM9, which strongly associates with the location of recombination events in several vertebrates, is inferred to be inactive in the dog genome. Here, we address several questions regarding the control of recombination and its influence on genome evolution in dogs. First, we address whether the association between CpG islands (CGIs) and recombination hotspots is generated by lack of methylation, GC-biased gene conversion (gBGC), or both. Using a genome-wide dog single nucleotide polymorphism data set and comparisons of the dog genome with related species, we show that recombination-associated CGIs have low CpG mutation rates, and that CpG mutation rate is negatively correlated with recombination rate genome wide, indicating that nonmethylation attracts the recombination machinery. We next use a neighbor-dependent model of nucleotide substitution to disentangle the effects of CpG mutability and gBGC and analyze the effects that loss of PRDM9 has on these rates. We infer that methylation patterns have been stable during canid genome evolution, but that dog CGIs have experienced a drastic increase in substitution rate due to gBGC, consistent with increased levels of recombination in these regions. We also show that gBGC is likely to have generated many new CGIs in the dog genome, but these mostly occur away from genes, whereas the number of CGIs in gene promoter regions has not increased greatly in recent evolutionary history. Recombination has a major impact on the distribution of CGIs that are detected in the dog genome due to the interaction between methylation and gBGC. The results indicate that germline methylation patterns are the main determinant of recombination rates in the absence of PRDM9.

Whole genome sequencing and methylome analysis of the wild guinea pig

Background

DNA methylation is a heritable mechanism that acts in response to environmental changes, lifestyle and diseases by influencing gene expression in eukaryotes. Epigenetic studies of wild organisms are mandatory to understand their role in e.g. adaptational processes in the great variety of ecological niches. However, strategies to address those questions on a methylome scale are widely missing. In this study we present such a strategy and describe a whole genome sequence and methylome analysis of the wild guinea pig.

Results

We generated a full Wild guinea pig (Cavia aperea) genome sequence with enhanced coverage of methylated regions, benefiting from the available sequence of the domesticated relative Cavia porcellus. This new genome sequence was then used as reference to map the sequence reads of bisulfite treated Wild guinea pig sequencing libraries to investigate DNA-methylation patterns at nucleotide-specific level, by using our here described method, named ‘DNA-enrichment-bisulfite-sequencing’ (MEBS). The results achieved using MEBS matched those of standard methods in other mammalian model species. The technique is cost efficient, and incorporates both methylation enrichment results and a nucleotide-specific resolution even without a whole genome sequence available. Thus MEBS can be easily applied to extend methylation enrichment studies to a nucleotide-specific level.

Conclusions

The approach is suited to study methylomes of not yet sequenced mammals at single nucleotide resolution. The strategy is transferable to other mammalian species by applying the nuclear genome sequence of a close relative. It is therefore of interest for studies on a variety of wild species trying to answer evolutionary, adaptational, ecological or medical questions by epigenetic mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1036) contains supplementary material, which is available to authorized users.

Orphan and gene related CpG Islands follow power-law-like distributions in several genomes: evidence of function-related and taxonomy-related modes of distribution

CpG Islands (CGIs) are compositionally defined short genomic stretches, which have been studied in the human, mouse, chicken and later in several other genomes. Initially, they were assigned the role of transcriptional regulation of protein-coding genes, especially the house-keeping ones, while more recently there is found evidence that they are involved in several other functions as well, which might include regulation of the expression of RNA genes, DNA replication etc. Here, an investigation of their distributional characteristics in a variety of genomes is undertaken for both whole CGI populations as well as for CGI subsets that lie away from known genes (gene-unrelated or "orphan" CGIs). In both cases power-law-like linearity in double logarithmic scale is found. An evolutionary model, initially put forward for the explanation of a similar pattern found in gene populations is implemented. It includes segmental duplication events and eliminations of most of the duplicated CGIs, while a moderate rate of non-duplicated CGI eliminations is also applied in some cases. Simulations reproduce all the main features of the observed inter-CGI chromosomal size distributions. Our results on power-law-like linearity found in orphan CGI populations suggest that the observed distributional pattern is independent of the analogous pattern that protein coding segments were reported to follow. The power-law-like patterns in the genomic distributions of CGIs described herein are found to be compatible with several other features of the composition, abundance or functional role of CGIs reported in the current literature across several genomes, on the basis of the proposed evolutionary model.

Mammalian epigenetic mechanisms

The mammalian genome is packaged into chromatin that is further compacted into three-dimensional structures consisting of distinct functional domains. The higher order structure of chromatin is in part dictated by enzymatic DNA methylation and histone modifications to establish epigenetic layers controlling gene expression and cellular functions, without altering the underlying DNA sequences. Apart from DNA and histone modifications, non-coding RNAs can also regulate the dynamics of the mammalian gene expression and various physiological functions including cell division, differentiation, and apoptosis. Aberrant epigenetic signatures are associated with abnormal developmental processes and diseases such as cancer. In this review, we will discuss the different layers of epigenetic regulation, including writer enzymes for DNA methylation, histone modifications, non-coding RNA, and chromatin conformation. We will highlight the combinatorial role of these structural and chemical modifications along with their partners in various cellular processes in mammalian cells. We will also address the cis and trans interacting "reader" proteins that recognize these modifications and "eraser" enzymes that remove these marks. Furthermore, an attempt will be made to discuss the interplay between various epigenetic writers, readers, and erasures in the establishment of mammalian epigenetic mechanisms.

Genetic recombination is targeted towards gene promoter regions in dogs

The identification of the H3K4 trimethylase, PRDM9, as the gene responsible for recombination hotspot localization has provided considerable insight into the mechanisms by which recombination is initiated in mammals. However, uniquely amongst mammals, canids appear to lack a functional version of PRDM9 and may therefore provide a model for understanding recombination that occurs in the absence of PRDM9, and thus how PRDM9 functions to shape the recombination landscape. We have constructed a fine-scale genetic map from patterns of linkage disequilibrium assessed using high-throughput sequence data from 51 free-ranging dogs, Canis lupus familiaris. While broad-scale properties of recombination appear similar to other mammalian species, our fine-scale estimates indicate that canine highly elevated recombination rates are observed in the vicinity of CpG rich regions including gene promoter regions, but show little association with H3K4 trimethylation marks identified in spermatocytes. By comparison to genomic data from the Andean fox, Lycalopex culpaeus, we show that biased gene conversion is a plausible mechanism by which the high CpG content of the dog genome could have occurred.

Recombination in mammalian genomes tends to occur within highly localized regions known as recombination hotspots. These hotspots appear to be a ubiquitous feature of mammalian genomes, but tend to not be shared between closely related species despite high levels of DNA sequence similarity. This disparity has been largely explained by the discovery of PRDM9 as the gene responsible for localizing recombination hotspots via recognition and binding to specific DNA motifs. Variation within PRDM9 can lead to changes to the recognized motif, and hence changes to the location of recombination hotspots thought the genome. Multiple studies have shown that PRDM9 is under strong selective pressure, apparently leading to a rapid turnover of hotspot locations between species. However, uniquely amongst mammals, PRDM9 appears to be dysfunctional in dogs and other canids. In this paper, we investigate how the loss of PRDM9 has affected the fine-scale recombination landscape in dogs and contrast this with patterns seen in other species.

Cigarette smoking and DNA methylation

DNA methylation is the most studied epigenetic modification, capable of controlling gene expression in the contexts of normal traits or diseases. It is highly dynamic during early embryogenesis and remains relatively stable throughout life, and such patterns are intricately related to human development. DNA methylation is a quantitative trait determined by a complex interplay of genetic and environmental factors. Genetic variants at a specific locus can influence both regional and distant DNA methylation. The environment can have varying effects on DNA methylation depending on when the exposure occurs, such as during prenatal life or during adulthood. In particular, cigarette smoking in the context of both current smoking and prenatal exposure is a strong modifier of DNA methylation. Epigenome-wide association studies have uncovered candidate genes associated with cigarette smoking that have biologically relevant functions in the etiology of smoking-related diseases. As such, DNA methylation is a potential mechanistic link between current smoking and cancer, as well as prenatal cigarette-smoke exposure and the development of adult chronic diseases.

Investigating the relationship of DNA methylation with mutation rate and allele frequency in the human genome

Background

DNA methylation, which mainly occurs at CpG dinucleotides, is a dynamic epigenetic regulation mechanism in most eukaryotic genomes. It is already known that methylated CpG dinucleotides can lead to a high rate of C to T mutation at these sites. However, less is known about whether and how the methylation level causes a different mutation rate, especially at the single-base resolution.

Results

In this study, we used genome-wide single-base resolution methylation data to perform a comprehensive analysis of the mutation rate of methylated cytosines from human embryonic stem cell. Through the analysis of the density of single nucleotide polymorphisms, we first confirmed that the mutation rate in methylated CpG sites is greater than that in unmethylated CpG sites. Then, we showed that among methylated CpG sites, the mutation rate is markedly increased in low-intermediately (20-40% methylation level) to intermediately methylated CpG sites (40-60% methylation level) of the human genome. This mutation pattern was observed regardless of DNA strand direction and the sequence coverage over the site on which the methylation level was calculated. Moreover, this highly non-random mutation pattern was found more apparent in intergenic and intronic regions than in promoter regions and CpG islands. Our investigation suggested this pattern appears primarily in autosomes rather than sex chromosomes. Further analysis based on human-chimpanzee divergence confirmed these observations. Finally, we observed a significant correlation between the methylation level and cytosine allele frequency.

Conclusions

Our results showed a high mutation rate in low-intermediately to intermediately methylated CpG sites at different scales, from the categorized genomic region, whole chromosome, to the whole genome level, thereby providing the first supporting evidence of mutation rate variation at human methylated CpG sites using the genome-wide sing-base resolution methylation data.

Evaluation of 5-methyl-2'-deoxycytidine stability in hydrolyzed and nonhydrolyzed DNA by HPLC-UV

BACKGROUND

Although the determination of 5-methyl-2'-deoxycytidine (5-MedC) in various biological samples is gaining increasing scientific interest, there are no data available regarding its stability.

RESULTS

We have currently evaluated the stability of 5-MedC and 2'-deoxycytidine (dC) at -20°C, both in hydrolyzed and nonhydrolyzed calf thymus DNA (CT DNA), as well as following repetitive freeze-thaw cycles. HPLC-UV was used for the accurate determination of the two 2'-deoxynucleosides. Statistical evaluation of the results revealed that 5-MedC and dC were stable in hydrolyzed CT DNA for at least 7 days and in nonhydrolyzed CT DNA for at least 65 days, when these were stored at -20°C. Furthermore, both 2'-deoxynucleosides were stable for at least three repetitive freeze-thaw cycles.

CONCLUSION

By using HPLC-UV, we have evaluated the stability of 5-MedC and dC under storage conditions and repetitive freeze-thaw cycles. Our results are informative about the way samples should be handled and stored in epigenetic studies.

Features of methylation and gene expression in the promoter-associated CpG islands using human methylome data

CpG islands are typically located in the 5′ end of genes and considered as gene markers because they play important roles in gene regulation via epigenetic change. In this study, we compared the features of CpG islands identified by several major algorithms by setting the parameter cutoff values in order to obtain a similar number of CpG islands in a genome. This approach allows us to systematically compare the methylation and gene expression patterns in the identified CpG islands. We found that Takai and Jones' algorithm tends to identify longer CpG islands but with weaker CpG island features (e.g., lower GC content and lower ratio of the observed over expected CpGs) and higher methylation level. Conversely, the CpG clusters identified by Hackenberg et al.'s algorithm using stringent criteria are shorter and have stronger features and lower methylation level. In addition, we used the genome-wide base-resolution methylation profile in two cell lines to show that genes with a lower methylation level at the promoter-associated CpG islands tend to express in more tissues and have stronger expression. Our results validated that the DNA methylation of promoter-associated CpG islands suppresses gene expression at the genome level.

Long-range autocorrelations of CpG islands in the human genome

In this paper, we use a statistical estimator developed in astrophysics to study the distribution and organization of features of the human genome. Using the human reference sequence we quantify the global distribution of CpG islands (CGI) in each chromosome and demonstrate that the organization of the CGI across a chromosome is non-random, exhibits surprisingly long range correlations (10 Mb) and varies significantly among chromosomes. These correlations of CGI summarize functional properties of the genome that are not captured when considering variation in any particular separate (and local) feature. The demonstration of the proposed methods to quantify the organization of CGI in the human genome forms the basis of future studies. The most illuminating of these will assess the potential impact on phenotypic variation of inter-individual variation in the organization of the functional features of the genome within and among chromosomes, and among individuals for particular chromosomes.

LINE-1 hypomethylation in familial and sporadic cancer

Increased and decreased methylation at specific sequences (hypermethylation and hypomethylation, respectively) is characteristic of tumor DNA compared to normal DNA and promotes carcinogenesis in multiple ways including genomic instability. Long interspersed element (LINE), an abundant class of retrotransposons, provides a surrogate marker for global hypomethylation. We developed methylation-specific multiplex ligation-dependent probe amplification assays to study LINE-1 methylation in cases of colorectal, gastric, and endometrial cancer (N = 276), stratified by patient category [sporadic; Lynch syndrome (LS); familial colorectal cancer type X (FCCX)] and microsatellite instability status. Within each patient group, LINE-1 showed lower methylation in tumor DNA relative to paired normal DNA and hypomethylation was statistically significant in most cases. Interestingly, normal colorectal mucosa samples from different patient groups displayed differences in LINE-1 methylation that mirrored differences between the respective tumor tissues, with a decreasing trend for LINE-1 methylation from patients with sporadic colorectal cancer to LS to FCCX. Despite the fact that the degree of LINE-1 methylation is generally tissue specific, normal colorectal mucosa, gastric mucosa, and endometrium from LS patients showed similar levels of LINE-1 methylation. Our results suggest that the degree of LINE-1 methylation may constitute a “field defect” that may predispose normal tissues for cancer development.

Functional complementation between transcriptional methylation regulation and post-transcriptional microRNA regulation in the human genome

BACKGROUND

DNA methylation in the 5' promoter regions of genes and microRNA (miRNA) regulation at the 3' untranslated regions (UTRs) are two major epigenetic regulation mechanisms in most eukaryotes. Both DNA methylation and miRNA regulation can suppress gene expression and their corresponding protein product; thus, they play critical roles in cellular processes. Although there have been numerous investigations of gene regulation by methylation changes and miRNAs, there is no systematic genome-wide examination of their coordinated effects in any organism.

RESULTS

In this study, we investigated the relationship between promoter methylation at the transcription level and miRNA regulation at the post-transcription level by taking advantage of recently released human methylome data and high quality miRNA and other gene annotation data. We found methylation level in the promoter regions and expression level was negatively correlated. Then, we showed that miRNAs tended to target the genes with a low DNA methylation level in their promoter regions. We further demonstrated that this observed pattern was not attributed to the gene expression level, expression broadness, or the number of transcription factor binding sites. Interestingly, we found miRNA target sites were significantly enriched in the genes located in differentially methylated regions or partially methylated domains. Finally, we explored the features of DNA methylation and miRNA regulation in cancer genes and found cancer genes tended to have low methylation level and more miRNA target sites.

CONCLUSION

This is the first genome-wide investigation of the combined regulation of gene expression. Our results supported a complementary regulation between DNA methylation (transcriptional level) and miRNA function (post-transcriptional level) in the human genome. The results were helpful for our understanding of the evolutionary forces towards organisms' complexity beyond traditional sequence level investigation.

Mutations and epimutations in the origin of cancer

Cancer is traditionally viewed as a disease of abnormal cell proliferation controlled by a series of mutations. Mutations typically affect oncogenes or tumor suppressor genes thereby conferring growth advantage. Genomic instability facilitates mutation accumulation. Recent findings demonstrate that activation of oncogenes and inactivation of tumor suppressor genes, as well as genomic instability, can be achieved by epigenetic mechanisms as well. Unlike genetic mutations, epimutations do not change the base sequence of DNA and are potentially reversible. Similar to genetic mutations, epimutations are associated with specific patterns of gene expression that are heritable through cell divisions. Knudson's hypothesis postulates that inactivation of tumor suppressor genes requires two hits, with the first hit occurring either in somatic cells (sporadic cancer) or in the germline (hereditary cancer) and the second one always being somatic. Studies on hereditary and sporadic forms of colorectal carcinoma have made it evident that, apart from genetic mutations, epimutations may serve as either hit or both. Furthermore, recent next-generation sequencing studies show that epigenetic genes, such as those encoding histone modifying enzymes and subunits for chromatin remodeling systems, are themselves frequent targets of somatic mutations in cancer and can act like tumor suppressor genes or oncogenes. This review discusses genetic vs. epigenetic origin of cancer, including cancer susceptibility, in light of recent discoveries. Situations in which mutations and epimutations occur to serve analogous purposes are highlighted.

Analysis of CpG methylation sites and CGI among human papillomavirus DNA genomes

Background

The Human Papillomavirus (HPV) genome is divided into early and late coding sequences, including 8 open reading frames (ORFs) and a regulatory region (LCR). Viral gene expression may be regulated through epigenetic mechanisms, including cytosine methylation at CpG dinucleotides. We have analyzed the distribution of CpG sites and CpG islands/clusters (CGI) among 92 different HPV genomes grouped in function of their preferential tropism: cutaneous or mucosal. We calculated the proportion of CpG sites (PCS) for each ORF and calculated the expected CpG values for each viral type.

Results

CpGs are underrepresented in viral genomes. We found a positive correlation between CpG observed and expected values, with mucosal high-risk (HR) virus types showing the smallest O/E ratios. The ranges of the PCS were similar for most genomic regions except E4, where the majority of CpGs are found within islands/clusters. At least one CGI belongs to each E2/E4 region. We found positive correlations between PCS for each viral ORF when compared with the others, except for the LCR against four ORFs and E6 against three other ORFs. The distribution of CpG islands/clusters among HPV groups is heterogeneous and mucosal HR-HPV types exhibit both lower number and shorter island sizes compared to cutaneous and mucosal Low-risk (LR) HPVs (all of them significantly different).

Conclusions

There is a difference between viral and cellular CpG underrepresentation. There are significant correlations between complete genome PCS and a lack of correlations between several genomic region pairs, especially those involving LCR and E6. L2 and L1 ORF behavior is opposite to that of oncogenes E6 and E7. The first pair possesses relatively low numbers of CpG sites clustered in CGIs while the oncogenes possess a relatively high number of CpG sites not associated to CGIs. In all HPVs, E2/E4 is the only region with at least one CGI and shows a higher content of CpG sites in every HPV type with an identified E4. The mucosal HR-HPVs show either the shortest CGI size, followed by the mucosal LR-HPVs and lastly by the cutaneous viral subgroup, and a trend to the lowest CGI number, followed by the cutaneous viral subgroup and lastly by the mucosal LR-HPVs.

Multifractal analysis of nonhyperbolic coupled map lattices: application to genomic sequences

Symbolic sequences generated by coupled map lattices (CMLs) can be used to model the chaotic-like structure of genomic sequences. In this study it is shown that diffusively coupled Chebyshev maps of order 4 (corresponding to a shift of four symbols) very closely reproduce the multifractal spectrum D(q) of human genomic sequences for coupling constant α = 0.35 ± 0.01 if q > 0. The presence of rare configurations causes deviations for q < 0, which disappear if the rare event statistics of the CML is modified. Such rare configurations are known to play specific functional roles in genomic sequences serving as promoters or regulatory elements.

Imprinting in the schizophrenia candidate gene GABRB2 encoding GABA(A) receptor β(2) subunit

Schizophrenia is a complex genetic disorder, the inheritance pattern of which is likely complicated by epigenetic factors yet to be elucidated. In this study, transmission disequilibrium tests with family trios yielded significant differences between paternal and maternal transmissions of the disease-associated single-nucleotide polymorphism (SNP) rs6556547 and its haplotypes. The minor allele (T) of rs6556547 was paternally undertransmitted to male schizophrenic offsprings, and this parent-of-origin effect strongly suggested that GABRB2 is imprinted. 'Flipping' of allelic expression in heterozygotes of SNP rs2229944 (C/T) in GABRB2 or rs2290732 (G/A) in the neighboring GABRA1 was compatible with imprinting effects on gene expression. Clustering analysis of GABRB2 mRNA expressions suggested that imprinting brought about the observed two-tiered distribution of expression levels in controls with heterozygous genotype at the disease-associated SNP rs1816071 (A/G). The deficit of upper-tiered expressions accounted for the lowered expression levels in the schizophrenic heterozygotes. The occurrence of a two-tiered distribution furnished support for imprinting, and also pointed to the necessity of differentiating between two kinds of heterozygotes of different parental origins in disease association studies on GABRB2. Bisulfite sequencing revealed hypermethylation in the neighborhood of SNP rs1816071, and methylation differences between controls and schizophrenia patients. Notably, the two schizophrenia-associated SNPs rs6556547 and rs1816071 overlapped with a CpG dinucleotide, thereby opening the possibility that CpG methylation status of these sites could have an impact on the risk of schizophrenia. Thus multiple lines of evidence pointed to the occurrence of imprinting in the GABRB2 gene and its possible role in the development of schizophrenia.

Anonymous nuclear loci in the white-faced storm-petrel Pelagodroma marina and their applicability to other Procellariiform seabirds

Procellariiform seabirds are among the avian species with the fastest rates of extinction due to interactions with fisheries and introduction of alien predators to the breeding colonies. Conservation and management policies targeting populations of these species must include information on colony demographics and levels of isolation and genetic markers go a long way toward providing reliable estimates of these parameters. To this end, we report isolation and characterization of 14 anonymous nuclear loci, with average length of 657 bp, in the pelagic seabird White-faced Storm-petrel Pelagodroma marina, a species for which there is virtually no genetic information available. These loci, initially isolated from a genomic library built from P. marina, were further tested, for a range of conditions, in 7 other species representing all Procellariiform families. We found high levels of cross-species amplification success, varying between 79% and 86% in representatives of Diomedeidae, Procellariidae, Pelecanoididae, and other Hydrobatidae. We also sequenced 11 loci for 22 P. marina individuals and report higher levels of anonymous genetic variation (π = 0.002), with an average of 1 single nucleotide polymorphism every 100 bp surveyed, relative to the levels found on a typically variable intron in avian species. These markers will be a valuable tool in future population genetics and phylogenetic studies, particularly of nonmodel seabird species.

Bovine leukemia virus integration site selection in cattle that develop leukemia

It is essential for efficient replication of retroviruses that the viral genome is integrated into the host genome after reverse transcription. Some retroviruses are preferentially integrated into certain genomic regions that may differ depending on the disease. In this study, we analyzed the integration site of bovine leukemia virus (BLV) in leukemic cells and 55 integration sites were determined. Although the integration sites were not located in a particular chromosome, the BLV provirus was integrated into transcription units at a frequency of 43.6% (24/55) and the transcriptional direction of the provirus was in accordance with that of the integrated host genes in 62.5% (15/24). The integration sites were located in introns of the host gene, excluding only one site, which was located in downstream from a stop codon. BLV provirus was never found in a protein coding sequence (CDS) in this study. Moreover, the BLV provirus did not favor integration near transcription start sites and CpG islands, or repetitive sequences such as transposons. Therefore, the possibility that the integration of the BLV provirus disrupts the host gene is very low. Although a hot spot was not found in the BLV provirus integration sites, the provirus favored the integration into regions disadvantageous for viral gene expression since no integration site was preferentially located into/near CDS, transcription start site or CpG island. It is suggested that the integration site of the BLV provirus in leukemic cells is related to the suppression of viral gene expression.

Reconstructing the ancestral germ line methylation state of young repeats

One of the key objectives of comparative genomics is the characterization of the forces that shape genomes over the course of evolution. In the last decades, evidence has been accumulated that for vertebrate genomes also epigenetic modifications have to be considered in this context. Especially, the elevated mutation frequency of 5-methylcytosine (5mC) is assumed to facilitate the depletion of CpG dinucleotides in species that exhibit global DNA methylation. For instance, the underrepresentation of CpG dinucleotides in many mammalian genomes is attributed to this effect, which is only neutralized in so-called CpG islands (CGIs) that are preferentially unmethylated and thus partially protected from rapid CpG decay. For primate-specific CpG-rich transposable elements from the ALU family, it is unclear whether their elevated CpG frequency is caused by their small age or by the absence of DNA methylation. In consequence, these elements are often misclassified in CGI annotations. We present a method for the estimation of germ line methylation from pairwise ancestral-descendant alignments. The approach is validated in a simulation study and tested on DNA repeats from the AluSx family. We conclude that a predicted unmethylated state in the germ line is highly correlated with epigenetic activity of the respective genomic region. Thus, CpG-rich repeats can be facilitated as in silico probes for the epigenetic potential of their genomic neighborhood.

Variation in crossover rates across a 3-Mb contig of bread wheat (Triticum aestivum) reveals the presence of a meiotic recombination hotspot

In bread wheat (Triticum aestivum L.), initial studies using deletion lines indicated that crossover (CO) events occur mainly in the telomeric regions of the chromosomes with a possible correlation with the presence of genes. However, little is known about the distribution of COs at the sequence level. To investigate this, we studied in detail the pattern of COs along a contig of 3.110 Mb using two F2 segregating populations (Chinese Spring × Renan (F2-CsRe) and Chinese Spring × Courtot (F2-CsCt)) each containing ~2,000 individuals. The availability of the sequence of the contig from Cs enabled the development of 318 markers among which 23 co-dominant polymorphic markers (11 SSRs and 12 SNPs) were selected for CO distribution analyses. The distribution of CO events was not homogeneous throughout the contig, ranging from 0.05 to 2.77 cM/Mb, but was conserved between the two populations despite very different contig recombination rate averages (0.82 cM/Mb in F2-CsRe vs 0.35 cM/Mb in F2-CsCt). The CO frequency was correlated with the percentage of coding sequence in Cs and with the polymorphism rate between Cs and Re or Ct in both populations, indicating an impact of these two factors on CO distribution. At a finer scale, COs were found in a region covering 2.38 kb, spanning a gene coding for a glycosyl transferase (Hga3), suggesting the presence of a CO hotspot. A non-crossover event covering at least 453 bp was also identified in the same interval. From these results, we can conclude that gene content could be one of the factors driving recombination in bread wheat.

Comparative analysis of CpG islands in four fish genomes

There has been much interest in CpG islands (CGIs), clusters of CpG dinucleotides in GC-rich regions, because they are considered gene markers and involved in gene regulation. To date, there has been no genome-wide analysis of CGIs in the fish genome. We first evaluated the performance of three popular CGI identification algorithms in four fish genomes (tetraodon, stickleback, medaka, and zebrafish). Our results suggest that Takai and Jones' (2002) algorithm is most suitable for comparative analysis of CGIs in the fish genome. Then, we performed a systematic analysis of CGIs in the four fish genomes using Takai and Jones' algorithm, compared to other vertebrate genomes. We found that both the number of CGIs and the CGI density vary greatly among these genomes. Remarkably, each fish genome presents a distinct distribution of CGI density with some genomic factors (e.g., chromosome size and chromosome GC content). These findings are helpful for understanding evolution of fish genomes and the features of fish CGIs.

The evolution of transcription-associated biases of mutations across vertebrates

Background

The interplay between transcription and mutational processes can lead to particular mutation patterns in transcribed regions of the genome. Transcription introduces several biases in mutational patterns; in particular it invokes strand specific mutations. In order to understand the forces that have shaped transcripts during evolution, one has to study mutation patterns associated with transcription across animals.

Results

Using multiple alignments of related species we estimated the regional single-nucleotide substitution patterns along genes in four vertebrate taxa: primates, rodents, laurasiatheria and bony fishes. Our analysis is focused on intronic and intergenic regions and reveals differences in the patterns of substitution asymmetries between mammals and fishes. In mammals, the levels of asymmetries are stronger for genes starting within CpG islands than in genes lacking this property. In contrast to all other species analyzed, we found a mutational pressure in dog and stickleback, promoting an increase of GC-contents in the proximity to transcriptional start sites.

Conclusions

We propose that the asymmetric patterns in transcribed regions are results of transcription associated mutagenic processes and transcription coupled repair, which both seem to evolve in a taxon related manner. We also discuss alternative mechanisms that can generate strand biases and involves error prone DNA polymerases and reverse transcription. A localized increase of the GC content near the transcription start site is a signature of biased gene conversion (BGC) that occurs during recombination and heteroduplex formation. Since dog and stickleback are known to be subject to rapid adaptations due to population bottlenecks and breeding, we further hypothesize that an increase in recombination rates near gene starts has been part of an adaptive process.

Hierarchical multifractal representation of symbolic sequences and application to human chromosomes

The two-dimensional density correlation matrix is constructed for symbolic sequences using contiguous segments of arbitrary size. The multifractal spectrum obtained from this matrix motif is shown to characterize the correlations in the symbolic sequences. This method is applied to entire human chromosomes, shuffled human chromosomes, reconstructed human genomic sequences and to artificial random sequences. It is shown that all human chromosomes have common characteristics in their multifractal spectrum and deviate substantially from random and uncorrelated sequences of the same size. Small deviations are observed between the longer and the shorter chromosomes, especially for the higher (in absolute values) statistical moments. The correlations are crucial for the form of the multifractal spectrum; surrogate shuffled chromosomes present randomlike spectrum, distinctly different from the actual chromosomes. Analytical approaches based on hierarchical superposition of tensor products show that retaining pair correlations in the sequences leads to a closer representation of the genomic multifractal spectra, especially in the region of negative exponents, due to the underrepresentation of various functional units (such as the cytosine-guanine CG combination and its complementary GC complex). Retaining higher-order correlations in the construction of the tensor products is a way to approach closer the structure of the multifractal spectra of the actual genomic sequences. This hierarchical approach is generic and is applicable to other correlated symbolic sequences.

Intergenic, gene terminal, and intragenic CpG islands in the human genome

Background

Recently, it has been discovered that the human genome contains many transcription start sites for non-coding RNA. Regulatory regions related to transcription of this non-coding RNAs are poorly studied. Some of these regulatory regions may be associated with CpG islands located far from transcription start-sites of any protein coding gene. The human genome contains many such CpG islands; however, until now their properties were not systematically studied.

Results

We studied CpG islands located in different regions of the human genome using methods of bioinformatics and comparative genomics. We have observed that CpG islands have a preference to overlap with exons, including exons located far from transcription start site, but usually extend well into introns. Synonymous substitution rate of CpG-containing codons becomes substantially reduced in regions where CpG islands overlap with protein-coding exons, even if they are located far downstream from transcription start site. CAGE tag analysis displayed frequent transcription start sites in all CpG islands, including those found far from transcription start sites of protein coding genes. Computational prediction and analysis of published ChIP-chip data revealed that CpG islands contain an increased number of sites recognized by Sp1 protein. CpG islands containing more CAGE tags usually also contain more Sp1 binding sites. This is especially relevant for CpG islands located in 3' gene regions. Various examples of transcription, confirmed by mRNAs or ESTs, but with no evidence of protein coding genes, were found in CAGE-enriched CpG islands located far from transcription start site of any known protein coding gene.

Conclusions

CpG islands located far from transcription start sites of protein coding genes have transcription initiation activity and display Sp1 binding properties. In exons, overlapping with these islands, the synonymous substitution rate of CpG containing codons is decreased. This suggests that these CpG islands are involved in transcription initiation, possibly of some non-coding RNAs.

CpG_MI: a novel approach for identifying functional CpG islands in mammalian genomes

CpG islands (CGIs) are CpG-rich regions compared to CpG-depleted bulk DNA of mammalian genomes and are generally regarded as the epigenetic regulatory regions in association with unmethylation, promoter activity and histone modifications. Accurate identification of CpG islands with epigenetic regulatory function in bulk genomes is of wide interest. Here, the common features of functional CGIs are identified using an average mutual information method to differentiate functional CGIs from the remaining CGIs. A new approach (CpG mutual information, CpG_MI) was further explored to identify functional CGIs based on the cumulative mutual information of physical distances between two neighboring CpGs. Compared to current approaches, CpG_MI achieved the highest prediction accuracy. This approach also identified new functional CGIs overlapping with gene promoter regions which were missed by other algorithms. Nearly all CGIs identified by CpG_MI overlapped with histone modification marks. CpG_MI could also be used to identify potential functional CGIs in other mammalian genomes, as the CpG dinucleotide contents and cumulative mutual information distributions are almost the same among six mammalian genomes in our analysis. It is a reliable quantitative tool for the identification of functional CGIs from bulk genomes and helps in understanding the relationships between genomic functional elements and epigenomic modifications.

Gene conversion causing human inherited disease: evidence for involvement of non-B-DNA-forming sequences and recombination-promoting motifs in DNA breakage and repair

A variety of DNA sequence motifs including inverted repeats, minisatellites, and the chi recombination hotspot, have been reported in association with gene conversion in human genes causing inherited disease. However, no methodical statistically based analysis has been performed to formalize these observations. We have performed an in silico analysis of the DNA sequence tracts involved in 27 nonoverlapping gene conversion events in 19 different genes reported in the context of inherited disease. We found that gene conversion events tend to occur within (C+G)- and CpG-rich regions and that sequences with the potential to form non-B-DNA structures, and which may be involved in the generation of double-strand breaks that could, in turn, serve to promote gene conversion, occur disproportionately within maximal converted tracts and/or short flanking regions. Maximal converted tracts were also found to be enriched (P<0.01) in a truncated version of the chi-element (a TGGTGG motif), immunoglobulin heavy chain class switch repeats, translin target sites and several novel motifs including (or overlapping) the classical meiotic recombination hotspot, CCTCCCCT. Finally, gene conversions tend to occur in genomic regions that have the potential to fold into stable hairpin conformations. These findings support the concept that recombination-inducing motifs, in association with alternative DNA conformations, can promote recombination in the human genome.

Contrast features of CpG islands in the promoter and other regions in the dog genome

The recent release of the domestic dog genome provides us with an ideal opportunity to investigate dog-specific genomic features. In this study, we performed a systematic analysis of CpG islands (CGIs), which are often considered gene markers, in the dog genome. Relative to the human and mouse genomes, the dog genome has a remarkably large number of CGIs and high CGI density, which is contributed by its noncoding sequences. Surprisingly, the dog genome has fewer CGIs associated with the promoter regions of genes than the human or the mouse. Further examination of functional features of dog-human-mouse homologous genes suggests that the dog might have undergone a faster erosion rate of promoter-associated CGIs than the human or mouse. Some genetic or genomic factors such as local recombination rate and karyotype may be related to the unique dog CGI features.

A database of vertebrate longevity records and their relation to other life-history traits

Longevity is a major characteristic of animals that has long fascinated scientists. In this work, we present a comprehensive database of animal longevity records and related life-history traits entitled AnAge, which we compiled and manually curated from an extensive literature. AnAge started as a collection of longevity records, but has since been expanded to include quantitative data for numerous other life-history traits, including body masses at different developmental stages, reproductive data such as age at sexual maturity and measurements of reproductive output, and physiological traits related to metabolism. AnAge features over 4000 vertebrate species and is a central resource for applying the comparative method to studies of longevity and life-history evolution across the tree of life. Moreover, by providing a reference value for longevity and other life-history traits, AnAge can prove valuable to a broad range of biologists working in evolutionary biology, ecology, zoology, physiology and conservation biology. AnAge is freely available online (http://genomics.senescence.info/species/).

Transcription initiation activity sets replication origin efficiency in mammalian cells

Genomic mapping of DNA replication origins (ORIs) in mammals provides a powerful means for understanding the regulatory complexity of our genome. Here we combine a genome-wide approach to identify preferential sites of DNA replication initiation at 0.4% of the mouse genome with detailed molecular analysis at distinct classes of ORIs according to their location relative to the genes. Our study reveals that 85% of the replication initiation sites in mouse embryonic stem (ES) cells are associated with transcriptional units. Nearly half of the identified ORIs map at promoter regions and, interestingly, ORI density strongly correlates with promoter density, reflecting the coordinated organisation of replication and transcription in the mouse genome. Detailed analysis of ORI activity showed that CpG island promoter-ORIs are the most efficient ORIs in ES cells and both ORI specification and firing efficiency are maintained across cell types. Remarkably, the distribution of replication initiation sites at promoter-ORIs exactly parallels that of transcription start sites (TSS), suggesting a co-evolution of the regulatory regions driving replication and transcription. Moreover, we found that promoter-ORIs are significantly enriched in CAGE tags derived from early embryos relative to all promoters. This association implies that transcription initiation early in development sets the probability of ORI activation, unveiling a new hallmark in ORI efficiency regulation in mammalian cells.

The duplication of the genetic information of a cell starts from specific sites on the chromosomes called DNA replication origins. Their number varies from a few hundred in yeast cells to several thousands in human cells, distributed along the genome at comparable distances in both systems. An important question in the field is to understand how origins of replication are specified and regulated in the mammalian genome, as neither their location nor their activity can be directly inferred from the DNA sequence. Previous studies at individual origins and, more recently, at large scale across 1% of the human genome, have revealed that most origins overlap with transcriptional regulatory elements, and specifically with gene promoters. To gain insight into the nature of the relationship between active transcription and origin specification we have combined a genomic mapping of origins at 0.4% of the mouse genome with detailed studies of activation efficiency. The data identify two types of origins with distinct regulatory properties: highly efficient origins map at CpG island-promoters and low efficient origins locate elsewhere in association with transcriptional units. We also find a remarkable parallel organisation of the replication initiation sites and transcription start sites at efficient promoter-origins that suggests a prominent role of transcription initiation in setting the efficiency of replication origin activation.

CpG islands: algorithms and applications in methylation studies

Methylation occurs frequently at 5'-cytosine of the CpG dinucleotides in vertebrate genomes; however, this epigenetic feature is rarely observed in CpG islands (CGIs) or CpG clusters in the promoter regions of genes. Aberrant methylation of the promoter-associated CGIs might influence gene expression and cause carcinogenesis. Because of the functional importance, multiple algorithms have been available for identifying CGIs in a genome or a sequence. They can be categorized into the traditional algorithms (e.g., Gardiner-Garden and Frommer (1987), Takai and Jones (2002), and CpGPRoD (2002)) or statistical property based algorithms (CpGcluster (2006) and CG cluster (2007)). We reviewed the features of these algorithms and evaluated their performance on identifying functional CGIs using genome-wide methylation data. Moreover, identification of CGIs is an initial step in many recent studies for predicting methylation status as well as in the design of methylation detection platforms. We reviewed the benchmarks and features used in these studies.

CpG islands or CpG clusters: how to identify functional GC-rich regions in a genome?

BACKGROUND

CpG islands (CGIs), clusters of CpG dinucleotides in GC-rich regions, are often located in the 5' end of genes and considered gene markers. Hackenberg et al. (2006) recently developed a new algorithm, CpGcluster, which uses a completely different mathematical approach from previous traditional algorithms. Their evaluation suggests that CpGcluster provides a much more efficient approach to detecting functional clusters or islands of CpGs.

RESULTS

We systematically compared CpGcluster with the traditional algorithm by Takai and Jones (2002). Our comparisons of (1) the number of islands versus the number of genes in a genome, (2) the distribution of islands in different genomic regions, (3) island length, (4) the distance between two neighboring islands, and (5) methylation status suggest that Takai and Jones' algorithm is overall more appropriate for identifying promoter-associated islands of CpGs in vertebrate genomes.

CONCLUSION

The generation of genome sequence and DNA methylation data is expected to accelerate greatly. The information in this study is important for its extensive utility in gene feature analysis and epigenomics including gene prediction and methylation chip design in different genomes.

SNP-based prediction of the human germ cell methylation landscape

Base substitution occurs at a high rate at CpG dinucleotides due to the frequent methylation of CpG and the deamination of methylated cytosine to thymine. If these substitutions occur in germ cells, they constitute a heritable mutation that may eventually rise to polymorphic frequencies, hence resulting in a SNP that is methylation associated. In this study, we sought to identify clusters of methylation associated SNPs as a basis for prediction of methylation landscapes of germ cell genomes. Genomic regions enriched with methylation associated SNPs, namely "methylation associated SNP clusters", were identified with an agglomerative hierarchical clustering algorithm. Repetitive elements, segmental duplications, and syntenic tandem DNA repeats were enriched in methylation associated SNP clusters. The frequency of methylation associated SNPs in Alu Y/S elements exhibited a gradient pattern suggestive of linear spreading, being higher in proximity to methylation associated SNP clusters and lower closer to CpG islands. Interestingly, methylation associated SNP clusters were over-represented near the transcriptional initiation sites of immune response genes. We propose a de novo DNA methylation model during germ cell development whereby a pattern is established by long-range chromatic interactions through syntenic repeats combined with regional methylation spreading from methylation associated SNP clusters.

Divergence of the vertebrate sp1A/ryanodine receptor domain and SOCS box-containing (Spsb) gene family and its expression and regulation within the mouse brain

The Spsb family of genes encode well-conserved proteins of unknown function. Mammalian Spsb genes are likely the result of three separate duplication and divergence events during vertebrate evolution. The phylogenetic relationship along with expression and regulation of Spsb genes may offer insight into the evolution and function of this gene family in vertebrates. We have established that Spsb genes are expressed in numerous tissues, however their pattern and level of expression is tissue-dependent. Further, only Spsb1 is responsive to stress caused by ethanol exposure in the mouse brain, which suggests that Spsb genes have acquired different regulatory mechanisms. Analysis of cis-regulatory elements supports this, but also reveals some common regulatory modules involved in cell proliferation and stress response. Our results contribute to the growing body of data on the expression and function of Spsb genes, which serve as a model for studies on the origin, divergence and specialization of eukaryotic gene families.