Gene Ontology Categories Research Articles

Comparative high-throughput analysis of sperm membrane proteins from crossbred bulls with contrasting fertility.

The aim of the present study was to identify fertility associated sperm membrane proteins in crossbred bulls. Sperm membrane proteins from high- and low-fertile Holstein Friesian crossbred bulls (n=3 each) were subjected to high-throughput liquid chromatography-mass spectrometry (LC-MS/MS) for comparative proteomic analysis. Proteomic profiling identified a total of 456 proteins in crossbred bull spermatozoa; it was found that 108 proteins were up regulated while 26 proteins were down regulated (>1.5-folds) in spermatozoa from low- compared to high-fertile bulls. Gene ontology classification revealed that upregulated proteins in low-fertile bulls were involved in biological process such as oxidation-reduction process (p=3.14E-06), fusion of sperm to egg plasma membrane (p=7.51E-04), sperm motility (p=0.03), and capacitation (p=0.09), while down regulated proteins were associated with transport (p=6.94E-04), superoxide metabolic process (p=0.02), and tricarboxylic acid cycle (p=0.04). KEGG pathway analysis revealed that oxidative phosphorylation and tricarboxylic acid cycle pathways are the most significantly affected pathway in low-fertile bulls. It was concluded that expression of proteins associated with oxidative phosphorylation and tricarboxylic acid cycle pathways were altered in low-fertile crossbred bulls, and expression levels of SPATA19, ELSPBP1, ACRBP, CLU, SUCLA2, and SPATC1 could aid in assessing potential fertility of crossbred bulls.

Transcriptome Analysis of the Cf-13-Mediated Hypersensitive Response of Tomato to Cladosporium fulvum Infection.

Tomato leaf mold disease caused by Cladosporium fulvum (C. fulvum) is one of the most common diseases affecting greenhouse tomato production. Cf proteins can recognize corresponding AVR proteins produced by C. fulvum, and Cf genes are associated with leaf mold resistance. Given that there are many physiological races of C. fulvum and that these races rapidly mutate, resistance to common Cf genes (such as Cf-2, Cf-4, Cf-5, and Cf-9) has decreased. In the field, Ont7813 plants (carrying the Cf-13 gene) show effective resistance to C. fulvum; thus, these plants could be used as new, disease-resistant materials. To explore the mechanism of the Cf-13-mediated resistance response, transcriptome sequencing was performed on three replicates each of Ont7813 (Cf-13) and Moneymaker (MM; carrying the Cf-0 gene) at 0, 9, and 15 days after inoculation (dai) for a total of 18 samples. In total, 943 genes were differentially expressed, specifically in the Ont7813 response process as compared to the Moneymaker response process. Gene ontology (GO) classification of these 943 differentially expressed genes (DEGs) showed that GO terms, including “hydrogen peroxide metabolic process (GO_Process)”, “secondary active transmembrane transporter activity (GO_Function)”, and “mismatch repair complex (GO_Component)”, which were the same as 11 other GO terms, were significantly enriched. An analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that many key regulatory genes of the Cf-13-mediated resistance response processes were involved in the “plant hormone signal transduction” pathway, the “plant–pathogen interaction” pathway, and the “MAPK signaling pathway–plant” pathway. Moreover, during C. fulvum infection, jasmonic acid (JA) and salicylic acid (SA) contents significantly increased in Ont7813 at the early stage. These results lay a vital foundation for further understanding the molecular mechanism of the Cf-13 gene in response to C. fulvum infection.

Brominated Carbazole with Antibiotic Adjuvant Activity Displays Pleiotropic Effects in MRSA's Transcriptome.

Methicillin-resistant Staphylococcus aureus (MRSA) is a major threat to human health, as the US mortality rate outweighs those from HIV, tuberculosis, and viral hepatitis combined. In the wake of the COVID-19 pandemic, antibiotic-resistant bacterial infections acquired during hospital stays have increased. Antibiotic adjuvants are a key strategy to combat these bacteria. We have evaluated several small molecule antibiotic adjuvants that have strong potentiation with β-lactam antibiotics and are likely inhibiting a master regulatory kinase, Stk1. Here, we investigated how the lead adjuvant (compound 8) exerts its effects in a more comprehensive manner. We hypothesized that the expression levels of key resistance genes would decrease once cotreated with oxacillin and the adjuvant. Furthermore, bioinformatic analyses would reveal biochemical pathways enriched in differentially expressed genes. RNA-seq analysis showed 176 and 233 genes significantly up- and downregulated, respectively, in response to cotreatment. Gene ontology categories and biochemical pathways that were significantly enriched with downregulated genes involved carbohydrate utilization, such as the citrate cycle and the phosphotransferase system. One of the most populated pathways was S. aureus infection. Results from an interaction network constructed with affected gene products supported the hypothesis that Stk1 is a target of compound 8. This study revealed a dramatic impact of our lead adjuvant on the transcriptome that is consistent with a pleiotropic effect due to Stk1 inhibition. These results point to this antibiotic adjuvant having potential broad therapeutic use in combatting MRSA.

Newly developed sugarbeet lines with altered postharvest respiration rates differ in transcription factor and glycolytic enzyme expression

AbstractRespiration is the principal cause for postharvest sucrose loss in sugarbeet (Beta vulgaris L.) roots. Although reductions in respiration rate could mitigate these losses, developing sugarbeet cultivars and storage procedures that reduce respiration are hindered by a lack of knowledge of the genetic and metabolic factors that control storage respiration rate. Research was conducted to identify genes and gene products that affect storage respiration rate by creating two sugarbeet lines that differ in respiration rate and characterizing gene expression differences between these lines. Sugarbeet lines F1056 and F1057, which differ by up to 42% in respiration rate, were created by divergent selection of a sugarbeet population using root respiration rate after 30 d in storage as the principal selection criterion. RNA sequencing identified 287 differentially expressed genes (DEGs) between these lines on the day of harvest and after 28‐d storage. Of these DEGs, nine encoded transcription factors and five encoded enzymes involved in the respiratory pathway. Other DEGs contributed to a variety of biological and molecular functions based on gene ontology classifications. Of respiratory pathway DEGS, genes for NAD+‐ and NADP+‐dependent forms of glyceraldehyde‐3‐phosphate dehydrogenase were of note due to their high upregulation in the high respiring line and for their established role in glycolysis, a pathway identified as a likely bottleneck in respiratory substrate production. Overall, lines F1056 and F1057 provide new tools for investigating genetic and physiological differences in storage respiration rate, and their DEGs identify candidates for genes affecting sugarbeet root respiration rate.

Identification of host cell proteins that interact with the M protein of Pigeon paramyxovirus type 1

Pigeon paramyxovirus type 1 (PPMV-1) belongs to the avian paramyxovirus type 1 group of viruses, which can cause tremors, torticollis, and respiratory signs in domestic and wild pigeons. The M protein of PPMV-1 is a multifunctional structural protein. It not only helps in the assembly, budding, and positioning of the virus but also inhibits the host's immune response and promotes replication of the virus in the host. In this study, the GST pull-down method was used to screen host proteins that interact with PPMV-1 M protein, and then mass spectrometry (MS) was used to analyse the screened host proteins. Enrichment analysis of the differentially expressed genes showed that the 77 screened proteins were highly associated with the gene ontology categories: protein synthesis, metabolism, and cell signalling pathway transduction. We selected NIMA-related kinase 7 (NEK7) as the candidate protein for co-localization analysis and co-immunoprecipitation verification. The results revealed that PPMV-1 M protein interacts with NEK7 of the host cell. This interactome study of PPMV-1 M protein will serve to clarify its function during viral replication and will provide a crucial theoretical basis for studying the pathogenic mechanism of PPMV-1.

MiRNA Profiling and Its Role in Multi-Omics Regulatory Networks Connected with Somaclonal Variation in Cucumber (Cucumis sativus L.).

The role of miRNAs in connection with the phenomenon of somaclonal variation, which occurs during plant in vitro culture, remains uncertain. This study aims to investigate the possible role of miRNAs in multi-omics regulatory pathways in cucumber somaclonal lines. For this purpose, we performed sRNA sequencing (sRNA-seq) from cucumber fruit samples identified 8, 10 and 44 miRNAs that are differentially expressed between somaclones (S1, S2, S3 lines) and the reference B10 line of Cucumis sativus. For miRNA identification, we use ShortStack software designed to filter miRNAs from sRNAs according to specific program criteria. The identification of predicted in-silico targets revealed 2,886 mRNAs encoded by 644 genes. The functional annotation of miRNA’s target genes and gene ontology classification revealed their association with metabolic processes, response to stress, multicellular organism development, biosynthetic process and catalytic activity. We checked with bioinformatic analyses for possible interactions at the level of target proteins, differentially expressed genes (DEGs) and genes affected by genomic polymorphisms. We assume that miRNAs can indirectly influence molecular networks and play a role in many different regulatory pathways, leading to somaclonal variation. This regulation is supposed to occur through the process of the target gene cleavage or translation inhibition, which in turn affects the proteome, as we have shown in the example of molecular networks. This is a new approach combining levels from DNA-seq through mRNA-seq, sRNA-seq and in silico PPI in the area of plants’ somaclonal variation.

New Insight on Vitality Differences for the Penaeid Shrimp, Fenneropenaeus chinensis, in Low Salinity Environment Through Transcriptomics

Excessive rainfall changes salinity in shrimp farming ponds in short period and exerts low salinity stress on the outdoor breeding shrimp under global warming. Fenneropenaeus chinensis can have different performance on vitality in low salinity environments. To reveal mechanisms of vitality difference in shrimp living in low saline environments. This study based on the normal and moribund F. chinensis in 10 ppt salinity environment using high-throughput sequencing identifies 1,429 differentially expressed genes (DEGs), 586 of which are upregulated, while 843 of which are downregulated in the normal group (FCN10) as compared to the moribund group (FCM10). Meanwhile, another transcriptomic analysis is conducted on the normal and moribund shrimp from 25 ppt (FCN25 vs. FCM25) salinity environment as the control, in which 1,311 DEGs (upregulated: 327 genes, downregulated: 984 genes) are identified. In this study, intersective pathways, GO (Gene Ontology) categories and DEGs from the two groups of comparative transcriptome are investigated. The two intersective pathways (Metabolism of xenobiotics by cytochrome P450, Pentose, and glucuronate interconversions) significantly enriched by DEGs are related to detoxification. In these two pathways, there is one vitality regulation-related gene (VRRG), the Dhdh (dihydrodiol dehydrogenase), which is upregulated in both the groups of FCN10 and FCN25 as compared to the groups of FCM10 and FCM25, respectively. Similarly, in the 25 top intersective GO categories, four VRRGs are revealed. Three of them are upregulated (Itgbl, kielin/chordin-like protein, Slc2a8, solute carrier family 2, facilitated glucose transporter member 8-like protein and Cyp3a30, cytochrome P450 3A30-like protein); one of them is downregulated (Slc6a9, sodium-dependent nutrient amino acid transporter 1-like protein isoform X2). These GO categories are related to transmembrane transporter activity of substance, enzyme inhibitor activity, monooxygenase activity. RT-qPCR analysis further verifies the VRRGs. The study gives new insight into understanding the vitality differences for F. chinensis, in low salinity environment. The pathways and DEGs in response to low salinity stress in modulating the vitality of F. chinensis that could serve as tools in future genetic studies and molecular breeding.

Predicting embryonic aneuploidy rate in IVF patients using whole-exome sequencing.

Infertility is a major reproductive health issue that affects about 12% of women of reproductive age in the United States. Aneuploidy in eggs accounts for a significant proportion of early miscarriage and in vitro fertilization failure. Recent studies have shown that genetic variants in several genes affect chromosome segregation fidelity and predispose women to a higher incidence of egg aneuploidy. However, the exact genetic causes of aneuploid egg production remain unclear, making it difficult to diagnose infertility based on individual genetic variants in mother's genome. In this study, we evaluated machine learning-based classifiers for predicting the embryonic aneuploidy risk in female IVF patients using whole-exome sequencing data. Using two exome datasets, we obtained an area under the receiver operating curve of 0.77 and 0.68, respectively. High precision could be traded off for high specificity in classifying patients by selecting different prediction score cutoffs. For example, a strict prediction score cutoff of 0.7 identified 29% of patients as high-risk with 94% precision. In addition, we identified MCM5, FGGY, and DDX60L as potential aneuploidy risk genes that contribute the most to the predictive power of the model. These candidate genes and their molecular interaction partners are enriched for meiotic-related gene ontology categories and pathways, such as microtubule organizing center and DNA recombination. In summary, we demonstrate that sequencing data can be mined to predict patients' aneuploidy risk thus improving clinical diagnosis. The candidate genes and pathways we identified are promising targets for future aneuploidy studies.

Transcriptomic analysis of methyl jasmonate treatment reveals gene networks involved in drought tolerance in pearl millet

Water deficit stress at the early stage of development is one of the main factors limiting pearl millet production. One practice to counteract this limitation would be to resort to the application of hormones to stimulate plant growth and development at critical stages. Exogenous methyl jasmonate (MeJA) can improve drought tolerance by modulating signaling, metabolism, and photosynthesis pathways, therefore, we assumed that can occur in pearl millet during the early stage of development. To decipher the molecular mechanisms controlling these pathways, RNAseq was conducted in two pearl millet genotypes, drought-sensitive SosatC88 and drought-tolerant Souna3, in response to 200 μM of MeJA. Pairwise comparison between the MeJA-treated and non-treated plants revealed 3270 differentially expressed genes (DEGs) among 20,783 transcripts in SosatC88 and 127 DEGs out of 20,496 transcripts in Souna3. Gene ontology (GO) classification assigned most regulated DEGs in SosatC88 to heme binding, oxidation–reduction process, response to oxidative stress and membrane, and in Souna3 to terpene synthase activity, lyase activity, magnesium ion binding, and thylakoid. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis reveals that DEGs in SosatC88 are related to the oxidation–reduction process, the biosynthesis of other secondary metabolites, the signal transduction, and the metabolism of terpenoids, while in Souna3, DEGs are related to the metabolism of terpenoids and the energy metabolism. Two genes encoding a diterpenoid biosynthesis-related (Pgl_GLEAN_10009413) and a Glutathione S transferase T3 (Pgl_GLEAN_10034098) were contra-regulated between SosatC88 and Souna3. Additionally, five random genes differentially expressed by RNAseq were validated using qPCR, therefore, they are potential targets for the development of novel strategies breeding schemes for plant growth under water deficit stress. These insights into the molecular mechanisms of pearl millet genotype tolerance at the early stage of development contribute to the understanding of the role of hormones in adaptation to drought-prone environments.

Identification and Integrated Analysis of the miRNA-mRNA Regulatory Network in Lens from an H2O2-Induced Zebrafish Cataract Model

Purpose This study aimed to explore the regulatory mechanisms of age-related cataract (ARC) formation. Methods Cataracts in zebrafish were induced by injecting hydrogen peroxide into the fish anterior chamber. The mRNA and miRNA expression profiles of the lens from H2O2-injected and PBS-injected zebrafishes were detected by RNA sequencing. The LIMMA package was applied to identify differentially expressed genes (DEGs). Gene Ontology categories were enriched by the R “cluster Profiler” package and Kyoto Encyclopedia of Genes and Genomes pathway enrichment was performed based on hypergeometric distribution using the R “phyper” function. The protein-protein interaction network of DEGs was built via the STRING. Target genes of differentially expressed miRNAs (DEmiRs) were predicted by miRanda. Furthermore, DEGs were selected as DEmiR targets and a DEmiR-DEG regulatory network was constructed via Cytoscape. Results In total, 3689 DEGs (such as opn1mw4, LOC103908930, si:dkeyp-1h4.8, crispld1b, cyp1a, and gdpd3a) including 2478 upregulated and 1211 downregulated genes were identified. 177 DEmiRs (such as dre-miR-96-3p, dre-miR-182-5p, dre-miR-9-7-3p, and dre-miR-124-4-5p) including 108 upregulated and 69 downregulated miRNAs were detected. The DEGs are involved in cell death, DNA repair, and cell development-related pathways. A protein-protein interaction network including 79 node genes was constructed to explore the interactions of DEGs. Furthermore, a DEmiR-DEG regulatory network focusing on the DNA repair process was constructed, including 21 hub DEGs and 15 hub DEmiRs. Conclusions We identified several DEGs and constructed a miRNA-mRNA regulatory network related to the DNA repair process in a zebrafish cataract model. These genes participate in the oxidative stress response of lens epithelium cells and finally contribute to the formation of zebrafish cataracts. The hub DEGs and hub DEmiRs could be potential therapeutic targets for ARC.

Network Theoretical Approach to Explore Factors Affecting Signal Propagation and Stability in Dementia's Protein-Protein Interaction Network.

Dementia—a syndrome affecting human cognition—is a major public health concern given to its rising prevalence worldwide. Though multiple research studies have analyzed disorders such as Alzheimer’s disease and Frontotemporal dementia using a systems biology approach, a similar approach to dementia syndrome as a whole is required. In this study, we try to find the high-impact core regulating processes and factors involved in dementia’s protein–protein interaction network. We also explore various aspects related to its stability and signal propagation. Using gene interaction databases such as STRING and GeneMANIA, a principal dementia network (PDN) consisting of 881 genes and 59,085 interactions was achieved. It was assortative in nature with hierarchical, scale-free topology enriched in various gene ontology (GO) categories and KEGG pathways, such as negative and positive regulation of apoptotic processes, macroautophagy, aging, response to drug, protein binding, etc. Using a clustering algorithm (Louvain method of modularity maximization) iteratively, we found a number of communities at different levels of hierarchy in PDN consisting of 95 “motif-localized hubs”, out of which, 7 were present at deepest level and hence were key regulators (KRs) of PDN (HSP90AA1, HSP90AB1, EGFR, FYN, JUN, CELF2 and CTNNA3). In order to explore aspects of network’s resilience, a knockout (of motif-localized hubs) experiment was carried out. It changed the network’s topology from a hierarchal scale-free topology to scale-free, where independent clusters exhibited greater control. Additionally, network experiments on interaction of druggable genome and motif-localized hubs were carried out where UBC, EGFR, APP, CTNNB1, NTRK1, FN1, HSP90AA1, MDM2, VCP, CTNNA1 and GRB2 were identified as hubs in the resultant network (RN). We finally concluded that stability and resilience of PDN highly relies on motif-localized hubs (especially those present at deeper levels), making them important therapeutic intervention candidates. HSP90AA1, involved in heat shock response (and its master regulator, i.e., HSF1), and EGFR are most important genes in pathology of dementia apart from KRs, given their presence as KRs as well as hubs in RN.

Non-negative matrix factorization and differential expression analyses identify hub genes linked to progression and prognosis of glioblastoma multiforme

One of the most prevailing primary brain tumors in adult human male is glioblastoma multiforme (GBM), which is categorized by rapid cellular growth. Even though the combination therapy comprises surgery, chemotherapy, and adjuvant therapies, the survival rate, on average, is 14.6 months. Glioma stem cells (GSCs) have key roles in tumorigenesis, progression, and defiance against chemotherapy and radiotherapy. In our study, firstly, the gene expression dataset GSE124145 was retrieved; the non-negative matrix factorization (NMF) method was applied on GBM dataset, and differentially expressed genes analysis (DEGs) was performed. After which, overlapping genes between metagenes and DEGs were detected to examine the Gene Ontology (GO) categories in the biological process (BP) in the stemness of GBM. The common hub genes were used to construct protein–protein interaction (PPI) network and further GO, while Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway was utilized to pinpoint the real hub genes. The analysis of hub genes particular for the same GO categories demonstrated that specific hub genes triggered distinct features of the same biological processes. After utilizing GSE124145 and The Cancer Genome Atlas (TCGA) dataset for survival analysis, we screened five real hub genes: GUCA1A, RFC2, GNG11, MMP19, and NRG1, which are strongly associated with the progression and prognosis of GBM. The DEGs analysis revealed that all real hub genes were overexpressed in GBM and TCGA datasets, which further validates our results. The constructed study of PPI, GO, and KEGG pathway on common hub genes was performed. Finally, the KEGG pathways performed on the top 15 candidate hub genes (including six real hub genes) of the PPI network in the GBM gene expression dataset study found mitogen-activated protein kinase (Mapk) signaling pathway to be the most significant pathway. The rest of the hub genes reviewed throughout the analysis might be favorable targets for diagnosing and treating GBM and lower-grade gliomas.

Abstract P5-12-07: Proteomic profiling of extracellular vesicles released from leptin-treated breast cancer cells: A potential role in cancer metabolism

Abstract Extracellular vesicles (EVs), nanosized particles enclosed within a phospholipid bilayer membrane, are recognized key determinants of cell-to-cell communication in breast cancer, able to modulate phenotypic behavior of recipient cells by transferring their molecular cargo. Recently, we identified the adipokine leptin, whose circulating levels are elevated in obesity, as an important regulator of EV biogenesis and release in mammary carcinoma cells. Here, we focused on the identification of a specific leptin-induced EV proteomic signature in attempt to find molecular effectors associated with breast cancer progression. To this end, we analyzed protein content of EVs, isolated by ultracentrifugation method from conditioned media of ER-α positive MCF-7 breast cancer cells grown in the presence of leptin. EVs were fully characterized using Transmission Electron Microscopy (TEM), Nanoparticle Tracking Analysis (NTA) and Immunoblotting. Proteomic analysis of EV cargo, obtained by liquid chromatography-tandem mass spectrometry analysis (LC-MS/MS), revealed an altered protein profile, identifying 274 up-regulated proteins (FC>1.5) in EVs from leptin-treated cells compared to EV from untreated ones. Gene set enrichment analysis (GSEA) was used to explore Gene Ontology categories enriched in the entire list of up-regulated proteins. We identified a number of significantly enriched biologic processes mainly involved in energy production and mitochondrial metabolism, including “generation of precursor metabolites and energy” (GO:0006091), “fatty acid metabolic process” (GO:0006631), “mitochondrial respiratory chain” (GO:0033108), “mitochondrial gene expression” (GO:0140053) and “mitochondrial transport” (GO:0006839). Expression of up-regulated proteins, in both cell and EV lysates, was validated by immunoblotting. Protein-protein interaction network analysis revealed CYC1 (Cytochrome c-1), NDUFV2 (NADH:ubiquinone oxidoreductase core subunit V2), and NDUFAB1 (NADH:ubiquinone oxidoreductase subunit AB1), important subunits of mitochondrial complexes and enhancer of mitochondrial metabolism, as the most interactive proteins among that up-regulated in EVs released by leptin-treated breast cancer cells. This suggests a role for leptin-treated EVs in enhancing oxidative phosphorylation in both epithelial breast cancer cells and their surrounding microenvironment. Collectively, our data indicate that the adipokine leptin could induce in breast cancer cells the release of EVs enriched in specific proteins mainly involved in mitochondrial metabolism. Since, mitochondrial function is essential for supplying energy to support breast cancer cell growth and progression, understanding the metabolic mechanisms mediated by EV released by breast cancer cells may provide important clues to develop novel therapeutic approaches for treatment of breast cancer especially in obesity setting. Citation Format: Luca Gelsomino, Giusi La Camera, Ines Barone, Salvatore Panza, Giovanna Morello, Amanda Caruso, Chiara Chiodo, Daniela Bonofiglio, Cinzia Giordano, Sebastiano Andò, Stefania Catalano. Proteomic profiling of extracellular vesicles released from leptin-treated breast cancer cells: A potential role in cancer metabolism [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P5-12-07.

SPDL1 Is an Independent Predictor of Patient Outcome in Colorectal Cancer.

Spindle Apparatus Coiled-Coil Protein 1 (SPDL1) is a relatively recently identified coiled-coil domain containing protein and an important determinant of DNA fidelity by ensuring faithful mitosis. Hence, SPDL1 is suspected to underlie genomic (in-)stability in human cancers, yet its exact roles in these diseases remain largely underexplored. Given that genomic instability (GIN) is a crucial feature in colorectal cancer (CRC), we primarily asked whether the expression of this protein may account for differences in clinicopathological features and survival rates of CRC patients. Protein expression was evaluated by immunohistochemistry in the institutional tissue microarray (TMA), and gene expression by the analysis of publicly available datasets. To place the prognostic relevance in a predicted biological context, gene co-expression set around SPDL1 identified by public data mining was annotated and assessed for enrichment in gene ontology (GO) categories, BRITE hierarchies, and Reactome pathways. The comparison with adjacent normal tissue revealed a high expression of SPDL1 protein in a subset of tumor cases (48.84%), and these had better prognosis than the SPDL1-low expression counterpart even after adjustment for multiple confounders. SPDL1-high expression within tumors was associated with a median 56-month survival advantage, but not with any clinicopathological characteristics of our cohort. In the TCGA cohort, SPDL1 was overexpressed in tumor tissue and positively associated with improved survival, chromosome instability phenotype, and various GIN markers. In addition to the genes critically involved in the cell cycle and mitosis, a gene set co-expressed with SPDL1 contained checkpoint members of both chromosome segregation and DNA replication, as well as those associated with defective DNA repair, and retrograde vesicle-mediated transport. In conclusion, SPDL1 is an independent predictor of CRC patient survival in a possible connection with chromosomal instability.

Berberine Inhibits FOXM1 Dependent Transcriptional Regulation of POLE2 and Interferes With the Survival of Lung Adenocarcinoma.

Background: Berberine is one of the most interesting and promising natural anticancer drugs. POLE2 is involved in many cellular functions such as DNA replication and is highly expressed in a variety of cancers. However, the specific molecular mechanism of berberine interfering with POLE2 expression in lung adenocarcinoma (LUAD) is still unknown to a great extent. Method: The KEGG database (Release 91.0) and Gene Ontology (GO) category database were used for functional annotation of differentially expressed genes after berberine treatment. Reproducibility assessment using TCGA dataset. The biological functions of berberine in LUAD were investigated by a series of in vitro and in vivo experiments: MTT, colony formation, mouse xenograft and plasmid transfection. The molecular mechanisms of berberine were demonstrated by plasmid transfection, quantitative RT-PCR and Western blotting. Result: The elevated expression of FOXM1 and the high enrichment of DNA replication pathway were confirmed in LUAD by microarray and TCGA analysis, and were positively correlated with poor prognosis. Functionally, berberine inhibited the proliferation and survival of LUAD cell lines in vitro and in vivo. Mechanistically, berberine treatment down regulated the expression of FOXM1which closely related to survival, survival related genes in Cell cycle and DNA replication pathway, and significantly down regulated the expression of survival related POLE2. Interestingly, we found that the transcription factor FOXM1 could act as a bridge between berberine and POLE2. Conclusion: Berberine significantly inhibited LUAD progression via the FOXM1/POLE2, and FOXM1/POLE2 may act as a clinical prognostic factor and a therapeutic target for LUAD. Berberine may be used as a promising therapeutic candidate for LUAD patients.

Population genomics provides insights into lineage divergence and local adaptation within the cotton bollworm.

The cotton bollworm Helicoverpa armigera is a cosmopolitan pest and its diverse habitats plausibly contribute to the formation of diverse lineages. Despite the significant threat it poses to economic crops worldwide, its evolutionary history and genetic basis of local adaptation are poorly understood. In this study, we de novo assembled a high-quality chromosome-level reference genome of H. a. armigera (contig N50=7.34Mb), with 99.13% of the HaSCD2 assembly assigned to 31 chromosomes (Z-chromosome + 30 autosomes). We constructed an ultradense variation map across 14 cotton bollworm populations and identified a novel lineage in northwestern China. Historical inference showed that effective population size changes coincided with global temperature fluctuation. We identified nine differentiated genes in the three H. armigera lineages (H. a. armigera, H. a. conferta and the new northwestern Chinese lineage), of which per and clk genes are involved in circadian rhythm. Selective sweep analyses identified a series of Gene Ontology categories related to climate adaptation, feeding behaviour and insecticide tolerance. Our findings reveal fundamental knowledge of the local adaptation of different cotton bollworm lineages and will guide the formulation of cotton bollworm management measures at different scales.

Construction and Analysis of the lncRNA-miRNA-mRNA Network Based on Competing Endogenous RNA in Atrial Fibrillation.

BackgroundAccumulated studies have revealed that long non-coding RNAs (lncRNAs) play critical roles in human diseases by acting as competing endogenous RNAs (ceRNAs). However, functional roles and regulatory mechanisms of lncRNA-mediated ceRNA in atrial fibrillation (AF) remain unknown. In the present study, we aimed to construct the lncRNA-miRNA-mRNA network based on ceRNA theory in AF by using bioinformatic analyses of public datasets.MethodsMicroarray data sets of GSE115574 and GSE79768 from the Gene Expression Omnibus database were downloaded. Twenty-one AF right atrial appendage (RAA) samples and 22 sinus rhythm (SR) subjects RAA samples were selected for subsequent analyses. After merging all microarray data and adjusting for batch effect, differentially expressed genes were identified. Gene Ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were carried out. A ceRNA network was constructed.ResultA total of 8 lncRNAs and 43 mRNAs were significantly differentially expressed with fold change >1.5 (p < 0.05) in RAA samples of AF patients when compared with SR. GO and KEGG pathway analysis showed that cardiac muscle contraction pathway were involved in AF development. The ceRNA was predicted by co-expressing LOC101928304/ LRRC2 from the constructional network analysis, which was competitively combined with miR-490-3p. The expression of LOC101928304 and LRRC were up-regulated in myocardial tissue of patients with AF, while miR-490-3p was down-regulated.ConclusionWe constructed the LOC101928304/miR-490-3p/LRRC2 network based on ceRNA theory in AF in the bioinformatic analyses of public datasets. The ceRNA network found from this study may help improve our understanding of lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of AF.

Analysis of Global Gene Expression in Maize (Zea mays) Vegetative and Reproductive Tissues That Differ in Accumulation of Starch and Sucrose.

Carbon allocation between vegetative and reproductive tissues impacts cereal grain production. Despite great agricultural importance, sink–source relationships have not been fully characterized at the early reproductive stages in maize. Here, we quantify the accumulation of non-structural carbohydrates and patterns of gene expression in the top internode of the stem and the female inflorescence of maize at the onset of grain filling (reproductive stage R1). Top internode stem and female inflorescence tissues of the Puma maize inbred line were collected at reproductive stage R1 (without pollination) and non-structural carbohydrates were quantified by spectrophotometry. The female inflorescence accumulated starch at higher levels than the top internode of the stem. Global mRNA transcript levels were then evaluated in both tissues by RNA sequencing. Gene expression analysis identified 491 genes differentially expressed between the female inflorescence and the top stem internode. Gene ontology classification of differentially expressed genes showed enrichment for sucrose synthesis, the light-dependent reactions of photosynthesis, and transmembrane transporters. Our results suggest that sugar transporters play a key role in sugar partitioning in the maize stem and reveal previously uncharacterized differences between the female inflorescence and the top internode of the stem at early reproductive stages.

Transcriptome analysis of grass carp (Ctenopharyngodon idella) and Holland's spinibarbel (Spinibarbus hollandi) infected with Ichthyophthirius multifiliis

Ichthyophthirius multifiliis is a protozoan ciliate that causes white spot disease (also known as ichthyophthiriasis) in freshwater fish. Holland's spinibarbel (Spinibarbus hollandi) was less susceptible to white spot disease than grass carp (Ctenopharyngodon Idella). In this study, grass carp and Holland's spinibarbel are infected by I. multifiliis and the amount of infection is 10,000 theronts per fish. All grass carp died within 12 days after infection, and the survival rate of Holland's spinibarbel was more than 80%. In order to study the difference in sensitivity of these two fish species to I. multifiliis, transcriptome analysis was conducted using gill, skin, liver, spleen and head kidney of Holland's spinibarbel and grass carp at 48 h post-infection with I. multifiliis. A total of 489,296,696 clean reads were obtained by sequencing. A total of 105 significantly up-regulated immune-related genes were obtained by Gene Ontology (GO) classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis in grass carp. Cluster of differentiation 40 (CD40), cluster of differentiation 80 (CD 80), tumor necrosis factor-alpha (TNF-α), toll-like receptor 4 (TLR-4), interleukin 1 beta (IL-1β) and other inflammatory-related genes in grass carp were enriched in the cytokine-cytokine receptor interaction pathway and toll-like receptor pathway. In Holland's spinibarbel, a total of 46 significantly up-regulated immune-related genes were obtained by GO classification and KEGG pathway enrichment analysis. Immune-related genes, such as Immunoglobin heavy chain (IgH), cathepsin S (CTSS), complement C1q A chain (C1qA), complement component 3 (C3) and complement component (C9) were enriched in phagosome pathway, lysosome pathway and complement and coagulation concatenation pathway. C3 was significantly up-regulated in gill and head kidney. Fluorescence in situ hybridization (FISH) showed that the C3 gene was highly expressed in gill tissue of Holland's spinibarbel infected with I. multifiliis. A small amount of C3 gene was expressed in the gill arch of grass carp after infected with I. multifiliis. In conclusion, the severe inflammatory response in vivo after infecting grass carp with I. multifiliis might be the main cause of the death of grass carp. The extrahepatic expression of the gene of Holland's spinibarbel might play an important role in the immune defense against I. multifiliis.

Transcriptional alterations in bladder epithelial cells in response to infection with different morphological states of uropathogenic Escherichia coli

Uropathogenic Escherichia coli (UPEC) may undergo a cyclic cascade of morphological alterations that are believed to enhance the potential of UPEC to evade host responses and re-infect host cell. However, knowledge on the pathogenic potential and host activation properties of UPEC during the morphological switch is limited. Microarray analysis was performed on mRNA isolated from human bladder epithelial cells (HBEP) after exposure to three different morphological states of UPEC (normal coliform, filamentous form and reverted form). Cells stimulated with filamentous bacteria showed the lowest number of significant gene alterations, although the number of enriched gene ontology classes was high suggesting diverse effects on many different classes of host genes. The normal coliform was in general superior in stimulating transcriptional activity in HBEP cells compared to the filamentous and reverted form. Top-scored gene entities activated by all three morphological states included IL17C, TNFAIP6, TNF, IL20, CXCL2, CXCL3, IL6 and CXCL8. The number of significantly changed canonical pathways was lower in HBEP cells stimulated with the reverted form (32 pathways), than in cells stimulated with the coliform (83 pathways) or filamentous bacteria (138 pathways). A host cell invasion assay showed that filamentous bacteria were unable to invade bladder cells, and that the number of intracellular bacteria was markedly lower in cells infected with the reverted form compared to the coliform. In conclusion, the morphological state of UPEC has major impact on the host bladder response both when evaluating the number and the identity of altered host genes and pathways.