Related Gene Ontology Terms Research Articles

Transcriptomic profiling of Rana [Lithobates] catesbeiana back skin during natural and thyroid hormone-induced metamorphosis under different temperature regimes with particular emphasis on innate immune system components

As amphibians undergo thyroid hormone (TH)-dependent metamorphosis from an aquatic tadpole to the terrestrial frog, their innate immune system must adapt to the new environment. Skin is a primary line of defense, yet this organ undergoes extensive remodelling during metamorphosis and how it responds to TH is poorly understood. Temperature modulation, which regulates metamorphic timing, is a unique way to uncover early TH-induced transcriptomic events. Metamorphosis of premetamorphic tadpoles is induced by exogenous TH administration at 24 °C but is paused at 5 °C. However, at 5 °C a “molecular memory” of TH exposure is retained that results in an accelerated metamorphosis upon shifting to 24 °C. We used RNA-sequencing to identify changes in Rana (Lithobates) catesbeiana back skin gene expression during natural and TH-induced metamorphosis. During natural metamorphosis, significant differential expression (DE) was observed in >6500 transcripts including classic TH-responsive transcripts (thrb and thibz), heat shock proteins, and innate immune system components: keratins, mucins, and antimicrobial peptides (AMPs). Premetamorphic tadpoles maintained at 5 °C showed 83 DE transcripts within 48 h after TH administration, including thibz which has previously been identified as a molecular memory component in other tissues. Over 3600 DE transcripts were detected in TH-treated tadpoles at 24 °C or when tadpoles held at 5 °C were shifted to 24 °C. Gene ontology (GO) terms related to transcription, RNA metabolic processes, and translation were enriched in both datasets and immune related GO terms were observed in the temperature-modulated experiment. Our findings have implications on survival as climate change affects amphibia worldwide.

Alterations in gene expressions of Caco-2 cell responses to LPS and ploy(I:C) stimulation.

The intestinal epithelium barrier serves as a highly dynamic immunologic frontier in the defense against invading pathogenic bacteria and viruses. Hence, understanding of the complicated underlying relationship between enteric pathogens and the intestinal epithelium barrier is vital for developing strategies to improve the intestinal health of farm animals. To this end, Caco-2 cells were stimulated by 1 µg/ml lipopolysaccharide (LPS) for 24 h and 5 µg/ml polyinosinic-polycytidylic acid (ploy(I:C)) for 4 h to imitate bacterial and viral infection processes, respectively. The specific alterations in gene expression of Caco-2 cells after stimulation were characterized by transcriptome sequencing. Seventy differentially expressed genes (DEGs) were identified under LPS exposure, and 17 DEGs were observed under ploy(I:C) exposure. We found that most DEGs were specific, and only one common DEG SPAG7 was observed. Gene Ontology (GO) annotation analysis indicated that all DEGs identified in the different treatments were mainly derived from GO terms related to the maintenance of cellular homeostasis. Moreover, specific DEGs such as SLC39A10, MT2A, and MT1E regulated by LPS treatment, while IFIT2 and RUNX2 mediated by ploy(I:C) treatment, which are derived from immune function modulation related GO terms, were confirmed by both transcriptome sequencing and qRT-PCR. In addition, both transcriptome sequencing and qRT-PCR results verified that LPS specifically down-regulated the DEGs INHBE and ARF6, which are involved in inflammation responses related to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway including the TGF-beta signaling pathways and the Ras signaling pathway. Ploy(I:C) uniquely suppressed the DEGs GABARAP and LAMTOR3, which participated in viral replication-associated pathways including autophagy and mTOR signaling pathway.

Mating-Induced Common and Sex-Specific Behavioral, Transcriptional Changes in the Moth Fall Armyworm (Spodoptera frugiperda, Noctuidae, Lepidoptera) in Laboratory.

The intermediate process between mating and postmating behavioral changes in insects is still poorly known. Here, we studied mating-induced common and sex-specific behavioral and transcriptional changes in both sexes of Spodoptera frugiperda and tested whether the transcriptional changes are linked to postmating behavioral changes in each sex. A behavioral study showed that mating caused a temporary suppression of female calling and male courting behavior, and females did not lay eggs until the next day after the first mating. The significant differences on daily fecundity under the presence of males or not, and the same or novel males, suggest that females may intentionally retain eggs to be fertilized by novel males or to be fertilized competitively by different males. RNA sequencing in females revealed that there are more reproduction related GO (gene ontology) terms and KEGG (Kyoto encyclopedia of genes and genomes) pathways (mainly related to egg and zygote development) enriched to upregulated DEGs (differentially expressed genes) than to downregulated DEGs at 0 and 24 h postmating. In males, however, mating induced DEGs did not enrich any reproduction related terms/pathways, which may be because male reproductive bioinformatics is relatively limited in moths. Mating also induced upregulation on soma maintenance (such as immune activity and stress reaction) related processes in females at 0, 6 and 24 h postmating. In males, mating also induced upregulation on soma maintenance related processes at 0 h postmating, but induced downregulation on these processes at 6 and 24 h postmating. In conclusion, this study demonstrated that mating induced sex-specific postmating behavioral and transcriptional changes in both sexes of S. frugiperda and suggested that the transcriptional changes are correlated with postmating physiological and behavioral changes in each sex.

KRTA6A and FA2H Are Hub Genes Associated With Cgas-STING-related Immunogenic Cell Death in Lung Adenocarcinoma.

The immunogenic cell death (ICD) pathway plays a crucial prognostic role in lung adenocarcinoma (LUAD) therapy. The cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway is an upstream mechanism that drives ICD activation, but the interaction of hub genes remains unclear. The present study aimed to investigate the hub genes involved in ICD and the cGAS-STING pathway. The prognostic performance for hub genes and related Gene Ontology (GO) terms were also investigated. Gene expression data of ICD induction and cGAS-STING pathway activation samples were extracted from the Gene Expression Omnibus (GEO) database, and gene expression as well as clinical data of LUAD patients who received pharmaceutical therapy were extracted from The Cancer Genome Atlas (TCGA). Differentially expressed genes (DEGs) were identified, and protein-protein interaction (PPI) analysis was used to identify hub genes. Hazard risk (HR) scores were identified using Kaplan-Meier (K-M) and COX analyses. Gene set enrichment analysis (GSEA) was performed to identify the related GO terms, and receiver operating characteristic (ROC) analysis was used to evaluate the prognosis performance of the related gene sets. A total of 22 DEGs were identified in the two GEO datasets, and six hub genes were identified by PPI analysis. Keratin 6A (KRTA6A) and fatty acid 2-hydroxylase (FA2H) were selected as the hub genes after survival analysis. GSEA and ROC analysis indicated that there was no difference between the KRTA6A and FA2H gene sets on prognosis performance. KRTA6A and FA2H are hub genes associated with the induction of cGAS-STING-related ICD in LUAD.

Ontological Analysis of Coronavirus Associated Human Genes at the COVID-19 Disease Portal.

The COVID-19 pandemic stemmed a parallel upsurge in the scientific literature about SARS-CoV-2 infection and its health burden. The Rat Genome Database (RGD) created a COVID-19 Disease Portal to leverage information from the scientific literature. In the COVID-19 Portal, gene-disease associations are established by manual curation of PubMed literature. The portal contains data for nine ontologies related to COVID-19, an embedded enrichment analysis tool, as well as links to a toolkit. Using these information and tools, we performed analyses on the curated COVID-19 disease genes. As expected, Disease Ontology enrichment analysis showed that the COVID-19 gene set is highly enriched with coronavirus infectious disease and related diseases. However, other less related diseases were also highly enriched, such as liver and rheumatic diseases. Using the comparison heatmap tool, we found nearly 60 percent of the COVID-19 genes were associated with nervous system disease and 40 percent were associated with gastrointestinal disease. Our analysis confirms the role of the immune system in COVID-19 pathogenesis as shown by substantial enrichment of immune system related Gene Ontology terms. The information in RGD's COVID-19 disease portal can generate new hypotheses to potentiate novel therapies and prevention of acute and long-term complications of COVID-19.

Integrative Analysis of Proteome-wide Association Studies and Functional Enrichment Analysis to Identify Genes and Chemicals Associated with Alcohol Dependence.

Alcohol dependence accounts for a large proportion of the global burden of disease and disability. This study aims to investigate the candidate genes and chemicals associated with alcohol dependence. Using data from published alcohol dependence genome-wide association studies, we first conducted a proteome-wide association study of alcohol dependence by integrating alcohol dependence genome-wide association studies with 2 human brain reference proteomes of dorsolateral prefrontal cortex from the Religious Order Study and Rush Memory and Aging Project and the Banner Sun Health Research Institute. Then, based on the identified genes in proteome-wide association study, we conducted functional enrichment analysis and chemical-related functional enrichment analysis to detect the related Gene Ontology terms and chemicals. Proteome-wide association study identified several potential candidate genes for alcohol dependence, such as GOT2 ( P = 7.59 × 10 -6 ) and C3orf33 ( P = 5.00 × 10 -3 ). Furthermore, functional enrichment analysis identified multiple candidate Gene Ontology terms associated with alcohol dependence, such as glyoxylate metabolic process (adjusted P = 2.99 × 10 -6 ) and oxoglutarate metabolic process (adjusted P = 9.95 × 10 -6 ). Chemical-related functional enrichment analysis detected several alcohol dependence-related candidate chemicals, such as pitavastatin ( P = 2.00 × 10 -4 ), cannabinoids ( P = 4.00 × 10 -4 ), 11-nor-Δ(9)-tetrahydrocannabinol-9-carboxylic acid ( P = 4.00 × 10 -4 ), and gabapentin ( P = 2.00 × 10 -3 ). Our study reports multiple candidate genes and chemicals associated with alcohol dependence, providing novel clues for understanding the biological mechanism of alcohol dependence.

Toxicity Effects and Mechanisms of MgO Nanoparticles on the Oomycete Pathogen Phytophthora infestans and Its Host Solanum tuberosum.

Engineered nanoparticles have recently been used for innovation in agricultural disease management. However, both the toxicity effects and mechanisms of nanoparticles in target pathogens and their host plants are still largely unknown. Here, we found that magnesium oxide nanoparticles (MgO NPs) could protect potatoes against Phytophthora infestans (P. infestans) at a low dosage (50 μg/mL). Through scanning electron microscopy observation, antioxidant enzymes activity measurement, and gene transcriptome analysis, we found that the cell surfaces of P. infestans were destroyed, endogenous superoxide dismutase continuously remained in a higher active state, oxidoreductase activity-related gene ontology (GO) terms were enriched with upregulation, and transporter-activity related GO terms and six essential metabolism-related pathways were enriched with downregulation in P. infestans after 30 min MgO NPs treatment, whereas only 89 genes were changed without enriched GO and pathways terms, and no change in antioxidant activities and phenylalnine ammonialyase in potato appeared at 6 h post-MgO NPs treatment. Only the “plant hormone signal transduction pathway” was enriched with upregulation under differential expression analysis in potatoes. In conclusion, cell surface distortion, continuous oxidative stress, and inhibitions of membrane transport activity and metabolic pathways were toxic mechanisms of Mg ONPs in P. infestans, and the “plant hormone signal transduction pathway” was potentially regulated by Mg-ONPs without obviously harmful effects on potato after Mg ONPs exposure.

The functional insight into the genetics of cardiovascular disease: results from the post-GWAS study.

Cardiovascular diseases (CVDs), the leading cause of death worldwide, generally refer to a range of pathological conditions with the involvement of the heart and the blood vessels. A sizable fraction of the susceptibility loci is known, but the underlying mechanisms have been established only for a small proportion. Therefore, there is an increasing need to explore the functional relevance of trait-associated variants and, moreover, to search for novel risk genetic variation. We have reported the bioinformatic approach allowing effective identif ication of functional non-coding variants by integrated analysis of genome-wide data. Here, the analysis of 1361 previously identif ied regulatory SNPs (rSNPs) was performed to provide new insights into cardiovascular risk. We found 773,471 coding co-segregating markers for input rSNPs using the 1000 Genomes Project. The intersection of GWAS-derived SNPs with a relevance to cardiovascular traits with these markers was analyzed within a window of 10 Kbp. The effects on the transcription factor (TF) binding sites were explored by DeFine models. Functional pathway enrichment and protein–protein interaction (PPI) network analyses were performed on the targets and the extended genes by STRING and DAVID. Eighteen rSNPs were functionally linked to cardiovascular risk. A signif icant impact on binding sites of thirteen TFs including those involved in blood cells formation, hematopoiesis, macrophage function, inf lammation, and vasoconstriction was found in K562 cells. 21 rSNP gene targets and 5 partners predicted by PPI were enriched for spliceosome and endocytosis KEGG pathways, endosome sorting complex and mRNA splicing REACTOME pathways. Related Gene Ontology terms included mRNA splicing and processing, endosome transport and protein catabolic processes. Together, the f indings provide further insight into the biological basis of CVDs and highlight the importance of the precise regulation of splicing and alternative splicing

Analysis of gene expression profiles of Enterococcus faecalis induced by type I collagen

Enterococcus faecalis is an etiological agent of endodontic infections. The present study was performed to investigate the gene profiles of E. faecalis induced by type I collagen stimulation. E. faecalis ATCC 19433 was cultivated with [collagen (+)] or without type I collagen [collagen (−)], and transcriptome analysis was performed using high-throughput sequencing technology. A total of 3.6 gb of information was obtained by sequence analysis and 77 differentially expressed genes (DEGs) between the two culture conditions were identified. Among the 77 DEGs, 35 genes were upregulated in collagen (+) E. faecalis, whereas 42 genes were downregulated. Gene Ontology (GO) enrichment analysis was performed and 11 GO terms, including metalloendopeptidase activity (GO:0004222) and two related GO terms (GO:0031012, GO:0044421), were significantly enriched in the set of upregulated genes. We focused on an upregulated DEG belonging to the matrixin metalloprotease gene family, and matrix metalloprotease (MMP) activities of the bacterial cell were examined. The generic MMP, MMP-8, and MMP-9 activities of collagen (+) E. faecalis were significantly higher than those of collagen (−) E. faecalis. These results suggested that contact with type I collagen may alter the gene expression profile of E. faecalis, and upregulation of metalloprotease genes may result in enhanced MMP activities in E. faecalis.

Identification of hub genes associated with human osteoarthritis cartilage: An in silico approach

Osteoarthritis is a common orthopedic disease among the greatest causes of morbidity and disability worldwide; however, research on its pathogenesis and diagnostic methods remains limited. The present study focuses on analyzing the microarray dataset to elucidate the hub genes and pathways related to the osteoarthritis. The gene expression data was retrieved from GEO database (GSE169077) and differentially expressed genes were determined using GEO2R tool based on adjusted P-value and log2FC values based on which 27 genes were found to be significant of which 9 genes were up-regulated while 18 genes were down-regulated. This was followed by gene enrichment analysis identify the related Gene Ontology terms and pathways. The Protein-Protein Network is constructed with 15 nodes and 22 edges using STRING database and then exported to Cytoscape 3.8 to predict the hub genes. The hub genes identified are POSTN, COL1A2, COL1A1, BMP1, MXRA5, MMP13 and SERPINF1. The hub genes identified were not only found to be associated with bone related disorders but also few were involved in other diseases. Therefore targeting these genes for disease management could be a viable option in osteoarthritis.

Prediction of potential prognostic biomarkers in metastatic prostate cancer based on a circular RNA-mediated competing endogenous RNA regulatory network.

Recently, studies on competing endogenous RNA (ceRNA) networks have become prevalent, and circular RNAs (circRNAs) have crucial implications for the development and progression of carcinoma. However, studies relevant to metastatic prostate cancer (mPCa) are scant. This study aims to discover potential ceRNAs that may be related to the prognosis of mPCa. RNA-Seq data were obtained from the MiOncoCirc database and Gene Expression Omnibus (GEO). Differential expression patterns of RNAs were examined using R packages. Circular RNA Interactome, miRTarBase, miRDB and TargetScan were applied to predict the corresponding relation between circRNAs, miRNAs and mRNAs. The Gene Ontology (GO) annotations were performed to present related GO terms, and Gene Set Enrichment Analysis (GSEA) tools were applied for pathway annotations. Moreover, survival analysis was conducted for the hub genes. We found 820 circRNAs, 81 miRNAs and 179 mRNAs that were distinguishingly expressed between primary prostate cancer (PCa) and mPCa samples. A ceRNA network including 45 circRNAs, 24 miRNAs and 56 mRNAs was constructed. In addition, the protein–protein interaction (PPI) network was built, and 10 hub genes were selected by using the CytoHubba application. Among the 10 hub genes, survival analysis showed that ITGA1, LMOD1, MYH11, MYLK, SORBS1 and TGFBR3 were significantly connected with disease-free survival (DFS). The circRNA-mediated ceRNA network provides potential prognostic biomarkers for metastatic prostate cancer.

Systematical Screening of Intracellular Protein Targets of Polyphemusin-I Using Escherichia coli Proteome Microarrays.

With their wide repertoire of mechanisms, antimicrobial peptides (AMPs) are promising alternatives to fight against varied pathogenic microorganisms (bacteria, fungi, viruses, parasites, etc.). AMPs, novel components of the innate immune defense system, are secreted by all organisms. The aquatic environment represents a huge population and an enormous source of varied AMPs. Polyphemusin-I, a marine AMP isolated from hemocytes of an American horseshoe crab, possesses high antimicrobial activities. Studies on polyphemusin-I have verified the intracellular mechanisms of action, however, its intracellular targets are not yet explored. In this study, we employed Escherichia coli proteome microarrays to systematically screen the entire intracellular protein targets of polyphemusin-I. A total of 97 protein targets of polyphemusin-I were statistically analyzed from the quadruplicate Escherichia coli proteome microarrays assays. Among these identified protein targets, 56 proteins had cellular location inside the cell (i.e., cytoplasm), one in the plasma membrane, one in the periplasm and the rest 39 proteins had no specified cellular location. The bioinformatics analysis of these identified protein targets of polyphemusin-I in gene ontology (GO) enrichment category of molecular function revealed significant enrichment in nucleic acid related GO terms i.e., “RNA binding”, “nucleotide binding”, “nuclease activities”, “uracil DNA N-glycosylase activities” and others. Moreover, enrichment in GO category of biological process also depicted enrichment in nucleic acid related GO terms, such as “nucleic acid phosphodiester bond hydrolysis”, “deoxyribonucleotide metabolism”, and others. In accordance to GO enrichment analysis, protein families (PFAM) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analysis also showed significant enrichment in nucleic acid terms. These enrichment results suggest that polyphemusin-I targets nucleic acid-associated proteins. Furthermore, to provide a comprehensive study, we compared the identified protein targets of polyphemusin-I with previously identified protein targets of four AMPs (P-Der, Lfcin B, PR-39, and Bac 7) using Escherichia coli proteome microarrays. The comparison study of five AMPs (polyhemusin-I, P-Der, Lfcin B, PR-39, and Bac 7) showed only nine common protein targets in all the five AMPs, whereas a total of 39 and 43 common protein targets were identified among the two marine AMPs (polyphemusin-I and P-Der) and three terrestrial AMPs (Lfcin B, PR-39 and Bac7), respectively. To further reveal the target pattern of marine and terrestrial AMPs, the enrichment results obtained from common protein targets of marine AMPs with terrestrial AMPs were compared. The comparison result indicated that AMPs have unique mechanism of action among marine or terrestrial AMPs. Hence, in this study, we have not only identified the intracellular protein targets of polyphemusin-I, but also revealed the protein target differences between marine AMPs and terrestrial AMPs.

Identification of Candidate Circular RNAs Underlying Intramuscular Fat Content in the Donkey.

Intramuscular fat (IMF) content is a crucial indicator of meat quality. Circular RNAs (circRNAs) are a large class of endogenous RNAs that are involved in many physiological processes. However, the expression and function of circRNA in IMF in the donkey remains unresolved. Here we performed an expression profiling of circRNAs in the donkey longissimus dorsi muscle and identified 12,727 candidate circRNAs. Among these, 70% were derived from the exons of protein genes. Furthermore, a total of 127 differentially expressed (DE) circRNAs were identified in high (H) and low (L) IMF content groups, including 63 upregulated and 64 downregulated circRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the host genes of the DE circRNAs showed that the host genes were enriched in lipid metabolism related GO terms (e.g., fatty acid beta-oxidation using acyl-CoA dehydrogenase and MLL3/4 complex), and signaling pathways (e.g., TGF-beta and lysine degradation signaling pathway). Further analyses indicated that 127 DE circRNAs were predicted to potentially interact with miRNAs, leading to the construction of circRNA-miRNA regulatory network. Multiple circRNAs can potentially function as sponges of miRNAs that regulate the differentiation of adipocytes. Our results provide valuable expression profile information for circRNA in the donkey and new insight into the regulatory mechanisms of circRNAs in the regulation of IMF content.

ROP18-Mediated Transcriptional Reprogramming of HEK293T Cell Reveals New Roles of ROP18 in the Interplay Between Toxoplasma gondii and the Host Cell.

Toxoplasma gondii secretes a number of virulence-related effector proteins, such as the rhoptry protein 18 (ROP18). To further broaden our understanding of the molecular functions of ROP18, we examined the transcriptional response of human embryonic kidney cells (HEK293T) to ROP18 of type I T. gondii RH strain. Using RNA-sequencing, we compared the transcriptome of ROP18-expressing HEK293T cells to control HEK293T cells. Our analysis revealed that ROP18 altered the expression of 750 genes (467 upregulated genes and 283 downregulated genes) in HEK293T cells. Gene ontology (GO) and pathway enrichment analyses showed that differentially expressed genes (DEGs) were significantly enriched in extracellular matrix– and immune–related GO terms and pathways. KEGG pathway enrichment analysis revealed that DEGs were involved in several disease-related pathways, such as nervous system diseases and eye disease. ROP18 significantly increased the alternative splicing pattern “retained intron” and altered the expression of 144 transcription factors (TFs). These results provide new insight into how ROP18 may influence biological processes in the host cells via altering the expression of genes, TFs, and pathways. More in vitro and in vivo studies are required to substantiate these findings.

Transcriptome analysis reveals potential immune function-related regulatory genes/pathways of female Lubo goat submandibular glands at different developmental stages.

BackgroundThe submandibular glands, as major salivary glands, participate in rumen digestion in goats. Sialic acid, lysozyme, immunoglobulin A (IgA), lactoferrin and other biologically active substances secreted in the submandibular glands were reported in succession, which suggests that the submandibular gland may have immune functions in addition to participating in digestion. The aim of this study was to map the expression profile of differentially expressed genes (DEGs) at three different stages by transcriptome sequencing, screen immune-related genes and pathways by bioinformatics methods, and predict the immune function of submandibular glands at different developmental stages.MethodsNine submandibular gland tissue samples were collected from groups of 1-month-old kids, 12-month-old adolescent goats and 24-month-old adult goats (3 samples from each group), and high-throughput transcriptome sequencing was conducted on these samples. The DEGs among the three stages were screened and analysed. Key genes and signalling pathways were selected via protein-protein interaction (PPI) network analysis.ResultsThe results revealed 2,706, 2,525 and 52 DEGs between 1-month-old and 12-month-old goats, between 1-month-old and 24-month-old goats, and between 12-month-old and 24-month-old goats, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated that most of the DEGs were enriched in immune- related GO terms and pathways. Based on functional enrichment analysis and network analysis, 10 genes (PTPRC, CD28, SELL, LCP2, MYC, LCK, ZAP70, ITGB2, SYK and CCR7), two signalling pathways (the T cell receptor signalling pathway and the NF-κβ signalling pathway) and eight GO terms (T cell receptor signalling pathway, neutrophil mediated immunity, B cell mediated immunity, regulation of alpha-beta T cell activation, positive regulation of T cell proliferation, regulation of leukocyte differentiation, positive regulation of antigen receptor-mediated signalling pathway, positive regulation of lymphocyte proliferation) that may play key roles in the immune functions of the goat submandibular glands at different developmental stages were identified. Moreover, we found that eight antibacterial peptide-encoding genes were downregulated in the tuberculosis and salivary secretion pathways, while all immunoglobulins were upregulated in 10 immune system pathways. These findings indicate that the submandibular glands may be important immunological organs during the growth process of goats and that the immune function of these glands gradually weakens with age up to 12 months but remains relatively stable after 12 months of age. Overall, this study will improve our understanding of transcriptional regulation related to goat submandibular gland immune function.

Integrated small RNA and mRNA expression profiles reveal miRNAs and their target genes in response to Aspergillus flavus growth in peanut seeds

BackgroundMicroRNAs are important gene expression regulators in plants immune system. Aspergillus flavus is the most common causal agents of aflatoxin contamination in peanuts, but information on the function of miRNA in peanut-A. flavus interaction is lacking. In this study, the resistant cultivar (GT-C20) and susceptible cultivar (Tifrunner) were used to investigate regulatory roles of miRNAs in response to A. flavus growth.ResultsA total of 30 miRNAs, 447 genes and 21 potential miRNA/mRNA pairs were differentially expressed significantly when treated with A. flavus. A total of 62 miRNAs, 451 genes and 44 potential miRNA/mRNA pairs exhibited differential expression profiles between two peanut varieties. Gene Ontology (GO) analysis showed that metabolic-process related GO terms were enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses further supported the GO results, in which many enriched pathways were related with biosynthesis and metabolism, such as biosynthesis of secondary metabolites and metabolic pathways. Correlation analysis of small RNA, transcriptome and degradome indicated that miR156/SPL pairs might regulate the accumulation of flavonoids in resistant and susceptible genotypes. The miR482/2118 family might regulate NBS-LRR gene which had the higher expression level in resistant genotype. These results provided useful information for further understanding the roles of miR156/157/SPL and miR482/2118/NBS-LRR pairs.ConclusionsIntegration analysis of the transcriptome, miRNAome and degradome of resistant and susceptible peanut varieties were performed in this study. The knowledge gained will help to understand the roles of miRNAs of peanut in response to A. flavus.

Enhanced Immune Responses with Serum Proteomic Analysis of Hu Sheep to Foot-and-Mouth Disease Vaccine Emulsified in a Vegetable Oil Adjuvant.

Our previous study demonstrated that a vegetable oil consisting of soybean oil, vitamin E, and ginseng saponins (SO-VE-GS) had an adjuvant effect on a foot-and-mouth disease (FMD) vaccine in a mouse model. The present study was to compare the adjuvant effects of SO-VE-GS and the conventional ISA 206 on an FMD vaccine in Hu sheep. Animals were intramuscularly (i.m.) immunized twice at a 3-week interval with 1 mL of an FMD vaccine adjuvanted with SO-VE-GS (n = 10) or ISA 206 (n = 9). Animals without immunization served as control (n = 10). Blood was sampled prior to vaccination and at 2, 4, 6, and 8 weeks post the booster immunization to detect FMD virus (FMDV)-specific IgG. Blood collected at 8 weeks after the booster was used for the analyses of IgG1 and IgG2, serum neutralizing (SN) antibody, IL-4 and IFN-γ production, and proteomic profiles. The results showed that IgG titers rose above the protection level (1:128) in SO-VE-GS and ISA 206 groups after 2 and 4 weeks post the booster immunization. At 6 weeks post the booster, the ISA 206 group had 1 animal with IgG titer less than 1:128 while all the animals in the SO-VE-GS group retained IgG titers of more than 1:128. At 8 weeks post the booster, 6 of 9 animals had IgG titers less than 1:128 with a protective rate of 33.3% in the ISA 206 group, while only 1 of 10 animals had IgG titer less than 1:128 with a protective rate of 90% in the SO-VE-GS group, with statistical significance. In addition, IgG1, IgG2, SN antibodies, IL-4, and IFN-γ in the SO-VE-GS group were significantly higher than those of the ISA 206 group. Different adjuvant effects of SO-VE-GS and ISA 206 may be explained by the different proteomic profiles in the two groups. There were 39 and 47 differentially expressed proteins (DEPs) identified in SO-VE-GS compared to the control or ISA 206 groups, respectively. In SO-VE-GS vs. control, 3 immune related gene ontology (GO) terms and 8 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were detected, while 2 immune related GO terms and 5 KEGG pathways were found in ISA 206 vs. control. GO and KEGG analyses indicated that ‘positive regulation of cytokine secretion’, ‘Th1/Th2 cell differentiation’, and ‘Toll-like receptor signaling pathways’, were obviously enriched in the SO-VE-GS group compared to the other groups. Coupled with protein–protein interaction (PPI) analysis, we found that B7TJ15 (MAPK14) was a key DEP for SO-VE-GS to activate the immune responses in Hu sheep. Therefore, SO-VE-GS might be a promising adjuvant for an FMD vaccine in Hu sheep.

A 17-gene expression-based prognostic signature associated with the prognosis of patients with breast cancer: A STROBE-compliant study.

Identification of reliable predictive biomarkers for patients with breast cancer (BC).Univariate Cox proportional hazards regression model was conducted to identify genes correlated with the overall survival (OS) of patients in the TCGA-BRCA cohort. Functional enrichment analysis was conducted to investigate the biological meaning of these survival related genes. Then, patients in TCGA-BCRA were randomly divided into training set and test. Least absolute shrinkage and selection operator (LASSO) penalized Cox regression model was performed and the risk score of BC patients in this model was used to build a prognostic signature. The prognostic performance of the signature was evaluated in the training set, test set, and an independent validation set GSE7390.2519 genes were demonstrated to be significantly associated with the OS of BC patients. Functional annotation of the 2519 genes suggested that these genes were associated with immune response and protein synthesis related gene ontology terms and pathways. 17 genes were identified in the LASSO Cox regression model and used to construct a 17-gene signature. Patients in the 17-gene signature low risk group have better OS and event-free survival compared with those in the 17-gene signature high risk group in the TCGA-BRCA cohort. The prognostic role of the 17-gene signature has been confirmed in the validation cohort. Multivariable Cox proportional hazards regression model suggested the 17-gene signature was an independent prognostic factor in BC.The 17-gene signature we developed could successfully classify patients into high- and low-risk groups, indicating that it might serve as candidate biomarker in BC.

A genome-wide multiphenotypic association analysis identified candidate genes and gene ontology shared by four common risky behaviors.

Background: Risky behaviors can lead to huge economic and health losses. However, limited efforts are paid to explore the genetic mechanisms of risky behaviors.Result: MASH analysis identified a group of target genes for risky behaviors, such as APBB2, MAPT and DCC. For GO enrichment analysis, FUMA detected multiple risky behaviors related GO terms and brain related diseases, such as regulation of neuron differentiation (adjusted P value = 2.84×10-5), autism spectrum disorder (adjusted P value =1.81×10-27) and intelligence (adjusted P value =5.89×10-15).Conclusion: We reported multiple candidate genes and GO terms shared by the four risky behaviors, providing novel clues for understanding the genetic mechanism of risky behaviors.Methods: Multivariate Adaptive Shrinkage (MASH) analysis was first applied to the GWAS data of four specific risky behaviors (automobile speeding, drinks per week, ever-smoker, number of sexual partners) to detect the common genetic variants shared by the four risky behaviors. Utilizing genomic functional annotation data of SNPs, the SNPs detected by MASH were then mapped to target genes. Finally, gene set enrichment analysis of the identified candidate genes were conducted by the FUMA platform to obtain risky behaviors related gene ontology (GO) terms as well as diseases and traits, respectively.

Identification of gene modules and hub genes in colon adenocarcinoma associated with pathological stage based on WGCNA analysis

Colorectal cancer (CRC) is one of the leading causes of cancer-related mortality in the world, in which colon adenocarcinoma (COAD) is the most common histological subtype of CRC. In this study, our aim is to identify gene modules and representative candidate biomarkers for clinical prognosis of patients with COAD, and help to predict prognosis and reveal the mechanisms of cancer progression. Weighted gene co-expression network analysis (WGCNA) was performed to construct a co-expression network and identify gene modules correlated with TNM clinical staging of COAD patients. The Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed with the module gene. Protein-protein interaction (PPI) network and hub gene identification were explored with Cytoscape software. Finally, the hub gene mRNA level was validated in Oncomine database. Five gene modules, related with the pathological TNM stage, were constructed, and the gene module was enriched in cell proliferation, invasion and migration related GO terms and metabolic related KEGG pathways. A total of top 10 hub genes was identified, and in which six of the hub genes show a significant up-regulation in COAD as compared to normal tissue, including IVL, KRT16, KRT6C, KRT6A, KRT78 and SBSN. In conclusion, we identified five gene modules and six candidate biomarkers correlated with the TNM staging of COAD patients. These findings may help us to understand the tumor progression of COAD and provide prognostic biomarkers as well as therapeutic targets.