Search

DNA methylation is a labile modification associated with gene expression control and environmental adaptations. High throughput, scalable and quantitative assessments of specific DNA methylation modifications in complex genomic regions for use in large population studies are needed. The performance of Droplet Digital™ PCR (ddPCR™) was investigated for DNA methylation detection against next-generation bisulfite sequencing (NGS) to demonstrate the ability of ddPCR to detect and validate DNA methylation levels and complex patterns among neighboring CpGs in regions associated with prenatal tobacco exposure. While both techniques are reproducible, ddPCR demonstrates a unique advantage for high-throughput DNA methylation analysis in large-scale population studies and provides the specificity to accurately measure DNA methylation of target CpGs in complex regions.

Genome-wide association studies have identified a growing number of single nucleotide polymorphisms (SNPs) associated with childhood acute lymphoblastic leukemia (ALL), yet the functional roles of most SNPs are unclear. Multiple lines of evidence suggest that epigenetic mechanisms may mediate the impact of heritable genetic variation on phenotypes. Here, we investigated whether DNA methylation mediates the effect of genetic risk loci for childhood ALL. We performed an epigenome-wide association study (EWAS) including 808 childhood ALL cases and 919 controls from California-based studies using neonatal blood DNA. For differentially methylated CpG positions (DMPs), we next conducted association analysis with 23 known ALL risk SNPs followed by causal mediation analyses addressing the significant SNP-DMP pairs. DNA methylation at CpG cg01139861, in the promoter region of IKZF1, mediated the effects of the intronic IKZF1 risk SNP rs78396808, with the average causal mediation effect (ACME) explaining ~30% of the total effect (ACME P = 0.0031). In analyses stratified by self-reported race/ethnicity, the mediation effect was only significant in Latinos, explaining ~41% of the total effect of rs78396808 on ALL risk (ACME P = 0.0037). Conditional analyses confirmed the presence of at least three independent genetic risk loci for childhood ALL at IKZF1, with rs78396808 unique to non-European populations. We also demonstrated that the most significant DMP in the EWAS, CpG cg13344587 at gene ARID5B (P = 8.61 × 10−10), was entirely confounded by the ARID5B ALL risk SNP rs7090445. Our findings provide new insights into the functional pathways of ALL risk SNPs and the DNA methylation differences associated with risk of childhood ALL.

Background: Parental smoking is implicated in the etiology of acute lymphoblastic leukemia (ALL), the most common childhood cancer. We recently reported an association between an epigenetic biomarker of early-life tobacco smoke exposure at the AHRR gene and increased frequency of somatic gene deletions among ALL cases. Methods: Here, we further assess this association using two epigenetic biomarkers for maternal smoking during pregnancy—DNA methylation at AHRR CpG cg05575921 and a recently established polyepigenetic smoking score—in an expanded set of 482 B-cell ALL (B-ALL) cases in the California Childhood Leukemia Study with available Illumina 450K or MethylationEPIC array data. Multivariable Poisson regression models were used to test the associations between the epigenetic biomarkers and gene deletion numbers. Results: We found an association between DNA methylation at AHRR CpG cg05575921 and deletion number among 284 childhood B-ALL cases with MethylationEPIC array data, with a ratio of means (RM) of 1.31 [95% confidence interval (CI), 1.02–1.69] for each 0.1 β value reduction in DNA methylation, an effect size similar to our previous report in an independent set of 198 B-ALL cases with 450K array data [meta-analysis summary RM (sRM) = 1.32; 95% CI, 1.10–1.57]. The polyepigenetic smoking score was positively associated with gene deletion frequency among all 482 B-ALL cases (sRM = 1.31 for each 4-unit increase in score; 95% CI, 1.09–1.57). Conclusions: We provide further evidence that prenatal tobacco-smoke exposure may influence the generation of somatic copy-number deletions in childhood B-ALL. Impact: Analyses of deletion breakpoint sequences are required to further understand the mutagenic effects of tobacco smoke in childhood ALL.

The activation of monocytes and macrophages plays a central role in chronic inflammation associated with type 2 diabetes mellitus (T2D) and cardiovascular disease. However, the precise epigenetic mechanisms involved are unclear. From RNA-seq analysis we found dysregulated expression of several genes, including downregulation of JARID2 in human CD14+ monocytes obtained from T2D patients versus controls. JARID2 protein functions as a transcription repressor by facilitating the recruitment of polycomb repressor complex 2 and associated histone H3K27-trimethylation to target genes. RT-qPCR and RNA-seq analysis showed that high glucose +TNF-α (HT) , key pathological factors in diabetes, downregulated JARID2 in human monocytes along with genes involved in inflammation and cytokine related signaling networks. Western blotting also showed HT reduces JARID2 protein levels in human monocytes. Furthermore, JARID2 knockdown with siRNAs increased proliferation and enhanced LPS and HT induced inflammatory genes in human monocytes. ATAC-seq revealed that HT increased chromatin-accessibility at promoters of inflammatory genes like IL1B that were upregulated by JARID2 knockdown. Interestingly, microRNA-seq analysis revealed that CD14+ monocytes from T2D patients showed increased expression of several microRNAs with complementarity to the JARID2 3’-UTR, indicating cross-talk between epigenetic factors and microRNA-dependent mechanisms. Mouse orthologous Jarid2 was downregulated in macrophages from diabetic db/db mice and in macrophages from normal C57BL/6 mice treated with HT, similar to human JARID2. These results demonstrate that JARID2 plays a key role in chromatin remodeling and restraining inflammatory genes in monocytes, and its downregulation in T2D can augment inflammation. Further, understanding of epigenetic regulators like JARID2 can lead to novel therapies for chronic inflammatory diabetic complications. Disclosure M.A.Reddy: None. V.S.Tanwar: None. L.Zhang: None. L.L.Lanting: None. M.Wang: None. P.Pandey: None. Z.Chen: None. X.Wu: None. R.Natarajan: None. Funding National Institutes of Health (DK065073) to RN AR-DMRI Endowment Pilot grant to MAR

Acute lymphoblastic leukemia (ALL) is the most common cancer in children. Thirdhand smoke (THS) is the residual tobacco contamination that remains after the smoke clears. We investigated the effects of THS exposure in utero and during early life in a transgenic Cdkn2a knockout mouse model that is vulnerable to the development of leukemia/lymphoma. Female mice, and their offspring, were exposed from the first day of pregnancy to weaning. Plasma cytokines, body weight and hematologic parameters were measured in the offspring. To investigate THS exposure effects on the development of leukemia/lymphoma, bone marrow (BM) was collected from control and THS-exposed mice and transplanted into BM-ablated recipient mice, which were followed for tumor development for 1 year. We found that in utero and early-life THS exposure caused significant changes in plasma cytokine concentrations and in immune cell populations; changes appeared more pronounced in male mice. Spleen (SP) and BM B-cell populations were significantly lower in THS-exposed mice. We furthermore observed that THS exposure increased the leukemia/lymphoma-free survival in BM transplantation recipient mice, potentially caused by THS-induced B-cell toxicity. A trend towards increased solid tumors in irradiated mice reconstituted with THS-exposed BM stimulates the hypothesis that the immunosuppressive effects of in utero and early-life THS exposure might contribute to carcinogenesis by lowering the host defense to other toxic exposures. Our study adds to expanding evidence that THS exposure alters the immune system and that in utero and early-life developmental periods represent vulnerable windows of susceptibility for these effects.

Chronic myelocytic leukemia cell line K562 undergoes differentiation by phorbol esters to megakaryocytes and we have used this system to understand miRNA processing leading to isomiR generation. PMA treatment significantly altered the production of miRNA in K562 cells. Expression of 24.4% of miRNAs were found to be stimulated whereas expression of 10% miRNAs were inhibited by PMA treatment. Our results suggest that miRNA precursors are processed into isomiRs in a deterministic manner. The relative levels of different isomiRs of a miRNA remained mainly unchanged even after PMA treatment irrespective of overall changes in expression (either up-regulation or down-regulation). However, not all miRNAs behave in the same way, about 7% showed a variation of isomiR profiles after PMA treatment. Most of the later class of miRNAs were found to be oncogenic miRNAs. Further, it was also found that number of isomiRs was independent of abundance of a miRNA. Functional importance of different isomiRs was demonstrated using three different isomiRs of miR-22. Our results showed that different isomiRs could inhibit expression of targets genes with different efficiencies. Our study suggests that the heterogeneity of a miRNA population generated during processing is in general regulated and that variation in the generation of an isomiR can be a functionally important regulatory feature.

Background: Down syndrome (DS) is associated with an up to 30-fold increased risk of B-cell acute lymphoblastic leukemia (ALL), and DS-ALL patients have worse overall survival and increased long-term treatment-related health conditions compared with non-DS ALL patients. In a recent genome-wide association study of DS-ALL, established ALL genetic risk loci were associated with DS-ALL, with several single nucleotide polymorphisms (SNPs) conferring a larger effect on ALL risk in the context of DS than in euploidy. We performed an epigenome-wide association study (EWAS) to elucidate whether epigenetic differences at birth are associated with risk of subsequent DS-ALL. Methods: The DS-ALL Discovery Study included 147 DS-ALL cases and 198 DS controls from the International Study of Down Syndrome Acute Leukemia, with newborn dried bloodspots (DBS) obtained from California (n=326) and Washington state (n=19) biobanks. The DS-ALL Replication Study included 24 DS-ALL cases and 24 DS controls with newborn DBS from the Michigan Neonatal Biobank. DNA was isolated from DBS, bisulfite converted, and assayed using Illumina Infinium MethylationEPIC Beadchip genome-wide DNA methylation arrays. Raw data were processed using "minfi" and "noob" packages in R. Reference-based deconvolution of blood cell proportions was performed using the Identifying Optimal DNA methylation Libraries (IDOL) algorithm, using DNA methylation data from cord blood reference samples, to estimate proportions of B cells, T cells (CD4+ and CD8+), monocytes, granulocytes, natural killer cells, and nucleated red blood cells. We compared each cell type proportion between DS-ALL cases and DS controls using linear regression adjusting for sex, plate, and principal components (PCs) to account for genetic ancestry. To identify single CpG probes associated with DS-ALL risk, we performed a multiethnic EWAS of DS-ALL in each study using linear regression adjusting for sex, plate, and PCs related to: 1) cell-type proportions and 2) genetic ancestry. Differentially methylated regions (DMRs) were identified using DMRcate and comb-p methods. In the Discovery Study, genome-wide SNP array data were available for 131 cases and 130 controls, and data from targeted sequencing of somatic mutations in exons 2/3 of GATA1 were available for 184/198 DS controls. Results: Deconvolution of blood cell proportions in the DS-ALL Discovery Study showed significantly higher B cell proportions in newborns with DS who later developed ALL (mean=0.0128, sd=0.0151) compared with DS controls (mean=0.00826, sd=0.0115) (P=6.4x10 -4, coefficient=0.0052). A significantly higher B cell proportion at birth was also found in DS-ALL cases in the independent Replication Study (cases mean=0.048, sd=0.024; controls mean=0.039, sd=0.028; P=0.03, coefficient=0.015). In the Discovery Study, the B cell difference remained significant (P=5.8x10 -3) with a similar effect size (coefficient=0.0045) after removal of GATA1 mutation-positive DS controls (n=30). We also investigated whether DS-ALL risk SNPs at ARID5B, IKZF1, GATA3, and CDKN2A may confound the association, but the increased B cell proportions in DS-ALL remained significant and effect estimates slightly increased in SNP genotype-adjusted models (coefficient range:0.0055-0.0059). In the EWAS of DS-ALL, 9 CpGs reached epigenome-wide significance (P Conclusions: Increased B cell proportions in newborns with DS may be a risk factor for development of DS-ALL in childhood. This finding, based on DNA methylation data, requires confirmation using conventional cell count measures, and should be explored as a novel biomarker for ALL risk in the non-DS population. Single CpGs and DMRs associated with DS-ALL risk in our Discovery Study require further investigation, including in additional ALL case-control studies in DS and non-DS populations. Disclosures Ma: Celgene/Bristol Myers Squibb: Consultancy, Research Funding.

Objective Recently, many tumor sequencing studies have inferred and reported on mutational signatures, short nucleotide patterns at which particular somatic base substitutions appear more often. A number of signatures reflect biological processes in the patient and factors associated with cancer risk. Our goal is to infer mutational signatures appearing in colon cancer, a cancer for which environmental risk factors vary by cancer subtype, and compare the signatures to those in adult stem cells from normal colon. We also compare the mutational signatures to others in the literature. Results We apply a probabilistic mutation signature model to somatic mutations previously reported for six adult normal colon stem cells and 431 colon adenocarcinomas. We infer six mutational signatures in colon cancer, four being specific to tumors with hypermutation. Just two signatures explained the majority of mutations in the small number of normal aging colon samples. All six signatures are independently identified in a series of 295 Chinese colorectal cancers.

We propose a hierarchical latent Dirichlet allocation model (HiLDA) for characterizing somatic mutation data in cancer. The method allows us to infer mutational patterns and their relative frequencies in a set of tumor mutational catalogs and to compare the estimated frequencies between tumor sets. We apply our method to two datasets, one containing somatic mutations in colon cancer by the time of occurrence, before or after tumor initiation, and the second containing somatic mutations in esophageal cancer by sex, age, smoking status, and tumor site. In colon cancer, the relative frequencies of mutational patterns were found significantly associated with the time of occurrence of mutations. In esophageal cancer, the relative frequencies were significantly associated with the tumor site. Our novel method provides higher statistical power for detecting differences in mutational signatures.

Additional file 1. Additional figures and tables.

Regulation of gene expression in higher eukaryotes has undergone a paradigm shift from an essentially protein factor-based regulatory system to a complex network of different systems involving, not only proteins, but also non-coding RNAs. This has brought in different layers of control, resulting in the fine-tuning of the final expressed product. The major class of regulatory RNAs is microRNA. MiRNAs function through regulation of expression of their cognate target genes and since target genes can also have regulatory role, overall effect of miRNAs can be quite significant. Computational methods have been at the core of miRNA research. Starting from prediction of miRNA genes, expression profiling and target gene identification computational methods have been very useful at every stage. Moreover, theoretical modelling of gene regulatory networks has enhanced our understanding of the mechanism of miRNA function in different cellular processes.

Children with Down syndrome (DS), caused by constitutive trisomy of chromosome 21, have a 20-fold increased risk of acute lymphoblastic leukemia (DS-ALL) and 500-fold increased risk of acute megakaryoblastic leukemia (DS-AMKL). At least 10-15% of DS neonates are born with the pre-leukemic syndrome transient abnormal myelopoiesis (TAM). Trisomy 21 affects hematopoiesis and leukemia risk; however, the underlying mechanisms are not fully understood. Previous small-scale studies (sample N DNA was extracted from newborn dried bloodspots from 196 children born with DS and 442 non-DS controls from the California Biobank Program, and assayed using Illumina Infinium MethylationEPIC arrays containing >850,000 CpG probes. Data preprocessing was performed using the SeSAMe R package, and the conumee package confirmed trisomy 21 or euploidy in DS and non-DS subjects. Cell type deconvolution was performed, and blood cell proportions compared between DS and non-DS newborns using linear regression adjusting for covariates. ReFACTor principal components (PCs) were used to adjust for cell type heterogeneity and EPISTRUCTURE PCs to adjust for ancestry. Epigenome-wide association studies (EWAS) were performed overall, and stratified by ethnicity, using linear regression models adjusting for sex, plate, first six ReFACTor PCs, and first six ancestry-related PCs. Differentially methylated regions (DMRs) were evaluated using DMRcate and comb-p, with a consensus list of significant DMRs using the overlap. Deconvolution of blood cell proportions revealed highly significant (P Intriguingly, PC analysis and hierarchical clustering of DNA methylation data identified a subset of DS subjects (N=34/196, 17%) that clustered separately: these all had significantly higher nRBC proportions than other subjects and possibly represent DS neonates with TAM, although further investigation is needed for confirmation. Removing this cluster did not affect our main findings in our EWAS or DMR analyses. Constitutive trisomy 21 has profound effects on DNA methylation across the genome, in particular resulting in repression of known regulators of hematopoiesis including RUNX1 and FLI1. Our findings highlight potential mechanisms for the increased risk of both lymphoid and myeloid malignancies in children with DS. Disclosures No relevant conflicts of interest to declare.

Repeat elements identified in Heterorhabditis indica transcripts. (DOCX 16 kb)

Repeat elements identified in Heterorhabditis indica transcripts. (DOCX 16 kb)

Background Nematodes are the most numerous animals in the soil. Insect parasitic nematodes of the genus Heterorhabditis are capable of selectively seeking, infecting and killing their insect-hosts in the soil. The infective juvenile (IJ) stage of the Heterorhabditis nematodes is analogous to Caenorhabditis elegans dauer juvenile stage, which remains in ‘arrested development’ till it finds and infects a new insect-host in the soil. H. indica is the most prevalent species of Heterorhabditis in India. To understand the genes and molecular processes that govern the biology of the IJ stage, and to create a resource to facilitate functional genomics and genetic exploration, we sequenced the transcriptome of H. indica IJs. Results The de-novo sequence assembly using Velvet-Oases pipeline resulted in 13,593 unique transcripts at N50 of 1,371 bp, of which 53 % were annotated by blastx. H. indica transcripts showed higher orthology with parasitic nematodes as compared to free living nematodes. In-silico expression analysis showed 30 % of transcripts expressing with ≥100 FPKM value. All the four canonical dauer formation pathways like cGMP-PKG, insulin, dafachronic acid and TGF-β were active in the IJ stage. Several other signaling pathways were highly represented in the transcriptome. Twenty-four orthologs of C. elegans RNAi pathway effector genes were discovered in H. indica, including nrde-3 that is reported for the first time in any of the parasitic nematodes. An ortholog of C. elegans tol-1 was also identified. Further, 272 kinases belonging to 137 groups, and several previously unidentified members of important gene classes were identified. Conclusions We generated high-quality transcriptome sequence data from H. indica IJs for the first time. The transcripts showed high similarity with the parasitic nematodes, M. hapla, and A. suum as opposed to C. elegans, a species to which H. indica is more closely related. The high representation of transcripts from several signaling pathways in the IJs indicates that despite being a developmentally arrested stage; IJs are a hotbed of signaling and are actively interacting with their environment. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2510-z) contains supplementary material, which is available to authorized users.

Background Hypoxia is a critical aspect of the glioma microenvironment and has been associated with poor prognosis and resistance to various therapies. However, the mechanisms responsible for hypoxic survival of glioma cells remain unclear. Recent studies strongly suggest that microRNAs act as critical mediators of the hypoxic response. We thus hypothesized their prominent role in hypoxia resistance in glioblastoma (GBM) and aimed to identify those. Results With this study, we present the first detailed analysis of small RNA transcriptome of cell line U87MG, a grade IV glioma cell line, and its alteration under hypoxic condition. Based on deep sequencing and microarray data, we identify a set of hypoxia regulated microRNAs, with the miR-210-3p and its isomiRs showing highest induction in GBM cell lines U87MG and U251MG. We show miR-210-3p, miR-1275, miR-376c-3p, miR-23b-3p, miR-193a-3p and miR-145-5p to be up-regulated, while miR-92b-3p, miR-20a-5p, miR-10b-5p, miR-181a-2-3p and miR-185-5p are down-regulated by hypoxia. Interestingly, certain hypoxia-induced miRNAs are also known to be over-expressed in GBM tumors, suggesting that hypoxia may be one of the factors involved in establishing the miRNA signature of GBM. Transcription factor binding sites for Hypoxia inducible factor 1 A (HIF1A) were identified in the promoter region (5 kb upstream) of 30 hypoxia-induced miRNAs. HIF-1A over-expression and silencing studies show regulation of specific miRNAs, including miR-210-3p, to be HIF1A dependent. On the other hand, miR-210-3p leads to an increase in transcriptional activity of HIF and its target genes vascular endothelial growth factor (VEGF) and carbonic anhydrase 9 (CA9). MiR-210-3p levels were found to be high in GBM patient samples and showed good correlation with the known hypoxia markers CA9 and VEGF. We show that miR-210-3p promotes hypoxic survival and chemoresistance in GBM cells and targets a negative regulator of hypoxic response, HIF3A. Additionally, a total of 139 novel miRNAs were discovered by the analysis of deep sequencing data and three of these were found to be differentially expressed under hypoxia. Conclusions Overall, our study reveals a novel miRNA signature of hypoxia in GBM and suggests miR-210-3p to be an oncogenic player and a novel potential intrinsic marker of hypoxia in glioblastoma. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-686) contains supplementary material, which is available to authorized users.

MicroRNAs control cellular processes by regulating expression of their target genes. Here we report that neuro-epithelial transforming gene 1 (NET1) is a target of tumor suppressor microRNA 22 (miR-22). miR-22 is downregulated in peripheral blood mononuclear cells derived from chronic myeloid leukemia (CML) patients and in CML cell line K562. NET1 was identified as one of the targets of miR-22 using both in vitro and in vivo experiments. Either mutations or naturally occurring single-nucleotide polymorphisms in NET1 3′-UTR that map at the miR-22 binding site were found to affect binding of miR-22 to NET1 mRNA. Over expression of NET1 in K562 cells resulted in increased proliferation. However decreased proliferation and alteration in cell cycle were observed on either overexpression of miR-22 or knockdown of NET1 expression respectively. We also found that overexpression of miR-22 or NET1 knockdown inhibits actin fiber formation, probably by downregulation of NET1 as NET1 knockdown also resulted in depletion of actin fiber formation. We suggest that the oncogenic properties of CML cells are probably due to deregulated expression of NET1 as a result of altered expression of miR-22.

Background MicroRNAs (miRNAs) have been recognized as one of the key regulatory non-coding RNAs that are involved in a number of basic cellular processes. miRNA expression profiling helps to identify miRNAs that could serve as biomarkers. Next generation sequencing (NGS) platforms provide the most effective way of miRNA profiling, particularly as expression of different isoforms of miRNA (IsomiRs) can be estimated by NGS. Therefore, it is now possible to discern the overall complexity of miRNA populations that participate in gene regulatory networks. It is thus important to consider different isoforms of miRNA as part of total profiling in order to understand all aspects of the biology of miRNAs. Results Here next generation sequencing data of small RNAs derived from normal peripheral blood mononuclear cells (PBMC) and Chronic myeloid leukemia (CML) patients has been used to generate miRNA profiles using a computation pipeline which can identify isomiRs that are natural variants of mature miRNAs. IsomiR profiles have been generated for all the 5p and 3p miRNAs (previously known as major mature miRNA and minor or miRNA*) and the data has been presented as a composite total miRNA transcriptome. The results indicated that the most abundant isomiR sequence of about 68% miRNAs, did not match the reference miRNA sequence as entered in the miRBase and that there is a definite pattern in relative concentration of different isomiRs derived from same precursors. Finally, a total of 17 potential novel miRNA sequences were identified suggesting that there are still some new miRNAs yet to be discovered. Conclusions Inclusion of different isoforms provides a detailed miRnome of a cell type or tissues. Availability of miRnome will be useful for finding biomarkers of different cell types and disease states. Our results also indicate that the relative expression levels of different isoforms of a miRNA are likely to be dynamic and may change with respect to changes in the cell or differentiation status.