Common Regulatory Network Research Articles

Transcriptional mechanisms of phenotypic plasticity in glioblastoma.

e14027 Background: Primary glioblastoma (GBM) tumor cells exhibit plasticity and can adopt different cellular states that exhibit dynamic transitions between states. This could be one of the mechanisms for treatment escape and leads to poor prognosis in GBM. In this regard, we sought to uncover transcriptional regulatory mechanisms responsible for this phenotypic plasticity in GBM. Methods: Standard of care (SOC) treatment (radiation and temozolomide) modulates transcriptional regulatory networks (TRNs) of the tumor cells and drives treatment resistance. Since TRNs could be key mechanistic drivers for state transitions, integrated single-cell RNA sequencing (scRNAseq) data of primary and relapse patient samples (11) were analyzed with scMINER. Our method uses a systems biology framework to identify cell states (St_1 to St_n) in terms of regulatory networks and their activities. Clustering of co-regulated gene modules (regulons), enabled the identification of transcriptional programs and regulatory mechanisms that specify tumor cell states. Results: Activity patterns of 38 transcriptional programs clustered tumor cells from the patients into 21 states. Preliminarily, we observed common regulatory network cell states (St_1, St_2, St_5, St_6) across patient tumors. Additionally, primary tumor cells of a patient (UW14) transitioned from initial cell state(s) observed across multiple patients (e.g., St_1) to a recurrent state (e.g., St_9) through a trajectory driven by transcription factors (TFs). Herein, scMINER regulatory network inference method also defines master TFs (TFs that regulate the majority of the other TFs in the network) that were used to characterize initial (St_1) or recurrent (St_9) cell states. Moreover, certain transcription factors (22) and their co-expressed target genes (regulons) also exhibit variation in expression across the pseudotime trajectory (from St_1 to St_9) including masterTFs such as SOX9 and MECOM. Conclusions: Therefore, we observed that shared cell states may exist across patients (primary and relapse) although in different proportions. Our analysis suggests a de-differentiation trajectory from neuronal cells to developing GBM stem cells and predicts transcription factors that could be playing an important role therein. Moreover, similar transcriptional regulatory mechanisms may promote transitions between primary and recurrent tumor cell states that exhibit phenotypic plasticity.

CHD7 and SOX2 act in a common gene regulatory network during mammalian semicircular canal and cochlear development

Inner ear morphogenesis requires tightly regulated epigenetic and transcriptional control of gene expression. CHD7, an ATP-dependent chromodomain helicase DNA-binding protein, and SOX2, an SRY-related HMG box pioneer transcription factor, are known to contribute to vestibular and auditory system development, but their genetic interactions in the ear have not been explored. Here, we analyzed inner ear development and the transcriptional regulatory landscapes in mice with variable dosages of Chd7 and/or Sox2. We show that combined haploinsufficiency for Chd7 and Sox2 results in reduced otic cell proliferation, severe malformations of semicircular canals, and shortened cochleae with ectopic hair cells. Examination of mice with conditional, inducible Chd7 loss by Sox2CreER reveals a critical period (~E9.5) of susceptibility in the inner ear to combined Chd7 and Sox2 loss. Data from genome-wide RNA-sequencing and CUT&Tag studies in the otocyst show that CHD7 regulates Sox2 expression and acts early in a gene regulatory network to control expression of key otic patterning genes, including Pax2 and Otx2. CHD7 and SOX2 directly bind independently and cooperatively at transcription start sites and enhancers to regulate otic progenitor cell gene expression. Together, our findings reveal essential roles for Chd7 and Sox2 in early inner ear development and may be applicable for syndromic and other forms of hearing or balance disorders.

Exploring large-scale gene coexpression networks in peach (Prunus persica L.): a new tool for predicting gene function.

Peach is a model for Prunus genetics and genomics, however, identifying and validating genes associated to peach breeding traits is a complex task. A gene coexpression network (GCN) capable of capturing stable gene-gene relationships would help researchers overcome the intrinsic limitations of peach genetics and genomics approaches and outline future research opportunities. In this study, we created four GCNs from 604 Illumina RNA-Seq libraries. We evaluated the performance of every GCN in predicting functional annotations using an algorithm based on the 'guilty-by-association' principle. The GCN with the best performance was COO300, encompassing 21 956 genes. To validate its performance predicting gene function, we performed two case studies. In case study 1, we used two genes involved in fruit flesh softening: the endopolygalacturonases PpPG21 and PpPG22. Genes coexpressing with both genes were extracted and referred to as melting flesh (MF) network. Finally, we performed an enrichment analysis of MF network and compared the results with the current knowledge regarding peach fruit softening. The MF network mostly included genes involved in cell wall expansion and remodeling, and with expressions triggered by ripening-related phytohormones, such as ethylene, auxin, and methyl jasmonate. In case study 2, we explored potential targets of the anthocyanin regulator PpMYB10.1 by comparing its gene-centered coexpression network with that of its grapevine orthologues, identifying a common regulatory network. These results validated COO300 as a powerful tool for peach and Prunus research. This network, renamed as PeachGCN v1.0, and the scripts required to perform a function prediction analysis are available at https://github.com/felipecobos/PeachGCN.

Common and divergent gene regulatory networks control injury-induced and developmental neurogenesis in zebrafish retina

Following acute retinal damage, zebrafish possess the ability to regenerate all neuronal subtypes through Müller glia (MG) reprogramming and asymmetric cell division that produces a multipotent Müller glia-derived neuronal progenitor cell (MGPC). This raises three key questions. First, do MG reprogram to a developmental retinal progenitor cell (RPC) state? Second, to what extent does regeneration recapitulate retinal development? And finally, does loss of different retinal cell subtypes induce unique MG regeneration responses? We examined these questions by performing single-nuclear and single-cell RNA-Seq and ATAC-Seq in both developing and regenerating retinas. Here we show that injury induces MG to reprogram to a state similar to late-stage RPCs. However, there are major transcriptional differences between MGPCs and RPCs, as well as major transcriptional differences between activated MG and MGPCs when different retinal cell subtypes are damaged. Validation of candidate genes confirmed that loss of different subtypes induces differences in transcription factor gene expression and regeneration outcomes.

Functional relevance of circRNA aberrant expression in pediatric acute leukemia with KMT2A::AFF1 fusion

AbstractCircular RNAs (circRNAs) are emerging molecular players in leukemogenesis and promising therapeutic targets. In KMT2A::AFF1 (MLL::AF4)-rearranged leukemia, an aggressive disease compared with other pediatric B-cell precursor (BCP) acute lymphoblastic leukemia (ALL), data about circRNAs are limited. Here, we disclose the circRNA landscape of infant patients with KMT2A::AFF1 translocated BCP-ALL showing dysregulated, mostly ectopically expressed, circRNAs in leukemia cells. Most of these circRNAs, apart from circHIPK3 and circZNF609, previously associated with oncogenic behavior in ALL, are still uncharacterized. An in vitro loss-of-function screening identified an oncogenic role of circFKBP5, circKLHL2, circNR3C1, and circPAN3 in KMT2A::AFF1 ALL, whose silencing affected cell proliferation and apoptosis. Further study in an extended cohort disclosed a significantly correlated expression of these oncogenic circRNAs and their putative involvement in common regulatory networks. Moreover, it showed that circAFF1 upregulation occurs in a subset of cases with HOXA KMT2A::AFF1 ALL. Collectively, functional analyses and patient data reveal oncogenic circRNA upregulation as a relevant mechanism that sustains the malignant cell phenotype in KMT2A::AFF1 ALL.

OR11-01 A New Model for Estrogen Receptor-alpha Antagonist Mechanism of Action

Abstract Disclosure: J.C. Nwachukwu: None. J.W. Njeri: None. T. Venables: None. C. Seath: None. M.E. Pipkin: None. Y. Hou: None. B.S. Katzenellenbogen: None. J.A. Katzenellenbogen: None. K.W. Nettles: None. We recently developed Dual Mechanism Estrogen Receptor-α (ERα) Inhibitors (DMERI) with more antagonist efficacy than selective ERα modulators (SERMs) such as tamoxifen, and selective ERα degraders (SERDs) such as fulvestrant in certain breast cancer models. Current models of anti-estrogen action are largely based on genomics studies in growth arrested cells, where the effects primarily reversed estradiol-mediated gene expression. In that context, most ER antagonists, especially the SERDs, show very little gene activation. We were interested in understanding antagonist mechanism of action in asynchronous cells that were growing in intact medium that had not been charcoal-stripped to remove estradiol and growth factors. RNA-seq of asynchronous MCF-7 cells treated with 4-hydroxytamoxifen, fulvestrant, a SERM-like DMERI, and a SERD-like DMERI showed that all four classes of ligands equally up-regulated and down-regulated differentially expressed genes (DEGs), which included a core set of regulatory networks, as well as ligand-selective gene networks. The number of DEGs predicted efficacy, with tamoxifen regulating few genes, and the SERD-like DMERI regulating the most genes. The SERD-like DMERI also inhibited expression of many more genes that encode cell cycle regulators and key cancer drivers, such as ERα, AKT1, HRAS, and several growth factor receptors. Landscape in silico analysis based on chromatin immunoprecipitation (ChIP)-seq data suggests that recruitment of the ERα MegaTrans complex in collaboration with the oxysterol receptor, LXRα, predicts up-regulation, while cooperation with the retinoblastoma-associated proteins, RB1 and RBL1, predicts down-regulation of target genes. APEX2-mediated proximity labeling in asynchronous cells revealed that fulvestrant and the SERD-like DMERI labeled both coactivators and corepressors. Both compounds also increase labeling of STAG1 and other cohesin subunits that coordinate chromatin looping, while also associating with the nuclear matrix protein Lamin B1 and heterochromatin-associated proteins. These data suggest that antagonists have dual roles in genome organization. An shRNAmiR screen targeting over 300 chromatin regulatory factors revealed networks of shared and exclusive coregulators used by antiestrogens, including the heterochromatin protein, CBX1, which was exclusively required for growth inhibition by fulvestrant, and subunits of histone acetylase and methylase complexes such as KMT2A (the catalytic subunit of the MLL1 complex), which were broadly required. Thus, antiestrogen-driven growth suppression involves both common and ligand-selective transcriptional regulatory networks through utilization of both coactivators and corepressors. Supported by NIH R01CA220284 and BCRF-083 and BCRF-084 grants. Presentation Date: Friday, June 16, 2023

Identifying cellular senescence associated genes involved in the progression of end-stage renal disease as new biomarkers

BackgroundCellular senescence plays an essential role in the development and progression of end-stage renal disease (ESRD). However, the detailed mechanisms phenomenon remains unclear.MethodsThe mRNA expression profiling dataset GSE37171 was taken from the Gene Expression Omnibus (GEO) database. The cell senescence-associated hub genes were selected by applying protein–protein interaction (PPI), followed by correlation analysis, gene interaction analysis, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. We next explored the relationships of hub genes with miRNAs, TFs, and diseases. The absolute abundance of eight immune cells and two stromal cells were calculated by MCPcount and the correlation of hub genes with these ten cells was analyzed. Lasso was used to selecting for trait genes. ROC curves and DCA decision curves were used to assess the accuracy and predictive power of the trait genes.ResultsA total of 65 cellular senescence signature genes were identified among patients and controls. The PPI network screened out ten hub genes. GO and KEGG indicated that ten hub genes were associated with ESRD progression. Transcription factor gene interactions and common regulatory networks of miRNAs were also identified in the datasets. The hub genes were significantly correlated with immune cells and stromal cells. Then the lasso model was constructed to screen out the five most relevant signature genes (FOS, FOXO3, SIRT1, TP53, SMARCA4). The area under the ROC curve (AUC) showed that these five characteristic genes have good resolving power for the diagnostic model.ConclusionsOur findings suggested that cellular senescence-associated genes played an important role in the development of ESRD and immune regulation.

The underlying molecular mechanisms and biomarkers between periodontitis and COVID-19

ObjectiveEmerging evidence shows the clinical consequences of patient with COVID-19 and periodontitis are not promising, and periodontitis is a risk factor. Periodontitis and COVID-19 probably have a relationship. Hence, this study aimed to identify the common molecular mechanism that may help to devise potential therapeutic strategies in the future.Material and methodsWe analyzed two RNA-seq datasets for differential expressed genes, enrichment of biological processes, transcription factors (TFs) and deconvolution-based immune cell types in periodontitis, COVID-19 and healthy controls. Relationships between TFs and mRNA were established by Pearson correlation analysis, and the common TFs-mRNA regulatory network and nine co-upregulated TFs of the two diseases was obtained. The RT-PCR detected the TFs.ResultsA total of 1616 and 10201 differentially expressed gene (DEGs) from periodontitis and COVID-19 are found. Moreover, nine shared TFs and common biological processes associated with lymphocyte activation involved in immune response were identified across periodontitis and COVID-19. The cell type enrichment revealed elevated plasma cells among two diseases. The RT-PCR further confirmed the nine TFs up-regulation in periodontitis.ConclusionThe pathogenesis of periodontitis and COVID-19 is closely related to the expression of TFs and lymphocyte activation, which can provide potential targets for treatment.

Conserved hierarchical gene regulatory networks for drought and cold stress response in Myrica rubra.

Stress response in plant is regulated by a large number of genes co-operating in diverse networks that serve multiple adaptive process. To understand how gene regulatory networks (GRNs) modulating abiotic stress responses, we compare the GRNs underlying drought and cold stresses using samples collected at 4 or 6h intervals within 48h in Chinese bayberry (Myrica rubra). We detected 7,583 and 8,840 differentially expressed genes (DEGs) under drought and cold stress respectively, which might be responsive to environmental stresses. Drought- and cold-responsive GRNs, which have been built according to the timing of transcription under both abiotic stresses, have a conserved trans-regulator and a common regulatory network. In both GRNs, basic helix-loop-helix family transcription factor (bHLH) serve as central nodes. MrbHLHp10 transcripts exhibited continuous increase in the two abiotic stresses and acts upstream regulator of ASCORBATE PEROXIDASE (APX) gene. To examine the potential biological functions of MrbHLH10, we generated a transgenic Arabidopsis plant that constitutively overexpresses the MrbHLH10 gene. Compared to wild-type (WT) plants, overexpressing transgenic Arabidopsis plants maintained higher APX activity and biomass accumulation under drought and cold stress. Consistently, RNAi plants had elevated susceptibility to both stresses. Taken together, these results suggested that MrbHLH10 mitigates abiotic stresses through the modulation of ROS scavenging.

Analysis of co-expression gene network associated with intracranial aneurysm and type 2 diabetes mellitus.

To screen for common target genes in intracranial aneurysms (IA) and type 2 diabetes mellitus (T2DM), construct a common transcriptional regulatory network to predict clusters of candidate genes involved in the pathogenesis of T2DM and IA, and identify the common neurovascular markers and pathways in T2DM causing IA. Microarray datasets (GSE55650, GSE25462, GSE26969, GSE75436, and GSE13353) from the GEO database were analyzed in this research. Screening of the IA and the T2DM datasets yielded a total of 126 DEGs, among which 78 were upregulated and 138 were downregulated. Functional enrichment analysis revealed that these DEGs were enriched for a total of 68 GO pathways, including extracellular matrix composition, coagulation regulation, hemostasis regulation, and collagen fiber composition pathways. We also constructed transcriptional regulatory networks, and identified key transcription factors involved in both the conditions. Univariate logistic regression analysis showed that ARNTL2 and STAT1 were significantly associated with the development of T2DM and IA, acting as the common neurovascular markers for both the diseases. In cellular experiments, hyperglycemic microenvironments exhibited upregulated STAT1 expression. STAT1 may be involved in the pathogenesis of IA in T2DM patients. Being the common neurovascular markers, STAT1 may acts as novel therapeutic targets for the treatment of IA and T2DM.

Transposable elements orchestrate subgenome-convergent and -divergent transcription in common wheat

The success of common wheat as a global staple crop was largely attributed to its genomic diversity and redundancy due to the merge of different genomes, giving rise to the major question how subgenome-divergent and -convergent transcription is mediated and harmonized in a single cell. Here, we create a catalog of genome-wide transcription factor-binding sites (TFBSs) to assemble a common wheat regulatory network on an unprecedented scale. A significant proportion of subgenome-divergent TFBSs are derived from differential expansions of particular transposable elements (TEs) in diploid progenitors, which contribute to subgenome-divergent transcription. Whereas subgenome-convergent transcription is associated with balanced TF binding at loci derived from TE expansions before diploid divergence. These TFBSs have retained in parallel during evolution of each diploid, despite extensive unbalanced turnover of the flanking TEs. Thus, the differential evolutionary selection of paleo- and neo-TEs contribute to subgenome-convergent and -divergent regulation in common wheat, highlighting the influence of TE repertory plasticity on transcriptional plasticity in polyploid.

Pathogen-induced biosynthetic pathways encode defense-related molecules in bread wheat

Wheat is a widely grown food crop that suffers major yield losses due to attack by pests and pathogens. A better understanding of biotic stress responses in wheat is thus of major importance. The recently assembled bread wheat genome coupled with extensive transcriptomic resources provides unprecedented new opportunities to investigate responses to pathogen challenge. Here, we analyze gene coexpression networks to identify modules showing consistent induction in response to pathogen exposure. Within the top pathogen-induced modules, we identify multiple clusters of physically adjacent genes that correspond to six pathogen-induced biosynthetic pathways that share a common regulatory network. Functional analysis reveals that these pathways, all of which are encoded by biosynthetic gene clusters, produce various different classes of compounds—namely, flavonoids, diterpenes, and triterpenes, including the defense-related compound ellarinacin. Through comparative genomics, we also identify associations with the known rice phytoalexins momilactones, as well as with a defense-related gene cluster in the grass model plant Brachypodium distachyon. Our results significantly advance the understanding of chemical defenses in wheat and open up avenues for enhancing disease resistance in this agriculturally important crop. They also exemplify the power of transcriptional networks to discover the biosynthesis of chemical defenses in plants with large, complex genomes.

Bioinformatics and Network Pharmacology Identify the Therapeutic Role and Potential Mechanism of Melatonin in AD and Rosacea.

Rosacea is significantly associated with dementia, particularly Alzheimer’s disease (AD). However, the common underlying molecular mechanism connecting these two diseases remains limited. This study aimed to reveal the common molecular regulatory networks and identify the potential therapeutic drugs for rosacea and AD. There were 747 overlapped DEGs (ol-DEGs) that were detected in AD and rosacea, enriched in inflammation-, metabolism-, and apoptosis-related pathways. Using the TF regulatory network analysis, 37 common TFs and target genes were identified as hub genes. They were used to predict the therapeutic drugs for rosacea and AD using the DGIdb/CMap database. Among the 113 predicted drugs, melatonin (MLT) was co-associated with both RORA and IFN-γ in AD and rosacea. Subsequently, network pharmacology analysis identified 19 pharmacological targets of MLT and demonstrated that MLT could help in treating AD/rosacea partly by modulating inflammatory and vascular signaling pathways. Finally, we verified the therapeutic role and mechanism of MLT on rosacea in vivo and in vitro. We found that MLT treatment significantly improved rosacea-like skin lesion by reducing keratinocyte-mediated inflammatory cytokine secretion and repressing the migration of HUVEC cells. In conclusion, this study contributes to common pathologies shared by rosacea and AD and identified MLT as an effective treatment strategy for rosacea and AD via regulating inflammation and angiogenesis.

Single-Cell Gene Expression Analysis Reveals Gene Regulatory Networks Driving Proliferation in Pituitary Stem and Endocrine Cells

A fundamental question for pituitary development and disease is to understand the mechanisms that regulate proliferation, quiescence, and differentiation of stem cells and endocrine lineage-committed precursor cells. Pituitary stem cells, marked by expression of SOX2, are highly proliferative during development and the early postnatal period. This pool becomes quiescent over time, but stem cells retain the ability to re-enter the cell cycle and differentiate into nascent endocrine cells of all lineages in response to physiological demands. The rodent pituitary gland undergoes substantial growth in the postnatal period, and much of this increase in size is due to proliferation of committed progenitors, such as Pou1f1-expressing cells. We performed single-cell RNA transcriptomics analyses of over 8,000 male and female 4-day-old mouse pituitary cells in order to assess stem cell heterogeneity and to identify novel pituitary stem cell biomarkers. We identified a number of factors enriched in pituitary stem cells relative to differentiating cells, including the transcription factors TGIF1 and NR4A3 as well as several members of the nuclear factor I family (NFIA, NFIB, NFIX). We also detected stem cell-specific expression of the cortisol synthesizing enzyme HSD11B1 and folate receptor FOLR1, suggesting novel roles for these pathways in pituitary stem cells. There were few transcriptomic differences between proliferating and non-proliferating stem cells. However, proliferating stem cells and proliferating committed progenitors shared expression of cell-cycle associated genes and novel transcription factors such as BRCA1 and E2F1. Furthermore, single-cell gene network inference and clustering (SCENIC) analyses demonstrated activation of common gene regulatory networks in both proliferating stem and proliferating endocrine populations, including both the E2f1 and Brca1 regulons. RNA velocity, trajectory, and phylogenetic analyses and find that proliferating stem and endocrine cells are likely independent cellular states. Instead, they support the idea that proliferating stem cells become quiescent and transition to committed progenitors that re-enter the cell cycle and then subsequently activate hormone gene expression. In conclusion, our single-cell gene expression analyses of early postnatal pituitary cells have shed light on the developmental trajectory from proliferating stem cell to quiescent differentiated hormone-producing cells.

Genome-wide analysis of European sea bass provides insights into the evolution and functions of single-exon genes.

Several studies have attempted to understand the origin and evolution of single‐exon genes (SEGs) in eukaryotic organisms, including fishes, but few have examined the functional and evolutionary relationships between SEGs and multiple‐exon gene (MEG) paralogs, in particular the conservation of promoter regions. Given that SEGs originate via the reverse transcription of mRNA from a “parental” MEGs, such comparisons may enable identifying evolutionarily‐related SEG/MEG paralogs, which might fulfill equivalent physiological functions. Here, the relationship of SEG proportion with MEG count, gene density, intron count, and chromosome size was assessed for the genome of the European sea bass, Dicentrarchus labrax. Then, SEGs with an MEG parent were identified, and promoter sequences of SEG/MEG paralogs were compared, to identify highly conserved functional motifs. The results revealed a total count of 1,585 (8.3% of total genes) SEGs in the European sea bass genome, which was correlated with MEG count but not with gene density. The significant correlation of SEG content with the number of MEGs suggests that SEGs were continuously and independently generated over evolutionary time following species divergence through retrotranscription events, followed by tandem duplications. Functional annotation showed that the majority of SEGs are functional, as is evident from their expression in RNA‐seq data used to support homology‐based genome annotation. Differences in 5′UTR and 3′UTR lengths between SEG/MEG paralogs observed in this study may contribute to gene expression divergence between them and therefore lead to the emergence of new SEG functions. The comparison of nonsynonymous to synonymous changes (Ka/Ks) between SEG/MEG parents showed that 74 of them are under positive selection (Ka/Ks > 1; p = .0447). An additional fifteen SEGs with an MEG parent have a common promoter, which implies that they are under the influence of common regulatory networks.

Gene Regulatory Network of ETS Domain Transcription Factors in Different Stages of Glioma.

The ETS domain family of transcription factors is involved in a number of biological processes, and is commonly misregulated in various forms of cancer. Using microarray datasets from patients with different grades of glioma, we have analyzed the expression profiles of various ETS genes, and have identified ETV1, ELK3, ETV4, ELF4, and ETV6 as novel biomarkers for the identification of different glioma grades. We have further analyzed the gene regulatory networks of ETS transcription factors and compared them to previous microarray studies, where Elk-1-VP16 or PEA3-VP16 were overexpressed in neuroblastoma cell lines, and we identify unique and common regulatory networks for these ETS proteins.

Identification of Differentially Expressed Long Noncoding RNAs in Mice Intestines After Severe Burns and a Preliminary Study into the Key Gene H19.

The intestine is considered the key organ in stress response to severe burns and injury to the intestine after severe burns can be fatal. However, the injury and subsequent repair of intestinal tissues after severe burns at the genetic level are poorly understood. Long noncoding RNAs (lncRNAs) have important functions in regulating many biological processes, including gene transcription and translation. Autophagy is a process of intracellular degradation and reutilization of cytoplasmic proteins and organelles. We herein analyzed the genome-wide expression profile of lncRNAs and mRNAs after severe burns in the intestines of mice by lncRNA microarray. Quantitative reverse transcription-polymerase chain reaction was performed to verify the reliability of microarray analysis results, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used for bioinformatics analysis of differentially expressed mRNAs. The common regulatory network between the top 10 differentially expressed lncRNAs and trans-related mRNAs were visualized by Cytoscape (v3.7.2). Next, we hypothesized that H19 is the key gene for intestinal mucosal repair. After H19 was overexpressed, the changes in downstream autophagy protein expression levels were observed. GO and KEGG analysis indicated that the differentially expressed mRNAs were mainly enriched in a cell cycle- and mitosis-related genes. Overexpression of lncRNA-H19 showed that the autophagy-related gene Trim21 was upregulated, while HIF1α was downregulated. LncRNA-H19 played a key role in repairing the intestinal mucosa, and overexpression of lncRNA-H19 activated autophagy and migration of intestinal epithelial cells (IEC-6).

Housing microbial symbionts: evolutionary origins and diversification of symbiotic organs in animals.

In many animal hosts, microbial symbionts are housed within specialized structures known as symbiotic organs, but the evolutionary origins of these structures have rarely been investigated. Here, I adopt an evolutionary developmental (evo-devo) approach, specifically to apply knowledge of the development of symbiotic organs to gain insights into their evolutionary origins and diversification. In particular, host genetic changes associated with evolution of symbiotic organs can be inferred from studies to identify the host genes that orchestrate the development of symbiotic organs, recognizing that microbial products may also play a key role in triggering the developmental programme in some associations. These studies may also reveal whether higher animal taxonomic groups (order, class, phylum, etc.) possess a common genetic regulatory network for symbiosis that is latent in taxa lacking symbiotic organs, and activated at the origination of symbiosis in different host lineages. In this way, apparent instances of convergent evolution of symbiotic organs may be homologous in terms of a common genetic blueprint for symbiosis. Advances in genetic technologies, including reverse genetic tools and genome editing, will facilitate the application of evo-devo approaches to investigate the evolution of symbiotic organs in animals. This article is part of the theme issue 'The role of the microbiome in host evolution'.

Analysis of Dynamic Global Transcriptional Atlas Reveals Common Regulatory Networks of Hormones and Photosynthesis Across Nicotiana Varieties in Response to Long-Term Drought.

Land plants evolve drought acclimation. Existing knowledge of gene regulation mainly comes from short-term drought treatment. However, common regulatory mechanism shared by multiple varieties under long-term drought is little explored. Here we investigated changes in physiology, hormones and transcriptomes in leaves of Nicotiana varieties K326 and Basma Xanthi with/without drought treatment at time courses spanning 1 month. Analyses of deep RNA-Seq data and further full-length Iso-Seq data revealed an atlas of dynamic changes of transcripts, spliced isoforms, gene expression, associated Gene Ontology, and metabolism pathways. Fewer differentially expressed genes (DEGs) were induced by drought in high tolerance variety than susceptible variety. Comparison among seven hormone signal pathways identified that genes in both abscisic acid and auxin signaling pathways were highly induced although specific genes were depended on the variety. Common hormone regulatory network analysis revealed that genes encoding clade A protein phosphatase 2C gene (PP2C) in abscisic acid pathway was the pivotal hub. Expressional regulation in photosynthesis was also common and variety specific. We conclude that long-term drought inducing gene regulatory networks of hormones and photosynthesis are variety dependent, and PP2C is the center of the common hormone regulatory network. Thus, this study improves our understanding of gene regulatory network in drought response.

Genetic Network between Leaf Senescence and Plant Immunity: Crucial Regulatory Nodes and New Insights.

Leaf senescence is an essential physiological process that is accompanied by the remobilization of nutrients from senescent leaves to young leaves or other developing organs. Although leaf senescence is a genetically programmed process, it can be induced by a wide variety of biotic and abiotic factors. Accumulating studies demonstrate that senescence-associated transcription factors (Sen-TFs) play key regulatory roles in controlling the initiation and progression of leaf senescence process. Interestingly, recent functional studies also reveal that a number of Sen-TFs function as positive or negative regulators of plant immunity. Moreover, the plant hormone salicylic acid (SA) and reactive oxygen species (ROS) have been demonstrated to be key signaling molecules in regulating leaf senescence and plant immunity, suggesting that these two processes share similar or common regulatory networks. However, the interactions between leaf senescence and plant immunity did not attract sufficient attention to plant scientists. Here, we review the regulatory roles of SA and ROS in biotic and abiotic stresses, as well as the cross-talks between SA/ROS and other hormones in leaf senescence and plant immunity, summarize the transcriptional controls of Sen-TFs on SA and ROS signal pathways, and analyze the cross-regulation between senescence and immunity through a broad literature survey. In-depth understandings of the cross-regulatory mechanisms between leaf senescence and plant immunity will facilitate the cultivation of high-yield and disease-resistant crops through a molecular breeding strategy.