J. Wade Davis

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.

Estrogen and Androgen Receptor Activities of Hydraulic Fracturing Chemicals and Surface and Ground Water in a Drilling-Dense Region

The rapid rise in natural gas extraction using hydraulic fracturing increases the potential for contamination of surface and ground water from chemicals used throughout the process. Hundreds of products containing more than 750 chemicals and components are potentially used throughout the extraction process, including more than 100 known or suspected endocrine-disrupting chemicals. We hypothesized that a selected subset of chemicals used in natural gas drilling operations and also surface and ground water samples collected in a drilling-dense region of Garfield County, Colorado, would exhibit estrogen and androgen receptor activities. Water samples were collected, solid-phase extracted, and measured for estrogen and androgen receptor activities using reporter gene assays in human cell lines. Of the 39 unique water samples, 89%, 41%, 12%, and 46% exhibited estrogenic, antiestrogenic, androgenic, and antiandrogenic activities, respectively. Testing of a subset of natural gas drilling chemicals revealed novel antiestrogenic, novel antiandrogenic, and limited estrogenic activities. The Colorado River, the drainage basin for this region, exhibited moderate levels of estrogenic, antiestrogenic, and antiandrogenic activities, suggesting that higher localized activity at sites with known natural gas-related spills surrounding the river might be contributing to the multiple receptor activities observed in this water source. The majority of water samples collected from sites in a drilling-dense region of Colorado exhibited more estrogenic, antiestrogenic, or antiandrogenic activities than reference sites with limited nearby drilling operations. Our data suggest that natural gas drilling operations may result in elevated endocrine-disrupting chemical activity in surface and ground water.

Effects of Vendor and Genetic Background on the Composition of the Fecal Microbiota of Inbred Mice

The commensal gut microbiota has been implicated as a determinant in several human diseases and conditions. There is mounting evidence that the gut microbiota of laboratory mice (Mus musculus) similarly modulates the phenotype of mouse models used to study human disease and development. While differing model phenotypes have been reported using mice purchased from different vendors, the composition and uniformity of the fecal microbiota in mice of various genetic backgrounds from different vendors is unclear. Using culture-independent methods and robust statistical analysis, we demonstrate significant differences in the richness and diversity of fecal microbial populations in mice purchased from two large commercial vendors. Moreover, the abundance of many operational taxonomic units, often identified to the species level, as well as several higher taxa, differed in vendor- and strain-dependent manners. Such differences were evident in the fecal microbiota of weanling mice and persisted throughout the study, to twenty-four weeks of age. These data provide the first in-depth analysis of the developmental trajectory of the fecal microbiota in mice from different vendors, and a starting point from which researchers may be able to refine animal models affected by differences in the gut microbiota and thus possibly reduce the number of animals required to perform studies with sufficient statistical power.

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.

Aberrant gene expression profile in a mouse model of endometriosis mirrors that observed in women

OBJECTIVE To define the altered gene expression profile of endometriotic lesions in a mouse model of surgically induced endometriosis. DESIGN Autologous experimental mouse model. SETTING Medical school department. ANIMAL(S) Adult C57Bl6 mice. INTERVENTION(S) Endometriosis was surgically induced by autotransplantation of uterine tissue to the intestinal mesentery. Endometriotic lesions and eutopic uteri were recovered at 3 or 29 days after induction. MAIN OUTCOME MEASURE(S) Altered gene expression was measured in the endometriotic lesion relative to the eutopic uterus by genome-wide complementary DNA microarray analysis and was confirmed by real-time reverse transcriptase-polymerase chain reaction for six genes. Relevant categories of altered genes were identified using gene ontology analysis to determine groups of genes enriched for altered expression. RESULT(S) The expression of 479 and 114 genes was altered in the endometriotic lesion compared with the eutopic uterus at 3 or 29 days after induction, respectively. Gene ontology enrichment analysis revealed that genes associated with the extracellular matrix, cell adhesions, immune function, cell growth, and angiogenesis were altered in the endometriotic lesion compared with the eutopic uterus. CONCLUSION(S) According to gene expression analysis, the mouse model of surgically induced endometriosis is a good model for studying the pathophysiology and treatment of endometriosis.

The RNA binding protein HuR differentially regulates unique subsets of mRNAs in estrogen receptor negative and estrogen receptor positive breast cancer

BackgroundThe discordance between steady-state levels of mRNAs and protein has been attributed to posttranscriptional control mechanisms affecting mRNA stability and translation. Traditional methods of genome wide microarray analysis, profiling steady-state levels of mRNA, may miss important mRNA targets owing to significant posttranscriptional gene regulation by RNA binding proteins (RBPs).MethodsThe ribonomic approach, utilizing RNA immunoprecipitation hybridized to microarray (RIP-Chip), provides global identification of putative endogenous mRNA targets of different RBPs. HuR is an RBP that binds to the AU-rich elements (ARE) of labile mRNAs, such as proto-oncogenes, facilitating their translation into protein. HuR has been shown to play a role in cancer progression and elevated levels of cytoplasmic HuR directly correlate with increased invasiveness and poor prognosis for many cancers, including those of the breast. HuR has been described to control genes in several of the acquired capabilities of cancer and has been hypothesized to be a tumor-maintenance gene, allowing for cancers to proliferate once they are established.ResultsWe used HuR RIP-Chip as a comprehensive and systematic method to survey breast cancer target genes in both MCF-7 (estrogen receptor positive, ER+) and MDA-MB-231 (estrogen receptor negative, ER-) breast cancer cell lines. We identified unique subsets of HuR-associated mRNAs found individually or in both cell types. Two novel HuR targets, CD9 and CALM2 mRNAs, were identified and validated by quantitative RT-PCR and biotin pull-down analysis.ConclusionThis is the first report of a side-by-side genome-wide comparison of HuR-associated targets in wild type ER+ and ER- breast cancer. We found distinct, differentially expressed subsets of cancer related genes in ER+ and ER- breast cancer cell lines, and noted that the differential regulation of two cancer-related genes by HuR was contingent upon the cellular environment.

Gastrotomy creation and closure for NOTES using a gastropexy technique (with video)

BACKGROUND Safe and efficient gastrotomy creation and closure is pivotal for natural orifice transluminal endoscopic surgery (NOTES). OBJECTIVE To test a method of transgastric access and closure with commercially available devices. DESIGN An animal survival study. SETTING University hospital. PATIENTS Fifteen pigs. INTERVENTIONS By using a surgical suture passer, under endoscopic guidance, 3 percutaneous stay sutures were placed, in a triangular fashion, through the gastric wall. A gastrotomy was created with a dilation balloon, which was introduced over a guidewire through the gastric wall in the center of the 3 sutures. After performing a NOTES procedure, the gastrotomy was closed by tying the sutures. Necropsies were performed after 2 to 4 weeks. MAIN OUTCOME MEASUREMENTS Success and time of gastrotomy creation and closure, and intraoperative and postoperative complications. RESULTS Gastrotomies were successfully created and closed in all the animals. The median time to create a gastrotomy was 19 minutes (range 11-85 minutes), and the median closure time was 1 minute (range 1-45 minutes). One pig died on postoperative day 1 because of peritonitis caused by a leaking gastrotomy site that extended beyond the stay sutures. There were no other gastrotomy-related complications. All gastrotomies were well healed at the necropsy. LIMITATION No control group. CONCLUSIONS We evaluated a simple method by using the principles of the PEG technique combined with a gastropexy, which is familiar to the majority of endoscopists. Strict attention to the gastrotomy site is needed, because one leak was from the gastrotomy site that extended beyond the stay sutures.

Vascular transcriptional alterations produced by juvenile obesity in Ossabaw swine

We adopted a transcriptome-wide microarray analysis approach to determine the extent to which vascular gene expression is altered as a result of juvenile obesity and identify obesity-responsive mRNAs. We examined transcriptional profiles in the left anterior descending coronary artery (LAD), perivascular fat adjacent to the LAD, and descending thoracic aorta between obese (n = 5) and lean (n = 6) juvenile Ossabaw pigs (age = 22 wk). Obesity was experimentally induced by feeding the animals a high-fat/high-fructose corn syrup/high-cholesterol diet for 16 wk. We found that expression of 189 vascular cell genes in the LAD and expression of 165 genes in the thoracic aorta were altered with juvenile obesity (false discovery rate ≤ 10%) with an overlap of only 28 genes between both arteries. Notably, a number of genes found to be markedly upregulated in the LAD of obese pigs are implicated in atherosclerosis, including ACP5, LYZ, CXCL14, APOE, PLA2G7, LGALS3, SPP1, ITGB2, CYBB, and P2RY12. Furthermore, pathway analysis revealed the induction of proinflammatory and pro-oxidant pathways with obesity primarily in the LAD. Gene expression in the LAD perivascular fat was minimally altered with juvenile obesity. Together, we provide new evidence that obesity produces artery-specific changes in pretranslational regulation with a clear upregulation of proatherogenic genes in the LAD. Our data may offer potential viable drug targets and mechanistic insights regarding the molecular precursors involved in the origins of overnutrition and obesity-associated vascular disease. In particular, our results suggest that the oxidized LDL/LOX-1/NF-κB signaling axis may be involved in the early initiation of a juvenile obesity-induced proatherogenic coronary artery phenotype.

Transcriptome-wide RNA sequencing analysis of rat skeletal muscle feed arteries. I. Impact of obesity

We employed next-generation RNA sequencing (RNA-Seq) technology to determine the influence of obesity on global gene expression in skeletal muscle feed arteries. Transcriptional profiles of the gastrocnemius and soleus muscle feed arteries (GFA and SFA, respectively) and aortic endothelial cell-enriched samples from obese Otsuka Long-Evans Tokushima Fatty (OLETF) and lean Long-Evans Tokushima Otsuka (LETO) rats were examined. Obesity produced 282 upregulated and 133 downregulated genes in SFA and 163 upregulated and 77 downregulated genes in GFA [false discovery rate (FDR) < 10%] with an overlap of 93 genes between the arteries. In LETO rats, there were 89 upregulated and 114 downregulated genes in the GFA compared with the SFA. There were 244 upregulated and 275 downregulated genes in OLETF rats (FDR < 10%) in the GFA compared with the SFA, with an overlap of 76 differentially expressed genes common to both LETO and OLETF rats in both the GFA and SFA. A total of 396 transcripts were found to be differentially expressed between LETO and OLETF in aortic endothelial cell-enriched samples. Overall, we found 1) the existence of heterogeneity in the transcriptional profile of the SFA and GFA within healthy LETO rats, 2) that this between-vessel heterogeneity was markedly exacerbated in the hyperphagic, obese OLETF rat, and 3) a greater number of genes whose expression was altered by obesity in the SFA compared with the GFA. Also, results indicate that in OLETF rats the GFA takes on a relatively more proatherogenic phenotype compared with the SFA.