Search

Post-traumatic stress disorder (PTSD) is a severe mental illness with grave social, political, economic, and humanitarian implications. To apply the principles of personalized omics-based medicine to this psychiatric problem, we implemented our previously introduced drug efficiency index (DEI) to the PTSD gene expression datasets. Generally, omics-based personalized medicine evaluates individual drug action using two classes of data: (1) gene expression, mutation, and Big Data profiles, and (2) molecular pathway graphs that reflect the protein–protein interaction. In the particular case of the DEI metric, we evaluate the drug action according to the drug’s ability to restore healthy (control) activation levels of molecular pathways. We have curated five PTSD and one TRD (treatment-resistant depression) cohorts of next-generation sequencing (NGS) and microarray hybridization (MH) gene expression profiles, which, in total, comprise 791 samples, including 379 cases and 413 controls. To check the applicability of our DEI metrics, we have performed three differential studies with gene expression and pathway activation data: (1) case samples vs. control samples, (2) case samples after treatment or/and observation vs. before treatment, and (3) samples from patients positively responding to the treatment vs. those responding negatively or non-responding patients. We found that the DEI values that use the signaling pathway impact activation (SPIA) metric were better than those that used the Oncobox pathway activation level (Oncobox PAL) approach. However, SPIA, Oncobox PAL, and DEI evaluations were reliable only if there were differential genes between case and control, or treated and untreated, samples.

Avian influenza is a severe viral infection that has the potential to cause human pandemics. In particular, chickens are susceptible to many highly pathogenic strains of the virus, resulting in significant losses. In contrast, ducks have been reported to exhibit rapid and effective innate immune responses to most avian influenza virus (AIV) infections. To explore the distinct genetic programs that potentially distinguish the susceptibility/resistance of both species to AIV, the investigation of coincident SNPs (coSNPs) and their differing causal effects on gene functions in both species is important to gain novel insight into the varying immune-related responses of chickens and ducks. By conducting a pairwise genome alignment between these species, we identified coSNPs and their respective effect on AIV-related differentially expressed genes (DEGs) in this study. The examination of these genes (e.g., CD74, RUBCN, and SHTN1 for chickens and ABCA3, MAP2K6, and VIPR2 for ducks) reveals their high relevance to AIV. Further analysis of these genes provides promising effector molecules (such as IκBα, STAT1/STAT3, GSK-3β, or p53) and related key signaling pathways (such as NF-κB, JAK/STAT, or Wnt) to elucidate the complex mechanisms of immune responses to AIV infections in both chickens and ducks.

During the cultivation of Marsupenaeus japonicus, there are often obvious differences in the growth within the same family under the same food, water quality, and environment, which greatly affects cultivation efficiency. To explore the molecular mechanism of this growth difference, this study used RNA-seq technology to compare the transcriptomes of M. japonicus individuals with significant growth differences from the same family. A total of 1375 differentially expressed genes were identified, of which 1109 were upregulated and 266 were downregulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed on the differentially expressed genes, and the results showed that growth-related processes, including chitin metabolism, chitin-binding amino sugar metabolism, and antioxidant processes, including response to oxidative stress, peroxidase activity, antioxidant activity, and peroxidase activity, showed significant differences between the large-size and small-size groups. The expression levels of some differentially expressed genes, such as cuticular protein, low-density lipoprotein receptor, ecdysteroid kinase, myosin heavy chain, and apoptosis inhibitor, were verified by quantitative PCR experiments. One cuticle gene was annotated, and phylogenetic analysis showed that this sequence clustered with the penaeid cuticle genes. This study provides valuable data and a scientific basis for understanding the mechanism of growth differences in M. japonicus at the molecular-genetic level.

Golden mahseer (Tor putitora) is a critically endangered fish with significant economic importance. However, its reproductive challenges in a captive environment pose a limitation to the successful domestication and aquaculture potential of this species. To understand the role of various genes in gonad maturation and reproduction in golden mahseer, we conducted an RNA-sequencing (RNA-Seq) study on the brains of mature male and female specimens. Altogether, 20.6 and 21.5 million reads were generated from the brains of the male and female fish, respectively. A total of 26,989 and 55,600 cDNA coding sequences (CDS) were identified from the male and female brains, respectively, among which 26,258 CDS from the male brain and 53,446 CDS from the female brain demonstrated homology to known protein database sequences. A comprehensive analysis revealed a total of 1187 distinct differentially expressed upregulated genes (DEGs), encompassing 953 DEGs anticipated to exhibit upregulation in the female brain and 234 DEGs in the male brain. Furthermore, in the brain of female and male golden mahseer, a significant downregulation was observed in 492 and 744 genes, respectively, resulting in a cumulative count of 1236 downregulated genes. Validation of the RNA-Seq results was performed by quantitative real-time PCR (qPCR) using 24 genes. Sixteen candidate genes with differential expression levels between the male and female fish were then selected and analyzed via qPCR. The results confirmed that amh, foxl3, dax1, kif20, and tkt were upregulated in the male golden mahseer brain, while cyp19a1a, dmrt2a, gdf9, sox9b, wt-1a, and aqp1 were upregulated in the female golden mahseer brain. Our study elucidates the distinct gene expression profiles in male and female golden mahseer brains and thus offers valuable insights for potential reproductive manipulation strategies in this fish species.

To further investigate the immune response of Macrobrachium rosenbergii against Aeromonas veronii, comparative transcriptomic analyses of the M. rosenbergii hepatopancreas were conducted on challenge and control groups at 6, 12, and 24 h post-infection (hpi), independently. A total of 51,707 high-quality unigenes were collected from the RNA-seq data, and 8060 differentially expressed genes (DEGs) were discovered through paired comparisons. Among the three comparison groups, a KEGG pathway enrichment analysis showed that 173 immune-related DEGs were considerably clustered into 28 immune-related pathways, including the lysosome, the phagosome, etc. Moreover, the expression levels of the four key immune-related genes (TOLL, PAK1, GSK3β, and IKKα) were evaluated at various stages following post-infection in the hepatopancreas, hemolymph, and gills. Both PAK1 and GSK3β genes were highly up-regulated in all three tissues at 6 hpi with A. veronii; TOLL was up-regulated in the hepatopancreas and hemolymph but down-regulated in the gill at 6 hpi, and IKKα was up-regulated in hemolymph and gill, but down-regulated in the hepatopancreas at 6 hpi. These findings lay the groundwork for understanding the immune mechanism of M. rosenbergii after contracting A. veronii.

For a half-century, the commercial wild silkworm, Antheraea pernyi, has been protected by coumaphos, which is an internal organophosphorus insecticide used to kill the potential parasitic fly larvae inside. Knowledge about the detoxification genes of A. pernyi as well as the detoxification mechanism for this species remains severely limited. In this study, we identified 281 detoxification genes (32 GSTs, 48 ABCs, 104 CYPs, and 97 COEs) in the genome of this insect, which are unevenly distributed over 46 chromosomes. When compared to the domesticated silkworm, Bombyx mori, a lepidopteran model species, A. pernyi has a similar number of ABCs, but a greater number of GSTs, CYPs, and COEs. By transcriptome-based expression analysis, we found that coumaphos at a safe concentration level significantly changed the pathways related to ATPase complex function and the transporter complex in A. pernyi. KEGG functional enrichment analysis indicated that protein processing in the endoplasmic reticulum was the most affected pathway after coumaphos treatment. Finally, we identified four significantly up-regulated detoxification genes (ABCB1, ABCB3, ABCG11, and ae43) and one significantly down-regulated detoxification gene (CYP6AE9) in response to coumaphos treatment, suggesting that these five genes may contribute to detoxification of coumaphos in A. pernyi. Our study provides the first set of detoxification genes for wild silkworms from Saturniidae and highlights the importance of detoxification gene repertoire in insect pesticide tolerance.

Mammals experience similar stages of embryonic development, birth, infancy, youth, adolescence, maturity, and senescence. While embryonic developmental processes have been extensively researched, many molecular mechanisms regulating the different life stages after birth, such as aging, remain unresolved. We investigated the conserved and global molecular transitions in transcriptional remodeling with age in dogs of 15 breeds, which revealed that genes underlying hormone level regulation and developmental programs were differentially regulated during aging. Subsequently, we show that the candidate genes associated with tumorigenesis also exhibit age-dependent DNA methylation patterns, which might have contributed to the tumor state through inhibiting the plasticity of cell differentiation processes during aging, and ultimately suggesting the molecular events that link the processes of aging and cancer. These results highlight that the rate of age-related transcriptional remodeling is influenced not only by the lifespan, but also by the timing of critical physiological milestones.

Trichomonas vaginalis is the most common non-viral sexually transmitted infection. 5-nitroimidazoles are the only FDA-approved medications for T. vaginalis treatment. However, 5-nitroimidazole resistance has been increasingly recognized and may occur in up to 10% of infections. We aimed to delineate mechanisms of T. vaginalis resistance using transcriptome profiling of metronidazole (MTZ)-resistant and sensitive T. vaginalis clinical isolates. In vitro, 5-nitroimidazole susceptibility testing was performed to determine minimum lethal concentrations (MLCs) for T. vaginalis isolates obtained from women who had failed treatment (n = 4) or were successfully cured (n = 4). RNA sequencing, bioinformatics, and biostatistical analyses were performed to identify differentially expressed genes (DEGs) in the MTZ-resistant vs. sensitive T. vaginalis isolates. RNA sequencing identified 304 DEGs, 134 upregulated genes and 170 downregulated genes in the resistant isolates. Future studies with more T. vaginalis isolates with a broad range of MLCs are needed to determine which genes may represent the best alternative targets in drug-resistant strains.

Breast cancer is one of the most prevalent types of cancer diagnosed globally and continues to have a significant impact on the global number of cancer deaths. Despite all efforts of epidemiological and experimental research, therapeutic concepts in cancer are still unsatisfactory. Gene expression datasets are widely used to discover the new biomarkers and molecular therapeutic targets in diseases. In the present study, we analyzed four datasets using R packages with accession number GSE29044, GSE42568, GSE89116, and GSE109169 retrieved from NCBI-GEO and differential expressed genes (DEGs) were identified. Protein–protein interaction (PPI) network was constructed to screen the key genes. Subsequently, the GO function and KEGG pathways were analyzed to determine the biological function of key genes. Expression profile of key genes was validated in MCF-7 and MDA-MB-231 human breast cancer cell lines using qRT-PCR. Overall expression level and stage wise expression pattern of key genes was determined by GEPIA. The bc-GenExMiner was used to compare expression level of genes among groups of patients with respect to age factor. OncoLnc was used to analyze the effect of expression levels of LAMA2, TIMP4, and TMTC1 on the survival of breast cancer patients. We identified nine key genes, of which COL11A1, MMP11, and COL10A1 were found up-regulated and PCOLCE2, LAMA2, TMTC1, ADAMTS5, TIMP4, and RSPO3 were found down-regulated. Similar expression pattern of seven among nine genes (except ADAMTS5 and RSPO3) was observed in MCF-7 and MDA-MB-231 cells. Further, we found that LAMA2, TMTC1, and TIMP4 were significantly expressed among different age groups of patients. LAMA2 and TIMP4 were found significantly associated and TMTC1 was found less correlated with breast cancer occurrence. We found that the expression level of LAMA2, TIMP4, and TMTC1 was abnormal in all TCGA tumors and significantly associated with poor survival. Indian Council of Medical Research | Ref. BMI/11/(35)/2020 MICINN | Ref. RYC-2017-22891

This study integrated four microarray datasets by Robust Rank Aggregation (RRA) method to identify the differentially expressed genes (DEG) in bovine mammary epithelial (BME) cells in response to Escherichia coli and Staphylococcus aureus infection. Furthermore, the GO function and KEGG pathway enrichment analysis of the integrated DEG were performed. Finally, the protein-protein interaction (PPI) network was constructed. A total of 72 integrated DEG were identified from the four datasets. The most significantly enriched terms within the integrated DEG were mainly involved in the immune response. The PPI network of DEG was constructed with 53 nodes. Seventeen genes, which constitute a significant module, were identified as hub genes. Among them, CD40, CXCL6, and NFKBIZ were further screened as the key genes and have the potential to become biomarkers of E. coli and S. aureus mastitis, considering the specificity of biomarkers for diseases. The identified key genes and pathways in this study can assist in the search for biomarkers for mastitis diagnosis and disease resistance breeding.

Ribosomal protein uL15 (RPL27a) carries a specific modification, hydroxylation, at the His39 residue, which neighbors the CCA terminus of the E-site-bound tRNA at the mammalian ribosome. Under hypoxia, the level of hydroxylation of this protein decreases. We transiently transfected HEK293T cells with constructs expressing wild-type uL15 or mutated uL15 (His39Ala) incapable of hydroxylation, and demonstrated that ribosomes containing both proteins are competent in translation. By applying RNA-seq to the total cellular and polysome-associated mRNAs, we identified differentially expressed genes (DEGs) in cells containing exogenous uL15 or its mutant form. Analyzing mRNA features of up- and down-regulated DEGs, we found an increase in the level of more abundant mRNAs and shorter CDSs in cells with uL15 mutant for both translated and total cellular mRNAs. The level of longer and rarer mRNAs, on the contrary, decreased. Our data show how ribosome heterogeneity can change the composition of the translatome and transcriptome, depending on the properties of the translated mRNAs.

This study estimated the effect of substituting fishmeal completely with cottonseed protein concentrate (CPC) in the diet of sturgeon (Acipenser schrenckii) on growth, digestive physiology, and hepatic gene expression. A control diet containing fishmeal and an experimental diet based on CPC was designed. The study was conducted for 56 days in indoor recirculating aquaculture systems. The results showed that weight gain, feed efficiency, and whole-body essential amino acids (EAAs) all decreased significantly in the experimental group, while whole-body non-essential amino acids (NEAAs) and serum transaminase activity increased (p < 0.05). The activity of digestive enzymes in the mid-intestine was significantly reduced (p < 0.05), and liver histology revealed fatty infiltration of hepatocytes. The hepatic transcriptome revealed an upregulation of genes linked to metabolism, including steroid biosynthesis, pyruvate metabolism, fatty acid metabolism, and amino acid biosynthesis. These findings indicate that fully replacing fishmeal with CPC harms A. schrenckii growth and physiology. This study provides valuable data for the development of improved aquafeeds and the use of molecular methods to evaluate the diet performance of sturgeon.

Pecan (Carya illinoinensis) is an economically important nut tree that is greatly affected by drought, limiting its production and distribution. Although the complete genome of the pecan was published several years ago, the molecular basis of the pecan’s response to drought remains unclear. In this study, we analyzed the high-throughput transcriptome data for pecans under 3, 6, 9, 12, and 15 days of drought stress compared with the controls. A total of 12,893 differentially expressed genes (DEGs) were identified under drought stress, with 11,684 of them showing significant changes after 15 d of drought treatment. Among these, 4448 genes were up-regulated while 7226 were down-regulated. The trend analysis revealed that DEGs could be classified into 20 clusters. Surprisingly, the majority of genes (6148) showed a gradual down-regulation, and 3683 genes showed a gradual up-regulation in response to drought. Gene ontology enrichment analysis showed that the DEGs were mainly enriched in biological processes. The KEGG pathway enrichment results indicated that the DEGs were mainly enriched in several pathways, including metabolic pathways, the biosynthesis of secondary metabolites, and plant hormone signal transduction processes. Among the DEGs, 457 protein kinase and 734 transcription factor genes were shown to be drought-responsive and may play key roles in the response to drought, and the expression patterns of selected candidate genes were further validated using quantitative real-time PCR. Collectively, these findings highlighted the multiple processes in pecans under drought stress and provided valuable insights into the further investigation of the functions of drought stress-responsive genes and the molecular basis of the pecan drought stress response.

Vedolizumab is efficacious in the treatment of Crohn’s disease (CD) and ulcerative colitis (UC). However, a significant proportion of patients present with a non-response. To investigate whether differences in the clinical response to vedolizumab is reflected in changes in gene expression levels in whole blood, samples were collected at baseline before treatment, and at follow-up after 10–12 weeks. Whole genome transcriptional profiles were established by RNA sequencing. Before treatment, no differentially expressed genes were noted between responders (n = 9, UC 4, CD 5) and non-responders (n = 11, UC 3, CD 8). At follow-up, compared with baseline, responders displayed 201 differentially expressed genes, and 51 upregulated (e.g., translation initiation, mitochondrial translation, and peroxisomal membrane protein import) and 221 downregulated (e.g., Toll-like receptor activating cascades, and phagocytosis related) pathways. Twenty-two of the upregulated pathways in responders were instead downregulated in non-responders. The results correspond with a dampening of inflammatory activity in responders. Although considered a gut-specific drug, our study shows a considerable gene regulation in the blood of patients responding to vedolizumab. It also suggests that whole blood is not optimal for identifying predictive pre-treatment biomarkers based on individual genes. However, treatment outcomes may depend on several interacting genes, and our results indicate a possible potential of pathway analysis in predicting response to treatment, which merits further investigation.

Although nitrogen (N) deficiency greatly affects N absorption and metabolism in barley, the global transcriptomic changes in morphological and physiological adaptation to altered N availability remains largely unclear. We conducted a comparative transcriptome analysis of roots in A9-29 (low N tolerant line of barley) and Hua 30 (low N-sensitive variety of barley) under low N conditions to elucidate the responses and the underlying molecular mechanism. The results demonstrated that the root architecture was strongly influenced and that the root morphological indexes (total root length, total root area surface, and root volume) were remarkably promoted in A9-29 compared to Hua30 under low N stress. The transcriptome analysis of roots identified 1779 upregulated differentially expressed genes (DEGs) and 1487 downregulated DEGs specifically expressed in A9-29 under low N stress. Specific DEGs in A9-29 were largely enriched in energy metabolism, lipid metabolism, and the metabolism of other amino acids. In addition, transcription factor genes ERFs and IAA-related genes were specifically expressed in A9-29. To conclude, this study could provide a foundation for improving low N tolerance in barley.

Through excellent absorption and transformation, the macrophyte Myriophyllum (M.) aquaticum can considerably remove phosphorus from wastewater. The results of changes in growth rate, chlorophyll content, and roots number and length showed that M. aquaticum could cope better with high phosphorus stress compared with low phosphorus stress. Transcriptome and differentially expressed genes (DEGs) analyses revealed that, when exposed to phosphorus stresses at various concentrations, the roots were more active than the leaves, with more DEGs regulated. M. aquaticum also showed different gene expression and pathway regulatory patterns when exposed to low phosphorus and high phosphorus stresses. M. aquaticum’s capacity to cope with phosphorus stress was maybe due to its improved ability to regulate metabolic pathways such as photosynthesis, oxidative stress reduction, phosphorus metabolism, signal transduction, secondary metabolites biosynthesis, and energy metabolism. In general, M. aquaticum has a complex and interconnected regulatory network that deals efficiently with phosphorus stress to varying degrees. This is the first time that the mechanisms of M. aquaticum in sustaining phosphorus stress have been fully examined at the transcriptome level using high-throughput sequencing analysis, which may indicate the direction of follow-up research and have some guiding value for its future applications.

Yue Wu,1,2,&ast; Zijun Zhao,3,&ast; Jinru Zhang,1,2 Yaye Wang,1,2 Xueqin Song1,2 1Department of Neurology, The Second Hospital of Hebei Medical University, Shijiazhuang, Hebei, 050000, People’s Republic of China; 2Neurological Laboratory of Hebei Province, Shijiazhuang, Hebei, People’s Republic of China; 3Department of Neurosurgery, The Second Hospital of Hebei Medical University, Shijiazhuang, Heibei, 050000, People’s Republic of China&ast;These authors contributed equally to this workCorrespondence: Xueqin SongDepartment of Neurology, The Second Hospital of Hebei Medical University, 215 Heping West Road, Shijiazhuang, Hebei, 050000, People’s Republic of China, Tel/Fax +86-318-2187209, Email shenjnksxq@163.comBackground: Inclusion body myositis (IBM) is a unique idiopathic inflammatory myopathy with unclear pathogenesis and poor prognosis. Although previous publications have identified some molecular biomarkers, the value of these biomarkers is unknown.Objective: To identify hub genes and signaling pathways related to IBM for understanding the IBM-related mechanisms and providing guidance for therapy development.Methods: Two microarray datasets (GSE3112 and GSE128470) were downloaded from the Gene Expression Omnibus (GEO) database. GEO2R was used to detect differentially expressed genes (DEGs) between IBM and normal muscle tissues. The hub genes were determined using protein–protein interaction (PPI) network in Cytoscape. The specific signaling pathways and biological functions of IBM were identified using GO, KEGG, and GSEA enrichment analyses. Moreover, CIBERSORT was applied to estimate the expression level of 22 immune cell types in IBM and normal muscle tissue. The relationship between the immune cell types and hub genes was then explored.Results: A total of 219 DEGs and 10 hub genes were identified. Enrichment analyses revealed that the chemokine signaling pathway, cellular response to interferon-gamma, and P53 pathway have crucial roles in IBM. Immune infiltration analyses showed that IBM was associated with high level of CD8 T cells, Tregs, and macrophages. Finally, five potential drugs were predicted for IBM patients through CMap (connectivity map) database.Conclusion: In this study, the underlying molecular mechanisms and immunological landscape of IBM were investigated, and thus may provide new directions for future research on IBM pathogenesis.Keywords: bioinformatic analysis, inclusion body myositis, differentially expressed genes, biological pathways, immune infiltrating cell

AIM:To profile gene expression changes induced by dexamethasone in human trabecular meshwork using bioinformatics analysis, and to elucidate the possible mechanisms involved in the pathogenesis of glucocorticoid-induced glaucoma(GIG). METHODS:The gene expression dataset GSE37474 was obtained from the Gene Expression Omnibus(GEO). GEO2R was utilized to identify differentially expressed genes(DEGs)in trabecular meshwork between the dexamethasone group and the control group. Gene Ontology(GO)function annotation and the Kyoto Encyclopedia of Genes and Genomes(KEGG)pathway enrichment were constructed using the DAVID database. STRING database and Cytoscape software were used to construct protein-protein interaction(PPI). The hub genes were screened by CytoHubba plug-in. Finally, the mRNA expression of key hub genes was verified by RT-PCR. RESULTS:Compared to normal trabecular meshwork, dexamethasone-treated trabecular meshwork had a total of 252 DEGs, with 141 genes up-regulated and 111 genes down-regulated. GO function annotation showed that DEGs were mostly located in the extracellular matrix, where they engaged in the biological processes of positive regulation of inflammation and extracellular matrix remodeling. KEGG pathway enrichment showed that DEGs were largely involved in vascular smooth muscle contraction, arachidonic acid metabolism, ether lipid metabolism and choline metabolism. The PPI network yielded seven hub genes, three of which were up-regulated(EDN1, FOS, and LPL)and four of which were down-regulated(CCL2, IGF1, PTGS2, CCL5). In RT-PCR, the mRNA expression levels of the seven hub genes matched those in the gene expression profile. CONCLUSION:Dexamethasone can cause dramatic changes in the gene expression profile in trabecular meshwork. The enriched pathways of DEGs and certain hub genes play an important role in the remodeling of the extracellular matrix and the regulation of aqueous humor outflow, providing a full knowledge of the molecular mechanism of GIG.

The underlying molecular mechanisms of colorectal cancer (CRC) has attracted great attention from the scholarly community. The aim of our study is to identify pivotal genes related to the pathogenesis and prognosis of CRC. We integrated five microarray datasets from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were analyzed with the limma package. DAVID and OmicShare tools were used for Gene Ontology (GO) and KEGG enrichment analysis. The protein–protein interaction (PPI) network of DEGs was constructed. The prognostic analysis of hub genes was performed through Kaplan Meier-plotter. Finally, potential drugs were predicted in the CMap database. Through five microarray datasets, a total of 90 DEGs were detected including 54 up-regulated and 36 down-regulated genes. Biological process analysis showed DEGs were mainly enriched in positive regulation of neutrophil chemotaxis, chemokine-mediated signaling pathway, and bicarbonate transport. Signaling pathway analysis indicated that DEGs played a vital in proximal tubule bicarbonate reclamation, cell cycle and progesterone-mediated oocyte maturation. The GEPIA database confirmed that overexpression levels of hub genes were significantly associated with better survival of patients. Finally, the 20 most significant small molecules were obtained based on the CMap database. Our study has identified novel candidate biomarkers, pathways, and kinases associated with CRC prognosis.

Defeng Pan,1,&ast; Yufei Zhou,2,&ast; Shengjue Xiao,1,&ast; Yue Hu,3,&ast; Chunyan Huan,1 Qi Wu,1 Xiaotong Wang,1 Qinyuan Pan,1 Jie Liu,1 Hong Zhu1 1Department of Cardiology, The Affiliated Hospital of Xuzhou Medical University, Xuzhou, Jiangsu, 221004, People’s Republic of China; 2Department of Cardiology, Shanghai Institute of Cardiovascular Diseases, Zhongshan Hospital and Institutes of Biomedical Sciences, Fudan University, Shanghai, 200032, People’s Republic of China; 3Department of General Practice, The Affiliated Hospital of Xuzhou Medical University, Xuzhou, Jiangsu, 221004, People’s Republic of China&ast;These authors contributed equally to this workCorrespondence: Defeng PanDepartment of Cardiology, The Affiliated Hospital of Xuzhou Medical University, 99 Huaihai West Road, Xuzhou, 221004, Jiangsu, People’s Republic of ChinaTel +86 516-83262017Email xzdefengpan@xzhmu.edu.cnIntroduction: Atrial fibrillation (AF) is the most prevalent sustained cardiac arrhythmia, but the molecular mechanisms underlying AF are not known. We aimed to identify the pivotal genes and pathways involved in AF pathogenesis because they could become potential biomarkers and therapeutic targets of AF.Methods: The microarray datasets of GSE31821 and GSE41177 were downloaded from the Gene Expression Omnibus database. After combining the two datasets, differentially expressed genes (DEGs) were screened by the Limma package. MicroRNAs (miRNAs) confirmed experimentally to have an interaction with AF were screened through the miRTarBase database. Target genes of miRNAs were predicted using the miRNet database, and the intersection between DEGs and target genes of miRNAs, which were defined as common genes (CGs), were analyzed. Functional and pathway-enrichment analyses of DEGs and CGs were performed using the databases DAVID and KOBAS. Protein–protein interaction (PPI) network, miRNA- messenger(m) RNA network, and drug-gene network was visualized. Finally, reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) was used to validate the expression of hub genes in the miRNA–mRNA network.Results: Thirty-three CGs were acquired from the intersection of 65 DEGs from the integrated dataset and 9777 target genes of miRNAs. Fifteen “hub” genes were selected from the PPI network, and the miRNA-mRNA network, including 82 miRNAs and 9 target mRNAs, was constructed. Furthermore, with the validation by RT-qPCR, macrophage migration inhibitory factor (MIF), MYC proto-oncogene, bHLH transcription factor (MYC), inhibitor of differentiation 1 (ID1), and C-X-C Motif Chemokine Receptor 4 (CXCR4) were upregulated and superoxide Dismutase 2 (SOD2) was downregulated in patients with AF compared with healthy controls. We also found MIF, MYC, and ID1 were enriched in the transforming growth factor (TGF)-β and Hippo signaling pathway.Conclusion: We identified several pivotal genes and pathways involved in AF pathogenesis. MIF, MYC, and ID1 might participate in AF progression through the TGF-β and Hippo signaling pathways. Our study provided new insights into the mechanisms of action of AF.Keywords: atrial fibrillation, bioinformatics analysis, differentially expressed genes, pathways, protein–protein interaction network

Yeletai Nurzati, Zhu Zhu, Heng Xu, Yixin Zhang Department of Plastic and Reconstructive Surgery, Shanghai Ninth People’s Hospital, Shanghai JiaoTong University School of Medicine, Shanghai, 201100, People’s Republic of ChinaCorrespondence: Heng Xu; Yixin Zhang Email xuh1990@gmail.com; zhangyx1576@163.comIntroduction: Keloid is a pathological scar type, which invades normal surrounding tissue without self-limiting to cause pain, itching, cosmetic disfigurement, etc. Knowledge of the molecular mechanisms underlying keloid remains unclear. This dilemma leads to no biomarker available for diagnosis. Thus, to seek accurate diagnosis, biomarkers are necessary for keloid diagnosis to help control its incidence.Methods: Gene Expression Omnibus (GEO) database was used to select differentially expressed miRNAs (DE-miRNAs) in GSE113620. miRTarBase miRNA–target tools were used to predict the interactions between miRNAs and their target mRNAs. Target mRNAs that were differentially expressed in keloid were selected by analyzing differentially expressed genes (DEGs) in GSE44270 and GSE92566. PPI network analysis, gene enrichment analysis, cell-specific and tissue-specific expression analyses of DE-target mRNAs were conducted. RT-PCR analysis was conducted to validate our results.Results: Three novel miRNAs (miR-30b-5p, miR-212-3p, miR-149-5p) and five target mRNAs (SIX1, CCNA2, CCNB1, FOXM1, RUNX2) were identified as potential biomarkers for keloid patients. Additionally, the potential functions of those miRNAs-mRNAs pathways were analyzed.Discussion: These findings of keloid-related miRNAs, mRNAs, and miRNA–mRNAs regulatory networks may provide insights into the underlying pathogenesis of keloid and serve as potential biomarkers for keloid diagnosis.Keywords: diagnostic biomarkers, miRNA–mRNA network, keloids, differentially expressed genes

Yang Guo,1– 4 Bobin Ning,5 Qunhui Zhang,1– 4 Jing Ma,4 Linlin Zhao,1– 4 QiQin Lu,1– 4 Dejun Zhang1,4 1Research Center for High Altitude Medicine, Medical College of Qinghai University, Xining, 810001, People’s Republic of China; 2Key Laboratory of Application and Foundation for High Altitude Medicine Research in Qinghai Province, Medical College of Qinghai University, Xining, 810001, People’s Republic of China; 3Qinghai-Utah Joint Research Key Lab for High Altitude Medicine, Medical College of Qinghai University, Xining, 810001, People’s Republic of China; 4Department of Eco-Environmental Engineering, Qinghai University, Xining, 810016, People’s Republic of China; 5Department of Medicine, The General Hospital of the People’s Liberation Army, Beijing, 100038, People’s Republic of ChinaCorrespondence: Dejun ZhangDepartment of Eco-Environmental Engineering, Qinghai University, Qinghai, 810016, People’s Republic of ChinaEmail djzhangqhu@163.comPurpose: The objective of this study was to identify the potential regulatory mechanisms, diagnostic biomarkers, and therapeutic drugs for heart failure (HF).Methods: Differentially expressed genes (DEGs) between HF and non-failing donors were screened from the GSE57345, GSE5406, and GSE3586 datasets. Database for Annotation Visualization and Integrated Discovery and Metascape were used for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses respectively. The GSE57345 dataset was used for weighted gene co-expression network analysis (WGCNA). The intersecting hub genes from the DEGs and WGCNA were identified and verified with the GSE5406 and GSE3586 datasets. The diagnostic value of the hub genes was calculated through receiver operating characteristic analysis and net reclassification index (NRI). Gene set enrichment analysis (GSEA) was used to filter out the signaling pathways associated with the hub genes. SYBYL 2.1 was used for molecular docking of hub targets and potential HF drugs obtained from the connection map.Results: Functional annotation of the DEGs showed enrichment of negative regulation of angiogenesis, endoplasmic reticulum stress response, and heart development. PTN, LUM, ISLR, and ASPN were identified as the hub genes of HF. GSEA showed that the key genes were related to the transforming growth factor-β (TGF-β) and Wnt signaling pathways. Sirolimus, LY-294002, and wortmannin have been confirmed as potential drugs for HF.Conclusion: We identified new hub genes and candidate therapeutic drugs for HF, which are potential diagnostic, therapeutic and prognostic targets and warrant further investigation.Keywords: differentially expressed genes, weighted gene co-expression network analysis, diagnostic biomarkers, therapeutic drugs, heart failure

Background: Papillary thyroid carcinoma (PTC) is one of the most common endocrine malignancies and has a favorable prognosis. However, optimal treatments and prognostic markers have not been clearly identified. Methods: Gene expression data from primary PTC were downloaded from the Gene Expression Omnibus database and subjected to two analyses of differentially expressed genes (DEGs), followed by intersecting individual and integrated DEGs analyses as well as gene set enrichment analysis. Analysis of data from Sequence Read Archive and The Cancer Genome Atlas, immunohistochemistry and qRT-PCR of TFF3 were performed to validate the results. Finally, the relationship between gene expression and disease-free survival as well as immune cell infiltration were investigated. Results: Six critical DEGs and several tumor-enriched signaling pathways were identified. Immunohistochemistry and qRT-PCR validated the low expression of TFF3 in PTC. TFF3 and FCGBP are coexpressed in PTC, and patients with lower gene expression had worse disease-free survival but higher immune cell infiltration. Conclusion: TFF3 was significantly underexpressed and may function with FCGBP synergistically in PTC.

The conditions for sex reversal in vertebrate species have been extensively studied, and the results highlighted numerous key factors involved in sex differentiation. However, the transcriptomes in hypothalamic and pituitary tissues from intersex goats have rarely been studied. The aim of this study was to screen candidate genes and signalling pathways related to sex reversal in Huai goats by analyzing gene expression in hypothalamic and pituitary tissues via transcriptome sequencing and bioinformatics analyses. In total, 612 and 139 differentially expressed genes (DEGs) were identified between the intersex and non‐intersex groups in the hypothalamus and pituitary, respectively. The DEGs in the hypothalamus and pituitary were significantly enriched in 41 and 16 signalling pathways, respectively, including the calcium signalling pathway, neuroactive ligand‐receptor interaction signalling pathway, and oestrogen signalling pathway, which might be related to intersex sex development disorders. A candidate gene from the tachykinin family (TACR1) was significantly enriched in the calcium signalling pathway. Thirty‐one DEGs were shared between these two comparisons and were enriched in several acetyl‐CoA‐related processes and the oestrogen signalling pathway. The results of the real‐time PCR analysis show that the transcriptome sequencing results were reliable. The transcriptome data indicate that the regulation of various physiological systems is involved in intersex goat development. Therefore, these results provide helpful data enhancing our understanding of the molecular mechanisms underlying intersex syndrome in goats., Our study compared the gene expression profiles of hypothalamus and pituitary tissues from intersex and normal goats. And 31 genes were shared between these two sets of comparisons. These genes were enriched in the oestrogen signalling pathway and acetyl‐CoA biosynthetic process, suggesting that they might exhibit important functions in the regulation of hormone biosynthesis related to gonad development.

Transcriptomic changes in plants during drought stress give insights into the mechanisms with which plants stabilise their metabolic processes in order to cope with the drought condition. In potato, understanding such drought-induced transcriptomic changes is critical because prolonged field drought interferes with tuber formation and bulking period of potato development, which eventually affects yield. We hypothesised that phenotypic drought responses of potato genotypes may be linked to differences in transcriptomic changes. Using an RNA sequencing approach, we investigated such transcriptomic changes in leaves of three cultivars (Biogold, Hansa and Lady Rosetta) under drought. We found more differentially expressed genes (DEGs) in the tolerant cultivars, Lady Rosetta and Biogold, than in the sensitive cultivar (Hansa). The differential gene expression trend reflected the phenotypic drought responses of the cultivars. For instance, we found in both Biogold and Lady Rosetta but not in Hansa, an upregulation of genes involved in carbohydrate metabolism (e.g., Alpha-glucosidase), flavonoid biosynthesis (e.g., Flavanone 3 beta-hydroxylase), lipid biosynthesis/transfer (e.g., nonspecific Lipid Transfer Proteins), heat shock proteins and secondary metabolites like phenolics and lignins. Furthermore, a prolonged drought stress resulted in reduced DEGs in Biogold and Hansa, but not in Lady Rosetta that also maintained its tuber yield under such prolonged stress suggesting a more robust drought tolerance. Our findings suggest that a synergistic expression of genes involved in several different aspects of drought response is required in order to obtain a robust tolerance.

Yang Guo,1– 4 Bobin Ning,5 Qunhui Zhang,1– 4 Jing Ma,4 Linlin Zhao,1– 4 QiQin Lu,1– 4 Dejun Zhang1,4 1Research Center for High Altitude Medicine, Medical College of Qinghai University, Xining, 810001, People’s Republic of China; 2Key Laboratory of Application and Foundation for High Altitude Medicine Research in Qinghai Province, Medical College of Qinghai University, Xining, 810001, People’s Republic of China; 3Qinghai-Utah Joint Research Key Lab for High Altitude Medicine, Medical College of Qinghai University, Xining, 810001, People’s Republic of China; 4Department of Eco-Environmental Engineering, Qinghai University, Xining, 810016, People’s Republic of China; 5Department of Medicine, The General Hospital of the People’s Liberation Army, Beijing, 100038, People’s Republic of ChinaCorrespondence: Dejun ZhangDepartment of Eco-Environmental Engineering, Qinghai University, Qinghai, 810016, People’s Republic of ChinaEmail djzhangqhu@163.comPurpose: The objective of this study was to identify the potential regulatory mechanisms, diagnostic biomarkers, and therapeutic drugs for heart failure (HF).Methods: Differentially expressed genes (DEGs) between HF and non-failing donors were screened from the GSE57345, GSE5406, and GSE3586 datasets. Database for Annotation Visualization and Integrated Discovery and Metascape were used for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses respectively. The GSE57345 dataset was used for weighted gene co-expression network analysis (WGCNA). The intersecting hub genes from the DEGs and WGCNA were identified and verified with the GSE5406 and GSE3586 datasets. The diagnostic value of the hub genes was calculated through receiver operating characteristic analysis and net reclassification index (NRI). Gene set enrichment analysis (GSEA) was used to filter out the signaling pathways associated with the hub genes. SYBYL 2.1 was used for molecular docking of hub targets and potential HF drugs obtained from the connection map.Results: Functional annotation of the DEGs showed enrichment of negative regulation of angiogenesis, endoplasmic reticulum stress response, and heart development. PTN, LUM, ISLR, and ASPN were identified as the hub genes of HF. GSEA showed that the key genes were related to the transforming growth factor-β (TGF-β) and Wnt signaling pathways. Sirolimus, LY-294002, and wortmannin have been confirmed as potential drugs for HF.Conclusion: We identified new hub genes and candidate therapeutic drugs for HF, which are potential diagnostic, therapeutic and prognostic targets and warrant further investigation.Keywords: differentially expressed genes, weighted gene co-expression network analysis, diagnostic biomarkers, therapeutic drugs, heart failure

Vaccines can be powerful tools, but for some diseases, safe and effective vaccines have been elusive. New developments in nucleic acid sequencing, bioinformatics, and protein modeling are facilitating the discovery of previously unknown antigens through reverse vaccinology approaches. Sequencing the complementarity- determining region of antibodies and T cell receptors allows detailed assessment of the immune repertoire and identification of paratopes shared by many individuals, supporting the selection of antigens that may be broadly protective. Systems vaccinology approaches to asses the global host response to vaccination by evaluation of differentially expressed genes in blood, cellular or tissue transcriptomes can reveal previously unknown pathways and interactions related to protective immunity. While it is important to remember that discoveries made through reverse vaccinology and systems vaccinology must still be confirmed with traditional challenge models and clinical trials, these approaches can provide new perspectives that may help solve longstanding problems in veterinary vaccinology.

Lodging is a major problem limiting maize yield worldwide. However, the mechanisms of lodging resistance remain incompletely understood for maize. Here, we evaluated 443 maize accessions for lodging resistance in the field. Five lodging-resistant accessions and five lodging-sensitive accessions were selected for further research. The leaf number, plant height, stem diameter, and rind penetrometer resistance were similar between lodging-resistant and -sensitive inbred lines. The average thickness of sclerenchymatous hypodermis layer was thicker and the vascular area was larger in the lodging-resistant lines compared with lodging-sensitive lines. Although total lignin content in stem tissue did not significantly differ between lodging-resistant and -sensitive lines, phloroglucinol staining revealed that the lignin content of the cell wall in the stem cortex and in the stem vascular tissue near the cortex was higher in the lodging-resistant lines than in the lodging-sensitive lines. Analysis of strand-specific RNA-seq transcriptome showed that a total of 793 genes were up-regulated and 713 genes were down-regulated in lodging-resistant lines relative to lodging-sensitive lines. The up-regulated genes in lodging-resistant lines were enriched in cell wall biogenesis. These results indicated that modification of cell wall biosynthesis would contribute to lodging resistance of maize.

Background The diameter of the abdominal aortic aneurysm (AAA) is the most commonly used parameter for the prediction of occurrence of AAA rupture. However, the most vulnerable region of the aortic wall may be different from the most dilated region of AAA under pressure. The present study is the first to use weighted gene coexpression network analysis (WGCNA) to detect the coexpressed genes that result in regional weakening of the aortic wall. Methods The GSE165470 raw microarray dataset was used in the present study. Differentially expressed genes (DEGs) were filtered using the “limma” R package. DEGs were assessed by Gene Ontology biological process (GO-BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. WGCNA was used to construct the coexpression networks in the samples with regional weakening of the AAA wall and in the control group to detect the gene modules. The hub genes were defined in the significant functional modules, and a hub differentially expressed gene (hDEG) coexpression network was constructed with the highest confidence based on protein–protein interactions (PPIs). Molecular compound detection (MCODE) was used to identify crucial genes in the hDEG coexpression network. Crucial genes in the hDEG coexpression network were validated using the GSE7084 and GSE57691 microarray gene expression datasets. Result A total of 350 DEGs were identified, including 62 upregulated and 288 downregulated DEGs. The pathways were involved in immune responses, vascular smooth muscle contraction and cell–matrix adhesion of DEGs in the samples with regional weakening in AAA. Antiquewhite3 was the most significant module and was used to identify downregulated hDEGs based on the result of the most significant modules negatively related to the trait of weakened aneurysm walls. Seven crucial genes were identified and validated: ACTG2, CALD1, LMOD1, MYH11, MYL9, MYLK, and TPM2. These crucial genes were associated with the mechanisms of AAA progression. Conclusion We identified crucial genes that may play a significant role in weakening of the AAA wall and may be potential targets for medical therapies and diagnostic biomarkers. Further studies are required to more comprehensively elucidate the functions of crucial genes in the pathogenesis of regional weakening in AAA.

Breast cancer (BC), an extremely aggressive malignant tumor, causes a large number of deaths worldwide. In this study, we pooled profile datasets from three cohorts to illuminate the underlying key genes and pathways of BC. Expression profiles GSE42568, GSE45827, and GSE124646, including 244 BC tissues and 28 normal breast tissues, were integrated and analyzed. Differentially expressed genes (DEGs) were screened out based on these three datasets. Functional analysis including Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway were performed using The Database for Annotation, Visualization and Integrated Discovery (DAVID). Moreover, Cytoscape with Search Tool for the Retrieval of Interacting Genes (STRING) and Molecular Complex Detection (MCODE) plugin were utilized to visualize protein protein interaction (PPI) of these DEGs. The module with the highest connectivity of gene interactions was selected for further analysis. All of these hub genes had a significantly worse prognosis in BC by survival analysis. Additionally, four genes (CDK1, CDC20, AURKA, and MCM4) dramatically were enriched in oocyte meiosis and cell cycle pathways through re-analysis of DAVID. Moreover, the mRNA and protein levels of CDK1, CDC20, AURKA, and MCM4 were significantly increased in BC patients. In addition, knockdown of CDK1 and CDC20 by small interfering RNA remarkably suppressed cell migration and invasion in MCF-7 and MDA-MB-231 cells. In conclusion, our results suggested that CDK1, CDC20, AURKA, and MCM4 were reliable biomarkers of BC via bioinformatics analysis and experimental validation and may act as prospective targets for BC diagnosis and treatment.

Background At the end of 2019, the world witnessed the emergence and ravages of a viral infection induced by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Also known as the coronavirus disease 2019 (COVID-19), it has been identified as a public health emergency of international concern (PHEIC) by the World Health Organization (WHO) because of its severity. Methods The gene data of 51 samples were extracted from the GSE150316 and GSE147507 data set and then processed by means of the programming language R, through which the differentially expressed genes (DEGs) that meet the standards were screened. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed on the selected DEGs to understand the functions and approaches of DEGs. The online tool STRING was employed to construct a protein–protein interaction (PPI) network of DEGs and, in turn, to identify hub genes. Results A total of 52 intersection genes were obtained through DEG identification. Through the GO analysis, we realized that the biological processes (BPs) that have the deepest impact on the human body after SARS-CoV-2 infection are various immune responses. By using STRING to construct a PPI network, 10 hub genes were identified, including IFIH1, DDX58, ISG15, EGR1, OASL, SAMD9, SAMD9L, XAF1, IFITM1, and TNFSF10. Conclusion The results of this study will hopefully provide guidance for future studies on the pathophysiological mechanism of SARS-CoV-2 infection.

Cancer is a global health problem. There are diverse types of cancers, but there are several common pathways which lead to the development of cancer. Changes in gene expression might be the most common similarity found in almost all cancers. An understanding of the underlying changes in gene expression during cancer progression could lay a valuable foundation for the development of cancer therapeutics and even cancer vaccines. In this study, a well-known carcinogen, N-methyl-N’- nitro-N-nitrosoguanidine (MNNG), was employed to induce changes in gene expression in normal stomach cells. MNNG is known to cause cancer by inducing damage to DNA in MNNG-treated mammalian cells and animals fed with this carcinogen. An analysis was performed by comparing the differentially expressed genes (DEGs) caused by MNNG treatment with DEGs in stomach cancer cell lines. To this end, methods of analysis for functional categorization and protein-protein interaction networks, such as gene ontology (GO), the database for annotation, visualization, and integrated discovery (DAVID), Kyoto encyclopedia of genes and genomics (KEGG) and search tool for the retrieval of interacting genes/proteins (STRING), were used. As a result of these analyses, MNNG-regulated specific genes and interaction networks of their protein products that contributed to stomach cancer were identified.

In recent years, Fusarium head blight (FHB) has developed into a global disease that seriously affects the yield and quality of wheat. Effective measures to solve this problem include exploring disease-resistant genes and breeding disease-resistant varieties. In this study, we conducted a comparative transcriptome analysis to identify the important genes that are differentially expressed in FHB medium-resistant (Nankang 1) and FHB medium-susceptible (Shannong 102) wheat varieties for various periods after Fusarium graminearum infection using RNA-seq technology. In total, 96,628 differentially expressed genes (DEGs) were identified, 42,767 from Shannong 102 and 53,861 from Nankang 1 (FDR < 0.05 and |log2FC| > 1). Of these, 5754 and 6841 genes were found to be shared among the three time points in Shannong 102 and Nankang 1, respectively. After inoculation for 48 h, the number of upregulated genes in Nankang 1 was significantly lower than that of Shannong 102, but at 96 h, the number of DEGs in Nankang 1 was higher than that in Shannong 102. This indicated that Shannong 102 and Nankang 1 had different defensive responses to F. graminearum in the early stages of infection. By comparing the DEGs, there were 2282 genes shared at the three time points between the two strains. GO and KEGG analyses of these DEGs showed that the following pathways were associated with disease resistance genes: response to stimulus pathway in GO, glutathione metabolism, phenylpropanoid biosynthesis, plant hormone signal transduction, and plant–pathogen interaction in KEGG. Among them, 16 upregulated genes were identified in the plant–pathogen interaction pathway. There were five upregulated genes, TraesCS5A02G439700, TraesCS5B02G442900, TraesCS5B02G443300, TraesCS5B02G443400, and TraesCS5D02G446900, with significantly higher expression levels in Nankang 1 than in Shannong 102, and these genes may have an important role in regulating the resistance of Nankang 1 to F. graminearum infection. The PR proteins they encode are PR protein 1-9, PR protein 1-6, PR protein 1-7, PR protein 1-7, and PR protein 1-like. In addition, the number of DEGs in Nankang 1 was higher than that in Shannong 102 on almost all chromosomes, except chromosomes 1A and 3D, but especially on chromosomes 6B, 4B, 3B, and 5A. These results indicate that gene expression and the genetic background must be considered for FHB resistance in wheat breeding.

WRKY-, PHD-, and MYB-like proteins are three important types of transcription factors in mungbeans, and play an important role in development and stress resistance. The genes’ structures and characteristics were clearly reported and were shown to contain the conservative WRKYGQK heptapeptide sequence, Cys4-His-cys3 zinc binding motif, and HTH (helix) tryptophan cluster W structure, respectively. Knowledge on the response of these genes to salt stress is largely unknown. To address this issue, 83 VrWRKYs, 47 VrPHDs, and 149 VrMYBs were identified by using comparative genomics, transcriptomics, and molecular biology methods in mungbeans. An intraspecific synteny analysis revealed that the three gene families had strong co-linearity and an interspecies synteny analysis showed that mungbean and Arabidopsis were relatively close in genetic relationship. Moreover, 20, 10, and 20 genes showed significantly different expression levels after 15 days of salt treatment (p < 0.05; Log2 FC > 0.5), respectively. Additionally, in the qRT-PCR analysis, VrPHD14 had varying degrees of response to NaCl and PEG treatments after 12 h. VrWRKY49 was upregulated by ABA treatment, especially in the beginning (within 24 h). VrMYB96 was significantly upregulated in the early stages of ABA, NaCl, and PEG stress treatments (during the first 4 h). VrWRKY38 was significantly upregulated by ABA and NaCl treatments, but downregulated by PEG treatment. We also constructed a gene network centered on the seven DEGs under NaCl treatment; the results showed that VrWRKY38 was in the center of the PPI network and most of the homologous Arabidopsis genes of the interacted genes were reported to have response to biological stress. Candidate genes identified in this study provide abundant gene resources for the study of salt tolerance in mungbeans.

In order to decipher the functional genes and reveal the molecular mechanism of gonadal development in Macrobrachium rosenbergii, a comparative transcriptome analysis was performed on the testes and ovaries at different developmental stages. A total of 146,537 unigenes with an N50 of 2008 bp and an average length of 1144 bp were obtained from the sequencing raw data via quality control and denovo assembly. Identification of differentially expressed genes (DEGs) showed that there were 339 and 468 DEGs among the different developmental stages of testes and ovaries, respectively, and 7993 DEGs between the testes and ovaries. The KEGG enrichment analysis identified 13 candidate pathways related to gonadal development, including insulin synthesis, oocyte maturation, and steroid biosynthesis, which were involved in biological processes such as regulation of hormone metabolism, sex cell proliferation and development, and amino acid metabolism. The DEGs related to the above pathways such as insulin-like growth factor 1 receptor (IGF1R), heat shock protein 90 (Hsp 90), and cyclooxygenase (COX) genes were highly expressed during yolk protein synthesis, indicating that these genes might be involved in yolk accumulation and oogenesis. Meanwhile, calmodulin (CaM) and other genes were highly expressed during spermatogenesis, suggesting that these genes might play an important role in spermatogenesis. Ten differentially expressed genes in the KEGG signaling pathway, including CRQ, COX, APP, Cdc42, Hsd17b12, Art-1, Hsp70, Hsp90, PRMT1, and GP, were selected for real-time quantitative PCR (RT- qPCR) to validate the transcriptome data, and the results showed that RT- qPCR obtained consistent results with the RNA-Seq data. The present findings provide new insights into the molecular regulation mechanism of gonadal development in M. rosenbergii.

Acanthamoeba keratitis (AK) is a blinding corneal infection caused by the protozoan Acanthamoeba. The long-term course of AK suggests the host immunity could not kill Acanthamoeba rapidly. The immune status is still unclear in the late stage of AK. The comparative transcriptome analysis was made based on the bulk RNA sequencing of cornea tissues from AK patients and donors. Differentially expressed genes and enriched signaling pathways were calculated. CIBERSORT algorithm was used for immune infiltration analysis of cornea tissue between AK and normal controls. A total of 2668 differentially expressed genes, including 1477 upregulated genes and 1191 downregulated genes, were detected. Gene Ontology analysis revealed that the pathways were significantly enriched in leukocyte migration, regulation of T-cell activation, the external side of plasma membrane, collagen-containing extracellular matrix, immune receptor activity, and cytokine binding. The Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the pathways were significantly enriched in the cytokine–cytokine receptor interaction, hematopoietic cell lineage, and Staphylococcus aureus infection pathway. The immune infiltration profiles varied little between AK and normal controls. Compared with normal tissue, cornea tissue of AK contained a higher proportion of M0 macrophages and CD8 T cells, while resting memory CD4 T cells contributed to a relatively lower portion (p < 0.05). Finally, the expression levels of cell markers and SLAMF7/STAT6 pathway were confirmed by histopathology examinations, RT-qPCR, and Western blot.

Severe coronavirus disease 2019 (COVID-19) has led to a rapid increase in death rates all over the world. Sepsis is a life-threatening disease associated with a dysregulated host immune response. It has been shown that COVID-19 shares many similarities with sepsis in many aspects. However, the molecular mechanisms underlying sepsis and COVID-19 are not well understood. The aim of this study was to identify common transcriptional signatures, regulators, and pathways between COVID-19 and sepsis, which may provide a new direction for the treatment of COVID-19 and sepsis. First, COVID-19 blood gene expression profile (GSE179850) data and sepsis blood expression profile (GSE134347) data were obtained from GEO. Then, we intersected the differentially expressed genes (DEG) from these two datasets to obtain common DEGs. Finally, the common DEGs were used for functional enrichment analysis, transcription factor and miRNA prediction, pathway analysis, and candidate drug analysis. A total of 307 common DEGs were identified between the sepsis and COVID-19 datasets. Protein–protein interactions (PPIs) were constructed using the STRING database. Subsequently, hub genes were identified based on PPI networks. In addition, we performed GO functional analysis and KEGG pathway analysis of common DEGs, and found a common association between sepsis and COVID-19. Finally, we identified transcription factor–gene interaction, DEGs-miRNA co-regulatory networks, and protein–drug interaction, respectively. Through ROC analysis, we identified 10 central hub genes as potential biomarkers. In this study, we identified SARS-CoV-2 infection as a high risk factor for sepsis. Our study may provide a potential therapeutic direction for the treatment of COVID-19 patients suffering from sepsis.

The marbled rockfish, Sebastiscus marmoratus, belongs to the euryhaline fish and is an oviparous scleractinian fish. There are few studies on the adaptation mechanism, functional genes, and related pathways of S. marmoratus and salinity. The results showed that a total of 72.1 GB of clean reads were obtained and all clean reads annotated a total of 25,278 Unigenes, of which 2,160 were novel genes. Compared to 20‰, 479 and 520 differential genes were obtained for 35‰ and 10‰, respectively. Gene ontology (GO) enrichment analysis revealed significant enrichment in protein binding, ion binding, ATP binding, and catalytic activity. Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that differentially expressed genes significantly expressed under salinity stress were mainly involved in the pathways of the cytochrome P450 metabolism of xenobiotics, tryptophan metabolism, cellular senescence, and calcium signaling pathways. Among them, pik3r6b, cPLA2γ-like, and WSB1 were differentially expressed in all three groups, and they were associated with apoptosis, inflammation, DNA damage, immune regulation, and other physiological processes. Six differentially expressed genes were randomly selected for qRT-PCR validation, and the results showed that the transcriptomic data were of high confidence.

The ichthyotoxic raphidophyte Heterosigma akashiwo is associated with harmful algal blooms (HABs). Although the influence of H. akashiwo on fish has been reported, little is known of its effects on other marine organisms, such as shellfish. In this study, the physiological and transcriptomic responses of the Asiatic hard clam Meretrix meretrix to H. akashiwo were investigated. The results showed that clam survival was not influenced by H. akashiwo (at a concentration of 2 × 104 cells/mL), and M. meretrix eliminated the alga effectively after 48 h in co-culture. H. akashiwo did not alter the antioxidant capacity of the clams. However, 544 unigenes were found to be differentially expressed between the control and H. akashiwo-treated clams. Interestingly, the significant up-regulation of lysosome-related genes and transport proteins suggested remodeling of the clam metabolism possibly for digestion of H. akashiwo. In summary, our results indicate that the clam M. meretrix can effectively feed on H. akashiwo, and intensification of shellfish farming may be helpful to control and mitigate H. akashiwo blooms.

Glioblastoma (GBM) is the most frequent primary brain tumor in adults and has a poor prognosis due to its resistance to Temozolomide (TMZ). However, there is limited research regarding the tumor microenvironment and genes related to the prognosis of TMZ-treated GBM patients. This study aimed to identify putative transcriptomic biomarkers with predictive value in patients with GBM who were treated with TMZ. Publicly available datasets from The Cancer Genome Atlas and Gene Expression Omnibus were analyzed using CIBERSORTx and Weighted Gene Co-expression Network Analysis (WGCNA) to obtain types of highly expressed cell types and gene clusters. Differentially Expressed Genes analysis was performed and was intersected with the WGCNA results to obtain a candidate gene list. Cox proportional-hazard survival analysis was performed to acquire genes related to the prognosis of TMZ-treated GBM patients. Inflammatory microglial cells, dendritic cells, myeloid cells, and glioma stem cells were highly expressed in GBM tissue, and ACP7, EPPK1, PCDHA8, RHOD, DRC1, ZIC3, and PRLR were significantly associated with survival. While the listed genes have been previously reported to be related to glioblastoma or other types of cancer, ACP7 was identified as a novel gene related to the prognosis of GBM. These findings may have potential implications for developing a diagnostic tool to predict GBM resistance and optimize treatment decisions.

Spermatogenesis is a complicated course of several rigorous restrained steps that spermatogonial stem cells undergo to develop into highly specialized spermatozoa; however, specific genes and signal pathways, which regulate the amplification, differentiation and maturation of these cells, remain unclear. We performed bioinformatics analyses to investigate the dynamic changes of the gene expression patterns at three time points in the course of the first wave of rat spermatogenesis. Differently expressed genes (DEGs) were identified, and the features of DEGs were further analyzed with GO (Gene Ontology), KEGG (Kyoto Encyclopedia of Genes and Genomes) and Short Time-series Expression Miner (STEM). A total of 2954 differentially expressed genes were identified. By using STEM, the top 10 key genes were selected in the profile according to the enrichment results, and the distinguishable biological functions encoded by these DEGs were automatically divided into three parts. Genes from 6, 8 and 10 days were related to biosynthesis, immune response and cell junction, and genes from 14, 15 and 16 days were related to energy metabolic pathways. The results also suggest that genes from 29, 31 and 35 days may shift metabolic to sperm motility, sperm flagellum and cilium movement.

The transcriptomic analysis of microarray and RNA-Seq datasets followed our own bioinformatic pipeline to identify a transcriptional regulatory network of lung cancer. Twenty-six transcription factors are dysregulated and co-expressed in most of the lung cancer and pulmonary arterial hypertension datasets, which makes them the most frequently dysregulated transcription factors. Co-expression, gene regulatory, coregulatory, and transcriptional regulatory networks, along with fibration symmetries, were constructed to identify common connection patterns, alignments, main regulators, and target genes in order to analyze transcription factor complex formation, as well as its synchronized co-expression patterns in every type of lung cancer. The regulatory function of the most frequently dysregulated transcription factors over lung cancer deregulated genes was validated with ChEA3 enrichment analysis. A Kaplan–Meier plotter analysis linked the dysregulation of the top transcription factors with lung cancer patients’ survival. Our results indicate that lung cancer has unique and common deregulated genes and transcription factors with pulmonary arterial hypertension, co-expressed and regulated in a coordinated and cooperative manner by the transcriptional regulatory network that might be associated with critical biological processes and signaling pathways related to the acquisition of the hallmarks of cancer, making them potentially relevant tumor biomarkers for lung cancer early diagnosis and targets for the development of personalized therapies against lung cancer.

Transcription factors (TF) are a wide class of genes in plants, and these can regulate the expression of other genes in response to various environmental stresses (biotic and abiotic). In the current study, transcription factor activity in sugarcane was examined during cold stress. Initially, RNA transcript reads of two sugarcane cultivars (ROC22 and GT08-1108) under cold stress were downloaded from SRA NCBI database. The reads were aligned into a reference genome and the differential expression analyses were performed with the R/Bioconductor edgeR package. Based on our analyses in the ROC22 cultivar, 963 TF genes were significantly upregulated under cold stress among a total of 5649 upregulated genes, while 293 TF genes were downregulated among a total of 3,289 downregulated genes. In the GT08-1108 cultivar, 974 TF genes were identified among 5,649 upregulated genes and 283 TF genes were found among 3,289 downregulated genes. Most transcription factors were annotated with GO categories related to protein binding, transcription factor binding, DNA-sequence-specific binding, transcription factor complex, transcription factor activity in RNA polymerase II, the activity of nucleic acid binding transcription factor, transcription corepressor activity, sequence-specific regulatory region, the activity of transcription factor of RNA polymerase II, transcription factor cofactor activity, transcription factor activity from plastid promoter, transcription factor activity from RNA polymerase I promoter, polymerase II and RNA polymerase III. The findings of above results will help to identify differentially expressed transcription factors during cold stress. It also provides a comprehensive analysis of the regulation of the transcription activity of many genes. Therefore, this study provides the molecular basis for improving cold tolerance in sugarcane and other economically important grasses. Resumo Fatores de transcrição (FT) são uma ampla classe de genes em plantas e podem regular a expressão de outros genes em resposta a vários estresses ambientais (estresses bióticos e abióticos). No presente estudo, a atividade do fator de transcrição na cana-de-açúcar foi examinada durante o estresse pelo frio. Inicialmente, as leituras de transcrição de RNA de duas cultivares de cana-de-açúcar (ROC22 e GT08-1108) sob estresse frio foram baixadas do banco de dados SRA NCBI. As leituras foram alinhadas em um genoma de referência e as análises de expressão diferencial foram realizadas com o pacote R / Bioconductor edgeR. Com base em nossas análises no cultivar ROC22, 963 genes TF foram significativamente regulados positivamente sob estresse pelo frio entre um total de 5.649 genes regulados positivamente, enquanto 293 genes TF foram regulados negativamente entre um total de 3.289 genes regulados negativamente. No cultivar GT08-1108, 974 genes TF foram identificados entre 5.649 genes regulados positivamente e 283 genes TF foram encontrados entre 3.289 genes regulados negativamente. Os fatores de transcrição, em sua maioria, foram anotados com categorias GO relacionadas à ligação de proteína, ligação de fator de transcrição, ligação específica de sequência de DNA, complexo de fator de transcrição, atividade de fator de transcrição em RNA polimerase II, atividade de fator de transcrição de ligação de ácido nucleico, atividade de corepressor de transcrição, sequência específica da região reguladora, atividade do fator de transcrição da RNA polimerase II, atividade do cofator do fator de transcrição, atividade do fator de transcrição do promotor do plastídio, atividade do fator de transcrição do promotor da RNA polimerase I, polimerase II e RNA polimerase III. As descobertas dos resultados acima ajudarão a identificar fatores de transcrição expressos diferencialmente durante o estresse pelo frio. Ele também fornece uma análise abrangente da regulação da atividade de transcrição de muitos genes. Portanto, este estudo fornece base molecular para melhorar a tolerância ao frio em cana-de-açúcar e outras gramíneas economicamente importantes.

Increases in chicken production are mainly due to specialised breeds. However, local breeds are of increasing importance, known for ability to adapt to the environment and unique products. Conventional poultry products contain lower levels of n-3 fatty acids (FAs) compared to those obtained from local breeds, therefore the aim of this study was to evaluate the modulation of expression of genes involved in long-chain polyunsaturated FA (PUFA) biosynthesis pathways according to genetic background, diet conditions, and sex. Animals from two local breeds and a commercial line were fed different diets: control and experimental diet (10% linseed supplementation). For each breed and diet group, both sexes were reared. The RNA was extracted from 36 liver samples and sequenced by RNAseq method. Bioinformatic analysis was carried out to find differentially expressed genes from comparisons between experimental groups. Results showed low impact of diet on differentially expressed genes related to FA biosynthesis, but linseed diet increased percentage of n-3 FAs of liver. Sex and genetic background determined the differential expression of genes related to long-chain PUFA biosynthesis. Specifically, females of local breeds shared 23 up-regulated genes when compared to their respective commercial line groups. Some of the shared genes had a role in de novo triglyceride biosynthesis (MTTPL and GPAM), and in de novo FA biosynthesis (ACACA and SCD) was detected. In conclusion, local breeds are able to better adapt to a diet rich in PUFA, by triggering certain transcriptomic shifts in the liver that allow birds to process the high PUFA content provided by diet.

Patients not responding to therapy as expected is one of the most pressing healthcare concerns of today. It causes economical, medical and societal issues along with suffering for patients. This project aimed to address this problem and evaluate how to find the best suited drug treatments for individual patients to treat Covid-19. This project was carried out in collaboration with the company AB Mavatar, that have two networks, one experimental and one predicted, which produce drug treatment rankings differently. Different methods are used to connect drug targets to disease associated genes and thus evaluate what drugs are best suited for specific patients to treat Covid-19. The aim of this project is to examine how network, method and drug category affect the ranking of a drug treatment for four mapped Covid-19 patients. Which drug category a drug belongs to did not seem to significantly affect the drug ranking. Yet, certain drug subcategories were closely correlated. However, these subcategories were not those that are typically associated with Covid-19. The method used to connect drug targets to disease associated genes heavily impacts the ranking of the drug treatment. The methods should be further evaluated to see if some should be excluded or weighted less in drug ranking calculations. The two networks are similar in how they rank different drugs, especially in severely ill patients. Through this project and the evaluation of the impact of method choice, one can start to figure out what should be prioritized among disease related changes. Also, important parameters for personalized treatment can be evaluated. Patienter som inte svarar på terapi som förväntat är en av de största utmaningarna inom hälso- och sjukvård idag. Det orsakar ekonomiska, medicinska och samhälleliga problem samt lidande för patienter. Det här projektet adresserade detta problem och evaluerade hur man kan hitta det bäst lämpade läkemedlet för specifika patienter för att behandla Covid-19. Projektet gjordes tillsammans med företaget AB Mavatar, som har två interaktom, en experimentell och en datadriven, som rangordnar läkemedelsrekommendationer på olika sätt. Olika metoder används för att koppla samman läkemedelsmål med sjukdomsrelaterade gener och således evaluera vilka läkemedel som är bäst lämpade för specifika patienter för behandling av Covid-19. Syftet med projektet var att undersöka hur nätverk, metod och läkemedelskategori påverkar hur läkemedel rangordnas för fyra kartlagda Covid-19-patienter. Vilken läkemedelskategori ett läkemedel tillhör tycks inte märkbart påverka läkemedelsrangordning. Trots detta var vissa läkemedelsunderkategorier nära korrelerade. Dock var dessa underkategorier inte typiskt associerade med Covid-19. Metoden för att koppla samman läkemedelsmål med sjukdomsassocierade gener påverkade läkemedelsrangordningen väsentligt. Metoderna borde dock evalueras ytterligare för att eventuellt exkludera eller vikta vissa mindre i uträkningar av läkemedelsrang. De två nätverken är lika i hur de rangordnar olika läkemedel, särskilt för svårt sjuka patienter. Genom detta projekt och genom evaluering av metodvalets påverkan kan man börja begripa hur man borde priorita bland sjukdomsrelaterade förändringar. Dessutom kunde viktiga parametrar inom personlig behandling evalueras.

Meat quality has a close relationship with fat and connective tissue; therefore, screening and identifying functional genes related to lipid metabolism is essential for the production of high-grade beef. The transcriptomes of the Longissimus dorsi muscle in Wagyu and Chinese Red Steppe cattle, breeds with significant differences in meat quality and intramuscular fat deposition, were analyzed using RNA-seq to screen for candidate genes associated with beef quality traits. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the 388 differentially expressed genes (DEGs) were involved in biological processes such as short-chain fatty acid metabolism, regulation of fatty acid transport and the peroxisome proliferator-activated receptor (PPAR) signaling pathway. In addition, crystallin alpha B (CRYAB), ankyrin repeat domain 2 (ANKRD2), aldehyde dehydrogenase 9 family member A1 (ALDH9A1) and enoyl-CoA hydratase and 3-hydroxyacyl CoA dehydrogenase (EHHADH) were investigated for their effects on intracellular triglyceride and fatty acid content and their regulatory effects on genes in lipogenesis and fatty acid metabolism pathways. This study generated a dataset from transcriptome profiling of two cattle breeds, with differing capacities for fat-deposition in the muscle, and revealed molecular evidence that CRYAB, ANKRD2, ALDH9A1 and EHHADH are related to fat metabolism in bovine fetal fibroblasts (BFFs). The results provide potential functional genes for maker-assisted selection and molecular breeding to improve meat quality traits in beef cattle.

Background We aim to discover some prognostic factors, provide a basis for discovering molecular markers, and provide a basis for molecular features of early lung adenocarcinoma (LUAD) to predict patient prognosis. Methods Sequence data of LUAD were downloaded from The Cancer Genome Atlas (TCGA) database to screen out differentially expressed lncRNAs, miRNAs, and mRNAs (DERs). DERs were identified using R software’s limma package. The competitive endogenous RNA (ceRNA) network was constructed based on these RNAs. Univariate and multivariate Cox regression analysis on the RNAs in the ceRNA screened out independent prognostic-related RNAs to construct a prognostic risk score (PS) model. Combined with clinical data, we can calculate the survival rate of patients with early LUAD. Results There were 2,701 differentially expressed mRNAs (DEmRNAs), 47 differentially expressed lncRNAs (DElncRNAs), and 161 differentially expressed miRNAs (DEmiRNAs) identified in early LUAD. Based on these RNAs, 32 lncRNAs, 87 miRNAs, and 174 mRNAs participated in the ceRNA network. Twelve independently prognostic-related RNAs form an optimized combination. We developed a PS model based on these RNAs. Age, tumor recurrence and PS model status were independent survival prognostic clinical factors. Nomogram was established to predict the 3-year and 5-year survival rates. Conclusions We successfully constructed a ceRNA regulatory network based on the DERs in early LUAD. It can help us clarify the molecular mechanism of early LUAD. Simultaneously, the prognostic-related RNAs in early LUAD were also screened out. This network could provide new bases for diagnoses and prognoses of patients with LUAD.

Background The prognosis of high grade serous ovarian cancer (HGSOC) patients is closely related to the immune microenvironment and immune response. Based on this, the purpose of this study was to construct a model to predict chemosensitivity and prognosis, and provide novel biomarkers for immunotherapy and prognosis evaluation of HGSOC. Methods GSE40595 (38 samples), GSE18520 (63 samples), GSE26712 (195 samples), TCGA (321 samples) and GTEx (88 samples) were integrated to screen differential expressed genes (DEGs) of HGSOC. The prognosis related DEGs (DEPGs) were screened through overall survival analysis. The DEGs-encoded protein–protein interaction network was constructed and hub genes of DEPGs (DEPHGs) were generated by STRING. Immune characteristics of the samples were judged by ssGSEA, ESTIMATE and CYBERSORT. TIMER was used to analyze the relationship between DEPHGs and tumor-infiltrating immunocytes, as well as the immune checkpoint genes, finally immune-related DEPHGs (IDEPHGs) were determined, and whose expression in 12 pairs of HGSOC tissues and tumor-adjacent tissues were analyzed by histological verification. Furthermore, the chemosensitivity genes in IDEPHGs were screened according to GSE15622 (n = 65). Finally, two prediction models of paclitaxel sensitivity score (PTX score) and carboplatin sensitivity score (CBP score) were constructed by lasso algorithm. The area under curve was calculated to estimate the accuracy of candidate gene models in evaluating chemotherapy sensitivity. Results 491 DEGs were screened and 37 DEGs were identified as DEPGs, and 11 DEPHGs were further identified. Among them, CXCL13, IDO1, PI3, SPP1 and TRIM22 were screened as IDEPHGs and verified in the human tissues. Further analysis showed that IDO1, PI3 and TRIM22 could independently affect the chemotherapy sensitivity of HGSOC patients. The PTX score was significantly better than TRIM22, PI3, SPP1, IDO1 and CXCL13 in predicting paclitaxel sensitivity, so was CBP score in predicting carboplatin sensitivity. What’s more, both of the HGSOC patients with high PTX score or high CBP score had longer survival time. Conclusions Five IDEPHGs identified through comprehensive bioinformatics analysis were closely related with the prognosis, immune microenvironment and chemotherapy sensitivity of HGSOC. Two prediction models based on IDEPHGs might have potential application of chemotherapy sensitivity and prognosis for patients with HGSOC.

BackgroundApplication of the wheat phytohormones abscisic acid (ABA) or gibberellic acid (GA) affect Fusarium head blight (FHB) disease severity; however, the molecular underpinnings of the elicited phenotypes remain unclear. Herein, the transcriptomic responses of an FHB-susceptible wheat cultivar ‘Fielder’ were characterized upon treatment with ABA, an ABA receptor antagonist (AS6), or GA in the presence or absence of Fusarium graminearum (Fg) challenge.ResultsA total of 30,876 differentially expressed genes (DEGs) where identified in ‘Fielder’ (26,004) and Fg (4,872). Fg challenge alone resulted in the most substantial wheat DEGs contributing to 57.2% of the total transcriptomic variation. Using a combination of topology overlap and correlation analyses, 9,689 Fg-related wheat DEGs were defined. Further enrichment analysis of the top 1% networked wheat DEGs identified critical expression changes within defense responses, cell structural metabolism, molecular transport, and membrane/lipid metabolism. Fg-challenged conditions also included the expression of a putative Fg ABA-biosynthetic cytochrome P450 and repression of wheat FUS3 for dysregulating ABA and GA crosstalk. ABA treatment alone elicited 4536 (32%) wheat DEGs common to those of the Fg-challenge, and Fg+ABA further enhanced 888 (12.5%) of them. These ABA elicited DEGs are involved in defense through both classical and non-classical phytohormone signaling and regulating cell wall structures including polyphenolic metabolism. Conversely, Fg+GA opposed 2239 (33%) Fg-elicited wheat DEGs, including modulating primary and secondary metabolism, defense responses, and flowering genes. ABA and jointly ABA⍰Fg⍰[Fg+ABA] treatments repressed, while Fg+GA induced an over-representation of wheat DEGs mapping to chromosome 6BL. Finally, compared to Fg+ABA, co-application of Fg+AS6 did not antagonize ABA biosynthesis or signal but rather elicited antagonistic Fg (557) and wheat (11) DEGs responses directly tied to stress responses, phytohormone transport, and FHB.ConclusionsComparative transcriptomics highlight the effects of wheat phytohormones on individual pathway and global metabolism simultaneously. Application of ABA may reduce FHB severity through misregulating defense mechanisms and cell wall fortification pathways. GA application may alter primary and secondary metabolism, creating a metabolic shift to ultimately reduce FHB severity. By comparing these findings to those previously reported for four additional plant genotypes, an additive model of the wheat-Fg interaction is proposed.

Kaidi Guo,1,2,&ast; Li Li1,2,&ast; 1Department of Gynecology and Oncology, Guangxi Medical University Cancer Hospital, Nanning, Guangxi, People’s Republic of China; 2Key Laboratory of Early Prevention and Treatment for Regional High Frequency Tumor, Ministry of Education, Nanning, Guangxi, People’s Republic of China&ast;These authors contributed equally to this workCorrespondence: Li LiDepartment of Gynecology and Oncology, Guangxi Medical University Cancer Hospital, Nanning, Guangxi, 530021, People’s Republic of ChinaTel +86 138 7811 3406Email lili@gxmu.edu.cnPurpose: Ovarian cancer is one of the common malignant tumors of female reproductive organs, which seriously threatens the life and health of women. Resistance to chemotherapeutic drugs for ovarian cancer is the root cause of recurrence in most patients. The purpose of this study is to determine the differentially expressed genes of platinum resistance in ovarian cancer, and to screen out molecular targets and diagnostic markers that could be used to treat ovarian cancer platinum resistance.Methods: We downloaded 5 gene microarray datasets GSE58470, GSE45553, GSE41499, GSE33482, and GSE15372 from the Gene Expression Omnibus database, all of which are associated with ovarian cancer platinum resistance. Subsequently, the intersection of the statistically significant differentially expressed genes in 5 gene chips was taken, and relevant bioinformatics and clinical parameters were performed on the screened differential genes. qRT-PCR was utilized to examine the mRNA expression levels in ovarian cancer sensitive and cisplatin-resistant cells.Results: Three differential genes, IFI27, JAG1, DNM3, may be closely related to platinum resistance of ovarian cancer, were screened by microarray datasets. According to the combined verification of bioinformatics, clinical case analyses and experiments, it was inferred that the increased expression of DNM3 was beneficial to patients with platinum resistance, but the high expression of IFI27 and JAG1 may lead to the risk of platinum resistance.Conclusion: IFI27, JAG1 and DNM3 screened by relevant gene chips may serve as new biomarkers of platinum resistance in ovarian cancer.Keywords: ovarian cancer, platinum resistance, bioinformatical analysis, differentially expressed genes