Histone Acetylome-wide Association Study of Autism Spectrum Disorder.

The association of histone modification changes with autism spectrum disorder (ASD) has not been systematically examined. We conducted a histone acetylome-wide association study (HAWAS) by performing H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq) on 257 postmortem samples from ASD and matched control brains. Despite etiological heterogeneity, ≥68% of syndromic and idiopathic ASD cases shared a common acetylome signature at >5,000 cis-regulatory elements in prefrontal and temporal cortex. Similarly, multiple genes associated with rare genetic mutations in ASD showed common "epimutations." Acetylome aberrations in ASD were not attributable to genetic differentiation at cis-SNPs and highlighted genes involved in synaptic transmission, ion transport, epilepsy, behavioral abnormality, chemokinesis, histone deacetylation, and immunity. By correlating histone acetylation with genotype, we discovered >2,000 histone acetylation quantitative trait loci (haQTLs) in human brain regions, including four candidate causal variants for psychiatric diseases. Due to the relative stability of histone modifications postmortem, we anticipate that the HAWAS approach will be applicable to multiple diseases.

The E-Id Protein Axis Specifies Adaptive Lymphoid Cell Identity and Suppresses Thymic Innate Lymphoid Cell Development.

Innate and adaptive lymphoid development is orchestrated by the activities of E proteins and their antagonist Id proteins, but how these factors regulate early T cell progenitor (ETP) and innate lymphoid cell (ILC) development remains unclear. Using multiple genetic strategies, we demonstrated that E proteins E2A and HEB acted in synergy in the thymus to establish T cell identity and to suppress the aberrant development of ILCs, including ILC2s and lymphoid-tissue-inducer-like cells. E2A and HEB orchestrated T cell fate and suppressed the ILC transcription signature by activating the expression of genes associated with Notch receptors, T cell receptor (TCR) assembly, and TCR-mediated signaling. E2A and HEB acted in ETPs to establish and maintain a T-cell-lineage-specific enhancer repertoire, including regulatory elements associated with the Notch1, Rag1, and Rag2 loci. On the basis of these and previous observations, we propose that the E-Id protein axis specifies innate and adaptive lymphoid cell fate.

Regulome-wide association study identifies genetically driven accessible regions associated with pancreatic cancer risk.

Genome-wide association studies (GWAS) have identified two dozen genetic variants that are associated with the risk of pancreatic ductal adenocarcinoma (PDAC), a deadly malignancy. However, a majority of these variants are located in noncoding regions of the genome, which limits the translation of GWAS findings into clinical applications. The regulome-wide association study (RWAS) is a recently developed method for identifying TF binding-induced accessibility regions for diseases. However, their potential connection to PDAC has yet to be fully explored. We evaluated the associations between genetically predicted levels of chromatin accessibility and risk of PDAC by using pan-cancer chromatin accessibility genetic prediction models. Our analysis included 8275 cases and 6723 controls from the PanScan (I, II, and III) and PanC4 consortia. To further refine our results, we also integrated genes associated to allele-specific accessibility quantitative trait loci (as-aQTL) and TF motifs located in the as-aQTL. We found that 50 chromatin accessibility features were associated with PDAC risk at a false discovery rate (FDR) of less than 0.05. A total of 28 RWAS peaks were identified as conditionally significant. By integrating the results from as-aQTL, motif analysis, and RWAS, we identified candidate causal regulatory elements for two potential chromatin accessibility regions (THCA_89956 and ESCA_89167) that are associated with PDAC risk. Our study identified chromatin accessibility features in noncoding genomic regions that are associated with PDAC risk. We also predicted the associated genes and disrupt motifs. Our findings provide new insights into the regulatory mechanisms of noncoding regions for pancreatic tumorigenesis.

De novo and inherited variants in coding and regulatory regions in genetic cardiomyopathies.

BACKGROUND: Cardiomyopathies are a leading cause of progressive heart failure and sudden cardiac death; however, their genetic aetiology remains poorly understood. We hypothesised that variants in noncoding regulatory regions and oligogenic inheritance mechanisms may help close the diagnostic gap. METHODS: We first analysed whole-genome sequencing data of 143 parent-offspring trios from Genomics England 100,000 Genomes Project. We used gene panel testing and a phenotype-based, variant prioritisation framework called Exomiser to identify candidate genes in trios. To assess the contribution of noncoding DNVs to cardiomyopathies, we intersected DNVs with open chromatin sequences from single-cell ATAC-seq data of cardiomyocytes. We also performed a case-control analysis in an exome-negative cohort, including 843 probands and 19,467 controls, to assess the association between noncoding variants in known cardiomyopathy genes and disease. RESULTS: In the trio analysis, a definite or probable genetic diagnosis was identified in 21 probands according to the American College of Medical Genetics guidelines. We identified novel DNVs in diagnostic-grade genes (RYR2, TNNT2, PTPN11, MYH7, LZR1, NKX2-5), and five cases harbouring a combination of prioritised variants, suggesting that oligogenic inheritance and genetic modifiers contribute to cardiomyopathies. Phenotype-based ranking of candidate genes identified in noncoding DNV analysis revealed JPH2 as the top candidate. Moreover, a case-control analysis revealed an enrichment of rare noncoding variants in regulatory elements of cardiomyopathy genes (p = .035, OR = 1.43, 95% Cl = 1.095-1.767) versus controls. Of the 25 variants associated with disease  (p< 0.5), 23 are novel and nine are predicted to disrupt transcription factor binding motifs. CONCLUSION: Our results highlight complex genetic mechanisms in cardiomyopathies and reveal novel genes for future investigations.

Pathway perturbations in signaling networks: Linking genotype to phenotype.

Genes and gene products interact with each other to form signal transduction networks in the cell. The interactome networks are under intricate regulation in physiological conditions, but could go awry upon genome instability caused by genetic mutations. In the past decade with next-generation sequencing technologies, an increasing number of genomic mutations have been identified in a variety of disease patients and healthy individuals. As functional and systematic studies on these mutations leap forward, they begin to reveal insights into cellular homeostasis and disease mechanisms. In this review, we discuss recent advances in the field of network biology and signaling pathway perturbations upon genomic changes, and highlight the success of various omics datasets in unraveling genotype-to-phenotype relationships.

The Developmental Gene Hypothesis for Punctuated Equilibrium: Combined Roles of Developmental Regulatory Genes and Transposable Elements.

Theories of the genetics underlying punctuated equilibrium (PE) have been vague to date. Here the developmental gene hypothesis is proposed, which states that: 1) developmental regulatory (DevReg) genes are responsible for the orchestration of metazoan morphogenesis and their extreme conservation and mutation intolerance generates the equilibrium or stasis present throughout much of the fossil record and 2) the accumulation of regulatory elements and recombination within these same genes-often derived from transposable elements-drives punctuated bursts of morphological divergence and speciation across metazoa. This two-part hypothesis helps to explain the features that characterize PE, providing a theoretical genetic basis for the once-controversial theory. Also see the video abstract here https://youtu.be/C-fu-ks5yDs.

[miRNA Regulome in Different Atherosclerosis Phenotypes].

Dysregulation of microRNA (miRNA) expression is associated with a susceptibility to many diseases, including atherosclerotic lesions of the coronary and carotid arteries and the development of clinical complications such as coronary heart disease, myocardial infarction, chronic cerebral ischemia, ischemic stroke. Recently, more and more studies analyze the miRNA regulome including a network of regulatory elements for the expression of miRNAs themselves and targets under their control. The review summarizes the data from articles concerned miRNA expression and changes in DNA methylation in the miRNA genes in human atherosclerotic arteries, as well as with the analysis of the association between single nucleotide polymorphisms and copy number variations in the miRNA genes with clinical complications of atherosclerosis.

The CHRNA5/CHRNA3/CHRNB4 Nicotinic Receptor Regulome: Genomic Architecture, Regulatory Variants, and Clinical Associations.

Functionally related genes often cluster into a genome region under coordinated regulation, forming a local regulome. To understand regulation of the CHRNA5/CHRNA3/CHRNB4 nicotinic receptor gene cluster, we integrate large-scale RNA expression data (brain and peripheral) from GTEx (Genotype Tissue Expression), clinical associations (GRASP), and linkage disequilibrium data (1000 Genomes) to find candidate SNPs representing independent regulatory variants. CHRNA3, CHRNA5, CHRNB4 mRNAs, and a well-expressed CHRNA5 antisense RNA (RP11-650L12.2) are co-expressed in many human tissues, suggesting common regulatory elements. The CHRNA5 enhancer haplotype tagged by rs880395 not only increases CHRNA5 mRNA expression in all tissues, but also enhances RP11-650L12.2 and CHRNA3 expression, suggesting DNA looping to multiple promoters. However, in nucleus accumbens and putamen, but not other brain regions, CHRNA3 expression associates uniquely with a haplotype tagged by rs1948 (located in the CHRNB4 3'UTR). Haplotype/diplotype analysis of rs880395 and rs1948 plus rs16969968 (a nonsynonymous CHRNA5 risk variant) in GWAS (COGEND, UW-TTURC, SAGE) yields a nicotine dependence risk profile only partially captured by rs16969968 alone. An example of local gene clusters, this nicotinic regulome is controlled by complex genetic variation, with broad implications for interpreting GWAS.

Understanding the response to endurance exercise using a systems biology approach: combining blood metabolomics, transcriptomics and miRNomics in horses.

BACKGROUND: Endurance exercise in horses requires adaptive processes involving physiological, biochemical, and cognitive-behavioral responses in an attempt to regain homeostasis. We hypothesized that the identification of the relationships between blood metabolome, transcriptome, and miRNome during endurance exercise in horses could provide significant insights into the molecular response to endurance exercise. For this reason, the serum metabolome and whole-blood transcriptome and miRNome data were obtained from ten horses before and after a 160 km endurance competition. RESULTS: We obtained a global regulatory network based on 11 unique metabolites, 263 metabolic genes and 5 miRNAs whose expression was significantly altered at T1 (post- endurance competition) relative to T0 (baseline, pre-endurance competition). This network provided new insights into the cross talk between the distinct molecular pathways (e.g. energy and oxygen sensing, oxidative stress, and inflammation) that were not detectable when analyzing single metabolites or transcripts alone. Single metabolites and transcripts were carrying out multiple roles and thus sharing several biochemical pathways. Using a regulatory impact factor metric analysis, this regulatory network was further confirmed at the transcription factor and miRNA levels. In an extended cohort of 31 independent animals, multiple factor analysis confirmed the strong associations between lactate, methylene derivatives, miR-21-5p, miR-16-5p, let-7 family and genes that coded proteins involved in metabolic reactions primarily related to energy, ubiquitin proteasome and lipopolysaccharide immune responses after the endurance competition. Multiple factor analysis also identified potential biomarkers at T0 for an increased likelihood for failure to finish an endurance competition. CONCLUSIONS: To the best of our knowledge, the present study is the first to provide a comprehensive and integrated overview of the metabolome, transcriptome, and miRNome co-regulatory networks that may have a key role in regulating the metabolic and immune response to endurance exercise in horses.

Pleiotropic genes for metabolic syndrome and inflammation.

Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular disease. Our study hypothesis is that additional to genes influencing individual MetS risk factors, genetic variants exist that influence MetS and inflammatory markers forming a predisposing MetS genetic network. To test this hypothesis a staged approach was undertaken. (a) We analyzed 17 metabolic and inflammatory traits in more than 85,500 participants from 14 large epidemiological studies within the Cross Consortia Pleiotropy Group. Individuals classified with MetS (NCEP definition), versus those without, showed on average significantly different levels for most inflammatory markers studied. (b) Paired average correlations between 8 metabolic traits and 9 inflammatory markers from the same studies as above, estimated with two methods, and factor analyses on large simulated data, helped in identifying 8 combinations of traits for follow-up in meta-analyses, out of 130,305 possible combinations between metabolic traits and inflammatory markers studied. (c) We performed correlated meta-analyses for 8 metabolic traits and 6 inflammatory markers by using existing GWAS published genetic summary results, with about 2.5 million SNPs from twelve predominantly largest GWAS consortia. These analyses yielded 130 unique SNPs/genes with pleiotropic associations (a SNP/gene associating at least one metabolic trait and one inflammatory marker). Of them twenty-five variants (seven loci newly reported) are proposed as MetS candidates. They map to genes MACF1, KIAA0754, GCKR, GRB14, COBLL1, LOC646736-IRS1, SLC39A8, NELFE, SKIV2L, STK19, TFAP2B, BAZ1B, BCL7B, TBL2, MLXIPL, LPL, TRIB1, ATXN2, HECTD4, PTPN11, ZNF664, PDXDC1, FTO, MC4R and TOMM40. Based on large data evidence, we conclude that inflammation is a feature of MetS and several gene variants show pleiotropic genetic associations across phenotypes and might explain a part of MetS correlated genetic architecture. These findings warrant further functional investigation.

Subolesin knockdown in tick cells provides insights into vaccine protective mechanisms.

Ticks as obligate blood-feeding arthropod vectors of pathogenic viruses, bacteria, protozoa and helminths associated with prevalent tick-borne diseases (TBDs) worldwide. These arthropods constitute the second vector after mosquitoes that transmit pathogens to humans and the first vector in domestic animals. Vaccines constitute the safest and more effective approach to control tick infestations and TBDs, but research is needed to identify new antigens and improve vaccine formulations. The tick protein Subolesin (Sub) is a well-known vaccine protective antigen with a highly conserved sequence at both gene and protein levels in the Ixodidae and among arthropods and vertebrates. In this study, transcriptomics and proteomics analyses were conducted together with graph theory data analysis in wild type and Sub knockdown (KD) tick ISE6 cells in order to identify and characterize the functional implications of Sub in tick cells. The results support a key role for Sub in the regulation of gene expression in ticks and the relevance of this antigen in vaccine development against ticks and TBDs. Proteins with differential representation in response to Sub KD provide insights into vaccine protective mechanisms and candidate tick protective antigens.

The ChIP-Hub Resource: Toward plantEncode.

Recent advances in sequencing technologies lead to the generation of an enormous amount of regulome and epigenome data in a variety of plant species. However, a comprehensive standardized resource is so far not available. In this chapter, we present ChIP-Hub, an integrative platform that has been developed based on the ENCODE standards by collecting and reanalyzing regulatory genomic datasets from 41 plant species. The ChIP-hub website is introduced in this chapter, including information on detailed steps of searching, data download, and online analyses, which facilitates users to explore ChIP-Hub. We also provide a cross-species comparison of chromatin accessibility information that gives a thorough view of evolutionary regulatory networks in plants.

YWHAG Deficiency Disrupts the EMT-Associated Network to Induce Oxidative Cell Death and Prevent Metastasis.

Metastasis involves epithelial-to-mesenchymal transition (EMT), a process that is regulated by complex gene networks, where their deliberate disruption may yield a promising outcome. However, little is known about mechanisms that coordinate these metastasis-associated networks. To address this gap, hub genes with broad engagement across various human cancers by analyzing the transcriptomes of different cancer cell types undergoing EMT are identified. The oncogenic signaling adaptor protein tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein gamma (YWHAG) is ranked top for its clinical relevance and impact. The cellular kinome and transcriptome data are surveyed to construct the regulome of YWHAG, revealing stress responses and metabolic processes during cancer EMT. It is demonstrated that a YWHAG-dependent cytoprotective mechanism in the regulome is embedded in EMT-associated networks to protect cancer cells from oxidative catastrophe through enhanced autophagy during EMT. YWHAG deficiency results in a rapid accumulation of reactive oxygen species (ROS), delayed EMT, and cell death. Tumor allografts show that metastasis potential and overall survival time are correlated with the YWHAG expression level of cancer cell lines. Metastasized tumors have higher expression of YWHAG and autophagy-related genes than primary tumors. Silencing YWHAG diminishes primary tumor volumes, prevents metastasis, and prolongs the median survival period of the mice.

A Differential Transcriptional Regulome Approach to Unpack Cancer Biology: Insights on Renal Cell Carcinoma Subtypes.

Cancer research calls for new approaches that account for the regulatory complexities of biology. We present, in this study, the differential transcriptional regulome (DIFFREG) approach for the identification and prioritization of key transcriptional regulators and apply it to the case of renal cell carcinoma (RCC) biology. Of note, RCC has a poor prognosis and the biomarker and drug discovery studies to date have tended to focus on gene expression independent from mutations and/or post-translational modifications. DIFFREG focuses on the differential regulation between transcription factors (TFs) and their target genes rather than differential gene expression and integrates transcriptome profiling with the human transcriptional regulatory network to analyze differential gene regulation between healthy and RCC cases. In this study, RNA-seq tissue samples (n = 1020) from the Cancer Genome Atlas (TCGA), including healthy and tumor subjects, were integrated with a comprehensive human TF-gene interactome dataset (1122603 interactions between 1289 TFs and 25177 genes). Comparative analysis of DIFFREG profiles, consisting of perturbed TF-gene interactions, from three common subtypes (clear cell RCC, papillary RCC and chromophobe RCC) revealed subtype-specific alterations, supporting the hypothesis that these signatures in the transcriptional regulome profiles may be considered potential biomarkers that may play an important role in elucidating the molecular mechanisms of RCC development and translating knowledge about the genetic basis of RCC into the clinic. In addition, these indicators may help oncologists make the best decisions for diagnosis and prognosis management.

Reverse-ChIP Techniques for Identifying Locus-Specific Proteomes: A Key Tool in Unlocking the Cancer Regulome.

A phenotypic hallmark of cancer is aberrant transcriptional regulation. Transcriptional regulation is controlled by a complicated array of molecular factors, including the presence of transcription factors, the deposition of histone post-translational modifications, and long-range DNA interactions. Determining the molecular identity and function of these various factors is necessary to understand specific aspects of cancer biology and reveal potential therapeutic targets. Regulation of the genome by specific factors is typically studied using chromatin immunoprecipitation followed by sequencing (ChIP-Seq) that identifies genome-wide binding interactions through the use of factor-specific antibodies. A long-standing goal in many laboratories has been the development of a 'reverse-ChIP' approach to identify unknown binding partners at loci of interest. A variety of strategies have been employed to enable the selective biochemical purification of sequence-defined chromatin regions, including single-copy loci, and the subsequent analytical detection of associated proteins. This review covers mass spectrometry techniques that enable quantitative proteomics before providing a survey of approaches toward the development of strategies for the purification of sequence-specific chromatin as a 'reverse-ChIP' technique. A fully realized reverse-ChIP technique holds great potential for identifying cancer-specific targets and the development of personalized therapeutic regimens.

Akirin/Subolesin regulatory mechanisms at host/tick-pathogen interactions.

Ticks and tick-borne pathogens such as Anaplasma phagocytophilum affect human and animal health worldwide and thus the characterization of host/tick-pathogen interactions is important for the control of tick-borne diseases. The vertebrate regulatory proteins Akirins and its tick ortholog, Subolesin, are conserved throughout the metazoan and involved in the regulation of different biological processes such as immune response to pathogen infection. Akirin/Subolesin have a key role in host/tick-pathogen interactions and exert its regulatory function primarily through interacting proteins such as transcription factors, chromatin remodelers and RNA-associated proteins. Recent results have provided evidence of akirin/subolesin genetic interactions and the interaction of Akirin/Subolesin with histones, thus suggesting a role in direct chromatin remodeling. Finally, and still to be proven, some models suggest the possibility of direct Akirin/Subolesin protein interactions with DNA. Future research should advance the characterization of Akirin/Subolesin interactome and its functional role at the host/tick-pathogen interface. These results have implications for translational biotechnology and medicine for the development of new effective interventions for the control of ticks and tick-borne diseases.

Comprehensive enhancer-target gene assignments improve gene set level interpretation of genome-wide regulatory data.

BACKGROUND: Revealing the gene targets of distal regulatory elements is challenging yet critical for interpreting regulome data. Experiment-derived enhancer-gene links are restricted to a small set of enhancers and/or cell types, while the accuracy of genome-wide approaches remains elusive due to the lack of a systematic evaluation. We combined multiple spatial and in silico approaches for defining enhancer locations and linking them to their target genes aggregated across >500 cell types, generating 1860 human genome-wide distal enhancer-to-target gene definitions (EnTDefs). To evaluate performance, we used gene set enrichment (GSE) testing on 87 independent ENCODE ChIP-seq datasets of 34 transcription factors (TFs) and assessed concordance of results with known TF Gene Ontology annotations, and other benchmarks. RESULTS: The top ranked 741 (40%) EnTDefs significantly outperform the common, naïve approach of linking distal regions to the nearest genes, and the top 10 EnTDefs perform well when applied to ChIP-seq data of other cell types. The GSE-based ranking of EnTDefs is highly concordant with ranking based on overlap with curated benchmarks of enhancer-gene interactions. Both our top general EnTDef and cell-type-specific EnTDefs significantly outperform seven independent computational and experiment-based enhancer-gene pair datasets. We show that using our top EnTDefs for GSE with either genome-wide DNA methylation or ATAC-seq data is able to better recapitulate the biological processes changed in gene expression data performed in parallel for the same experiment than our lower-ranked EnTDefs. CONCLUSIONS: Our findings illustrate the power of our approach to provide genome-wide interpretation regardless of cell type.

Bone Marrow Stroma-Induced Transcriptome and Regulome Signatures of Multiple Myeloma.

Multiple myeloma (MM) is a hematological cancer with inevitable drug resistance. MM cells interacting with bone marrow stromal cells (BMSCs) undergo substantial changes in the transcriptome and develop de novo multi-drug resistance. As a critical component in transcriptional regulation, how the chromatin landscape is transformed in MM cells exposed to BMSCs and contributes to the transcriptional response to BMSCs remains elusive. We profiled the transcriptome and regulome for MM cells using a transwell coculture system with BMSCs. The transcriptome and regulome of MM cells from the upper transwell resembled MM cells that coexisted with BMSCs from the lower chamber but were distinctive to monoculture. BMSC-induced genes were enriched in the JAK2/STAT3 signaling pathway, unfolded protein stress, signatures of early plasma cells, and response to proteasome inhibitors. Genes with increasing accessibility at multiple regulatory sites were preferentially induced by BMSCs; these genes were enriched in functions linked to responses to drugs and unfavorable clinic outcomes. We proposed JUNB and ATF4::CEBPß as candidate transcription factors (TFs) that modulate the BMSC-induced transformation of the regulome linked to the transcriptional response. Together, we characterized the BMSC-induced transcriptome and regulome signatures of MM cells to facilitate research on epigenetic mechanisms of BMSC-induced multi-drug resistance in MM.

Tick Importin-α Is Implicated in the Interactome and Regulome of the Cofactor Subolesin.

Ticks and tick-borne diseases (TBDs) represent a burden for human and animal health worldwide. Currently, vaccines constitute the safest and most effective approach to control ticks and TBDs. Subolesin (SUB) has been identified as a vaccine antigen for the control of tick infestations and pathogen infection and transmission. The characterization of the molecular function of SUB and the identification of tick proteins interacting with SUB may provide the basis for the discovery of novel antigens and for the rational design of novel anti-tick vaccines. In the present study, we used the yeast two-hybrid system (Y2H) as an unbiased approach to identify tick SUB-interacting proteins in an Ixodes ricinus cDNA library, and studied the possible role of SUB as a chromatin remodeler through direct interaction with histones. The Y2H screening identified Importin-α as a potential SUB-interacting protein, which was confirmed in vitro in a protein pull-down assay. The sub gene expression levels in tick midgut and fat body were significantly higher in unfed than fed female ticks, however, the importin-α expression levels did not vary between unfed and fed ticks but tended to be higher in the ovary when compared to those in other organs. The effect of importin-α RNAi was characterized in I. ricinus under artificial feeding conditions. Both sub and importin-α gene knockdown was observed in all tick tissues and, while tick weight was significantly lower in sub RNAi-treated ticks than in controls, importin-α RNAi did not affect tick feeding or oviposition, suggesting that SUB is able to exert its function in the absence of Importin-α. Furthermore, SUB was shown to physically interact with histone 4, which was corroborated by protein pull-down and western blot analysis. These results confirm that by interacting with numerous tick proteins, SUB is a key cofactor of the tick interactome and regulome. Further studies are needed to elucidate the nature of the SUB-Importin-α interaction and the biological processes and functional implications that this interaction may have.