Transcriptome and Weighted Gene Co-Expression Network Analysis for Feather Follicle Density in a Chinese Indigenous Breed.

Feather follicle density plays an important role in appealing to consumers' first impressions when making purchasing decisions. However, the molecular network that contributes to this trait remains largely unknown. The aim of this study was to perform transcriptome and weighted gene co-expression network analyses to determine the candidate genes relating to feather follicle density in Wannan male chickens. In total, five hundred one-day-old Wannan male chickens were kept in a conventional cage system. Feather follicle density was recorded for each bird at 12 weeks of age. At 12 weeks, fifteen skin tissue samples were selected for weighted gene co-expression network analysis, of which six skin tissue samples (three birds in the H group and three birds in the L group) were selected for transcriptome analysis. The results showed that, in total, 95 DEGs were identified, and 56 genes were upregulated and 39 genes were downregulated in the high-feather-follicle-density group when compared with the low-feather-follicle-density group. Thirteen co-expression gene modules were identified. The red module was highly significantly negatively correlated with feather follicle density (p < 0.01), with a significant negative correlation coefficient of -0.72. In total, 103 hub genes from the red module were screened. Upon comparing the 103 hub genes with differentially expressed genes (DEGs), it was observed that 13 genes were common to both sets, including MELK, GTSE1, CDK1, HMMR, and CENPE. From the red module, FOXM1, GTSE1, MELK, CDK1, ECT2, and NEK2 were selected as the most important genes. These genes were enriched in the DNA binding pathway, the heterocyclic compound binding pathway, the cell cycle pathway, and the oocyte meiosis pathway. This study suggests that FOXM1, GTSE1, MELK, CDK1, ECT2, and NEK2 may be involved in regulating the development of feather follicle density in Wannan male chickens. The results of this study reveal the genetic structure and molecular regulatory network of feather follicle density in Wannan male chickens, and provide a basis for further elucidating the genetic regulatory mechanism and identifying molecular markers with breeding value.

Prenatal PM2.5 exposure impairs spatial learning and memory in male mice offspring: from transcriptional regulation to neuronal morphogenesis.

BACKGROUND: As one of the environmental risk factors for human health, atmospheric fine particulate matter (PM2.5) contributes to cognitive deterioration in addition to respiratory and cardiovascular injuries. Recently, increasing evidence implicates that PM2.5 inhalation can affect neurological functions in offspring, but the sex-specific outcomes and the underlying biological processes are largely unknown. OBJECTIVES: To observe the influence of prenatal PM2.5 exposure on cognitive performance in offspring, to elucidate the neuronal morphological alterations and possible transcriptional regulation based on mRNA-sequencing (mRNA-Seq) data after birth, and to determine the key components of PM2.5 contributing to the adverse effects. METHODS: Pregnant C57BL/6J mice were exposed to sterile saline or PM2.5 suspension. Morris water maze test was used to assess the cognitive function in weanling offspring. Microscopic observation was applied to detect neuronal morphogenesis in vivo and in vitro. The cortex tissues from male offspring were collected on postnatal days (PNDs) 1, 7, and 21 for mRNA-Seq analysis. The organic and inorganic components of PM2.5 were separated to assess their contributions using primary cultured neurons. RESULTS: Prenatal PM2.5 exposure impaired spatial learning and memory in weanling male mice, but not female mice. The sex-specific outcomes were associated with mRNA expression profiles of the cortex during postnatal critical windows, and the annotations in Gene Ontology (GO) of differentially expressed genes (DEGs) revealed that the exposure persistently disrupted the expression of genes involved in neuronal features in male offspring. Consistently, axonal growth impairment and dendritic complexity reduction were observed. Importantly, Homeobox A5 (Hoxa5), a critical transcription factor regulating all of the neuronal morphogenesis-associated hub genes on PNDs 1, 7, and 21, significantly decreased in the cortex of male offspring following PM2.5 exposure. In addition, both inorganic and organic components were harmful to axonal and dendritic growth, with organic components exhibiting stronger inhibition than inorganic ones. CONCLUSION: Prenatal PM2.5 exposure affected spatial learning and memory in male mice by disrupting Hoxa5-mediated neuronal morphogenesis, and the organic components, including polycyclic aromatic hydrocarbons (PAHs), posed more adverse effects than the inorganic components.

Based on mRNA Sequencing Techniques to Explore the Molecular Mechanism of Buzhong Yiqi Decoction for Autoimmune Thyroiditis.

OBJECTIVE: Autoimmune diseases (AD) account for a high percentage of the population. One of the most prevalent is autoimmune thyroiditis (AIT). However, the therapeutic effects of Buzhong Yiqi (BZYQ) decoction on AIT have not been studied yet. The majority of the present study was conducted on NOD.H-2h4 mice in an attempt to ascertain the therapeutic effects of BZYQ decoction on AIT. METHODS: The 0.05% sodium iodide water (NaI)-induced AIT mice model was established. A total of nine NOD.H-2h4 mice were randomly divided into three groups: the normal group provided with regular water, the model group drinking freely 0.05% NaI, and the treatment group treated with BZYQ decoction (9.56 g/kg) after NaI supplementation (NaI + BZYQ). BZYQ decoction was administered orally once daily for eight weeks. The thyroid histopathology test was used to measure the severity of lymphocytic infiltration. An enzyme-linked immunosorbent assay (ELISA) was used to determine the levels of anti-thyroglobulin antibody (TgAb), interleukin (IL)-1ß, IL-6, and IL-17. The Illumina HiSeq X sequencing platform was utilized to analyze the thyroid tissue by mRNA expression profiles. Bioinformatics analysis was used to investigate the biological function of the differentially expressed mRNAs. In addition, the expression of Carbonyl Reductase 1 (CBR1), 6-Pyruvoyltetrahydropterin Synthase (PTS), Major Histocompatibility Complex, Class II (H2-EB1), Interleukin 23 Subunit Alpha (IL-23A), Interleukin 6 Receptor (IL-6RA), and Janus Kinase 1 (JAK1) was measured by quantitative real-time PCR (qRT-PCR). RESULTS: The treatment group exhibited significantly lower rates of thyroiditis and lymphocyte infiltration compared to the model group. Serum levels of TgAb, IL-1ß, IL-6, and IL-17 were significantly higher in the model group, but they fell dramatically after BZYQ decoction administration. According to our results, 495 genes showed differential expression in the model group compared to the control group. Six hundred twenty-five genes were significantly deregulated in the treatment group compared to the model group. Bioinformatic analysis showed that most mRNAs were associated with immune-inflammatory responses and were involved in multiple signaling pathways, including folate biosynthesis and the Th17 cell differentiation pathway. CBR1, PTS, H2-EB1, IL23A, IL-6RA and JAK1 mRNA participated in folate biosynthesis and the Th17 cell differentiation pathway. The qRT-PCR analysis confirmed that the above mRNAs were regulated in the model group compared to the treatment group Conclusion: The results of this investigation have revealed novel insights into the molecular mechanism of action of BZYQ decoction against AIT. The mechanism may be partially attributed to the regulation of mRNA expression and pathways.

Low expression of ESR1 correlates with ascending aortic dilation and acute type A aortic dissection.

BACKGROUND: The common genetic signatures of ascending aortic dilation and acute type A aortic dissection (AAAD) remained unclear. This study aims to investigate the role of ESR1 in ascending aortic dilation and AAAD via a multi-omics approach combining transcriptome-wide association study (TWAS) and mRNA expression profiles. METHODS: The TWAS analysis was performed by integrating expression quantitative trait loci (eQTL) data of aorta and the GWAS dataset of ascending aortic diameter using the FUSION software. Joint/Conditional Analysis was used to screen conditionally independent genes from TWAS significant regions. mRNA expression profiles were used to confirm the differential expression of ESR1 in ascending aortic dilation and AAAD. An independent mRNA expression profile dataset was used to validate the diagnostic efficiency of ESR1 expression on thoracic aortic aneurysm (TAA). Gene set enrichment analysis (GSEA) was utilized to explore the potential molecular function of ESR1 in mouse aorta. Immune infiltration analysis was performed to evaluate the association between ESR1 and immune infiltration in AAAD. RESULTS: ESR1 was among the top 10 most significant genes identified by TWAS of ascending aortic diameter (Z score = -7.3, PTWAS = 3.30 × 10-13). mRNA expression profiles confirmed the low expression of ESR1 both in ascending aortic dilation (PmRNA = 0.0361) and AAAD (PmRNA = 0.0142). The diagnostic efficiency of peripheral blood ESR1 expression on TAA was further validated by an independent mRNA expression profile dataset. GSEA showed that GO terms and KEGG pathways mainly involved in metabolism and oxidoreductase activity were enriched in aortic tissue of ERα knockout mice. The immune infiltration ration of naive CD8 + T cells was significantly lower in AAAD compared with normal aorta, and was positively correlated with ESR1 expression. CONCLUSION: Low expression of ESR1 is associated with ascending aortic dilation and AAAD. Peripheral blood ESR1 expression may be a novel diagnostic biomarker of TAA. ESR1 expression is positively correlated with immune infiltration ration of naive CD8 + T cells in AAAD.

A Cross-Tissue Transcriptome-Wide Association Study Identifies Novel Susceptibility Genes for Juvenile Idiopathic Arthritis in Asia and Europe.

Background: Juvenile idiopathic arthritis (JIA) is the most common rheumatic disease in children, and its pathogenesis is still unclear. Genome-wide association studies (GWASs) of JIA have identified hundreds of risk factors, but few of them implicated specific biological mechanisms. Methods: A cross-tissue transcriptome-wide association study (TWAS) was performed with the functional summary-based imputation software (FUSION) tool based on GWAS summary datasets (898 JIA patients and 346,102 controls from BioBank Japan (BBJ)/FinnGen). The gene expression reference weights of skeletal muscle and the whole blood were obtained from the Genotype-Tissue Expression (GTExv8) project. JIA-related genes identified by TWAS findings genes were further compared with the differentially expressed genes (DEGs) identified by the mRNA expression profile of JIA from the Gene Expression Omnibus (GEO) database (accession number: GSE1402). Last, candidate genes were analyzed using functional enrichment and annotation analysis by Metascape to examine JIA-related gene sets. Results: The TWAS identified 535 significant genes with P < 0.05 and contains 350 for Asian and 195 for European (including 10 genes both expressed in Asian and European), such as CDC16 (P = 1.72E-03) and PSMD5-AS1 (P = 3.65E-02). Eight overlapping genes were identified based on TWAS results and DEGs of JIA patients, such as SIRPB1 (PTWAS = 4.21E-03, PDEG = 1.50E-04) and FRAT2 (PTWAS = 2.82E-02, PDEG = 1.43E-02). Pathway enrichment analysis of TWAS identified 183 pathways such as cytokine signaling in the immune system and cell adhesion molecules. By integrating the results of DEGs pathway and process enrichment analyses, 19 terms were identified such as positive regulation of T-cell activation. Conclusion: By conducting two populations TWAS, we identified a group of JIA-associated genes and pathways, which may provide novel clues to uncover the pathogenesis of JIA.

Bioinformatics Screening of Potential Biomarkers from mRNA Expression Profiles to Discover Drug Targets and Agents for Cervical Cancer.

Bioinformatics analysis has been playing a vital role in identifying potential genomic biomarkers more accurately from an enormous number of candidates by reducing time and cost compared to the wet-lab-based experimental procedures for disease diagnosis, prognosis, and therapies. Cervical cancer (CC) is one of the most malignant diseases seen in women worldwide. This study aimed at identifying potential key genes (KGs), highlighting their functions, signaling pathways, and candidate drugs for CC diagnosis and targeting therapies. Four publicly available microarray datasets of CC were analyzed for identifying differentially expressed genes (DEGs) by the LIMMA approach through GEO2R online tool. We identified 116 common DEGs (cDEGs) that were utilized to identify seven KGs (AURKA, BRCA1, CCNB1, CDK1, MCM2, NCAPG2, and TOP2A) by the protein-protein interaction (PPI) network analysis. The GO functional and KEGG pathway enrichment analyses of KGs revealed some important functions and signaling pathways that were significantly associated with CC infections. The interaction network analysis identified four TFs proteins and two miRNAs as the key transcriptional and post-transcriptional regulators of KGs. Considering seven KGs-based proteins, four key TFs proteins, and already published top-ranked seven KGs-based proteins (where five KGs were common with our proposed seven KGs) as drug target receptors, we performed their docking analysis with the 80 meta-drug agents that were already published by different reputed journals as CC drugs. We found Paclitaxel, Vinorelbine, Vincristine, Docetaxel, Everolimus, Temsirolimus, and Cabazitaxel as the top-ranked seven candidate drugs. Finally, we investigated the binding stability of the top-ranked three drugs (Paclitaxel, Vincristine, Vinorelbine) by using 100 ns MD-based MM-PBSA simulations with the three top-ranked proposed receptors (AURKA, CDK1, TOP2A) and observed their stable performance. Therefore, the proposed drugs might play a vital role in the treatment against CC.

Integrative Analysis of Transcriptome-Wide Association Study and mRNA Expression Profiles Identified Candidate Genes and Pathways Associated With Acute Myocardial Infarction.

BACKGROUND: Acute myocardial infarction (AMI), characterized by an event of myocardial necrosis, is a common cardiac emergency worldwide. However, the genetic mechanisms of AMI remain largely elusive. METHODS: A genome-wide association study dataset of AMI was obtained from the CARDIoGRAMplusC4D project. A transcriptome-wide association study (TWAS) was conducted using the FUSION tool with gene expression references of the left ventricle and whole blood. Significant genes detected by TWAS were subjected to Gene Ontology (GO) enrichment analysis. Then the TWAS results of AMI were integrated with mRNA expression profiling to identify common genes and biological processes. Finally, the identified common genes were validated by RT-qPCR analysis. RESULTS: TWAS identified 1,050 genes for the left ventricle and 1,079 genes for whole blood. Upon comparison with the mRNA expression profile, 4 common genes were detected, including HP (PTWAS = 1.22 × 10-3, PGEO = 4.98 × 10-2); CAMP (PTWAS = 2.48 × 10-2, PGEO = 2.36 × 10-5); TNFAIP6 (PTWAS = 1.90 × 10-2, PGEO = 3.46 × 10-2); and ARG1 (PTWAS = 8.35 × 10-3, PGEO = 4.93 × 10-2). Functional enrichment analysis of the genes identified by TWAS detected multiple AMI-associated biological processes, including autophagy of mitochondrion (GO: 0000422) and mitochondrion disassembly (GO: 0061726). CONCLUSION: This integrative study of TWAS and mRNA expression profiling identified multiple candidate genes and biological processes for AMI. Our results may provide a fundamental clue for understanding the genetic mechanisms of AMI.

Differentiation potential and mRNA profiles of human dedifferentiated adipose cells and adipose­derived stem cells from young donors.

Dedifferentiated adipose cells (DAs) and adipose­derived stem cells (ADSCs) are two of the primary types of stem cells derived from adipose tissue, which have been reported to possess similar characteristics, but also exhibit unique phenotypic and functional advantages. However, several reports have described inconsistent results regarding their differences in multilineage differentiation function. Moreover, to the best of our knowledge, there are no studies assessing their myogenic ability, or the differences in the transcriptome between the two cell types derived from lipoaspirates via tumescent liposuction from the same donors. The aim of the present study was to compare the properties and expression profiles of these cell types. Subcutaneous adipose tissue of three female patients (aged 23­30 years) with a physiological BMI (19.1­23.9 kg/m2) were obtained during tumescent liposuction of the abdomen or the thigh. The stromal vascular fraction and mature adipocytes were obtained via collagenase digestion, and ADSCs and DAs were cultured successively. To determine the differences between DAs and ADSCs after 6­7 passages, cell proliferation assays, phenotypic assessment, differentiation assays and high­throughput RNA sequencing (seq) were used. Similar cell morphologies, proliferation dynamics, surface markers and transcriptome expression profiles were observed between the DAs and ADSCs. Whilst there were notable individual differences in the osteogenic, lipogenic, chondrogenic and myogenic abilities of the DAs and ADSCs, it was difficult to determine their differentiation potential based only on the cell source. Interestingly, the myogenic ability was relatively stronger in cells with relatively weaker lipogenic ability. Only 186 differentially expressed genes between the two groups were identified using RNAseq. Several of these genes were involved in biological functions such as transcription regulation, protein translation regulation, cytokine interactions and energy metabolism regulation. The results of the present study suggested a similar functional potential of DAs and ADSCs from young donors undergoing tumescent liposuction operation in regeneration areas and the balance of the differentiative ability of the same cell populations. These data may provide a foundation for further clinical administration of stem cells derived from adipose tissues in therapy.

Integrative analysis of miRNA and mRNA expression profiles reveals a novel mRNA/miRNA signature to improve risk classification for patients with gastric cancer.

Gastric cancer (GC) is one of the most common types of malignant cancer and is associated with poor prognosis. Although the prognosis of patients with GC is associated with grade, stage and lymph node metastases, these traditional clinical features are inadequate to predict the outcome of GC. Therefore, there has been an increased focus on identifying novel molecular biomarkers for early diagnosis and prognosis, in order to improve outcomes in GC. In the present study, an integrative analysis of microRNA (miRNA) expression profiles, mRNA expression profiles and clinical characteristics was performed in a large cohort of patients with GC in order to identify an integrative prognostic model for improving postoperative risk classification. An integrative mRNA/miRNA signature (IMMIS), comprised of three miRNAs and one mRNA, was identified from a large number of differentially expressed miRNAs and mRNAs using univariate and multivariate Cox regression analysis. The prognostic value of the IMMIS was validated in the discovery cohort, testing cohort and The Cancer Genome Atlas (TCGA) cohort. The present results suggested that the identified signature had a reliable predictive performance and could classify the patients into high- and low-risk groups with significantly different overall survival times. In the discovery cohort, the hazard ratio (HR) was 2.805 with a 95% CI=1.722-4.567 (P<0.001). The median overall survival time as 1.49 vs. 3.85 years. In the testing cohort, the HR was 1.625 with a 95% CI=1.004-2.638 (P=0.039) and the median overall survival time was 2.17 vs. 4.62 years. In the TCGA cohort, the HR was 2.139 with a 95% CI=1.519-3.012 (P<0.001) and the median overall survival time was 1.53 vs. 4.62 years. The IMMIS constituted a reliable independent prognostic factor compared with clinical covariates, including age, sex, grade and stage, as indicated by multivariate and stratified analyses. Furthermore, comparative analysis revealed that the predictive value of the IMMIS was superior to the mRNA-based signature alone. The present results suggested the potential value of the IMMIS as a promising novel biomarker for improving the clinical management of patients with GC.

Comparison of mRNA expression profiles of drug-metabolizing enzymes and transporters in fresh and cryopreserved cynomolgus monkey hepatocytes.

In this study, freshly isolated and cryopreserved cynomolgus monkey hepatocytes were seeded on Cell-able® plates with feeder cells to form spheroids and were cultured for 28 days. As a control, hepatocytes were also cultured with or without feeder cells on collagen-coated plates. We verified the mRNA expression levels of drug-metabolizing enzyme-related genes and the leakage of enzymes (AST, ALT, LDH, and γ-GTP) as indicators of cell survival. As a result, the patterns of target mRNA expression in fresh and cryopreserved hepatocytes were very similar during the culture period between culture methods. mRNA expression levels were highly maintained at day 28 using the 3D spheroid and co-culture methods, demonstrating that these methods are useful for maintenance of liver function. Leakage of AST and ALT was higher at day 3 but decreased at day 14. LDH was not detected, suggesting that the cell viability was also maintained during the culture period. Furthermore, the functional differences between fresh and cryopreserved hepatocytes were not clearly detected. The co-culture method was useful for long-term culture not requiring 3D structure, and the 3D spheroid culture method was effective as well. With these techniques, cynomolgus monkey hepatocytes are expected to exhibit smaller individual differences and high reproducibility.

Identification of global mRNA expression profiles and comprehensive bioinformatic analyses of abnormally expressed genes in cholestatic liver disease.

BACKGROUND: Cholestatic liver disease (CLD) is a highly heterogeneous hepatobiliary disease with various causes. The purpose of this research was to explore the gene expression changes throughout the course of CLD revealing potential causative molecular mechanisms and therapeutic targets. METHODS: We established two animal models of cholestasis: 3,5-diethoxycarbonyl-1,4-dihydrocollidine feeding for 2, 4 and 6 weeks and bile duct ligation for 14 days. Using these two models, we identified differentially expressed genes (DEGs) by RNA-Seq analysis and used the newly-found knowledge of DEGs in comprehensive bioinformatic analyses to investigate key molecular events. Sequencing results were confirmed by experimental verification. RESULTS: Our study detected overlapping DEGs in the two models, of these 568 genes were upregulated and 117 genes were downregulated. Gene Ontology analysis demonstrated that the upregulated genes were associated with the biological processes of cell adhesion, cell migration and cell motility, while the metabolic processes of various substances were enriched for the downregulated genes. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the upregulated pathways were mainly distributed in focal adhesion, ECM-receptor interaction and amoebiasis, while downregulated pathways focused on peroxisome proliferator-activated receptor signaling pathway, metabolic pathways and primary bile acid biosynthesis. These findings were further confirmed by protein-protein interaction network modeling. Hub genes Src, Pdgfb, Col15a1, Mmp9, Egfr were selected using centralities analyses and verified by qRT-PCR. CONCLUSION: We profiled a global mRNA landscape in CLD to promote a complete understanding of transcriptomic events of this disease, offering candidate biomarkers and therapeutic targets for the clinic.

Comparative analysis of mRNA expression profiles in Type 1 and Type 2 diabetes mellitus.

Aim: We aimed to understand the individual and shared features of Type 1 diabetes (T1D) and Type 2 diabetes (T2D) by analyzing the gene expression profile. Materials & methods: An integrated analysis was performed with microarray datasets for T1D and T2D. Compared with normal control, shared and specific differentially expressed genes (DEGs) in T1D and T2D were obtained. Functional annotation, further validation and receiver operating characteristic curve analysis were performed. Results: Five and three datasets for T1D and T2D were downloaded, respectively. In total, 141 (85 T1D vs 56 normal controls) and 70 (29 T2D vs 41 normal controls) peripheral blood samples were included in T1D and T2D group, respectively. Compared with normal controls, 119 and 146 DEGs were found in T1D and T2D, respectively. PNP and CCR1 have great diagnostic value for both T1D and T2D. MGAM and NAMPT had great diagnostic value for T2D. Conclusion: Our finding provided clues for developing biomarkers for diabetes.

Integrative analysis of transcriptome-wide association study and mRNA expression profiles identifies candidate genes associated with autism spectrum disorders.

Autism spectrum disorders (ASD) are a group of highly heritable psychiatric syndromes with high prevalence. The genetic mechanism of ASD remains elusive now. Here we conducted a transcriptome-wide association study (TWAS) of ASD. The GWAS summary data of ASD was driven from the Psychiatric Genomics Consortium (PGC) portal, totally involving 5,305 ASD cases and 5,305 controls. FUSION software was applied to the GWAS summary data for tissue-related TWAS of ASD considering brain and blood. The ASD associated genes identified by TWAS were further validated by mRNA expression profiling of ASD and the Simons Foundation for Autism Research (SFARI) Gene tool. DAVID 6.8 was used to perform gene ontology (GO) enrichment analysis of ASD associated genes identified by TWAS. TWAS identified 85 genes with TWAS P value <0.05 for ASD. Further comparing the 85 genes with the differentially expressed genes identified by mRNA expression profiling of ASD patients found 5 overlapped genes, including MUTYH (PTWAS = 0.0460, PmRNA = 0.0040), ARHGAP27 (PTWAS = 0.0100, PmRNA = 0.0016), GCA (PTWAS = 0.0480, PmRNA = 0.0063), CCDC14 (PTWAS = 0.0067, PmRNA = 0.0035), and MED15 (PTWAS = 0.0324, PmRNA = 0.0092). Gene Ontology (GO) enrichment analysis of the genes identified by TWAS detected 10 significant GO terms, such as mitochondrion (P = 0.0051), NAD or NADH binding (P = 0.0169), mitochondrial part (P = 0.0386) and 2-oxoglutarate metabolic process (P = 0.0399). In conclusion, this study identified multiple ASD associated genes and gene sets, providing novel clues for revealing the pathogenesis of ASD. Autism Research 2019, 12: 33-38. © 2018 International Society for Autism Research, Wiley Periodicals, Inc. LAY SUMMARY: Recent genetic studies of autism spectrum disorders (ASD) have found multiple ASD related genes. However, the results of these studies were hardly replicated with each other, providing limited clues for exploring the genetic mechanism of ASD. This study detected a group of candidate genes showing transcriptome-wide associations with ASD. These results may provide novel clues for revealing the pathogenesis of ASD.

Use of integrative epigenetic and mRNA expression analyses to identify significantly changed genes and functional pathways in osteoarthritic cartilage.

AIM: Osteoarthritis (OA) is caused by complex interactions between genetic and environmental factors. Epigenetic mechanisms control the expression of genes and are likely to regulate the OA transcriptome. We performed integrative genomic analyses to define methylation-gene expression relationships in osteoarthritic cartilage. PATIENTS AND METHODS: Genome-wide DNA methylation profiling of articular cartilage from five patients with OA of the knee and five healthy controls was conducted using the Illumina Infinium HumanMethylation450 BeadChip (Illumina, San Diego, California). Other independent genome-wide mRNA expression profiles of articular cartilage from three patients with OA and three healthy controls were obtained from the Gene Expression Omnibus (GEO) database. Integrative pathway enrichment analysis of DNA methylation and mRNA expression profiles was performed using integrated analysis of cross-platform microarray and pathway software. Gene ontology (GO) analysis was conducted using the Database for Annotation, Visualization and Integrated Discovery (DAVID). RESULTS: We identified 1265 differentially methylated genes, of which 145 are associated with significant changes in gene expression, such as DLX5, NCOR2 and AXIN2 (all p-values of both DNA methylation and mRNA expression < 0.05). Pathway enrichment analysis identified 26 OA-associated pathways, such as mitogen-activated protein kinase (MAPK) signalling pathway (p = 6.25 × 10-4), phosphatidylinositol (PI) signalling system (p = 4.38 × 10-3), hypoxia-inducible factor 1 (HIF-1) signalling pathway (p = 8.63 × 10-3 pantothenate and coenzyme A (CoA) biosynthesis (p = 0.017), ErbB signalling pathway (p = 0.024), inositol phosphate (IP) metabolism (p = 0.025), and calcium signalling pathway (p = 0.032). CONCLUSION: We identified a group of genes and biological pathwayswhich were significantly different in both DNA methylation and mRNA expression profiles between patients with OA and controls. These results may provide new clues for clarifying the mechanisms involved in the development of OA.Cite this article: A. He, Y. Ning, Y. Wen, Y. Cai, K. Xu, Y. Cai, J. Han, L. Liu, Y. Du, X. Liang, P. Li, Q. Fan, J. Hao, X. Wang, X. Guo, T. Ma, F. Zhang. Use of integrative epigenetic and mRNA expression analyses to identify significantly changed genes and functional pathways in osteoarthritic cartilage. Bone Joint Res 2018;7:343-350. DOI: 10.1302/2046-3758.75.BJR-2017-0284.R1.

Identification and characterization of atypical 2-cysteine peroxiredoxins from mud crab Scylla paramamosain: The first evidence of two peroxiredoxin 5 genes in non-primate species and their involvement in immune defense against pathogen infection.

Peroxiredoxin 5 (Prx5) belongs to a novel family of evolutionarily conserved antioxidant proteins that protect cells against various oxidative stresses. Generally, no more than one Prx5 transcript had been reported in non-primate species. In this study, two Prx5 genes (coined as SpPrx5-1 and SpPrx5-2) were firstly isolated from the mud crab, Scylla paramamosain, through RT-PCR and RACE methods. The open reading frame of SpPrx5-1 and SpPrx5-2 were 561 bp and 429 bp in length, encoding 186 and 142 amino acids polypeptide, respectively. Both the conserved signatures of peroxiredoxin catalytic center and Prx5-specific domain were identified in SpPrx5-1 and SpPrx5-2. Phylogenetic analysis indicated that both SpPrx5 clustered together with other animal Prx proteins and were classified into Prx5 subfamily. Tissue-specific expression analysis revealed that both SpPrx5-1 and SpPrx5-2 were ubiquitously expressed, highest in hepatopancreas, and showed remarkably similar transcription patterns. Quantitative RT-PCR analysis exhibited that both SpPrx5 genes changed dramatically in hepatopancreas, although showing different expression profiles, after virus-analog poly (I:C) or Vibrio alginolyticus challenge. The expression levels of both SpPrx5s were significantly enhanced in hepatopancreas after poly (I:C) stimulation, while SpPrx5-2 exhibited a more prompt response than SpPrx5-1. Nevertheless, the expression levels of both SpPrx5s were significantly reduced in hepatopancreas after Vibrio alginolyticus challenge in which SpPrx5-1 showed a more prompt response than SpPrx5-2. These results suggested the involvement of SpPrx5s in responses against viral and bacterial infections and further highlighted their functional importance in the immune system of Scylla paramamosain.

Predicting Functional MicroRNA-mRNA Interactions.

MicroRNAs (miRNAs) are small RNA molecules that play key regulatory roles in general biological processes and disease pathogenesis. These small RNA molecules interact with their target mRNAs to induce mRNA degradation and/or inhibit the translation of mRNAs into proteins. Therefore, identifying miRNA targets is an essential step to fully understand the regulatory effects of miRNAs. Here, we describe a regularized regression approach that integrates the sequence information with the miRNA and mRNA expression profiles for detecting miRNA targets. This method takes into account the full spectrum of gene sequence features of miRNA targets, including the thermodynamic stability, the accessibility energy, and the context features of the target sites,. Given these sequence features for each putative miRNA-mRNA interaction and their expression values, this model is able to quantify the down-regulation effect of each miRNA on its targets while simultaneously estimating the contribution of each sequence feature for predicting functional miRNA-mRNA interactions.

Ultrasonic irradiation of low intensity with a mode of sweeping frequency enhances the membrane permeability and cell growth rate of Candida tropicalis.

Here we report the enhancement of both cellular permeability and cell growth rate of Candida tropicalis after treatment with the ultrasonic irradiation of low intensity using a mode of sweeping frequency (UILS) generated by a self-developed ultrasonic device in our lab. After the ultrasonic treatment, remarkable biomass enhancement of the yeast was observed; the hyphae became significantly longer; the seeped cellular protein and nucleic acid from the yeast increased and the cellular Ca2+ content became lower. Illumina transcriptome sequencing showed that the ultrasonic treatment affected the expression of genes involved in diverse cellular components, biological processes and molecular functions. RT-PCR and Western blotting further confirmed the up-/down-regulation of genes in the ultrasound-treated yeasts. The optimal conditions of the ultrasonic treatment for the maximum biomass addition were determined as follows: the yeast was treated for 1h at the mid logarithmic phase, the frequency was 28±2kHz and the power density was 120W/L. Under these conditions, the Candida tropicalis biomass increased by 142.5% compared with the untreated yeast.

The detection of risk pathways, regulated by miRNAs, via the integration of sample-matched miRNA-mRNA profiles and pathway structure.

The use of genome-wide, sample-matched miRNA (miRNAs)-mRNA expression data provides a powerful tool for the investigation of miRNAs and genes involved in diseases. The identification of miRNA-regulated pathways has been crucial for analysis of the role of miRNAs. However, the classical identification method fails to consider the structural information of pathways and the regulation of miRNAs simultaneously. We proposed a method that simultaneously integrated the change in gene expression and structural information in order to identify pathways. Our method used fold changes in miRNAs and gene products, along with the quantification of the regulatory effect on target genes, to measure the change in gene expression. Topological characteristics were investigated to measure the influence of gene products on entire pathways. Through the analysis of multiple myeloma and prostate cancer expression data, our method was proven to be effective and reliable in identifying disease risk pathways that are regulated by miRNAs. Further analysis showed that the structure of a pathway plays a crucial role in the recognition of the pathway as a factor in disease risk.

Identifying dysfunctional miRNA-mRNA regulatory modules by inverse activation, cofunction, and high interconnection of target genes: a case study of glioblastoma.

BACKGROUND: Accumulating evidence demonstrates that complex diseases may arise from cooperative effects of multiple dysfunctional miRNAs. Thus, identifying abnormal functions cooperatively regulated by multiple miRNAs is useful for understanding the pathogenesis of complex diseases. METHODS: In this study, we proposed a multistep method to identify dysfunctional miRNA-mRNA regulatory modules (dMiMRMs) in a specific disease, in which a group of miRNAs cooperatively regulate a group of target genes involved in a specific function. We identified dysfunctional miRNAs, which were differentially expressed and inversely regulated most of their target genes, by integrating paired miRNA and mRNA expression profiles and miRNA target information. Then, we identified cooperative functional units, in each of which a pair of miRNAs cooperatively repressed function-enriched and highly interconnected target genes. Finally, the cooperative functional units were assembled into dMiMRMs. RESULTS: We applied our method to glioblastoma (GBM) and identified GBM-associated dMiMRMs at the population, subtype, and individual levels. We identified 5 common dMiMRMs using all GBM samples, 3 of which were associated with the prognosis in patients with GBM and were better predictors of prognosis than were miRNAs or mRNAs alone. By applying our approach to GBM subtypes, we found consistent dMiMRMs across GBM subtypes, and some subtype-specific dMiMRMs were observed. Furthermore, personalized dMiMRMs were identified, suggesting significant individual differences in different patients with GBM. CONCLUSIONS: Our method provides the capability to identify miRNA-mediated dysfunctional mechanisms underlying complex diseases.