TRANSFAC: transcriptional regulation, from patterns to profiles

Notch1 blockade by a novel, selective anti-Notch1 neutralizing antibody improves immunotherapy efficacy in melanoma by promoting an inflamed TME

BACKGROUND

Immune checkpoint inhibitors (ICI) have dramatically improved the life expectancy of patients with metastatic melanoma. However, about half of the patient population still present resistance to these treatments. We have previously shown Notch1 contributes to a non-inflamed TME in melanoma that reduces the response to ICI. Here, we addressed the therapeutic effects of a novel anti-Notch1 neutralizing antibody we produced, alone and in combination with immune checkpoint inhibition in melanoma models.

METHODS

Anti-Notch1 was designed to interfere with ligand binding. Mice were immunized with a peptide encompassing EGF-like repeats 11-15 of human Notch1, the minimal required region that allows ligand binding and Notch1 activation. Positive clones were expanded and tested for neutralizing capabilities. Anti-Notch1-NIC was used to determine whether anti-Notch1 was able to reduce Notch1 cleavage; while anti-SNAP23 and BCAT2 were used as downstream Notch1 and Notch2 targets, respectively. K457 human melanoma cells and the YUMM2.1 and 1.7 syngeneic mouse melanoma cells were used. Cell death after anti-Notch1 treatment was determined by trypan blue staining and compared to the effects of the gamma-secretase inhibitor DBZ. 10 mg/kg anti-Notch1 was used for in vivo tumor growth of YUMM2.1 and 1.7 cells. Tumors were measured and processed for flow cytometry using antibodies against major immune cell populations.

RESULTS

Anti-Notch1 selectively inhibited Notch1 but not Notch2; caused significant melanoma cell death in vitro but did not affect normal melanocytes. In vivo, it delayed tumor growth without evident signs of gastro-intestinal toxicities; and importantly promoted an inflamed TME by increasing the cytotoxic CD8+ T cells while reducing the tolerogenic Tregs and MDSCs, resulting in enhanced efficacy of anti-PD-1.

CONCLUSIONS

Anti-Notch1 safely exerts anti-melanoma effects and improves immune checkpoint inhibitor efficacy. Thus, anti-Notch1 could represent a novel addition to the immunotherapy repertoire for melanoma.

Deletion upstream of MAB21L2 highlights the importance of evolutionarily conserved non-coding sequences for eye development

Anophthalmia, microphthalmia and coloboma (AMC) comprise a spectrum of developmental eye disorders, accounting for approximately 20% of childhood visual impairment. While non-coding regulatory sequences are increasingly recognised as contributing to disease burden, characterising their impact on gene function and phenotype remains challenging. Furthermore, little is known of the nature and extent of their contribution to AMC phenotypes. We report two families with variants in or near MAB21L2, a gene where genetic variants are known to cause AMC in humans and animal models. The first proband, presenting with microphthalmia and coloboma, has a likely pathogenic missense variant (c.338 G > C; p.[Trp113Ser]), segregating within the family. The second individual, presenting with microphthalmia, carries an ~ 113.5 kb homozygous deletion 19.38 kb upstream of MAB21L2. Modelling of the deletion results in transient small lens and coloboma as well as midbrain anomalies in zebrafish, and microphthalmia and coloboma in Xenopus tropicalis. Using conservation analysis, we identify 15 non-coding conserved elements (CEs) within the deleted region, while ChIP-seq data from mouse embryonic stem cells demonstrates that two of these (CE13 and 14) bind Otx2, a protein with an established role in eye development. Targeted disruption of CE14 in Xenopus tropicalis recapitulates an ocular coloboma phenotype, supporting its role in eye development. Together, our data provides insights into regulatory mechanisms underlying eye development and highlights the importance of non-coding sequences as a source of genetic diagnoses in AMC.

Notch1 signaling pathway promotes growth and metastasis of gastric cancer via modulating CDH5

OBJECTIVE

To explore the underlying molecular mechanism of Notch1/cadherin 5 (CDH5) pathway in modulating in cell malignant behaviors of gastric cancer (GC).

METHODS

We performed bioinformatic analyses to screen the potential target genes of Notch1 from cadherins in GC. Western blot and RT-PCR were conducted to detect CDH5 expression in GC tissues and cells. We utilized chromatin immunoprecipitation (CHIP) assays to assess the interaction of Notch1 with CDH5 gene. The effects of Notch1/CDH5 axis on the proliferation, invasion, migration and vasculogenic mimicry in GC cells were evaluated by EdU, wound healing, transwell, and tubule formation assays.

RESULTS

Significantly increased CDH5 expression was found in GC tissues compared with paracancerous tissues and associated to clinical stage and poor overall survival (OS) in patients with GC. Notch1 positively regulate the expression of CDH5 in GC cells. CHIP assays validated that CDH5 was a direct target of Notch1. In addition, Notch1 upregulation enhanced the proliferation, migration, invasion and vasculogenic mimicry capacity of GC cells, which could be attenuated by CDH5 silencing.

CONCLUSIONS

These results indicated Notch1 upregulation enhanced GC malignant behaviors by triggering CDH5, suggesting that targeting Notch1/CDH5 axis could be a potential therapeutic strategy for GC progression.

Regulation of the Lewis Blood Group Antigen Expression: A Literature Review Supplemented with Computational Analysis

Background

The Lewis (Le) blood group system, unlike most other blood groups, is not defined by antigens produced internally to the erythrocytes and their precursors but rather by glycan antigens adsorbed on to the erythrocyte membrane from the plasma. These oligosaccharides are synthesized by the two fucosyltransferases FUT2 and FUT3 mainly in epithelial cells of the digestive tract and transferred to the plasma. At their place of synthesis, some Lewis blood group carbohydrate antigen variants also seem to be involved in various gastrointestinal malignancies. However, relatively little is known about the transcriptional regulation of FUT2 and FUT3.

Summary

To address this question, we screened existing literature and additionally used in silico prediction tools to identify novel candidate regulators for FUT2 and FUT3 and combine these findings with already known data on their regulation. With this approach, we were able to describe a variety of transcription factors, RNA binding proteins and microRNAs, which increase FUT2 and FUT3 transcription and translation upon interaction.

Key Messages

Understanding the regulation of FUT2 and FUT3 is crucial to fully understand the blood group system Lewis (ISBT 007 LE) phenotypes, to shed light on the role of the different Lewis antigens in various pathologies, and to identify potential new diagnostic targets for these diseases.

Identifying cell type-specific transcription factor-mediated activity immune modules reveal implications for immunotherapy and molecular classification of pan-cancer

Systematic investigation of tumor-infiltrating immune (TII) cells is important to the development of immunotherapies, and the clinical response prediction in cancers. There exists complex transcriptional regulation within TII cells, and different immune cell types display specific regulation patterns. To dissect transcriptional regulation in TII cells, we first integrated the gene expression profiles from single-cell datasets, and proposed a computational pipeline to identify TII cell type-specific transcription factor (TF) mediated activity immune modules (TF-AIMs). Our analysis revealed key TFs, such as BACH2 and NFKB1 play important roles in B and NK cells, respectively. We also found some of these TF-AIMs may contribute to tumor pathogenesis. Based on TII cell type-specific TF-AIMs, we identified eight CD8+ T cell subtypes. In particular, we found the PD1 + CD8+ T cell subset and its specific TF-AIMs associated with immunotherapy response. Furthermore, the TII cell type-specific TF-AIMs displayed the potential to be used as predictive markers for immunotherapy response of cancer patients. At the pan-cancer level, we also identified and characterized six molecular subtypes across 9680 samples based on the activation status of TII cell type-specific TF-AIMs. Finally, we constructed a user-friendly web interface CellTF-AIMs (http://bio-bigdata.hrbmu.edu.cn/CellTF-AIMs/) for exploring transcriptional regulatory pattern in various TII cell types. Our study provides valuable implications and a rich resource for understanding the mechanisms involved in cancer microenvironment and immunotherapy.

Prostate cancer detection through unbiased capture of methylated cell-free DNA

Prostate cancer screening using prostate-specific antigen (PSA) has been shown to reduce mortality but with substantial overdiagnosis, leading to unnecessary biopsies. The identification of a highly specific biomarker using liquid biopsies, represents an unmet need in the diagnostic pathway for prostate cancer. In this study, we employed a method that enriches for methylated cell-free DNA fragments coupled with a machine learning algorithm which enabled the detection of metastatic and localized cancers with AUCs of 0.96 and 0.74, respectively. The model also detected 51.8% (14/27) of localized and 88.7% (79/89) of patients with metastatic cancer in an external dataset. Furthermore, we show that the differentially methylated regions reflect epigenetic and transcriptomic changes at the tissue level. Notably, these regions are significantly enriched for biologically relevant pathways associated with the regulation of cellular proliferation and TGF-beta signaling. This demonstrates the potential of circulating tumor DNA methylation for prostate cancer detection and prognostication.

Identification of the New GmJAG1 Transcription Factor Binding Motifs Using DAP-Seq

Interaction between transcription factors (TFs) and motifs is essential for gene regulation and the subsequent phenotype formation. Soybean (Glycine max) JAGGEED 1 (GmJAG1) is a key TF that controls leaf shape, seed number and flower size. To understand the GmJAG1 binding motifs, in this study, we performed the GmJAG1 DNA affinity purification sequencing (DAP-seq) experiment, which is a powerful tool for the de novo motif prediction method. Two new significant GmJAG1 binding motifs were predicted and the EMSA experiments further verified the ability of GmJAG1 bound to these motifs. The potential binding sites in the downstream gene promoter were identified through motif scanning and a potential regulatory network mediated by GmJAG1 was constructed. These results served as important genomic resources for further understanding the regulatory mechanism of GmJAG1.

Metabolic regulator ERRγ governs gastric stem cell differentiation into acid-secreting parietal cells

Parietal cells (PCs) produce gastric acid to kill pathogens and aid digestion. Dysregulated PC census is common in disease, yet how PCs differentiate is unclear. Here, we identify the PC progenitors arising from isthmal stem cells, using mouse models and human gastric cells, and show that they preferentially express cell-metabolism regulator and orphan nuclear receptor Estrogen-related receptor gamma (Esrrg, encoding ERRγ). Esrrg expression facilitated the tracking of stepwise molecular, cellular, and ultrastructural stages of PC differentiation. EsrrgP2ACreERT2 lineage tracing revealed that Esrrg expression commits progenitors to differentiate into mature PCs. scRNA-seq indicated the earliest Esrrg+ PC progenitors preferentially express SMAD4 and SP1 transcriptional targets and the GTPases regulating acid-secretion signal transduction. As progenitors matured, ERRγ-dependent metabolic transcripts predominated. Organoid and mouse studies validated the requirement of ERRγ for PC differentiation. Our work chronicles stem cell differentiation along a single lineage in vivo and suggests ERRγ as a therapeutic target for PC-related disorders.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

CACIMAR: cross-species analysis of cell identities, markers, regulations, and interactions using single-cell RNA sequencing data

Transcriptomic analysis across species is increasingly used to reveal conserved gene regulations which implicate crucial regulators. Cross-species analysis of single-cell RNA sequencing (scRNA-seq) data provides new opportunities to identify the cellular and molecular conservations, especially for cell types and cell type-specific gene regulations. However, few methods have been developed to analyze cross-species scRNA-seq data to uncover both molecular and cellular conservations. Here, we built a tool called CACIMAR, which can perform cross-species analysis of cell identities, markers, regulations, and interactions using scRNA-seq profiles. Based on the weighted sum models of the conserved features, we developed different conservation scores to measure the conservation of cell types, regulatory networks, and intercellular interactions. Using publicly available scRNA-seq data on retinal regeneration in mice, zebrafish, and chick, we demonstrated four main functions of CACIMAR. First, CACIMAR allows to identify conserved cell types even in evolutionarily distant species. Second, the tool facilitates the identification of evolutionarily conserved or species-specific marker genes. Third, CACIMAR enables the identification of conserved intracellular regulations, including cell type-specific regulatory subnetworks and regulators. Lastly, CACIMAR provides a unique feature for identifying conserved intercellular interactions. Overall, CACIMAR facilitates the identification of evolutionarily conserved cell types, marker genes, intracellular regulations, and intercellular interactions, providing insights into the cellular and molecular mechanisms of species evolution.

Gene expression supports a single origin of horns and antlers in hoofed mammals

Horns, antlers, and other bony cranial appendages of even-toed hoofed mammals (ruminant artiodactyls) challenge traditional morphological homology assessments. Cranial appendages all share a permanent bone portion with family-specific integument coverings, but homology determination depends on whether the integument covering is an essential component or a secondary elaboration of each structure. To enhance morphological homology assessments, we tested whether juvenile cattle horn bud transcriptomes share homologous gene expression patterns with deer antlers relative to pig outgroup tissues, treating the integument covering as a secondary elaboration. We uncovered differentially expressed genes that support horn and antler homology, potentially distinguish them from non-cranial-appendage bone and other tissues, and highlight the importance of phylogenetic outgroups in homology assessments. Furthermore, we found differentially expressed genes that could support a shared cranial neural crest origin for horns and antlers and expression patterns that refine our understanding of the timing of horn and antler differentiation.

Exploring the regulatory interactions between mutated genes and homeobox genes in the head and neck cancer progression

OBJECTIVE

Understanding the regulatory role of homeobox (HOX) and mutated genes in the progression of head and neck cancers is essential, although their interaction remains elusive. This study aims to decipher the critical regulation of mutation driven effects on homeobox genes to enhance our understanding of head and neck cancer progression.

METHODS

Genomic mutation data from The Cancer Genome Atlas-Head and Neck Squamous Cell Carcinoma were analyzed using VarScan2 for somatic variant detection. Mutational clustering, driver mutation identification, and cancer signaling pathway analysis were performed using the OncodriveCLUST method. Harmonizome datasets were retrieved to identify critical cancer driver genes affecting HOX genes. The effects of HPV infection on HOX and mutated genes were assessed using the oncoviral database. Altered pathway activity due to the effects of cancer drivers on HOX genes was analyzed with Gene Set Cancer Analysis. Functional enrichment analysis of gene ontology biological processes and molecular functions was conducted using the ClusterProfiler R package.

RESULTS

Significant alterations in HOX genes were observed in head and neck cancer cohorts with mutated TP53, FAT1, and CDKN2A. HOX genes were identified as functionally downstream targets of TP53, signifying transcriptionally mediated regulation. The interaction between HOX genes and mutated TP53, FAT1, and CDKN2A dysregulated the epithelial-to-mesenchymal transition, cell cycle, and apoptosis pathways in head and neck cancer progression.

CONCLUSION

The interplay between cancer driver genes and HOX genes is pivotal in regulating the oncogenic processes underlying the pathogenesis of head and neck squamous cell carcinoma.

Cis-regulatory control of mammalian Trps1 gene expression

TRPS1 serves as the causative gene for tricho-rhino phalangeal syndrome, known for its craniofacial and skeletal abnormalities. The Trps1 gene encodes a protein that represses Wnt signaling through strong interactions with Wnt signaling inhibitors. The identification of genomic cis-acting regulatory sequences governing Trps1 expression is crucial for understanding its role in embryogenesis. Nevertheless, to date, no investigations have been conducted concerning these aspects of Trps1. To identify deeply conserved noncoding elements (CNEs) within the Trps1 locus, we employed a comparative genomics approach, utilizing slowly evolving fish such as coelacanth and spotted gar. These analyses resulted in the identification of eight CNEs in the intronic region of the Trps1 gene. Functional characterization of these CNEs in zebrafish revealed their regulatory potential in various tissues, including pectoral fins, heart, and pharyngeal arches. RNA in-situ hybridization experiments revealed concordance between the reporter expression pattern induced by the identified set of CNEs and the spatial expression pattern of the trps1 gene in zebrafish. Comparative in vivo data from zebrafish and mice for CNE7/hs919 revealed conserved functions of these enhancers. Each of these eight CNEs was further investigated in cell line-based reporter assays, revealing their repressive potential. Taken together, in vivo and in vitro assays suggest a context-dependent dual functionality for the identified set of Trps1-associated CNE enhancers. This functionally characterized set of CNE-enhancers will contribute to a more comprehensive understanding of the developmental roles of Trps1 and can aid in the identification of noncoding DNA variants associated with human diseases.

Gene expression analysis reveals diabetes-related gene signatures

BACKGROUND

Diabetes is a spectrum of metabolic diseases affecting millions of people worldwide. The loss of pancreatic β-cell mass by either autoimmune destruction or apoptosis, in type 1-diabetes (T1D) and type 2-diabetes (T2D), respectively, represents a pathophysiological process leading to insulin deficiency. Therefore, therapeutic strategies focusing on restoring β-cell mass and β-cell insulin secretory capacity may impact disease management. This study took advantage of powerful integrative bioinformatic tools to scrutinize publicly available diabetes-associated gene expression data to unveil novel potential molecular targets associated with β-cell dysfunction.

METHODS

A comprehensive literature search for human studies on gene expression alterations in the pancreas associated with T1D and T2D was performed. A total of 6 studies were selected for data extraction and for bioinformatic analysis. Pathway enrichment analyses of differentially expressed genes (DEGs) were conducted, together with protein-protein interaction networks and the identification of potential transcription factors (TFs). For noncoding differentially expressed RNAs, microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), which exert regulatory activities associated with diabetes, identifying target genes and pathways regulated by these RNAs is fundamental for establishing a robust regulatory network.

RESULTS

Comparisons of DEGs among the 6 studies showed 59 genes in common among 4 or more studies. Besides alterations in mRNA, it was possible to identify differentially expressed miRNA and lncRNA. Among the top transcription factors (TFs), HIPK2, KLF5, STAT1 and STAT3 emerged as potential regulators of the altered gene expression. Integrated analysis of protein-coding genes, miRNAs, and lncRNAs pointed out several pathways involved in metabolism, cell signaling, the immune system, cell adhesion, and interactions. Interestingly, the GABAergic synapse pathway emerged as the only common pathway to all datasets.

CONCLUSIONS

This study demonstrated the power of bioinformatics tools in scrutinizing publicly available gene expression data, thereby revealing potential therapeutic targets like the GABAergic synapse pathway, which holds promise in modulating α-cells transdifferentiation into β-cells.

A meta-analysis of Boolean network models reveals design principles of gene regulatory networks

Gene regulatory networks (GRNs) play a central role in cellular decision-making. Understanding their structure and how it impacts their dynamics constitutes thus a fundamental biological question. GRNs are frequently modeled as Boolean networks, which are intuitive, simple to describe, and can yield qualitative results even when data are sparse. We assembled the largest repository of expert-curated Boolean GRN models. A meta-analysis of this diverse set of models reveals several design principles. GRNs exhibit more canalization, redundancy, and stable dynamics than expected. Moreover, they are enriched for certain recurring network motifs. This raises the important question why evolution favors these design mechanisms.

Targeting class A GPCRs for hard tissue regeneration

G protein-coupled receptors (GPCRs) play important roles in various pathogeneses and physiological regulations. Owing to their functional diversity, GPCRs are considered one of the primary pharmaceutical targets. However, drugs targeting GPCRs have not been developed yet to regenerate hard tissues such as teeth and bones. Mesenchymal stromal cells (MSCs) have high proliferation and multi-lineage differentiation potential, which are essential for hard tissue regeneration. Here, we present a strategy for targeting class A GPCRs for hard tissue regeneration by promoting the differentiation of endogenous MSCs into osteogenic and odontogenic progenitor cells. Through in vitro screening targeted at class A GPCRs, we identified six target receptors (LPAR1, F2R, F2RL1, F2RL2, S1PR1, and ADORA2A) and candidate drugs with potent biomineralization effects. Through a combination of profiling whole transcriptome and accessible chromatin regions, we identified that p53 acts as a key transcriptional activator of genes that modulate the biomineralization process. Moreover, the therapeutic potential of class A GPCR-targeting drugs was demonstrated in tooth pulpotomy and calvarial defect models. The selected drugs revealed potent regenerative effects in both tooth and bone defects, represented by newly formed highly mineralized regions. Consequently, this study provides translational evidence for a new regenerative strategy for damaged hard tissue.

Design principles for synthetic control systems to engineer plants

KEY MESSAGE

Synthetic control systems have led to significant advancement in the study and engineering of unicellular organisms, but it has been challenging to apply these tools to multicellular organisms like plants. The ability to predictably engineer plants will enable the development of novel traits capable of alleviating global problems, such as climate change and food insecurity. Engineering predictable multicellular phenotypes will require the development of synthetic control systems that can precisely regulate how the information encoded in genomes is translated into phenotypes. Many efficient control systems have been developed for unicellular organisms. However, it remains challenging to use such tools to study or engineer multicellular organisms. Plants are a good chassis within which to develop strategies to overcome these challenges, thanks to their capacity to withstand large-scale reprogramming without lethality. Additionally, engineered plants have great potential for solving major societal problems. Here we briefly review the progress of control system development in unicellular organisms, and how that information can be leveraged to characterize control systems in plants. Further, we discuss strategies for developing control systems designed to regulate the expression of transgenes or endogenous loci and generate dosage-dependent or discrete traits. Finally, we discuss the utility that mathematical models of biological processes have for control system deployment.

Biocaiv: an integrative webserver for motif-based clustering analysis and interactive visualization of biological networks

BACKGROUND

As an important task in bioinformatics, clustering analysis plays a critical role in understanding the functional mechanisms of many complex biological systems, which can be modeled as biological networks. The purpose of clustering analysis in biological networks is to identify functional modules of interest, but there is a lack of online clustering tools that visualize biological networks and provide in-depth biological analysis for discovered clusters.

RESULTS

Here we present BioCAIV, a novel webserver dedicated to maximize its accessibility and applicability on the clustering analysis of biological networks. This, together with its user-friendly interface, assists biological researchers to perform an accurate clustering analysis for biological networks and identify functionally significant modules for further assessment.

CONCLUSIONS

BioCAIV is an efficient clustering analysis webserver designed for a variety of biological networks. BioCAIV is freely available without registration requirements at http://bioinformatics.tianshanzw.cn:8888/BioCAIV/ .

Pooled genetic screens to identify vulnerabilities in TERT-promoter-mutant glioblastoma

Pooled genetic screens represent a powerful approach to identify vulnerabilities in cancer. Here we used pooled CRISPR/Cas9-based approaches to identify vulnerabilities associated with telomerase reverse transcriptase (TERT) promoter mutations (TPMs) found in >80% of glioblastomas. We first developed a platform to detect perturbations that cause long-term growth defects in a TPM-mutated glioblastoma cell line. However, we could not detect dependencies on either TERT itself or on an E-twenty six transcription (ETS) factor known to activate TPMs. To explore this finding, we cataloged TPM status for 441 cell lines and correlated this with genome-wide screening data. We found that TPM status was not associated with differential dependency on TERT, but that E-twenty six (ETS) transcription factors represent key dependencies in both TPM+ and TPM- lines. Further, we found that TPMs are associated with expression of gene programs regulated by a wide array of ETS-factors in both cell lines and primary glioblastoma tissues. This work contributes a unique TPM cell line reagent, establishes TPM status for many deeply-profiled cell lines, and catalogs TPM-associated vulnerabilities. The results highlight challenges in executing genetic screens to detect TPM-specific vulnerabilities, and suggest redundancy in the genetic network that regulates TPM function with therapeutic implications.

Towards a comprehensive regulatory map of Mammalian Genomes

Genome mapping studies have generated a nearly complete collection of genes for the human genome, but we still lack an equivalently vetted inventory of human regulatory sequences. Cis-regulatory modules (CRMs) play important roles in controlling when, where, and how much a gene is expressed. We developed a training data-free CRM-prediction algorithm, the Mammalian Regulatory MOdule Detector (MrMOD) for accurate CRM prediction in mammalian genomes. MrMOD provides genome position-fixed CRM models similar to the fixed gene models for the mouse and human genomes using only genomic sequences as the inputs with one adjustable parameter - the significance p-value. Importantly, MrMOD predicts a comprehensive set of high-resolution CRMs in the mouse and human genomes including all types of regulatory modules not limited to any tissue, cell type, developmental stage, or condition. We computationally validated MrMOD predictions used a compendium of 21 orthogonal experimental data sets including thousands of experimentally defined CRMs and millions of putative regulatory elements derived from hundreds of different tissues, cell types, and stimulus conditions obtained from multiple databases. In ovo transgenic reporter assay demonstrates the power of our prediction in guiding experimental design. We analyzed CRMs located in the chromosome 17 using unsupervised machine learning and identified groups of CRMs with multiple lines of evidence supporting their functionality, linking CRMs with upstream binding transcription factors and downstream target genes. Our work provides a comprehensive base pair resolution annotation of the functional regulatory elements and non-functional regions in the mammalian genomes.

Plasma cortisol-linked gene networks in hepatic and adipose tissues implicate corticosteroid-binding globulin in modulating tissue glucocorticoid action and cardiovascular risk

Genome-wide association meta-analysis (GWAMA) by the Cortisol Network (CORNET) consortium identified genetic variants spanning the SERPINA6/SERPINA1 locus on chromosome 14 associated with morning plasma cortisol, cardiovascular disease (CVD), and SERPINA6 mRNA expression encoding corticosteroid-binding globulin (CBG) in the liver. These and other findings indicate that higher plasma cortisol levels are causally associated with CVD; however, the mechanisms by which variations in CBG lead to CVD are undetermined. Using genomic and transcriptomic data from The Stockholm Tartu Atherosclerosis Reverse Networks Engineering Task (STARNET) study, we identified plasma cortisol-linked single-nucleotide polymorphisms (SNPs) that are trans-associated with genes from seven different vascular and metabolic tissues, finding the highest representation of trans-genes in the liver, subcutaneous fat, and visceral abdominal fat, [false discovery rate (FDR) = 15%]. We identified a subset of cortisol-associated trans-genes that are putatively regulated by the glucocorticoid receptor (GR), the primary transcription factor activated by cortisol. Using causal inference, we identified GR-regulated trans-genes that are responsible for the regulation of tissue-specific gene networks. Cis-expression Quantitative Trait Loci (eQTLs) were used as genetic instruments for identification of pairwise causal relationships from which gene networks could be reconstructed. Gene networks were identified in the liver, subcutaneous fat, and visceral abdominal fat, including a high confidence gene network specific to subcutaneous adipose (FDR = 10%) under the regulation of the interferon regulatory transcription factor, IRF2. These data identify a plausible pathway through which variation in the liver CBG production perturbs cortisol-regulated gene networks in peripheral tissues and thereby promote CVD.

Pan-cancer Analysis Reveals Cancer-dependent Expression of SOX17 and Associated Clinical Outcomes

BACKGROUND/AIM

SRY-box containing gene 17 (SOX17) plays a pivotal role in cancer onset and progression and is considered a potential target for cancer diagnosis and treatment. However, the expression pattern of SOX17 in cancer and its clinical relevance remains unknown. Here, we explored the relationship between the expression of SOX17 and drug response by examining SOX17 expression patterns across multiple cancer types.

MATERIALS AND METHODS

Single-cell and bulk RNA-seq analyses were used to explore the expression profile of SOX17. Analysis results were verified with qPCR and immunohistochemistry. Survival, drug response, and co-expression analyses were performed to illustrate its correlation with clinical outcomes.

RESULTS

The results revealed that abnormal expression of SOX17 is highly heterogenous across multiple cancer types, indicating that SOX17 manifests as a cancer type-dependent feature. Furthermore, the expression pattern of SOX17 is also associated with cancer prognosis in certain cancer types. Strong SOX17 expression correlates with the potency of small molecule drugs that affect PI3K/mTOR signaling. FGF18, a gene highly relevant to SOX17, is involved in the PI3K-AKT signaling pathway. Single-cell RNA-seq analysis demonstrated that SOX17 is mainly expressed in endothelial cells and barely expressed in other cells but spreads to other cell types during the development of ovarian cancer.

CONCLUSION

Our study revealed the expression pattern of SOX17 in pan-cancer through bulk and single-cell RNA-seq analyses and determined that SOX17 is related to the diagnosis, staging, and prognosis of some tumors. These findings have clinical implications and may help identify mechanistic pathways amenable to pharmacological interventions.

ACE2 in chronic disease and COVID-19: gene regulation and post-translational modification

Angiotensin-converting enzyme 2 (ACE2), a counter regulator of the renin-angiotensin system, provides protection against several chronic diseases. Besides chronic diseases, ACE2 is the host receptor for SARS-CoV or SARS-CoV-2 virus, mediating the first step of virus infection. ACE2 levels are regulated by transcriptional, post-transcriptional, and post-translational regulation or modification. ACE2 transcription is enhanced by transcription factors including Ikaros, HNFs, GATA6, STAT3 or SIRT1, whereas ACE2 transcription is reduced by the transcription factor Brg1-FoxM1 complex or ERRα. ACE2 levels are also regulated by histone modification or miRNA-induced destabilization. The protein kinase AMPK, CK1α, or MAP4K3 phosphorylates ACE2 protein and induces ACE2 protein levels by decreasing its ubiquitination. The ubiquitination of ACE2 is induced by the E3 ubiquitin ligase MDM2 or UBR4 and decreased by the deubiquitinase UCHL1 or USP50. ACE2 protein levels are also increased by the E3 ligase PIAS4-mediated SUMOylation or the methyltransferase PRMT5-mediated ACE2 methylation, whereas ACE2 protein levels are decreased by AP2-mediated lysosomal degradation. ACE2 is downregulated in several human chronic diseases like diabetes, hypertension, or lung injury. In contrast, SARS-CoV-2 upregulates ACE2 levels, enhancing host cell susceptibility to virus infection. Moreover, soluble ACE2 protein and exosomal ACE2 protein facilitate SARS-CoV-2 infection into host cells. In this review, we summarize the gene regulation and post-translational modification of ACE2 in chronic disease and COVID-19. Understanding the regulation and modification of ACE2 may help to develop prevention or treatment strategies for ACE2-mediated diseases.

Flame (v2.0): advanced integration and interpretation of functional enrichment results from multiple sources

Functional enrichment is the process of identifying implicated functional terms from a given input list of genes or proteins. In this article, we present Flame (v2.0), a web tool which offers a combinatorial approach through merging and visualizing results from widely used functional enrichment applications while also allowing various flexible input options. In this version, Flame utilizes the aGOtool, g: Profiler, WebGestalt, and Enrichr pipelines and presents their outputs separately or in combination following a visual analytics approach. For intuitive representations and easier interpretation, it uses interactive plots such as parameterizable networks, heatmaps, barcharts, and scatter plots. Users can also: (i) handle multiple protein/gene lists and analyse union and intersection sets simultaneously through interactive UpSet plots, (ii) automatically extract genes and proteins from free text through text-mining and Named Entity Recognition (NER) techniques, (iii) upload single nucleotide polymorphisms (SNPs) and extract their relative genes, or (iv) analyse multiple lists of differentially expressed proteins/genes after selecting them interactively from a parameterizable volcano plot. Compared to the previous version of 197 supported organisms, Flame (v2.0) currently allows enrichment for 14 436 organisms.

AVAILABILITY AND IMPLEMENTATION

Web Application: http://flame.pavlopouloslab.info. Code: https://github.com/PavlopoulosLab/Flame. Docker: https://hub.docker.com/r/pavlopouloslab/flame.

Proteomics analysis of plasma from middle-aged adults identifies protein markers of dementia risk in later life

A diverse set of biological processes have been implicated in the pathophysiology of Alzheimer's disease (AD) and related dementias. However, there is limited understanding of the peripheral biological mechanisms relevant in the earliest phases of the disease. Here, we used a large-scale proteomics platform to examine the association of 4877 plasma proteins with 25-year dementia risk in 10,981 middle-aged adults. We found 32 dementia-associated plasma proteins that were involved in proteostasis, immunity, synaptic function, and extracellular matrix organization. We then replicated the association between 15 of these proteins and clinically relevant neurocognitive outcomes in two independent cohorts. We demonstrated that 12 of these 32 dementia-associated proteins were associated with cerebrospinal fluid (CSF) biomarkers of AD, neurodegeneration, or neuroinflammation. We found that eight of these candidate protein markers were abnormally expressed in human postmortem brain tissue from patients with AD, although some of the proteins that were most strongly associated with dementia risk, such as GDF15, were not detected in these brain tissue samples. Using network analyses, we found a protein signature for dementia risk that was characterized by dysregulation of specific immune and proteostasis/autophagy pathways in adults in midlife ~20 years before dementia onset, as well as abnormal coagulation and complement signaling ~10 years before dementia onset. Bidirectional two-sample Mendelian randomization genetically validated nine of our candidate proteins as markers of AD in midlife and inferred causality of SERPINA3 in AD pathogenesis. Last, we prioritized a set of candidate markers for AD and dementia risk prediction in midlife.

Signaling mechanisms in renal compensatory hypertrophy revealed by multi-omics

Loss of a kidney results in compensatory growth of the remaining kidney, a phenomenon of considerable clinical importance. However, the mechanisms involved are largely unknown. Here, we use a multi-omic approach in a unilateral nephrectomy model in male mice to identify signaling processes associated with renal compensatory hypertrophy, demonstrating that the lipid-activated transcription factor peroxisome proliferator-activated receptor alpha (PPARα) is an important determinant of proximal tubule cell size and is a likely mediator of compensatory proximal tubule hypertrophy.

Erythroid lineage chromatin accessibility maps facilitate identification and validation of NFIX as a fetal hemoglobin repressor

Human genetics has validated de-repression of fetal gamma globin (HBG) in adult erythroblasts as a powerful therapeutic paradigm in diseases involving defective adult beta globin (HBB)¹. To identify factors involved in the switch from HBG to HBB expression, we performed Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq)² on sorted erythroid lineage cells derived from bone marrow (BM) or cord blood (CB), representing adult and fetal states, respectively. BM to CB cell ATAC-seq profile comparisons revealed genome-wide enrichment of NFI DNA binding motifs and increased NFIX promoter chromatin accessibility, suggesting that NFIX may repress HBG. NFIX knockdown in BM cells increased HBG mRNA and fetal hemoglobin (HbF) protein levels, coincident with increased chromatin accessibility and decreased DNA methylation at the HBG promoter. Conversely, overexpression of NFIX in CB cells reduced HbF levels. Identification and validation of NFIX as a new target for HbF activation has implications in the development of therapeutics for hemoglobinopathies.

Potential biomarkers in Japanese encephalitis from different hosts and geographical locations

Japanese encephalitis (JE) is a single-stranded, mosquito-borne, positive-sense RNA flavivirus that causes one of the most severe encephalitides. There are treatments available for those who contact this illness; however, there are no known cures. This disease has a 30% fatality rate, and of the people who survive, 30-50% develops neurologic and psychiatric sequelae. The JE virus genome size is 10.98 kb and contains two coding DNA sequences (CDS), two genes, and 15 mature peptides; the CDS polyprotein is 10.3 kb. In this study, we used 29 genomics sequences of the JE virus reported from different countries and infecting different animals and analysed vast dimensions of the genomic annotation of JE comparatively to understand its evolutionary aspects. The extensive SNPs analysis revealed that KF907505.1, reported from Taiwan, has only three SNPs, similar to sequences reported from India. Repeat and polymorphism analyses revealed that the genome tends to be similar in most JE sequences.

Positive selection of somatically mutated clones identifies adaptive pathways in metabolic liver disease

Somatic mutations in nonmalignant tissues accumulate with age and injury, but whether these mutations are adaptive on the cellular or organismal levels is unclear. To interrogate genes in human metabolic disease, we performed lineage tracing in mice harboring somatic mosaicism subjected to nonalcoholic steatohepatitis (NASH). Proof-of-concept studies with mosaic loss of Mboat7, a membrane lipid acyltransferase, showed that increased steatosis accelerated clonal disappearance. Next, we induced pooled mosaicism in 63 known NASH genes, allowing us to trace mutant clones side by side. This in vivo tracing platform, which we coined MOSAICS, selected for mutations that ameliorate lipotoxicity, including mutant genes identified in human NASH. To prioritize new genes, additional screening of 472 candidates identified 23 somatic perturbations that promoted clonal expansion. In validation studies, liver-wide deletion of Tbx3, Bcl6, or Smyd2 resulted in protection against hepatic steatosis. Selection for clonal fitness in mouse and human livers identifies pathways that regulate metabolic disease.

Identification and characterization of the genome of a papillomavirus from skin lesions of four-toed hedgehogs (Atelerix albiventris)

The present study describes the clinical and pathological characteristics of skin lesions in two four-toed hedgehogs (Atelerix albiventris). We performed inverse PCR to identify the genome of papillomavirus (PV) in the skin lesions and subsequently sequenced the full genome of the virus, which was tentatively named Atelerix albiventris papillomavirus 1 (AalbPV1). The overall sequences of the viral genomes of both four-toed hedgehogs were identical. This study first identified the presence of a novel PV in Japanese four-toed hedgehogs and provided genetic information about this virus.

RegVar: Tissue-specific Prioritization of Non-coding Regulatory Variants

Non-coding genomic variants constitute the majority of trait-associated genome variations; however, the identification of functional non-coding variants is still a challenge in human genetics, and a method for systematically assessing the impact of regulatory variants on gene expression and linking these regulatory variants to potential target genes is still lacking. Here, we introduce a deep neural network (DNN)-based computational framework, RegVar, which can accurately predict the tissue-specific impact of non-coding regulatory variants on target genes. We show that by robustly learning the genomic characteristics of massive variant-gene expression associations in a variety of human tissues, RegVar vastly surpasses all current non-coding variant prioritization methods in predicting regulatory variants under different circumstances. The unique features of RegVar make it an excellent framework for assessing the regulatory impact of any variant on its putative target genes in a variety of tissues. RegVar is available as a web server at https://regvar.omic.tech/.

Positive selection of somatically mutated clones identifies adaptive pathways in metabolic liver disease

Somatic mutations in non-malignant tissues accumulate with age and insult, but whether these mutations are adaptive on the cellular or organismal levels is unclear. To interrogate mutations found in human metabolic disease, we performed lineage tracing in mice harboring somatic mosaicism subjected to non-alcoholic steatohepatitis (NASH). Proof-of-concept studies with mosaic loss of Mboat7 , a membrane lipid acyltransferase, showed that increased steatosis accelerated clonal disappearance. Next, we induced pooled mosaicism in 63 known NASH genes, allowing us to trace mutant clones side-by-side. This in vivo tracing platform, which we coined MOSAICS, selected for mutations that ameliorate lipotoxicity, including mutant genes identified in human NASH. To prioritize new genes, additional screening of 472 candidates identified 23 somatic perturbations that promoted clonal expansion. In validation studies, liver-wide deletion of Bcl6, Tbx3, or Smyd2 resulted in protection against NASH. Selection for clonal fitness in mouse and human livers identifies pathways that regulate metabolic disease.

Highlights

Mosaic Mboat7 mutations that increase lipotoxicity lead to clonal disappearance in NASH. In vivo screening can identify genes that alter hepatocyte fitness in NASH. Mosaic Gpam mutations are positively selected due to reduced lipogenesis. In vivo screening of transcription factors and epifactors identified new therapeutic targets in NASH.

A Database of Lung Cancer-Related Genes for the Identification of Subtype-Specific Prognostic Biomarkers

The molecular subtype is critical for accurate treatment and follow-up in patients with lung cancer; however, information regarding subtype-associated genes is dispersed among thousands of published studies. Systematic curation and cross-validation of the scientific literature would provide a solid foundation for comparative genetic studies of the major molecular subtypes of lung cancer. Here, we constructed a literature-based lung cancer gene database (LCGene). In the current release, we collected and curated 2507 unique human genes, including 2267 protein-coding and 240 non-coding genes from comprehensive manual examination of 10,960 PubMed article abstracts. Extensive annotations were added to aid identification of differentially expressed genes, potential gene editing sites, and non-coding gene regulation. For instance, we prepared 607 curated genes with CRISPR knockout information in 43 lung cancer cell lines. Further comparison of these implicated genes among different subtypes identified several subtype-specific genes with high mutational frequencies. Common tumor suppressors and oncogenes shared by lung adenocarcinoma and lung squamous cell carcinoma, for example, exhibited different mutational frequencies and prognostic features, suggesting the presence of subtype-specific biomarkers. Our retrospective analysis revealed 43 small cell lung cancer-specific genes. Moreover, 52 tumor suppressors and oncogenes shared by lung adenocarcinoma and squamous cell carcinoma confirmed the different molecular mechanisms of these two cancer subtypes. The subtype-based genetic differences, when combined, may provide insight into subtype-specific biomarkers for genetic testing.

Investigation of autism-related transcription factors underlying sex differences in the effects of bisphenol A on transcriptome profiles and synaptogenesis in the offspring hippocampus

BACKGROUND

Bisphenol A (BPA) has been linked to susceptibility to autism spectrum disorder (ASD). Our recent studies have shown that prenatal BPA exposure disrupted ASD-related gene expression in the hippocampus, neurological functions, and behaviors associated with ASD in a sex-specific pattern. However, the molecular mechanisms underlying the effects of BPA are still unclear.

METHODS

Transcriptome data mining and molecular docking analyses were performed to identify ASD-related transcription factors (TFs) and their target genes underlying the sex-specific effects of prenatal BPA exposure. Gene ontology analysis was conducted to predict biological functions associated with these genes. The expression levels of ASD-related TFs and targets in the hippocampus of rat pups prenatally exposed to BPA were measured using qRT-PCR analysis. The role of the androgen receptor (AR) in BPA-mediated regulation of ASD candidate genes was investigated using a human neuronal cell line stably transfected with AR-expression or control plasmid. Synaptogenesis, which is a function associated with genes transcriptionally regulated by ASD-related TFs, was assessed using primary hippocampal neurons isolated from male and female rat pups prenatally exposed to BPA.

RESULTS

We found that there was a sex difference in ASD-related TFs underlying the effects of prenatal BPA exposure on the transcriptome profiles of the offspring hippocampus. In addition to the known BPA targets AR and ESR1, BPA could directly interact with novel targets (i.e., KDM5B, SMAD4, and TCF7L2). The targets of these TFs were also associated with ASD. Prenatal BPA exposure disrupted the expression of ASD-related TFs and targets in the offspring hippocampus in a sex-dependent manner. Moreover, AR was involved in the BPA-mediated dysregulation of AUTS2, KMT2C, and SMARCC2. Prenatal BPA exposure altered synaptogenesis by increasing synaptic protein levels in males but not in females, but the number of excitatory synapses was increased in female primary neurons only.

CONCLUSIONS

Our findings suggest that AR and other ASD-related TFs are involved in sex differences in the effects of prenatal BPA exposure on transcriptome profiles and synaptogenesis in the offspring hippocampus. These TFs may play an essential role in an increased ASD susceptibility associated with endocrine-disrupting chemicals, particularly BPA, and the male bias of ASD.

Integration of single sample and population analysis for understanding immune evasion mechanisms of lung cancer

A deep understanding of the complex interaction mechanism between the various cellular components in tumor microenvironment (TME) of lung adenocarcinoma (LUAD) is a prerequisite for understanding its drug resistance, recurrence, and metastasis. In this study, we proposed two complementary computational frameworks for integrating multi-source and multi-omics data, namely ImmuCycReg framework (single sample level) and L0Reg framework (population or subtype level), to carry out difference analysis between the normal population and different LUAD subtypes. Then, we aimed to identify the possible immune escape pathways adopted by patients with different LUAD subtypes, resulting in immune deficiency which may occur at different stages of the immune cycle. More importantly, combining the research results of the single sample level and population level can improve the credibility of the regulatory network analysis results. In addition, we also established a prognostic scoring model based on the risk factors identified by Lasso-Cox method to predict survival of LUAD patients. The experimental results showed that our frameworks could reliably identify transcription factor (TF) regulating immune-related genes and could analyze the dominant immune escape pathways adopted by each LUAD subtype or even a single sample. Note that the proposed computational framework may be also applicable to the immune escape mechanism analysis of pan-cancer.

Building a knowledge graph to enable precision medicine

Developing personalized diagnostic strategies and targeted treatments requires a deep understanding of disease biology and the ability to dissect the relationship between molecular and genetic factors and their phenotypic consequences. However, such knowledge is fragmented across publications, non-standardized repositories, and evolving ontologies describing various scales of biological organization between genotypes and clinical phenotypes. Here, we present PrimeKG, a multimodal knowledge graph for precision medicine analyses. PrimeKG integrates 20 high-quality resources to describe 17,080 diseases with 4,050,249 relationships representing ten major biological scales, including disease-associated protein perturbations, biological processes and pathways, anatomical and phenotypic scales, and the entire range of approved drugs with their therapeutic action, considerably expanding previous efforts in disease-rooted knowledge graphs. PrimeKG contains an abundance of 'indications', 'contradictions', and 'off-label use' drug-disease edges that lack in other knowledge graphs and can support AI analyses of how drugs affect disease-associated networks. We supplement PrimeKG's graph structure with language descriptions of clinical guidelines to enable multimodal analyses and provide instructions for continual updates of PrimeKG as new data become available.

Nucleus Pulposus Cells from Calcified Discs Promote the Degradation of the Extracellular Matrix through Upregulation of the GATA3 Expression

OBJECTIVE

Disc calcification is strongly associated with disc degeneration; however, the underlying mechanisms driving its pathogenesis are poorly understood. This study aimed to provide a gene expression profile of nucleus pulposus cells (NPCs) from calcified discs, and clarify the potential mechanism in disc degeneration.

METHODS

Primary NPCs were isolated from calcified and control discs (CAL-NPC and CON-NPC), respectively. The proliferation and extracellular matrix (ECM) metabolism capacities of the cells were evaluated using MTT and Western blotting, respectively. RNA sequencing was used to identify differentially expressed genes (DEGs) in the CAL-NPCs. The biological functions of the DEGs were analyzed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The transcription factor database and Cytoscape software were used to construct the transcription factor-DEGs regulatory network. The role of the verified transcription factor in NPC proliferation and ECM metabolism was also investigated.

RESULTS

The CAL-NPCs exhibited a lower proliferation rate and higher ECM degradation capacity than the CON-NPCs. In total, 375 DEGs were identified in the CAL-NPCs. The GO and KEGG analyses showed that the DEGs were primarily involved in the regulation of ribonuclease activity and NF-kappa B and p53 signaling pathways. GATA-binding protein 3 (GATA3) with the highest verified levels was selected for further studies. Overexpression of GATA3 in the CON-NPCs significantly inhibited their proliferation and promoted their ECM degradation function, while the knockdown of GATA3 in the CAL-NPCs resulted in the opposite phenotypes.

CONCLUSION

This study provided a comprehensive gene expression profile of the NPCs from the calcified discs and supported that GATA3 could be a potential target for reversing calcification-associated disc degeneration.

Transcriptional and post-transcriptional controls for tuning gene expression in plants

Plant biotechnologists seek to modify plants through genetic reprogramming, but our ability to precisely control gene expression in plants is still limited. Here, we review transcription and translation in the model plants Arabidopsis thaliana and Nicotiana benthamiana with an eye toward control points that may be used to predictably modify gene expression. We highlight differences in gene expression requirements between these plants and other species, and discuss the ways in which our understanding of gene expression has been used to engineer plants. This review is intended to serve as a resource for plant scientists looking to achieve precise control over gene expression.

Downregulation of hepatic lncRNA Gm19619 improves gluconeogenesis and lipogenesis following vertical sleeve gastrectomy in mice

Long non-coding RNAs (lncRNAs) are emerging important epigenetic regulators in metabolic processes. Whether they contribute to the metabolic effects of vertical sleeve gastrectomy (VSG), one of the most effective treatments for sustainable weight loss and metabolic improvement, is unknown. Herein, we identify a hepatic lncRNA Gm19619, which is strongly repressed by VSG but highly up-regulated by diet-induced obesity and overnight-fasting in mice. Forced transcription of Gm19619 in the mouse liver significantly promotes hepatic gluconeogenesis with the elevated expression of G6pc and Pck1. In contrast, AAV-CasRx mediated knockdown of Gm19619 in high-fat diet-fed mice significantly improves hepatic glucose and lipid metabolism. Mechanistically, Gm19619 is enriched along genomic regions encoding leptin receptor (Lepr) and transcription factor Foxo1, as revealed in chromatin isolation by RNA purification (ChIRP) assay and is confirmed to modulate their transcription in the mouse liver. In conclusion, Gm19619 may enhance gluconeogenesis and lipid accumulation in the liver.

HGCA2.0: An RNA-Seq Based Webtool for Gene Coexpression Analysis in Homo sapiens

Genes with similar expression patterns in a set of diverse samples may be considered coexpressed. Human Gene Coexpression Analysis 2.0 (HGCA2.0) is a webtool which studies the global coexpression landscape of human genes. The website is based on the hierarchical clustering of 55,431 Homo sapiens genes based on a large-scale coexpression analysis of 3500 GTEx bulk RNA-Seq samples of healthy individuals, which were selected as the best representative samples of each tissue type. HGCA2.0 presents subclades of coexpressed genes to a gene of interest, and performs various built-in gene term enrichment analyses on the coexpressed genes, including gene ontologies, biological pathways, protein families, and diseases, while also being unique in revealing enriched transcription factors driving coexpression. HGCA2.0 has been successful in identifying not only genes with ubiquitous expression patterns, but also tissue-specific genes. Benchmarking showed that HGCA2.0 belongs to the top performing coexpression webtools, as shown by STRING analysis. HGCA2.0 creates working hypotheses for the discovery of gene partners or common biological processes that can be experimentally validated. It offers a simple and intuitive website design and user interface, as well as an API endpoint.

NF-YB family transcription factors in Arabidopsis: Structure, phylogeny, and expression analysis in biotic and abiotic stresses

Nuclear factor-Y (NF-Y) transcription factors (TFs) are conserved heterotrimeric complexes present and widespread across eukaryotes. Three main subunits make up the structural and functional aspect of the NF-Y TFs: NF-YA, NF-YB and NF-YC, which bind to the conserved CCAAT- box of the promoter region of specific genes, while also interacting with each other, thereby forming myriad combinations. The NF-YBs are expressed differentially in various tissues and plant development stages, likely impacting many of the cellular processes constitutively and under stress conditions. In this study, ten members of NF-YB family from Arabidopsis thaliana were identified and expression profiles were mined from microarray data under different biotic and abiotic conditions, revealing key insights into the involvement of this class of proteins in the cellular and biological processes in Arabidopsis. Analysis of cis-acting regulatory elements (CAREs) indicated the presence of abiotic and biotic stress-related transcription factor binding sites (TFBs), shedding light on the multifaceted roles of these TFs. Microarray data analysis inferred distinct patterns of expression in various tissues under differing treatments such as drought, cold and heat stress as well as bacterial, fungal, and viral stress, indicating their likelihood of having an expansive range of regulatory functions under native and stressed conditions; while quantitative real-time PCR (qRT-PCR) based expression analysis revealed that these TFs get real-time-modulated in a stress dependent manner. This study, overall, provides an understanding of the AtNF-YB family of TFs in their regulation and participation in various morphogenetic and defense- related pathways and can provide insights for development of transgenic plants for trait dependent studies.

Sexually concordant and dimorphic transcriptional responses to maternal trichloroethylene and/or N-acetyl cysteine exposure in Wistar rat placental tissue

Numerous Superfund sites are contaminated with the volatile organic chemical trichloroethylene (TCE). In women, exposure to TCE in pregnancy is associated with reduced birth weight. Our previous study reported that TCE exposure in pregnant rats decreased fetal weight and elevated oxidative stress biomarkers in placentae, suggesting placental injury as a potential mechanism of TCE-induced adverse birth outcomes. In this study, we investigated if co-exposure with the antioxidant N-acetylcysteine (NAC) attenuates TCE exposure effects on RNA expression. Timed-pregnant Wistar rats were exposed orally to 480 mg TCE/kg/day on gestation days 6-16. Exposure of 200 mg NAC/kg/day alone or as a pre/co-exposure with TCE occurred on gestation days 5-16 to stimulate antioxidant genes prior to TCE exposure. Tissue was collected on gestation day 16. In male and female placentae, we evaluated TCE- and/or NAC-induced changes to gene expression and pathway enrichment analyses using false discovery rate (FDR) and fold-change criteria. In female placentae, exposure to TCE caused significant differential expression 129 genes while the TCE+NAC altered 125 genes, compared with controls (FDR< 0.05 + fold-change >1). In contrast, in male placentae TCE exposure differentially expressed 9 genes and TCE+NAC differentially expressed 35 genes, compared with controls (FDR< 0.05 + fold-change >1). NAC alone did not significantly alter gene expression in either sex. Differentially expressed genes observed with TCE exposure were enriched in mitochondrial biogenesis and oxidative phosphorylation pathways in females whereas immune system pathways and endoplasmic reticulum stress pathways were differentially expressed in both sexes (FDR<0.05). TCE treatment was differentially enriched for genes regulated by the transcription factors ATF6 (both sexes) and ATF4 (males only), indicating a cellular condition triggered by misfolded proteins during endoplasmic reticulum stress. This study demonstrates novel genes and pathways involved in TCE-induced placental injury and showed antioxidant co-treatment largely did not attenuate TCE exposure effects.

Quantitative Proteomics of Nervous System Regeneration: From Sample Preparation to Functional Data Analyses

Mammals have a limited regenerative capacity, especially of the central nervous system. Consequently, any traumatic injury or neurodegenerative disease results in irreversible damage. An important approach to finding strategies to promote regeneration in mammals has been the study of regenerative organisms like Xenopus, the axolotl, and teleost fish. High-throughput technologies like RNA-Seq and quantitative proteomics are starting to provide valuable insight into the molecular mechanisms that drive nervous system regeneration in these organisms. In this chapter, we present a detailed protocol for performing iTRAQ proteomics that can be applied to the analysis of nervous system samples, using Xenopus laevis as an example. The quantitative proteomics protocol and directions for performing functional enrichment data analyses of gene lists (e.g., differentially abundant proteins from a proteomic study, or any type of high-throughput analysis) are aimed at the general bench biologist and do not require previous programming knowledge.

A pituitary gene network linking vgll3 to regulators of sexual maturation in male Atlantic salmon

Age at maturity is a key life history trait and a significant contributor to life history strategy variation. The maturation process is influenced by genetic and environmental factors, but specific causes of variation in maturation timing remain elusive. In many species, the increase in the regulatory gonadotropin-releasing hormone 1 (GnRH1) marks the onset of puberty. Atlantic salmon, however, lacks the gnrh1 gene, suggesting gnrh3 and/or other regulatory factors are involved in the maturation process. Earlier research in Atlantic salmon has found a strong association between alternative alleles of vgll3 and maturation timing. Recently we reported strong induction of gonadotropin genes (fshb and lhb) in the pituitary of Atlantic salmon homozygous for the early maturation allele (E) of vgll3. The induction of gonadotropins was accompanied by increased expression of their direct upstream regulators, c-jun and sf1 (nr5a1b) but the regulatory connection between vgll3 and these regulators has never been investigated in any organism. In this study, we investigated the potential regulatory connection between vgll3 genotypes and these regulators through a stepwise approach of identifying a gene regulatory network (GRN) containing c-jun and sf1, and transcription factor motif enrichment analysis. We found a GRN containing c-jun with predicted upstream regulators, e2f1, egr1, foxj1 and klf4, to be differentially expressed in the pituitary. Finally, we suggest a vgll3 and Hippo pathway -dependent model for transcriptional regulation of c-jun and sf1 in the pituitary, which may have broader implications across vertebrates.

Gene expression patterns associated with caudal fin shape in the cichlid Lamprologus tigripictilis

Variation in fin shape is one of the most prominent features of morphological diversity among fish. Regulation of fin growth has mainly been studied in zebrafish, and it is not clear whether the molecular mechanisms underlying shape variation are equally diverse or rather conserved across species. In the present study, expression levels of 37 candidate genes were tested for association with fin shape in the cichlid fish Lamprologus tigripictilis. The tested genes included members of a fin shape-associated gene regulatory network identified in a previous study and novel candidates selected within this study. Using both intact and regenerating fin tissue, we tested for expression differences between the elongated and the short regions of the spade-shaped caudal fin and identified 20 genes and transcription factors (including angptl5, cd63, csrp1a, cx43, esco2, gbf1, and rbpj), whose expression patterns were consistent with a role in fin growth. Collated with available gene expression data of two other cichlid species, our study not only highlights several genes that were correlated with fin growth in all three species (e.g., angptl5, cd63, cx43, and mmp9), but also reveals species-specific gene expression and correlation patterns, which indicate considerable divergence in the regulatory mechanisms of fin growth across cichlids.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10750-022-05068-4.

Detection and validation of cis-regulatory motifs in osmotic stress-inducible synthetic gene switches via computational and experimental approaches

Synthetic cis -regulatory modules can improve our understanding of gene regulatory networks. We applied an ensemble approach for de novo cis motif discovery among the promoters of 181 drought inducible differentially expressed soybean (Glycine max L.) genes. A total of 43 cis motifs were identified in promoter regions of all gene sets using the binding site estimation suite of tools (BEST). Comparative analysis of these motifs revealed similarities with known cis -elements found in PLACE database and led to the discovery of cis -regulatory motifs that were not yet implicated in drought response. Compiled with the proposed synthetic promoter design rationale, three synthetic assemblies were constructed by concatenating multiple copies of drought-inducible cis motifs in a specific order with inter-motif spacing using random bases and placed upstream of 35s minimal core promoter. Each synthetic module substituted 35S promoter in pBI121 and pCAMBIA3301 to drive glucuronidase expression in soybean hairy roots and Arabidopsis thaliana L. Chimeric soybean seedlings and 3-week-old transgenic Arabidopsis plants were treated with simulated with different levels of osmotic stress. Histochemical staining of transgenic soybean hairy roots and Arabidopsis displayed drought-inducible GUS activity of synthetic promoters. Fluorometric assay and expression analysis revealed that SP2 is the better manual combination of cis -elements for stress-inducible expression. qRT-PCR results further demonstrated that designed synthetic promoters are not tissue-specific and thus active in different parts upon treatment with osmotic stress in Arabidopsis plants. This study provides tools for transcriptional upgradation of valuable crops against drought stress and adds to the current knowledge of synthetic biology.

Non-coding de novo mutations in chromatin interactions are implicated in autism spectrum disorder

Three-dimensional chromatin interactions regulate gene expressions. The significance of de novo mutations (DNMs) in chromatin interactions remains poorly understood for autism spectrum disorder (ASD). We generated 813 whole-genome sequences from 242 Korean simplex families to detect DNMs, and identified target genes which were putatively affected by non-coding DNMs in chromatin interactions. Non-coding DNMs in chromatin interactions were significantly involved in transcriptional dysregulations related to ASD risk. Correspondingly, target genes showed spatiotemporal expressions relevant to ASD in developing brains and enrichment in biological pathways implicated in ASD, such as histone modification. Regarding clinical features of ASD, non-coding DNMs in chromatin interactions particularly contributed to low intelligence quotient levels in ASD probands. We further validated our findings using two replication cohorts, Simons Simplex Collection (SSC) and MSSNG, and showed the consistent enrichment of non-coding DNM-disrupted chromatin interactions in ASD probands. Generating human induced pluripotent stem cells in two ASD families, we were able to demonstrate that non-coding DNMs in chromatin interactions alter the expression of target genes at the stage of early neural development. Taken together, our findings indicate that non-coding DNMs in ASD probands lead to early neurodevelopmental disruption implicated in ASD risk via chromatin interactions.

Evolution of the ErbB gene family and analysis of regulators of Egfr expression during development of the rat spinal cord

Egfr, a member of the ErbB gene family, plays a critical role in tissue development and homeostasis, wound healing, and disease. However, expression and regulators of Egfr during spinal cord development remain poorly understood. In this study, we investigated ErbB evolution and analyzed co-expression modules, miRNAs, and transcription factors that may regulate Egfr expression in rats. We found that ErbB family members formed via Egfr duplication in the ancient vertebrates but diverged after speciation of gnathostomes. We identified a module that was co-expressed with Egfr, which involved cell proliferation and blood vessel development. We predicted 25 miRNAs and nine transcription factors that may regulate Egfr expression. Dual-luciferase reporter assays showed six out of nine transcription factors significantly affected Egfr promoter reporter activity. Two of these transcription factors (KLF1 and STAT3) inhibited the Egfr promoter reporter, whereas four transcription factors (including FOXA2) activated the Egfr promoter reporter. Real-time PCR and immunofluorescence experiments showed high expression of FOXA2 during the embryonic period and FOXA2 was expressed in the floor plate of the spinal cord, suggesting the importance of FOXA2 during embryonic spinal cord development. Considering the importance of Egfr in embryonic spinal cord development, wound healing, and disease (specifically in cancer), regulatory elements identified in this study may provide candidate targets for nerve regeneration and disease treatment in the future.

Gender-related differentially expressed genes in pancreatic cancer: possible culprits or accomplices?

Pancreatic cancer (PC) is one of the leading causes of cancer mortality worldwide, and its incidence and mortality rate in several regions is higher in male patients. Although numerous efforts have been made to enhance the clinical outcomes of existing therapeutic regimens, their efficiency is still low, and drug resistance usually occurs in many patients. In addition, the exact underlying molecular basis that makes PC slightly more prevalent among males remains unknown. Providing information regarding the possible association between gender and PC tumorigenesis may offer important clues for how certain molecular cross-talks can affect PC initiation and/or progression. In this study, we used several microarray expression data to identify the common up- and downregulated genes within one specific gender, which were also specified to have binding sites for androgen and/or estrogen receptors. Using functional enrichment analysis among the others, for all the gene sets found in this study, we have shed light on the plausible importance of the androgenic effectors in tumorigenesis, such as the androgen-regulated expression of the GLI transcription factor and the potential role of testosterone in the extracellular matrix (ECM)–cell interaction, which are known for their importance in tumorigenesis. Moreover, we demonstrated that the biological process axon guidance was highlighted regarding the upregulated genes in male patients. Overall, identification of gene candidates as the possible link between gender and PC progression or survival rates may help in developing strategies to reduce the incidence of this cancer.

Motif and conserved module analysis in DNA (promoters, enhancers) and RNA (lncRNA, mRNA) using AlModules

Nucleic acid motifs consist of conserved and variable nucleotide regions. For functional action, several motifs are combined to modules. The tool AIModules allows identification of such motifs including combinations of them and conservation in several nucleic acid stretches. AIModules recognizes conserved motifs and combinations of motifs (modules) allowing a number of interesting biological applications such as analysis of promoter and transcription factor binding sites (TFBS), identification of conserved modules shared between several gene families, e.g. promoter regions, but also analysis of shared and conserved other DNA motifs such as enhancers and silencers, in mRNA (motifs or regulatory elements e.g. for polyadenylation) and lncRNAs. The tool AIModules presented here is an integrated solution for motif analysis, offered as a Web service as well as downloadable software. Several nucleotide sequences are queried for TFBSs using predefined matrices from the JASPAR DB or by using one's own matrices for diverse types of DNA or RNA motif discovery. Furthermore, AIModules can find TFBSs common to two or more sequences. Demanding high or low conservation, AIModules outperforms other solutions in speed and finds more modules (specific combinations of TFBS) than alternative available software. The application also searches RNA motifs such as polyadenylation site or RNA-protein binding motifs as well as DNA motifs such as enhancers as well as user-specified motif combinations ( https://bioinfo-wuerz.de/aimodules/ ; alternative entry pages: https://aimodules.heinzelab.de or https://www.biozentrum.uni-wuerzburg.de/bioinfo/computing/aimodules ). The application is free and open source whether used online, on-site, or locally.