PASTAA: identifying transcription factors associated with sets of co-regulated genes

Her2 amplification, Rel-A, and Bach1 can influence APOBEC3A expression in breast cancer cells

APOBEC-induced mutations occur in 50% of sequenced human tumors, with APOBEC3A (A3A) being a major contributor to mutagenesis in breast cancer cells. The mechanisms that cause A3A activation and mutagenesis in breast cancers are still unknown. Here, we describe factors that influence basal A3A mRNA transcript levels in breast cancer cells. We found that basal A3A mRNA correlates with A3A protein levels and predicts the amount of APOBEC signature mutations in a panel of breast cancer cell lines, indicating that increased basal transcription may be one mechanism leading to breast cancer mutagenesis. We also show that alteration of ERBB2 expression can drive A3A mRNA levels, suggesting the enrichment of the APOBEC mutation signature in Her2-enriched breast cancer could in part result from elevated A3A transcription. Hierarchical clustering of transcripts in primary breast cancers determined that A3A mRNA was co-expressed with other genes functioning in viral restriction and interferon responses. However, reduction of STAT signaling via inhibitors or shRNA in breast cancer cell lines had only minor impact on A3A abundance. Analysis of single cell RNA-seq from primary tumors indicated that A3A mRNA was highest in infiltrating immune cells within the tumor, indicating that correlations of A3A with STAT signaling in primary tumors may be result from higher immune infiltrates and are not reflective of STAT signaling controlling A3A expression in breast cancer cells. Analysis of ATAC-seq data in multiple breast cancer cell lines identified two transcription factor sites in the APOBEC3A promoter region that could promote A3A transcription. We determined that Rel-A, and Bach1, which have binding sites in these peaks, elevated basal A3A expression. Our findings highlight a complex and variable set of transcriptional activators for A3A in breast cancer cells.

Construction of ferroptosis-related prediction model for pathogenesis, diagnosis and treatment of ruptured abdominal aortic aneurysm

Abdominal aortic aneurysm (AAA) is a dangerous cardiovascular disease, which often brings great psychological burden and economic pressure to patients. If AAA rupture occurs, it is a serious threat to patients' lives. Therefore, it is of clinical value to actively explore the pathogenesis of ruptured AAA and prevent its occurrence. Ferroptosis is a new type of cell death dependent on lipid peroxidation, which plays an important role in many cardiovascular diseases. In this study, we used online data and analysis of ferroptosis-related genes to uncover the formation of ruptured AAA and potential therapeutic targets. We obtained ferroptosis-related differentially expressed genes (Fe-DEGs) from GSE98278 dataset and 259 known ferroptosis-related genes from FerrDb website. Enrichment analysis of differentially expressed genes (DEGs) was performed by gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG). Receiver Operating characteristic (ROC) curve was employed to evaluate the diagnostic abilities of Fe-DEGs. Transcription factors and miRNAs of Fe-DEGs were identified through PASTAA and miRDB, miRWalk, TargetScan respectively. Single-sample gene set enrichment analysis (ssGSEA) was used to observe immune infiltration between the stable group and the rupture group. DGIdb database was performed to find potential targeted drugs of DEGs. GO and KEGG enrichment analysis found that DEGs mainly enriched in "cellular divalent inorganic cation homeostasis," "cellular zinc ion homeostasis," "divalent inorganic cation homeostasis," "Mineral absorption," "Cytokine - cytokine receptor interaction," "Coronavirus disease - COVID-19." Two up-regulated Fe-DEGs MT1G and DDIT4 were found to further analysis. Both single and combined applications of MT1G and DDIT4 showed good diagnostic efficacy (AUC = 0.8254, 0.8548, 0.8577, respectively). Transcription factors STAT1 and PU1 of MT1G and ARNT and MAX of DDIT4 were identified. Meanwhile, has_miR-548p-MT1G pairs, has_miR-53-3p/has_miR-181b-5p/ has_miR-664a-3p-DDIT4 pairs were found. B cells, NK cells, Th2 cells were high expression in the rupture group compared with the stable group, while DCs, Th1 cells were low expression in the rupture group. Targeted drugs against immunity, GEMCITABINE and INDOMETHACIN were discovered. We preliminarily explored the clinical significance of Fe-DEGs MT1G and DDIT4 in the diagnosis of ruptured AAA, and proposed possible upstream regulatory transcription factors and miRNAs. In addition, we also analyzed the immune infiltration of stable and rupture groups, and found possible targeted drugs for immunotherapy.

mir-98-5p regulates gluconeogenesis and lipogenesis by targeting PPP1R15B in hepatocytes

Several reports suggest that circulatory miRNAs are deregulated in diverse diseases and used as markers for disease diagnosis and prognosis. Here we show that miR-98-5p, that is down-regulated in the circulation during diabetes, regulates hepatic gluconeogenesis and lipogenesis by targeting PPP1R15B. miR-98-5p overexpression significantly decreased the transcript and protein levels of PPP1R15B in hepatic HepG2 cells and increased p-eIF2α expression and these were prevented in the presence of its inhibitor. Two major hepatic hallmarks during diabetes i.e. hepatic lipid accumulation and glucose output were explored towards physiological relevance. As compared to scramble, overexpression of miR-98-5p decreased the transcript levels of both gluconeogenic and lipogenic genes together with a significant reduction in hepatic glucose production and fat accumulation in HepG2 cells. Using PASTAA to detect common transcription factors regulating these altered genes, CREB emerged as the most significantly enriched transcription factor. While miR-98-5p overexpression did not change the transcript levels of CREB, there was a significant change in its protein levels. While similar effects on gluconeogenic and lipogenic gene expression were detected using the PPP1R15B siRNA, the opposite was observed in the presence of miR-98-5p inhibitor alone. All these suggest that by targeting PPP1R15B, miR-98-5p regulates hepatic steatosis and glucose output; deregulation of which are characteristic hepatic features during diabetes. Therapeutic intervention of the miR-98/PPP1R15B axis might offer a potential strategy to target aberrant hepatic metabolism during diabetes.

Circular RNA circPLOD2 regulates pericyte function by targeting the transcription factor KLF4

Circular RNAs are generated by backsplicing and control cellular signaling and phenotypes. Pericytes stabilize capillary structures and play important roles in the formation and maintenance of blood vessels. Here, we characterize hypoxia-regulated circular RNAs (circRNAs) in human pericytes and show that the circular RNA of procollagen-lysine,2-oxoglutarate 5-dioxygenase-2 (circPLOD2) is induced by hypoxia and regulates pericyte functions. Silencing of circPLOD2 affects pericytes and increases proliferation, migration, and secretion of soluble angiogenic proteins, thereby enhancing endothelial migration and network capability. Transcriptional and epigenomic profiling of circPLOD2-depleted cells reveals widespread changes in gene expression and identifies the transcription factor krüppel-like factor 4 (KLF4) as a key effector of the circPLOD2-mediated changes. KLF4 depletion mimics circPLOD2 silencing, whereas KLF4 overexpression reverses the effects of circPLOD2 depletion on proliferation and endothelial-pericyte interactions. Together, these data reveal an important function of circPLOD2 in controlling pericyte proliferation and capillary formation and show that the circPLOD2-mediated regulation of KLF4 significantly contributes to the transcriptional response to hypoxia.

Patterns of transcription factor binding and epigenome at promoters allow interpretable predictability of multiple functions of non-coding and coding genes

Understanding the biological roles of all genes only through experimental methods is challenging. A computational approach with reliable interpretability is needed to infer the function of genes, particularly for non-coding RNAs. We have analyzed genomic features that are present across both coding and non-coding genes like transcription factor (TF) and cofactor ChIP-seq (823), histone modifications ChIP-seq (n = 621), cap analysis gene expression (CAGE) tags (n = 255), and DNase hypersensitivity profiles (n = 255) to predict ontology-based functions of genes. Our approach for gene function prediction was reliable (>90% balanced accuracy) for 486 gene-sets. PubMed abstract mining and CRISPR screens supported the inferred association of genes with biological functions, for which our method had high accuracy. Further analysis revealed that TF-binding patterns at promoters have high predictive strength for multiple functions. TF-binding patterns at the promoter add an unexplored dimension of explainable regulatory aspects of genes and their functions. Therefore, we performed a comprehensive analysis for the functional-specificity of TF-binding patterns at promoters and used them for clustering functions to reveal many latent groups of gene-sets involved in common major cellular processes. We also showed how our approach could be used to infer the functions of non-coding genes using the CRISPR screens of coding genes, which were validated using a long non-coding RNA CRISPR screen. Thus our results demonstrated the generality of our approach by using gene-sets from CRISPR screens. Overall, our approach opens an avenue for predicting the involvement of non-coding genes in various functions.

Relationships between some transcription factors and concordantly expressed drought stress-related genes in bread wheat

The challenge of climate change makes it mandatory to improve tolerance to drought stress in bread wheat (Triticum aestivum) via biotechnological approaches. Drought stress experiment was conducted followed by RNA-Seq analysis for leaves of two wheat cultivars namely Giza 168 and Gemmiza 10 with contrasting genotypes. Expression patterns of the regulated stress-related genes and concordantly expressed TFs were detected, then, validated via qPCR for two loss-of-function mutants in Arabidopsis background harboring mutated genes analogue to those in wheat. Drought-stress related genes were searched for concordantly expressed TFs and a total of eight TFs were shown to coexpress with 14 stress-related genes. Among these genes, one TF belongs to the zinc finger protein CONSTANS family and proved via qPCR to drive expression of a gene encoding a speculative TF namely zinc transporter 3-like and two other stress related genes encoding tryptophan synthase alpha chain and asparagine synthetase. Known functions of the two TFs under drought stress complement those of the two concordantly expressed stress-related genes, thus, it is likely that they are related. This study highlights the possibility to utilize metabolic engineering approaches to decipher and incorporate existing regulatory frameworks under drought stress in future breeding programs of bread wheat.

HAPSTR1 localizes HUWE1 to the nucleus to limit stress signaling pathways

HUWE1 is a large, enigmatic HECT-domain ubiquitin ligase implicated in the regulation of diverse pathways, including DNA repair, apoptosis, and differentiation. How HUWE1 engages its structurally diverse substrates and how HUWE1 activity is regulated are unknown. Using unbiased quantitative proteomics, we find that HUWE1 targets substrates in a largely cell-type-specific manner. However, we identify C16orf72/HAPSTR1 as a robust HUWE1 substrate in multiple cell lines. Previously established physical and genetic interactions between HUWE1 and HAPSTR1 suggest that HAPSTR1 positively regulates HUWE1 function. Here, we show that HAPSTR1 is required for HUWE1 nuclear localization and nuclear substrate targeting. Nuclear HUWE1 is required for both cell proliferation and modulation of stress signaling pathways, including p53 and nuclear factor κB (NF-κB)-mediated signaling. Combined, our results define a role for HAPSTR1 in gating critical nuclear HUWE1 functions.

Deep characterization of human γδ T cell subsets defines shared and lineage-specific traits

Under non-pathological conditions, human γδ T cells represent a small fraction of CD3+ T cells in peripheral blood (1-10%). They constitute a unique subset of T lymphocytes that recognize stress ligands or non-peptide antigens through MHC-independent presentation. Major human γδ T cell subsets, Vδ1 and Vδ2, expand in response to microbial infection or malignancy, but possess distinct tissue localization, antigen recognition, and effector responses. We hypothesized that differences at the gene, phenotypic, and functional level would provide evidence that γδ T cell subpopulations belong to distinct lineages. Comparisons between each subset and the identification of the molecular determinants that underpin their differences has been hampered by experimental challenges in obtaining sufficient numbers of purified cells. By utilizing a stringent FACS-based isolation method, we compared highly purified human Vδ1 and Vδ2 cells in terms of phenotype, gene expression profile, and functional responses. We found distinct genetic and phenotypic signatures that define functional differences in γδ T cell populations. Differences in TCR components, repertoire, and responses to calcium-dependent pathways suggest that Vδ1 and Vδ2 T cells are different lineages. These findings will facilitate further investigation into the ligand specificity and unique role of Vδ1 and Vδ2 cells in early immune responses.

Targeting HIC1/TGF-β axis-shaped prostate cancer microenvironment restrains its progression

Prostate cancer (PCa) is a malignant tumor that seriously threatens men's health worldwide. Recently, stromal cells in the tumor microenvironment (TME) have been reported to contribute to the progression of PCa. However, the role and mechanism of how PCa cells interact with stromal cells to reshape the TME remain largely unknown. Here, using a spontaneous prostate adenocarcinoma (PRAD) model driven by the loss of Pten and Hic1, we found that M2 macrophages markedly infiltrated the stroma of Pten and Hic1 double conditional knockout (dCKO) mice compared with those in control (Ctrl) mice due to higher TGF-β levels secreted by HIC1-deleted PCa cells. Mechanistically, TGF-β in TME promoted the polarization of macrophages into "M2" status by activating the STAT3 pathway and modulating c-Myc to upregulate CXCR4 expression. Meanwhile, TGF-β activated the fibroblasts to form cancer-associated fibroblasts (CAFs) that secrete higher CXCL12 levels, which bound to its cognate receptor CXCR4 on M2 macrophages. Upon interaction with CAFs, M2 macrophages secreted more CXCL5, which promoted the epithelial-mesenchymal transition (EMT) of PCa via CXCR2. Moreover, using the TGF-β receptor I antagonist, galunisertib, significantly inhibited the tumor growth and progression of the TRAMP-C1 cell line-derived subcutaneous tumor model. Finally, we confirmed that the stromal microenvironment was shaped by TGF-β in HIC1-deficient PCa and was associated with the progression of PCa.

Exome Sequencing Reveals Genetic Variability and Identifies Chronic Prognostic Loci in Chinese Sarcoidosis Patients

Background

Sarcoidosis is an inflammatory disease characterized by non-caseating granuloma formation in various organs, with several recognized genetic and environmental risk factors. Despite substantial progress, the genetic determinants associated with its prognosis remain largely unknown.

Objectives

This study aimed to identify the genetic changes involved in sarcoidosis and evaluate their clinical relevance.

Methods

We performed whole-exome sequencing (WES) in 116 sporadic sarcoidosis patients (acute sarcoidosis patients, n=58; chronic sarcoidosis patients, n=58). In addition, 208 healthy controls were selected from 1000 G East Asian population data. To identify genes enriched in sarcoidosis, Fisher exact tests were performed. The identified genes were included for further pathway analysis using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Additionally, we used the STRING database to construct a protein network of rare variants and Cytoscape to identify hub genes of signaling pathways.

Results

WES and Fisher’s exact test identified 1,311 variants in 439 protein-coding genes. A total of 135 single nucleotide polymorphisms (SNPs) on 30 protein-coding genes involved in the immunological process based on the GO and KEGG enrichment analysis. Pathway enrichment analysis showed osteoclast differentiation and cytokine–cytokine receptor interactions. Three missense mutations (rs76740888, rs149664918, and rs78251590) in two genes (PRSS3 and CNN2) of immune-related genes showed significantly different mutation frequencies between the disease group and healthy controls. The correlation of genetic abnormalities with clinical outcomes using multivariate analysis of the clinical features and mutation loci showed that the missense variant (rs76740888, Chr9:33796673 G>A) of PRSS3 [p=0.04, odds ratio (OR) = 2.49] was significantly associated with chronic disease prognosis. Additionally, the top two hub genes were CCL4 and CXCR4 based on protein–protein interaction (PPI) network analysis.

Conclusion

Our study provides new insights into the molecular pathogenesis of sarcoidosis and identifies novel genetic alterations in this disease, especially PRSS3, which may be promising targets for future therapeutic strategies for chronic sarcoidosis.

β-Klotho promotes glycolysis and glucose-stimulated insulin secretion via GP130

Impaired glucose-stimulated insulin secretion (GSIS) is a hallmark of type-2 diabetes. However, cellular signaling machineries that control GSIS remain incompletely understood. Here, we report that β-klotho (KLB), a single-pass transmembrane protein known as a co-receptor for fibroblast growth factor 21 (FGF21), fine tunes GSIS via modulation of glycolysis in pancreatic β-cells independent of the actions of FGF21. β-cell-specific deletion of Klb but not Fgf21 deletion causes defective GSIS and glucose intolerance in mice and defective GSIS in islets of type-2 diabetic mice is mitigated by adenovirus-mediated restoration of KLB. Mechanistically, KLB interacts with and stabilizes the cytokine receptor subunit GP130 by blockage of ubiquitin-dependent lysosomal degradation, thereby facilitating interleukin-6-evoked STAT3-HIF1α signaling, which in turn transactivates a cluster of glycolytic genes for adenosine triphosphate production and GSIS. The defective glycolysis and GSIS in Klb-deficient islets are rescued by adenovirus-mediated replenishment of STAT3 or HIF1α. Thus, KLB functions as a key cell-surface regulator of GSIS by coupling the GP130 receptor signaling to glucose catabolism in β-cells and represents a promising therapeutic target for diabetes.

A HIF-1 signature dominates the attenuation in the human skeletal muscle transcriptional response to high-intensity interval training

High-intensity interval training (HIIT) generates profound metabolic adaptations in skeletal muscle. These responses mirror performance improvements but follow a non-linear pattern comprised of an initial fast phase followed by a gradual plateau effect. The complete time-dependent molecular sequelae that regulates this plateau effect remains unknown. We hypothesize that the plateau effect during HIIT is restricted to specific pathways with communal upstream transcriptional regulation. To investigate this, eleven healthy men performed nine sessions of HIIT (10x4 minutes of cycling at 91 % of HR_max) over a 3-week period. Before and 3h after the 1^st and 9^th exercise bout, skeletal muscle biopsies were obtained, and RNA sequencing performed. Almost 2,000 genes across 84 pathways were differentially expressed in response to a single HIIT session. The overall transcriptional response to acute exercise was strikingly similar at 3 weeks, 83 % (n=1650) of the genes regulated after the 1^st bout of exercise were similarly regulated by the 9^th bout, albeit with a smaller effect size, and the response attenuated to on average 70 % of the 1^st bout. The attenuation differed substantially between pathways and was very pronounced for glycolysis and cellular adhesion but more preserved for MAPK and VEGF-A signalling. The attenuation was driven by a combination of changes in steady-state expression and specific transcriptional regulation. Given that the exercise intensity was progressively increased, and that the attenuation was pathway specific, we suggest that moderation of muscular adaptation after a period of training stems from targeted regulation rather than a diminished exercise stimulus.

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.

PTPN2 Regulates the Interferon Signaling and Endoplasmic Reticulum Stress Response in Pancreatic β-Cells in Autoimmune Diabetes

Type 1 diabetes (T1D) results from autoimmune destruction of β-cells in the pancreas. Protein tyrosine phosphatases (PTPs) are candidate genes for T1D and play a key role in autoimmune disease development and β-cell dysfunction. Here, we assessed the global protein and individual PTP profiles in the pancreas from nonobese mice with early-onset diabetes (NOD) mice treated with an anti-CD3 monoclonal antibody and interleukin-1 receptor antagonist. The treatment reversed hyperglycemia, and we observed enhanced expression of PTPN2, a PTP family member and T1D candidate gene, and endoplasmic reticulum (ER) chaperones in the pancreatic islets. To address the functional role of PTPN2 in β-cells, we generated PTPN2-deficient human stem cell-derived β-like and EndoC-βH1 cells. Mechanistically, we demonstrated that PTPN2 inactivation in β-cells exacerbates type I and type II interferon signaling networks and the potential progression toward autoimmunity. Moreover, we established the capacity of PTPN2 to positively modulate the Ca2+-dependent unfolded protein response and ER stress outcome in β-cells. Adenovirus-induced overexpression of PTPN2 partially protected from ER stress-induced β-cell death. Our results postulate PTPN2 as a key protective factor in β-cells during inflammation and ER stress in autoimmune diabetes.

Synaptic disruption and CREB-regulated transcription are restored by K+ channel blockers in ALS

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease, which is still missing effective therapeutic strategies. Although manipulation of neuronal excitability has been tested in murine and human ALS models, it is still under debate whether neuronal activity might represent a valid target for efficient therapies. In this study, we exploited a combination of transcriptomics, proteomics, optogenetics and pharmacological approaches to investigate the activity-related pathological features of iPSC-derived C9orf72-mutant motoneurons (MN). We found that human ALSC9orf72 MN are characterized by accumulation of aberrant aggresomes, reduced expression of synaptic genes, loss of synaptic contacts and a dynamic "malactivation" of the transcription factor CREB. A similar phenotype was also found in TBK1-mutant MN and upon overexpression of poly(GA) aggregates in primary neurons, indicating a strong convergence of pathological phenotypes on synaptic dysregulation. Notably, these alterations, along with neuronal survival, could be rescued by treating ALS-related neurons with the K+ channel blockers Apamin and XE991, which, respectively, target the SK and the Kv7 channels. Thus, our study shows that restoring the activity-dependent transcriptional programme and synaptic composition exerts a neuroprotective effect on ALS disease progression.

Impact of Interleukin 10 Deficiency on Intestinal Epithelium Responses to Inflammatory Signals

Interleukin 10 (IL-10) is a pleiotropic, anti-inflammatory cytokine that has a major protective role in the intestine. Although its production by cells of the innate and adaptive immune system has been extensively studied, its intrinsic role in intestinal epithelial cells is poorly understood. In this study, we utilised both ATAC sequencing and RNA sequencing to define the transcriptional response of murine enteroids to tumour necrosis factor (TNF). We identified that the key early phase drivers of the transcriptional response to TNF within intestinal epithelium were NFκB transcription factor dependent. Using wild-type and Il10-/- enteroid cultures, we showed an intrinsic, intestinal epithelium specific effect of IL-10 deficiency on TNF-induced gene transcription, with significant downregulation of identified NFκB target genes Tnf, Ccl20, and Cxcl10, and delayed overexpression of NFκB inhibitor encoding genes, Nfkbia and Tnfaip3. IL-10 deficiency, or immunoblockade of IL-10 receptor, impacted on TNF-induced endogenous NFκB activity and downstream NFκB target gene transcription. Intestinal epithelium-derived IL-10 appears to play a crucial role as a positive regulator of the canonical NFκB pathway, contributing to maintenance of intestinal homeostasis. This is particularly important in the context of an inflammatory environment and highlights the potential for future tissue-targeted IL-10 therapeutic intervention.

Meta-Analysis of Heifer Traits Identified Reproductive Pathways in Bos indicus Cattle

Fertility traits measured early in life define the reproductive potential of heifers. Knowledge of genetics and biology can help devise genomic selection methods to improve heifer fertility. In this study, we used ~2400 Brahman cattle to perform GWAS and multi-trait meta-analysis to determine genomic regions associated with heifer fertility. Heifer traits measured were pregnancy at first mating opportunity (PREG1, a binary trait), first conception score (FCS, score 1 to 3) and rebreeding score (REB, score 1 to 3.5). The heritability estimates were 0.17 (0.03) for PREG1, 0.11 (0.05) for FCS and 0.28 (0.05) for REB. The three traits were highly genetically correlated (0.75-0.83) as expected. Meta-analysis was performed using SNP effects estimated for each of the three traits, adjusted for standard error. We identified 1359 significant SNPs (p-value < 9.9 × 10-6 at FDR < 0.0001) in the multi-trait meta-analysis. Genomic regions of 0.5 Mb around each significant SNP from the meta-analysis were annotated to create a list of 2560 positional candidate genes. The most significant SNP was in the vicinity of a genomic region on chromosome 8, encompassing the genes SLC44A1, FSD1L, FKTN, TAL2 and TMEM38B. The genomic region in humans that contains homologs of these genes is associated with age at puberty in girls. Top significant SNPs pointed to additional fertility-related genes, again within a 0.5 Mb region, including ESR2, ITPR1, GNG2, RGS9BP, ANKRD27, TDRD12, GRM1, MTHFD1, PTGDR and NTNG1. Functional pathway enrichment analysis resulted in many positional candidate genes relating to known fertility pathways, including GnRH signaling, estrogen signaling, progesterone mediated oocyte maturation, cAMP signaling, calcium signaling, glutamatergic signaling, focal adhesion, PI3K-AKT signaling and ovarian steroidogenesis pathway. The comparison of results from this study with previous transcriptomics and proteomics studies on puberty of the same cattle breed (Brahman) but in a different population identified 392 genes in common from which some genes-BRAF, GABRA2, GABR1B, GAD1, FSHR, CNGA3, PDE10A, SNAP25, ESR2, GRIA2, ORAI1, EGFR, CHRNA5, VDAC2, ACVR2B, ORAI3, CYP11A1, GRIN2A, ATP2B3, CAMK2A, PLA2G, CAMK2D and MAPK3-are also part of the above-mentioned pathways. The biological functions of the positional candidate genes and their annotation to known pathways allowed integrating the results into a bigger picture of molecular mechanisms related to puberty in the hypothalamus-pituitary-ovarian axis. A reasonable number of genes, common between previous puberty studies and this study on early reproductive traits, corroborates the proposed molecular mechanisms. This study identified the polymorphism associated with early reproductive traits, and candidate genes that provided a visualization of the proposed mechanisms, coordinating the hypothalamic, pituitary, and ovarian functions for reproductive performance in Brahman cattle.

Disruption of innate defense responses by endoglycosidase HPSE promotes cell survival

The drive to withstand environmental stresses and defend against invasion is a universal trait extant in all forms of life. While numerous canonical signaling cascades have been characterized in detail, it remains unclear how these pathways interface to generate coordinated responses to diverse stimuli. To dissect these connections, we followed heparanase (HPSE), a protein best known for its endoglycosidic activity at the extracellular matrix but recently recognized to drive various forms of late-stage disease through unknown mechanisms. Using herpes simplex virus-1 (HSV-1) infection as a model cellular perturbation, we demonstrate that HPSE acts beyond its established enzymatic role to restrict multiple forms of cell-intrinsic defense and facilitate host cell reprogramming by the invading pathogen. We reveal that cells devoid of HPSE are innately resistant to infection and counteract viral takeover through multiple amplified defense mechanisms. With a unique grasp of the fundamental processes of transcriptional regulation and cell death, HPSE represents a potent cellular intersection with broad therapeutic potential.

Computational prediction of CRISPR-impaired non-coding regulatory regions

Genome-wide CRISPR screens are becoming more widespread and allow the simultaneous interrogation of thousands of genomic regions. Although recent progress has been made in the analysis of CRISPR screens, it is still an open problem how to interpret CRISPR mutations in non-coding regions of the genome. Most of the tools concentrate on the interpretation of mutations introduced in gene coding regions. We introduce a computational pipeline that uses epigenomic information about regulatory elements for the interpretation of CRISPR mutations in non-coding regions. We illustrate our analysis protocol on the analysis of a genome-wide CRISPR screen in hTERT-RPE1 cells and reveal novel regulatory elements that mediate chemoresistance against doxorubicin in these cells. We infer links to established and to novel chemoresistance genes. Our analysis protocol is general and can be applied on any cell type and with different CRISPR enzymes.

Dissociation of DNA damage sensing by endoglycosidase HPSE

Balance between cell proliferation and elimination is critical in handling threats both exogenous and of internal dysfunction. Recent work has implicated a conserved but poorly understood endoglycosidase heparanase (HPSE) in the restriction of innate defense responses, yet biochemical mediators of these key functions remained unclear. Here, an unbiased immunopurification proteomics strategy is employed to identify and rank uncharacterized interactions between HPSE and mediators of canonical signaling pathways linking cell cycle and stress responses. We demonstrate with models of genotoxic stress including herpes simplex virus infection and chemotherapeutic treatment that HPSE dampens innate responses to double-stranded DNA breakage by interfering with signal transduction between initial sensors and downstream mediators. Given the long-standing recognition of HPSE in driving late-stage inflammatory disease exemplified by tissue destruction and cancer metastasis, modulation of this protein with control over the DNA damage response imparts a unique strategy in the development of unconventional multivalent therapy.

•

HPSE binds key proteins at interface of DNA damage signaling and IFN responses

•

Nuclear translocation of DNA damage transducer ATM is enhanced in absence of HPSE

•

Cells lacking HPSE display enhanced sensitivity to DNA damage-induced death

•

HPSE interfaces with regulators of DNA damage response to influence cell fate

Single cell RNA sequencing of human microglia uncovers a subset associated with Alzheimer's disease

The extent of microglial heterogeneity in humans remains a central yet poorly explored question in light of the development of therapies targeting this cell type. Here, we investigate the population structure of live microglia purified from human cerebral cortex samples obtained at autopsy and during neurosurgical procedures. Using single cell RNA sequencing, we find that some subsets are enriched for disease-related genes and RNA signatures. We confirm the presence of four of these microglial subpopulations histologically and illustrate the utility of our data by characterizing further microglial cluster 7, enriched for genes depleted in the cortex of individuals with Alzheimer's disease (AD). Histologically, these cluster 7 microglia are reduced in frequency in AD tissue, and we validate this observation in an independent set of single nucleus data. Thus, our live human microglia identify a range of subtypes, and we prioritize one of these as being altered in AD.

Application of the RBBP9 Serine Hydrolase Inhibitor, ML114, Decouples Human Pluripotent Stem Cell Proliferation and Differentiation

Retinoblastoma binding protein 9 (RBBP9) is required for maintaining the expression of both pluripotency and cell cycle genes in human pluripotent stem cells (hPSCs). An siRNA-based study from our group showed it does so by influencing cell cycle progression through the RB/E2F pathway. In non-pluripotent cells, RBBP9 is also known to have serine hydrolase (SH) activity, acting on currently undefined target proteins. The role of RBBP9 SH activity in hPSCs, and during normal development, is currently unknown. To begin assessing whether RBBP9 SH activity might contribute to hPSC maintenance, hPSCs were treated with ML114—a selective chemical inhibitor of RBBP9 SH activity. Stem cells treated with ML114 showed significantly reduced population growth rate, colony size and progression through the cell cycle, with no observable change in cell morphology or decrease in pluripotency antigen expression—suggesting no initiation of hPSC differentiation. Consistent with this, hPSCs treated with ML114 retained the capacity for tri-lineage differentiation, as seen through teratoma formation. Subsequent microarray and Western blot analyses of ML114-treated hPSCs suggest the nuclear transcription factor Y subunit A (NFYA) may be a candidate effector of RBBP9 SH activity in hPSCs. These data support a role for RBBP9 in regulating hPSC proliferation independent of differentiation, whereby inhibition of RBBP9 SH activity de-couples decreased hPSC proliferation from initiation of differentiation.

Gene regulation could be attributed to TCF3 and other key transcription factors in the muscle of pubertal heifers

Puberty is a whole-body event, driven by the hypothalamic integration of peripheral signals such as leptin or IGF-1. In the process of puberty, reproductive development is simultaneous to growth, including muscle growth. To enhance our understanding of muscle function related to puberty, we performed transcriptome analyses of muscle samples from six pre- and six post-pubertal Brahman heifers (Bos indicus). Our aims were to perform differential expression analyses and co-expression analyses to derive a regulatory gene network associate with puberty. As a result, we identified 431 differentially expressed (DEx) transcripts (genes and non-coding RNAs) when comparing pre- to post-pubertal average gene expression. The DEx transcripts were compared with all expressed transcripts in our samples (over 14,000 transcripts) for functional enrichment analyses. The DEx transcripts were associated with "extracellular region," "inflammatory response" and "hormone activity" (adjusted p < .05). Inflammatory response for muscle regeneration is a necessary aspect of muscle growth, which is accelerated during puberty. The term "hormone activity" may signal genes that respond to progesterone signalling in the muscle, as the presence of this hormone is an important difference between pre- and post-pubertal heifers in our experimental design. The DEx transcript with the highest average expression difference was a mitochondrial gene, ENSBTAG00000043574 that might be another important link between energy metabolism and puberty. In the derived co-expression gene network, we identified six hub genes: CDC5L, MYC, TCF3, RUNX2, ATF2 and CREB1. In the same network, 48 key regulators of DEx transcripts were identified, using a regulatory impact factor metric. The hub gene TCF3 was also a key regulator. The majority of the key regulators (22 genes) are members of the zinc finger family, which has been implicated in bovine puberty in other tissues. In conclusion, we described how puberty may affect muscle gene expression in cattle.

BATF acts as an essential regulator of IL-25-responsive migratory ILC2 cell fate and function

A transitory, interleukin-25 (IL-25)-responsive, group 2 innate lymphoid cell (ILC2) subset induced during type 2 inflammation was recently identified as iILC2s. This study focuses on understanding the significance of this population in relation to tissue-resident nILC2s in the lung and intestine. RNA-sequencing and pathway analysis revealed the AP-1 superfamily transcription factor BATF (basic leucine zipper transcription factor, activating transcription factor-like) as a potential modulator of ILC2 cell fate. Infection of BATF-deficient mice with Nippostrongylus brasiliensis showed a selective defect in IL-25-mediated helminth clearance and a corresponding loss of iILC2s in the lung characterized as IL-17RB^high, KLRG1^high, BATF^high, and Arg1^low BATF deficiency selectively impaired iILC2s because it had no impact on tissue-resident nILC2 frequency or function. Pulmonary-associated iILC2s migrated to the lung after infection, where they represented an early source of IL-4 and IL-13. Although the composition of ILC2s in the small intestine was distinct from those in the lung, their frequency and IL-13 expression remained dependent on BATF, which was also required for optimal goblet and tuft cell hyperplasia. Findings support IL-25-responsive ILC2s as early sentinels of mucosal barrier integrity.

MEIRLOP: improving score-based motif enrichment by incorporating sequence bias covariates

BACKGROUND

Motif enrichment analysis (MEA) identifies over-represented transcription factor binding (TF) motifs in the DNA sequence of regulatory regions, enabling researchers to infer which transcription factors can regulate transcriptional response to a stimulus, or identify sequence features found near a target protein in a ChIP-seq experiment. Score-based MEA determines motifs enriched in regions exhibiting extreme differences in regulatory activity, but existing methods do not control for biases in GC content or dinucleotide composition. This lack of control for sequence bias, such as those often found in CpG islands, can obscure the enrichment of biologically relevant motifs.

RESULTS

We developed Motif Enrichment In Ranked Lists of Peaks (MEIRLOP), a novel MEA method that determines enrichment of TF binding motifs in a list of scored regulatory regions, while controlling for sequence bias. In this study, we compare MEIRLOP against other MEA methods in identifying binding motifs found enriched in differentially active regulatory regions after interferon-beta stimulus, finding that using logistic regression and covariates improves the ability to call enrichment of ISGF3 binding motifs from differential acetylation ChIP-seq data compared to other methods. Our method achieves similar or better performance compared to other methods when quantifying the enrichment of TF binding motifs from ENCODE TF ChIP-seq datasets. We also demonstrate how MEIRLOP is broadly applicable to the analysis of numerous types of NGS assays and experimental designs.

CONCLUSIONS

Our results demonstrate the importance of controlling for sequence bias when accurately identifying enriched DNA sequence motifs using score-based MEA. MEIRLOP is available for download from https://github.com/npdeloss/meirlop under the MIT license.

Neuronal Apolipoprotein E4 Expression Results in Proteome-Wide Alterations and Compromises Bioenergetic Capacity by Disrupting Mitochondrial Function

Apolipoprotein (apo) E4, the major genetic risk factor for Alzheimer's disease (AD), alters mitochondrial function and metabolism early in AD pathogenesis. When injured or stressed, neurons increase apoE synthesis. Because of its structural difference from apoE3, apoE4 undergoes neuron-specific proteolysis, generating fragments that enter the cytosol, interact with mitochondria, and cause neurotoxicity. However, apoE4's effect on mitochondrial respiration and metabolism is not understood in detail. Here we used biochemical assays and proteomic profiling to more completely characterize the effects of apoE4 on mitochondrial function and cellular metabolism in Neuro-2a neuronal cells stably expressing apoE4 or apoE3. Under basal conditions, apoE4 impaired respiration and increased glycolysis, but when challenged or stressed, apoE4-expressing neurons had 50% less reserve capacity to generate ATP to meet energy requirements than apoE3-expressing neurons. ApoE4 expression also decreased the NAD+/NADH ratio and increased the levels of reactive oxygen species and mitochondrial calcium. Global proteomic profiling revealed widespread changes in mitochondrial processes in apoE4 cells, including reduced levels of numerous respiratory complex subunits and major disruptions to all detected subunits in complex V (ATP synthase). Also altered in apoE4 cells were levels of proteins related to mitochondrial endoplasmic reticulum-associated membranes, mitochondrial fusion/fission, mitochondrial protein translocation, proteases, and mitochondrial ribosomal proteins. ApoE4-induced bioenergetic deficits led to extensive metabolic rewiring, but despite numerous cellular adaptations, apoE4-expressing neurons remained vulnerable to metabolic stress. Our results provide insights into potential molecular targets of therapies to correct apoE4-associated mitochondrial dysfunction and altered cellular metabolism.

MiR-302b as a Combinatorial Therapeutic Approach to Improve Cisplatin Chemotherapy Efficacy in Human Triple-Negative Breast Cancer

Introduction: Chemotherapy is still the standard of care for triple-negative breast cancers (TNBCs). Here, we investigated miR-302b as a therapeutic tool to enhance cisplatin sensitivity in vivo and unraveled the molecular mechanism. Materials and Methods: TNBC-xenografted mice were treated with miR-302b or control, alone or with cisplatin. Genome-wide transcriptome analysis and independent-validation of Integrin Subunit Alpha 6 (ITGA6) expression was assessed on mice tumor samples. Silencing of ITGA6 was performed to evaluate cisplatin response in vitro. Further, potential transcription factors of ITGA6 (E2F transcription facor 1 (E2F1), E2F transcription factor 2 (E2F2), and Yin Yang 1 (YY1)) were explored to define the miRNA molecular mechanism. The miR-302b expression was also assessed in TNBC patients treated with chemotherapy. Results: The miR–302b-cisplatin combination significantly impaired tumor growth versus the control through indirect ITGA6 downregulation. Indeed, ITGA6 was downmodulated in mice treated with miR-302b–cisplatin, and ITGA6 silencing increased drug sensitivity in TNBC cells. In silico analyses and preclinical assays pointed out the regulatory role of the E2F family and YY1 on ITGA6 expression under miR-302b–cisplatin treatment. Finally, miR-302b enrichment correlated with better overall survival in 118 TNBC patients. Conclusion: MiR-302b can be exploited as a new therapeutic tool to improve the response to chemotherapy, modulating the E2F family, YY1, and ITGA6 expression. Moreover, miR-302b could be defined as a new prognostic factor in TNBC patients.

Mapping Gene Circuits Essential for Germ Layer Differentiation via Loss-of-Function Screens in Haploid Human Embryonic Stem Cells

Pluripotent stem cells can differentiate into all embryonic germ layers, yet the genes essential for these cell fate transitions in human remain elusive. Here, we mapped the essential genes for the differentiation of human pluripotent stem cells (hPSCs) into the three germ layers by using a genome-wide loss-of-function library established in haploid hPSCs. Strikingly, we observed a high fraction of essential genes associated with plasma membrane, highlighting signaling pathways needed for each lineage differentiation. Interestingly, analysis of all hereditary neurological disorders uncovered high essentiality among microcephaly-causing genes. Furthermore, we demonstrated lineage-specific hierarchies among essential transcription factors and a set of Golgi- and endoplasmic reticulum-related genes needed for the differentiation into all germ layers. Our work sheds light on the gene networks regulating early gastrulation events in human by defining essential drivers of specific embryonic germ layer fates and essential genes for the exit from pluripotency.

MS2-TRIBE Evaluates Both Protein-RNA Interactions and Nuclear Organization of Transcription by RNA Editing

Both UV-cross-linking and immunoprecipitation (CLIP) and RNA editing (TRIBE) can identify the targets of RNA-binding proteins. To evaluate false-positives of CLIP and TRIBE, endogenous β-actin mRNA was tagged with MS2 stem loops, making it the only bona fide target mRNA for the MS2 capsid protein (MCP). CLIP and TRIBE detected β-actin, albeit with false-positives. False-positive CLIP signals were attributed to nonspecific antibody interactions. In contrast, putative false-positive TRIBE targets were genes spatially proximal to the β-actin gene. MCP-ADAR edited nearby nascent transcripts consistent with interchromosomal contacts observed in Hi-C. The identification of nascent contacts implies RNA regulatory proteins (e.g., splicing factors) associated with multiple nascent transcripts, forming domains of post-transcriptional activity. Repeating these results with an integrated inducible MS2 reporter indicated that MS2-TRIBE can be applied to a broad array of cells and transcripts to study spatial organization and nuclear RNA regulation.

Identification of Target Genes in Hypertension and Left Ventricular Remodeling

INTRODUCTION

Hypertension occurs profoundly in the world, and left ventricular (LV) remodeling containing functional, structural, and mechanical changes induced by uncontrolled blood pressure is a well-known complication, however the underlying mechanism is still obscure.

METHODS

To determine differences in gene expression profiles of hypertension and LV remodeling consequence to hypertension, Gene Expression Omnibus 2R online tool was used to identify differently expressed genes. Publicly available databases including GeneMANIA, database for annotation, visualization and integrated discovery, search tool for the retrieva predicting associated transcription factors (TF) from annotated affinities interacting genes, Predicting Associated TF from Annotated Affinities, JASPAR and Comparative Toxicogenomics Database (CTD) were accessed to perform an integrated bioinformatic analysis.

RESULTS

Twenty-one genes (SEC14L3, EML7, PSMD7, PSMA1, GLRX, CNOT10, NBR1, DUSP12, STRAP, SMIM14, RBM8A, TMEM59, TMEM87A,PSMC1, CASP4, ITGB8, DNAJA1, PINK1, PRNP, SAP30L, and EIF3M) were found overexpression in both hypertension and hypertensive LV remodeling. Biological process analysis first revealed that enrichment of these target genes correlated with regulation of cellular amino acid metabolic process, antigen processing and presentation of exogenous peptide antigen via MHC class I, TAP-dependent and proteasome complex, 3 different expression genes (DEGs) participate significantly enriched in NFκB, WNT, and MAPK pathways, meanwhile, 47% DEGs displayed similar co-expression characteristics. Furthermore, the transcription factors associated with key DEGs were identified. Finally, the TF (HAND1, E4BP4, ESR1, VBP, ELK-1, POU3F2) associated with LV remodeling in hypertension were confirmed to act a crucial role in correlated heart diseases.

CONCLUSION

The present study reveals the targeted genes probably associated with LV remodeling in hypertension by bioinformatics-based analyses, which provides clues for prognosis judgement and pharmacological therapies.

Characterization of a Novel Compound That Stimulates STING-Mediated Innate Immune Activity in an Allele-Specific Manner

The innate immune response to cytosolic DNA involves transcriptional activation of type I interferons (IFN-I) and proinflammatory cytokines. This represents the culmination of intracellular signaling pathways that are initiated by pattern recognition receptors that engage DNA and require the adaptor protein Stimulator of Interferon Genes (STING). These responses lead to the generation of cellular and tissue states that impair microbial replication and facilitate the establishment of long-lived, antigen-specific adaptive immunity. Ultimately this can lead to immune-mediated protection from infection but also to the cytotoxic T cell-mediated clearance of tumor cells. Intriguingly, pharmacologic activation of STING-dependent phenotypes is known to enhance both vaccine-associated immunogenicity and immune-based anti-tumor therapies. Unfortunately, the STING protein exists as multiple variant forms in the human population that exhibit differences in their reactivity to chemical stimuli and in the intensity of molecular signaling they induce. In light of this, STING-targeting drug discovery efforts require an accounting of protein variant-specific activity. Herein we describe a small molecule termed M04 that behaves as a novel agonist of human STING. Importantly, we find that the molecule exhibits a differential ability to activate STING based on the allelic variant examined. Furthermore, while M04 is inactive in mice, expression of human STING in mouse cells rescues reactivity to the compound. Using primary human cells in ex vivo assays we were also able to show that M04 is capable of simulating innate responses important for adaptive immune activation such as cytokine secretion, dendritic cell maturation, and T cell cross-priming. Collectively, this work demonstrates the conceivable utility of a novel agonist of human STING both as a research tool for exploring STING biology and as an immune potentiating molecule.

EpiRegio: analysis and retrieval of regulatory elements linked to genes

A current challenge in genomics is to interpret non-coding regions and their role in transcriptional regulation of possibly distant target genes. Genome-wide association studies show that a large part of genomic variants are found in those non-coding regions, but their mechanisms of gene regulation are often unknown. An additional challenge is to reliably identify the target genes of the regulatory regions, which is an essential step in understanding their impact on gene expression. Here we present the EpiRegio web server, a resource of regulatory elements (REMs). REMs are genomic regions that exhibit variations in their chromatin accessibility profile associated with changes in expression of their target genes. EpiRegio incorporates both epigenomic and gene expression data for various human primary cell types and tissues, providing an integrated view of REMs in the genome. Our web server allows the analysis of genes and their associated REMs, including the REM's activity and its estimated cell type-specific contribution to its target gene's expression. Further, it is possible to explore genomic regions for their regulatory potential, investigate overlapping REMs and by that the dissection of regions of large epigenomic complexity. EpiRegio allows programmatic access through a REST API and is freely available at https://epiregio.de/.

IGAP-integrative genome analysis pipeline reveals new gene regulatory model associated with nonspecific TF-DNA binding affinity

The human genome is regulated in a multi-dimensional way. While biophysical factors like Non-specific Transcription factor Binding Affinity (nTBA) act at DNA sequence level, other factors act above sequence levels such as histone modifications and 3-D chromosomal interactions. This multidimensionality of regulation requires many of these factors for a proper understanding of the regulatory landscape of the human genome. Here, we propose a new biophysical model for estimating nTBA. Integration of nTBA with chromatin modifications and chromosomal interactions, using a new Integrative Genome Analysis Pipeline (IGAP), reveals additive effects of nTBA to regulatory DNA sequences and identifies three types of genomic zones in the human genome (Inactive Genomic Zones, Poised Genomic Zones, and Active Genomic Zones). It also unveils a novel long distance gene regulatory model: chromosomal interactions reduce the physical distance between the high occupancy target (HOT) regions that results in high nTBA to DNA in the area, which in turn attract TFs to such regions with higher binding potential. These findings will help to elucidate the three-dimensional diffusion process that TFs use during their search for the right targets.

Seq-ing answers: Current data integration approaches to uncover mechanisms of transcriptional regulation

Advancements in the field of next generation sequencing lead to the generation of ever-more data, with the challenge often being how to combine and reconcile results from different OMICs studies such as genome, epigenome and transcriptome. Here we provide an overview of the standard processing pipelines for ChIP-seq and RNA-seq as well as common downstream analyses. We describe popular multi-omics data integration approaches used to identify target genes and co-factors, and we discuss how machine learning techniques may predict transcriptional regulators and gene expression.

Functional effects of variation in transcription factor binding highlight long-range gene regulation by epromoters

Identifying DNA cis-regulatory modules (CRMs) that control the expression of specific genes is crucial for deciphering the logic of transcriptional control. Natural genetic variation can point to the possible gene regulatory function of specific sequences through their allelic associations with gene expression. However, comprehensive identification of causal regulatory sequences in brute-force association testing without incorporating prior knowledge is challenging due to limited statistical power and effects of linkage disequilibrium. Sequence variants affecting transcription factor (TF) binding at CRMs have a strong potential to influence gene regulatory function, which provides a motivation for prioritizing such variants in association testing. Here, we generate an atlas of CRMs showing predicted allelic variation in TF binding affinity in human lymphoblastoid cell lines and test their association with the expression of their putative target genes inferred from Promoter Capture Hi-C and immediate linear proximity. We reveal >1300 CRM TF-binding variants associated with target gene expression, the majority of them undetected with standard association testing. A large proportion of CRMs showing associations with the expression of genes they contact in 3D localize to the promoter regions of other genes, supporting the notion of 'epromoters': dual-action CRMs with promoter and distal enhancer activity.

Testing Proximity of Genomic Regions to Transcription Start Sites and Enhancers Complements Gene Set Enrichment Testing

Large sets of genomic regions are generated by the initial analysis of various genome-wide sequencing data, such as ChIP-seq and ATAC-seq experiments. Gene set enrichment (GSE) methods are commonly employed to determine the pathways associated with them. Given the pathways and other gene sets (e.g., GO terms) of significance, it is of great interest to know the extent to which each is driven by binding near transcription start sites (TSS) or near enhancers. Currently, no tool performs such an analysis. Here, we present a method that addresses this question to complement GSE methods for genomic regions. Specifically, the new method tests whether the genomic regions in a gene set are significantly closer to a TSS (or to an enhancer) than expected by chance given the total list of genomic regions, using a non-parametric test. Combining the results from a GSE test with our novel method provides additional information regarding the mode of regulation of each pathway, and additional evidence that the pathway is truly enriched. We illustrate our new method with a large set of ENCODE ChIP-seq data, using the chipenrich Bioconductor package. The results show that our method is a powerful complementary approach to help researchers interpret large sets of genomic regions.

Improved linking of motifs to their TFs using domain information

Motivation

A central aim of molecular biology is to identify mechanisms of transcriptional regulation. Transcription factors (TFs), which are DNA-binding proteins, are highly involved in these processes, thus a crucial information is to know where TFs interact with DNA and to be aware of the TFs’ DNA-binding motifs. For that reason, computational tools exist that link DNA-binding motifs to TFs either without sequence information or based on TF-associated sequences, e.g. identified via a chromatin immunoprecipitation followed by sequencing (ChIP-seq) experiment.

In this paper, we present MASSIF, a novel method to improve the performance of existing tools that link motifs to TFs relying on TF-associated sequences. MASSIF is based on the idea that a DNA-binding motif, which is correctly linked to a TF, should be assigned to a DNA-binding domain (DBD) similar to that of the mapped TF. Because DNA-binding motifs are in general not linked to DBDs, it is not possible to compare the DBD of a TF and the motif directly. Instead we created a DBD collection, which consist of TFs with a known DBD and an associated motif. This collection enables us to evaluate how likely it is that a linked motif and a TF of interest are associated to the same DBD. We named this similarity measure domain score, and represent it as a P-value. We developed two different ways to improve the performance of existing tools that link motifs to TFs based on TF-associated sequences: (i) using meta-analysis to combine P-values from one or several of these tools with the P-value of the domain score and (ii) filter unlikely motifs based on the domain score.

Results

We demonstrate the functionality of MASSIF on several human ChIP-seq datasets, using either motifs from the HOCOMOCO database or de novo identified ones as input motifs. In addition, we show that both variants of our method improve the performance of tools that link motifs to TFs based on TF-associated sequences significantly independent of the considered DBD type.

Availability and implementation

MASSIF is freely available online at https://github.com/SchulzLab/MASSIF.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome-wide DNA copy number profiling and bioinformatics analysis of ovarian cancer reveals key genes and pathways associated with distinct invasive/migratory capabilities

Ovarian cancer (OC) metastasis presents major hurdles that must be overcome to improve patient outcomes. Recent studies have demonstrated copy number variations (CNVs) frequently contribute to alterations in oncogenic drivers. The present study used a CytoScan HD Array to analyse CNVs and loss of heterozygosity (LOH) in the entire genomes of 6 OC patients and human OC cell lines to determine the genetic target events leading to the distinct invasive/migratory capacities of OC. The results showed that LOH at Xq11.1 and Xp21.1 and gains at 8q21.13 were novel, specific CNVs. Ovarian cancer-related CNVs were then screened by bioinformatics analysis. In addition, transcription factors-target gene interactions were predicted with information from PASTAA analysis. As a result, six genes (i.e., GAB2, AKT1, EGFR, COL6A3, UGT1A1 and UGT1A8) were identified as strong candidates by integrating the above data with gene expression and clinical outcome data. In the transcriptional regulatory network, 4 known cancer-related transcription factors (TFs) interacted with 6 CNV-driven genes. The protein/DNA arrays revealed 3 of these 4 TFs as potential candidate gene-related transcription factors in OC. We then demonstrated that these six genes can serve as potential biomarkers for OC. Further studies are required to elucidate the pathogenesis of OC.

Multifunctional transcriptional coactivator PC4 is a global co-regulator of p53-dependent stress response and gene regulation

Human positive coactivator 4 (PC4), a multifunctional chromatin-associated protein, is known to directly interact with p53 and modulate expressions of a few p53-dependent genes. However, the role of PC4 in p53's myriad of other regulatory functions is not known. The p53-PC4 interaction was selectively perturbed by a small peptide which led to abrogation of genotoxic stress-induced up-regulation of many p53-dependent genes and reduction of apoptosis in A549 cells. Over-expression of a PC4 point mutant, incapable of binding p53, recapitulated many of the effects of the peptide. Global gene expression profiling in A549 cells, upon peptide treatment, revealed PC4's involvement in the regulation of many p53-dependent pathways, including the Hippo pathway. Introduction of the peptide in neuronal cells significantly reduced its amyloid-β-induced death. Thus, PC4 emerges as a global co-regulator of p53 and a therapeutic target against pathogeneses where the p53-dependent cell death process plays a crucial role.

Hierarchical cooperation of transcription factors from integration analysis of DNA sequences, ChIP-Seq and ChIA-PET data

Background

Chromosomal architecture, which is constituted by chromatin loops, plays an important role in cellular functions. Gene expression and cell identity can be regulated by the chromatin loop, which is formed by proximal or distal enhancers and promoters in linear DNA (1D). Enhancers and promoters are fundamental non-coding elements enriched with transcription factors (TFs) to form chromatin loops. However, the specific cooperation of TFs involved in forming chromatin loops is not fully understood.

Results

Here, we proposed a method for investigating the cooperation of TFs in four cell lines by the integrative analysis of DNA sequences, ChIP-Seq and ChIA-PET data. Results demonstrate that the interaction of enhancers and promoters is a hierarchical and dynamic complex process with cooperative interactions of different TFs synergistically regulating gene expression and chromatin structure. The TF cooperation involved in maintaining and regulating the chromatin loop of cells can be regulated by epigenetic factors, such as other TFs and DNA methylation.

Conclusions

Such cooperation among TFs provides the potential features that can affect chromatin’s 3D architecture in cells. The regulation of chromatin 3D organization and gene expression is a complex process associated with the hierarchical and dynamic prosperities of TFs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5535-2) contains supplementary material, which is available to authorized users.

Stems cells, big data and compendium-based analyses for identifying cell types, signalling pathways and gene regulatory networks

Identification of new drug and cell therapy targets for disease treatment will be facilitated by a detailed molecular understanding of normal and disease development. Human pluripotent stem cells can provide a large in vitro source of human cell types and, in a growing number of instances, also three-dimensional multicellular tissues called organoids. The application of stem cell technology to discovery and development of new therapies will be aided by detailed molecular characterisation of cell identity, cell signalling pathways and target gene networks. Big data or 'omics' techniques-particularly transcriptomics and proteomics-facilitate cell and tissue characterisation using thousands to tens-of-thousands of genes or proteins. These gene and protein profiles are analysed using existing and/or emergent bioinformatics methods, including a growing number of methods that compare sample profiles against compendia of reference samples. This review assesses how compendium-based analyses can aid the application of stem cell technology for new therapy development. This includes via robust definition of differentiated stem cell identity, as well as elucidation of complex signalling pathways and target gene networks involved in normal and diseased states.

Identification of Potential Gene Interactions in Heart Failure Caused by Idiopathic Dilated Cardiomyopathy

Background

Many heart failure (HF) cases are caused by idiopathic dilated cardiomyopathy (iDCM). This study explored the mechanisms of the development and progression of HF caused by iDCM.

Material/Methods

The gene expression profiles of 102 samples were downloaded from the GEO database (GSE5406). Differentially expressed genes (DEGs) were identified through GO analysis and a KEGG pathway analysis, respectively. A protein–protein interaction (PPI) network was constructed and analyzed to screen potential regulatory proteins. In addition, MCODE and a cytoHubba plugin were used to identify the module and hub genes of DEGs. Finally, transcription factors (TFs) were predicted using PASTAA. We did not perform whole-exome sequencing (WES) for detecting mitochondrial DNA (mtDNA).

Results

A total of 197 DEGs were screened, and 3 modules, and 4 upregulated and 11 downregulated hub genes were screened. The GO analysis focused on the terms and 12 KEGG pathways were enriched. The FOS, TIMP1, and SERPINE1 hub genes, as well as some key TFs, demonstrated important roles in the progression of HF caused by iDCM. CEBPD, CEBOB, CDC37L1, and SRGN may be new targets for HF in iDCM patients.

Conclusions

The identified DEGs and their enriched pathways provide references for exploring the mechanisms of the development and progression of HF patients with iDCM. Moreover, modules, hub genes, and TFs may be useful in the treatment and diagnosis of HF patients with iDCM. However, mtDNA was not investigated.

Lung eosinophilia induced by house dust mites or ovalbumin is modulated by nicotinic receptor α7 and inhibited by cigarette smoke

Eosinophilia (EOS) is an important component of airway inflammation and hyperresponsiveness in allergic reactions including those leading to asthma. Although cigarette smoking (CS) is a significant contributor to long-term adverse outcomes in these lung disorders, there are also the curious reports of its ability to produce acute suppression of inflammatory responses including EOS through poorly understood mechanisms. One possibility is that proinflammatory processes are suppressed by nicotine in CS acting through nicotinic receptor α7 (α7). Here we addressed the role of α7 in modulating EOS with two mouse models of an allergic response: house dust mites (HDM; Dermatophagoides sp.) and ovalbumin (OVA). The influence of α7 on EOS was experimentally resolved in wild-type mice or in mice in which a point mutation of the α7 receptor (α7E260A:G) selectively restricts normal signaling of cellular responses. RNA analysis of alveolar macrophages and the distal lung epithelium indicates that normal α7 function robustly impacts gene expression in the epithelium to HDM and OVA but to different degrees. Notable was allergen-specific α7 modulation of Ccl11 and Ccl24 (eotaxins) expression, which was enhanced in HDM but suppressed in OVA EOS. CS suppressed EOS induced by both OVA and HDM, as well as the inflammatory genes involved, regardless of α7 genotype. These results suggest that EOS in response to HDM or OVA is through signaling pathways that are modulated in a cell-specific manner by α7 and are distinct from CS suppression.

Identification of target genes in cardiomyopathy with fibrosis and cardiac remodeling

Background

Identify genes probably associated with chronic heart failure and predict potential target genes for dilated cardiomyopathy using bioinformatics analyses.

Methods

Gene expression profiles (series number GSE3585 and GSE42955) of cardiomyopathy patients and healthy controls were downloaded from the Expression Omnibus Gene (GEO) database. Differential expression of genes (DEGS) between the two groups of total 14 cardiomyopathy patients and 10 healthy controls were subsequently identified by limma package of R. Database for Annotation, Visualization, and Integrated Discovery (DAVID Tool), which is an analysis of enriched biological processes. Search Tool for the Retrieval Interacting Genes (STRING) was used as well for the analysis of protein-protein interaction network (PPI). Prediction of the potential drugs was suggested based on the preliminarily identified genes using Connectivity Map (CMap).

Results

Eighty-nine DEGs were identified (57 up-regulated and 32 down-regulated). The most enrichment Gene Ontology (GO) terms (P < 0.05) contain genes involved in extracellular matrix (ECM) and biological adhesion signal pathways (P < 0.05, ES > 1.5) such as ECM-receptors, focal adhesion and transforming growth factor beta (TGF-β), etc. Fifty-one differentially expressed genes were found to encode interacting proteins. Eleven key genes along with related transcription factors were identified including CTGF, POSTN, CORIN, FIGF, etc.

Conclusion

Bioinformatics-based analyses reveal the targeted genes probably associated with cardiomyopathy, which provide clues for pharmacological therapies aiming at the targets.

An improved compound Poisson model for the number of motif hits in DNA sequences

Motivation

Transcription factors play a crucial role in gene regulation by binding to specific regulatory sequences. The sequence motifs recognized by a transcription factor can be described in terms of position frequency matrices. When scanning a sequence for matches to a position frequency matrix, one needs to determine a cut-off, which then in turn results in a certain number of hits. In this paper we describe how to compute the distribution of match scores and of the number of motif hits, which are the prerequisites to perform motif hit enrichment analysis.

Results

We put forward an improved compound Poisson model that supports general order-d Markov background models and which computes the number of motif-hits more accurately than earlier models. We compared the accuracy of the improved compound Poisson model with previously proposed models across a range of parameters and motifs, demonstrating the improvement. The importance of the order-d model is supported in a case study using CpG-island sequences.

Availability and implementation

The method is available as a Bioconductor package named ’motifcounter’ https://bioconductor.org/packages/motifcounter.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genomic approaches for the elucidation of genes and gene networks underlying cardiovascular traits

Genome-wide association studies have shed light on the association between natural genetic variation and cardiovascular traits. However, linking a cardiovascular trait associated locus to a candidate gene or set of candidate genes for prioritization for follow-up mechanistic studies is all but straightforward. Genomic technologies based on next-generation sequencing technology nowadays offer multiple opportunities to dissect gene regulatory networks underlying genetic cardiovascular trait associations, thereby aiding in the identification of candidate genes at unprecedented scale. RNA sequencing in particular becomes a powerful tool when combined with genotyping to identify loci that modulate transcript abundance, known as expression quantitative trait loci (eQTL), or loci modulating transcript splicing known as splicing quantitative trait loci (sQTL). Additionally, the allele-specific resolution of RNA-sequencing technology enables estimation of allelic imbalance, a state where the two alleles of a gene are expressed at a ratio differing from the expected 1:1 ratio. When multiple high-throughput approaches are combined with deep phenotyping in a single study, a comprehensive elucidation of the relationship between genotype and phenotype comes into view, an approach known as systems genetics. In this review, we cover key applications of systems genetics in the broad cardiovascular field.

Candidate genes and pathways downstream of PAX8 involved in ovarian high-grade serous carcinoma

Understanding the biology and molecular pathogenesis of ovarian epithelial cancer (EOC) is key to developing improved diagnostic and prognostic indicators and effective therapies. Although research has traditionally focused on the hypothesis that high-grade serous carcinoma (HGSC) arises from the ovarian surface epithelium (OSE), recent studies suggest that additional sites of origin exist and a substantial proportion of cases may arise from precursor lesions located in the Fallopian tubal epithelium (FTE). In FTE cells, the transcription factor PAX8 is a marker of the secretory cell lineage and its expression is retained in 96% of EOC. We have recently reported that PAX8 is involved in the tumorigenic phenotype of ovarian cancer cells. In this study, to uncover genes and pathways downstream of PAX8 involved in ovarian carcinoma we have determined the molecular profiles of ovarian cancer cells and in parallel of Fallopian tube epithelial cells by means of a silencing approach followed by an RNA-seq analysis. Interestingly, we highlighted the involvement of pathways like WNT signaling, epithelial-mesenchymal transition, p53 and apoptosis. We believe that our analysis has led to the identification of candidate genes and pathways regulated by PAX8 that could be additional targets for the therapy of ovarian carcinoma.

Profiling of the transcriptional response to all-trans retinoic acid in breast cancer cells reveals RARE-independent mechanisms of gene expression

Retinoids, derivatives of vitamin A, are key physiological molecules with regulatory effects on cell differentiation, proliferation and apoptosis. As a result, they are of interest for cancer therapy. Specifically, models of breast cancer have varied responses to manipulations of retinoid signaling. This study characterizes the transcriptional response of MDA-MB-231 and MDA-MB-468 breast cancer cells to retinaldehyde dehydrogenase 1A3 (ALDH1A3) and all-trans retinoic acid (atRA). We demonstrate limited overlap between ALDH1A3-induced gene expression and atRA-induced gene expression in both cell lines, suggesting that the function of ALDH1A3 in breast cancer progression extends beyond its role as a retinaldehyde dehydrogenase. Our data reveals divergent transcriptional responses to atRA, which are largely independent of genomic retinoic acid response elements (RAREs) and consistent with the opposing responses of MDA-MB-231 and MDA-MB-468 to in vivo atRA treatment. We identify transcription factors associated with each gene set. Manipulation of the IRF1 transcription factor demonstrates that it is the level of atRA-inducible and epigenetically regulated transcription factors that determine expression of target genes (e.g. CTSS, cathepsin S). This study provides a paradigm for complex responses of breast cancer models to atRA treatment, and illustrates the need to characterize RARE-independent responses to atRA in a variety of models.

Benzotriazoles Reactivate Latent HIV-1 through Inactivation of STAT5 SUMOylation

The presence of latent HIV-1 in infected individuals represents a major barrier preventingviral eradication. For that reason, reactivation of latent viruses in the presence of antiretroviral regimens has been proposed as a therapeutic strategy to achieve remission. We screened for small molecules and identified several benzotriazole derivatives with the ability to reactivate latent HIV-1. In the presence of IL-2, benzotriazoles reactivated and reduced the latent reservoir in primary cells, and, remarkably, viral reactivation was achieved without inducing cell proliferation, T cell activation, or cytokine release. Mechanistic studies showed that benzotriazoles block SUMOylation of phosphorylated STAT5, increasing STAT5’s activity and occupancy of the HIV-1 LTR. Our results identify benzotriazoles as latency reversing agents and STAT5 signaling and SUMOylation as targets for HIV-1 eradication strategies. These compounds represent a different direction in the search for “shock and kill” therapies.