The Encyclopedia of DNA elements (ENCODE): data portal update

Replication timing networks reveal a link between transcription regulatory circuits and replication timing control

DNA replication occurs in a defined temporal order known as the replication timing (RT) program and is regulated during development, coordinated with 3D genome organization and transcriptional activity. However, transcription and RT are not sufficiently coordinated to predict each other, suggesting an indirect relationship. Here, we exploit genome-wide RT profiles from 15 human cell types and intermediate differentiation stages derived from human embryonic stem cells to construct different types of RT regulatory networks. First, we constructed networks based on the coordinated RT changes during cell fate commitment to create highly complex RT networks composed of thousands of interactions that form specific functional subnetwork communities. We also constructed directional regulatory networks based on the order of RT changes within cell lineages, and identified master regulators of differentiation pathways. Finally, we explored relationships between RT networks and transcriptional regulatory networks (TRNs) by combining them into more complex circuitries of composite and bipartite networks. Results identified novel trans interactions linking transcription factors that are core to the regulatory circuitry of each cell type to RT changes occurring in those cell types. These core transcription factors were found to bind cooperatively to sites in the affected replication domains, providing provocative evidence that they constitute biologically significant directional interactions. Our findings suggest a regulatory link between the establishment of cell-type-specific TRNs and RT control during lineage specification.

The lncRNA Locus Handsdown Regulates Cardiac Gene Programs and Is Essential for Early Mouse Development

Precisely controlled gene regulatory networks are required during embryonic development to give rise to various structures, including those of the cardiovascular system. Long non-coding RNA (lncRNA) loci are known to be important regulators of these genetic programs. We have identified a novel and essential lncRNA locus Handsdown (Hdn), active in early heart cells, and show by genetic inactivation that it is essential for murine development. Hdn displays haploinsufficiency for cardiac development as Hdn-heterozygous adult mice exhibit hyperplasia in the right ventricular wall. Transcriptional activity of the Hdn locus, independent of its RNA, suppresses its neighboring gene Hand2. We reveal a switch in a topologically associated domain in differentiation of the cardiac lineage, allowing the Hdn locus to directly interact with regulatory elements of the Hand2 locus.

DNA methylation fingerprint of monozygotic twins and their singleton sibling with intellectual disability carrying a novel KDM5C mutation

Mutations in KDM5C (lysine (K)-specific demethylase 5C) were causally associated with up to 3% of X-linked intellectual disability (ID) in males. By exome and Sanger sequencing, a novel frameshift KDM5C variant, predicted to eliminate the JmjC catalytic domain from the protein, was identified in two monozygotic twins and their older brother, which was inherited from their clinically normal mother, who had completely skewed X-inactivation. DNA methylation (DNAm) data were evaluated using the Illumina 450 K Methylation Beadchip arrays. Comparison of methylation levels between the three patients and male controls identified 399 differentially methylated CpG sites, which were enriched among those CpG sites modulated during brain development. Most of them were hypomethylated (72%), and located mainly in shores, whereas the hypermethylated CpGs were more represented in open sea regions. The DNAm changes did not differ between the monozygotic twins nor between them and their older sibling, all presenting a global hypomethylation, similar to other studies that associated DNA methylation changes to different KDM5C mutations. The 38 differentially methylated regions (DMRs) were enriched for H3K4me3 marks identified in developing brains. The remarkable similarity between the methylation changes in the monozygotic twins and their older brother is indicative that these epigenetic changes were mostly driven by the KDM5C mutation.

RNA sequencing for research and diagnostics in clinical oncology

Molecular diagnostics is becoming one of the major drivers of personalized oncology. With hundreds of different approved anticancer drugs and regimens of their administration, selecting the proper treatment for a patient is at least nontrivial task. This is especially sound for the cases of recurrent and metastatic cancers where the standard lines of therapy failed. Recent trials demonstrated that mutation assays have a strong limitation in personalized selection of therapeutics, consequently, most of the drugs cannot be ranked and only a small percentage of patients can benefit from the screening. Other approaches are, therefore, needed to address a problem of finding proper targeted therapies. The analysis of RNA expression (transcriptomic) profiles presents a reasonable solution because transcriptomics stands a few steps closer to tumor phenotype than the genome analysis. Several recent studies pioneered using transcriptomics for practical oncology and showed truly encouraging clinical results. The possibility of directly measuring of expression levels of molecular drugs' targets and profiling activation of the relevant molecular pathways enables personalized prioritizing for all types of molecular-targeted therapies. RNA sequencing is the most robust tool for the high throughput quantitative transcriptomics. Its use, potentials, and limitations for the clinical oncology will be reviewed here along with the technical aspects such as optimal types of biosamples, RNA sequencing profile normalization, quality controls and several levels of data analysis.

GWAS4D: multidimensional analysis of context-specific regulatory variant for human complex diseases and traits

Genome-wide association studies have generated over thousands of susceptibility loci for many human complex traits, and yet for most of these associations the true causal variants remain unknown. Tissue/cell type-specific prediction and prioritization of non-coding regulatory variants will facilitate the identification of causal variants and underlying pathogenic mechanisms for particular complex diseases and traits. By leveraging recent large-scale functional genomics/epigenomics data, we develop an intuitive web server, GWAS4D (http://mulinlab.tmu.edu.cn/gwas4d or http://mulinlab.org/gwas4d), that systematically evaluates GWAS signals and identifies context-specific regulatory variants. The updated web server includes six major features: (i) updates the regulatory variant prioritization method with our new algorithm; (ii) incorporates 127 tissue/cell type-specific epigenomes data; (iii) integrates motifs of 1480 transcriptional regulators from 13 public resources; (iv) uniformly processes Hi-C data and generates significant interactions at 5 kb resolution across 60 tissues/cell types; (v) adds comprehensive non-coding variant functional annotations; (vi) equips a highly interactive visualization function for SNP-target interaction. Using a GWAS fine-mapped set for 161 coronary artery disease risk loci, we demonstrate that GWAS4D is able to efficiently prioritize disease-causal regulatory variants.

Glial response to hypoxia in mutants of NPAS1/3 homolog Trachealess through Wg signaling to modulate synaptic bouton organization

Synaptic structure and activity are sensitive to environmental alterations. Modulation of synaptic morphology and function is often induced by signals from glia. However, the process by which glia mediate synaptic responses to environmental perturbations such as hypoxia remains unknown. Here, we report that, in the mutant for Trachealess (Trh), the Drosophila homolog for NPAS1 and NPAS3, smaller synaptic boutons form clusters named satellite boutons appear at larval neuromuscular junctions (NMJs), which is induced by the reduction of internal oxygen levels due to defective tracheal branches. Thus, the satellite bouton phenotype in the trh mutant is suppressed by hyperoxia, and recapitulated in wild-type larvae raised under hypoxia. We further show that hypoxia-inducible factor (HIF)-1α/Similar (Sima) is critical in mediating hypoxia-induced satellite bouton formation. Sima upregulates the level of the Wnt/Wingless (Wg) signal in glia, leading to reorganized microtubule structures within presynaptic sites. Finally, hypoxia-induced satellite boutons maintain normal synaptic transmission at the NMJs, which is crucial for coordinated larval locomotion.

Disruption of chromatin organisation causes MEF2C gene overexpression in intellectual disability: a case report

BACKGROUND

Balanced structural variants are mostly described in disease with gene disruption or subtle rearrangement at breakpoints.

CASE PRESENTATION

Here we report a patient with mild intellectual deficiency who carries a de novo balanced translocation t(3;5). Breakpoints were fully explored by microarray, Array Painting and Sanger sequencing. No gene disruption was found but the chromosome 5 breakpoint was localized 228-kb upstream of the MEF2C gene. The predicted Topologically Associated Domains analysis shows that it contains only the MEF2C gene and a long non-coding RNA LINC01226. RNA studies looking for MEF2C gene expression revealed an overexpression of MEF2C in the lymphoblastoid cell line of the patient.

CONCLUSIONS

Pathogenicity of MEF2C overexpression is still unclear as only four patients with mild intellectual deficiency carrying 5q14.3 microduplications containing MEF2C are described in the literature. The microduplications in these individuals also contain other genes expressed in the brain. The patient presented the same phenotype as 5q14.3 microduplication patients. We report the first case of a balanced translocation leading to an overexpression of MEF2C similar to a functional duplication.

Lung Cancer Risk in Never-Smokers of European Descent is Associated With Genetic Variation in the 5p15.33 TERT-CLPTM1Ll Region

INTRODUCTION

Inherited susceptibility to lung cancer risk in never-smokers is poorly understood. The major reason for this gap in knowledge is that this disease is relatively uncommon (except in Asians), making it difficult to assemble an adequate study sample. In this study we conducted a genome-wide association study on the largest, to date, set of European-descent never-smokers with lung cancer.

METHODS

We conducted a two-phase (discovery and replication) genome-wide association study in never-smokers of European descent. We further augmented the sample by performing a meta-analysis with never-smokers from the recent OncoArray study, which resulted in a total of 3636 cases and 6295 controls. We also compare our findings with those in smokers with lung cancer.

RESULTS

We detected three genome-wide statistically significant single nucleotide polymorphisms rs31490 (odds ratio [OR]: 0.769, 95% confidence interval [CI]: 0.722-0.820; p value 5.31 × 10-16), rs380286 (OR: 0.770, 95% CI: 0.723-0.820; p value 4.32 × 10-16), and rs4975616 (OR: 0.778, 95% CI: 0.730-0.829; p value 1.04 × 10-14). All three mapped to Chromosome 5 CLPTM1L-TERT region, previously shown to be associated with lung cancer risk in smokers and in never-smoker Asian women, and risk of other cancers including breast, ovarian, colorectal, and prostate.

CONCLUSIONS

We found that genetic susceptibility to lung cancer in never-smokers is associated to genetic variants with pan-cancer risk effects. The comparison with smokers shows that top variants previously shown to be associated with lung cancer risk only confer risk in the presence of tobacco exposure, underscoring the importance of gene-environment interactions in the etiology of this disease.

U2AF65 assemblies drive sequence-specific splice site recognition

The essential splicing factor U2AF65 is known to help anchoring U2 snRNP at the branch site. Its C-terminal UHM domain interacts with ULM motifs of SF3b155, an U2 snRNP protein. Here, we report a cooperative binding of U2AF65 and the related protein CAPERα to the multi-ULM domain of SF3b155. In addition, we show that the RS domain of U2AF65 drives a liquid-liquid phase separation that is amplified by intronic RNA with repeated pyrimidine tracts. In cells, knockdown of either U2AF65 or CAPERα improves the inclusion of cassette exons that are preceded by such repeated pyrimidine-rich motifs. These results support a model in which liquid-like assemblies of U2AF65 and CAPERα on repetitive pyrimidine-rich RNA sequences are driven by their RS domains, and facilitate the recruitment of the multi-ULM domain of SF3b155. We anticipate that posttranslational modifications and proteins recruited in dynamical U2AF65 and CAPERα condensates may further contribute to the complex mechanisms leading to specific splice site choice that occurs in cells.

Mapping cis-regulatory chromatin contacts in neural cells links neuropsychiatric disorder risk variants to target genes

Mutations in gene regulatory elements have been associated with a wide range of complex neuropsychiatric disorders. However, due to their cell-type specificity and difficulties in characterizing their regulatory targets, the ability to identify causal genetic variants has remained limited. To address these constraints, we perform an integrative analysis of chromatin interactions, open chromatin regions and transcriptomes using promoter capture Hi-C, assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and RNA sequencing, respectively, in four functionally distinct neural cell types: induced pluripotent stem cell (iPSC)-induced excitatory neurons and lower motor neurons, iPSC-derived hippocampal dentate gyrus-like neurons and primary astrocytes. We identify hundreds of thousands of long-range cis-interactions between promoters and distal promoter-interacting regions, enabling us to link regulatory elements to their target genes and reveal putative processes that are dysregulated in disease. Finally, we validate several promoter-interacting regions by using clustered regularly interspaced short palindromic repeats (CRISPR) techniques in human excitatory neurons, demonstrating that CDK5RAP3, STRAP and DRD2 are transcriptionally regulated by physically linked enhancers.

Genome-Wide Association Study of Metabolic Syndrome Reveals Primary Genetic Variants at CETP Locus in Indians

Indians, a rapidly growing population, constitute vast genetic heterogeneity to that of Western population; however they have become a sedentary population in past decades due to rapid urbanization ensuing in the amplified prevalence of metabolic syndrome (MetS). We performed a genome-wide association study (GWAS) of MetS in 10,093 Indian individuals (6,617 MetS and 3,476 controls) of Indo-European origin, that belong to our previous biorepository of The Indian Diabetes Consortium (INDICO). The study was conducted in two stages-discovery phase (N = 2,158) and replication phase (N = 7,935). We discovered two variants within/near the CETP gene-rs1800775 and rs3816117-associated with MetS at genome-wide significance level during replication phase in Indians. Additional CETP loci rs7205804, rs1532624, rs3764261, rs247617, and rs173539 also cropped up as modest signals in Indians. Haplotype association analysis revealed GCCCAGC as the strongest haplotype within the CETP locus constituting all seven CETP signals. In combined analysis, we perceived a novel and functionally relevant sub-GWAS significant locus-rs16890462 in the vicinity of SFRP1 gene. Overlaying gene regulatory data from ENCODE database revealed that single nucleotide polymorphism (SNP) rs16890462 resides in repressive chromatin in human subcutaneous adipose tissue as characterized by the enrichment of H3K27me3 and CTCF marks (repressive gene marks) and diminished H3K36me3 marks (activation gene marks). The variant displayed active DNA methylation marks in adipose tissue, suggesting its likely regulatory activity. Further, the variant also disrupts a potential binding site of a key transcription factor, NRF2, which is known for involvement in obesity and metabolic syndrome.

Hidden Markov models lead to higher resolution maps of mutation signature activity in cancer

Knowing the activity of the mutational processes shaping a cancer genome may provide insight into tumorigenesis and personalized therapy. It is thus important to characterize the signatures of active mutational processes in patients from their patterns of single base substitutions. However, mutational processes do not act uniformly on the genome, leading to statistical dependencies among neighboring mutations. To account for such dependencies, we develop the first sequence-dependent model, SigMa, for mutation signatures. We apply SigMa to characterize genomic and other factors that influence the activity of mutation signatures in breast cancer. We show that SigMa outperforms previous approaches, revealing novel insights on signature etiology. The source code for SigMa is publicly available at https://github.com/lrgr/sigma.

Connectivity Mapping: Methods and Applications

Connectivity mapping resources consist of signatures representing changes in cellular state following systematic small-molecule, disease, gene, or other form of perturbations. Such resources enable the characterization of signatures from novel perturbations based on similarity; provide a global view of the space of many themed perturbations; and allow the ability to predict cellular, tissue, and organismal phenotypes for perturbagens. A signature search engine enables hypothesis generation by finding connections between query signatures and the database of signatures. This framework has been used to identify connections between small molecules and their targets, to discover cell-specific responses to perturbations and ways to reverse disease expression states with small molecules, and to predict small-molecule mimickers for existing drugs. This review provides a historical perspective and the current state of connectivity mapping resources with a focus on both methodology and community implementations.

Best practice for CRISPR design using current tools and resources

Users facing the task of designing gRNAs for a CRISPR-mutagenesis experiment are typically confronted with a large variety of possible tools and existing libraries. Here we examine the design principles for such resources, and suggest a best practice which allows a user to evaluate and effectively use any of the existing CRISPR design tools or genome-wide libraries.

ChEA3: transcription factor enrichment analysis by orthogonal omics integration

Identifying the transcription factors (TFs) responsible for observed changes in gene expression is an important step in understanding gene regulatory networks. ChIP-X Enrichment Analysis 3 (ChEA3) is a transcription factor enrichment analysis tool that ranks TFs associated with user-submitted gene sets. The ChEA3 background database contains a collection of gene set libraries generated from multiple sources including TF-gene co-expression from RNA-seq studies, TF-target associations from ChIP-seq experiments, and TF-gene co-occurrence computed from crowd-submitted gene lists. Enrichment results from these distinct sources are integrated to generate a composite rank that improves the prediction of the correct upstream TF compared to ranks produced by individual libraries. We compare ChEA3 with existing TF prediction tools and show that ChEA3 performs better. By integrating the ChEA3 libraries, we illuminate general transcription factor properties such as whether the TF behaves as an activator or a repressor. The ChEA3 web-server is available from https://amp.pharm.mssm.edu/ChEA3.

EpiAlignment: alignment with both DNA sequence and epigenomic data

Comparative epigenomics, which subjects both epigenome and genome to interspecies comparison, has become a powerful approach to reveal regulatory features of the genome. Thus elucidated regulatory features surpass the information derived from comparison of genomic sequences alone. Here, we present EpiAlignment, a web-based tool to align genomic regions with both DNA sequence and epigenomic data. EpiAlignment takes DNA sequence and epigenomic profiles derived by ChIP-seq from two species as input data, and outputs the best semi-global alignments. These alignments are based on EpiAlignment scores, computed by a dynamic programming algorithm that accounts for both sequence alignment and epigenome similarity. For timely response, the EpiAlignment web server automatically initiates up to 140 computing threads depending on the size of user input data. For users’ convenience, we have pre-compiled the comparable human and mouse epigenome datasets in matched cell types and tissues from the Roadmap Epigenomics and ENCODE consortia. Users can either upload their own data or select pre-compiled datasets as inputs for EpiAlignment analyses. Results are presented in graphical and tabular formats where the entries can be interactively expanded to visualize additional features of these aligned regions. EpiAlignment is available at https://epialign.ucsd.edu/.

Reversible Regulation of Polyubiquitin Gene UBC via Modified Inducible CRISPR/Cas9 System

Ubiquitin is required under both normal and stress conditions. Under stress conditions, upregulation of the polyubiquitin gene UBC is essential to meet the requirement of increased ubiquitin levels to confer stress resistance. However, UBC upregulation is usually observed only under stress conditions and not under normal conditions. Therefore, it has not been possible to upregulate UBC under normal conditions to study the effect of excess ubiquitin on cellular machinery. Recently, the CRISPR/Cas9 system has been widely used in biological research as a useful tool to study gene disruption effects. In this study, using an inducible CRISPR/Cas9 variant, a dCas9-VP64 fusion protein, combined with a single guide RNA (sgRNA) containing MS2 aptamer loops and MS2-p65-HSF1, we developed a system to increase the ubiquitin pool via upregulation of UBC. Although it is challenging to upregulate the expression of a gene that is already expressed at high levels, the significance of our system is that UBC upregulation can be induced in an efficient, reversible manner that is compatible with cellular processes, even under normal conditions. This system can be used to study ubiquitin pool dynamics and it will be a useful tool in identifying the role of ubiquitin under normal and stress conditions.

Hit and Run Transcriptional Repressors Are Difficult to Catch in the Act

Transcriptional silencing may not necessarily depend on the continuous residence of a sequence-specific repressor at a control element and may act via a "hit and run" mechanism. Due to limitations in assays that detect transcription factor (TF) binding, such as chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq), this phenomenon may be challenging to detect and therefore its prevalence may be underappreciated. To explore this possibility, erythroid gene promoters that are regulated directly by GATA1 in an inducible system are analyzed. It is found that many regulated genes are bound immediately after induction of GATA1 but the residency of GATA1 decreases over time, particularly at repressed genes. Furthermore, it is shown that the repressive mark H3K27me3 is seldom associated with bound repressors, whereas, in contrast, the active (H3K4me3) histone mark is overwhelmingly associated with TF binding. It is hypothesized that during cellular differentiation and development, certain genes are silenced by repressive TFs that subsequently vacate the region. Catching such repressor TFs in the act of silencing via assays such as ChIP-seq is thus a temporally challenging prospect. The use of inducible systems, epitope tags, and alternative techniques may provide opportunities for detecting elusive "hit and run" transcriptional silencing. Also see the video abstract here https://youtu.be/vgrsoP_HF3g.

Integrating proteomics and transcriptomics for the identification of potential targets in early colorectal cancer

Colorectal cancer (CRC) is one of the most common malignancies worldwide. At present, CRC can often be treated upon diagnosis at stage I or II, or when dysplasia is detected; however, 60-70% of cases are not diagnosed until they have developed into late stages of the disease or until the malignancy is identified. Diagnosis of CRC at an early stage remains a challenge due to the absence of early-stage-specific biomarkers. To identify potential targets of early stage CRC, label-free proteomics analysis was applied to paired tumor-benign tissue samples from patients with stage II CRC (n=21). A total of 2,968 proteins were identified; corresponding RNA-Sequencing data were retrieved from The Cancer Genome Atlas-colon adenocarcinoma. Numerous bioinformatics methods, including differential expression analysis, weighted correlation network analysis, Gene Ontology and protein-protein interaction analyses, were applied to the proteomics and transcriptomics data. A total of 111 key proteins, which appeared as both differentially expressed proteins and mRNAs in the hub module, were identified as key candidates. Among these, three potential targets [protein-arginine deiminase type-2 (PADI2), Fc fragment of IgG binding protein (FCGBP) and phosphoserine aminotransferase 1] were identified from the pathological data. Furthermore, the survival analysis indicated that PADI2 and FCGBP were associated with the prognosis of CRC. The findings of the present study suggested potential targets for the identification of early stage CRC, and may improve understanding of the mechanism underlying the occurrence of CRC.

METTL1 Promotes let-7 MicroRNA Processing via m7G Methylation

7-methylguanosine (m7G) is present at mRNA caps and at defined internal positions within tRNAs and rRNAs. However, its detection within low-abundance mRNAs and microRNAs (miRNAs) has been hampered by a lack of sensitive detection strategies. Here, we adapt a chemical reactivity assay to detect internal m7G in miRNAs. Using this technique (Borohydride Reduction sequencing [BoRed-seq]) alongside RNA immunoprecipitation, we identify m7G within a subset of miRNAs that inhibit cell migration. We show that the METTL1 methyltransferase mediates m7G methylation within miRNAs and that this enzyme regulates cell migration via its catalytic activity. Using refined mass spectrometry methods, we map m7G to a single guanosine within the let-7e-5p miRNA. We show that METTL1-mediated methylation augments let-7 miRNA processing by disrupting an inhibitory secondary structure within the primary miRNA transcript (pri-miRNA). These results identify METTL1-dependent N7-methylation of guanosine as a new RNA modification pathway that regulates miRNA structure, biogenesis, and cell migration.

Preformed chromatin topology assists transcriptional robustness of Shh during limb development

Long-range gene regulation involves physical proximity between enhancers and promoters to generate precise patterns of gene expression in space and time. However, in some cases, proximity coincides with gene activation, whereas, in others, preformed topologies already exist before activation. In this study, we investigate the preformed configuration underlying the regulation of the Shh gene by its unique limb enhancer, the ZRS, in vivo during mouse development. Abrogating the constitutive transcription covering the ZRS region led to a shift within the Shh-ZRS contacts and a moderate reduction in Shh transcription. Deletion of the CTCF binding sites around the ZRS resulted in the loss of the Shh-ZRS preformed interaction and a 50% decrease in Shh expression but no phenotype, suggesting an additional, CTCF-independent mechanism of promoter-enhancer communication. This residual activity, however, was diminished by combining the loss of CTCF binding with a hypomorphic ZRS allele, resulting in severe Shh loss of function and digit agenesis. Our results indicate that the preformed chromatin structure of the Shh locus is sustained by multiple components and acts to reinforce enhancer-promoter communication for robust transcription.

Can Epigenetics of Endothelial Dysfunction Represent the Key to Precision Medicine in Type 2 Diabetes Mellitus?

In both developing and industrialized Countries, the growing prevalence of Type 2 Diabetes Mellitus (T2DM) and the severity of its related complications make T2DM one of the most challenging metabolic diseases worldwide. The close relationship between genetic and environmental factors suggests that eating habits and unhealthy lifestyles may significantly affect metabolic pathways, resulting in dynamic modifications of chromatin-associated proteins and homeostatic transcriptional responses involved in the progression of T2DM. Epigenetic mechanisms may be implicated in the complex processes linking environmental factors to genetic predisposition to metabolic disturbances, leading to obesity and type 2 diabetes mellitus (T2DM). Endothelial dysfunction represents an earlier marker and an important player in the development of this disease. Dysregulation of the endothelial ability to produce and release vasoactive mediators is recognized as the initial feature of impaired vascular activity under obesity and other insulin resistance conditions and undoubtedly concurs to the accelerated progression of atherosclerotic lesions and overall cardiovascular risk in T2DM patients. This review aims to summarize the most current knowledge regarding the involvement of epigenetic changes associated with endothelial dysfunction in T2DM, in order to identify potential targets that might contribute to pursuing “precision medicine” in the context of diabetic illness.

Growth disrupting mutations in epigenetic regulatory molecules are associated with abnormalities of epigenetic aging

Germline mutations in fundamental epigenetic regulatory molecules including DNA methyltransferase 3 alpha (DNMT3A) are commonly associated with growth disorders, whereas somatic mutations are often associated with malignancy. We profiled genome-wide DNA methylation patterns in DNMT3A c.2312G > A; p.(Arg771Gln) carriers in a large Amish sibship with Tatton-Brown–Rahman syndrome (TBRS), their mosaic father, and 15 TBRS patients with distinct pathogenic de novo DNMT3A variants. This defined widespread DNA hypomethylation at specific genomic sites enriched at locations annotated as genes involved in morphogenesis, development, differentiation, and malignancy predisposition pathways. TBRS patients also displayed highly accelerated DNA methylation aging. These findings were most marked in a carrier of the AML-associated driver mutation p.Arg882Cys. Our studies additionally defined phenotype-related accelerated and decelerated epigenetic aging in two histone methyltransferase disorders: NSD1 Sotos syndrome overgrowth disorder and KMT2D Kabuki syndrome growth impairment. Together, our findings provide fundamental new insights into aberrant epigenetic mechanisms, the role of epigenetic machinery maintenance, and determinants of biological aging in these growth disorders.

Bridging the gap between reference and real transcriptomes

Genetic, transcriptional, and post-transcriptional variations shape the transcriptome of individual cells, rendering establishing an exhaustive set of reference RNAs a complicated matter. Current reference transcriptomes, which are based on carefully curated transcripts, are lagging behind the extensive RNA variation revealed by massively parallel sequencing. Much may be missed by ignoring this unreferenced RNA diversity. There is plentiful evidence for non-reference transcripts with important phenotypic effects. Although reference transcriptomes are inestimable for gene expression analysis, they may turn limiting in important medical applications. We discuss computational strategies for retrieving hidden transcript diversity.

Towards a deeper annotation of human lncRNAs

A substantial fraction of the human transcriptome is composed of the so-called long noncoding RNAs (lncRNAs), yet the available catalogs of known lncRNAs are far from complete. Moreover, functional studies of these RNAs are challenged by several factors, such as their tissue-specific expression and functional heterogeneity, resulting in only ca. 1% of them being well characterized. Here, we describe a set of 41,400 novel lncRNAs discovered with RNA-Seq data from 1463 samples encompassing diverse tissues and cell lines. We utilized publicly available transcriptomic and genomic data to provide their characteristics, such as tissue specificity, cellular abundance, polyA status, cellular localization, evolutionary conservation and transcript stability, which allowed us to speculate on their possible biological roles. We also pinpointed 24 novel lncRNAs as candidates for breast cancer biomarkers. The results bring us closer to a comprehensive annotation of human lncRNAs, though vast amounts of further work are needed to validate the predictions and fully decipher their biology. This article is part of a Special Issue entitled: ncRNA in control of gene expression edited by Kotb Abdelmohsen.

Securing the future of research computing in the biosciences

Improvements in technology often drive scientific discovery. Therefore, research requires sustained investment in the latest equipment and training for the researchers who are going to use it. Prioritising and administering infrastructure investment is challenging because future needs are difficult to predict. In the past, highly computationally demanding research was associated primarily with particle physics and astronomy experiments. However, as biology becomes more quantitative and bioscientists generate more and more data, their computational requirements may ultimately exceed those of physical scientists. Computation has always been central to bioinformatics, but now imaging experiments have rapidly growing data processing and storage requirements. There is also an urgent need for new modelling and simulation tools to provide insight and understanding of these biophysical experiments. Bioscience communities must work together to provide the software and skills training needed in their areas. Research-active institutions need to recognise that computation is now vital in many more areas of discovery and create an environment where it can be embraced. The public must also become aware of both the power and limitations of computing, particularly with respect to their health and personal data.

An in-silico approach for discovery of microRNA-TF regulation of DISC1 interactome mediating neuronal migration

Neuronal migration constitutes an important step in corticogenesis; dysregulation of the molecular mechanisms mediating this crucial step in neurodevelopment may result in various neuropsychiatric disorders. By curating experimental data from published literature, we identified eight functional modules involving Disrupted-in-schizophrenia 1 (DISC1) and its interacting proteins that regulate neuronal migration. We then identified miRNAs and transcription factors (TFs) that form functional feedback loops and regulate gene expression of the DISC1 interactome. Using this curated data, we conducted in-silico modeling of the DISC1 interactome involved in neuronal migration and identified the proteins that either facilitate or inhibit neuronal migrational processes. We also studied the effect of perturbation of miRNAs and TFs in feedback loops on the DISC1 interactome. From these analyses, we discovered that STAT3, TCF3, and TAL1 (through feedback loop with miRNAs) play a critical role in the transcriptional control of DISC1 interactome thereby regulating neuronal migration. To the best of our knowledge, regulation of the DISC1 interactome mediating neuronal migration by these TFs has not been previously reported. These potentially important TFs can serve as targets for undertaking validation studies, which in turn can reveal the molecular processes that cause neuronal migration defects underlying neurodevelopmental disorders. This underscores the importance of the use of in-silico techniques in aiding the discovery of mechanistic evidence governing important molecular and cellular processes. The present work is one such step towards the discovery of regulatory factors of the DISC1 interactome that mediates neuronal migration.

Retrotransposon elements among initial sites of hepatitis B virus integration into human genome in the HepG2-NTCP cell infection model

Integration of hepatitis B virus (HBV) DNA into host's genome is evident in all stages and models of HBV infection. Investigations of the initial virus-host junctions have been just recently initiated since their nature may promote liver oncogenesis immediately following infection. We examined the time-frame and host sites at which HBV integrates in HepG2 cells overexpressing sodium taurocholate co-transporting polypeptide (NTCP) receptor mediating HBV entry. HepG2-NTCP cells were analyzed from 15 min to 13 days post-infection (p.i.). The results showed that except for 15 min p.i., HBV-host integrations were detected at all time points thereafter. At 30 min p.i., virus junctions with retrotransposon SINE and with neuroblastoma breakpoint family member 1 gene were detected. At one-hour p.i., HBV integration with retrotransposon THE-1B-LTR was identified, while virus insertions into proline-rich protein and protein kinase cGMP-dependent type 1 encoding genes were found at 3 h p.i. Fusion with runt-related transcription factor 1 was detected at 24 h p.i. and merges with 9 different genes at 13 day p.i. The data showed that retrotransposon elements are frequent among first-hit sites of HBV insertion. This may suggest a mechanism by which HBV DNA may spread across host's genome from earliest stages of infection.

Atlas of RNA sequencing profiles for normal human tissues

Comprehensive analysis of molecular pathology requires a collection of reference samples representing normal tissues from healthy donors. For the available limited collections of normal tissues from postmortal donors, there is a problem of data incompatibility, as different datasets generated using different experimental platforms often cannot be merged in a single panel. Here, we constructed and deposited the gene expression database of normal human tissues based on uniformly screened original sequencing data. In total, 142 solid tissue samples representing 20 organs were taken from post-mortal human healthy donors of different age killed in road accidents no later than 36 hours after death. Blood samples were taken from 17 healthy volunteers. We then compared them with the 758 transcriptomic profiles taken from the other databases. We found that overall 463 biosamples showed tissue-specific rather than platform- or database-specific clustering and could be aggregated in a single database termed Oncobox Atlas of Normal Tissue Expression (ANTE). Our data will be useful to all those working with the analysis of human gene expression.

ETV6-RUNX1 interacts with a region in SPIB intron 1 to regulate gene expression in pre-B-cell acute lymphoblastic leukemia

The most frequently occurring genetic abnormality in pediatric B-lymphocyte-lineage acute lymphoblastic leukemia is the t(12;21) chromosomal translocation that results in a ETV6-RUNX1 (also known as TEL-AML1) fusion gene. Expression of ETV6-RUNX1 induces a preleukemic condition leading to acquisition of secondary driver mutations, but the mechanism is poorly understood. SPI-B (encoded by SPIB) is an important transcriptional activator of B-cell development and differentiation. We hypothesized that SPIB is directly transcriptionally repressed by ETV6-RUNX1. Using chromatin immunoprecipitation, we identified a regulatory region in the first intron of SPIB that interacts with ETV6-RUNX1. Mutation of the RUNX1 binding site in SPIB intron 1 prevented transcriptional repression in transient transfection assays. Next, we sought to determine to what extent gene expression in REH cells can be altered by ectopic SPI-B expression. SPI-B expression was forced using CRISPR-mediated gene activation and also using a retroviral vector. Forced expression of SPI-B resulted in altered gene expression and, at high levels, impaired cell proliferation and induced apoptosis. Finally, we identified CARD11 and CDKN1A (encoding p21) as transcriptional targets of SPI-B involved in regulation of proliferation and apoptosis. Taken together, this study identifies SPIB as an important target of ETV6-RUNX1 in regulation of B-cell gene expression in t(12;21) leukemia.

The Promoter Regions of Intellectual Disability-Associated Genes Are Uniquely Enriched in LTR Sequences of the MER41 Primate-Specific Endogenous Retrovirus: An Evolutionary Connection Between Immunity and Cognition

Social behavior and neuronal connectivity in rodents have been shown to be shaped by the prototypical T lymphocyte-derived pro-inflammatory cytokine Interferon-gamma (IFNγ). It has also been demonstrated that STAT1 (Signal Transducer And Activator Of Transcription 1), a transcription factor (TF) crucially involved in the IFNγ pathway, binds consensus sequences that, in humans, are located with a high frequency in the LTRs (Long Terminal Repeats) of the MER41 family of primate-specific HERVs (Human Endogenous Retroviruses). However, the putative role of an IFNγ/STAT1/MER41 pathway in human cognition and/or behavior is still poorly documented. Here, we present evidence that the promoter regions of intellectual disability-associated genes are uniquely enriched in LTR sequences of the MER41 HERVs. This observation is specific to MER41 among more than 130 HERVs examined. Moreover, we have not found such a significant enrichment in the promoter regions of genes that associate with autism spectrum disorder (ASD) or schizophrenia. Interestingly, ID-associated genes exhibit promoter-localized MER41 LTRs that harbor TF binding sites (TFBSs) for not only STAT1 but also other immune TFs such as, in particular, NFKB1 (Nuclear Factor Kappa B Subunit 1) and STAT3 (Signal Transducer And Activator Of Transcription 3). Moreover, IL-6 (Interleukin 6) rather than IFNγ, is identified as the main candidate cytokine regulating such an immune/MER41/cognition pathway. Of note, differences between humans and chimpanzees are observed regarding the insertion sites of MER41 LTRs in the promoter regions of ID-associated genes. Finally, a survey of the human proteome has allowed us to map a protein-protein network which links the identified immune/MER41/cognition pathway to FOXP2 (Forkhead Box P2), a key TF involved in the emergence of human speech. Our work suggests that together with the evolution of immune genes, the stepped self-domestication of MER41 in the genomes of primates could have contributed to cognitive evolution. We further propose that non-inherited forms of ID might result from the untimely or quantitatively inappropriate expression of immune signals, notably IL-6, that putatively regulate cognition-associated genes via promoter-localized MER41 LTRs.

Genome-wide association study of blood lipids in Indians confirms universality of established variants

Lipids foster energy production and their altered levels have been coupled with metabolic ailments. Indians feature high prevalence of metabolic diseases, yet uncharacterized for genes regulating lipid homeostasis. We performed first GWAS for quantitative lipids (total cholesterol, LDL, HDL, and triglycerides) exclusively in 5271 Indians. Further to corroborate our genetic findings, we investigated DNA methylation marks in peripheral blood in Indians at the identified loci (N = 233) and retrieved gene regulatory features from public domains. Recurrent GWAS loci-CELSR2, CETP, LPL, ZNF259, and BUD13 cropped up as lead signals in Indians, reflecting their universal applicability. Besides established variants, we found certain unreported variants at sub-genome-wide level-QKI, REEP3, TMCC2, FAM129C, FAM241B, and LOC100506207. These variants though failed to attain GWAS significance in Indians, but largely turned out to be active CpG sites in human subcutaneous adipose tissue and showed robust association to two or more lipid traits. Of which, QKI variants showed significant association to all four lipid traits and their designated region was observed to be a key gene regulatory segment denoting active transcription particularly in human subcutaneous adipose tissue. Both established and novel loci were observed to be significantly associated with altered DNA methylation in Indians for specific CpGs that resided in key regulatory elements. Further, gene-based association analysis pinpointed novel GWAS loci-LINC01340 and IQCJ-SCHIP1 for TC; IFT27, IFT88, and LINC02141 for HDL; and TEX26 for TG. Present study ascertains universality of selected known genes and also identifies certain novel loci for lipids in Indians by integrating data from various levels of gene regulation.

Computational analysis of the evolutionarily conserved Missing In Metastasis/Metastasis Suppressor 1 gene predicts novel interactions, regulatory regions and transcriptional control

Missing in Metastasis (MIM), or Metastasis Suppressor 1 (MTSS1), is a highly conserved protein, which links the plasma membrane to the actin cytoskeleton. MIM has been implicated in various cancers, however, its modes of action remain largely enigmatic. Here, we performed an extensive in silico characterisation of MIM to gain better understanding of its function. We detected previously unappreciated functional motifs including adaptor protein (AP) complex interaction site and a C-helix, pointing to a role in endocytosis and regulation of actin dynamics, respectively. We also identified new functional regions, characterised with phosphorylation sites or distinct hydrophilic properties. Strong negative selection during evolution, yielding high conservation of MIM, has been combined with positive selection at key sites. Interestingly, our analysis of intra-molecular co-evolution revealed potential regulatory hotspots that coincided with reduced potentially pathogenic polymorphisms. We explored databases for the mutations and expression levels of MIM in cancer. Experimentally, we focused on chronic lymphocytic leukaemia (CLL), where MIM showed high overall expression, however, downregulation on poor prognosis samples. Finally, we propose strong conservation of MTSS1 also on the transcriptional level and predict novel transcriptional regulators. Our data highlight important targets for future studies on the role of MIM in different tissues and cancers.

Making Sense of the Epigenome Using Data Integration Approaches

Epigenetic research involves examining the mitotically heritable processes that regulate gene expression, independent of changes in the DNA sequence. Recent technical advances such as whole-genome bisulfite sequencing and affordable epigenomic array-based technologies, allow researchers to measure epigenetic profiles of large cohorts at a genome-wide level, generating comprehensive high-dimensional datasets that may contain important information for disease development and treatment opportunities. The epigenomic profile for a certain disease is often a result of the complex interplay between multiple genetic and environmental factors, which poses an enormous challenge to visualize and interpret these data. Furthermore, due to the dynamic nature of the epigenome, it is critical to determine causal relationships from the many correlated associations. In this review we provide an overview of recent data analysis approaches to integrate various omics layers to understand epigenetic mechanisms of complex diseases, such as obesity and cancer. We discuss the following topics: (i) advantages and limitations of major epigenetic profiling techniques, (ii) resources for standardization, annotation and harmonization of epigenetic data, and (iii) statistical methods and machine learning methods for establishing data-driven hypotheses of key regulatory mechanisms. Finally, we discuss the future directions for data integration that shall facilitate the discovery of epigenetic-based biomarkers and therapies.

CDK4/6 inhibitors target SMARCA4-determined cyclin D1 deficiency in hypercalcemic small cell carcinoma of the ovary

Inactivating mutations in SMARCA4 (BRG1), a key SWI/SNF chromatin remodelling gene, underlie small cell carcinoma of the ovary, hypercalcemic type (SCCOHT). To reveal its druggable vulnerabilities, we perform kinase-focused RNAi screens and uncover that SMARCA4-deficient SCCOHT cells are highly sensitive to the inhibition of cyclin-dependent kinase 4/6 (CDK4/6). SMARCA4 loss causes profound downregulation of cyclin D1, which limits CDK4/6 kinase activity in SCCOHT cells and leads to in vitro and in vivo susceptibility to CDK4/6 inhibitors. SCCOHT patient tumors are deficient in cyclin D1 yet retain the retinoblastoma-proficient/p16INK4a-deficient profile associated with positive responses to CDK4/6 inhibitors. Thus, our findings indicate that CDK4/6 inhibitors, approved for a breast cancer subtype addicted to CDK4/6 activation, could be repurposed to treat SCCOHT. Moreover, our study suggests a novel paradigm whereby critically low oncogene levels, caused by loss of a driver tumor suppressor, may also be exploited therapeutically.

Evolution of Immune Systems From Viruses and Transposable Elements

Virus-derived sequences and transposable elements constitute a substantial portion of many cellular genomes. Recent insights reveal the intimate evolutionary relationship between these sequences and various cellular immune pathways. At the most basic level, superinfection exclusion may be considered a prototypical virus-mediated immune system that has been described in both prokaryotes and eukaryotes. More complex immune mechanisms fully or partially derived from mobile genetic elements include CRISPR-Cas of prokaryotes and the RAG1/2 system of vertebrates, which provide immunological memory of foreign genetic elements and generate antibody and T cell receptor diversity, respectively. In this review, we summarize the current knowledge on the contribution of mobile genetic elements to the evolution of cellular immune pathways. A picture is emerging in which the various cellular immune systems originate from and are spread by viruses and transposable elements. Immune systems likely evolved from simple superinfection exclusion to highly complex defense strategies.

Systems Biology Approach to Identify Novel Genomic Determinants for Pancreatic Cancer Pathogenesis

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy with a 5-year survival rate of <8%. Its dismal prognosis stems from inefficient therapeutic modalities owing to the lack of understanding about pancreatic cancer pathogenesis. Considering the molecular complexity and heterogeneity of PDAC, identification of novel molecular contributors involved in PDAC onset and progression using global "omics" analysis will pave the way to improved strategies for disease prevention and therapeutic targeting. Meta-analysis of multiple miRNA microarray datasets containing healthy controls (HC), chronic pancreatitis (CP) and PDAC cases, identified 13 miRNAs involved in the progression of PDAC. These miRNAs showed dysregulation in both tissue as well as blood samples, along with progressive decrease in expression from HC to CP to PDAC. Gene-miRNA interaction analysis further elucidated 5 miRNAs (29a/b, 27a, 130b and 148a) that are significantly downregulated in conjunction with concomitant upregulation of their target genes throughout PDAC progression. Among these, miRNA-29a/b targeted genes were found to be most significantly altered in comparative profiling of HC, CP and PDAC, indicating its involvement in malignant evolution. Further, pathway analysis suggested direct involvement of miRNA-29a/b in downregulating the key pathways associated with PDAC development and metastasis including focal adhesion signaling and extracellular matrix organization. Our systems biology data analysis, in combination with real-time PCR validation indicates direct functional involvement of miRNA-29a in PDAC progression and is a potential prognostic marker and therapeutic candidate for patients with progressive disease.

Cistrome Data Browser: expanded datasets and new tools for gene regulatory analysis

The Cistrome Data Browser (DB) is a resource of human and mouse cis-regulatory information derived from ChIP-seq, DNase-seq and ATAC-seq chromatin profiling assays, which map the genome-wide locations of transcription factor binding sites, histone post-translational modifications and regions of chromatin accessible to endonuclease activity. Currently, the Cistrome DB contains approximately 47,000 human and mouse samples with about 24,000 newly collected datasets compared to the previous release two years ago. Furthermore, the Cistrome DB has a new Toolkit module with several features that allow users to better utilize the large-scale ChIP-seq, DNase-seq, and ATAC-seq data. First, users can query the factors which are likely to regulate a specific gene of interest. Second, the Cistrome DB Toolkit facilitates searches for factor binding, histone modifications, and chromatin accessibility in any given genomic interval shorter than 2Mb. Third, the Toolkit can determine the most similar ChIP-seq, DNase-seq, and ATAC-seq samples in terms of genomic interval overlaps with user-provided genomic interval sets. The Cistrome DB is a user-friendly, up-to-date, and well maintained resource, and the new tools will greatly benefit the biomedical research community. The database is freely available at http://cistrome.org/db, and the Toolkit is at http://dbtoolkit.cistrome.org.

Long noncoding RNA H19 is a critical oncogenic driver and contributes to epithelial-mesenchymal transition in papillary thyroid carcinoma

BACKGROUND

Growing evidence has indicated that the long noncoding RNA H19 (lncRNA H19), frequently deregulated in almost all tumor types tested, acted as a pivotal contributor to both cancer initiation and progression. However, the role of lncRNA H19 in human papillary thyroid carcinoma (PTC) remains controversial. The aim of the study was to investigate the expression and potential function of lncRNA H19 in human PTC.

PATIENTS AND METHODS

The lncRNA H19 level was determined by quantitative real-time (RT)-PCR analyses in 58 PTC tissue samples and their paired paracancerous tissue samples. RNA interference, RT-PCR analysis, and Western blot assay were used to determine the impact of lncRNA H19 on epithelial-mesenchymal transition (EMT) markers in human PTC cells. The migratory and invasive capacities of PTC cells were determined by wound-healing and transwell migration and invasion assays.

RESULTS

lncRNA H19 expression was 2.417-fold higher in PTC tissues than their paired paracancerous tissue (95% CI: 1.898-2.935, P<0.0001). Higher level of lncRNA H19 was correlated to elevated expression of Vimentin, ZEB2, Twist, and Snail2. Inhibition of lncRNA H19 resulted in upregulation of E-cadherin and downregulation of Vimentin both at mRNA and protein levels. Conversely, enforced expression of the exogenous lncRNA H19 led to E-cadherin mRNA and protein downregulation and relative upregulation of Vimentin. Moreover, wound-healing and transwell migration and invasion assays showed that lncRNA H19 could promote the migratory and invasive abilities of PTC cells.

CONCLUSION

The level of lncRNA H19 was significantly higher in PTC tissues than paired paracancerous tissue or normal tissues. Overexpression of lncRNA H19 was correlated with higher tumor burden of PTC. It also contributes to EMT process, as well as promotes migration and invasion of PTC cells.

Hi-C analysis: from data generation to integration

In the epigenetics field, large-scale functional genomics datasets of ever-increasing size and complexity have been produced using experimental techniques based on high-throughput sequencing. In particular, the study of the 3D organization of chromatin has raised increasing interest, thanks to the development of advanced experimental techniques. In this context, Hi-C has been widely adopted as a high-throughput method to measure pairwise contacts between virtually any pair of genomic loci, thus yielding unprecedented challenges for analyzing and handling the resulting complex datasets. In this review, we focus on the increasing complexity of available Hi-C datasets, which parallels the adoption of novel protocol variants. We also review the complexity of the multiple data analysis steps required to preprocess Hi-C sequencing reads and extract biologically meaningful information. Finally, we discuss solutions for handling and visualizing such large genomics datasets.

CRISPR-mediated deletion of prostate cancer risk-associated CTCF loop anchors identifies repressive chromatin loops

BACKGROUND

Recent genome-wide association studies (GWAS) have identified more than 100 loci associated with increased risk of prostate cancer, most of which are in non-coding regions of the genome. Understanding the function of these non-coding risk loci is critical to elucidate the genetic susceptibility to prostate cancer.

RESULTS

We generate genome-wide regulatory element maps and performed genome-wide chromosome confirmation capture assays (in situ Hi-C) in normal and tumorigenic prostate cells. Using this information, we annotate the regulatory potential of 2,181 fine-mapped prostate cancer risk-associated SNPs and predict a set of target genes that are regulated by prostate cancer risk-related H3K27Ac-mediated loops. We next identify prostate cancer risk-associated CTCF sites involved in long-range chromatin loops. We use CRISPR-mediated deletion to remove prostate cancer risk-associated CTCF anchor regions and the CTCF anchor regions looped to the prostate cancer risk-associated CTCF sites, and we observe up to 100-fold increases in expression of genes within the loops when the prostate cancer risk-associated CTCF anchor regions are deleted.

CONCLUSIONS

We identify GWAS risk loci involved in long-range loops that function to repress gene expression within chromatin loops. Our studies provide new insights into the genetic susceptibility to prostate cancer.

Developmental Chromatin Restriction of Pro-Growth Gene Networks Acts as an Epigenetic Barrier to Axon Regeneration in Cortical Neurons

Axon regeneration in the central nervous system is prevented in part by a developmental decline in the intrinsic regenerative ability of maturing neurons. This loss of axon growth ability likely reflects widespread changes in gene expression, but the mechanisms that drive this shift remain unclear. Chromatin accessibility has emerged as a key regulatory mechanism in other cellular contexts, raising the possibility that chromatin structure may contribute to the age-dependent loss of regenerative potential. Here we establish an integrated bioinformatic pipeline that combines analysis of developmentally dynamic gene networks with transcription factor regulation and genome-wide maps of chromatin accessibility. When applied to the developing cortex, this pipeline detected overall closure of chromatin in sub-networks of genes associated with axon growth. We next analyzed mature CNS neurons that were supplied with various pro-regenerative transcription factors. Unlike prior results with SOX11 and KLF7, here we found that neither JUN nor an activated form of STAT3 promoted substantial corticospinal tract regeneration. Correspondingly, chromatin accessibility in JUN or STAT3 target genes was substantially lower than in predicted targets of SOX11 and KLF7. Finally, we used the pipeline to predict pioneer factors that could potentially relieve chromatin constraints at growth-associated loci. Overall this integrated analysis substantiates the hypothesis that dynamic chromatin accessibility contributes to the developmental decline in axon growth ability and influences the efficacy of pro-regenerative interventions in the adult, while also pointing toward selected pioneer factors as high-priority candidates for future combinatorial experiments. © 2018 Wiley Periodicals, Inc. Develop Neurobiol 00: 000-000, 2018.

KLF6 and STAT3 co-occupy regulatory DNA and functionally synergize to promote axon growth in CNS neurons

The failure of axon regeneration in the CNS limits recovery from damage and disease. Members of the KLF family of transcription factors can exert both positive and negative effects on axon regeneration, but the underlying mechanisms are unclear. Here we show that forced expression of KLF6 promotes axon regeneration by corticospinal tract neurons in the injured spinal cord. RNA sequencing identified 454 genes whose expression changed upon forced KLF6 expression in vitro, including sub-networks that were highly enriched for functions relevant to axon extension including cytoskeleton remodeling, lipid synthesis, and bioenergetics. In addition, promoter analysis predicted a functional interaction between KLF6 and a second transcription factor, STAT3, and genome-wide footprinting using ATAC-Seq data confirmed frequent co-occupancy. Co-expression of the two factors yielded a synergistic elevation of neurite growth in vitro. These data clarify the transcriptional control of axon growth and point the way toward novel interventions to promote CNS regeneration.

Patterns and mechanisms of structural variations in human cancer

Next-generation sequencing technology has enabled the comprehensive detection of genomic alterations in human somatic cells, including point mutations, chromosomal rearrangements, and structural variations (SVs). Using sophisticated bioinformatics algorithms, unbiased catalogs of SVs are emerging from thousands of human cancer genomes for the first time. Via careful examination of SV breakpoints at single-nucleotide resolution as well as local DNA copy number changes, diverse patterns of genomic rearrangements are being revealed. These "SV signatures" provide deep insight into the mutational processes that have shaped genome changes in human somatic cells. This review summarizes the characteristics of recently identified complex SVs, including chromothripsis, chromoplexy, microhomology-mediated breakage-induced replication (MMBIR), and others, to provide a holistic snapshot of the current knowledge on genomic rearrangements in somatic cells.

Non-coding genetic variation shaping mental health

•

Gene expression determined by the genome mediating a response to cell environment.

•

Genetic variation results in distinct individual response in gene expression.

•

Non-coding DNA is an important site for such functional genetic variation.

•

Gene expression is a major modulator of brain chemistry and thus behavior.

Over 98% of our genome is non-coding and is now recognised to have a major role in orchestrating the tissue specific and stimulus inducible gene expression pattern which underpins our wellbeing and mental health. The non-coding genome responds functionally to our environment at all levels, encompassing the span from psychological to physiological challenge. The gene expression pattern, termed the transcriptome, ultimately gives us our neurochemistry. Therefore a major modulator of mental wellbeing is how our genes are regulated in response to life experiences. Superimposed on the aforementioned non-coding DNA framework is a vast body of genetic variation in the elements that control response to challenges. These differences, termed polymorphisms, allow for a differential response from a specific DNA element to the same challenge thus potentially allowing ‘individuality’ in the modulation of our transcriptome. This review will focus on a fundamental mechanism defining our psychological and psychiatric wellbeing, namely how genetic variation can be correlated with differential gene expression in response to specific challenges, thus resulting in altered neurochemistry which consequently may shape behaviour.