Accelerated epigenetic aging in newborns with Down syndrome.

Accelerated aging is a hallmark of Down syndrome (DS), with adults experiencing early-onset Alzheimer's disease and premature aging of the skin, hair, and immune and endocrine systems. Accelerated epigenetic aging has been found in the blood and brain tissue of adults with DS but when premature aging in DS begins remains unknown. We investigated whether accelerated aging in DS is already detectable in blood at birth. We assessed the association between age acceleration and DS using five epigenetic clocks in 346 newborns with DS and 567 newborns without DS using Illumina MethylationEPIC DNA methylation array data. We compared two epigenetic aging clocks (DNAmSkinBloodClock and pan-tissue DNAmAge) and three epigenetic gestational age clocks (Haftorn, Knight, and Bohlin) between DS and non-DS newborns using linear regression adjusting for observed age, sex, batch, deconvoluted blood cell proportions, and genetic ancestry. Targeted sequencing of GATA1 was performed in a subset of 184 newborns with DS to identify somatic mutations associated with transient abnormal myelopoiesis. DS was significantly associated with increased DNAmSkinBloodClock (effect estimate = 0.2442, p < 0.0001), with an epigenetic age acceleration of 244 days in newborns with DS after adjusting for potential confounding factors (95% confidence interval: 196-292 days). We also found evidence of epigenetic age acceleration associated with somatic GATA1 mutations among newborns with DS (p = 0.015). DS was not associated with epigenetic gestational age acceleration. We demonstrate that accelerated epigenetic aging in the blood of DS patients begins prenatally, with implications for the pathophysiology of immunosenescence and other aging-related traits in DS.

Investigating DNA methylation as a mediator of genetic risk in childhood acute lymphoblastic leukemia.

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.

Development of a Droplet Digital™ PCR DNA methylation detection and quantification assay of prenatal tobacco exposure.

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.

Lysine acetylation regulates the interaction between proteins and membranes.

Lysine acetylation regulates the function of soluble proteins in vivo, yet it remains largely unexplored whether lysine acetylation regulates membrane protein function. Here, we use bioinformatics, biophysical analysis of recombinant proteins, live-cell fluorescent imaging and genetic manipulation of Drosophila to explore lysine acetylation in peripheral membrane proteins. Analysis of 50 peripheral membrane proteins harboring BAR, PX, C2, or EHD membrane-binding domains reveals that lysine acetylation predominates in membrane-interaction regions. Acetylation and acetylation-mimicking mutations in three test proteins, amphiphysin, EHD2, and synaptotagmin1, strongly reduce membrane binding affinity, attenuate membrane remodeling in vitro and alter subcellular localization. This effect is likely due to the loss of positive charge, which weakens interactions with negatively charged membranes. In Drosophila, acetylation-mimicking mutations of amphiphysin cause severe disruption of T-tubule organization and yield a flightless phenotype. Our data provide mechanistic insights into how lysine acetylation regulates membrane protein function, potentially impacting a plethora of membrane-related processes.

Epigenetic Biomarkers of Prenatal Tobacco Smoke Exposure Are Associated with Gene Deletions in Childhood Acute Lymphoblastic Leukemia.

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.

In utero and early-life exposure to thirdhand smoke causes profound changes to the immune system.

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.

The genome-wide impact of trisomy 21 on DNA methylation and its implications for hematopoiesis.

Down syndrome is associated with genome-wide perturbation of gene expression, which may be mediated by epigenetic changes. We perform an epigenome-wide association study on neonatal bloodspots comparing 196 newborns with Down syndrome and 439 newborns without Down syndrome, adjusting for cell-type heterogeneity, which identifies 652 epigenome-wide significant CpGs (P < 7.67 × 10-8) and 1,052 differentially methylated regions. Differential methylation at promoter/enhancer regions correlates with gene expression changes in Down syndrome versus non-Down syndrome fetal liver hematopoietic stem/progenitor cells (P < 0.0001). The top two differentially methylated regions overlap RUNX1 and FLI1, both important regulators of megakaryopoiesis and hematopoietic development, with significant hypermethylation at promoter regions of these two genes. Excluding Down syndrome newborns harboring preleukemic GATA1 mutations (N = 30), identified by targeted sequencing, has minimal impact on the epigenome-wide association study results. Down syndrome has profound, genome-wide effects on DNA methylation in hematopoietic cells in early life, which may contribute to the high frequency of hematological problems, including leukemia, in children with Down syndrome.

iMutSig: a web application to identify the most similar mutational signature using shiny.

There are two frameworks for characterizing mutational signatures which are commonly used to describe the nucleotide patterns that arise from mutational processes. Estimated mutational signatures from fitting these two methods in human cancer can be found online, in the Catalogue Of Somatic Mutations In Cancer (COSMIC) website or a GitHub repository. The two frameworks make differing assumptions regarding independence of base pairs and for that reason may produce different results. Consequently, there is a need to compare and contrast the results of the two methods, but no such tool currently exists. In this paper, we provide a simple and intuitive interface that allows such comparisons to be easily performed. When using our software, the user may download published mutational signatures of either type. Mutational signatures from the pmsignature data source are expanded to probabilistic vectors of 96-possible mutation types, the same model specification used by COSMIC, and then compared to COSMIC signatures. Cosine similarity measures the extent of signature similarity. iMutSig provides a simple and user-friendly web application allowing researchers to compare signatures from COSMIC to those from pmsignature, and vice versa. Furthermore, iMutSig allows users to input a self-defined mutational signature and examine its similarity to published signatures from both data sources. iMutSig is accessible online and source code is available for download on GitHub.

Mutational signatures in colon cancer.

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.

HiLDA: a statistical approach to investigate differences in mutational signatures.

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.

De novo genome and transcriptome analyses provide insights into the biology of the trematode human parasite Fasciolopsis buski.

Many trematode parasites cause infection in humans and are thought to be a major public health problem. Their ecological diversity in different regions provides challenging questions on evolution of these organisms. In this report, we perform transcriptome analysis of the giant intestinal fluke, Fasciolopsis buski, using next generation sequencing technology. Short read sequences derived from polyA containing RNA of this organism were assembled into 30,677 unigenes that led to the annotation of 12,380 genes. Annotation of the assembled transcripts enabled insight into processes and pathways in the intestinal fluke, such as RNAi pathway and energy metabolism. The expressed kinome of the organism was characterized by identifying all protein kinases. A rough draft genome assembly for Fasciolopsis buski is also reported herewith with SRA accessions for crosschecking the findings in the analyzed transcriptome data. Transcriptome data also helped us to identify some of the expressed transposable elements. Though many Long Interspersed elements (LINEs) were identified, only two Short Interspersed Elements (SINEs) were visible. Overall transcriptome and draft genome analysis of F. buski helped us to characterize some of its important biological characteristics and provided enormous resources for development of a suitable diagnostic system and anti-parasitic therapeutic molecules.

The 17-residue-long N terminus in huntingtin controls stepwise aggregation in solution and on membranes via different mechanisms.

Aggregation of huntingtin protein arising from expanded polyglutamine (polyQ) sequences in the exon-1 region of mutant huntingtin plays a central role in the pathogenesis of Huntington's disease. The huntingtin aggregation pathways are of therapeutic and diagnostic interest, but obtaining critical information from the physiologically relevant htt exon-1 (Httex1) protein has been challenging. Using biophysical techniques and an expression and purification protocol that generates clean, monomeric Httex1, we identified and mapped three distinct aggregation pathways: 1) unseeded in solution; 2) seeded in solution; and 3) membrane-mediated. In solution, aggregation proceeded in a highly stepwise manner, in which the individual domains (N terminus containing 17 amino acids (N17), polyQ, and proline-rich domain (PRD)) become ordered at very different rates. The aggregation was initiated by an early oligomer requiring a pathogenic, expanded Gln length and N17 α-helix formation. In the second phase, ß-sheet forms in the polyQ. The slowest step is the final structural maturation of the PRD. This stepwise mechanism could be bypassed by seeding, which potently accelerated aggregation and was a prerequisite for prion-like spreading in vivo Remarkably, membranes could catalyze aggregation even more potently than seeds, in a process that caused significant membrane damage. The N17 governed membrane-mediated aggregation by anchoring Httex1 to the membrane, enhancing local concentration and promoting collision via two-dimensional diffusion. Considering its central roles in solution and in membrane-mediated aggregation, the N17 represents an attractive target for inhibiting multiple pathways. Our approach should help evaluate such inhibitors and identify diagnostic markers for the misfolded forms identified here.

IsomiR processing during differentiation of myelogenous leukemic cell line K562 by phorbol ester PMA.

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.

A transcriptomic insight into the infective juvenile stage of the insect parasitic nematode, Heterorhabditis indica.

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.

Hypoxic signature of microRNAs in glioblastoma: insights from small RNA deep sequencing.

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.

miR-22 regulates expression of oncogenic neuro-epithelial transforming gene 1, NET1.

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.

Analysis of the microRNA transcriptome and expression of different isomiRs in human peripheral blood mononuclear cells.

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.