Nutritional factors and gender influence age-related DNA methylation in the human rectal mucosa

Aberrant methylation of CpG islands (CGI) occurs in many genes expressed in colonic epithelial cells, and may contribute to the dysregulation of signalling pathways associated with carcinogenesis. This cross-sectional study assessed the relative importance of age, nutritional exposures and other environmental factors in the development of CGI methylation. Rectal biopsies were obtained from 185 individuals (84 male, 101 female) shown to be free of colorectal disease, and for whom measurements of age, body size, nutritional status and blood cell counts were available. We used quantitative DNA methylation analysis combined with multivariate modelling to investigate the relationships between nutritional, anthropometric and metabolic factors and the CGI methylation of 11 genes, together with LINE-1 as an index of global DNA methylation. Age was a consistent predictor of CGI methylation for 9/11 genes but significant positive associations with folate status and negative associations with vitamin D and selenium status were also identified for several genes. There was evidence for positive associations with blood monocyte levels and anthropometric factors for some genes. In general, CGI methylation was higher in males than in females and differential effects of age and other factors on methylation in males and females were identified. In conclusion, levels of age-related CGI methylation in the healthy human rectal mucosa are influenced by gender, the availability of folate, vitamin D and selenium, and perhaps by factors related to systemic inflammation

PREDICTION OF DNA METHYLATION BASED ON GENOMIC ARCHITECTURE AND APPLICATIONS OF POSITIONAL WEIGHT MATRICES

Gene silencing due to epigenetic mechanisms shows evidence of significant contributions to cancer development. We hypothesis that the genetic architecture based on retrotransposon elements surrounding the transcription start site, plays an important role in the suppression and promotion of DNA methylation. In our investigation we found a high rate of SINE and LINEs retrotransposon elements near the transcription start site of unmethylated genes when compared to methylated genes. The presence of these elements were positively associated with promoter methylation, contrary to logical expectations, due to the malicious effects of retrotransposon elements which insert themselves randomly into the genome causing possible loss of gene function. In our genome wide analysis of human genes, results suggested that 22% of the genes in cancer were predicted to be methylation-prone; in cancer these genes are generally down-regulated and function in the development process. In summary, our investigation validated our hypothesis and showed that these widespread genomic elements in cancer are highly associated with promoter DNA methylation and may further participate in influencing epigenetic regulation.

IDENTIFICATION OF FACTORS INVOLVED IN DNA METHYLATION OF CPG-ISLAND-PROMOTERS

Repression of many tumor suppressor genes (TSGs) in cancer is mediated by aberrantly increased DNA methylation levels at promoter CpG islands (CGI). About one-fourth of empirically defined human promoters are surrounded by or contain clustered repetitive elements. It was previously observed that a sharp transition of methylation occurs between highly methylated repetitive elements (SINE or LINE) and unmethylated CGI-promoters (e.g. P16, VHL, CDH and RIL) in normal tissues. The functions that lead to increased CGI methylation in cancer remain poorly understood. We propose that CGI-promoters contain cis-elements for triggering de novo DNA methylation. In the first part of our project, we established a site-specific integration system with enforced local transcriptional repression in colorectal cancer cells and monitored the occurrence of de novo DNA methylation in exogenous fragments containing a CGI-promoter and repetitive elements. Initial de novo methylation was seeded at specific CG sites in a repetitive element, and accelerated by persistent binding of a KRAB-containing transcriptional repressor. Furthermore, additional repetitive elements (LINE and SINE) located adjacent to the promoter could confer DNA methylation spreading into the CGI particularly in the setting of KRAB-factor binding. However, a repressive chromatin alone was not sufficient to initiate DNA methylation, which required specific DNA sequences and was integration-site (and/or cell-line) specific. In addition, all the methylation observed showed slow and gradual accumulation over several months of culture. Overall, these results demonstrate a requirement for specific DNA sequences to trigger de novo DNA methylation, and repetitive elements as cis-regulatory factors to cooperate with strengthened transcriptional repression in promoting methylation spreading. In the second part, we re-introduced disrupted DNMT3B or DNMT1 into HCT116 DKO cells and mapped the remethylation pattern through a profiling method (DREAM). Moderate remethylation occurred when DNMT3B was re-expressed with a preference toward non-CGI and non-promoter regions. Hence, there exists a set of genomic regions with priority to be targets for DNMT3B in somatic cells.

Selective association of the methyl-CpG binding protein MBD2 with the silent p14/p16 locus in human neoplasia

DNA methylation of tumor suppressor genes is a common feature of human cancer. The cyclin-dependent kinase inhibitor gene p16/Ink4A is hypermethylated in a wide range of malignant tissues and the p14/ARF gene located 20 kb upstream on chromosome 9p21 is also methylated in carcinomas. p14/ARF (ARF, alternative reading frame) does not inhibit the activities of cyclins or cyclin-dependent kinase complexes; however, the importance of the two gene products in the etiology of cancer resides in their involvement in two major cell cycle regulatory pathways: p53 and the retinoblastoma protein, Rb, respectively. Distinct first exons driven from separate promoters are spliced onto the common exons 2 and 3 and the resulting proteins are translated in different reading frames. Both genes are expressed in normal cells but can be alternatively or coordinately silenced when their CpG islands are hypermethylated. Herein, we examined the presence of methyl-CpG binding proteins associated with aberrantly methylated promoters, the distribution of acetylated histones H3 and H4 by chromatin immunoprecipitation assays, and the effect of chemical treatment with 5-aza-2′-deoxycytidine (5aza-dC) and trichostatin A on gene induction in colon cell lines by quantitative reverse transcriptase–PCR. We observed that the methyl-CpG binding protein MBD2 is targeted to methylated regulatory regions and excludes the acetylated histones H3 and H4, resulting in a localized inactive chromatin configuration. When methylated, the genes can be induced by 5aza-dC but the combined action of 5aza-dC and trichostatin A results in robust gene expression. Thus, methyl-CpG binding proteins and histone deacetylases appear to cooperate in vivo, with a dominant effect of DNA methylation toward histone acetylation, and repress expression of tumor suppressor genes hypermethylated in cancers.

Computational analysis of base composition pattern and promoter elements in the putative promoter regions in relation to expression profiles of 682 human genes on chromosome 22

The base composition pattern (BCP) in the putative promoter region (PPRs) up to 5 Kb lengths of 682 human genes on Chromosome 22 (Chr22) was examined. Two-dimensional (2D) and three-dimensional (3D) functions were designed to delineate the DNA base composition, with four major patterns identified. It is found that 17.6% genes include TATA box, 28.0% GC box, 18.9% CAAT box and 38.4% CpG islands, and approximately 10% genes have one of four putative initiator (Inr) motifs. The occurrence of the promoter elements is tightly associated with the base composition features in the promoter regions, and the associations of the base composition features with occurrence of the promoter elements in the promoter regions mediate tissue-wide expression of the genes in human. The occurrence of two or more promoter elements in the promoter regions is required for the medium- and wide-range expression profiles of the human genes on Chr22. Thus, the reported data shed light on the characteristics of the PPRs of the human genes on Chr22, which may improve our understanding of regulatory roles of the PPRs with occurrence of the promoter elements in gene expression.

Integration pattern of HIV-1 based lentiviral vector carrying recombinant coagulation factor VIII in Sk-Hep and 293T cells

293T and Sk-Hep-1 cells were transduced with a replication-defective self-inactivating HIV-1 derived vector carrying FVIII cDNA. The genomic DNA was sequenced to reveal LTR/human genome junctions and integration sites. One hundred and thirty-two sequences matched human sequences, with an identity of at least 98%. The integration sites in 293T-FVIIIDB and in Sk-Hep-FVIIIDB cells were preferentially located in gene regions. The integrations in both cell lines were distant from the CpG islands and from the transcription start sites. A comparison between the two cell lines showed that the lentiviral-transduced DNA had the same preferred regions in the two different cell lines.

Phenotypic consequences of copy number variation: insights from Smith-Magenis and Potocki-Lupski syndrome mouse models.

A large fraction of genome variation between individuals is comprised of submicroscopic copy number variation of genomic DNA segments. We assessed the relative contribution of structural changes and gene dosage alterations on phenotypic outcomes with mouse models of Smith-Magenis and Potocki-Lupski syndromes. We phenotyped mice with 1n (Deletion/+), 2n (+/+), 3n (Duplication/+), and balanced 2n compound heterozygous (Deletion/Duplication) copies of the same region. Parallel to the observations made in humans, such variation in gene copy number was sufficient to generate phenotypic consequences: in a number of cases diametrically opposing phenotypes were associated with gain versus loss of gene content. Surprisingly, some neurobehavioral traits were not rescued by restoration of the normal gene copy number. Transcriptome profiling showed that a highly significant propensity of transcriptional changes map to the engineered interval in the five assessed tissues. A statistically significant overrepresentation of the genes mapping to the entire length of the engineered chromosome was also found in the top-ranked differentially expressed genes in the mice containing rearranged chromosomes, regardless of the nature of the rearrangement, an observation robust across different cell lineages of the central nervous system. Our data indicate that a structural change at a given position of the human genome may affect not only locus and adjacent gene expression but also "genome regulation." Furthermore, structural change can cause the same perturbation in particular pathways regardless of gene dosage. Thus, the presence of a genomic structural change, as well as gene dosage imbalance, contributes to the ultimate phenotype.

Nineteen additional unpredicted transcripts from human chromosome 21.

The identification of all human chromosome 21 (HC21) genes is a necessary step in understanding the molecular pathogenesis of trisomy 21 (Down syndrome). The first analysis of the sequence of 21q included 127 previously characterized genes and predicted an additional 98 novel anonymous genes. Recently we evaluated the quality of this annotation by characterizing a set of HC21 open reading frames (C21orfs) identified by mapping spliced expressed sequence tags (ESTs) and predicted genes (PREDs), identified only in silico. This study underscored the limitations of in silico-only gene prediction, as many PREDs were incorrectly predicted. To refine the HC21 annotation, we have developed a reliable algorithm to extract and stringently map sequences that contain bona fide 3' transcript ends to the genome. We then created a specific 21q graphical display allowing an integrated view of the data that incorporates new ESTs as well as features such as CpG islands, repeats, and gene predictions. Using these tools we identified 27 new putative genes. To validate these, we sequenced previously cloned cDNAs and carried out RT-PCR, 5'- and 3'-RACE procedures, and comparative mapping. These approaches substantiated 19 new transcripts, thus increasing the HC21 gene count by 9.5%. These transcripts were likely not previously identified because they are small and encode small proteins. We also identified four transcriptional units that are spliced but contain no obvious open reading frame. The HC21 data presented here further emphasize that current gene prediction algorithms miss a substantial number of transcripts that nevertheless can be identified using a combination of experimental approaches and multiple refined algorithms.

Epigenetic features of human mesenchymal stem cells determine their permissiveness for induction of relevant transcriptional changes by SYT-SSX1.

BACKGROUND: A characteristic SYT-SSX fusion gene resulting from the chromosomal translocation t(X;18)(p11;q11) is detectable in almost all synovial sarcomas, a malignant soft tissue tumor widely believed to originate from as yet unidentified pluripotent stem cells. The resulting fusion protein has no DNA binding motifs but possesses protein-protein interaction domains that are believed to mediate association with chromatin remodeling complexes. Despite recent advances in the identification of molecules that interact with SYT-SSX and with the corresponding wild type SYT and SSX proteins, the mechanisms whereby the SYT-SSX might contribute to neoplastic transformation remain unclear. Epigenetic deregulation has been suggested to be one possible mechanism. METHODOLOGY/PRINCIPAL FINDINGS: We addressed the effect of SYT/SSX expression on the transcriptome of four independent isolates of primary human bone marrow mesenchymal stem cells (hMSC). We observed transcriptional changes similar to the gene expression signature of synovial sarcoma, principally involving genes whose regulation is linked to epigenetic factors, including imprinted genes, genes with transcription start sites within a CpG island and chromatin related genes. Single population analysis revealed hMSC isolate-specific transcriptional changes involving genes that are important for biological functions of stem cells as well as genes that are considered to be molecular markers of synovial sarcoma including IGF2, EPHRINS, and BCL2. Methylation status analysis of sequences at the H19/IGF2 imprinted locus indicated that distinct epigenetic features characterize hMSC populations and condition the transcriptional effects of SYT-SSX expression. CONCLUSIONS/SIGNIFICANCE: Our observations suggest that epigenetic features may define the cellular microenvironment in which SYT-SSX displays its functional effects.

Ultra high throughput sequencing in human DNA variation detection: a comparative study on the NDUFA3-PRPF31 region.

BACKGROUND: Ultra high throughput sequencing (UHTS) technologies find an important application in targeted resequencing of candidate genes or of genomic intervals from genetic association studies. Despite the extraordinary power of these new methods, they are still rarely used in routine analysis of human genomic variants, in part because of the absence of specific standard procedures. The aim of this work is to provide human molecular geneticists with a tool to evaluate the best UHTS methodology for efficiently detecting DNA changes, from common SNPs to rare mutations. METHODOLOGY/PRINCIPAL FINDINGS: We tested the three most widespread UHTS platforms (Roche/454 GS FLX Titanium, Illumina/Solexa Genome Analyzer II and Applied Biosystems/SOLiD System 3) on a well-studied region of the human genome containing many polymorphisms and a very rare heterozygous mutation located within an intronic repetitive DNA element. We identify the qualities and the limitations of each platform and describe some peculiarities of UHTS in resequencing projects. CONCLUSIONS/SIGNIFICANCE: When appropriate filtering and mapping procedures are applied UHTS technology can be safely and efficiently used as a tool for targeted human DNA variations detection. Unless particular and platform-dependent characteristics are needed for specific projects, the most relevant parameter to consider in mainstream human genome resequencing procedures is the cost per sequenced base-pair associated to each machine.

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.

Islands of euchromatin-like sequence and expressed polymorphic sequences within the short arm of human chromosome 21

The goals of the human genome project did not include sequencing of the heterochromatic regions. We describe here an initial sequence of 1.1 Mb of the short arm of human chromosome 21 (HSA21p), estimated to be 10% of 21p. This region contains extensive euchromatic-like sequence and includes on average one transcript every 100 kb. These transcripts show multiple inter- and intrachromosomal copies, and extensive copy number and sequence variability. The sequencing of the "heterochromatic" regions of the human genome is likely to reveal many additional functional elements and provide important evolutionary information.

EGASP: the human ENCODE Genome Annotation Assessment Project

Background: We present the results of EGASP, a community experiment to assess the state-ofthe-art in genome annotation within the ENCODE regions, which span 1% of the human genomesequence. The experiment had two major goals: the assessment of the accuracy of computationalmethods to predict protein coding genes; and the overall assessment of the completeness of thecurrent human genome annotations as represented in the ENCODE regions. For thecomputational prediction assessment, eighteen groups contributed gene predictions. Weevaluated these submissions against each other based on a ‘reference set’ of annotationsgenerated as part of the GENCODE project. These annotations were not available to theprediction groups prior to the submission deadline, so that their predictions were blind and anexternal advisory committee could perform a fair assessment.Results: The best methods had at least one gene transcript correctly predicted for close to 70%of the annotated genes. Nevertheless, the multiple transcript accuracy, taking into accountalternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotidelevel, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programsrelying on mRNA and protein sequences were the most accurate in reproducing the manuallycurated annotations. Experimental validation shows that only a very small percentage (3.2%) of the selected 221 computationally predicted exons outside of the existing annotation could beverified.Conclusions: This is the first such experiment in human DNA, and we have followed thestandards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe theresults presented here contribute to the value of ongoing large-scale annotation projects and shouldguide further experimental methods when being scaled up to the entire human genome sequence.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

MGMT promoter methylation is prognostic but not predictive for outcome to adjuvant PCV chemotherapy in anaplastic oligodendroglial tumors: a report from EORTC Brain Tumor Group Study 26951.

PURPOSE: O6-methylguanine-methyltransferase (MGMT) promoter methylation has been shown to predict survival of patients with glioblastomas if temozolomide is added to radiotherapy (RT). It is unknown if MGMT promoter methylation is also predictive to outcome to RT followed by adjuvant procarbazine, lomustine, and vincristine (PCV) chemotherapy in patients with anaplastic oligodendroglial tumors (AOT). PATIENTS AND METHODS: In the European Organisation for the Research and Treatment of Cancer study 26951, 368 patients with AOT were randomly assigned to either RT alone or to RT followed by adjuvant PCV. From 165 patients of this study, formalin-fixed, paraffin-embedded tumor tissue was available for MGMT promoter methylation analysis. This was investigated with methylation specific multiplex ligation-dependent probe amplification. RESULTS: In 152 cases, an MGMT result was obtained, in 121 (80%) cases MGMT promoter methylation was observed. Methylation strongly correlated with combined loss of chromosome 1p and 19q loss (P = .00043). In multivariate analysis, MGMT promoter methylation, 1p/19q codeletion, tumor necrosis, and extent of resection were independent prognostic factors. The prognostic significance of MGMT promoter methylation was equally strong in the RT arm and the RT/PCV arm for both progression-free survival and overall survival. In tumors diagnosed at central pathology review as glioblastoma, no prognostic effect of MGMT promoter methylation was observed. CONCLUSION: In this study, on patients with AOT MGMT promoter methylation was of prognostic significance and did not have predictive significance for outcome to adjuvant PCV chemotherapy. The biologic effect of MGMT promoter methylation or pathogenetic features associated with MGMT promoter methylation may be different for AOT compared with glioblastoma.