A user's guide to the encyclopedia of DNA elements (ENCODE).

The mission of the Encyclopedia of DNA Elements (ENCODE) Project is to enable the scientific and medical communities to interpret the human genome sequence and apply it to understand human biology and improve health. The ENCODE Consortium is integrating multiple technologies and approaches in a collective effort to discover and define the functional elements encoded in the human genome, including genes, transcripts, and transcriptional regulatory regions, together with their attendant chromatin states and DNA methylation patterns. In the process, standards to ensure high-quality data have been implemented, and novel algorithms have been developed to facilitate analysis. Data and derived results are made available through a freely accessible database. Here we provide an overview of the project and the resources it is generating and illustrate the application of ENCODE data to interpret the human genome.

EGASP: the human ENCODE Genome Annotation Assessment Project.

BACKGROUND: We present the results of EGASP, a community experiment to assess the state-of-the-art in genome annotation within the ENCODE regions, which span 1% of the human genome sequence. The experiment had two major goals: the assessment of the accuracy of computational methods to predict protein coding genes; and the overall assessment of the completeness of the current human genome annotations as represented in the ENCODE regions. For the computational prediction assessment, eighteen groups contributed gene predictions. We evaluated these submissions against each other based on a 'reference set' of annotations generated as part of the GENCODE project. These annotations were not available to the prediction groups prior to the submission deadline, so that their predictions were blind and an external 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 account alternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotide level, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programs relying on mRNA and protein sequences were the most accurate in reproducing the manually curated 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 be verified. CONCLUSION: This is the first such experiment in human DNA, and we have followed the standards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe the results presented here contribute to the value of ongoing large-scale annotation projects and should guide further experimental methods when being scaled up to the entire human genome sequence.

EGASP: the human ENCODE Genome Annotation Assessment Project

Background: Non-long terminal repeat (non-LTR) retrotransposons have contributed to shaping the structure and function of genomes. In silico and experimental approaches have been used to identify the non-LTR elements of the urochordate Ciona intestinalis. Knowledge of the types and abundance of non-LTR elements in urochordates is a key step in understanding their contribution to the structure and function of vertebrate genomes. Results: Consensus elements phylogenetically related to the I, LINE1, LINE2, LOA and R2 elements of the 14 eukaryotic non-LTR clades are described from C. intestinalis. The ascidian elements showed conservation of both the reverse transcriptase coding sequence and the overall structural organization seen in each clade. The apurinic/apyrimidinic endonuclease and nucleic-acid-binding domains encoded upstream of the reverse transcriptase, and the RNase H and the restriction enzyme-like endonuclease motifs encoded downstream of the reverse transcriptase were identified in the corresponding Ciona families. Conclusions: The genome of C. intestinalis harbors representatives of at least five clades of non-LTR retrotransposons. The copy number per haploid genome of each element is low, less than 100, far below the values reported for vertebrate counterparts but within the range for protostomes. Genomic and sequence analysis shows that the ascidian non-LTR elements are unmethylated and flanked by genomic segments with a gene density lower than average for the genome. The analysis provides valuable data for understanding the evolution of early chordate genomes and enlarges the view on the distribution of the non-LTR retrotransposons in eukaryotes.

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.

Hypermethylation-mediated regulation of CD44 gene expression in human neuroblastoma

The CD44 adhesion receptor is silenced in highly malignant neuroblastomas (NBs) with MYCN amplification. Because its functional expression is associated with decreased tumorigenic properties, CD44 behaves as a tumor suppressor gene in NB and other cancers. Given that the precise mechanisms responsible for CD44 silencing are not elucidated, we investigated whether CD44 expression could be regulated by DNA hypermethylation. The methylation status of CD44 gene promoter and exon 1 regions was analyzed in 12 NB cell lines and 21 clinical samples after bisulfite genomic modification, followed by PCR and single-strand conformation polymorphism analysis and genomic sequencing. The results showed that almost all CD44-negative cell lines displayed hypermethylation in both regions, whereas all CD44-expressing cell lines were unmethylated. These observations correlated with the ability to restore CD44 mRNA and protein expression by treatment of CD44-negative cells with the 5-aza-2'-deoxycytidine demethylating agent. In contrast, no CD44 gene hypermethylation could be detected in 21 NB clinical samples of different stages, irrespective of CD44 expression. Although our results suggest that aberrant methylation of promoter and exon 1 regions is involved in CD44 silencing in NB cell lines, they also indicate that methylation of unidentified regulatory sequences or methylation-independent mechanisms also control the expression of CD44 in primary NB tumors and cell lines. We therefore conclude that CD44 silencing is controlled by complex and tumor cell-specific processes, including gene hypermethylation. Further investigation of other mechanisms and genes involved in CD44 regulation will be needed before demethylation-mediated reactivation of the CD44 gene can be considered as therapeutic strategy for neuroblastoma and perhaps other related cancers.

Epigenetic alterations in human brain tumors in a Brazilian population

Aberrant methylation of CpG islands located in promoter regions represents one of the major mechanisms for silencing cancer-related genes in tumor cells. We determined the frequency of aberrant CpG island methylation for several tumor-associated genes: DAPK, MGMT, p14^ARF, p16^INK4a, TP73, RB1 and TIMP-3 in 55 brain tumors, consisting of 26 neuroepithelial tumors, 6 peripheral nerve tumors, 13 meningeal tumors and 10 metastatic brain tumors. Aberrant methylation of at least one of the seven genes studied was detected in 83.6% of the cases. The frequencies of aberrant methylation were: 40% for p14^ARF, 38.2% for MGMT, 30.9% for, p16^INK4a, 14.6% for TP73 and for TIMP-3, 12.7% for DAPK and 1.8% for RB1. These data suggest that the hypermethylation observed in the genes p14^ARF, MGMT and p16^INK4a is a very important event in the formation or progression of brain tumors, since the inactivation of these genes directly interferes with the cell cycle or DNA repair. The altered methylation rate of the other genes has already been reported to be related to tumorigenesis, but the low methylation rate of RB1 found in tumors in our sample is different from that so far reported in the literature, suggesting that perhaps hypermethylation of the promoter is not the main event in the inactivation of this gene. Our results suggest that hypermethylation of the promoter region is a very common event in nervous system tumors.

Comprehensive genome methylation and whole genome expression analysis in penile carcinoma: Uncovering new molecular markers

Background: Penile carcinoma (PeCa) is frequently associated with high morbidity rates. Unlikely of the vast majority of tumors, there is no molecular markers described that are able to assist in diagnosis and prognosis or with potential to be therapeutic targets in PeCa. Patients and methods: DNA methylation status (244K Human DNA Methylation Microarray platform, Agilent Technologies) and large-scale expression analysis (4x44K Whole Human Genome Microarray, Agilent Technologies) were performed in 35 and 37 PeCa, respectively. Quantitative bisulfite pyrosequencing (qBP) and RT-qPCR were used to validate the findings in 93 samples. HPV status was assessed using the Linear Array HPV Genotyping kit (Roche Molecular Diagnostics, CA, USA). Results: Methylome analysis revealed 171 hypermethylated and 449 hypomethylated CpGs sites and the transcriptome profiling showed 2986 down- and 2817 over-expressed genes. HPV positivity was found in 32.7% of the cases, mainly the HPV16. The integrative analysis in 32 PeCa revealed a panel of 96 genes with inverse correlation between methylation and gene expression levels. The CpG hypermetlylation and gene downexpression, was confirmed for TWIST1, RSOP2, SOX3, SOX17, CD133, OTX2, HOXA3 and MEIS. In addition, BIRC5, DNMT1 and DNMT3B presented low levels of methylation and overexpression. The comparison of the results with clinical findings revealed that LIN28A, NKX2.2, NKX2.3, LHX5, BDNF, FOXA1 and CDX2 were associated with poor prognosis features. Conclusion: Putative prognostic markers were detected revealing that DNA methylation modulates the expression of several genes in PeCa. These data may prove instrumental for biomarker discovery in clinics and molecular epidemiology of PeCa.

DNA repair mechanisms protect our genome from carcinogenesis

Human cells are constantly exposed to DNA damage. Without repair, damage can result in genetic instability and eventually cancer. The strong association between the lack of DNA damage repair, mutations and cancer is dramatically demonstrated by a number of cancer-prone human syndromes, such as xeroderma pigmentosum (XP), ataxia-telangiectasia (AT) and Fanconi anemia (FA). This review focuses on the historical discoveries related with these three diseases and describes their impact on the understanding of DNA repair mechanisms and the causes of human cancer. As deficiencies in DNA repair are also often related with progeria symptoms, unrepaired damage and aging are somehow related. Several other pathologies associated with DNA repair defects, genetic instability and increased cancer risk are also discussed. In fact, studies with cells from these many syndromes have helped in understanding important levels of protection against cancer and aging, although little help has actually been conferred to the patients in terms of therapy. Finally, the recent advances in combined basic and translational research on DNA repair and chemotherapy are presented.

DNA methylation changes associated with risk factors in tumors of the upper aerodigestive tract

Cancers of the upper aerodigestive tract (UADT) are common forms of malignancy associated with tobacco and alcohol exposures, although human papillomavirus and nutritional deficiency are also important risk factors. While somatically acquired DNA methylation changes have been associated with UADT cancers, what triggers these events and precise epigenetic targets are poorly understood. In this study, we applied quantitative profiling of DNA methylation states in a panel of cancer-associated genes to a case-control study of UADT cancers. Our analyses revealed a high frequency of aberrant hypermethylation of several genes, including MYOD1, CHRNA3 and MTHFR in UADT tumors, whereas CDKN2A was moderately hypermethylated. Among differentially methylated genes, we identified a new gene (the nicotinic acetycholine receptor gene) as target of aberrant hypermethylation in UADT cancers, suggesting that epigenetic deregulation of nicotinic acetycholine receptors in non-neuronal tissues may promote the development of UADT cancers. Importantly, we found that sex and age is strongly associated with the methylation states, whereas tobacco smoking and alcohol intake may also influence the methylation levels in specific genes. This study identifies aberrant DNA methylation patterns in UADT cancers and suggests a potential mechanism by which environmental factors may deregulate key cellular genes involved in tumor suppression and contribute to UADT cancers.

Epigenetic regulation of somatic angiotensin-converting enzyme by DNA methylation and histone acetylation

Somatic angiotensin-converting enzyme (sACE) is crucial in cardiovascular homeostasis and displays a tissue-specific profile. Epigenetic patterns modulate genes expression and their alterations were implied in pathologies including hypertension. However, the influence of DNA methylation and chromatin condensation state on the expression of sACE is unknown. We examined whether such epigenetic mechanisms could participate in the control of sACE expression in vitro and in vivo. We identified two CpG islands in the human ace-1 gene 3 kb proximal promoter region. Their methylation abolished the luciferase activity of ace-1 promoter/reporter constructs transfected into human liver (HepG2), colon (HT29), microvascular endothelial (HMEC-1) and lung (SUT) cell lines (p < 0.001). Bisulphite sequencing revealed a cell-type specific basal methylation pattern of the ace-1 gene -1,466/+25 region. As assessed by RT-qPCR, inhibition of DNA methylation by 5-aza-2'-deoxycytidine and/or of histone deacetylation by trichostatin A highly stimulated sACE mRNA expression cell-type specifically (p < 0.001 vs. vehicle treated cells). In the rat, in vivo 5-aza-cytidine injections demethylated the ace-1 promoter and increased sACE mRNA expression in the lungs and liver (p = 0.05), but not in the kidney. In conclusion, the expression level of somatic ACE is modulated by CpG-methylation and histone deacetylases inhibition. The basal methylation pattern of the promoter of the ace-1 gene is cell-type specific and correlates to sACE transcription. DNMT inhibition is associated with altered methylation of the ace-1 promoter and a cell-type and tissue-specific increase of sACE mRNA levels. This study indicates a strong influence of epigenetic mechanisms on sACE expression.

Evolutionary forces shaping genomic islands of population differentiation in humans

Background Levels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure. Methods Individual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation. Results Under the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process. Conclusions Even though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.

A HIDDEN MARKOV MODEL FOR JOINT ESTIMATION OF GENOTYPE AND COPY NUMBER IN HIGH-THROUGHPUT SNP CHIPS

Amplifications and deletions of chromosomal DNA, as well as copy-neutral loss of heterozygosity have been associated with diseases processes. High-throughput single nucleotide polymorphism (SNP) arrays are useful for making genome-wide estimates of copy number and genotype calls. Because neighboring SNPs in high throughput SNP arrays are likely to have dependent copy number and genotype due to the underlying haplotype structure and linkage disequilibrium, hidden Markov models (HMM) may be useful for improving genotype calls and copy number estimates that do not incorporate information from nearby SNPs. We improve previous approaches that utilize a HMM framework for inference in high throughput SNP arrays by integrating copy number, genotype calls, and the corresponding confidence scores when available. Using simulated data, we demonstrate how confidence scores control smoothing in a probabilistic framework. Software for fitting HMMs to SNP array data is available in the R package ICE.

Whole genome methylation array analysis reveals new aspects in Balkan endemic nephropathy etiology

Background Balkan endemic nephropathy (BEN) represents a chronic progressive interstitial nephritis in striking correlation with uroepithelial tumours of the upper urinary tract. The disease has endemic distribution in the Danube river regions in several Balkan countries. DNA methylation is a primary epigenetic modification that is involved in major processes such as cancer, genomic imprinting, gene silencing, etc. The significance of CpG island methylation status in normal development, cell differentiation and gene expression is widely recognized, although still stays poorly understood. Methods We performed whole genome DNA methylation array analysis on DNA pool samples from peripheral blood from 159 affected individuals and 170 healthy individuals. This technique allowed us to determine the methylation status of 27 627 CpG islands throughout the whole genome in healthy controls and BEN patients. Thus we obtained the methylation profile of BEN patients from Bulgarian and Serbian endemic regions. Results Using specifically developed software we compared the methylation profiles of BEN patients and corresponding controls and revealed the differently methylated regions. We then compared the DMRs between all patient-control pairs to determine common changes in the epigenetic profiles. SEC61G, IL17RA, HDAC11 proved to be differently methylated throughout all patient-control pairs. The CpG islands of all 3 genes were hypomethylated compared to controls. This suggests that dysregulation of these genes involved in immunological response could be a common mechanism in BEN pathogenesis in both endemic regions and in both genders. Conclusion Our data propose a new hypothesis that immunologic dysregulation has a place in BEN etiopathogenesis. Keywords: Epigenetics; Whole genome array analysis; Balkan endemic nephropathy

Whole genome amplification using a degenerate oligonucleotide primer allows hundreds of genotypes to be performed on less than one nanogram of genomic DNA

Genetic analysis of limiting quantities of genomic DNA play an important role in DNA forensics, paleoarcheology, genetic disease diagnosis, genetic linkage analysis, and genetic diversity studies. We have tested the ability of degenerate oligonucleotide primed polymerase chain reaction (DOP-PCR) to amplify picogram quantities of human genomic DNA for the purpose of increasing the amount of template for genotyping with microsatellite repeat markers. DNA was uniformly amplified at a large number of typable loci throughout the human genome with starting template DNAs from as little as 15 pg to as much as 400 ng. A much greater-fold enrichment was seen for the smaller genomic DOP-PCRs. All markers tested were amplified from starting genomic DNAs in the range of 0.6–40 ng with amplifications of 200- to 600-fold. The DOP-PCR-amplified genomic DNA was an excellent and reliable template for genotyping with microsatellites, which give distinct bands with no increase in stutter artifact on di-, tri-, and tetranucleotide repeats. There appears to be equal amplification of genomic DNA from 55 of 55 tested discrete microsatellites implying near complete coverage of the human genome. Thus, DOP-PCR appears to allow unbiased, hundreds-fold whole genome amplification of human genomic DNA for genotypic analysis.

CpG DNA can induce strong Th1 humoral and cell-mediated immune responses against hepatitis B surface antigen in young mice

Successful neonatal immunization of humans has proven difficult. We have evaluated CpG-containing oligonucleotides as an adjuvant for immunization of young mice (1–14 days old) against hepatitis B virus surface antigen. The protein-alum-CpG formulation, like the DNA vaccine, produced seroconversion of the majority of mice immunized at 3 or 7 days of age, compared with 0–10% with the protein-alum or protein-CpG formulations. All animals, from neonates to adults, immunized with the protein-alum vaccine exhibited strong T helper (Th)2-like responses [predominantly IgG1, weak or absent cytotoxic T lymphocytes (CTL)]. Th2-type responses also were induced in young mice with protein-CpG (in 1-, 3-, and 7-day-old mice) and protein-alum-CpG (in 1- and 3-day-old mice) but immunization carried out at older ages gave mixed Th1/Th2 (Th0) responses. DNA vaccines gave Th0-like responses when administered at 1 and 7 days of age and Th1-like (predominantly IgG2a and CTL) responses with 14-day-old or adult mice. Surprisingly, the protein-alum-CpG formulation was better than the DNA vaccine for percentage of seroconversion, speed of appearance, and peak titer of the antibody response, as well as prevalence and strength of CTL. These findings may have important implications for immunization of human infants.