Pankaj Chopra

Age-associated DNA methylation in pediatric populations

DNA methylation (DNAm) plays diverse roles in human biology, but this dynamic epigenetic mark remains far from fully characterized. Although earlier studies uncovered loci that undergo age-associated DNAm changes in adults, little is known about such changes during childhood. Despite profound DNAm plasticity during embryogenesis, monozygotic twins show indistinguishable childhood methylation, suggesting that DNAm is highly coordinated throughout early development. Here we examine the methylation of 27,578 CpG dinucleotides in peripheral blood DNA from a cross-sectional study of 398 boys, aged 3-17 yr, and find significant age-associated changes in DNAm at 2078 loci. These findings correspond well with pyrosequencing data and replicate in a second pediatric population (N = 78). Moreover, we report a deficit of age-related loci on the X chromosome, a preference for specific nucleotides immediately surrounding the interrogated CpG dinucleotide, and a primary association with developmental and immune ontological functions. Meta-analysis (N = 1158) with two adult populations reveals that despite a significant overlap of age-associated loci, most methylation changes do not follow a lifelong linear pattern due to a threefold to fourfold higher rate of change in children compared with adults; consequently, the vast majority of changes are more accurately modeled as a function of logarithmic age. We therefore conclude that age-related DNAm changes in peripheral blood occur more rapidly during childhood and are imperfectly accounted for by statistical corrections that are linear in age, further suggesting that future DNAm studies should be matched closely for age.

Analysis of FMRP mRNA target datasets reveals highly associated mRNAs mediated by G-quadruplex structures formed via clustered WGGA sequences

Fragile X syndrome, a common cause of intellectual disability and a well-known cause of autism spectrum disorder, is the result of loss or dysfunction of fragile X mental retardation protein (FMRP), a highly selective RNA-binding protein and translation regulator. A major research priority has been the identification of the mRNA targets of FMRP, particularly as recent studies suggest an excess of FMRP targets among genes implicated in idiopathic autism and schizophrenia. Several large-scale studies have attempted to identify mRNAs bound by FMRP through several methods, each generating a list of putative target genes, leading to distinct hypotheses by which FMRP recognizes its targets; namely, by RNA structure or sequence. However, no in depth analyses have been performed to identify the level of consensus among the studies. Here, we analyze four large FMRP target datasets to generate high-confidence consensus lists, and examine all datasets for sequence elements within the target RNAs to validate reported FMRP binding motifs (GACR, ACUK and WGGA). We found GACR to be highly enriched in FMRP datasets, while ACUK was not. The WGGA pattern was modestly enriched in several, but not all datasets. The previous association between FMRP and G-quadruplexes prompted the analysis of the distribution of WGGA in the target genes. Consistent with the requirements for G-quadruplex formation, we observed highly clustered WGGA motifs in FMRP targets compared with other genes, implicating both RNA structure and sequence in the recognition motif of FMRP. In addition, we generate a list of the top 40 FMRP targets associated with FXS-related phenotypes.

Improving Cancer Classification Accuracy Using Gene Pairs

Recent studies suggest that the deregulation of pathways, rather than individual genes, may be critical in triggering carcinogenesis. The pathway deregulation is often caused by the simultaneous deregulation of more than one gene in the pathway. This suggests that robust gene pair combinations may exploit the underlying bio-molecular reactions that are relevant to the pathway deregulation and thus they could provide better biomarkers for cancer, as compared to individual genes. In order to validate this hypothesis, in this paper, we used gene pair combinations, called doublets, as input to the cancer classification algorithms, instead of the original expression values, and we showed that the classification accuracy was consistently improved across different datasets and classification algorithms. We validated the proposed approach using nine cancer datasets and five classification algorithms including Prediction Analysis for Microarrays (PAM), C4.5 Decision Trees (DT), Naive Bayesian (NB), Support Vector Machine (SVM), and k-Nearest Neighbor (k-NN).

Genome-wide analysis validates aberrant methylation in fragile X syndrome is specific to the FMR1 locus.

BackgroundFragile X syndrome (FXS) is a common form of inherited intellectual disability caused by an expansion of CGG repeats located in the 5′ untranslated region (UTR) of the FMR1 gene, which leads to hypermethylation and silencing of this locus. Although a dramatic increase in DNA methylation of the FMR1 full mutation allele is well documented, the extent to which these changes affect DNA methylation throughout the rest of the genome has gone unexplored.MethodsHere we examined genome-wide methylation in both peripheral blood (N = 62) and induced pluripotent stem cells (iPSCs; N = 10) from FXS individuals and controls.ResultsWe not only found the expected significant DNA methylation differences in the FMR1 promoter and 5′ UTR, we also saw that these changes inverse in the FMR1 gene body. Importantly, we found no other differentially methylated loci throughout the remainder of the genome, indicating the aberrant methylation of FMR1 in FXS is locus-specific.ConclusionsThis study provides a comprehensive methylation profile of FXS and helps refine our understanding of the mechanisms behind FMR1 silencing.

Differentially methylated plasticity genes in the amygdala of young primates are linked to anxious temperament, an at risk phenotype for anxiety and depressive disorders.

Children with an anxious temperament (AT) are at a substantially increased risk to develop anxiety and depression. The young rhesus monkey is ideal for studying the origin of human AT because it shares with humans the genetic, neural, and phenotypic underpinnings of complex social and emotional functioning. Heritability, functional imaging, and gene expression studies of AT in young monkeys revealed that the central nucleus of the amygdala (Ce) is a key environmentally sensitive substrate of this at risk phenotype. Because epigenetic marks (e.g., DNA methylation) can be modulated by environmental stimuli, these data led us to hypothesize a role for DNA methylation in the development of AT. To test this hypothesis, we used reduced representation bisulfite sequencing to examine the cross-sectional genome-wide methylation levels in the Ce of 23 age-matched monkeys (1.3 ± 0.2 years) phenotyped for AT. Because AT reflects a continuous trait-like variable, we used an analytical approach that is consistent with this biology to identify genes in the Ce with methylation patterns that predict AT. Expression data from the Ce of these same monkeys were then used to find differentially methylated candidates linked to altered gene regulation. Two genes particularly relevant to the AT phenotype were BCL11A and JAG1. These transcripts have well-defined roles in neurodevelopmental processes, including neurite arborization and the regulation of neurogenesis. Together, these findings represent a critical step toward understanding the effects of early environment on the neuromolecular mechanisms that underlie the risk to develop anxiety and depressive disorders.

Identification of consensus binding sites clarifies FMRP binding determinants

Fragile X mental retardation protein (FMRP) is a multifunctional RNA-binding protein with crucial roles in neuronal development and function. Efforts aimed at elucidating how FMRP target mRNAs are selected have produced divergent sets of target mRNA and putative FMRP-bound motifs, and a clear understanding of FMRPs binding determinants has been lacking. To clarify FMRPs binding to its target mRNAs, we produced a shared dataset of FMRP consensus binding sequences (FCBS), which were reproducibly identified in two published FMRP CLIP sequencing datasets. This comparative dataset revealed that of the various sequence and structural motifs that have been proposed to specify FMRP binding, the short sequence motifs TGGA and GAC were corroborated, and a novel TAY motif was identified. In addition, the distribution of the FCBS set demonstrates that FMRP preferentially binds to the coding region of its targets but also revealed binding along 3′ UTRs in a subset of target mRNAs. Beyond probing these putative motifs, the FCBS dataset of reproducibly identified FMRP binding sites is a valuable tool for investigating FMRP targets and function.

Genome-wide alterations in hippocampal 5-hydroxymethylcytosine links plasticity genes to acute stress

Environmental stress is among the most important contributors to increased susceptibility to develop psychiatric disorders, including anxiety and post-traumatic stress disorder. While even acute stress alters gene expression, the molecular mechanisms underlying these changes remain largely unknown. 5-hydroxymethylcytosine (5hmC) is a novel environmentally sensitive DNA modification that is highly enriched in post-mitotic neurons and is associated with active transcription of neuronal genes. Recently, we found a hippocampal increase of 5hmC in the glucocorticoid receptor gene (Nr3c1) following acute stress, warranting a deeper investigation of stress-related 5hmC levels. Here we used an established chemical labeling and affinity purification method coupled with high-throughput sequencing technology to generate the first genome-wide profile of hippocampal 5hmC following exposure to acute restraint stress and a one-hour recovery. This approach found a genome-wide disruption in 5hmC associated with acute stress response, primarily in genic regions, and identified known and potentially novel stress-related targets that have a significant enrichment for neuronal ontological functions. Integration of these data with hippocampal gene expression data from these same mice found stress-related hydroxymethylation correlated to altered transcript levels and sequence motif predictions indicated that 5hmC may function by mediating transcription factor binding to these transcripts. Together, these data reveal an environmental impact on this newly discovered epigenetic mark in the brain and represent a critical step toward understanding stress-related epigenetic mechanisms that alter gene expression and can lead to the development of psychiatric disorders.

Reactivation of FMR1 by CRISPR/Cas9-Mediated Deletion of the Expanded CGG-Repeat of the Fragile X Chromosome

Fragile X syndrome (FXS) is a common cause of intellectual disability that is most often due to a CGG-repeat expansion mutation in the FMR1 gene that triggers epigenetic gene silencing. Epigenetic modifying drugs can only transiently and modestly induce FMR1 reactivation in the presence of the elongated CGG repeat. As a proof-of-principle, we excised the expanded CGG-repeat in both somatic cell hybrids containing the human fragile X chromosome and human FXS iPS cells using the CRISPR/Cas9 genome editing. We observed transcriptional reactivation in approximately 67% of the CRISPR cut hybrid colonies and in 20% of isolated human FXS iPSC colonies. The reactivated cells produced FMRP and exhibited a decline in DNA methylation at the FMR1 locus. These data demonstrate the excision of the expanded CGG-repeat from the fragile X chromosome can result in FMR1 reactivation.

PEMapper and PECaller provide a simplified approach to whole-genome sequencing

Significance PEMapper and PECaller are paired software programs that simplify mapping and variant calling for whole-genome datasets. Whole-genome sequencing data are fast becoming the most natural dataset for all genetic studies. Analysis tools for data at this scale are essential. This manuscript describes tools, which solve the challenges of data analysis at whole-genome scale, using an approach involving 16-mer mapping and SNP calling based on a Pólya–Eggenberger distribution for SNP genotypes. We show that our software package is faster (cheaper to run), uses much less disk space (cheaper to store results), requires no previous knowledge of existing genetic variation (easier to deploy to nonhuman species), and achieves calling results that are as good as Genome Analysis Toolkit best practices. The analysis of human whole-genome sequencing data presents significant computational challenges. The sheer size of datasets places an enormous burden on computational, disk array, and network resources. Here, we present an integrated computational package, PEMapper/PECaller, that was designed specifically to minimize the burden on networks and disk arrays, create output files that are minimal in size, and run in a highly computationally efficient way, with the single goal of enabling whole-genome sequencing at scale. In addition to improved computational efficiency, we implement a statistical framework that allows for a base by base error model, allowing this package to perform as well or better than the widely used Genome Analysis Toolkit (GATK) in all key measures of performance on human whole-genome sequences.

Sex-specific hippocampal 5-hydroxymethylcytosine is disrupted in response to acute stress.

Environmental stress is among the most important contributors to increased susceptibility to develop psychiatric disorders. While it is well known that acute environmental stress alters gene expression, the molecular mechanisms underlying these changes remain largely unknown. 5-hydroxymethylcytosine (5hmC) is a novel environmentally sensitive epigenetic modification that is highly enriched in neurons and is associated with active neuronal transcription. Recently, we reported a genome-wide disruption of hippocampal 5hmC in male mice following acute stress that was correlated to altered transcript levels of genes in known stress related pathways. Since sex-specific endocrine mechanisms respond to environmental stimulus by altering the neuronal epigenome, we examined the genome-wide profile of hippocampal 5hmC in female mice following exposure to acute stress and identified 363 differentially hydroxymethylated regions (DhMRs) linked to known (e.g., Nr3c1 and Ntrk2) and potentially novel genes associated with stress response and psychiatric disorders. Integration of hippocampal expression data from the same female mice found stress-related hydroxymethylation correlated to altered transcript levels. Finally, characterization of stress-induced sex-specific 5hmC profiles in the hippocampus revealed 778 sex-specific acute stress-induced DhMRs some of which were correlated to altered transcript levels that produce sex-specific isoforms in response to stress. Together, the alterations in 5hmC presented here provide a possible molecular mechanism for the adaptive sex-specific response to stress that may augment the design of novel therapeutic agents that will have optimal effectiveness in each sex.