Search

Nitrogen is an essential nutrient for plants and a major determinant of plant growth and crop yield. Plants acquire nitrogen mainly in the form of nitrate and ammonium. Both nitrogen sources affect plant responses and signaling pathways in a different way, but these signaling pathways interact, complicating the study of nitrogen responses. Extensive transcriptome analyses and the construction of gene regulatory networks, mainly in response to nitrate, have significantly advanced our understanding of nitrogen signaling and responses in model plants and crops. In this study, we aimed to generate a more comprehensive gene regulatory network for the major crop, rice, by incorporating the interactions between ammonium and nitrate. To achieve this, we assessed transcriptome changes in rice roots and shoots over an extensive time course under single or combined applications of the two nitrogen sources. This dataset enabled us to construct a holistic co-expression network and identify potential key regulators of nitrogen responses. Next to known transcription factors, we identified multiple new candidates, including the transcription factors OsRLI and OsEIL1, which we demonstrated to induce the primary nitrate-responsive genes OsNRT1.1b and OsNIR1. Our network thus serves as a valuable resource to obtain novel insights in nitrogen signaling. [ABSTRACT FROM AUTHOR]

Introduction: Hemibiotrophic Phytophthora are a group of agriculturally and ecologically important pathogenic oomycetes causing severe decline in plant growth and fitness. The lifestyle of these pathogens consists of an initial biotrophic phase followed by a switch to a necrotrophic phase in the latter stages of infection. Between these two phases is the biotrophic to necrotrophic switch (BNS) phase, the timing and controls of which are not well understood particularly in Phytophthora spp. where host resistance has a purely quantitative genetic basis. Methods: To investigate this we sequenced and annotated the genome of Phytophthora medicaginis, causal agent of root rot and substantial yield losses to Fabaceae hosts. We analyzed the transcriptome of P. medicaginis across three phases of colonization of a susceptible chickpea host (Cicer arietinum) and performed co-regulatory analysis to identify putative small secreted protein (SSP) effectors that influence timing of the BNS in a quantitative pathosystem. Results: The genome of P. medicaginis is ~78 Mb, comparable to P. fragariae and P. rubi which also cause root rot. Despite this, it encodes the second smallest number of RxLR (arginine-any amino acid-leucine-arginine) containing proteins of currently sequenced Phytophthora species. Only quantitative resistance is known in chickpea to P. medicaginis, however, we found that many RxLR, Crinkler (CRN), and Nep1-like protein (NLP) proteins and carbohydrate active enzymes (CAZymes) were regulated during infection. Characterization of one of these, Phytmed_10271, which encodes an RxLR effector demonstrates that it plays a role in the timing of the BNS phase and root cell death. Discussion: These findings provide an important framework and resource for understanding the role of pathogenicity factors in purely quantitative Phytophthora pathosystems and their implications to the timing of the BNS phase. [ABSTRACT FROM AUTHOR]

N6-methyladenosine (m6A) RNA modification is a conserved mechanism that regulates the fate of RNA across eukaryotic organisms. Despite its significance, a comprehensive analysis of m6A-related genes in non-model plants, such as kiwifruit, is lacking. Here, we identified 36 m6Arelated genes in the kiwifruit genome according to homology and phylogenetic inference. We performed bioinformatics and evolutionary analyses of the writer, eraser, and reader families of m6A modification. Reanalysis of public RNA-seq data collected from samples under various biotic and abiotic stresses indicated that most m6A-related genes were remarkably expressed under different conditions. Through construction of gene co-expression networks, we found significant correlations between several m6A-related genes and transcription factors (TFs) as well as receptor-like genes during the development and ripening of kiwifruit. Furthermore, we performed ATAC-seq assays on diverse kiwifruit tissues to investigate the regulatory mechanisms of m6A-related genes. We identified 10 common open chromatin regions that were present in at least two tissues, and these regions might serve as potential binding sites for MADS protein, C2H2 protein, and other predicted TFs. Our study offers comprehensive insights into the gene family of m6A-related components in kiwifruit, which will lay foundation for exploring mechanisms of post-transcriptional regulation involved in development and adaptation of kiwifruit. [ABSTRACT FROM AUTHOR]

Background: Although rheumatoid arthritis (RA) is a chronic systemic tissue disease often accompanied by osteoporosis (OP), the molecular mechanisms underlying this association remain unclear. This study aimed to elucidate the pathogenesis of RA and OP by identifying differentially expressed mRNAs (DEmRNAs) and long non-coding RNAs (lncRNAs) using a bioinformatics approach. Methods: Expression profiles of individuals diagnosed with OP and RA were retrieved from the Gene Expression Omnibus database. Differential expression analysis was conducted. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) pathway enrichment analyses were performed to gain insights into the functional categories and molecular/biochemical pathways associated with DEmRNAs. We identified the intersection of common DEmRNAs and lncRNAs and constructed a protein-protein interaction (PPI) network. Correlation analysis between the common DEmRNAs and lncRNAs facilitated the construction of a coding-non-coding network. Lastly, serum peripheral blood mononuclear cells (PBMCs) from patients with RA and OP, as well as healthy controls, were obtained for TRAP staining and qRT-PCR to validate the findings obtained from the online dataset assessments. Results: A total of 28 DEmRNAs and 2 DElncRNAs were identified in individuals with both RA and OP. Chromosomal distribution analysis of the consensus DEmRNAs revealed that chromosome 1 had the highest number of differential expression genes. GO and KEGG analyses indicated that these DEmRNAs were primarily associated with " platelets (PLTs) degranulation", "platelet alpha granules", "platelet activation", "tight junctions" and "leukocyte transendothelial migration", with many genes functionally related to PLTs. In the PPI network, MT-ATP6 and PTGS1 emerged as potential hub genes, with MT-ATP6 originating from mitochondrial DNA. Co-expression analysis identified two key lncRNA-mRNA pairs: RP11 − 815J21.2 with MT − ATP6 and RP11 − 815J21.2 with PTGS1. Experimental validation confirmed significant differential expression of RP11-815J21.2, MT-ATP6 and PTGS1 between the healthy controls and the RA + OP groups. Notably, knockdown of RP11-815J21.2 attenuated TNF + IL-6-induced osteoclastogenesis. Conclusions: This study successfully identified shared dysregulated genes and potential therapeutic targets in individuals with RA and OP, highlighting their molecular similarities. These findings provide new insights into the pathogenesis of RA and OP and suggest potential avenues for further research and targeted therapies. [ABSTRACT FROM AUTHOR]

The evolutionary analysis showed that the GF14 family was conserved, however, there was limited evidence linking GF14s to plant height. In our investigations, we discovered a co-expression relationship between ZmGF14s and functionally characterized genes linked to plant height. In the co-expression network, we identified ZmGF14-3, a gene expression exhibiting a positive correlation with plant height in three maize varieties, we postulated that this gene could be intimately linked to plant height development. Subsequently, we cloned ZmGF14-3 from the maize B73 inbred line and overexpressed it in Arabidopsis, resulting in markedly dwarfed transgenic phenotypes. Measurements of endogenous phytohormones disclosed a significant reduction in concentrations of Gibberellic Acid 7 (GA7) and Indole-3-Acetic Acid (IAA) in the overexpressed Arabidopsis, furthermore, qPCR results highlighted a pronounced decrease in the expression levels of plant height-related genes when compared to the wild type, therefore, it is plausible to posit that ZmGF14-3 plays a pivotal role in regulating the growth and development of maize through interactions with various phytohormone-related genes. Thus, delving into the potential interactions between ZmGF14-3 and these genes holds the promise of yielding valuable insights into the molecular mechanisms underpinning plant height development in maize. [ABSTRACT FROM AUTHOR]

Background: Early blight and brown leaf spot are often cited as the most problematic pathogens of tomato in many agricultural regions. Their causal agents are Alternaria spp., a genus of Ascomycota containing numerous necrotrophic pathogens. Breeding programs have yielded quantitatively resistant commercial cultivars, but fungicide application remains necessary to mitigate the yield losses. A major hindrance to resistance breeding is the complexity of the genetic determinants of resistance and susceptibility. In the absence of sufficiently resistant germplasm, we sequenced the transcriptomes of Heinz 1706 tomatoes treated with strongly virulent and weakly virulent isolates of Alternaria spp. 3 h post infection. We expanded existing functional gene annotations in tomato and using network statistics, we analyzed the transcriptional modules associated with defense and susceptibility. Results: The induced responses are very distinct. The weakly virulent isolate induced a defense response of calcium-signaling, hormone responses, and transcription factors. These defense-associated processes were found in a single transcriptional module alongside secondary metabolite biosynthesis genes, and other defense responses. Co-expression and gene regulatory networks independently predicted several D clade ethylene response factors to be early regulators of the defense transcriptional module, as well as other transcription factors both known and novel in pathogen defense, including several JA-associated genes. In contrast, the strongly virulent isolate elicited a much weaker response, and a separate transcriptional module bereft of hormone signaling. Conclusions: Our findings have predicted major defense regulators and several targets for downstream functional analyses. Combined with our improved gene functional annotation, they suggest that defense is achieved through induction of Alternaria-specific immune pathways, and susceptibility is mediated by modulating hormone responses. The implication of multiple specific clade D ethylene response factors and upregulation of JA-associated genes suggests that host defense in this pathosystem involves ethylene response factors to modulate jasmonic acid signaling. [ABSTRACT FROM AUTHOR]

Abstract Background Although rheumatoid arthritis (RA) is a chronic systemic tissue disease often accompanied by osteoporosis (OP), the molecular mechanisms underlying this association remain unclear. This study aimed to elucidate the pathogenesis of RA and OP by identifying differentially expressed mRNAs (DEmRNAs) and long non-coding RNAs (lncRNAs) using a bioinformatics approach. Methods Expression profiles of individuals diagnosed with OP and RA were retrieved from the Gene Expression Omnibus database. Differential expression analysis was conducted. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) pathway enrichment analyses were performed to gain insights into the functional categories and molecular/biochemical pathways associated with DEmRNAs. We identified the intersection of common DEmRNAs and lncRNAs and constructed a protein-protein interaction (PPI) network. Correlation analysis between the common DEmRNAs and lncRNAs facilitated the construction of a coding-non-coding network. Lastly, serum peripheral blood mononuclear cells (PBMCs) from patients with RA and OP, as well as healthy controls, were obtained for TRAP staining and qRT-PCR to validate the findings obtained from the online dataset assessments. Results A total of 28 DEmRNAs and 2 DElncRNAs were identified in individuals with both RA and OP. Chromosomal distribution analysis of the consensus DEmRNAs revealed that chromosome 1 had the highest number of differential expression genes. GO and KEGG analyses indicated that these DEmRNAs were primarily associated with " platelets (PLTs) degranulation”, “platelet alpha granules”, “platelet activation”, “tight junctions” and “leukocyte transendothelial migration”, with many genes functionally related to PLTs. In the PPI network, MT-ATP6 and PTGS1 emerged as potential hub genes, with MT-ATP6 originating from mitochondrial DNA. Co-expression analysis identified two key lncRNA-mRNA pairs: RP11 − 815J21.2 with MT − ATP6 and RP11 − 815J21.2 with PTGS1. Experimental validation confirmed significant differential expression of RP11-815J21.2, MT-ATP6 and PTGS1 between the healthy controls and the RA + OP groups. Notably, knockdown of RP11-815J21.2 attenuated TNF + IL-6-induced osteoclastogenesis. Conclusions This study successfully identified shared dysregulated genes and potential therapeutic targets in individuals with RA and OP, highlighting their molecular similarities. These findings provide new insights into the pathogenesis of RA and OP and suggest potential avenues for further research and targeted therapies.

Sesuvium portulacastrum L., a perennial facultative halophyte, is extensively distributed across tropical and subtropical coastal regions. Its limited cold tolerance significantly impacts both the productivity and the geographical distribution of this species in higher-latitude areas. In this study, we employed RNA-Seq technology to delineate the transcriptomic alterations in Sesuvium plants exposed to low temperatures, thus advancing our comprehension of the molecular underpinnings of this physiological adaptation and root formation. Our findings demonstrated differential expression of 10,805, 16,389, and 10,503 genes in the low versus moderate temperature (LT vs. MT), moderate versus high temperature (MT vs. HT), and low versus high temperature (LT vs. HT) comparative analyses, respectively. Notably, the gene categories "structural molecule activity", "ribosome biogenesis", and "ribosome" were particularly enriched among the LT vs. HT-specific differentially expressed genes (DEGs). When synthesizing the insights from these three comparative studies, the principal pathways associated with the cold response mechanism were identified as "carbon fixation in photosynthetic organisms", "starch and sucrose metabolism", "plant hormone signal transduction", "glycolysis/gluconeogenesis", and "photosynthesis". In addition, we elucidated the involvement of auxin signaling pathways, adventitious root formation (ARF), lateral root formation (LRF), and novel genes associated with shoot system development in root formation. Subsequently, we constructed a network diagram to investigate the interplay between hormone levels and pivotal genes, thereby clarifying the regulatory pathways of plant root formation under low-temperature stress and isolating key genes instrumental in root development. This study has provided critical insights into the molecular mechanisms that facilitate the adaptation to cold stress and root formation in S. portulacastrum. [ABSTRACT FROM AUTHOR]

Abstract Background Soil salinization is one of the vital factors threatening the world’s food security. To reveal the biological mechanism of response to salt stress in wheat, this study was conducted to resolve the transcription level difference to salt stress between CM6005 (salt-tolerant) and KN9204 (salt-sensitive) at the germination and seedling stage. Results To investigate the molecular mechanism underlying salt tolerance in wheat, we conducted comprehensive transcriptome analyses at the seedling and germination stages. Two wheat cultivars, CM6005 (salt-tolerant) and KN9204 (salt-sensitive) were subjected to salt treatment, resulting in a total of 24 transcriptomes. Through expression-network analysis, we identified 17 modules, 16 and 13 of which highly correlate with salt tolerance-related phenotypes in the germination and seedling stages, respectively. Moreover, we identified candidate Hub genes associated with specific modules and explored their regulatory relationships using co-expression data. Enrichment analysis revealed specific enrichment of gibberellin-related terms and pathways in CM6005, highlighting the potential importance of gibberellin regulation in enhancing salt tolerance. In contrast, KN9204 exhibited specific enrichment in glutathione-related terms and activities, suggesting the involvement of glutathione-mediated antioxidant mechanisms in conferring resistance to salt stress. Additionally, glucose transport was found to be a fundamental mechanism for salt tolerance during wheat seedling and germination stages, indicating its potential universality in wheat. Wheat plants improve their resilience and productivity by utilizing adaptive mechanisms like adjusting osmotic balance, bolstering antioxidant defenses, accumulating compatible solutes, altering root morphology, and regulating hormones, enabling them to better withstand extended periods of salt stress. Conclusion Through utilizing transcriptome-level analysis employing WGCNA, we have revealed a potential regulatory mechanism that governs the response to salt stress and recovery in wheat cultivars. Furthermore, we have identified key candidate central genes that play a crucial role in this mechanism. These central genes are likely to be vital components within the gene expression network associated with salt tolerance. The findings of this study strongly support the molecular breeding of salt-tolerant wheat, particularly by utilizing the genetic advancements based on CM6005 and KN9204.

Key message: This study not only provide a theoretical basis for screening drought-resistant yellowhorn rootstocks but also lay the foundation for the research on the drought resistance mechanism of yellowhorn. Drought is a major limiting factor in the production and cultivation of yellowhorn, an important woody oil species in China. It is important to select excellent drought-resistant rootstocks, and the current reproduction of rootstocks of yellowhorn is mainly through seed propagation. In this study, we analyzed the phenotypes, physiology and anatomy of the F1 generation of five yellowhorn varieties ('Zhongshi 1', 'Zhongshi 4', 'Zhongshi 7', 'Zhongshi 9' and 'Yuanda') which are currently the main promoted varieties of yellowhorn. The changes in physiological traits were used to comprehensively rank drought resistance by the membership function method. The results showed that the drought resistance of different varieties was in the order 'Zhongshi 9' > 'Zhongshi 4' > 'Yuanda' > 'Zhongshi 7' > 'Zhongshi 1'. In addition, we investigated the molecular signature of two yellowhorn varieties exhibiting a large difference in levels of drought resistance using transcriptome data. A total of 119 differentially expressed genes (DEGs) were found between the 2 yellowhorn varieties under drought stress. Based on KEGG analysis, DEGs were mostly enriched in three pathways. According to the correlation analysis between weighted gene co-expression network analysis (WGCNA) and physiological traits, a total of four co-expression networks with high correlation were constructed, and some hub genes of yellowhorn in response to drought stress were found. We carried out a comprehensive evaluation of yellowhorn drought tolerance based on parental lines. It solves the current problems faced in the breeding of yellowhorn and lays the foundation for studying the molecular mechanism of drought resistance in yellowhorn. [ABSTRACT FROM AUTHOR]

Key message: We constructed a gene expression atlas and co-expression network for potatoes and identified several novel genes associated with various agronomic traits. This resource will accelerate potato genetics and genomics research. Potato (Solanum tuberosum L.) is the world's most crucial non-cereal food crop and ranks third in food production after wheat and rice. Despite the availability of several potato transcriptome datasets at public databases like NCBI SRA, an effort has yet to be put into developing a global transcriptome atlas and a co-expression network for potatoes. The objectives of our study were to construct a global expression atlas for potatoes using publicly available transcriptome datasets, identify housekeeping and tissue-specific genes, construct a global co-expression network and identify co-expression clusters, investigate the transcriptional complexity of genes involved in various essential biological processes related to agronomic traits, and provide a web server (StCoExpNet) to easily access the newly constructed expression atlas and co-expression network to investigate the expression and co-expression of genes of interest. In this study, we used data from 2299 publicly available potato transcriptome samples obtained from 15 different tissues to construct a global transcriptome atlas. We found that roughly 87% of the annotated genes exhibited detectable expression in at least one sample. Among these, we identified 281 genes with consistent and stable expression levels, indicating their role as housekeeping genes. Conversely, 308 genes exhibited marked tissue-specific expression patterns. We exemplarily linked some co-expression clusters to important agronomic traits of potatoes, such as self-incompatibility, anthocyanin biosynthesis, tuberization, and defense responses against multiple pathogens. The dataset compiled here constitutes a new resource (StCoExpNet), which can be accessed at . This transcriptome atlas and the co-expression network will accelerate potato genetics and genomics research. [ABSTRACT FROM AUTHOR]

After the loss of a trait, theory predicts that the molecular machinery underlying its phenotypic expression should decay. Yet, empirical evidence is contrasting. Here, we test the hypotheses that (i) the molecular ground plan of a lost trait could persist due to pleiotropic effects on other traits and (ii) that gene co-expression network architecture could constrain individual gene expression. Our testing ground has been the Bacillus stick insect species complex, which contains close relatives that are either bisexual or parthenogenetic. After the identification of genes expressed in male reproductive tissues in a bisexual species, we investigated their gene co-expression network structure in two parthenogenetic species. We found that gene co-expression within the male gonads was partially preserved in parthenogens. Furthermore, parthenogens did not show relaxed selection on genes upregulated in male gonads in the bisexual species. As these genes were mostly expressed in female gonads, this preservation could be driven by pleiotropic interactions and an ongoing role in female reproduction. Connectivity within the network also played a key role, with highly connected—and more pleiotropic—genes within male gonad also having a gonad-biased expression in parthenogens. Our findings provide novel insight into the mechanisms which could underlie the production of rare males in parthenogenetic lineages; more generally, they provide an example of the cryptic persistence of a lost trait molecular architecture, driven by gene pleiotropy on other traits and within their co-expression network. [ABSTRACT FROM AUTHOR]

Background: Gastric cancer (GC) is a significant health issue globally, ranking as the fifth most common cancer with over 10,000 new cases reported annually. Long non-coding RNA (lncRNA) has emerged as a critical player in cellular functions, influencing GC's development, growth, metastasis, and prognosis. However, our understanding of lncRNA's role in the pathogenesis of GC remains limited. Therefore, it is particularly important to explore the relationship between lncRNA and gastric cancer. Methods: we conducted a comprehensive analysis of RNA sequencing data from the GEO database and stomach adenocarcinoma (STAD) data from the TCGA database to identify lncRNAs that exhibit altered expression levels in GC and the mechanisms underlying lncRNA-mediated transcription and post-transcriptional regulation were explored. Results: This study uncovered 94 lncRNAs with differential expression and, through co-expression analysis, linked these to 1508 differentially expressed genes (DEGs). GO functional enrichment analysis highlighted that these DEGs are involved in critical pathways, such as cell adhesion and the positive regulation of cell migration. By establishing a lncRNA-miRNA-mRNA regulatory network, we found that the ceRNA mechanism, particularly involving RP11-357H14.17 and CTD-2377D24.4, could play a role in GC progression. Experimental validation of selected differentially expressed lncRNAs and mRNAs (including RP11-357H14.17-CLDN1, BBOX1, TRPM2-AS, CLDN1, PLAU, HOXB7) confirmed the RNA-seq results. Conclusions: Overall, our findings highlight the critical role of the lncRNA-mRNA regulatory network in the development and progression of GC, offering potential biomarkers for diagnosis and targets for innovative treatment strategies. [ABSTRACT FROM AUTHOR]

Background: Soil salinization is one of the vital factors threatening the world's food security. To reveal the biological mechanism of response to salt stress in wheat, this study was conducted to resolve the transcription level difference to salt stress between CM6005 (salt-tolerant) and KN9204 (salt-sensitive) at the germination and seedling stage. Results: To investigate the molecular mechanism underlying salt tolerance in wheat, we conducted comprehensive transcriptome analyses at the seedling and germination stages. Two wheat cultivars, CM6005 (salt-tolerant) and KN9204 (salt-sensitive) were subjected to salt treatment, resulting in a total of 24 transcriptomes. Through expression-network analysis, we identified 17 modules, 16 and 13 of which highly correlate with salt tolerance-related phenotypes in the germination and seedling stages, respectively. Moreover, we identified candidate Hub genes associated with specific modules and explored their regulatory relationships using co-expression data. Enrichment analysis revealed specific enrichment of gibberellin-related terms and pathways in CM6005, highlighting the potential importance of gibberellin regulation in enhancing salt tolerance. In contrast, KN9204 exhibited specific enrichment in glutathione-related terms and activities, suggesting the involvement of glutathione-mediated antioxidant mechanisms in conferring resistance to salt stress. Additionally, glucose transport was found to be a fundamental mechanism for salt tolerance during wheat seedling and germination stages, indicating its potential universality in wheat. Wheat plants improve their resilience and productivity by utilizing adaptive mechanisms like adjusting osmotic balance, bolstering antioxidant defenses, accumulating compatible solutes, altering root morphology, and regulating hormones, enabling them to better withstand extended periods of salt stress. Conclusion: Through utilizing transcriptome-level analysis employing WGCNA, we have revealed a potential regulatory mechanism that governs the response to salt stress and recovery in wheat cultivars. Furthermore, we have identified key candidate central genes that play a crucial role in this mechanism. These central genes are likely to be vital components within the gene expression network associated with salt tolerance. The findings of this study strongly support the molecular breeding of salt-tolerant wheat, particularly by utilizing the genetic advancements based on CM6005 and KN9204. [ABSTRACT FROM AUTHOR]

Nitrogen is an essential nutrient for plants and a major determinant of plant growth and crop yield. Plants acquire nitrogen mainly in the form of nitrate and ammonium. Both nitrogen sources affect plant responses and signaling pathways in a different way, but these signaling pathways interact, complicating the study of nitrogen responses. Extensive transcriptome analyses and the construction of gene regulatory networks, mainly in response to nitrate, have significantly advanced our understanding of nitrogen signaling and responses in model plants and crops. In this study, we aimed to generate a more comprehensive gene regulatory network for the major crop, rice, by incorporating the interactions between ammonium and nitrate. To achieve this, we assessed transcriptome changes in rice roots and shoots over an extensive time course under single or combined applications of the two nitrogen sources. This dataset enabled us to construct a holistic co-expression network and identify potential key regulators of nitrogen responses. Next to known transcription factors, we identified multiple new candidates, including the transcription factors OsRLI and OsEIL1, which we demonstrated to induce the primary nitrate-responsive genes OsNRT1.1b and OsNIR1. Our network thus serves as a valuable resource to obtain novel insights in nitrogen signaling.

IntroductionHemibiotrophic Phytophthora are a group of agriculturally and ecologically important pathogenic oomycetes causing severe decline in plant growth and fitness. The lifestyle of these pathogens consists of an initial biotrophic phase followed by a switch to a necrotrophic phase in the latter stages of infection. Between these two phases is the biotrophic to necrotrophic switch (BNS) phase, the timing and controls of which are not well understood particularly in Phytophthora spp. where host resistance has a purely quantitative genetic basis.MethodsTo investigate this we sequenced and annotated the genome of Phytophthora medicaginis, causal agent of root rot and substantial yield losses to Fabaceae hosts. We analyzed the transcriptome of P. medicaginis across three phases of colonization of a susceptible chickpea host (Cicer arietinum) and performed co-regulatory analysis to identify putative small secreted protein (SSP) effectors that influence timing of the BNS in a quantitative pathosystem.ResultsThe genome of P. medicaginis is ~78 Mb, comparable to P. fragariae and P. rubi which also cause root rot. Despite this, it encodes the second smallest number of RxLR (arginine-any amino acid-leucine-arginine) containing proteins of currently sequenced Phytophthora species. Only quantitative resistance is known in chickpea to P. medicaginis, however, we found that many RxLR, Crinkler (CRN), and Nep1-like protein (NLP) proteins and carbohydrate active enzymes (CAZymes) were regulated during infection. Characterization of one of these, Phytmed_10271, which encodes an RxLR effector demonstrates that it plays a role in the timing of the BNS phase and root cell death.DiscussionThese findings provide an important framework and resource for understanding the role of pathogenicity factors in purely quantitative Phytophthora pathosystems and their implications to the timing of the BNS phase.

Abstract Heart transplantation is the gold standard for treating patients with advanced heart failure. Although improvements in immunosuppressive therapies have significantly reduced the frequency of cardiac graft rejection, the incidences of T cell-mediated rejection (TCMR) and antibody-mediated rejection remain almost unchanged. A four-archetype analysis (4AA) model, developed by Philip F. Halloran, illustrated this problem well. It provided a new dimension to improve the accuracy of diagnoses and an independent system for recalibrating the histology guidelines. However, this model was based on the invasive method of endocardial biopsy, which undoubtedly increased the postoperative risk of heart transplant patients. Currently, little is known regarding the associated genes and specific functions of the different phenotypes. We performed bioinformatics analysis (using machine-learning methods and the WGCNA algorithm) to screen for hub-specific genes related to different phenotypes, based Gene Expression Omnibus accession number GSE124897. More immune cell infiltration was observed with the ABMR, TCMR, and injury phenotypes than with the stable phenotype. Hub-specific genes for each of the four archetypes were verified successfully using an external test set (accession number GSE2596). Logistic-regression models based on TCMR-specific hub genes and common hub genes were constructed with accurate diagnostic utility (area under the curve > 0.95). RELA, NFKB1, and SOX14 were identified as transcription factors important for TCMR/injury phenotypes and common genes, respectively. Additionally, 11 Food and Drug Administration-approved drugs were chosen from the DrugBank Database for each four-archetype model. Tyrosine kinase inhibitors may be a promising new option for transplant rejection treatment. KRAS signaling in cardiac transplant rejection is worth further investigation. Our results showed that heart transplant rejection subtypes can be accurately diagnosed by detecting expression of the corresponding specific genes, thereby enabling precise treatment or medication.

The double flower developmental process is regulated via a complex transcriptional regulatory network. To understand this highly dynamic and complex developmental process of Dianthus spp., we performed a comparative analysis of floral morphology and transcriptome dynamics in simple flowers and double flowers. We found that the primordium of double flowers of ‘X’ was larger in size compared to that of simple flowers of ‘L’ in Dianthus chinensis. RNA-seq and Weighted Gene Co-expression Network Analysis (WGCNA) during flower development, identified stage-specific gene network modules. Expression analysis by RNA-seq indicated that a group of genes related to floral meristem identity, primordia position and polarity were highly expressed in double flowers genotypes compared to simple flowers genotypes, suggesting their roles in double-petal formation. A total of 21 DEGs related to petal number were identified between simple and double flowers. The experiments of in situ hybridization revealed that DcaAP2L, DcaLFY and DcaUFO genes were expressed in the intra-sepal boundary and petal boundary. We proposed a potential transcriptional regulatory network for simple and double flower development. This study provides novel insights into the molecular mechanism underlying double flower formation in Dianthus spp.

Myricaria laxiflora is an endangered shrub plant with remarkable tolerance to waterlogging stress, however, little attention has been paid to understanding the underlying mechanisms. Here, physiological and transcriptomic approaches were applied to uncover the physiological and molecular reconfigurations in the stem of M. laxiflora in response to waterlogging stress. The accumulation of the contents of H2O2 and malonaldehyde (MDA) alongside increased activities of enzymes for scavenging the reactive oxygen species (ROS) in the stem of M. laxiflora were observed under waterlogging stress. The principal component analysis (PCA) of transcriptomes from five different timepoints uncovered PC1 counted for 17.3 % of total variations and separated the treated and non-treated samples. A total of 8714 genes in the stem of M. laxiflora were identified as differentially expressed genes (DEGs) under waterlogging stress, which could be assigned into two different subgroups with distinct gene expression patterns and biological functions. The DEGs involved in glycolysis were generally upregulated, whereas opposite results were observed for nitrogen uptake and the assimilation pathway. The contents of abscisic acid (ABA) and jasmonic acid (JA) were sharply decreased alongside the decreased mRNA levels of the genes involved in corresponding synthesis pathways upon waterlogging stress. A network centered by eight key transcription factors has been constructed, which uncovered the inhibited cell division processes in the stem of M. laxiflora upon waterlogging stress. Taken together, the obtained results showed that glycolysis, nitrogen metabolism and meristem activities played an important role in the stem of M. laxiflora in response to waterlogging stress.

Heart transplantation is the gold standard for treating patients with advanced heart failure. Although improvements in immunosuppressive therapies have significantly reduced the frequency of cardiac graft rejection, the incidences of T cell-mediated rejection (TCMR) and antibody-mediated rejection remain almost unchanged. A four-archetype analysis (4AA) model, developed by Philip F. Halloran, illustrated this problem well. It provided a new dimension to improve the accuracy of diagnoses and an independent system for recalibrating the histology guidelines. However, this model was based on the invasive method of endocardial biopsy, which undoubtedly increased the postoperative risk of heart transplant patients. Currently, little is known regarding the associated genes and specific functions of the different phenotypes. We performed bioinformatics analysis (using machine-learning methods and the WGCNA algorithm) to screen for hub-specific genes related to different phenotypes, based Gene Expression Omnibus accession number GSE124897. More immune cell infiltration was observed with the ABMR, TCMR, and injury phenotypes than with the stable phenotype. Hub-specific genes for each of the four archetypes were verified successfully using an external test set (accession number GSE2596). Logistic-regression models based on TCMR-specific hub genes and common hub genes were constructed with accurate diagnostic utility (area under the curve > 0.95). RELA, NFKB1, and SOX14 were identified as transcription factors important for TCMR/injury phenotypes and common genes, respectively. Additionally, 11 Food and Drug Administration-approved drugs were chosen from the DrugBank Database for each four-archetype model. Tyrosine kinase inhibitors may be a promising new option for transplant rejection treatment. KRAS signaling in cardiac transplant rejection is worth further investigation. Our results showed that heart transplant rejection subtypes can be accurately diagnosed by detecting expression of the corresponding specific genes, thereby enabling precise treatment or medication. [ABSTRACT FROM AUTHOR]

WUSCHEL-related homeobox (WOX) genes are a class of plant-specific transcription factors, regulating the development of multiple tissues. However, the genomic characterizations and expression patterns of WOX genes have not been analyzed in lotus. In this study, 15 NnWOX genes were identified based on the well-annotated reference genome of lotus. According to the phylogenetic analysis, the NnWOX genes were clustered into three clades, i.e., ancient clade, intermediate clade, and WUS clade. Except for the conserved homeobox motif, we further found specific motifs of NnWOX genes in different clades and divergence gene structures, suggesting their distinct functions. In addition, two NnWOX genes in the ancient clade have conserved expression patterns and other NnWOX genes exhibit different expression patterns in lotus tissues, suggesting a low level of functional redundancy in lotus WOX genes. Furthermore, we constructed the gene co-expression networks for each NnWOX gene. Based on weighted gene co-expression network analysis (WGCNA), ten NnWOX genes and their co-expressed genes were assigned to the modules that were significantly related to the cotyledon and seed coat. We further performed RT-qPCR experiments, validating the expression levels of ten NnWOX genes in the co-expression networks. Our study reveals comprehensive genomic features of NnWOX genes in lotus, providing a solid basis for further function studies. [ABSTRACT FROM AUTHOR]

Introduction: Fruit cracking not only affects the appearance of netted melons (Cucumis melo L. var. reticulatus Naud.) but also decreases their marketability. Methods: Herein, to comprehensively understand the role of expansin (EXP) proteins in netted melon, bioinformatics methods were employed to discover the EXP gene family in the melon genome and analyze its characteristic features. Furthermore, transcriptomics analysis was performed to determine the expression patterns of melon EXP (CmEXP) genes in crack-tolerant and cracksusceptible netted melon varieties. Discussion: Thirty-three CmEXP genes were identified. Chromosomal location analysis revealed that CmEXP gene distribution was uneven on 12 chromosomes. In addition, phylogenetic tree analysis revealed that CmEXP genes could be categorized into four subgroups, among which the EXPA subgroup had the most members. The same subgroup members shared similar protein motifs and gene structures. Thirteen duplicate events were identified in the 33 CmEXP genes. Collinearity analysis revealed that the CmEXP genes had 50, 50, and 44 orthologous genes with EXP genes in cucumber, watermelon, and Arabidopsis, respectively. However, only nine orthologous EXP genes were observed in rice. Promoter cis-acting element analysis demonstrated that numerous cis-acting elements in the upstream promoter region of CmEXP genes participate in plant growth, development, and environmental stress responses. Transcriptomics analysis revealed 14 differentially expressed genes (DEGs) in the non-cracked fruit peels between the crack-tolerant variety 'Xizhoumi 17' (N17) and the cracksusceptible variety 'Xizhoumi 25' (N25). Among the 14 genes, 11 were upregulated, whereas the remaining three were downregulated in N17. In the non-cracked (N25) and cracked (C25) fruit peels of 'Xizhoumi 25', 24DEGswere identified, and 4 of them were upregulated, whereas the remaining 20 were downregulated in N25. In the two datasets, only CmEXPB1 exhibited consistently upregulated expression, indicating its importance in the fruit peel crack resistance of netted melon. Transcription factor prediction revealed 56 potential transcription factors that regulate CmEXPB1 expression. Results: Our study findings enrich the understanding of the CmEXP gene family and present candidate genes for the molecular breeding of fruit peel crack resistance of netted melon. [ABSTRACT FROM AUTHOR]

Phoebe bournei is one of the main afforestation tree species in subtropical regions of China and is famous for its timber. Its distribution and growth are significantly impaired by water conditions. Thus, it is essential to understand the mechanism of the stress response in P. bournei. Here, we analyzed the phenotypic changes and transcriptomic rearrangement in the leaves and roots of P. bournei seedlings grown for 0 h, 1 h, 24 h, and 72 h under simulated drought conditions (10% PEG 6000). The results showed that drought stress inhibited plant photosynthesis and increased oxidoreductase activity and abscisic acid (ABA) accumulation. Spatio-temporal transcriptomic analysis identified 2836 and 3704 differentially expressed genes (DEGs) in leaves and roots, respectively. The responsive genes in different organs presented various expression profiles at different times. Gene co-expression network analysis identified two core transcription factors, TGA4 and APRR2, from two modules that showed a strong positive correlation with ABA accumulation. Our study investigated the different responses of aboveground and belowground organs of P. bournei to drought stress and provides critical information for improving the drought resistance of this timber species. [ABSTRACT FROM AUTHOR]

DUSP4 is a member of the DUSP (dual-specificity phosphatase) subfamily that is selective to the mitogen-activated protein kinases (MAPK) and has been implicated in a range of biological processes and functions in Alzheimer's disease (AD). In this study, we utilized the stereotactic delivery of adeno-associated virus (AAV)-DUSP4 to overexpress DUSP4 in the dorsal hippocampus of 5xFAD and wildtype (WT) mice, then used mass spectrometry (MS)-based proteomics along with the label-free quantification to profile the proteome and phosphoproteome in the hippocampus. We identified protein expression and phosphorylation patterns modulated in 5xFAD mice and examined the sex-specific impact of DUSP4 overexpression on the 5xFAD proteome/phosphoproteome. In 5xFAD mice, a substantial number of proteins were up- or down-regulated in both male and female mice in comparison to age and sex-matched WT mice, many of which are involved in AD-related biological processes, such as activated immune response or suppressed synaptic activities. Many proteins in pathways, such as immune response were found to be suppressed in response to DUSP4 overexpression in male 5xFAD mice. In contrast, such a shift was absent in female mice. For the phosphoproteome, we detected an array of phosphorylation sites regulated in 5xFAD compared to WT and modulated via DUSP4 overexpression in each sex. Interestingly, 5xFAD- and DUSP4-associated phosphorylation changes occurred in opposite directions. Strikingly, both the 5xFAD- and DUSP4-associated phosphorylation changes were found to be mostly in neurons and play key roles in neuronal processes and synaptic functions. Site-centric pathway analysis revealed that both the 5xFAD- and DUSP4-associated phosphorylation sites were enriched for a number of kinase sets in females but only a limited number of sets of kinases in male mice. Taken together, our results suggest that male and female 5xFAD mice responded to DUSP4 overexpression via shared and sex-specific molecular mechanisms, which might underly similar reductions in amyloid pathology in both sexes while learning deficits were reduced in only females with DUSP4 overexpression. Finally, we validated our findings with the sex-specific AD-associated proteomes in human cohorts and further developed DUSP4-centric proteomic network models and signaling maps for each sex. [ABSTRACT FROM AUTHOR]

Abstract Background The origin recognition complex 1 (ORC1) is a large subunit of the origin recognition complex and acts as the master subunit of the precoding complex. Objective To explore potential function and clinical significance of ORC1 in cancers. Methods The expression level of ORC1 in different types of tumor tissues and matched normal tissues were detected by The Cancer Genome Atlas (TCGA) and validated by datasets from the gene expression omnibus (GEO) database. The association between ORC1 expression and infiltration levels of immune cell was analyzed. ORC1 and its co-expression genes were subjected to enrichment analysis to explore potential mechanisms in cancers, and the protein-protein interaction (PPI) network was constructed. Finally, the expression of ORC1 in tumor tissue and adjacent tissue was verified by immunohistochemistry (IHC). Results ORC1 was highly expressed in the majority of tumors, and the expression level of ORC1 was associated with the pathological stages of ACC, LUAD, OV and SKCM. ORC1 was closely related with poor prognosis in ACC, LIHC, PAAD, READ and THCA. ORC1 in ACC and KICH was positively correlated with the infiltration level of immune cells while it was negatively correlated with the infiltration level of immune cells in THYM. Co-expression network analysis showed that CDCA3, GSG2, KIF2C, NCAPH and PLK1 were positively correlated with ORC1 in cancer, and enrichment analysis showed a correlation with cytosol, ATP binding and cell division. The expression of ORC1 in UCEC and KICH was higher than that in the adjacent tissues. Conclusion ORC1 over-expressed in most tumors and could be severed as a novel biomarker for diagnosis. This study revealed that ORC1 might inhibit tumor immunity and might be a potential therapeutic target in cancers.

One of the most destructive diseases, Gibberella stalk rot (GSR), caused by Fusarium graminearum, reduces maize yields significantly. An induced resistance response is a potent and cost-effective plant defense against pathogen attack. The functional counterpart of JAs, coronatine (COR), has attracted a lot of interest recently due to its ability to control plant growth and stimulate secondary metabolism. Although several studies have focused on COR as a plant immune elicitor to improve plant resistance to pathogens, the effectiveness and underlying mechanisms of the suppressive ability against COR to F. graminearum in maize have been limited. We investigated the potential physiological and molecular mechanisms of COR in modulating maize resistance to F. graminearum. COR treatment strongly enhanced disease resistance and promoted stomatal closure with H2O2 accumulation, and 10 μg/mL was confirmed as the best concentration. COR treatment increased defense-related enzyme activity and decreased the malondialdehyde content with enhanced antioxidant enzyme activity. To identify candidate resistance genes and gain insight into the molecular mechanism of GSR resistance associated with COR, we integrated transcriptomic and metabolomic data to systemically explore the defense mechanisms of COR, and multiple hub genes were pinpointed using weighted gene correlation network analysis (WGCNA). We discovered 6 significant modules containing 10 candidate genes: WRKY transcription factor (LOC100279570), calcium-binding protein (LOC100382070), NBR1-like protein (LOC100275089), amino acid permease (LOC100382244), glutathione S-transferase (LOC541830), HXXXD-type acyl-transferase (LOC100191608), prolin-rich extensin-like receptor protein kinase (LOC100501564), AP2-like ethylene-responsive transcription factor (LOC100384380), basic leucine zipper (LOC100275351), and glycosyltransferase (LOC606486), which are highly correlated with the jasmonic acid–ethylene signaling pathway and antioxidants. In addition, a core set of metabolites, including alpha-linolenic acid metabolism and flavonoids biosynthesis linked to the hub genes, were identified. Taken together, our research revealed differentially expressed key genes and metabolites, as well as co-expression networks, associated with COR treatment of maize stems after F. graminearum infection. In addition, COR-treated maize had higher JA (JA-Ile and Me-JA) levels. We postulated that COR plays a positive role in maize resistance to F. graminearum by regulating antioxidant levels and the JA signaling pathway, and the flavonoid biosynthesis pathway is also involved in the resistance response against GSR. [ABSTRACT FROM AUTHOR]

Aquaporins play a crucial role in helping water molecules move across cell membranes. While some studies have examined the role of AQPs in model plants like Arabidopsis, their impact on the ability of non-model plants to withstand environmental stress is largely unknown. In this study, we have explored the functions of the PaPIP1-2 gene, which encodes a protein called PIP, in apricot kernels (Prunus armeniaca L.). Our findings reveal how the PaPIP1-2 gene behaves during both dormancy and sprouting phases. Using a network analysis, we identified its interaction with six genes related to cold resistance. The mRNA levels of PaAQP genes, which co-express with cold resistance genes, remain consistent throughout different stages of P. armeniaca flower bud development, including physiological dormancy (PD), ecological dormancy (ED), sprouting period (SP), and germination stage (GS). Furthermore, our investigation of the location of the GFP-tagged PaPIP1-2 protein showed that it is mainly found in the cell membrane. Yeast strains with overexpressed PaPIP1-2 exhibited improved cold resistance and higher protein content. Similarly, when we overexpressed PaPIP1-2 in Arabidopsis, it enhanced the growth of these transgenic plants under cold stress. This improvement was associated with reduced levels of MDA (malondialdehyde); decreased ion leakage; increased proline accumulation; superoxide dismutase (SOD) activity; and the expression of cold resistance genes like AtPUB26, AtBTF3L, AtEBF1-1, and AtRAV1, compared with the wild-type plants. In summary, our results highlight the role of the P. armeniaca PaPIP1-2 gene in enhancing cold resistance and its importance in the dormancy and germination stages. [ABSTRACT FROM AUTHOR]

The induction of mutants using a variety of mutagens provides new genetic resources for crop breeding. In the current study, we generated mutants in sorghum using gamma-ray treatment and conducted transcriptome analyses between a mutant line and its origin. Based on the data generated by the analysis, we predicted genes that were differentially expressed in the mutant compared to its origin and identified their functions. Gene Ontology and KEGG pathway were performed to annotate DEGs functions. Furthermore, we constructed gene co-expression networks to demonstrate the association of the identified transcription factors. We also searched for single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDel) which can possibly cause the alteration of gene functions. The SNP data discovered in this transcript analysis found 11 genes that appeared in common with co-expression networks and examined the associations for salt-responsiveness. This study provides genetic backgrounds that can emerge when sorghum is subjected to gamma irradiation to progress the mutant. [ABSTRACT FROM AUTHOR]

Ear traits are key contributors to grain yield in maize; therefore, exploring their genetic basis facilitates the improvement of grain yield. However, the underlying molecular mechanisms of ear traits remain obscure in both inbred lines and hybrids. Here, two association panels, respectively, comprising 203 inbred lines (IP) and 246 F1 hybrids (HP) were employed to identify candidate genes for six ear traits. The IP showed higher phenotypic variation and lower phenotypic mean than the HP for all traits, except ear tip-barrenness length. By conducting a genome-wide association study (GWAS) across multiple environments, 101 and 228 significant single-nucleotide polymorphisms (SNPs) associated with six ear traits were identified in the IP and HP, respectively. Of these significant SNPs identified in the HP, most showed complete–incomplete dominance and over-dominance effects for each ear trait. Combining a gene co-expression network with GWAS results, 186 and 440 candidate genes were predicted in the IP and HP, respectively, including known ear development genes ids1 and sid1. Of these, nine candidate genes were detected in both populations and expressed in maize ear and spikelet tissues. Furthermore, two key shared genes (GRMZM2G143330 and GRMZM2G171139) in both populations were found to be significantly associated with ear traits in the maize Goodman diversity panel with high-density variations. These findings advance our knowledge of the genetic architecture of ear traits between inbred lines and hybrids and provide a valuable resource for the genetic improvement of ear traits in maize. [ABSTRACT FROM AUTHOR]

Introduction: Alzheimer's disease (AD) is heterogeneous, both clinically and neuropathologically. We investigated whether polygenic risk scores (PRSs) integrated with transcriptome profiles from AD brains can explain AD clinical heterogeneity. Methods: We conducted co‐expression network analysis and identified gene sets (modules) that were preserved in three AD transcriptome datasets and associated with AD‐related neuropathological traits including neuritic plaques (NPs) and neurofibrillary tangles (NFTs). We computed the module‐based PRSs (mbPRSs) for each module and tested associations with mbPRSs for cognitive test scores, cognitively defined AD subgroups, and brain imaging data. Results: Of the modules significantly associated with NPs and/or NFTs, the mbPRSs from two modules (M6 and M9) showed distinct associations with language and visuospatial functioning, respectively. They matched clinical subtypes and brain atrophy at specific regions. Discussion: Our findings demonstrate that polygenic profiling based on co‐expressed gene sets can explain heterogeneity in AD patients, enabling genetically informed patient stratification and precision medicine in AD. HIGHLIGHTS: Co‐expression gene‐network analysis in Alzheimer's disease (AD) brains identified gene sets (modules) associated with AD heterogeneity.AD‐associated modules were selected when genes in each module were enriched for neuritic plaques and neurofibrillary tangles.Polygenic risk scores from two selected modules were linked to the matching cognitively defined AD subgroups (language and visuospatial subgroups).Polygenic risk scores from the two modules were associated with cognitive performance in language and visuospatial domains and the associations were confirmed in regional‐specific brain atrophy data. [ABSTRACT FROM AUTHOR]

Background: The origin recognition complex 1 (ORC1) is a large subunit of the origin recognition complex and acts as the master subunit of the precoding complex. Objective: To explore potential function and clinical significance of ORC1 in cancers. Methods: The expression level of ORC1 in different types of tumor tissues and matched normal tissues were detected by The Cancer Genome Atlas (TCGA) and validated by datasets from the gene expression omnibus (GEO) database. The association between ORC1 expression and infiltration levels of immune cell was analyzed. ORC1 and its co-expression genes were subjected to enrichment analysis to explore potential mechanisms in cancers, and the protein-protein interaction (PPI) network was constructed. Finally, the expression of ORC1 in tumor tissue and adjacent tissue was verified by immunohistochemistry (IHC). Results: ORC1 was highly expressed in the majority of tumors, and the expression level of ORC1 was associated with the pathological stages of ACC, LUAD, OV and SKCM. ORC1 was closely related with poor prognosis in ACC, LIHC, PAAD, READ and THCA. ORC1 in ACC and KICH was positively correlated with the infiltration level of immune cells while it was negatively correlated with the infiltration level of immune cells in THYM. Co-expression network analysis showed that CDCA3, GSG2, KIF2C, NCAPH and PLK1 were positively correlated with ORC1 in cancer, and enrichment analysis showed a correlation with cytosol, ATP binding and cell division. The expression of ORC1 in UCEC and KICH was higher than that in the adjacent tissues. Conclusion: ORC1 over-expressed in most tumors and could be severed as a novel biomarker for diagnosis. This study revealed that ORC1 might inhibit tumor immunity and might be a potential therapeutic target in cancers. [ABSTRACT FROM AUTHOR]

BackgroundWith the advances in the high throughput next generation sequencing technologies, genome-wide association studies (GWAS) have identified a large set of variants associated with complex phenotypic traits at a very fine scale. Despite the progress in GWAS, identification of genotype-phenotype relationship remains challenging in maize due to its nature with dozens of variants controlling the same trait. As the causal variations results in the change in expression, gene expression analyses carry a pivotal role in unraveling the transcriptional regulatory mechanisms behind the phenotypes.ResultsTo address these challenges, we incorporated the gene expression and GWAS-driven traits to extend the knowledge of genotype-phenotype relationships and transcriptional regulatory mechanisms behind the phenotypes. We constructed a large collection of gene co-expression networks and identified more than 2 million co-expressing gene pairs in the GWAS-driven pan-network which contains all the gene-pairs in individual genomes of the nested association mapping (NAM) population. We defined four sub-categories for the pan-network: (1) core-network contains the highest represented ~ 1% of the gene-pairs, (2) near-core network contains the next highest represented 1-5% of the gene-pairs, (3) private-network contains ~ 50% of the gene pairs that are unique to individual genomes, and (4) the dispensable-network contains the remaining 50-95% of the gene-pairs in the maize pan-genome. Strikingly, the private-network contained almost all the genes in the pan-network but lacked half of the interactions. We performed gene ontology (GO) enrichment analysis for the pan-, core-, and private- networks and compared the contributions of variants overlapping with genes and promoters to the GWAS-driven pan-network.ConclusionsGene co-expression networks revealed meaningful information about groups of co-regulated genes that play a central role in regulatory processes. Pan-network approach enabled us to visualize the global view of the gene regulatory network for the studied system that could not be well inferred by the core-network alone.

Abstract Background Song performed in flocks by European starlings (Sturnus vulgaris), referred to here as gregarious song, is a non-sexual, social behavior performed by adult birds. Gregarious song is thought to be an intrinsically reinforced behavior facilitated by a low-stress, positive affective state that increases social cohesion within a flock. The medial preoptic area (mPOA) is a region known to have a role in the production of gregarious song. However, the neurochemical systems that potentially act within this region to regulate song remain largely unexplored. In this study, we used RNA sequencing to characterize patterns of gene expression in the mPOA of male and female starlings singing gregarious song to identify possibly novel neurotransmitter, neuromodulator, and hormonal pathways that may be involved in the production of gregarious song. Results Differential gene expression analysis and rank rank hypergeometric analysis indicated that dopaminergic, cholinergic, and GABAergic systems were associated with the production of gregarious song, with multiple receptor genes (e.g., DRD2, DRD5, CHRM4, GABRD) upregulated in the mPOA of starlings who sang at high rates. Additionally, co-expression network analyses identified co-expressing gene clusters of glutamate signaling-related genes associated with song. One of these clusters contained five glutamate receptor genes and two glutamate scaffolding genes and was significantly enriched for genetic pathways involved in neurodevelopmental disorders associated with social deficits in humans. Two of these genes, GRIN1 and SHANK2, were positively correlated with performance of gregarious song. Conclusions This work provides new insights into the role of the mPOA in non-sexual, gregarious song in starlings and highlights candidate genes that may play a role in gregarious social interactions across vertebrates. The provided data will also allow other researchers to compare across species to identify conserved systems that regulate social behavior.

Abstract Sepsis is a life-threatening condition characterized by a harmful host response to infection with organ dysfunction. Annually about 20 million people are dead owing to sepsis and its mortality rates is as high as 20%. However, no studies have been carried out to investigate sepsis from the system biology point of view, as previous research predominantly focused on individual genes without considering their interactions and associations. Here, we conducted a comprehensive exploration of genome-wide expression alterations in both mRNAs and long non-coding RNAs (lncRNAs) in sepsis, using six microarray datasets. Co-expression networks were conducted to identify mRNA and lncRNA modules, respectively. Comparing these sepsis modules with normal modules, we observed a homogeneous expression pattern within the mRNA/lncRNA members, with the majority of them displaying consistent expression direction. Moreover, we identified consistent modules across diverse datasets, consisting of 20 common mRNA members and two lncRNAs, namely CHRM3-AS2 and PRKCQ-AS1, which are potential regulators of sepsis. Our results reveal that the up-regulated common mRNAs are mainly involved in the processes of neutrophil mediated immunity, while the down-regulated mRNAs and lncRNAs are significantly overrepresented in T-cell mediated immunity functions. This study sheds light on the co-expression patterns of mRNAs and lncRNAs in sepsis, providing a novel perspective and insight into the sepsis transcriptome, which may facilitate the exploration of candidate therapeutic targets and molecular biomarkers for sepsis.

Abstract Background Lignin is a major restriction factor for the industrial production of biomass resources, such as pulp and bioenergy. Eucalyptus is one of the most important sources of pulp and bioenergy. After polyploidization, the lignin content of forest trees is generally reduced, which is considered a beneficial genetic improvement. However, the differences in the lignin content between triploid and diploid Eucalyptus and the underlying regulatory mechanism are still unclear. Results We conducted a comprehensive analysis at the phenotypic, transcriptional and metabolite levels between Eucalyptus urophylla triploids and diploids to reveal the effects of polyploidization on the lignin content and lignin metabolic pathway. The results showed that the lignin content of Eucalyptus urophylla triploid stems was significantly lower than that of diploids. Lignin-related metabolites were differentially accumulated between triploids and diploids, among which coniferaldehyde, p-coumaryl alcohol, sinapaldehyde and coniferyl alcohol had significant positive correlations with lignin content, indicating that they might be primarily contributing metabolites. Most lignin biosynthetic genes were significantly downregulated, among which 11 genes were significantly positively correlated with the lignin content and above metabolites. Furthermore, we constructed a co-expression network between lignin biosynthetic genes and transcription factors based on weighted gene co-expression network analysis. The network identified some putative orthologues of secondary cell wall (SCW)-related transcription factors, among which MYB52, MYB42, NAC076, and LBD15 were significantly downregulated in Eucalyptus urophylla triploids. In addition, potential important transcription factors, including HSL1, BEE3, HHO 3, and NAC046, also had high degrees of connectivity and high edge weights with lignin biosynthetic genes, indicating that they might also be involved in the variation of lignin accumulation between triploid and diploid Eucalyptus urophylla. Conclusions The results demonstrated that some lignin-related metabolites, lignin biosynthetic genes and transcription factors in Eucalyptus urophylla triploids may be relatively sensitive in response to the polyploidization effect, significantly changing their expression levels, which ultimately correlated with the varied lignin content. The analysis of the underlying formation mechanism could provide beneficial information for the development and utilization of polyploid biomass resources, which will be also valuable for genetic improvement in other bioenergy plants.

Abstract Background Network analysis is a powerful tool for studying gene regulation and identifying biological processes associated with gene function. However, constructing gene co-expression networks can be a challenging task, particularly when dealing with a large number of missing values. Results We introduce GeCoNet-Tool, an integrated gene co-expression network construction and analysis tool. The tool comprises two main parts: network construction and network analysis. In the network construction part, GeCoNet-Tool offers users various options for processing gene co-expression data derived from diverse technologies. The output of the tool is an edge list with the option of weights associated with each link. In network analysis part, the user can produce a table that includes several network properties such as communities, cores, and centrality measures. With GeCoNet-Tool, users can explore and gain insights into the complex interactions between genes.

Marigold black spot caused by Alternaria tagetica is a major disease that can decrease marigold production by 40%, resulting in serious economic losses. In this study, we identified many genes responsive to A. tagetica in the resistant and susceptible marigold genotypes. Analyses of differentially expressed genes, expression trends, and a weighted gene co-expression network revealed a series of hub genes with key roles in different A. tagetica infection stages. Additionally, 1 216 unigenes encoding transcription factors from eight families were differentially expressed between Ts and Ma. Moreover, R genes from various families (e.g., N, NL, RLP, and TNL) were differentially expressed in the two marigold genotypes before and after the inoculation with A. tagetica. Pathway diagrams were used to visualize the leaf transcriptional changes in the two marigold genotypes infected by A. tagetica to clarify the effects of A. tagetica on the expression patterns of genes involved in phosphatidylinositol signaling, plant–pathogen interactions, and plant hormone signal transduction. We identified candidate genes related to disease resistance and generated valuable resources for analyzing the candidate gene functions related to black spot resistance in marigold. The study data may be useful for the molecular marker-assisted screening and breeding of marigold lines with increased disease resistance.

Despite advancements in treatment, high‐grade serous ovarian cancer (HGSOC) is still characterized by poor patient outcomes. To understand the molecular heterogeneity of this disease, which underlies the challenge in selecting optimal treatments for HGSOC patients, we have integrated genomic, transcriptomic, and epigenetic information to identify seven new HGSOC subtypes using a multiscale clustering method. These subtypes not only have significantly distinct overall survival, but also exhibit unique patterns of gene expression, microRNA expression, DNA methylation, and copy number alterations. As determined by our analysis, patients with similar clinical outcomes have distinct profiles of activated or repressed cellular processes, including cell cycle, epithelial‐to‐mesenchymal transition, immune activation, interferon response, and cilium organization. Furthermore, we performed a multiscale gene co‐expression network analysis to identify subtype‐specific key regulators and predicted optimal targeted therapies based on subtype‐specific gene expression. In summary, this study provides new insights into the cellular heterogeneity of the HGSOC genomic, epigenetic, and transcriptomic landscapes and provides a basis for future studies into precision medicine for HGSOC patients.

The BBAA subgenomes of hexaploid common wheat can be 'extracted' to constitute a viable and self-reproducing novel tetraploid wheat, termed extracted tetraploid wheat (ETW). Prior studies have shown ETW manifesting phenotypic abnormalities and alteration in gene expression and epigenetic modifications. No population level investigation has been conducted, leaving the issue unclear regarding whether developmental stability, an essential property evolved in all natural organisms, might have been undermined in ETW. Here, we measured variations in five morphological traits and somatic chromosomal stability in populations of ETW and of its hexaploid donor, a resynthesized hexaploid and a natural tetraploid wheat. We observed phenotypic defects in ETW. Meanwhile, we documented much greater within-population variations in ETW than in the other wheat genotypes, most probably due to disrupted developmental stability in ETW. Also, somatic structural chromosome variations were detected only in ETW. Comparative transcriptome analyses indicated that the disrupted developmental stability of ETW is likely linked to massive dysregulation of genome-wide gene expression rather than to genetic mutations. Population network analysis of gene expression implicated intrinsic connectivity among the variable traits, while gene set enrichment analysis provided possible links between dysregulated gene expression and interlaced trait variation. [ABSTRACT FROM AUTHOR]

Gastric cancer is one of the most commonly occurring malignancies. While researches have led to the discovery of an increasing number of pathogenic factors, the pathogenesis is not yet clearly known. The purpose of this study is to summarize the existing studies on gastric cancer, and further explore the key role of non-coding genes in the progression and prognosis of gastric cancer. Taking into consideration the expression data of coding and non-coding genes of gastric cancer patients and known gastric cancer pathogenic genes， co-expression calculation of coding and non-coding genes was performed to identify lncRNAs mediated by miRNA that interact with known gastric cancer genes. The three interacting modules (mRNA-miRNA-lncRNA) were mined， the modules were then screened, and the genes of significant modules related to disease were subject to survival analysis. In addition to known gastric cancer-related genes， differentially expressed gastric cancer genes were also employed in study， which are significantly enriched in biological processes closely related to the basic functions of the stomach and development of gastric cancer， such as cell proliferation， cell adhesion， muscle contraction， vascular remodeling， cell division， and chromosome separation. The gene BGN in the triplet module with the highest score was not only significantly overexpressed in gastric cancer patients， but also significantly correlated with the prognosis of gastric cancer patients. The relations of the has-miRNA-153-5p and the has-miRNA-5001-5p in this module to gastric cancer were both confirmed. We then infer that the abnormal expression of mRNA and miRNA in the module may be caused by abnormal expression of lncRNA that is significantly related to them. Conclusion (s): This study has found a new breakthrough for the study of abnormal expression of known pathogenic genes in gastric cancer. The discovery of potential non-coding genes related to gastric cancer provides new targets for the prevention and treatment of gastric cancer, and provides basis for future clinical applications. [ABSTRACT FROM AUTHOR]

Background. Neurogenic bladder (NB) following suprasacral spinal cord injury (SSCI) is an interstitial disease with the structural remodeling of bladder tissue and matrix over-deposition. Circular RNAs (circRNAs) are involved in fibrotic disease development through their post-transcriptional regulatory functions. This study aimed to use transcriptome high-throughput sequencing to investigate the process of NB and bladder fibrosis after SSCI. Methods. Spinal cord transection at the T10≥T11 level was used to construct the SSCI model in rats (10≥week≥old female Wistar rats, weighing 200±20 g). The bladders were collected without (sham group) and with (SSCI 1≥3 groups) NB status. Morphological examination was conducted to assess the extent of bladder fibrosis. Additionally, RNA sequencing was utilized to determine mRNAs and circRNAs expression patterns. The dynamic changes of differentially expressed mRNAs (DEMs) and circRNAs (DECs) in different periods of SSCI were further analyzed. Results. Bladder weight, smooth muscle cell hypertrophy, and extracellular matrix gradually increased after SSCI. Compared with the sham group, 3,255 DEMs and 1,339 DECs, 3,449 DEMs and 1,324 DECs, 884 DEMs, and 1,151 DECs were detected in the SSCI 1≥3 groups, respectively. Specifically, circRNA3621, circRNA0617, circRNA0586, and circRNA4426 were significant DECs common to SSCI 1≥3 groups compared with the sham group. Moreover, Gene Ontology (GO) enrichment suggested that inflammatory and chronic inflammatory responses were the key events in NB progression following SSCI. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment associated with the "Chemokine signaling pathway", the "IL-17 signaling pathway", and the "TGF-beta signaling pathway" suggests their potential involvement in regulating biological processes. The circRNA≥miRNA≥mRNA interaction networks of DECs revealed rno-circ-2239 (micu2) as the largest node, indicating that the rno-circ-2239≥miRNA≥mRNA≥mediated network may play a critical role in the pathogenesis of SSCI-induced NB. Conclusions. This study offers a comprehensive outlook on the possible roles of DEMs and DECs in bladder fibrosis and NB progression following SSCI. These findings have the potential to serve as novel biomarkers and therapeutic targets. [ABSTRACT FROM AUTHOR]

Colorectal cancer (CRC) is the third most diagnosed cancer in the world. A better understanding of the molecular mechanism of CRC is essential for making novel strategies for the CRC management and its prevention. The present study aims to explore the molecular mechanism through integrated bioinformatics analysis by analyzing genes and their co-expression pattern in normal and CRC states. GSE110223, GSE110224 and GSE113513 gene expression profiles were analyzed in this study. The co-expression networks for normal and tumor samples were constructed separately and analyzed to identify the modules, sub-networks and key genes. Gene regulatory network analysis was done to understand the regulatory mechanism of selected genes. Survival analysis was performed for the identified sub-networks and key genes to understand their role in CRC progression. A total of seven modules were detected and the KEGG pathway analysis revealed these modules were mainly enriched with cell cycle, metabolism and signaling-related pathways. E2F6 and ETV4 transcription factors regulating the activity of multiple genes of identified modules were found to be up-regulated in CRC. Six Sub-networks and seven key genes, BORA, CCT7, DTL, RUVBL1, RUVBL2, THEM6 and TMEM97 associated with the CRC progression were identified. Disease-gene association analysis identified a novel association of the BORA gene with CRC that activates and regulates the AURORA-PLK1 cascades in the cell cycle. Survival analysis indicates that the overexpressed BORA is associated with unfavourable overall survival in CRC. The mechanistic role of BORA in the regulation of cell cycle progression suggests that BORA might act as a potential therapeutic target for CRC. [ABSTRACT FROM AUTHOR]

Cultivated peanut (Arachis hypogaea L.) is an important economic and oilseed crop worldwide, providing high-quality edible oil and high protein content. Seed size/weight and oil content are two important determinants of yield and quality in peanut breeding. To identify key regulators controlling these two traits, two peanut cultivars with contrasting phenotypes were compared to each other, one having a larger seed size and higher oil content (Zhonghua16, ZH16 for short), while the second cultivar had smaller-sized seeds and lower oil content (Zhonghua6, ZH6). Whole transcriptome analyses were performed on these two cultivars at four stages of seed development. The results showed that ~40% of the expressed genes were stage-specific in each cultivar during seed development, especially at the early stage of development. In addition, we identified a total of 5356 differentially expressed genes (DEGs) between ZH16 and ZH6 across four development stages. Weighted gene co-expression network analysis (WGCNA) based on DEGs revealed multiple hub genes with potential roles in seed size/weight and/or oil content. These hub genes were mainly involved in transcription factors (TFs), phytohormones, the ubiquitin–proteasome pathway, and fatty acid synthesis. Overall, the candidate genes and co-expression networks detected in this study could be a valuable resource for genetic breeding to improve seed yield and quality traits in peanut. [ABSTRACT FROM AUTHOR]

Enhancing the absorption and utilization of phosphorus by crops is an important aim for ensuring food security worldwide. However, the gene regulatory network underlying phosphorus use in foxtail millet remains unclear. In this study, the molecular mechanism underlying low-phosphorus (LP) responsiveness in foxtail millet was evaluated using a comparative transcriptome analysis. LP reduced the chlorophyll content in shoots, increased the anthocyanin content in roots, and up-regulated purple acid phosphatase and phytase activities as well as antioxidant systems (CAT, POD, and SOD). Finally, 13 differentially expressed genes related to LP response were identified and verified using transcriptomic data and qRT-PCR. Two gene co-expression network modules related to phosphorus responsiveness were positively correlated with POD, CAT, and PAPs. Of these, SiPHR1, functionally annotated as PHOSPHATE STARVATION RESPONSE 1, was identified as an MYB transcription factor related to phosphate responsiveness. SiPHR1 overexpression in Arabidopsis significantly modified the root architecture. LP stress caused cellular, physiological, and phenotypic changes in seedlings. SiPHR1 functioned as a positive regulator by activating downstream genes related to LP tolerance. These results improve our understanding of the molecular mechanism underlying responsiveness to LP stress, thereby laying a theoretical foundation for the genetic modification and breeding of new LP-tolerant foxtail millet varieties. [ABSTRACT FROM AUTHOR]

Background: Song performed in flocks by European starlings (Sturnus vulgaris), referred to here as gregarious song, is a non-sexual, social behavior performed by adult birds. Gregarious song is thought to be an intrinsically reinforced behavior facilitated by a low-stress, positive affective state that increases social cohesion within a flock. The medial preoptic area (mPOA) is a region known to have a role in the production of gregarious song. However, the neurochemical systems that potentially act within this region to regulate song remain largely unexplored. In this study, we used RNA sequencing to characterize patterns of gene expression in the mPOA of male and female starlings singing gregarious song to identify possibly novel neurotransmitter, neuromodulator, and hormonal pathways that may be involved in the production of gregarious song. Results: Differential gene expression analysis and rank rank hypergeometric analysis indicated that dopaminergic, cholinergic, and GABAergic systems were associated with the production of gregarious song, with multiple receptor genes (e.g., DRD2, DRD5, CHRM4, GABRD) upregulated in the mPOA of starlings who sang at high rates. Additionally, co-expression network analyses identified co-expressing gene clusters of glutamate signaling-related genes associated with song. One of these clusters contained five glutamate receptor genes and two glutamate scaffolding genes and was significantly enriched for genetic pathways involved in neurodevelopmental disorders associated with social deficits in humans. Two of these genes, GRIN1 and SHANK2, were positively correlated with performance of gregarious song. Conclusions: This work provides new insights into the role of the mPOA in non-sexual, gregarious song in starlings and highlights candidate genes that may play a role in gregarious social interactions across vertebrates. The provided data will also allow other researchers to compare across species to identify conserved systems that regulate social behavior. [ABSTRACT FROM AUTHOR]

Objective: RDH5 (Retinol dehydrogenase 5) is one of the member of the shortchain dehydrogenase/reductase (SDR) family and plays a critical role in a variety of tumor processes. In this study, we analyzed the expression, diagnostic significance and gene function of RDH5 in breast cancer. Methods: The Gene Expression Profiling Interactive Analysis (GEPIA) database was used to determine the diagnostic and prognostic value of RDH5 and the University of Alabama at Birmingham (UALCAN) database analyzed its expression in different subtypes, clinical stages and altered signal pathways in breast cancer. Then we explored the co-expression gene of RDH5 in breast cancer and constructed its corresponding network through the cBioPortal database. The String database was used to determine its interactions with other proteins. Results: RDH5 expression was differentiated in various tumours compared to normal tissue, and was down-regulated in breast cancer tissue and different subtypes of breast cancer. There was no correlation between RDH5 and the survival of breast cancer. RDH5 expression was significantly associated with different clinical parameters. Additionally, the protein expression of RDH5 in breast cancer with altered pathways indicated it was closely associated with the SWI-SNF (switching/sucrose non-fermenting) complex status and the mammalian Target of Rapamycin (mTOR), WNT (Wingless/Integrated), MYC/MYCN (Myelocytomatosis oncogene/euroblastoma derived MYC), RTK (receptor Tyrosine kinases), and p53/Rb (tumor protein p53/Retinoblastoma) pathways. Interaction analysis with RDH5 and its co-expression genes showed that its main related functions were ethanol oxidation, regulation of heart contraction and negative regulation of tumor necrosis factor. Lastly, the Protein-protein interaction (PPI) of RDH5 showed that the most relevant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were retinol metabolism, vitamin digestion and absorption and metabolic pathways. Conclusion: These results implied that RDH5 might be a larvaceous diagnostic biomarker for breast cancer and was closely associated with its metabolism. However, more studies are needed to confirm our findings and support the clinical importance of RDH5 in breast cancer. [ABSTRACT FROM AUTHOR]

INTRODUCTION: Recent studies revealed the association of abnormal methylomic changes with Alzheimer's disease (AD) but there is a lack of systematic study of the impact of methylomic alterations over the molecular networks underlying AD. METHODS: We profiled genome‐wide methylomic variations in the parahippocampal gyrus from 201 post mortem control, mild cognitive impaired, and AD brains. RESULTS: We identified 270 distinct differentially methylated regions (DMRs) associated with AD. We quantified the impact of these DMRs on each gene and each protein as well as gene and protein co‐expression networks. DNA methylation had a profound impact on both AD‐associated gene/protein modules and their key regulators. We further integrated the matched multi‐omics data to show the impact of DNA methylation on chromatin accessibility, which further modulates gene and protein expression. DISCUSSION: The quantified impact of DNA methylation on gene and protein networks underlying AD identified potential upstream epigenetic regulators of AD. Highlights: A cohort of DNA methylation data in the parahippocampal gyrus was developed from 201 post mortem control, mild cognitive impaired, and Alzheimer's disease (AD) brains.Two hundred seventy distinct differentially methylated regions (DMRs) were found to be associated with AD compared to normal control.A metric was developed to quantify methylation impact on each gene and each protein.DNA methylation was found to have a profound impact on not only the AD‐associated gene modules but also key regulators of the gene and protein networks.Key findings were validated in an independent multi‐omics cohort in AD.The impact of DNA methylation on chromatin accessibility was also investigated by integrating the matched methylomic, epigenomic, transcriptomic, and proteomic data. [ABSTRACT FROM AUTHOR]

Aflatoxin (AF) contamination, caused by Aspergillus flavus, compromises the food safety and marketability of commodities, such as maize, cotton, peanuts, and tree nuts. Multigenic inheritance of AF resistance impedes conventional introgression of resistance traits into high-yielding commercial maize varieties. Several AF resistance-associated quantitative trait loci (QTLs) and markers have been reported from multiple biparental mapping and genome-wide association studies (GWAS) in maize. However, QTLs with large confidence intervals (CI) explaining inconsistent phenotypic variance limit their use in marker-assisted selection. Meta-analysis of published QTLs can identify significant meta-QTLs (MQTLs) with a narrower CI for reliable identification of genes and linked markers for AF resistance. Using 276 out of 356 reported QTLs controlling resistance to A. flavus infection and AF contamination in maize, we identified 58 MQTLs on all 10 chromosomes with a 66.5% reduction in the average CI. Similarly, a metaanalysis of maize genes differentially expressed in response to (a)biotic stresses from the to-date published literature identified 591 genes putatively responding to only A. flavus infection, of which 14 were significantly differentially expressed (-1.0 ≤ Log2Fc ≥ 1.0; p ≤ 0.05). Eight MQTLs were validated by their colocalization with 14 A. flavus resistance-associated SNPs identified from GWAS in maize. A total of 15 genes were physically close between the MQTL intervals and SNPs. Assessment of 12 MQTL-linked SSR markers identified three markers that could discriminate 14 and eight cultivars with resistance and susceptible responses, respectively. A comprehensive meta-analysis of QTLs and differentially expressed genes led to the identification of genes and makers for their potential application in marker-assisted breeding of A. flavus-resistant maize varieties. [ABSTRACT FROM AUTHOR]

Sepsis is a life-threatening condition characterized by a harmful host response to infection with organ dysfunction. Annually about 20 million people are dead owing to sepsis and its mortality rates is as high as 20%. However, no studies have been carried out to investigate sepsis from the system biology point of view, as previous research predominantly focused on individual genes without considering their interactions and associations. Here, we conducted a comprehensive exploration of genome-wide expression alterations in both mRNAs and long non-coding RNAs (lncRNAs) in sepsis, using six microarray datasets. Co-expression networks were conducted to identify mRNA and lncRNA modules, respectively. Comparing these sepsis modules with normal modules, we observed a homogeneous expression pattern within the mRNA/lncRNA members, with the majority of them displaying consistent expression direction. Moreover, we identified consistent modules across diverse datasets, consisting of 20 common mRNA members and two lncRNAs, namely CHRM3-AS2 and PRKCQ-AS1, which are potential regulators of sepsis. Our results reveal that the up-regulated common mRNAs are mainly involved in the processes of neutrophil mediated immunity, while the down-regulated mRNAs and lncRNAs are significantly overrepresented in T-cell mediated immunity functions. This study sheds light on the co-expression patterns of mRNAs and lncRNAs in sepsis, providing a novel perspective and insight into the sepsis transcriptome, which may facilitate the exploration of candidate therapeutic targets and molecular biomarkers for sepsis. [ABSTRACT FROM AUTHOR]