Kristen H. Taylor

Ultradeep bisulfite sequencing analysis of DNA methylation patterns in multiple gene promoters by 454 sequencing

We developed a novel approach for conducting multisample, multigene, ultradeep bisulfite sequencing analysis of DNA methylation patterns in clinical samples. A massively parallel sequencing-by-synthesis method (454 sequencing) was used to directly sequence >100 bisulfite PCR products in a single sequencing run without subcloning. We showed the utility, robustness, and superiority of this approach by analyzing methylation in 25 gene-related CpG rich regions from >40 cases of primary cells, including normal peripheral blood lymphocytes, acute lymphoblastic leukemia (ALL), chronic lymphocytic leukemia (CLL), follicular lymphoma (FL), and mantle cell lymphoma (MCL). A total of 294,631 sequences was generated with an average read length of 131 bp. On average, >1,600 individual sequences were generated for each PCR amplicon far beyond the few clones (<20) typically analyzed by traditional bisulfite sequencing. Comprehensive analysis of CpG methylation patterns at a single DNA molecule level using clustering algorithms revealed differential methylation patterns between diseases. A significant increase in methylation was detected in ALL and FL samples compared with CLL and MCL. Furthermore, a progressive spreading of methylation was detected from the periphery toward the center of select CpG islands in the ALL and FL samples. The ultradeep sequencing also allowed simultaneous analysis of genetic and epigenetic data and revealed an association between a single nucleotide polymorphism and the methylation present in the LRP1B promoter. This new generation of methylome sequencing will provide digital profiles of aberrant DNA methylation for individual human cancers and offers a robust method for the epigenetic classification of tumor subtypes.

Large-scale CpG methylation analysis identifies novel candidate genes and reveals methylation hotspots in acute lymphoblastic leukemia

This study examined DNA methylation associated with acute lymphoblastic leukemia (ALL) and showed that selected molecular targets can be pharmacologically modulated to reverse gene silencing. A CpG island (CGI) microarray containing more than 3,400 unique clones that span all human chromosomes was used for large-scale discovery experiments and led to 262 unique CGI loci being statistically identified as methylated in ALL lymphoblasts. The methylation status of 10 clones encompassing 11 genes (DCC, DLC-1, DDX51, KCNK2, LRP1B, NKX6-1, NOPE, PCDHGA12, RPIB9, ABCB1, and SLC2A14) identified as differentially methylated between ALL patients and controls was independently verified. Consequently, the methylation status of DDX51 was found to differentiate patients with B- and T-ALL subtypes (P = 0.011, Fishers exact test). Next, the relationship between methylation and expression of these genes was examined in ALL cell lines (NALM-6 and Jurkat) before and after treatments with 5-aza-2-deoxycytidine and trichostatin A. More than a 10-fold increase in mRNA expression was observed for two previously identified tumor suppressor genes (DLC-1 and DCC) and also for RPIB9 and PCDHGA12. Although the mechanisms that lead to the CGI methylation of these genes are unknown, bisulfite sequencing of the promoter of RPIB9 suggests that expression is inhibited by methylation within SP1 and AP2 transcription factor binding motifs. Finally, specific chromosomal methylation hotspots were found to be associated with ALL. This study sets the stage for acquiring a better biological understanding of ALL and for the identification of epigenetic biomarkers useful for differential diagnosis, therapeutic monitoring, and the detection of leukemic relapse.

Differential DNA methylation patterns of small B-cell lymphoma subclasses with different clinical behavior.

Non-Hodgkins lymphoma (NHL) is a group of malignancies of the immune system with variable clinical behaviors and diverse molecular features. Despite the progress made in classification of NHLs based on classical methods, molecular classifications are a work in progress. Toward this goal, we used an array-based technique called differential methylation hybridization (DMH) to study small B-cell lymphoma (SBCL) subtypes. A total of 43 genomic DMH experiments were performed. From these results, several statistical methods were used to generate a set of differentially methylated genes for further validation. Methylation of LHX2, POU3F3, HOXC10, NRP2, PRKCE, RAMP, MLLT2, NKX6.1, LRP1B and ARF4 was validated in cell lines and patient samples and demonstrated subtype-related preferential methylation patterns. For LHX2 and LRP1B, bisulfite sequencing, real-time reverse transcriptase-polymerase chain reaction and induction of gene expression following treatment with the demethylating agent, 5′-aza-2′-deoxycytidine, were confirmed. This new epigenetic information is helping to define molecular portraits of distinct subtypes of SBCL that are not recognized by current classification systems and provides valuable potential insights into the biology of these tumors.

Genome-wide DNA methylation analysis reveals novel epigenetic changes in chronic lymphocytic leukemia

We conducted a genome-wide DNA methylation analysis in CD19+ B-cells from chronic lymphocytic leukemia (CLL) patients and normal control samples using reduced representation bisulfite sequencing (RRBS). The methylation status of 1.8–2.3 million CpGs in the CLL genome was determined; about 45% of these CpGs were located in more than 23,000 CpG islands (CGIs). While global CpG methylation was similar between CLL and normal B-cells, 1764 gene promoters were identified as being differentially methylated in at least one CLL sample when compared with normal B-cell samples. Nineteen percent of the differentially methylated genes were involved in transcriptional regulation. Aberrant hypermethylation was found in all HOX gene clusters and a significant number of WNT signaling pathway genes. Hypomethylation occurred more frequently in the gene body including introns, exons, and 3′-UTRs in CLL. The NFATc1 P2 promoter and first intron was found to be hypomethylated and correlated with upregulation of both NFATc1 RNA and protein expression levels in CLL suggesting that an epigenetic mechanism is involved in the constitutive activation of NFAT activity in CLL cells. This comprehensive DNA methylation analysis will further our understanding of the epigenetic contribution to cellular dysfunction in CLL.

DNA hypermethylation accompanied by transcriptional repression in follicular lymphoma

High‐throughput microarray technologies were used to study DNA methylation accompanied by transcriptional changes in follicular lymphoma (FL). Using Methylated CpG Island Amplification with Microarrays to study CpG Island DNA methylation in FL, we discovered widespread hypermethylation of homeobox genes and previously identified targets of polycomb repressive complex 2 (PRC2) in cell lines and primary tumors, but not in benign follicular hyperplasia (BFH). DNA methylation for HOXA11, HOXD10, HOXB7, HOXC12, PAX6, LHX9, SFMBT2, EN2, and PAX7 was independently validated in the RL cell line and HOXA11, HOXD10, PAX6, and EN2 in primary tumors. Combined Bisulfite Restriction Analysis (COBRA) also established DNA methylation for the previously identified PRC2 targets DCC, DES, GAD2, AQP5, GPR61, GRIA4, GJD2, and AMPH in FL but not in BFH. Gene expression analyses revealed 411 genes that were hypermethylated and transcriptionally repressed in RL, 74% of which were reactivated by the demethylating agent 5‐aza‐2′‐deoxycytidine (5‐azaD) plus or minus the histone deacetylase inhibitor trichostatin A (TSA). Forty genes were also downregulated in primary FL. Our results suggest that extensive hypermethylation in promoters of polycomb target genes is a characteristic of FL and that loss of expression of certain SUZ12 target genes could be functionally relevant for lymphomagenesis.

Large-scale analysis of DNA methylation in chronic lymphocytic leukemia

AIMS B-cell chronic lymphocytic leukemia (CLL) is a heterogeneous malignancy that clinically ranges from indolent to rapidly progressive. CLL, like other cancers, can be affected by epigenetic alterations. MATERIALS & METHODS A microarray discovery-based study was initiated to determine DNA methylation in CLL cases with a range of CD38 expression (1–92%). RESULTS Many loci were either methylated or unmethylated across all CD38 levels, but differential methylation was also observed for some genes. Genomic sequencing of DLEU7 confirmed extensive cytosine methylation preferentially in patient samples with low CD38 expression, whereas NRP2, SFRP2 and ADAM12 were more commonly methylated in those with high CD38 expression. CONCLUSION This study demonstrates that CLL is affected by CpG island methylation in some genes that segregate with CD38 expression levels, while most others show similar methylation patterns across all levels. The CpG island methylation in certain functional gene groups and pathway-associated genes that are known to be deregulated in CLL provides additional insights into the CLL methylome and epigenetic contribution to cellular dysfunction. It will now be useful to investigate the effectiveness of epigenetic therapeutic reversal of these alterations to develop effective treatments for the disease.

Genome-wide DNA methylation maps in follicular lymphoma cells determined by methylation-enriched bisulfite sequencing

Background Follicular lymphoma (FL) is a form of non-Hodgkins lymphoma (NHL) that arises from germinal center (GC) B-cells. Despite the significant advances in immunotherapy, FL is still not curable. Beyond transcriptional profiling and genomics datasets, there currently is no epigenome-scale dataset or integrative biology approach that can adequately model this disease and therefore identify novel mechanisms and targets for successful prevention and treatment of FL. Methodology/Principal Findings We performed methylation-enriched genome-wide bisulfite sequencing of FL cells and normal CD19+ B-cells using 454 sequencing technology. The methylated DNA fragments were enriched with methyl-binding proteins, treated with bisulfite, and sequenced using the Roche-454 GS FLX sequencer. The total number of bases covered in the human genome was 18.2 and 49.3 million including 726,003 and 1.3 million CpGs in FL and CD19+ B-cells, respectively. 11,971 and 7,882 methylated regions of interest (MRIs) were identified respectively. The genome-wide distribution of these MRIs displayed significant differences between FL and normal B-cells. A reverse trend in the distribution of MRIs between the promoter and the gene body was observed in FL and CD19+ B-cells. The MRIs identified in FL cells also correlated well with transcriptomic data and ChIP-on-Chip analyses of genome-wide histone modifications such as tri-methyl-H3K27, and tri-methyl-H3K4, indicating a concerted epigenetic alteration in FL cells. Conclusions/Significance This study is the first to provide a large scale and comprehensive analysis of the DNA methylation sequence composition and distribution in the FL epigenome. These integrated approaches have led to the discovery of novel and frequent targets of aberrant epigenetic alterations. The genome-wide bisulfite sequencing approach developed here can be a useful tool for profiling DNA methylation in clinical samples.

The Properties of Genome Conformation and Spatial Gene Interaction and Regulation Networks of Normal and Malignant Human Cell Types

The spatial conformation of a genome plays an important role in the long-range regulation of genome-wide gene expression and methylation, but has not been extensively studied due to lack of genome conformation data. The recently developed chromosome conformation capturing techniques such as the Hi-C method empowered by next generation sequencing can generate unbiased, large-scale, high-resolution chromosomal interaction (contact) data, providing an unprecedented opportunity to investigate the spatial structure of a genome and its applications in gene regulation, genomics, epigenetics, and cell biology. In this work, we conducted a comprehensive, large-scale computational analysis of this new stream of genome conformation data generated for three different human leukemia cells or cell lines by the Hi-C technique. We developed and applied a set of bioinformatics methods to reliably generate spatial chromosomal contacts from high-throughput sequencing data and to effectively use them to study the properties of the genome structures in one-dimension (1D) and two-dimension (2D). Our analysis demonstrates that Hi-C data can be effectively applied to study tissue-specific genome conformation, chromosome-chromosome interaction, chromosomal translocations, and spatial gene-gene interaction and regulation in a three-dimensional genome of primary tumor cells. Particularly, for the first time, we constructed genome-scale spatial gene-gene interaction network, transcription factor binding site (TFBS) – TFBS interaction network, and TFBS-gene interaction network from chromosomal contact information. Remarkably, all these networks possess the properties of scale-free modular networks.

Aberrant Epigenetic Gene Regulation in Lymphoid Malignancies

In lymphoid malignancies, aberrant epigenetic mechanisms such as DNA methylation and histone modifications influence chromatin architecture and can result in altered gene expression. These alterations commonly affect genes that play important roles in the cell cycle, apoptosis, and DNA repair in non-Hodgkin lymphoma (NHL). The ability to identify epigenetic modifications to these important genes has increased exponentially due to advances in technology. As a result, there are well-defined, gene-specific epigenetic aberrations associated with NHL comprising follicular lymphoma (FL), mantle cell lymphoma (MCL), chronic lymphocytic leukemia (CLL), and diffuse large B-cell lymphoma (DLBCL). The identification of these genes is important because they may be used as biomarkers for prognosis, diagnosis and in developing improved treatment strategies. Also important, in the control of gene expression, is the packaging of DNA within the nucleus of a cell. This packaging can be distorted by epigenetic alterations and may alter the accessibility of certain regions of the genome in cancer cells. This review discusses the impact of known epigenetic aberration on the regulation of gene expression in NHL and provides insight into the spatial conformation of the genome (DNA packaging) in acute lymphoblastic leukemia.

Integrated methylome and transcriptome analysis reveals novel regulatory elements in pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common cancer diagnosed in children under the age of 15. In addition to genetic aberrations, epigenetic modifications such as DNA methylation are altered in cancer and impact gene expression. To identify epigenetic alterations in ALL, genome-wide methylation profiles were generated using the methylated CpG island recovery assay followed by next-generation sequencing. More than 25,000 differentially methylated regions (DMR) were observed in ALL patients with ∼90% present within intronic or intergenic regions. To determine the regulatory potential of the DMR, whole-transcriptome analysis was performed and integrated with methylation data. Aberrant promoter methylation was associated with the altered expression of genes involved in transcriptional regulation, apoptosis, and proliferation. Novel enhancer-like sequences were identified within intronic and intergenic DMR. Aberrant methylation in these regions was associated with the altered expression of neighboring genes involved in cell cycle processes, lymphocyte activation and apoptosis. These genes include potential epi-driver genes, such as SYNE1, PTPRS, PAWR, HDAC9, RGCC, MCOLN2, LYN, TRAF3, FLT1, and MELK, which may provide a selective advantage to leukemic cells. In addition, the differential expression of epigenetic modifier genes, pseudogenes, and non-coding RNAs was also observed accentuating the role of erroneous epigenetic gene regulation in ALL.