Statistical methods for the analysis of DNA sequences: application to dinucleotide distribution in the human genome

Questa tesi si inserisce nell'ambito delle analisi statistiche e dei metodi stocastici applicati all'analisi delle sequenze di DNA. Nello specifico il nostro lavoro è incentrato sullo studio del dinucleotide CG (CpG) all'interno del genoma umano, che si trova raggruppato in zone specifiche denominate CpG islands. Queste sono legate alla metilazione del DNA, un processo che riveste un ruolo fondamentale nella regolazione genica. La prima parte dello studio è dedicata a una caratterizzazione globale del contenuto e della distribuzione dei 16 diversi dinucleotidi all'interno del genoma umano: in particolare viene studiata la distribuzione delle distanze tra occorrenze successive dello stesso dinucleotide lungo la sequenza. I risultati vengono confrontati con diversi modelli nulli: sequenze random generate con catene di Markov di ordine zero (basate sulle frequenze relative dei nucleotidi) e uno (basate sulle probabilità di transizione tra diversi nucleotidi) e la distribuzione geometrica per le distanze. Da questa analisi le proprietà caratteristiche del dinucleotide CpG emergono chiaramente, sia dal confronto con gli altri dinucleotidi che con i modelli random. A seguito di questa prima parte abbiamo scelto di concentrare le successive analisi in zone di interesse biologico, studiando l’abbondanza e la distribuzione di CpG al loro interno (CpG islands, promotori e Lamina Associated Domains). Nei primi due casi si osserva un forte arricchimento nel contenuto di CpG, e la distribuzione delle distanze è spostata verso valori inferiori, indicando che questo dinucleotide è clusterizzato. All’interno delle LADs si trovano mediamente meno CpG e questi presentano distanze maggiori. Infine abbiamo adottato una rappresentazione a random walk del DNA, costruita in base al posizionamento dei dinucleotidi: il walk ottenuto presenta caratteristiche drasticamente diverse all’interno e all’esterno di zone annotate come CpG island. Riteniamo pertanto che metodi basati su questo approccio potrebbero essere sfruttati per migliorare l’individuazione di queste aree di interesse nel genoma umano e di altri organismi.

Preferential isolation of DNA fragments associated with CpG islands.

We describe a procedure for preferential isolation of DNA fragments with G+C-rich portions. Such fragments occur in known genes within or adjacent to CpG islands. Since about 56% of human genes are associated with CpG islands, isolation of these fragments permits detection and probing of many genes within much larger segments of DNA, such as cosmids or yeast artificial chromosomes, which have not been sequenced. Cloned DNA fragments digested with four restriction endonucleases were subjected to denaturing gradient gel electrophoresis. Long G+C-rich sections in fragments inhibit strand dissociation after the fragments reach retardation level in the gradient; such fragments are retained in the gel after most others disappear. Nucleotide sequences of the retained fragments show that about half of these fragments appear to be derived from CpG islands. Northern analysis indicated the presence of RNA complementary to most of the retained fragments. A heuristic approach to the relation between base sequence and the kinetics of strand dissociation of partly melted molecules appears to account for retention and nonretention. The expectation that CpG island fragments will be enriched among fragments retained in a denaturing gradient is supported by rate estimates based on melting theory applied to known sequences. This method, designated SPM for segregation of partly melted molecules, is expected to provide a means for convenient and efficient isolation of genes from unsequenced DNA.

The properties of CpG islands in the putative promoter regions of human immunoglobulin (Ig) genes

CpG island is a GC-rich motif occurred in gene promoter region, which can play important roles in gene silencing and imprinting. Here, we present a set of discriminant functions that can recognize the structural and compositional features of CpG islands in the putative promoter regions (PPRs) of human and mouse immunoglobulin (Ig) genes. We showed that the PPRs of both human and mouse Ig genes irrespective of gene chromosomal localization are apparently CpG island poor, with a low percentage of the CpG islands overlapped with the transcription start site (TSS). The human Ig genes that have CpG islands in the PPRs show a very narrow range of CpG densities. 47% of the Ig genes fall in the range of 3.5-4 CpGs/100 bp. In contrast, the non-Ig genes examined have a wide range of the density of CpG island, with 10.5% having the density of 8.1-15 CpGs/100 bp. Meantime, five patterns of the CpG distributions within the CpG islands have been classified: Pat A, B, C, D, and E. 21.6% and 10.8% of the Ig genes fall into the Pat B and Pat D groups, respectively, which were significantly higher than the non-Ig genes examined (8.2% and 3.8%). Moreover, the length of CpG islands is shorter in human Ig genes than in non-Ig genes but is much longer than in mouse orthologues. These findings provide a clear picture of non-neutral and nonrandom occurrence of the CpG islands in the PPRs of human and mouse Ig genes, which facilitate rational recommendations regarding their nomenclature. (C) 2005 Elsevier B.V. All rights reserved.

Features and trend of loss of promoter-associated CpG islands in the human and mouse genomes

CpG islands (CGIs) are often considered as gene markers, but the number of CGIs varies among mammalian genomes that have similar numbers of genes. In this study, we investigated the distribution of CGIs in the promoter regions of 3,197 human-mouse ortholo

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

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 diffe

DNA methylation subgroups and the CpG island methylator phenotype in gastric cancer: a comprehensive profiling approach

BACKGROUND: Methylation-induced silencing of promoter CpG islands in tumor suppressor genes plays an important role in human carcinogenesis. In colorectal cancer, the CpG island methylator phenotype (CIMP) is defined as widespread and elevated levels of DNA methylation and CIMP+ tumors have distinctive clinicopathological and molecular features. In contrast, the existence of a comparable CIMP subtype in gastric cancer (GC) has not been clearly established. To further investigate this issue, in the present study we performed comprehensive DNA methylation profiling of a well-characterised series of primary GC.

METHODS: The methylation status of 1,421 autosomal CpG sites located within 768 cancer-related genes was investigated using the Illumina GoldenGate Methylation Panel I assay on DNA extracted from 60 gastric tumors and matched tumor-adjacent gastric tissue pairs. Methylation data was analysed using a recursively partitioned mixture model and investigated for associations with clinicopathological and molecular features including age, Helicobacter pylori status, tumor site, patient survival, microsatellite instability and BRAF and KRAS mutations.

RESULTS: A total of 147 genes were differentially methylated between tumor and matched tumor-adjacent gastric tissue, with HOXA5 and hedgehog signalling being the top-ranked gene and signalling pathway, respectively. Unsupervised clustering of methylation data revealed the existence of 6 subgroups under two main clusters, referred to as L (low methylation; 28% of cases) and H (high methylation; 72%). Female patients were over-represented in the H tumor group compared to L group (36% vs 6%; P = 0.024), however no other significant differences in clinicopathological or molecular features were apparent. CpG sites that were hypermethylated in group H were more frequently located in CpG islands and marked for polycomb occupancy.

CONCLUSIONS: High-throughput methylation analysis implicates genes involved in embryonic development and hedgehog signaling in gastric tumorigenesis. GC is comprised of two major methylation subtypes, with the highly methylated group showing some features consistent with a CpG island methylator phenotype.

Detecting recombination in 4-taxa DNA sequence alignments with Bayesian hidden Markov models and Markov chain Monte Carlo

This article presents a statistical method for detecting recombination in DNA sequence alignments, which is based on combining two probabilistic graphical models: (1) a taxon graph (phylogenetic tree) representing the relationship between the taxa, and (2) a site graph (hidden Markov model) representing interactions between different sites in the DNA sequence alignments. We adopt a Bayesian approach and sample the parameters of the model from the posterior distribution with Markov chain Monte Carlo, using a Metropolis-Hastings and Gibbs-within-Gibbs scheme. The proposed method is tested on various synthetic and real-world DNA sequence alignments, and we compare its performance with the established detection methods RECPARS, PLATO, and TOPAL, as well as with two alternative parameter estimation schemes.

MED1, a novel human methyl-CpG-binding endonuclease, interacts with DNA mismatch repair protein MLH1

The DNA mismatch repair (MMR) is a specialized system, highly conserved throughout evolution, involved in the maintenance of genomic integrity. To identify novel human genes that may function in MMR, we employed the yeast interaction trap. Using the MMR protein MLH1 as bait, we cloned MED1. The MED1 protein forms a complex with MLH1, binds to methyl-CpG-containing DNA, has homology to bacterial DNA repair glycosylases/lyases, and displays endonuclease activity. Transfection of a MED1 mutant lacking the methyl-CpG-binding domain (MBD) is associated with microsatellite instability (MSI). These findings suggest that MED1 is a novel human DNA repair protein that may be involved in MMR and, as such, may be a candidate eukaryotic homologue of the bacterial MMR endonuclease, MutH. In addition, these results suggest that cytosine methylation may play a role in human DNA repair.

Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands.

Precise mapping of DNA methylation patterns in CpG islands has become essential for understanding diverse biological processes such as the regulation of imprinted genes, X chromosome inactivation, and tumor suppressor gene silencing in human cancer. We describe a new method, MSP (methylation-specific PCR), which can rapidly assess the methylation status of virtually any group of CpG sites within a CpG island, independent of the use of methylation-sensitive restriction enzymes. This assay entails initial modification of DNA by sodium bisulfite, converting all unmethylated, but not methylated, cytosines to uracil, and subsequent amplification with primers specific for methylated versus unmethylated DNA. MSP requires only small quantities of DNA, is sensitive to 0.1% methylated alleles of a given CpG island locus, and can be performed on DNA extracted from paraffin-embedded samples. MSP eliminates the false positive results inherent to previous PCR-based approaches which relied on differential restriction enzyme cleavage to distinguish methylated from unmethylated DNA. In this study, we demonstrate the use of MSP to identify promoter region hypermethylation changes associated with transcriptional inactivation in four important tumor suppressor genes (p16, p15, E-cadherin, and von Hippel-Lindau) in human cancer.

Tiggers and DNA transposon fossils in the human genome.

We report several classes of human interspersed repeats that resemble fossils of DNA transposons, elements that move by excision and reintegration in the genome, whereas previously characterized mammalian repeats all appear to have accumulated by retrotransposition, which involves an RNA intermediate. The human genome contains at least 14 families and > 100,000 degenerate copies of short (180-1200 bp) elements that have 14- to 25-bp terminal inverted repeats and are flanked by either 8 bp or TA target site duplications. We describe two ancient 2.5-kb elements with coding capacity, Tigger1 and -2, that closely resemble pogo, a DNA transposon in Drosophila, and probably were responsible for the distribution of some of the short elements. The deduced pogo and Tigger proteins are related to products of five DNA transposons found in fungi and nematodes, and more distantly, to the Tc1 and mariner transposases. They also are very similar to the major mammalian centromere protein CENP-B, suggesting that this may have a transposase origin. We further identified relatively low-copy-number mariner elements in both human and sheep DNA. These belong to two subfamilies previously identified in insect genomes, suggesting lateral transfer between diverse species.

Methylenetetrahydrofolate reductase CpG islands: Epigenotyping

Background Determination of the differential DNA methylation patterns of methylenetetrahydrofolate reductase (MTHFR) that are associated with differential MTHFR activity is important to understand the pathogenesis of ischemic stroke. However, to date, no data are available on the differential DNA methylation profiles of Kelantanese Malays. Therefore, we developed a rapid and efficient serial pyrosequencing assay to determine differential DNA methylation profiles of MTHFR, which help to further our understanding of the pathogenesis of ischemic stroke. The developed assay also served as the validation platform for our previous computational epigenetic research on MTHFR. Methods Polymerase chain reaction primers were designed and validated to specifically amplify the cytosine that is followed by guanine residues (CpGs) A and B regions. Prior epigenotyping on 110 Kelantanese Malays, the serial pyrosequencing assays for the CpGs A and B regions were validated using five validation controls. The mean values of the DNA methylation profiles of CpGs A and B were calculated. Results The mean DNA methylation levels for CpGs A and B were 0.984 ± 0.582 and 2.456 ± 1.406, respectively. The CpGs 8 and 20 showed the highest (5.581 ± 4.497) and the lowest (0.414 ± 2.814) levels of DNA methylation at a single-base resolution. Conclusion We have successfully developed and validated a pyrosequencing assay that is fast and can yield high-quality pyrograms for DNA methylation analysis and is therefore applicable to high throughput study. Using this newly developed pyrosequencing assay, the MTHFR DNA methylation profiles of 110 Kelantanese Malays were successfully determined. It also validated our computational epigenetic research on MTHFR.

DNA methylation analysis of phenotype specific stratified Indian population

Background: DNA methylation and its perturbations are an established attribute to a wide spectrum of phenotypic variations and disease conditions. Indian traditional system practices personalized medicine through indigenous concept of distinctly descriptive physiological, psychological and anatomical features known as prakriti. Here we attempted to establish DNA methylation differences in these three prakriti phenotypes. Methods: Following structured and objective measurement of 3416 subjects, whole blood DNA of 147 healthy male individuals belonging to defined prakriti (Vata, Pitta and Kapha) between the age group of 20-30years were subjected to methylated DNA immunoprecipitation (MeDIP) and microarray analysis. After data analysis, prakriti specific signatures were validated through bisulfite DNA sequencing. Results: Differentially methylated regions in CpG islands and shores were significantly enriched in promoters/UTRs and gene body regions. Phenotypes characterized by higher metabolism (Pitta prakriti) in individuals showed distinct promoter (34) and gene body methylation (204), followed by Vata prakriti which correlates to motion showed DNA methylation in 52 promoters and 139 CpG islands and finally individuals with structural attributes (Kapha prakriti) with 23 and 19 promoters and CpG islands respectively. Bisulfite DNA sequencing of prakriti specific multiple CpG sites in promoters and 5'-UTR such as; LHX1 (Vata prakriti), SOX11 (Pitta prakriti) and CDH22 (Kapha prakriti) were validated. Kapha prakriti specific CDH22 5'-UTR CpG methylation was also found to be associated with higher body mass index (BMI). Conclusion: Differential DNA methylation signatures in three distinct prakriti phenotypes demonstrate the epigenetic basis of Indian traditional human classification which may have relevance to personalized medicine.

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

Sequence context analysis of 8.2 million single nucleotide polymorphisms in the human genome

We analyzed n-mers (n=3-8) in the local environment of 8,249,446 human SNPs and compared their distribution with that in the genome reference sequences. The results revealed that the short sequences, which contained at least one CpG dinucleotide, occurred