Receptor Interaction Pathway Research Articles

Transcriptional Regulatory Network Analysis for Gastric Cancer Based on mRNA Microarray.

We aimed to screen the differential expressed genes (DEGs) and transcriptional factors (TFs) related to gastric cancer. GSE19826 microarray data downloaded from Gene Expression Omnibus was used to identify the differentially expressed genes (DEGs) and PPI network of DEGs were constructed by the Retrieval of Interacting Genes database. Pathway enrichment analysis of DEGs were performed by Gene Set Enrichment Analysis. Then, the transcriptional regulatory network was constructed based on TRANSFAC database. Finally, regulatory impact factor (RIF) of TF was calculated. We identified 446 DEGs including 209 up- and 237 down-regulated genes. These DEGs were mainly significantly enriched in 5 pathways including ECM receptor interaction (p=0.013899), spliceosome (p=0.025591), bladder cancer (p=0.026316), focal adhesion (p=0.047809) and WNT signaling pathway (p=0.048077). PPI network with 247 nodes and 913 edges were constructed and COL5A2 was the hub node. Transcriptional regulatory network with 6 differently expressed TFs, 58 non-differently expressed TFs, 44 DEGs and 735 non-DEGs was constructed. Finally, top 5 TFs including CRX, TFAP4, NKX2-1, MYB and RARG with higher ZRIF were screened. The identified DEGs such as COL5A2 and TOP2A, and TFs including EGR2, FOXM1, NKX2-1 and TFAP4 might be the critical genes and TFs for gastric cancer.

Transcriptome analysis reveals dynamic changes in coxsackievirus A16 infected HEK 293T cells

BackgroundCoxsackievirus A16 (CVA16) and enterovirus 71 (EV71) are two of the major causes of hand, foot and mouth disease (HFMD) world-wide. Although many studies have focused on infection and pathogenic mechanisms, the transcriptome profile of the host cell upon CVA16 infection is still largely unknown.ResultsIn this study, we compared the mRNA and miRNA expression profiles of human embryonic kidney 293T cells infected and non-infected with CVA16. We highlighted that the transcription of SCARB2, a cellular receptor for both CVA16 and EV71, was up-regulated by nearly 10-fold in infected cells compared to non-infected cells. The up-regulation of SCARB2 transcription induced by CVA16 may increase the possibility of subsequent infection of CVA16/EV71, resulting in the co-infection with two viruses in a single cell. This explanation would partly account for the co-circulation and genetic recombination of a great number of EV71 and CVA16 viruses. Based on correlation analysis of miRNAs and genes, we speculated that the high expression of SCARB2 is modulated by down-regulation of miRNA has-miR-3605-5p. At the same time, we found that differentially expressed miRNA target genes were mainly reflected in the extracellular membrane (ECM)-receptor interaction and circadian rhythm pathways, which may be related to clinical symptoms of patients infected with CVA16, such as aphthous ulcers, cough, myocarditis, somnolence and potentially meningoencephalitis. The miRNAs hsa-miR-149-3p and hsa-miR-5001-5p may result in up-regulation of genes in these morbigenous pathways related to CVA16 and further cause clinical symptoms.ConclusionsThe present study elucidated the changes in 293T cells upon CVA16 infection at transcriptome level, containing highly up-regulated SCARB2 and genes in ECM-receptor interaction and circadian rhythm pathways, and key miRNAs in gene expression regulation. These results provided novel insight into the pathogenesis of HFMD induced by CVA16 infection.

Integrative Analysis of miRNA and mRNA Expression Profiles in Calcium Oxalate Nephrolithiasis Rat Model.

The microRNA (miRNA) expression profiles and their biological functions in calcium oxalate nephrolithiasis remain unclear. In this study, we investigate the miRNA and mRNA expression profiles of kidney tissues in calcium oxalate stone rats. 16 Sprague Dawley rats were divided into control group and stone-forming group. 24-hour urine samples and kidney tissues were collected for biochemical and histological determination after 4 weeks. MiRNA and mRNA microarray were applied to evaluate the miRNA and mRNA expression profiles. To validate the microarray results, the quantitative real-time PCR (qRT-PCR) was performed. A total of 38 miRNAs and 2728 mRNAs were significantly and differentially expressed in kidney tissues of stone-forming group versus control group. Gene Ontology (GO) analysis revealed that most of the target genes were enriched in terms of oxidation reduction, ion transport, inflammatory response, and response to wounding. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of these targets highlights their critical role in cytokine-cytokine receptor interaction, gap junction, and chemokine signaling pathway. Furthermore, the reliability of the microarray-based results was confirmed by using qRT-PCR determination. The miRNA and mRNA expressions in calcium oxalate stone rat kidneys might provide a basis for further research on urolithiasis mechanism.

Detecting discordance enrichment among a series of two-sample genome-wide expression data sets

BackgroundWith the current microarray and RNA-seq technologies, two-sample genome-wide expression data have been widely collected in biological and medical studies. The related differential expression analysis and gene set enrichment analysis have been frequently conducted. Integrative analysis can be conducted when multiple data sets are available. In practice, discordant molecular behaviors among a series of data sets can be of biological and clinical interest.MethodsIn this study, a statistical method is proposed for detecting discordance gene set enrichment. Our method is based on a two-level multivariate normal mixture model. It is statistically efficient with linearly increased parameter space when the number of data sets is increased. The model-based probability of discordance enrichment can be calculated for gene set detection.ResultsWe apply our method to a microarray expression data set collected from forty-five matched tumor/non-tumor pairs of tissues for studying pancreatic cancer. We divided the data set into a series of non-overlapping subsets according to the tumor/non-tumor paired expression ratio of gene PNLIP (pancreatic lipase, recently shown it association with pancreatic cancer). The log-ratio ranges from a negative value (e.g. more expressed in non-tumor tissue) to a positive value (e.g. more expressed in tumor tissue). Our purpose is to understand whether any gene sets are enriched in discordant behaviors among these subsets (when the log-ratio is increased from negative to positive). We focus on KEGG pathways. The detected pathways will be useful for our further understanding of the role of gene PNLIP in pancreatic cancer research. Among the top list of detected pathways, the neuroactive ligand receptor interaction and olfactory transduction pathways are the most significant two. Then, we consider gene TP53 that is well-known for its role as tumor suppressor in cancer research. The log-ratio also ranges from a negative value (e.g. more expressed in non-tumor tissue) to a positive value (e.g. more expressed in tumor tissue). We divided the microarray data set again according to the expression ratio of gene TP53. After the discordance enrichment analysis, we observed overall similar results and the above two pathways are still the most significant detections. More interestingly, only these two pathways have been identified for their association with pancreatic cancer in a pathway analysis of genome-wide association study (GWAS) data.ConclusionsThis study illustrates that some disease-related pathways can be enriched in discordant molecular behaviors when an important disease-related gene changes its expression. Our proposed statistical method is useful in the detection of these pathways. Furthermore, our method can also be applied to genome-wide expression data collected by the recent RNA-seq technology.

Analysis of differentially co-expressed genes based on microarray data of hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is the third leading cause of cancer related death worldwide. Although great progress in diagnosis and management of HCC have been made, the exact molecular mechanisms remain poorly understood. The study aims to identify potential biomarkers for HCC progression, mainly at transcription level. In this study, chip data GSE 29721 was utilized, which contains 10 HCC samples and 10 normal adjacent tissue samples. Differentially expressed genes (DEGs) between two sample types were selected by t-test method. Following, the differentially co-expressed genes (DCGs) and differentially co-expressed Links (DCLs) were identified by DCGL package in Rwith the threshold of q< 0.25. Afterwards, pathway enrichment analysis of the DCGs was carried out by DAVID. Then, DCLs were mapped to TRANSFAC database to reveal associations between relevant transcriptional factors (TFs) and their target genes. Quantitative real-time RT-PCR was performed for TFs or genes of interest. As aresult, atotal of 388 DCGs and 35,771 DCLs were obtained. The predominant pathways enriched by these genes were Cytokine-cytokine receptor interaction, ECM-receptor interaction and TGF-β signaling pathway. Three TF-target interactions, LEF1-NCAM1, EGR1-FN1 and FOS-MT2A were predicted. Compared with control, expressions of the TF genes EGR1, FOS and ETS2 were all up-regulated in the HCC cell line, HepG2; while LEF1 was down-regulated. Except NCAM1, all the target genes were up-regulated in HepG2. Our findings suggest these TFs and genes might play important roles in the pathogenesis of HCC and may be used as therapeutic targets for HCC management.

CpG Oligodeoxynucleotides Induce Differential Cytokine and Chemokine Gene Expression Profiles in Dapulian and Landrace Pigs.

Oligodeoxynucleotides containing unmethylated CpG motifs (CpG ODN) mimic the immunostimulatory activity of microbial DNA by interacting with Toll-like receptor 9 (TLR9) to activate both the innate and adaptive immune responses in different species. However, few studies have been published to compare the effects of CpG ODN on different pig breeds. Therefore, in this study, whole blood gene expression profiles of DPL and Landrace pigs treated with CpG ODN were studied using RNA-seq technology. Five Hundred differentially expressed genes (DEGs) were identified between the two breeds. DPL pigs had significantly higher number of immune-relevant DEGs than the Landrace pigs after CpG ODN treatment. Pathway analysis showed that cytokine-cytokine receptor interaction and chemokine signaling pathway were the major enriched pathways of the immune-relevant DEGs. Further in vitro experiments showed that PBMCs of the DPL pigs had significantly higher levels of TLR9 mRNA than those of the Landrace pigs, both before and after CpG ODN stimulation. Cytokine and chemokine induction in the PBMCs of both breeds were also measured after CpG ODN stimulation. Our data showed that mRNA levels of cytokines (IFNα, IL8, IL12 p40) and chemokines (CXCL9, CXCL13) were significantly higher in the PBMCs of the DPL pigs than those of the Landrace pigs. Taken together, our data provide new information regarding the pig breed difference in response to CpG ODN stimulation and that higher levels of TLR9 mRNA in DPL pigs may be a major contributor for disease resistance.

Bioinformatics analysis of key genes and latent pathway interactions based on the anaplastic thyroid carcinoma gene expression profile.

Anaplastic thyroid carcinoma (ATC) is an aggressive malignant disease in older adults with a high mortality rate. The present study aimed to examine several key genes and pathways, which are associated with ATC. The GSE33630 gene expression profile was downloaded from the Gene Expression Omnibus database, which included 11 ATC and 45 normal thyroid samples. The differentially expressed genes (DEGs) in ATC were identified using the Limma package in R. The Gene Ontology functions and Kyoto Encyclopedia of Genes and Genomes pathways of the selected DEGs were enriched using the Database for Annotation, Visualization and Integrated Discovery. A protein-protein interaction (PPI) network of the DEGs was constructed to select significant modules. Furthermore, a latent pathway interactive network was constructed to select the significant pathways associated with ATC. A total of 665 DEGs in the ATC samples were screened, and four significant modules were selected from the PPI network. The DEGs in the four modules were enriched in several functions and pathways. In addition, 29 significant pathways associated with ATC were selected, and he Toll-like receptor (TLR) signaling pathway, extracellular matrix (ECM)-receptor interaction and cytokine-cytokine interaction pathway were identified as important pathways. FBJ murine osteosarcoma viral oncogene homolog (FOS), chemokine C-X-C motif ligand 10 (CXCL10), collagen type V α1 (COL5A1) and chemokine (C-C motif) ligand 28 (CCL28) were the key DEGs involved in these significant pathways. The data obtained in the present study revealed that the TLR signaling pathway, ECM-receptor interaction and cytokine-cytokine receptor interaction pathway, and the FOS, CXCL10, COL5A1, COL11A1 and CCL28 genes have different roles in the progression of ATC, and these may be used as therapeutic targets for ATC.

Investigation of candidate genes for osteoarthritis based on gene expression profiles

ObjectiveTo explore the mechanism of osteoarthritis (OA) and provide valid biological information for further investigation. MethodsGene expression profile of GSE46750 was downloaded from Gene Expression Omnibus database. The Linear Models for Microarray Data (limma) package (Bioconductor project, http://www.bioconductor.org/packages/release/bioc/html/limma.html) was used to identify differentially expressed genes (DEGs) in inflamed OA samples. Gene Ontology function enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analysis of DEGs were performed based on Database for Annotation, Visualization and Integrated Discovery data, and protein–protein interaction (PPI) network was constructed based on the Search Tool for the Retrieval of Interacting Genes/Proteins database. Regulatory network was screened based on Encyclopedia of DNA Elements. Molecular Complex Detection was used for sub-network screening. Two sub-networks with highest node degree were integrated with transcriptional regulatory network and KEGG functional enrichment analysis was processed for 2 modules. ResultsIn total, 401 up- and 196 down-regulated DEGs were obtained. Up-regulated DEGs were involved in inflammatory response, while down-regulated DEGs were involved in cell cycle. PPI network with 2392 protein interactions was constructed. Moreover, 10 genes including Interleukin 6 (IL6) and Aurora B kinase (AURKB) were found to be outstanding in PPI network. There are 214 up- and 8 down-regulated transcription factor (TF)-target pairs in the TF regulatory network. Module 1 had TFs including SPI1, PRDM1, and FOS, while module 2 contained FOSL1. The nodes in module 1 were enriched in chemokine signaling pathway, while the nodes in module 2 were mainly enriched in cell cycle. ConclusionThe screened DEGs including IL6, AGT, and AURKB might be potential biomarkers for gene therapy for OA by being regulated by TFs such as FOS and SPI1, and participating in the cell cycle and cytokine–cytokine receptor interaction pathway.

Comparative Proteomic Analysis of Bleomycin-induced Pulmonary Fibrosis Based on Isobaric Tag for Quantitation

BackgroundPulmonary fibrosis (PF) is a destructive pulmonary disease and the molecular mechanisms underlying PF are unclear. This study investigated differentially expressed proteins associated with the occurrence and development of PF in rat lung tissue with bleomycin-induced PF. MethodsSixteen Sprague-Dawley rats were randomly divided into 2 groups: the PF model group (n = 8) and the control group (n = 8). After successfully establishing the rat PF model induced by bleomycin, the differentially expressed proteins in the 2 groups were identified through isobaric tag for relative and absolute quantitation coupled with liquid chromatography-mass spectrometry and bioinformatics analysis. ResultsA total of 146 differentially expressed proteins were identified; 88 of which displayed increased abundance and 58 were downregulated in the PF rat model group. Most functional proteins were associated with extracellular matrix, inflammation, damage response, vitamin A synthesis and metabolism. Critical proteins related to PF development and progression was identified, such as type V collagen-3, arachidonic acid 12-lipoxygenase, arachidonic acid 15-lipoxygenase and cytochrome P4501A1. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that these differentially expressed proteins were enriched in extracellular matrix receptor interaction pathway, renin-angiotensin system and metabolic pathway of retinol. ConclusionsThe proteins expressed in bleomycin-induced PF rat model provide important data for further functional analysis of proteins involved in PF.

Somatostatin receptor 5 is a prominent regulator of signaling pathways in cells with coexpression of Cannabinoid receptors 1

Endocannabinoids and somatostatin (SST) play critical roles in several pathophysiological conditions via binding to different receptor subtypes. Cannabinoid receptor 1 (CB1R) and somatostatin receptors (SSTRs) are expressed in several brain regions and share overlapping functions. Whether these two prominent members of G-protein-coupled receptor (GPCR) family interact with each other and constitute a functional receptor complex is not known. In the present study, we investigated the colocalization of CB1R and SSTR5 in rat brain, and studied receptor internalization, interaction and signal transduction pathways in HEK-293 cells cotransfected with human cannabinoid receptor 1 (hCB1R) and hSSTR5. Our results showed that CB1R and SSTR5 colocalized in rat brain cortex, striatum, and hippocampus. CB1R was expressed in SSTR5 immunoprecipitate prepared from the brain tissue lysate, indicating their association in a system where these receptors are endogenously expressed. In cotransfected HEK-293 cells, SSTR5 and CB1R existed in a constitutive heteromeric complex under basal condition, which was disrupted upon agonist treatments. Furthermore, concurrent receptor activation led to preferential formation of SSTR5 homodimer and dissociation of CB1R homodimer. We also discovered that second messenger cyclic adenosine monophosphate and downstream signaling pathways were modulated in a SSTR5-dominant and concentration-dependent manner in the presence of receptor-specific agonist. In conclusion, with predominant role of SSTR5, the functional consequences of crosstalk between SSTR5 and CB1R resulting in the regulation of receptor trafficking and signal transduction pathways open new therapeutic avenue in cancer biology and excitotoxicity.

Dysregulated Wnt signalling and recurrent mutations of the tumour suppressor RNF43 in early gastric carcinogenesis

Several recurrent mutations and epigenetic changes have been identified in advanced gastric cancer, but the genetic alterations associated with early gastric carcinogenesis and malignant transformation remain unclear. We investigated the genomic and transcriptomic landscape of adenomas with low-grade dysplasia (LGD) and high-grade dysplasia (HGD), and intestinal-type early gastric cancer (EGC). The results were validated in an independent cohort that included EGCs directly adjacent to adenoma (EGC-adenomas) that were in the process of malignant transformation, and de novo EGCs that do not seem to have been derived from adenoma. The expression patterns clearly divided into normal, LGD, and EGC, whereas those of HGD overlapped with LGD or EGC. These results suggest that HGD is the critical stage determining malignant transformation. We found that genes related to focal adhesion and extracellular matrix receptor interaction pathways were upregulated as LGD progressed to EGC, whereas canonical Wnt signalling and peroxisome proliferator-activated receptor (PPAR) signalling pathway genes were downregulated in EGC. Genomic alterations such as somatic mutation, gene fusion and copy number variation increased gradually from LGD to EGC. APC mutations were present in 67% of LGDs, 58% of HGDs, and 18% of EGCs. RNF43 mutations were present only in HGD and EGC, and TP53 mutations were present only in EGC. In a validation cohort, RNF43 mutations were present in 35.2% of EGC-adenomas, but in only 8.6% of de novo EGCs. This is the first study to investigate the genomic and transcriptomic landscape of multistep gastric carcinogenesis. We investigated important alterations and their related pathways in each step as tumours progressed from LGD to HGD and eventually to EGC. We suggest that mutations and downregulation of RNF43 may play a critical role in the transition from adenoma to carcinoma. Given these findings and Wnt dependency in tumours with RNF43 mutation, intestinal-type gastric cancer or adenoma with RNF43 mutation might represent a promising indication for Wnt-targeted agents. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Functional analysis of keratinocyte and fibroblast gene expression in skin and keloid scar tissue based on deviation analysis of dynamic capabilities.

The aim of the present study was to select key genes that are associated with fibroblasts and keratinocytes during keloid scar progression and development. The gene expression profile of GSE44270, which includes 32 samples, was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) in case samples compared with control samples were screened using the Limma R package followed by hierarchical clustering analysis. Protein-protein interaction (PPI) networks of the total selected DEGs were constructed using Cytoscape. Moreover, the Gene Ontology biological processes and significant Kyoto Encyclopedia of Genes and Genomes pathways of the total selected DEGs were enriched using the Database for Annotation, Visualization and Integrated Discovery. Significant pathways that may be associated with keloid scar were analyzed using deviation analysis of dynamic capabilities. There were 658 DEGs in fibroblast keloid vs. normal, 112 DEGs in fibroblast non-lesion vs. normal, 439 DEGs in fibroblast keloid vs. non-lesion, 523 DEGs in keratocyte keloid vs. normal, 186 DEGs in keratocyte non-lesion vs. normal, and 963 DEGs in keratocyte keloid vs. non-lesion groups. HOXA9, BMP4, CDKN1A and SMAD2 in fibroblasts, and HOXA7, MCM8, PSMA4 and PSMB2 in keratinocytes were key genes in the PPI networks. Moreover, the amino sugar and nucleotide sugar metabolism pathway, cell cycle, and extracellular matrix (ECM)-receptor interaction pathway were significant pathways. This study suggests that several key genes (BMP4, HOXA9, SMAD2, CDKN1A, HOXA7, PSMA4 and PSMB2) that participate in some significant pathways (cell cycle and ECM-receptor interaction pathways) may be potential therapeutic targets for keloid scars.

Co-expressed differentially expressed genes and long non-coding RNAs involved in the celecoxib treatment of gastric cancer: An RNA sequencing analysis.

The aim of the present study was to investigate the mechanisms of long non-coding RNAs (lncRNAs) in a gastric cancer cell line treated with celecoxib. The human gastric carcinoma cell line NCI-N87 was treated with 15 µM celecoxib for 72 h (celecoxib group) and an equal volume of dimethylsulfoxide (control group), respectively. Libraries were constructed by NEBNext Ultra RNA Library Prep kit for Illumina. Paired-end RNA sequencing reads were aligned to a human hg19 reference genome using TopHat2. Differentially expressed genes (DEGs) and lncRNAs were identified using Cuffdiff. Enrichment analysis was performed using GO-function package and KEGG profile in Bioconductor. A protein-protein interaction network was constructed using STRING database and module analysis was performed using ClusterONE plugin of Cytoscape. ATP5G1, ATP5G3, COX8A, CYC1, NDUFS3, UQCRC1, UQCRC2 and UQCRFS1 were enriched in the oxidative phosphorylation pathway. CXCL1, CXCL3, CXCL5 and CXCL8 were enriched in the chemokine signaling and cytokine-cytokine receptor interaction pathways. ITGA3, ITGA6, ITGB4, ITGB5, ITGB6 and ITGB8 were enriched in the integrin-mediated signaling pathway. DEGs co-expressed with lnc-SCD-1:13, lnc-LRR1-1:2, lnc-PTMS-1:3, lnc-S100P-3:1, lnc-AP000974.1-1:1 and lnc-RAB3IL1-2:1 were enriched in the pathways associated with cancer, such as the basal cell carcinoma pathway in cancer. In conclusion, these DEGs and differentially expressed lncRNAs may be important in the celecoxib treatment of gastric cancer.

Identification of distinct genes associated with seawater aspiration‑induced acute lung injury by gene expression profile analysis.

Seawater aspiration-induced acute lung injury (ALI) is a syndrome associated with a high mortality rate, which is characterized by severe hypoxemia, pulmonary edema and inflammation. The present study is the first, to the best of our knowledge, to analyze gene expression profiles from a rat model of seawater aspiration-induced ALI. Adult male Sprague-Dawley rats were instilled with seawater (4 ml/kg) in the seawater aspiration-induced ALI group (S group) or with distilled water (4 ml/kg) in the distilled water negative control group (D group). In the blank control group (C group) the rats' tracheae were exposed without instillation. Subsequently, lung samples were examined by histopathology; total protein concentration was detected in bronchoalveolar lavage fluid (BALF); lung wet/dry weight ratios were determined; and transcript expression was detected by gene sequencing analysis. The results demonstrated that histopathological alterations, pulmonary edema and total protein concentrations in BALF were increased in the S group compared with in the D group. Analysis of differential gene expression identified up and downregulated genes in the S group compared with in the D and C groups. A gene ontology analysis of the differential gene expression revealed enrichment of genes in the functional pathways associated with neutrophil chemotaxis, immune and defense responses, and cytokine activity. Kyoto Encyclopedia of Genes and Genomes analysis revealed that the cytokine-cytokine receptor interaction pathway was one of the most important pathways involved in seawater aspiration-induced ALI. In conclusion, activation of the cytokine-cytokine receptor interaction pathway may have an essential role in the progression of seawater aspiration-induced ALI, and the downregulation of tumor necrosis factor superfamily member 10 may enhance inflammation. Furthermore, IL-6 may be considered a biomarker in seawater aspiration-induced ALI.

The effects of Qi Teng Xiao Zhuo granules, traditional Chinese medicine, on the expression of genes in chronic glomerulonephritis rats

Background and aimChronic glomerulonephritis (CGN) is a primary glomerular disease that is related to immune-mediated inflammatory diseases. Qi Teng Xiao Zhuo granules have been proposed as a prescription of traditional Chinese medicine for treatment of CGN, but the comprehensive molecular mechanism underlying this therapeutic effect is not clear to date. The aim of this study was to evaluate and analyze the possible roles and molecular mechanisms of Qi Teng Xiao Zhuo granule-mediated treatment of CGN induced by adriamycin in rats. MethodsFor gene expression analysis, four samples of glomerular tissue from rats in the Qi Teng Xiao Zhuo granule group and four samples each from the adriamycin treated and control groups were hybridized with Agilent Rat 4×44K whole genome microarrays. KEGG and Gene Ontology (GO) analyses and LIMMA, String and Cytoscape software were used to analyze the functional microarray data and screen differentially expressed genes. Hub genes were identified using Pathway Studio software. Real-time PCR was performed to verify the selected genes. ResultsMicroarray gene expression analysis showed that Pnoc, Cacfd1, Fos, Igll1, Lcn2, and Syk were among the most downregulated genes in the Qi Teng Xiao Zhuo granule group compared with the adriamycin treated group, whereas Cyp2c7, Hsd3b6, Acsm5, and Ugt2b15 were significantly upregulated. Functional analysis demonstrated that metabolism of xenobiotics by cytochrome P450, the B cell receptor signaling pathway, and cytokine-cytokine receptor interaction pathways were significantly downregulated in the Qi Teng Xiao Zhuo granule group and that GO terms related to positive regulation of immune response, immune response-activating signal transduction, cell differentiation, cell cycle, proliferation, and adhesion were significantly affected. Fos and Syk were considered to be potential hub genes. ConclusionsIn the adriamycin-induced CGN rat model, comprehensive molecular mechanisms were involved with complex gene expression alterations containing many altered pathways and GO terms. However, how Qi Teng Xiao Zhuo granules regulate these events warrants further investigation.

Deciphering Transcriptome and Complex Alternative Splicing Transcripts in Mammary Gland Tissues from Cows Naturally Infected with Staphylococcus aureus Mastitis.

Alternative splicing (AS) contributes to the complexity of the mammalian proteome and plays an important role in diseases, including infectious diseases. The differential AS patterns of these transcript sequences between the healthy (HS3A) and mastitic (HS8A) cows naturally infected by Staphylococcus aureus were compared to understand the molecular mechanisms underlying mastitis resistance and susceptibility. In this study, using the Illumina paired-end RNA sequencing method, 1352 differentially expressed genes (DEGs) with higher than twofold changes were found in the HS3A and HS8A mammary gland tissues. Gene ontology and KEGG pathway analyses revealed that the cytokine–cytokine receptor interaction pathway is the most significantly enriched pathway. Approximately 16k annotated unigenes were respectively identified in two libraries, based on the bovine Bos taurus UMD3.1 sequence assembly and search. A total of 52.62% and 51.24% annotated unigenes were alternatively spliced in term of exon skipping, intron retention, alternative 5′ splicing and alternative 3ʹ splicing. Additionally, 1,317 AS unigenes were HS3A-specific, whereas 1,093 AS unigenes were HS8A-specific. Some immune-related genes, such as ITGB6, MYD88, ADA, ACKR1, and TNFRSF1B, and their potential relationships with mastitis were highlighted. From Chromosome 2, 4, 6, 7, 10, 13, 14, 17, and 20, 3.66% (HS3A) and 5.4% (HS8A) novel transcripts, which harbor known quantitative trait locus associated with clinical mastitis, were identified. Many DEGs in the healthy and mastitic mammary glands are involved in immune, defense, and inflammation responses. These DEGs, which exhibit diverse and specific splicing patterns and events, can endow dairy cattle with the potential complex genetic resistance against mastitis.

High expression of WISP1 in colon cancer is associated with apoptosis, invasion and poor prognosis.

Colon cancer (CC) likes many epithelial-derived cancers, resulting from a complex tumorigenic process. However, the exactly mechanisms of development and progression of CC are still unknown. In this study, integrated analysis in the GSE33113 and Fudan University Shanghai Cancer Center Hospital datasets revealed that WISP1 expression was significantly increased in CC cases, positivity correlated with the advanced pathologic stage and a poor prognosis was more likely in CC patients with higher levels of WISP1. Downregulation of WISP1 inhibited cell proliferation and invasion through increasing apoptosis and blocking cell cycle at G1 phase in CC LOVO and RKO cells. Besides, Gene set enrichment analysis (GSEA) revealed that relative genes involved in the Cell adhesion molecules and Cytokine-cytokine receptor interaction pathways were enriched in WISP1-higher expression patients. Western blot analysis showed that Cell adhesion molecules pathway associated genes (ICAM- 1, VCAM-1, SDC2 and CDH2) and Cytokine-cytokine receptor interaction pathway associated genes (VEGFC, CCL18, CXCR4 and TGFBR1) were also modulated by WISP1 downregulation. Then, we found that the protein β-catenin was identified as a binding partner of WISP1 and mediated the functions of WISP1 through promoting cell proliferation and invasion in LOVO and RKO cells. Further in vivo tumor formation study in nude mice indicated that inhibition of WISP1 delayed the progress of tumor formation and inhibited PCNA expression. These results indicate that WISP1 could act as an oncogene and may serve as a promising therapeutic strategy for colon cancer.

Activation of the Met kinase confers acquired drug resistance in FGFR-targeted lung cancer therapy.

Aberrant fibroblast growth factor receptor (FGFR) activation/expression is a common feature in lung cancer (LC). In this study, we evaluated the antitumor activity of and the mechanisms underlying acquired resistance to two potent selective FGFR inhibitors, AZD4547 and BAY116387, in LC cell lines. The antitumor activity of AZD4547 and BAY1163877 was screened in 24 LC cell lines, including 5 with FGFR1 amplification. Two cell lines containing FGFR1 amplifications, H1581 and DMS114, were sensitive to FGFR inhibitors (IC50<250 nm). Clones of FGFR1-amplified H1581 cells resistant to AZD4547 or BAY116387 (H1581AR and H1581BR cells, respectively) were established. Receptor tyrosine kinase (RTK) array and immunoblotting analyses showed strong overexpression and activation of Met in H1581AR/BR cells, compared with that in the parental cells. Gene set enrichment analysis against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database showed that cytokine–cytokine receptor interaction pathways were significantly enriched in H1581AR/BR cells, with Met contributing significantly to the core enrichment. Genomic DNA quantitative PCR and fluorescent in situ hybridization analyses showed MET amplification in H1581AR, but not in H1581BR, cells. Met amplification drives acquired resistance to AZD4547 in H1581AR cells by activating ErbB3. Combination treatment with FGFR inhibitors and an anaplastic lymphoma kinase (ALK)/Met inhibitor, crizotinib, or Met-specific short interfering RNA (siRNA) synergistically inhibited cell proliferation in both H1581AR and H1581BR cells. Conversely, ectopic expression of Met in H1581 cells conferred resistance to AZD4547 and BAY1163877. Acquired resistance to FGFR inhibitors not only altered cellular morphology, but also promoted migration and invasion of resistant clones, in part by inducing epithelial-to-mesenchymal transition. Taken together, our data suggest that Met activation is sufficient to bypass dependency on FGFR signaling. Concurrent inhibition of the Met and FGFR pathways may have synergistic clinical benefits when targeting FGFR-dependent LC.

Screening of gene signatures for rheumatoid arthritis and osteoarthritis based on bioinformatics analysis.

The current study aimed to identify gene signatures during rheumatoid arthritis (RA) and osteoarthritis (OA), and used these to elucidate the underlying modular mechanisms. Using the Gene Expression Omnibus database, the present study obtained the GSE7669 mRNA expression microarray data from RA and OA synovial fibroblasts (n=6 each). The differentially expressed genes (DEGs) in RA synovial samples compared with OA samples were identified using the Linear Models for Microarray Analysis package. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed using the Database for Annotation Visualization and Integrated Discovery. A protein-protein interaction network was constructed and the modules were further analyzed using the Molecular Complex Detection plugin of Cytoscape. A total of 181 DEGs were identified by comparing RA and OA synovial samples (96 up- and 85 downregulated genes). The significant DEGs in module 1, including collagen, type I, α 1 (COL1A1), COL3A1, COL4A1 and COL11A1, were predominantly enriched in the extracellular matrix (ECM)-receptor interaction and focal adhesion pathways. Additionally, significant DEGs in module 2, including radical S-adenosyl methionine domain containing 2 (RSAD2), 2′-5′-oligoadenylate synthetase 2 (OAS2), myxovirus (influenza virus) resistance 1 (MX1) and ISG15 ubiquitin-like modifier (ISG15), were predominantly associated with immune function pathways. In conclusion, the present study indicated that RSAD2, OAS2, MX1 and ISG15 may be notable gene signatures in RA development via regulation of the immune response. COL3A1, COL4A1, COL1A1 and COL11A1 may be important gene signatures in OA development via involvement in the pathways of ECM-receptor interactions and focal adhesions.

Analysis of the cell-free amniotic fluid transcriptome expressed during the euploid mid-trimester of pregnancy

ObjectiveThe amniotic fluid (AF) contains cell-free RNAs (cfRNAs), which are considered to reflect the fetal status in utero. However, there are limited numbers of data to examine the AF cell-free transcriptome because amniocentesis is an invasive procedure. In this study, the AF transcriptome expressed during the euploid mid-trimester of pregnancy was characterized. Study designFourteen AF samples were collected. RNA was extracted from AF supernatant, hybridized to Affymetrix GeneChip Human arrays, and the transcriptome was analyzed by using the DAVID toolkit. ResultWe detected 1069 genes in the 14 AF samples. The GNF atlas mapping showed that genes present in the AF were annotated with endocrine organs and blood components, including the pancreas, adrenal gland, thyroid, ovary and monocytes. The proteins encoded by the transcriptome were localized to several organs, which are directly in contact with the AF, including the placenta, lung, skin, epithelium, and kidney. During the early fetal period, there is a bi-directional diffusion between the fetus and AF. Therefore, the AF composition is similar to that of the fetal plasma. In addition, fetal urine, swallowing, pulmonary secretion, and diffusion across the placenta contribute to produce amniotic fluid by directly excreting fluid. The KEGG pathway analysis with placenta specific genes revealed that focal adhesion and extracellular matrix receptor interaction pathways were enriched. These pathways are important for the placental development. ConclusioncfRNA in the amniotic fluid originates from placenta and fetal organs directly contacting the amniotic fluid as well as from diffusion of the fetal plasma across the placenta. AF transcriptome may reflect not only fetal development, but also placental development.