Multiple Gene Ontology Terms Research Articles

Unlocking the molecular modifications of plasma-activated water-induced oxidation through redox proteomics: In the case of duck myofibrillar protein (Anas platyrhynchos)

Plasma-activated water (PAW) contains multiple active species that alter the structure of myofibrillar protein (MP) to enhance their gel properties. This work investigated the impact of PAW on the oxidation of cysteine in MP by label-free quantitative proteomics. PAW treatment caused the oxidation of 8241 cysteine sites on 2815 proteins, and structural proteins such as nebulin, myosin XVIIIB, myosin XVIIIA, and myosin heavy chain were susceptible to oxidation by PAW. Bioinformatics analysis, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, subcellular localization, and STRING analysis, indicated that these proteins with differential oxidation sites were mainly derived from the cytoplasm and membrane, and were involved in multiple GO terms and KEGG pathways. This is one of the first reports of the redox proteomic changes induced by PAW treatment, and the results are useful for understanding the possible mechanism of PAW-induced oxidation of MP.

The transcriptome of acute dehydration in myeloid leukemic cells.

Human myeloid leukemia cells (HL-60/S4) exposed to hyperosmotic stress with sucrose undergo dehydration and cell shrinkage. Interphase chromatin and mitotic chromosomes congeal, exhibiting altered phase separation (demixing) of chromatin proteins. To investigate changes in the transcriptome, we exposed HL-60/S4 cells to hyperosmotic sucrose stress (~600 milliOsmolar) for 30 and 60 minutes. We employed RNA-Seq of polyA mRNA to identify genes with increased or decreased transcript levels relative to untreated control cells (i.e., differential gene expression). These genes were examined for over-representation of Gene Ontology (GO) terms. In stressed cells, multiple GO terms associated with transcription, translation, mitochondrial function and proteosome activity, as well as "replication-dependent histones", were over-represented among genes with increased transcript levels; whereas, genes with decreased transcript levels were over-represented with transcription repressors. The transcriptome profiles of hyperosmotically-stressed cells suggest acquisition of cellular rebuilding, a futile homeostatic response, as these cells are ultimately doomed to a dehydrated death.

Comparative analysis of whole-transcriptome RNA expression of lung tissue of Chinese soft-shell turtle infected by Trionyx sinensis Hemorrhagic Syndrome Virus

Trionyx sinensis Hemorrhagic Syndrome Virus (TSHSV), the first aquatic arterivirus identified in China, causes severe mortality to T. sinensis. In this study, we sought to determine the functions of T. sinensis mRNAs and non-coding RNAs (ncRNAs) that were differentially expressed (DE) over different periods of TSHSV infection of T. sinensis lung. We used RT-qPCR to validate the sequencing results of select RNAs, confirming their reliable and referable nature. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were used to predict multiple biological functions and signaling pathways in various comparison groups (1-day versus mock, 3-day versus 1-day, and 5-day versus 3-day). Multiple types of differentially expressed RNA, including mRNA, lncRNA, circRNA, and miRNA, were associated with cardiac dysfunction, coagulation abnormalities, and arachidonic acid metabolism at day 1. Pre-inflammatory cytokines and inflammatory factors such as PLA2G4A, cPLA2, γ-GGT1, TNFRSF14, TCP11L2, PTER CYP2J2 and LTC4S, were noticeably regulated at the same time. On day 3, multiple GO terms and KEGG pathways were implicated, including those related to virus defense, apoptosis, pyroptosis, and inflammatory response. Notably, key genes such as RSAD2, TRIM39, STAT4, CASP1, CASP14, MYD88, CXCL3, CARD11, ZBP1, and ROBO4 exhibited significant regulation. The lncRNAs and circRNAs that targeted the genes involved in viral recognition (TLR5), apoptosis (CARD11), pyroptosis (ZBP1), inflammatory processes (NEK7, RASGRP4, and SELE) and angiogenesis (ROBO4) exhibited significant regulation. Significantly regulated miRNAs were primarily linked to genes involved in apoptosis (Let-7f-3p, miR-1260a, miR-455–3p), and inflammation (miR-146a, miR-125a, miR-17a, miR-301b, and miR-30a-3p). The findings could advance our understanding of the host immunological response to TSHSV and offer new ideas for developing effective strategies to prevent infection of T. sinensis.

Abstract 6606: Biomarker discovery in non-small-cell lung cancer enabled by deep multi-omics profiling of proteins, metabolites, transcripts, and genes in blood

Abstract Lung cancer is the leading cause of cancer-related deaths in the United States, with estimates of 236,740 new cases and 118,830 deaths in 2022 secondary to the disease. Blood-based liquid biopsies hold promise to reduce morbidity and mortality from lung cancer by enabling early detection to downstage disease at diagnosis, theragnostic identification of patients most likely to be helped or harmed by therapy, monitoring of therapeutic efficacy, and detection of residual disease. PrognomiQ’s multi-omics platform comprehensively profiles proteins, metabolites, lipids, mRNA, and cfDNA in blood samples which can be used for the development of liquid biopsy tests with high sensitivity and specificity for lung cancer. We conducted a case-control study comprising 1031 subjects: 361 subjects with untreated non-small-cell lung cancer (NSCLC) and 670 matched controls which included 340 subjects with salient pulmonary and gastrointestinal co-morbidities. Blood samples from each subject were processed to provide 7 different `omics readouts. LCMS was used to detect and quantify proteins, metabolites, and lipids. In addition, cfDNA and mRNA were assayed using next-generation sequencing. cfDNA reads were analyzed to estimate fragment-lengths, copy-number variation, and CpG site methylation. All molecular data were normalized using standard methods specific to each assay. Univariate analyses of cases vs controls were performed to identify differentially abundant features on all available samples per assay. We detected 9,868 proteins, 605 lipids, 329 metabolites, and 109,070 mRNA transcripts. Of these, 3,098 proteins, 210 lipids, 57 metabolites, and 30,236 mRNA transcripts were significantly different (FWER < 0.05) in cases versus controls. Gene set enrichment analysis on statistically significant transcripts and proteins identified multiple gene-ontology terms associated with cancer including the Wnt signaling process and IgA immunoglobulin complex, respectively. From cfDNA data, we identified 234 non-contiguous genomic regions associated with the fragment-length disorder, 4,790 with copy-number variation, and 74 differentially methylated genomic regions spanning 184 CpG sites (FWER < 0.05). With the premise that deviations from copy number neutrality are more likely to indicate a tumor contribution, we then focused our examination on those differentially expressed proteins that overlap with differentially expressed mRNA transcripts as well as CNV genomic regions. We identified 52 protein coding genes including E-cadherin (associated with EMT) and related binding proteins such as RAB11B, CAPZB, EPS15, FLNB, MYH9, STK24 and YWHAE. Ongoing machine-learning-based classifier training to distinguish between cancer and non-cancer can serve as the basis for the development of high-sensitivity liquid-biopsy tests for lung cancer. Citation Format: Jinlyung Choi, Ajinkya Kokate, Ehdieh Khaledian, Manway Liu, Preethi Prasad, John Blume, Jessica Chan, Rea Cuaresma, Kevin Dai, Manoj Khadka, Thidar Khin, Yuya Kodama, Joon-Yong Lee, Hoda Malekpour, Megan Mora, Nithya Mudaliar, Sara Nouri Golmaei, Madhuvanthi Ramaiah, Saividya Ramaswamy, Peter Spiro, Dijana Vitko, Kavya Swaminathan, James Yee, Brian Young, Chinmay Belthangady, Bruce Wilcox, Brian Koh, Philip Ma. Biomarker discovery in non-small-cell lung cancer enabled by deep multi-omics profiling of proteins, metabolites, transcripts, and genes in blood. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6606.

High-throughput transcriptome sequencing reveals the key stages of cardiovascular development in zebrafish embryos

BackgroundThe cardiovascular developmental process is a tightly regulated network involving multiple genes. The current understanding of the molecular mechanism behind cardiovascular development is insufficient and requires further research.ResultsTranscriptome sequencing of three developmental stages in zebrafish embryos was performed and revealed three key cardiovascular developmental stages. Then, the differentially expressed genes (DEGs) involved in cardiovascular development were screened out. The three developmental stages were 18 (T1), 24 (T2), and 42 h post fertilization (hpf) (T3), and the three stages were confirmed by detecting differences in expression between cardiomyocyte and endothelial marker genes (cmlc2, fli1) using in situ hybridization, which represents the characteristics of cardiovascular development. Thousands of DEGs were identified using transcriptome analysis. Of them, 2605 DEGs were in T1-vs-T2, including 2003 up-regulated and 602 down-regulated genes, 6446 DEGs were in T1-vs-T3, consisting of 4608 up-regulated and 1838 down-regulated genes, and 3275 DEGs were in T2-vs-T3, including 2420 up-regulated and 855 down-regulated genes. There were 644 common DEGs and 167 common five-fold higher differentially expressed genes (HDEGs) identified, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID). Significant differences was observed in the levels of gene expression among different developmental stages in multiple GO terms and KEGG pathways, such as cell migration to the midline involved in heart development, cardiovascular system development, circulatory system process for biological processes of GO terms; and cardiac muscle contraction, adrenergic signaling in cardiomyocytes for KEGG pathways. These results demonstrated that these three stages were important period for the development of the cardiovascular system. Lastly, we used quantitative real-time PCR (qPCR) to validate the reliability of RNA-sequencing by selecting 21 DEGs.ConclusionsThese results demonstrated that these three stages represented the important periods for cardiovascular system development of zebrafish and some candidate genes was obtained and provided a solid foundation for additional functional studies of the DEGs.

Protease-activated receptor-2 dependent and independent responses of bone cells to prostate cancer cell secretory products.

Metastatic prostate cancer lesions in the skeleton are frequently characterized by excessive formation of bone. Prostate cancer cells secrete factors, including serine proteases, that are capable of influencing the behavior of surrounding cells. Some of these proteases activate protease-activated receptor-2 (PAR2 ), which is expressed by osteoblasts (bone-forming cells) and precursors of osteoclasts (bone-resorbing cells). The aim of the current study was to investigate a possible role for PAR2 in regulating the behavior of bone cells exposed to metastatic prostate cancer cells. The effect of medium conditioned by the PC3, DU145, and MDA-PCa-2b prostate cancer cell lines was investigated in assays of bone cell function using cells isolated from wildtype and PAR2 -null mice. Osteoclast differentiation was assessed by counting tartrate-resistant acid phosphatase-positive multinucleate cells in bone marrow cultured in osteoclastogenic medium. Osteoblasts were isolated from calvariae of neonatal mice, and BrdU incorporation was used to assess their proliferation. Assays of alkaline phosphatase activity and quantitative PCR analysis of osteoblastic gene expression were used to assess osteoblast differentiation. Responses of osteoblasts to medium conditioned by MDA-PCa-2b cells were analyzed by RNAseq. Conditioned medium (CM) from all three cell lines inhibited osteoclast differentiation independently of PAR2 . Media from PC3 and DU145 cells had no effect on assays of osteoblast function. Medium conditioned by MDA-PCa-2b cells stimulated BrdU incorporation in both wildtype and PAR2 -null osteoblasts but increased alkaline phosphatase activity and Runx2 and Col1a1 expression in wildtype but not PAR2 -null cells. Functional enrichment analysis of RNAseq data identified enrichment of multiple gene ontology terms associated with lysosomal function in both wildtype and PAR2 -null cells in response to MDA-PCa-2b-CM. Analysis of individual genes identified osteogenesis-associated genes that were either upregulated by MDA-PCa-2b-CM selectively in wildtype cells or downregulated selectively in PAR2 -null cells. Factors secreted by prostate cancer cells influence bone cell behavior through both PAR2 -dependent and -independent mechanisms. Both PAR2 -independent suppression of osteoclast differentiation and PAR2 -dependent stimulation of osteogenesis are likely to determine the nature of prostate cancer metastases in bone.

Investigating the Transcriptome of Candida albicans in a Dual-Species Staphylococcus aureus Biofilm Model.

Candida albicans is an opportunistic pathogen found throughout multiple body sites and is frequently co-isolated from infections of the respiratory tract and oral cavity with Staphylococcus aureus. Herein we present the first report of the effects that S. aureus elicits on the C. albicans transcriptome. Dual-species biofilms containing S. aureus and C. albicans mutants defective in ALS3 or ECE1 were optimised and characterised, followed by transcriptional profiling of C. albicans by RNA-sequencing (RNA-seq). Altered phenotypes in C. albicans mutants revealed specific interaction profiles between fungus and bacteria. The major adhesion and virulence proteins Als3 and Ece1, respectively, were found to have substantial effects on the Candida transcriptome in early and mature biofilms. Despite this, deletion of ECE1 did not adversely affect biofilm formation or the ability of S. aureus to interact with C. albicans hyphae. Upregulated genes in dual-species biofilms corresponded to multiple gene ontology terms, including those attributed to virulence, biofilm formation and protein binding such as ACE2 and multiple heat-shock protein genes. This shows that S. aureus pushes C. albicans towards a more virulent genotype, helping us to understand the driving forces behind the increased severity of C. albicans-S. aureus infections.

The effect and mechanism of cypermethrin-induced hippocampal neurotoxicity as determined by network pharmacology analysis and experimental validation

ABSTRACT Cypermethrin (CMN) is a widely used artificial synthetic pesticide that causes neurotoxicity in the hippocampus. However, the underlying toxicological targets and mechanisms remain unclear. In this study, network pharmacology analysis and in vitro models were integrated to investigate the effect and mechanism of CMN-induced hippocampal neurotoxicity. A total of 88 targets of CMN-induced hippocampal neurotoxicity were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes enrichment (KEGG) analyses suggested that these targets were related to multiple GO terms and signaling pathways. To further investigate underlying mechanism, the top 10 hub targets (Akt1, Tnf, Ptgs2, Casp3, Igf1, Sirt1, Jun, Cat, Il10, and Bcl2l1) were screened. Furthermore, cell viability and lactate dehydrogenase (LDH) assays demonstrated that CMN was toxic to HT22 cells in a time- and dose-dependent manner. Terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) staining revealed that treatment with CMN increased the proportion of apoptotic cells. In addition, the real-time quantitative polymerase chain reaction (RT-qPCR) results indicated that CMN altered the mRNA expression levels of most of the hub targets, with the exceptions of Igf1 and Jun. The results demonstrated that multiple targets and signaling pathways were involved in CMN-induced hippocampal neurotoxicity. These findings provided reference values for subsequent studies of the toxicological mechanism of CMN.

Integrating network pharmacology and in vitro model to investigate hippocampal neurotoxicity induced by atrazine

Atrazine (ATR), a commonly applied herbicide in agriculture, has been found to cause hippocampal injury in rodents. However, the underlying toxicological targets and mechanisms are unclear. In this study, network pharmacology analysis and in vitro model were integrated to investigate the effect and mechanism of ATR-induced hippocampal neurotoxicity. In total, 71 targets of hippocampal neurotoxicity induced by ATR were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes enrichment (KEGG) enrichment analysis suggested that these targets were related to multiple GO terms and signaling pathways. To further investigate the underlying mechanisms, the top 10 hub targets were screened and included tumor protein p53 (Tp53), caspase 3 (Casp3), prostaglandin-endoperoxide synthase 2 (Ptgs2), cAMP-responsive element-binding protein 1 (Creb1), estrogen receptor 1 (Esr1), Jun proto-oncogene (Jun), brain-derived neurotrophic factor (Bdnf), catalase (Cat), sirtuin 1 (Sirt1) and Fos proto-oncogene (Fos). Moreover, the cell counting kit-8 (CCK8) and lactate dehydrogenase (LDH) assay showed that ATR had time and dose-dependent cytotoxicity on H19-7 cells. TUNEL staining revealed that ATR increased the apoptotic ratio. In addition, Real-time quantitative polymerase chain reaction (RT-qPCR) results indicated that the mRNA expression levels of all hub targets showed significant changes, except Esr1 and Jun. Our study demonstrated that ATR mainly acted on multiple targets and signaling pathways to exert its hippocampal neurotoxicity. These results provided initial evidence for the further exploration of the toxicological mechanism of ATR.

Cadherin and Wnt signaling pathways as key regulators in diabetic nephropathy.

AimA recent meta-analysis of genome-wide linkage studies (GWLS) has identified multiple genetic regions suggestive of linkage with DN harboring hundreds of genes. Moving this number of genetic loci forward into biological insight is truly the next step. Here, we approach this challenge with a gene ontology (GO) analysis in order to yield biological and functional role to the genes, an over-representation test to find which GO terms are enriched in the gene list, pathway analysis, as well as protein network analysis.MethodGO analysis was performed using protein analysis through evolutionary relationships (PANTHER) version 14.0 software and P-values less than 0.05 were considered statistically significant. GO analysis was followed by over-representation test for the identification of enriched terms. Statistical significance was calculated by Fisher’s exact test and adjusted using the false discovery rate (FDR) for correction of multiple tests. Cytoscape with the relevant plugins was used for the construction of the protein network and clustering analysis.ResultsThe GO analysis assign multiple GO terms to the genes regarding the molecular function, the biological process and the cellular component, protein class and pathway analysis. The findings of the over-representation test highlight the contribution of cell adhesion regarding the biological process, integral components of plasma membrane regarding the cellular component, chemokines and cytokines with regard to protein class, while the pathway analysis emphasizes the contribution of Wnt and cadherin signaling pathways.ConclusionsOur results suggest that a core feature of the pathogenesis of DN may be a disturbance in Wnt and cadherin signaling pathways, whereas the contribution of chemokines and cytokines need to be studied in additional studies.

A Metagenome-Wide Association Study of Gut Microbiome in Patients With Multiple Sclerosis Revealed Novel Disease Pathology.

While microbiome plays key roles in the etiology of multiple sclerosis (MS), its mechanism remains elusive. Here, we conducted a comprehensive metagenome-wide association study (MWAS) of the relapsing-remitting MS gut microbiome (ncase = 26, ncontrol = 77) in the Japanese population, by using whole-genome shotgun sequencing. Our MWAS consisted of three major bioinformatic analytic pipelines (phylogenetic analysis, functional gene analysis, and pathway analysis). Phylogenetic case-control association tests showed discrepancies of eight clades, most of which were related to the immune system (false discovery rate [FDR] < 0.10; e.g., Erysipelatoclostridium_sp. and Gemella morbillorum). Gene association tests found an increased abundance of one putative dehydrogenase gene (Clo1100_2356) and one ABC transporter related gene (Mahau_1952) in the MS metagenome compared with controls (FDR < 0.1). Molecular pathway analysis of the microbiome gene case-control comparisons identified enrichment of multiple Gene Ontology terms, with the most significant enrichment on cell outer membrane (P = 1.5 × 10−7). Interaction between the metagenome and host genome was identified by comparing biological pathway enrichment between the MS MWAS and the MS genome-wide association study (GWAS) results (i.e., MWAS-GWAS interaction). No apparent discrepancies in alpha or beta diversities of metagenome were found between MS cases and controls. Our shotgun sequencing-based MWAS highlights novel characteristics of the MS gut microbiome and its interaction with host genome, which contributes to our understanding of the microbiome’s role in MS pathophysiology.

Cardiac Regeneration and Tumor Growth-What Do They Have in Common?

Acute myocardial infarction is a leading cause of death. Unlike most adult mammals, zebrafish have the capability to almost fully regenerate their hearts after injury. In contrast, ischemic damage in adult human and mouse hearts usually results in scar tissue. mRNA-Sequencing (Seq) and miRNA-Seq analyses of heart regeneration in zebrafish over time showed that the process can be divided into three phases: the first phase represents dedifferentiation and proliferation of cells, the second phase is characterized by migration, and in the third phase cell signals indicate heart development and differentiation. The first two phases seem to share major similarities with tumor development and growth. To gain more insight into these similarities between cardiac regeneration and tumor development and growth, we used patient matched tumor normal (“healthy”) RNA-Seq data for several tumor entities from The Cancer Genome Atlas (TCGA). Subsequently, RNA data were processed using the same pipeline for both the zebrafish samples and tumor datasets. Functional analysis showed that multiple Gene Ontology terms (GO terms) are involved in both early stage cardiac regeneration and tumor development/growth across multiple tumor entities. These GO terms are mostly associated with cell cycle processes. Further analysis showed that orthologous genes are the same key players that regulated these changes in both diseases. We also observed that GO terms associated with heart development in the third late phase of cardiac regeneration are downregulated in the tumor entities. Taken together, our analysis illustrates similarities between cardiac remodeling and tumor progression.

Progesterone-induced RNA Hand2os1 regulates decidualization in mice uteri.

Decidualization is a critical process for successful embryo implantation and subsequent placenta formation. The characterization and physiological function of lncRNA during decidualization remain largely unknown. In the present study, we conducted RNA-sequencing analysis to compare gene expression between decidua of days 6 and 8, and normal pregnant endometrium (day 4). A total of 2332 high-confidence putative lncRNA transcripts were expressed. Functional clustering analysis of cis and trans lncRNA targets showed that differentially expressed lncRNAs may regulate multiple gene ontology terms and pathways that have important functions in decidualization. Subsequent analyses using qRT-PCR validated that eight of all lncRNAs were differentially regulated in mice uteri during decidualization, both in vivo and in vitro. Furthermore, we showed that differentially expressed lncRNA of Hand2os1 was specifically detected in stromal cells on days 2 to 5 of pregnancy and was strongly upregulated in decidual cells on days 6-8 of pregnancy. Similarly, Hand2os1 expression was also strongly expressed in decidualized cells following artificial decidualization, both in vivo and in vitro. In uterine stromal cells, P4 was able to significantly upregulate the expression of Hand2os1, but upregulation was impeded by RU486, whereas E2 appeared to have no regulating effect on Hand2os1 expression. Concurrently, Hand2os1 significantly promoted the decidual process in vitro and dramatically increased decidualization markers Prl8a2 and Prl3c1. Our results provide a valuable catalog for better understanding of the functional roles of lncRNAs in pregnant mouse uteri, as it relates to decidualization.

Detection and characterization of microRNA expression profiling and its target genes in response to canine parvovirus in Crandell Reese Feline Kidney cells.

BackgroundMicroRNAs (miRNAs) play an essential role in gene regulators in many biological and molecular phenomena. Unraveling the involvement of miRNA as a key cellular factor during in vitro canine parvovirus (CPV) infection may facilitate the discovery of potential intervention candidates. However, the examination of miRNA expression profiles in CPV in tissue culture systems has not been fully elucidated.MethodIn the present study, we utilized high-throughput small RNA-seq (sRNA-seq) technology to investigate the altered miRNA profiling in miRNA libraries from uninfected (Control) and CPV-2c infected Crandell Reese Feline Kidney cells.ResultsWe identified five of known miRNAs (miR-222-5p, miR-365-2-5p, miR-1247-3p, miR-322-5p and miR-361-3p) and three novel miRNAs (Novel 137, Novel 141 and Novel 102) by sRNA-seq with differentially expressed genes in the miRNA repertoire of CPV-infected cells over control. We further predicted the potential target genes of the aforementioned miRNAs using sequence homology algorithms. Notably, the targets of miR-1247-3p exhibited a potential function associated with cellular defense and humoral response to CPV. To extend the probing scheme for gene targets of miR-1247-3p, we explored and performed Gene Ontology (GO) enrichment analysis of its target genes. We discovered 229 putative targets from a total of 38 enriched GO terms. The top over-represented GO enrichment in biological process were lymphocyte activation and differentiation, marginal zone B cell differentiation, negative regulation of cytokine production, negative regulation of programed cell death, and negative regulation of signaling. We next constructed a GO biological process network composed of 28 target genes of miR-1247-3p, of which, some genes, namely BCL6, DLL1, GATA3, IL6, LEF1, LFNG and WNT1 were among the genes with obviously intersected in multiple GO terms.ConclusionThe miRNA-1247-3p and its cognate target genes suggested their great potential as novel therapeutic targets or diagnostic biomarkers of CPV or other related viruses.

NetGO: improving large-scale protein function prediction with massive network information.

Automated function prediction (AFP) of proteins is of great significance in biology. AFP can be regarded as a problem of the large-scale multi-label classification where a protein can be associated with multiple gene ontology terms as its labels. Based on our GOLabeler—a state-of-the-art method for the third critical assessment of functional annotation (CAFA3), in this paper we propose NetGO, a web server that is able to further improve the performance of the large-scale AFP by incorporating massive protein-protein network information. Specifically, the advantages of NetGO are threefold in using network information: (i) NetGO relies on a powerful learning to rank framework from machine learning to effectively integrate both sequence and network information of proteins; (ii) NetGO uses the massive network information of all species (>2000) in STRING (other than only some specific species) and (iii) NetGO still can use network information to annotate a protein by homology transfer, even if it is not contained in STRING. Separating training and testing data with the same time-delayed settings of CAFA, we comprehensively examined the performance of NetGO. Experimental results have clearly demonstrated that NetGO significantly outperforms GOLabeler and other competing methods. The NetGO web server is freely available at http://issubmission.sjtu.edu.cn/netgo/.

Identification of more objective biomarkers for Blood-Stasis syndrome diagnosis.

BackgroundBlood-stasis syndrome (BSS) is one of the Traditional Chinese medicine (TCM) syndrome differentiations that are commonly seen in stroke and ischemic heart diseases; however, the BSS differentiation criterion is not standardized. More objective biomarkers for BSS diagnosis are needed.MethodsAcute ischemic stroke (AIS) or unstable angina (UA) patients with BSS and healthy controls were enrolled. The miRNA and mRNA expression profiles of UA patients and AIS patients were compared to those of healthy controls to identify the differentially expressed miRNA and mRNA of BSS. Bioinformatics analysis was used to identify significantly deregulated miRNAs and mRNAs correlated to BSS. QRT-PCR was performed to validate the bioinformatics analysis results.ResultsApproximately 401 mRNAs and 11 miRNAs were differentially expressed in both UA and AIS patients compared to healthy controls. Gene ontology (GO) functional analysis was performed, and multiple GO terms were enriched. Among the overlapping DE miRNAs and mRNAs, miR-146b-5p, -199a-5p and 23 targeted mRNAs were pivotal genes in the BSS genomic characteristics. These 2 miRNAs and 23 mRNAs formed network-type biomarkers for BSS.ConclusionsThe genomic characteristics of BSS were shown in this study. miR-146b-5p, -199a-5p and the 23 targeted mRNAs formed a diagnostic network for BSS. Further improvement and validation of this diagnostic network might lead to more objective diagnostic criteria for BSS.

Identification of candidate causal genes and their associated pathogenic mechanisms underlying teratozoospermia based on the spermatozoa transcript profiles.

Teratozoospermia with unclear pathomechanism is one of the common causes for failed fertilisation. This study aimed to further explore the pathological mechanism for teratozoospermia. Spermatozoal transcript profiles generated from 13 normal fertile men and eight infertile males with a consistent severe heterogeneous teratozoospermia were used. These data were pre-processed, and differentially expressed genes were screened. Besides, gene ontology and pathway enrichment analysis were performed, and then, protein-protein interaction (PPI) network was constructed, and spermatogenesis-related genes in the PPI network were extracted. As a result, 366 up-regulated and 2158 down-regulated genes were identified. Multiple gene ontology terms and pathways including cell-cell signalling and reproduction enriched by differentially expressed genes were obtained. Moreover, four clusters including cluster 1 associated with RNA catabolic process were identified from the PPI network. In addition, genes including cyclin B1, proteasome (prosome, macropain) activator subunit 4, Rac GTPase-activating protein 1 and pituitary tumour-transforming 1 were received. In conclusion, abnormal expression of cyclin B1 and Rac GTPase-activating protein 1, still proteasome (prosome, macropain) activator subunit 4 and pituitary tumour-transforming 1 would impede cell cycle progression during sperm development and maturation, which may contribute to the occurrence and development of teratozoospermia.

Sequencing and Characterization of Divergent Marbling Levels in the Beef Cattle (Longissimus dorsi Muscle) Transcriptome.

Marbling is an important trait regarding the quality of beef. Analysis of beef cattle transcriptome and its expression profile data are essential to extend the genetic information resources and would support further studies on beef cattle. RNA sequencing was performed in beef cattle using the Illumina High-Seq2000 platform. Approximately 251.58 million clean reads were generated from a high marbling (H) group and low marbling (L) group. Approximately 80.12% of the 19,994 bovine genes (protein coding) were detected in all samples, and 749 genes exhibited differential expression between the H and L groups based on fold change (>1.5-fold, p<0.05). Multiple gene ontology terms and biological pathways were found significantly enriched among the differentially expressed genes. The transcriptome data will facilitate future functional studies on marbling formation in beef cattle and may be applied to improve breeding programs for cattle and closely related mammals.

Characterization of transcriptional complexity during adipose tissue development in bovines of different ages and sexes.

BackgroundAdipose tissue has long been recognized to play an extremely important role in development. In bovines, it not only serves a fundamental function but also plays a key role in the quality of beef and, consequently, has drawn much public attention. Age and sex are two key factors that affect the development of adipose tissue, and there has not yet been a global study detailing the effects of these two factors on expressional differences of adipose tissues.ResultsIn this study, total RNA from the back fat of fetal bovines, adult bulls, adult heifers and adult steers were used to construct libraries for Illumina next-generation sequencing. We detected the expression levels of 12,233 genes, with over 3,000 differently expressed genes when comparing fetal and adult patterns and an average of 1000 differently expressed genes when comparing adult patterns. Multiple Gene Ontology terms and pathways were found to be significantly enriched for these differentially expressed genes. Of the 12,233 detected genes, a total of 4,753 genes (38.85%) underwent alternative splicing events, and over 50% were specifically expressed in each library. Over 4,000 novel transcript units were discovered for one library, whereas only approximately 30% were considered to have coding ability, which supplied a large amount of information for the lncRNA study. Additionally, we detected 56,564 (fetal bovine), 65,154 (adult bull), 78,061 (adult heifer) and 86,965 (adult steer) putative single nucleotide polymorphisms located in coding regions of the four pooled libraries.ConclusionHere, we present, for the first time, a complete dataset involving the spatial and temporal transcriptome of bovine adipose tissue using RNA-seq. These data will facilitate the understanding of the effects of age and sex on the development of adipose tissue and supply essential information towards further studies on the genomes of beef cattle and other related mammals.

Phylogenomics and the Dynamic Genome Evolution of the Genus Streptococcus

The genus Streptococcus comprises important pathogens that have a severe impact on human health and are responsible for substantial economic losses to agriculture. Here, we utilize 46 Streptococcus genome sequences (44 species), including eight species sequenced here, to provide the first genomic level insight into the evolutionary history and genetic basis underlying the functional diversity of all major groups of this genus. Gene gain/loss analysis revealed a dynamic pattern of genome evolution characterized by an initial period of gene gain followed by a period of loss, as the major groups within the genus diversified. This was followed by a period of genome expansion associated with the origins of the present extant species. The pattern is concordant with an emerging view that genomes evolve through a dynamic process of expansion and streamlining. A large proportion of the pan-genome has experienced lateral gene transfer (LGT) with causative factors, such as relatedness and shared environment, operating over different evolutionary scales. Multiple gene ontology terms were significantly enriched for each group, and mapping terms onto the phylogeny showed that those corresponding to genes born on branches leading to the major groups represented approximately one-fifth of those enriched. Furthermore, despite the extensive LGT, several biochemical characteristics have been retained since group formation, suggesting genomic cohesiveness through time, and that these characteristics may be fundamental to each group. For example, proteolysis: mitis group; urea metabolism: salivarius group; carbohydrate metabolism: pyogenic group; and transcription regulation: bovis group.