Optimal gene expression analysis by microarrays

Tailoring Escherichia coli BL21 (DE3) for preferential xylose utilization via metabolic and regulatory engineering

Xylose is the most abundant pentose in nature. However, it is usually obtained in mixtures with glucose, leading to carbon catabolite repression in many microorganisms. Among E. coli lineages, significant metabolic and regulatory differences exist, requiring distinct metabolic engineering strategies to develop a xylose-selective phenotype in the strains W, K-12, and C. In this study, strain ES02 was engineered from Escherichia coli BL21 (DE3) as a xylose-selective strain by deleting the glk, ptsG, and manZ genes. However, when grown in a mixture of xylose and glucose, this strain's specific growth rate and xylose consumption rate decreased by about 50% compared to cultures with only xylose. A modified version of the xylose-responsive transcriptional activator XylRQ31K was utilized to overcome this issue. The resulting strain ES04 (BL21 (DE3) Δglk, ΔmanZ, ΔptsG, xylR::Kmr, lacZ::xylRC91A-Gmr) efficiently used xylose as carbon source either alone or in a mixture with glucose, with a specific xylose consumption rate 75% higher than that of the wild-type strain BL21(DE3). Unexpectedly, strain ES04 partially recovers the ability to grow and consume glucose at a low rate, preferentially consuming xylose over glucose in sugar mixtures, revealing an altered carbon catabolite repression phenotype. Transcriptomics analysis suggested that glucose assimilation in this strain was related to the overexpression of the galactitol operon gatDCBAZY. Further inactivation of this operon confirmed its participation in glucose assimilation. KEY POINTS: • XylRQ31K alleviates carbon catabolite repression in the xylose-selective strain ES04. • Galactitol operon overexpression in ES04 links to partial glucose utilization. • ES04 strain preferentially uses xylose over glucose, revealing altered CCR.

Inferring gene regulatory networks of ALS from blood transcriptome profiles

One of the most robust approaches to the prediction of causal driver genes of complex diseases is to apply reverse engineering methods to infer a gene regulatory network (GRN) from gene expression profiles (GEPs). In this work, we analysed 794 GEPs of 1117 human whole-blood samples from Amyotrophic Lateral Sclerosis (ALS) patients and healthy subjects reported in the GSE112681 dataset. GRNs for ALS and healthy individuals were reconstructed by ARACNe-AP (Algorithm for the Reconstruction of Accurate Cellular Networks - Adaptive Partitioning). In order to examine phenotypic differences in the ALS population surveyed, several datasets were built by arranging GEPs according to sex, spinal or bulbar onset, and survival time. The designed reverse engineering methodology identified a significant number of potential ALS-promoting mechanisms and putative transcriptional biomarkers that were previously unknown. In particular, the characterization of ALS phenotypic networks by pathway enrichment analysis has identified a gender-specific disease signature, namely network activation related to the radiation damage response, reported in the networks of bulbar and female ALS patients. Also, focusing on a smaller interaction network, we selected some hub genes to investigate their inferred pathological and healthy subnetworks. The inferred GRNs revealed the interconnection of the four selected hub genes (TP53, SOD1, ALS2, VDAC3) with p53-mediated pathways, suggesting the potential neurovascular response to ALS neuroinflammation. In addition to being well consistent with literature data, our results provide a novel integrated view of ALS transcriptional regulators, expanding information on the possible mechanisms underlying ALS and also offering important insights for diagnostic purposes and for developing possible therapies for a disease yet incurable.

Temporal gene expression profiling during early-stage traumatic temporomandibular joint bony ankylosis in a sheep model

BACKGROUND

Investigating the molecular biology underpinning the early-stage of traumatic temporomandibular joint (TMJ) ankylosis is crucial for discovering new ways to prevent the disease. This study aimed to explore the dynamic changes of transcriptome from the intra-articular hematoma or the newly generated ankylosed callus during the onset and early progression of TMJ ankylosis.

METHODS

Based on a well-established sheep model of TMJ bony ankylosis, the genome-wide microarray data were obtained from samples at postoperative Days 1, 4, 7, 9, 11, 14 and 28, with intra-articular hematoma at Day 1 serving as controls. Fold changes in gene expression values were measured, and genes were identified via clustering based on time series analysis and further categorised into three major temporal classes: increased, variable and decreased expression groups. The genes in these three temporal groups were further analysed to reveal pathways and establish their biological significance.

RESULTS

Osteoblastic and angiogenetic genes were found to be significantly expressed in the increased expression group. Genes linked to inflammation and osteoclasts were found in the decreased expression group. The various biological processes and pathways related to each temporal expression group were identified, and the increased expression group comprised genes exclusively involved in the following pathways: Hippo signaling pathway, Wnt signaling pathway and Rap 1 signaling pathway. The decreased expression group comprised genes exclusively involved in immune-related pathways and osteoclast differentiation. The variable expression group consisted of genes associated with DNA replication, DNA repair and DNA recombination. Significant biological pathways and transcription factors expressed at each time point postoperatively were also identified.

CONCLUSIONS

These data, for the first time, presented the temporal gene expression profiling and reveal the important process of molecular biology in the early-stage of traumatic TMJ bony ankylosis. The findings might contributed to identifying potential targets for the treatment of TMJ ankylosis.

Genes related to allergen exposure in allergic rhinitis: a gene-chip-based study in a mouse model

BACKGROUND

The typical clinical symptoms of allergic rhinitis (AR) are known to be associated with allergen exposure; however, the underlying mechanisms are not fully understood. We wanted to gain a comprehensive view of the molecular mechanisms related to allergen exposure in a well-controlled mouse model of AR.

METHODS

An OVA-induced AR model was developed. Two hours and 4 weeks after the last OVA challenge, AR symptoms and local immune responses were assessed. At the same time, differentially expressed genes (DEG) in nasal mucosa were identified by gene expression microarray and further analyzed by bioinformatics methods. Verification of DEG was done by quantitative RT-PCR and immunohistochemistry.

RESULTS

The number of nasal rubbings and sneezes, serum OVA-specific IgE concentrations, and the number of neutrophils and eosinophils in the nasal mucosa were significantly increased at 2 h and decreased at 4 weeks after the last allergen challenge compared to controls. A total of 2119 DEG were identified, and their expression dynamics were clustered into 8 profiles. Enriched functions in Profile 5, which had a similar trend to clinical features, were mainly related to inflammatory and immune response to environmental factors, eosinophils and neutrophils chemotaxis, and cell migration. Gene co-expression Network for genes from profile 5 identified BCL3, NFKB2, SOCS3, and CD53 having a higher degree. Profile 6 showed persistence of inflammatory and immune response at 4 weeks after the last allergen challenge. Olfactory and coagulation functions were enriched mainly in profiles with downward trends.

CONCLUSIONS

A wide range of genes with sequential cooperative action were identified to be associated with allergen exposure in AR. BCL3 may be the most vital in symptoms manifestation. Moreover, some inflammatory responses persisted for a period after allergen exposure, supporting a new treatment strategy of targeting inflammation out of season. This study may contribute to a better understanding of AR pathogenesis and provide potential therapeutic targets for AR patients.

Identification of Competing Endogenous RNAs (ceRNAs) Network Associated with Drought Tolerance in Medicago truncatula with Rhizobium Symbiosis

Drought, bringing the risks of agricultural production losses, is becoming a globally environmental stress. Previous results suggested that legumes with nodules exhibited superior drought tolerance compared with the non-nodule group. To investigate the molecular mechanism of rhizobium symbiosis impacting drought tolerance, transcriptome and sRNAome sequencing were performed to identify the potential mRNA-miRNA-ncRNA dynamic network. Our results revealed that seedlings with active nodules exhibited enhanced drought tolerance by reserving energy, synthesizing N-glycans, and medicating systemic acquired resistance due to the early effects of symbiotic nitrogen fixation (SNF) triggered in contrast to the drought susceptible with inactive nodules. The improved drought tolerance might be involved in the decreased expression levels of miRNA such as mtr_miR169l-5p, mtr_miR398b, and mtr_miR398c and its target genes in seedlings with active nodules. Based on the negative expression pattern between miRNA and its target genes, we constructed an mRNA-miR169l-ncRNA ceRNA network. During severe drought stress, the lncRNA alternative splicings TCONS_00049507 and TCONS_00049510 competitively interacted with mtr_miR169l-5p, which upregulated the expression of NUCLEAR FACTOR-Y (NF-Y) transcription factor subfamily NF-YA genes MtNF-YA2 and MtNF-YA3 to regulate their downstream drought-response genes. Our results emphasized the importance of SNF plants affecting drought tolerance. In conclusion, our work provides insight into ceRNA involvement in rhizobium symbiosis contributing to drought tolerance and provides molecular evidence for future study.

Prognosis of lasso-like penalized Cox models with tumor profiling improves prediction over clinical data alone and benefits from bi-dimensional pre-screening

BACKGROUND

Prediction of patient survival from tumor molecular '-omics' data is a key step toward personalized medicine. Cox models performed on RNA profiling datasets are popular for clinical outcome predictions. But these models are applied in the context of "high dimension", as the number p of covariates (gene expressions) greatly exceeds the number n of patients and e of events. Thus, pre-screening together with penalization methods are widely used for dimensional reduction.

METHODS

In the present paper, (i) we benchmark the performance of the lasso penalization and three variants (i.e., ridge, elastic net, adaptive elastic net) on 16 cancers from TCGA after pre-screening, (ii) we propose a bi-dimensional pre-screening procedure based on both gene variability and p-values from single variable Cox models to predict survival, and (iii) we compare our results with iterative sure independence screening (ISIS).

RESULTS

First, we show that integration of mRNA-seq data with clinical data improves predictions over clinical data alone. Second, our bi-dimensional pre-screening procedure can only improve, in moderation, the C-index and/or the integrated Brier score, while excluding irrelevant genes for prediction. We demonstrate that the different penalization methods reached comparable prediction performances, with slight differences among datasets. Finally, we provide advice in the case of multi-omics data integration.

CONCLUSIONS

Tumor profiles convey more prognostic information than clinical variables such as stage for many cancer subtypes. Lasso and Ridge penalizations perform similarly than Elastic Net penalizations for Cox models in high-dimension. Pre-screening of the top 200 genes in term of single variable Cox model p-values is a practical way to reduce dimension, which may be particularly useful when integrating multi-omics.

FGF21 at physiological concentrations regulates vascular endothelial cell function through multiple pathways

Cardiovascular diseases are closely associated with dysfunction of vascular endothelial cells (VECs), which can be influenced by various intrinsic and extrinsic factors, including fibroblast growth factor 21 (FGF21), but the effects of serum FGF21 on VECs remain unclear. We performed a cross-sectional study nested within a prospective cohort to assess the range of physiological concentrations of fasting serum FGF21 in 212 healthy individuals. We also treated human umbilical VECs (HUVECs) with recombinant FGF21 at different concentrations. The effects of FGF21 treatment on glycolysis, nitric oxide release and reduction of intracellular reactive oxygen species were assessed. The cells were also collected for RNA transcriptomic sequencing to investigate the potential mechanisms induced by FGF21 treatment. In addition, the roles of SIRT1 in the regulation of FGF21 were evaluated by SIRT1 knockdown. The results showed that the serum FGF21 concentration in healthy individuals ranged from 15.70 to 499.96 pg/mL and was positively correlated with age and pulse wave velocity. FGF21 at 400 pg/mL was sufficient to enhance glycolysis, increase nitric oxide release and protect cells from H₂O₂-induced oxidative damage. The upregulated genes after FGF21 treatment were mostly enriched in metabolic pathways, whereas the downregulated genes were mostly enriched in inflammation and apoptosis signaling pathways. Moreover, SIRT1 may be involved in the regulation of some genes by FGF21. In conclusion, our data indicate that FGF21 at a level within the physiological concentration range has a beneficial effect on HUVECs and that this effect may partly depend on the regulation of SIRT1.

Guanylate-binding protein 2b regulates the AMPK/mTOR/ULK1 signalling pathway to induce autophagy during Mycobacterium bovis infection

Autophagic isolation and degradation of intracellular pathogens are employed by host cells as primary innate immune defense mechanisms to control intercellular M. bovis infection. In this study, RNA-Seq technology was used to obtain the total mRNA from bone marrow-derived macrophages (BMDMs) infected with M. bovis at 6 and 24 h after infection. One of the differential genes, GBP2b, was also investigated. Analysis of the significant pathway involved in GBP2b-coexpressed mRNA demonstrated that GBP2b was associated with autophagy and autophagy-related mammalian target of rapamycin (mTOR) signaling and AMP-activated protein kinase (AMPK) signaling. The results of in vivo and in vitro experiments showed significant up-regulation of GBP2b during M. bovis infection. For in vitro validation, small interfering RNA-GBP2b plasmids were transfected into BMDMs and RAW264.7 cells lines to down-regulate the expression of GBP2b. The results showed that the down-regulation of GBP2b impaired autophagy via the AMPK/mTOR/ULK1 pathway, thereby promoting the intracellular survival of M. bovis. Further studies revealed that the activation of AMPK signaling was essential for the regulation of autophagy during M. bovis infection. These findings expand the understanding of how GBP2b regulates autophagy and suggest that GBP2b may be a potential target for the treatment of diseases caused by M. bovis.

Estrogen-Related Receptor γ Agonist DY131 Ameliorates Lipopolysaccharide-Induced Acute Liver Injury

Sepsis-associated liver dysfunction remains a challenge in clinical practice with high mortality and limited specific therapies. DY131 is a pharmacological agonist of the orphan receptor estrogen-related receptor (ERR) γ which plays a crucial role in regulating energy generation, oxidative metabolism, cell apoptosis, inflammatory responses, etc. However, its role in acute liver injury is unknown. In this study, we evaluated the effect of DY131 on lipopolysaccharide (LPS)-induced liver injury. Mice were pretreated with DY131 through intraperitoneal injection at a dose of 5 mg/kg/day for 3 days prior to LPS challenge (10 mg/kg). 24 h later, they were anesthetized and sacrificed. Blood and liver tissues were collected for further studies. In a separate experiment, mice were treated with saline (vehicle) or DY131 for 3 days to evaluate the toxicity of DY131. We found that ERRγ was downregulated in the liver tissues from LPS-treated mice. Pretreatment with DY131 ameliorated LPS-induced liver injury as demonstrated by reduced liver enzyme release (ALT, AST, and LDH), improved liver morphological damage, and attenuated oxidative stress, inflammation and apoptosis. Meanwhile, DY131 had no significant side effects on hepatic and renal functions in mice. Finally, transcriptomics analysis revealed that the dysregulated pathways associated with inflammation and metabolism were significantly reversed by DY131 in LPS-treated mice, providing more evidence in favor of the protective effect of DY131 against LPS-induced liver injury. Altogether, these findings highlighted the protective effect of DY131 on LPS-induced hepatotoxicity possibly via suppressing oxidative stress, inflammation, and apoptosis.

Integrated analysis of competing endogenous RNA in esophageal carcinoma

Background

The Competing endogenous RNA (CeRNA) network plays important roles in the development and progression of multiple human cancers. Increasing attention has been paid to CeRNA in esophageal carcinoma (ESCA).

Methods

We explored The Cancer Genome Atlas (TCGA) database and then analyzed the RNAs of 142 samples to obtain long non-coding RNAs (lncRNAs), micro RNAs (miRNAs), and messenger RNAs (mRNAs) with different expression trends alongside the progress of ESCA. A series test of cluster (STC) analysis was carried out to identify a set of unique model expression tendencies. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used to validate the function of key genes that were obtained from the STC analysis.

Results

Through our analysis, 272 lncRNAs, 87 miRNAs, and 692 mRNAs showed upward expression or downward expression trends, and these molecules were tightly involved in cell cycle, pathways in cancer, metabolic processes, and protein phosphorylation, among others. Ultimately, we constructed a CeRNA network containing a total of 71 lncRNAs, 56 miRNAs, and 125 mRNAs. The overall survival (OS) was analyzed using univariate Cox regression analysis to clarify the relationship between these key molecules from the CeRNA network and the prognosis of ESCA patients. Through survival analysis, we finally screened out two lncRNAs (DLEU2, RP11-890B15.3), three miRNAs (miR-26b-3p, miR-92a-3p, miR-324-5p), and one mRNA (SIK2) as crucial prognostic factors for ESCA.

Conclusions

The novel CeRNA network that we constructed will provide new novel prognostic biomarkers and therapeutic targets for patients with ESCA.

Integrated transcriptomics and metabolomics analysis of catechins, caffeine and theanine biosynthesis in tea plant (Camellia sinensis) over the course of seasons

BACKGROUND

Catechins, caffeine, and theanine as three important metabolites in the tea leaves play essential roles in the formation of specific taste and shows potential health benefits to humans. However, the knowledge on the dynamic changes of these metabolites content over seasons, as well as the candidate regulatory factors, remains largely undetermined.

RESULTS

An integrated transcriptomic and metabolomic approach was used to analyze the dynamic changes of three mainly metabolites including catechins, caffeine, and theanine, and to explore the potential influencing factors associated with these dynamic changes over the course of seasons. We found that the catechins abundance was higher in Summer than that in Spring and Autumn, and the theanine abundance was significantly higher in Spring than that in Summer and Autumn, whereas caffeine exhibited no significant changes over three seasons. Transcriptomics analysis suggested that genes in photosynthesis pathway were significantly down-regulated which might in linkage to the formation of different phenotypes and metabolites content in the tea leaves of varied seasons. Fifty-six copies of nine genes in catechins biosynthesis, 30 copies of 10 genes in caffeine biosynthesis, and 12 copies of six genes in theanine biosynthesis were detected. The correlative analysis further presented that eight genes can be regulated by transcription factors, and highly correlated with the changes of metabolites abundance in tea-leaves.

CONCLUSION

Sunshine intensity as a key factor can affect photosynthesis of tea plants, further affect the expression of major Transcription factors (TFs) and structural genes in, and finally resulted in the various amounts of catechins, caffeine and theaine in tea-leaves over three seasons. These findings provide new insights into abundance and influencing factors of metabolites of tea in different seasons, and further our understanding in the formation of flavor, nutrition and medicinal function.

Drp1 regulates mitochondrial dysfunction and dysregulated metabolism in ischemic injury via Clec16a-, BAX-, and GSH- pathways

The adaptation of mitochondrial homeostasis to ischemic injury is not fully understood. Here, we studied the role of dynamin-related protein 1 (Drp1) in this process. We found that mitochondrial morphology was altered in the early stage of ischemic injury while mitochondrial dysfunction occurred in the late stage of ischemia. Drp1 appeared to inhibit mitophagy by upregulating mito-Clec16a, which suppressed mito-Parkin recruitment and subsequently impaired the formation of autophagosomes in vascular tissues after ischemic injury. Moreover, ischemia-induced Drp1 activation enhanced apoptosis through inducing mitochondrial translocation of BAX and thereby increasing release of Cytochrome C to activate caspase-3/-9 signalling. Furthermore, Drp1 mediated metabolic disorders and inhibited the levels of mitochondrial glutathione to impair free radical scavenging, leading to further increases in ROS and the exacerbation of mitochondrial dysfunction after ischemic injury. Together, our data suggest a critical role for Drp1 in ischemic injury.

Dynamic transcriptome profiling exploring cold tolerance in forensically important blow fly, Aldrichina grahami (Diptera: Calliphoridae)

Background

Aldrichina grahami (Diptera: Calliphoridae) is a forensically important fly, which has been widely applied to practical legal investigations. Unlike other necrophagous flies, A. grahami exhibits cold tolerance which helps to maintain its activity during low-temperature months, when other species are usually not active. Hence, A. grahami is considered an important forensic insect especially in cold seasons. In this study, we aim to explore the molecular mechanisms of cold tolerance of A. grahami through transcriptome.

Results

We collected eggs and larvae (first-instar, second-instar and third-instar) at three different temperatures (4 °C, 12 °C and 20 °C) and performed RNA-seq analyses. The differentially expressed genes (DEGs) associated with the cold-tolerance were screened out. The Venn analysis of DEGs from egg to third-instar larvae at three different temperatures showed there were 9 common genes. Candidate biological processes and genes were identified which refer to growth, and development of different temperatures, especially the chitin and cuticle metabolic process. The series-clusters showed crucial and unique trends when the temperature changed. Moreover, by comparing the results of growth and developmental transcriptomes from different temperatures, we found that DEGs belonging to the family of larval cuticle proteins (LCP), pupal cuticle protein (CUP), and heat shock proteins (HSP) have certain differences.

Conclusions

This study identified functional genes and showed differences in the expression pattern of diverse temperatures. The DEGs series-clusters with increasing or decreasing trends were analyzed which may play an important role in cold-tolerance. Moreover, the findings in LCP, CUP and HSP showed more possible modulations in a cold environment. This work will provide valuable information for the future investigation of the molecular mechanism of cold tolerance in A. grahami.

Measuring similarity between gene interaction profiles

Background

Gene and protein interaction data are often represented as interaction networks, where nodes stand for genes or gene products and each edge stands for a relationship between a pair of gene nodes. Commonly, that relationship within a pair is specified by high similarity between profiles (vectors) of experimentally defined interactions of each of the two genes with all other genes in the genome; only gene pairs that interact with similar sets of genes are linked by an edge in the network. The tight groups of genes/gene products that work together in a cell can be discovered by the analysis of those complex networks.

Results

We show that the choice of the similarity measure between pairs of gene vectors impacts the properties of networks and of gene modules detected within them. We re-analyzed well-studied data on yeast genetic interactions, constructed four genetic networks using four different similarity measures, and detected gene modules in each network using the same algorithm. The four networks induced different numbers of putative functional gene modules, and each similarity measure induced some unique modules. In an example of a putative functional connection suggested by comparing genetic interaction vectors, we predict a link between SUN-domain proteins and protein glycosylation in the endoplasmic reticulum.

Conclusions

The discovery of molecular modules in genetic networks is sensitive to the way of measuring similarity between profiles of gene interactions in a cell. In the absence of a formal way to choose the “best” measure, it is advisable to explore the measures with different mathematical properties, which may identify different sets of connections between genes.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3024-x) contains supplementary material, which is available to authorized users.

lncRNA Gm10451 regulates PTIP to facilitate iPSCs-derived β-like cell differentiation by targeting miR-338-3p as a ceRNA

iPSCs-derived insulin-producing cell transplantation is a promising strategy for diabetes therapy. Although there have been many protocols of mature, glucose-responsive β cells induced in vitro over the past few years, many underlying problems remain to be resolved. As a crucial regulator, long noncoding RNAs (lncRNAs) participate in numerous biological processes, including the maintenance of pluripotency, and stem cell differentiation. In this study, we identified a novel lncRNA Gm10451 as a functional regulator for β-like cell differentiation. Localized to the cytoplasm, Gm10451 regulates histone H3K4 methyltransferase complex PTIP to facilitate Insulin⁺/Nkx6.1⁺ β-like cell differentiation by targeting miR-338-3p as a competing endogenous RNA (ceRNA). miR-338-3p has also been shown to suppress Nkx6.1⁺ early-stage β-like cell differentiation by targeting PTIP. Following transplantation into streptozotocin (STZ)-mice, Gm10451 loss in β-like cells prevented the expression of mature β-cell makers, such as Insulin, Nkx6.1, and Mafa. Accordingly, hyperglycemia in the mice was not resolved. Taken together, this study provides an efficient epigenetic target for generating more mature and functional iPSCs-derived β-like cells. We anticipate that pancreatic organoids, which are generated from human stem cells, biological materials, and epigenetic modifications, can be used in the future as a novel diabetes treatment option.

Role of vimentin in modulating immune cell apoptosis and inflammatory responses in sepsis

New diagnostic biomarkers or therapeutic targets for sepsis have substantial significance for critical care medicine. In this study, 192 differentially expressed proteins were selected through iTRAQ. Based on cluster analysis of protein expression dynamics and protein-protein interactions, hemopexin, vimentin, and heat shock protein 90 were selected for further investigation. It was demonstrated that serum vimentin (VIM) levels were significantly increased in patients with sepsis and septic shock compared to controls and that VIM expression was significantly increased in lymphocytes isolated from septic shock and sepsis patients compared to controls. Moreover, a nonsurvivor group had higher serum VIM levels and VIM expression in lymphocytes. Caspase-3 was significantly upregulated in Jurkat T cells lacking VIM and when exposed to LPS compared to control cells. In contrast, caspase-3 was reduced nearly 40% in cells over-expressing VIM. IL-2, IL-10 and IFN-α levels were significantly decreased in cells lacking VIM compared to control cells, whereas they were not significantly altered in cells over-expressing VIM. These findings suggest that VIM modulates lymphocyte apoptosis and inflammatory responses and that VIM could be a new target for the diagnosis and prognostic prediction of patients with sepsis or septic shock.

Heterotypic CAF-tumor spheroids promote early peritoneal metastatis of ovarian cancer

The study provides insights in HGSOC by identifying that ascitic CAFs selectively recruit ITGA5^high ascitic tumor cells to form heterotypic spheroids named metastatic units (MUs), which actively engage in peritoneal metastasis, discriminates HGSOC from LGSOC, and act as therapeutic targets in hampering OC metastasis.

High-grade serous ovarian cancer (HGSOC) is hallmarked by early onset of peritoneal dissemination, which distinguishes it from low-grade serous ovarian cancer (LGSOC). Here, we describe the aggressive nature of HGSOC ascitic tumor cells (ATCs) characterized by integrin α5^high (ITGA5^high) ATCs, which are prone to forming heterotypic spheroids with fibroblasts. We term these aggregates as metastatic units (MUs) in HGSOC for their advantageous metastatic capacity and active involvement in early peritoneal dissemination. Intriguingly, fibroblasts inside MUs support ATC survival and guide their peritoneal invasion before becoming essential components of the tumor stroma in newly formed metastases. Cancer-associated fibroblasts (CAFs) recruit ITGA5^high ATCs to form MUs, which further sustain ATC ITGA5 expression by EGF secretion. Notably, LGSOC is largely devoid of CAFs and the resultant MUs, which might explain its metastatic delay. These findings identify a specialized MU architecture that amplifies the tumor–stroma interaction and promotes transcoelomic metastasis in HGSOC, providing the basis for stromal fibroblast-oriented interventions in hampering OC peritoneal propagation.

Identification of novel lncRNAs involved in the pathogenesis of childhood acute lymphoblastic leukemia

This study aimed to explore novel long non-coding RNAs (lncRNAs) and the underlying mechanisms involved in childhood acute lymphoblastic leukemia (cALL). The GSE67684 dataset was downloaded from the Gene Expression Omnibus. Differentially expressed genes (DEGs) and lncRNAs (DELs) between Days 0, 8, 15 and 33 were isolated using random variance model corrective analysis of variance. Overlapping DEGs and DELs were clustered using Cluster 3.0. Bio-functional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Interactions between lncRNAs and mRNAs were calculated using dynamic simulations, and interactions among mRNAs were predicted using the STRING database. lncRNA-mRNA and protein-protein interaction (PPI) networks were visualized using Cytoscape. Subsequently, the expression levels of lncRNAs in biological samples from children with or without cALL were validated using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). A total of 593 overlapping DEGs and 21 DELs were identified. After clustering, Profile 26 exhibited a continuously increasing temporal trend, whereas Profile 1 exhibited a continuous decreasing trend. Upregulated DEGs were significantly enriched in 1,825 GO terms and 166 KEGG pathways, whereas downregulated DEGs were significantly enriched in 196 GO terms and 90 KEGG pathways. The lncRNAs NONHSAT027612.2 and NONHSAT134556.2 were the top two regulators in the lncRNA-mRNA network. Toll-like receptor 4, cathepsin G, nucleotide-binding oligomerization domain containing 2 and cathepsin S may be considered the hub genes of the PPI network. RT-qPCR results indicated that the expression levels of the lncRNAs NONHSAT027612.2 and NONHSAT134556.2 were significantly elevated in the blood and bone marrow of patients with cALL compared with the controls. In conclusion, the lncRNAs NONHSAT027612.2 and NONHSAT134556.2 may serve important roles in the pathogenesis of cALL via regulating immune response-associated pathways.

Dynamic transcriptome profiling towards understanding the morphogenesis and development of diverse feather in domestic duck

Background

Feathers with complex and fine structure are hallmark avian integument appendages, which have contributed significantly to the survival and breeding for birds. Here, we aimed to explore the differentiation, morphogenesis and development of diverse feathers in the domestic duck.

Results

Transcriptome profiles of skin owing feather follicle from two body parts at three physiological stages were constructed to understand the molecular network and excavate the candidate genes associated with the development of plumulaceous and flight feather structures. The venn analysis of differentially expressed genes (DEGs) between abdomen and wing skin tissues at three developmental stages showed that 38 genes owing identical differentially expression pattern. Together, our data suggest that feather morphological and structural diversity can be possibly related to the homeobox proteins. The key series-clusters, many candidate biological processes and genes were identified for the morphogenesis, growth and development of two feather types. Through comparing the results of developmental transcriptomes from plumulaceous and flight feather, we found that DEGs belonging to the family of WNT, FGF and BMP have certain differences; even the consistent DEGs of skin and feather follicle transcriptomes from abdomen and wing have the different expression patterns.

Conclusions

Overall, this study detected many functional genes and showed differences in the molecular mechanisms of diverse feather developments. The findings in WNT, FGF and BMP, which were consistent with biological experiments, showed more possible complex modulations. A correlative role of HOX genes was also suggested but future biological verification experiments are required. This work provided valuable information for subsequent research on the morphogenesis of feathers.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4778-7) contains supplementary material, which is available to authorized users.

Integrated analysis of microRNA and mRNA expression profiles in rats with selenium deficiency and identification of associated miRNA-mRNA network

Selenium deficiency is closely related with various type of cardiovascular disease. However, the miRNA-mRNA regulatory network in Selenium deficiency related cardiac change remains to be understand. In the present study, a reliable Selenium deficiency rat model was established and confirmed by pathological and biochemical examination. The mRNA and miRNA expression profiles were conducted by microarray technology. Gene Ontology (GO) Analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Analysis was performed to investigate the function of targeted genes, and the relationship between miRNA and mRNA was studied by network analysis. A total of 4931 mRNAs and 119 miRNAs was differentially expressed between any two groups (control group, low-selenium group and selenium supplementation group). GO and KEGG pathway analysis of selected miRNAs target genes found that selenium deficiency was related to several different biological processes. Furthermore, a miRNA-mRNA regulatory network was conducted to illustrate the interaction of miRNAs and these targeted genes. In conclusion, our present study provides a new insight that potential molecular mechanism of Selenium deficiency was a multiply miRNAs and mRNA caused biological change.

Comprehensive gene and microRNA expression profiling on cardiovascular system in zebrafish co-exposured of SiNPs and MeHg

Air pollution has been shown to increase cardiovascular diseases. However, little attention has been paid to the combined effects of PM and air pollutants on the cardiovascular system. To explore this, a high-throughput sequencing technology was used to determine combined effects of silica nanoparticles (SiNPs) and MeHg in zebrafish. Our study demonstrated that SiNPs and MeHg co-exposure could cause significant changes in mRNA and miRNA expression patterns in zebrafish. The differentially expressed (DE) genes in profiles 17 and 26 of STC analysis suggest that SiNPs and MeHg co-exposure had more proinflammatory and cardiovascular toxicity in zebrafish than single exposure. Major gene functions associated with cardiovascular system in the co-exposed zebrafish were discerned from the dynamic-gene-network, including stxbp1a, celf4, ahr1b and bai2. In addition, the prominently expressed pathway of cardiac muscle contraction was targeted by 3 DE miRNAs identified by the miRNA-pathway-network (dre-miR-7147, dre-miR-26a and dre-miR-375), which included 23 DE genes. This study presents a global view of the combined SiNPs and MeHg toxicity on the dynamic expression of both mRNAs and miRNAs in zebrafish, and could serve as fundamental research clues for future studies, especially on cardiovascular system toxicity.

Genome-Wide Search for Competing Endogenous RNAs Responsible for the Effects Induced by Ebola Virus Replication and Transcription Using a trVLP System

Understanding how infected cells respond to Ebola virus (EBOV) and how this response changes during the process of viral replication and transcription are very important for establishing effective antiviral strategies. In this study, we conducted a genome-wide screen to identify long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), micro RNAs (miRNAs), and mRNAs differentially expressed during replication and transcription using a tetracistronic transcription and replication-competent virus-like particle (trVLP) system that models the life cycle of EBOV in 293T cells. To characterize the expression patterns of these differentially expressed RNAs, we performed a series cluster analysis, and up- or down-regulated genes were selected to establish a gene co-expression network. Competing endogenous RNA (ceRNA) networks based on the RNAs responsible for the effects induced by EBOV replication and transcription in human cells, including circRNAs, lncRNAs, miRNAs, and mRNAs, were constructed for the first time. Based on these networks, the interaction details of circRNA-chr19 were explored. Our results demonstrated that circRNA-chr19 targeting miR-30b-3p regulated CLDN18 expression by functioning as a ceRNA. These findings may have important implications for further studies of the mechanisms of EBOV replication and transcription. These RNAs potentially have important functions and may be promising targets for EBOV therapy.

Global gene expression analysis combined with a genomics approach for the identification of signal transduction networks involved in postnatal mouse myocardial proliferation and development

Mammalian cardiomyocytes may permanently lose their ability to proliferate after birth. Therefore, studying the proliferation and growth arrest of cardiomyocytes during the postnatal period may enhance the current understanding regarding this molecular mechanism. The present study identified the differentially expressed genes in hearts obtained from 24 h-old mice, which contain proliferative cardiomyocytes; 7-day-old mice, in which the cardiomyocytes are undergoing a proliferative burst; and 10-week-old mice, which contain growth-arrested cardiomyocytes, using global gene expression analysis. Furthermore, myocardial proliferation and growth arrest were analyzed from numerous perspectives, including Gene Ontology annotation, cluster analysis, pathway enrichment and network construction. The results of a Gene Ontology analysis indicated that, with increasing age, enriched gene function was not only associated with cell cycle, cell division and mitosis, but was also associated with metabolic processes and protein synthesis. In the pathway analysis, 'cell cycle', proliferation pathways, such as the 'PI3K-AKT signaling pathway', and 'metabolic pathways' were well represented. Notably, the cluster analysis revealed that bone morphogenetic protein (BMP)1, BMP10, cyclin E2, E2F transcription factor 1 and insulin like growth factor 1 exhibited increased expression in hearts obtained from 7-day-old mice. In addition, the signal transduction pathway associated with the cell cycle was identified. The present study primarily focused on genes with altered expression, including downregulated anaphase promoting complex subunit 1, cell division cycle (CDC20), cyclin dependent kinase 1, MYC proto-oncogene, bHLH transcription factor and CDC25C, and upregulated growth arrest and DNA damage inducible α in 10-week group, which may serve important roles in postnatal myocardial cell cycle arrest. In conclusion, these data may provide important information regarding myocardial proliferation and development.

Long non-coding RNA expression profiles in different severity EV71-infected hand foot and mouth disease patients

Enterovirus 71 (EV71) is associated with the severe hand foot and mouth disease (HFMD) outcomes, however the host-virus interaction mechanism and the pathogenesis remain poorly understood. Long non-coding RNAs (lncRNAs) are involved in variety physiological and pathological processes, but the functions of lncRNAs in EV71 infection remain elusive. Here we profiled the expression of lncRNAs in peripheral blood mononuclear cells (PBMCs) from EV71-infected mild patients, severe patients as well as the healthy controls, and identified 8541 lncRNAs were differentially expressed. Focused on the dynamic changed lncRNAs, we performed systematic bioinformatics analysis with Series Test of Cluster (STC) algorithm, Gene Ontology (GO) analysis, pathway analysis and lncRNA-mRNA co-expression network analysis, and revealed the potential functions and related pathways of these lncRNAs were associated with immunity and inflammation during the clinical process of EV71-infected HFMD. Among the significant dynamic changed lncRNAs, ten lncRNAs were screened whose expression were further validated in EV71-infected mild patients, severe patients and healthy control. These results shed light on the potential roles of lncRNAs in EV71-infected HFMD, especially in distinguishing the mild and severe cases for early diagnose and treatment, moreover, provide deeper insight into the mechanism of EV71-induced immune and inflammatory responses, as well as the pathogenesis of the imbalanced inflammation in severe EV71 infection.

•

LncRNAs expression profile were analyzed with EV71-infected HFMD patients' PBMCs.

•

Dynamic changed lncRNAs were analyzed in different severity EV71-infected patients.

•

Most of the lncRNAs were related to innate immune and inflammatory responses.

•

The candidate lncRNAs may serve as potential markers for EV71-induced severe HFMD.

miR-181d-5p promotes neurite outgrowth in PC12 Cells via PI3K/Akt pathway

INTRODUCTION

miRNAs dysregulate in spinal cord injury (SCI) and have been demonstrated to play a crucial role in neurite outgrowth. However, the underlying mechanism remains elusive. In this study, we constructed a mouse model of SCI, extracted RNA from injured spinal cord tissue for the use of microarray assay. miR-181d-5p which is one of the most significantly expressed miRNAs in miRNA-mRNA network, abundantly expressed in center system and highly conserved across different spices, was chosen for our further study.

AIMS

To demonstrate whether miR-181d-5p can promote neurite outgrowth in PC12 cells via PI3K/Akt signaling pathway, we performed function analysis of miR-181d-5p with LV-miR-181d-5p and LV-sh-GFP to infect PC12 cells.

RESULTS

Through microarray assay analysis, we totally found 262 significantly expressed miRNAs and 2973 target genes in SCI and observed that their expression dynamically changed postinjury. Here, we provided enough evidences that the overexpression of miR181d-5p significantly decreased the expression of PTEN, upregulated p-Akt expression, increased neurite outgrowth-related proteins (GAP-43 and NF-200) and synaptic vesicle-related proteins (Synapsin and PSD95), and then promoted neurite outgrowth in PC12 cells. Furthermore, we confirmed that miR-181d-5p could directly target to the 3'-UTR of PTEN mRNA through dual-luciferase report assay.

CONCLUSIONS

Our study supports that aberrant expression of miRNAs is involved in the pathogenesis of SCI, miR-181d-5p plays an important role in neurite growth in PC12 cells via PI3K/Akt signaling pathway and may be a candidate target for the treatment of SCI in the future.

Dynamic gene expression analysis in a H1N1 influenza virus mouse pneumonia model

H1N1, a major pathogenic subtype of influenza A virus, causes a respiratory infection in humans and livestock that can range from a mild infection to more severe pneumonia associated with acute respiratory distress syndrome. Understanding the dynamic changes in the genome and the related functional changes induced by H1N1 influenza virus infection is essential to elucidating the pathogenesis of this virus and thereby determining strategies to prevent future outbreaks. In this study, we filtered the significantly expressed genes in mouse pneumonia using mRNA microarray analysis. Using STC analysis, seven significant gene clusters were revealed, and using STC-GO analysis, we explored the significant functions of these seven gene clusters. The results revealed GOs related to H1N1 virus-induced inflammatory and immune functions, including innate immune response, inflammatory response, specific immune response, and cellular response to interferon-beta. Furthermore, the dynamic regulation relationships of the key genes in mouse pneumonia were revealed by dynamic gene network analysis, and the most important genes were filtered, including Dhx58, Cxcl10, Cxcl11, Zbp1, Ifit1, Ifih1, Trim25, Mx2, Oas2, Cd274, Irgm1, and Irf7. These results suggested that during mouse pneumonia, changes in the expression of gene clusters and the complex interactions among genes lead to significant changes in function. Dynamic gene expression analysis revealed key genes that performed important functions. These results are a prelude to advancements in mouse H1N1 influenza virus infection biology, as well as the use of mice as a model organism for human H1N1 influenza virus infection studies.

Transcriptome analysis of endometrial tissues following GnRH agonist treatment in a mouse adenomyosis model

Purpose

Adenomyosis is a common, benign gynecological condition of the female reproductive tract characterized by heavy menstrual bleeding and dysmenorrhea. Gonadotropin-releasing hormone (GnRH) agonists are one of the medications used in adenomyosis treatment; however, their underlying mechanisms are poorly understood. Moreover, it is difficult to obtain endometrial samples from women undergoing such treatment. To overcome this, we generated an adenomyosis mouse model, which we treated with an GnRH agonist to determine its effect on pregnancy outcomes. We also analyzed endometrial gene expression following GnRH agonist treatment to determine the mechanisms that may affect pregnancy outcome in individuals with adenomyosis.

Methods

Neonatal female mice were divided into a control group, an untreated adenomyosis group, and an adenomyosis group treated with a GnRH agonist (n=6 each). The pregnancy outcome was observed and compared among the groups. Then, three randomly chosen transcriptomes from endometrial tissues from day 4 of pregnancy were analyzed between the adenomyosis group and the GnRH agonist treatment group by RNA sequencing and quantitative reverse transcription polymerase chain reaction (PCR).

Results

The litter size was significantly smaller in the adenomyosis group than in the control group (7±0.28 vs 11±0.26; P<0.05). However, the average live litter size was increased (10±0.28 vs 7±0.28; P<0.05) after GnRH agonist treatment. Three hundred and fifty-nine genes were differentially expressed in the GnRH agonist-treated group compared with the untreated group (218 were downregulated and 141 were upregulated). Differentially expressed genes were related to diverse biological processes, including estrogen metabolism, cell cycle, and metabolite biosynthesis.

Conclusion

GnRH agonist treatment appears to improve the pregnancy outcome of adenomyosis in a mouse model. Besides pituitary down-regulation, other possible mechanisms such as the regulation of cell proliferation may play a role in this. These new insights into GnRH agonist mechanisms will be useful for future adenomyosis treatment.

Ectopic overexpression of filamin C scaffolds MEK1/2 and ERK1/2 to promote the progression of human hepatocellular carcinoma

Hepatocellular carcinoma (HCC) invasion and metastasis are mediated by a complicated signal transduction network and downstream cytoskeletal and adhesion molecules. In this study, a microarray-based analysis revealed a dramatic increase in filamin C (FLNC), which is commonly expressed in muscle rather than in liver cells, in the two metastatic HCC cell lines MHCC97L and HCCLM3. Clinicopathological studies showed that increased FLNC expression was associated with microvascular invasion and poor prognosis. Specific hypomethylation was identified within the FLNC promoter region in HCC cell lines and patient tumor samples, which might contribute to the ectopic overexpression of FLNC. FLNC downregulation inhibited cell migration and impaired cell proliferation and promoted apoptosis. Mechanistic studies suggested that FLNC interacts with mitogen-activated extracellular signal-regulated kinase 1/2 (MEK1/2) and extracellular signal-regulated kinase 1/2 (ERK1/2) and that FLNC downregulation inhibited MEK1/2 and ERK1/2 activation. Xenographic tumor transplantation experiments in nude mice further confirmed the role of FLNC in HCC progression and metastasis. Our results reveal a novel mechanism by which the cytoskeletal protein FLNC enhances the mitogen-activated protein kinase signaling pathway during tumorigenesis.

miRNA targeted signaling pathway in the early stage of denervated fast and slow muscle atrophy

Denervation often results in skeletal muscle atrophy. Different mechanisms seem to be involved in the determination between denervated slow and fast skeletal muscle atrophy. At the epigenetic level, miRNAs are thought to be highly involved in the pathophysiological progress of denervated muscles. We used miRNA microarrays to determine miRNA expression profiles from a typical slow muscle (soleus muscle) and a typical fast muscle (tibialis anterior muscle) at an early denervation stage in a rat model. Results showed that miR-206, miR-195, miR-23a, and miR-30e might be key factors in the transformation process from slow to fast muscle in denervated slow muscles. Additionally, certain miRNA molecules (miR-214, miR-221, miR-222, miR-152, miR-320, and Let-7e) could be key regulatory factors in the denervated atrophy process involved in fast muscle. Analysis of signaling pathway networks revealed the miRNA molecules that were responsible for regulating certain signaling pathways, which were the final targets (e.g., p38 MAPK pathway; Pax3/Pax7 regulates Utrophin and follistatin by HDAC4; IGF1/PI3K/Akt/mTOR pathway regulates atrogin-1 and MuRF1 expression via FoxO phosphorylation). Our results provide a better understanding of the mechanisms of denervated skeletal muscle pathophysiology.

Aberrant LncRNA Expression Profile in a Contusion Spinal Cord Injury Mouse Model

Long noncoding RNAs (LncRNAs) play a crucial role in cell growth, development, and various diseases related to the central nervous system. However, LncRNA differential expression profiles in spinal cord injury are yet to be reported. In this study, we profiled the expression pattern of LncRNAs using a microarray method in a contusion spinal cord injury (SCI) mouse model. Compared with a spinal cord without injury, few changes in LncRNA expression levels were noted 1 day after injury. The differential changes in LncRNA expression peaked 1 week after SCI and subsequently declined until 3 weeks after injury. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate the reliability of the microarray, demonstrating that the results were reliable. Gene ontology (GO) analysis indicated that differentially expressed mRNAs were involved in transport, cell adhesion, ion transport, and metabolic processes, among others. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the neuroactive ligand-receptor interaction, the PI3K-Akt signaling pathway, and focal adhesions were potentially implicated in SCI pathology. We constructed a dynamic LncRNA-mRNA network containing 264 LncRNAs and 949 mRNAs to elucidate the interactions between the LncRNAs and mRNAs. Overall, the results from this study indicate for the first time that LncRNAs are differentially expressed in a contusion SCI mouse model.

Non-monotonic dose-response effect of bisphenol A on rare minnow Gobiocypris rarus ovarian development

Bisphenol A (BPA) is widely spread in the environment, and can cause various reproductive disrupting effects on different organisms, including fish. Our previous published study showed that BPA has non-monotonic (inverted U-shaped) dose-response effect on rare minnow Gobiocypris rarus ovarian weight at different concentrations. To investigate the potential mechanism, we exposed female rare minnow to 1, 15 and 225 µg L(-1) BPA for 7 days in the present study. The levels of vitellogenin (Vtg), sex hormones, hydrogen peroxide (H2O2), glutathione (GSH) and triglyceride (TG) were measured. RNA-seq of ovary tissues was also performed. Result showed that Vtg, sex hormone and TG levels showed an inverted U-shaped increased response, while H2O2 and GSH levels showed a U-shaped inhibited response. RNA-seq data showed that many genes involved in lipid metabolism, oxidative stress, and proteolysis processes were altered. The change of Vtg, H2O2, GSH and TG levels was possibly related to the altered sex hormone levels. Sex hormone's direct effect, Vtg accumulation, TG accumulation and oxidative stress induced proteolysis may contribute to the change of ovary weight.

MicroRNA expression profiles and networks in mouse lung infected with H1N1 influenza virus

Influenza A viruses can cause localized outbreaks and worldwide pandemics, owing to their high transmissibility and wide host range. As such, they are among the major diseases that cause human death. However, the molecular changes induced by influenza A virus infection in lung tissue are not entirely clear. Changes in microRNA (miRNA) expression occur in many pathological and physiological processes, and influenza A virus infection has been shown to alter miRNA expression in cultured cells and animal models. In this study, we mined key miRNAs closely related to influenza A virus infection and explored cellular regulatory mechanisms against influenza A virus infection, by building networks among miRNAs and genes, gene ontologies (GOs), and pathways. In this study, miRNAs and mRNAs induced by H1N1 influenza virus infection were measured by gene chips, and we found that 82 miRNAs and 3371 mRNAs were differentially expressed. The 82 miRNAs were further analyzed with the series test of cluster (STC) analysis. Three of the 16 cluster profiles identified by STC, which include 46 miRNAs in the three profiles, changed significantly. Using potential target genes of the 46 miRNAs, we looked for intersections of these genes with 3371 differentially expressed mRNAs; 719 intersection genes were identified. Based on the GO or KEGG databases, we attained GOs or pathways for all of the above intersection genes. Fisher's and χ (2) test were used to calculate p value and false discovery rate (FDR), and according to the standard of p < 0.001, 241 GOs and 76 pathways were filtered. Based on these data, miRNA-gene, miRNA-GO, and miRNA-pathway networks were built. We then extracted three classes of GOs (related to inflammatory and immune response, cell cycle, proliferation and apoptosis, and signal transduction) to build three subgraphs, and pathways strictly related with H1N1 influenza virus infection were filtered to extract a subgraph of the miRNA-pathway network. Last, according to the pathway analysis and miRNA-pathway network analysis, 17 miRNAs were found to be associated with the "influenza A" pathway. This study provides the most complete miRNAome profiles, and the most detailed miRNA regulatory networks to date, and is the first to report the most important 17 miRNAs closely related with the pathway of influenza A. These results are a prelude to advancements in mouse H1N1 influenza virus infection biology and the use of mice as a model for human H1N1 influenza virus infection studies.

PSTK is a novel gene associated with early lung injury in Paraquat Poisoning

AIMS

Paraquat Poisoning (PQ) can cause illness and death, and its main causes of mortality are acute respiratory failure and lung fibrosis. Early recognition of this condition and early treatment are vital. Thus, it is of importance to target the key genes controlling pathogenesis in the early stage of PQ.

MAIN METHODS

C57BL/6 mice were used for Paraquat intragastric administration as a model of PQ. Following a gene chip-based screening, the change of gene expression in the lung was further validated by bioinformatic analyses, co-expression network construction and real-time RT-PCR, Western blot and immunofluorescence assays.

KEY FINDINGS

2287 genes with differential expression were identified at the very early stage of PQ. From these, 76 genes that were linked to mitochondrion function were further pursued. Among these genes, PSTK was a phosphorylase kinase which serves a protective role in oxidative stress lung damage. PSTK was the central gene in a 30-gene network that is important for mitochondrial complex I assembly, mitochondrial apoptosis and mitochondrial fatty acid beta-oxidation, suggesting that they could conceivably be related to the pathogenesis of PQ induced lung damage. Lastly, we confirmed that PSTK was lowered in rodent lungs following PQ.

SIGNIFICANCE

PSTK emerges as a central gene in a network of mitochondrial function genes in PQ exposed mice. The functional role of PSTK in PQ induced lung injury warrants further examination.

Distance-based classifiers as potential diagnostic and prediction tools for human diseases

Typically, gene expression biomarkers are being discovered in course of high-throughput experiments, for example, RNAseq or microarray profiling. Analytic pipelines that extract so-called signatures suffer from the "Dimensionality curse": the number of genes expressed exceeds the number of patients we can enroll in the study and use to train the discriminator algorithm. Hence, problems with the reproducibility of gene signatures are more common than not; when the algorithm is executed using a different training set, the resulting diagnostic signature may turn out to be completely different.

In this paper we propose an alternative novel approach which takes into account quantifiable expression levels of all genes assayed. In our analysis, the cumulative gene expression pattern of an individual patient is represented as a point in the multidimensional space formed by all gene expression profiles assayed in given system, where the clusters of "normal samples" and "affected samples" and defined. The degree of separation of the given sample from the space occupied by "normal samples" reflects the drift of the sample away from homeostasis in the course of development of the pathophysiological process that underly the disease. The outlined approach was validated using the publicly available glioma dataset deposited in Rembrandt and associated with survival data. Additionally, the applicability of the distance analysis to the classification of non-malignant sampled was tested using psoriatic lesions and non-lesional matched controls as a model.

Keywords: biomarkers; clustering; human diseases; RNA

Protein expression profiling during wallerian degeneration after rat sciatic nerve injury

INTRODUCTION

Wallerian degeneration (WD) is an important area of research in modern neuroscience. Many protein expressions are regulated by differentially expressed genes in WD, but the precise mechanisms are elusive.

METHODS

In this study, we profiled differentially expressed proteins in WD after rat sciatic nerve injury using an antibody array.

RESULTS

Functional analysis positively identified cell proliferation, regulation of cell proliferation, and immune system processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed molecular networks related mainly to cytokine-cytokine receptor interaction, the mitogen-activated proteinkinase (MAPK) signaling pathway, apoptosis, the toll-like receptor (TLR) signaling pathway, and the Janus kinase (Jak) - signal transducer and activator of transcription (STAT) signaling pathway. Interactions between these differential proteins were well established and regulated by the key factors transforming growth factor beta 1 (TGF-β1), toll-like receptor 4 (TLR4), Fas ligand (FasL), and 5'-AMP-activated protein kinase catalytic subunit alpha-1 (PRKAA1).

CONCLUSIONS

These results provide information related to functional analysis of differentially expressed genes during WD.

Differential gene expression profiling analysis in workers occupationally exposed to benzene

Benzene is an important industrial chemical and an environmental contaminant, but the pathogenesis of hematotoxicity induced by chronic occupational benzene exposure remains to be elucidated. To gain an insight into the molecular mechanisms and new biomarkers, microarray analysis was used to identify the differentially expressed mRNA critical for benzene hematotoxicity. RNA was extracted from four chronic benzene poisoning patients occupationally exposed to benzene, three benzene-exposed workers without clinical symptoms and three health controls without benzene exposure and mRNA expression profiling was performed using Gene Chip Human Gene 2.0ST Arrays. Analysis of mRNA expression profiles were conducted to identify key genes, biological processes and pathways by the series test of cluster (STC), STC-Gene Ontology analysis (STC-GO), pathway analysis and Signal-net.

PRINCIPAL FINDINGS

1) 1661 differentially expressed mRNAs were identified and assigned to sixteen model profiles. 2) Profiles 14, 10, 11, 1 and 15 are the predominant expression profiles which were involved in immune response, inflammatory response, chemotaxis, defense response, anti-apoptosis and signal transduction. 3) The importance of immune response at benzene hematotoxicity is highlighted by several immune-related signaling pathways such as B/T cell receptor signaling pathway, acute myeloid leukemia, hematopoietic cell lineage and natural killer cell mediated cytotoxicity. 4) Signet analysis showed that PIK3R1, PIK3CG, PIK3R2, GNAI3, KRAS, NRAS, NFKB1, HLA-DMA, and HLA-DMB were key genes involved in benzene hematotoxicity. These genes might be new biomarkers for the prevention and early diagnosis of benzene poisoning. This is a preliminary study that paves the way to further functional study to understand the underlying regulatory mechanisms.

Preparation and use of reverse protein microarrays

Reverse-phase protein array (RPPA) is a multiplex, high-throughput proteomic technique for profiling the activation status of signal transduction pathways involved in cancer survival and progression, potentially allowing for identification of new biomarkers and drug targets. On RPPA, the entire patient proteome is immobilized on a spot and single proteins can be quantified across a set of samples, spotted on the same array, with high specificity and sensitivity. Array immunostaining and signal amplification systems are used to generate a signal proportional to the concentration of the analyte. Dedicated scanners and software are used to detect spots, measure intensity, subtract background, normalize signal, and generate a numeric value as output. The generated output file is then analyzed using several different bioinformatic and biostatistical tools. In this unit, the RPPA procedure is described in depth, from sample handling and preparation to data analysis, with particular emphasis on tissue sample analysis.

Clinical Study of Effects of Jian Ji Ning, a Chinese Herbal Medicine Compound Preparation, in Treating Patients with Myasthenia Gravis via the Regulation of Differential MicroRNAs Expression in Serum

Myasthenia gravis (MG) is an autoimmune disease, of which the pathogenesis has remained unclear. At present, MG does not have any effective treatment with minor side effects. Jian Ji Ning (JJN), a traditional Chinese medicine formula consisting of 11 medicinal plants, has been used in the treatment of MG for many years. The present study aims to determine if the Chinese herbal medicine JJN could lighten the clinical symptoms of patients with MG via the regulation of differential microRNAs (miRNAs) expression in serum. JJN should be orally administered twice a day for 6 months. In the efficacy evaluation adopting the Quantitative Myasthenia Gravis Score (QMG), we found that JJN could improve the clinical symptoms of patients with MG more effectively. Besides, we found that JJN could regulate differential miRNAs expression in serum of patients with MG. Accordingly, we speculate that the effects of JJN on improving clinical symptoms and blood test indicators of patients with MG may be due to its inhibition of apoptotic pathways of some immune cells and its connection with the regulation of serum miRNAs of some patients. In conclusion, we believe that JJN has a reliable curative effect on patients with MG-induced neuropathologic changes.

miR-429 identified by dynamic transcriptome analysis is a new candidate biomarker for colorectal cancer prognosis

Colorectal cancer (CRC) is a common malignant gastrointestinal cancer. Efforts for preventive and personalized medicine have intensified in the last decade with attention to novel forms of biomarkers. In the present study, microRNA and genetic analyses were performed in tandem for differential transcriptome profiling between primary tumors with or without nodes or distant metastases. Serial Test Cluster (STC) analysis demonstrated that 20 genes and two microRNAs showed distinctive expression patterns associated with the tumor, node, and metastasis (TNM) stage. The selected target genes were characterized by GO and Pathway analysis. A microRNA-target gene network analysis showed that miR-429 resided in the center of the network, indicating that miR-429 might serve important roles in the development of CRC. Real-time PCR and tissue microarrays showed that miR-429 had a dynamic expression pattern during the CRC progression stage, and was significantly downregulated in stage II and stage III clinical progression. The low expression of miR-429 was correlated with poor prognosis for CRC. Taken together, miR-429 warrant further clinical translation research as a candidate biomarker for CRC prognosis. Additional downstream targets and attendant gene function also need to be discerned to design a sound critical path to personalized medicine for persons susceptible to, or diagnosed with CRC.

Transcriptome changes associated with delayed flower senescence on transgenic petunia by inducing expression of etr1-1, a mutant ethylene receptor

Flowers of ethylene-sensitive ornamental plants transformed with ethylene-insensitive 1-1(etr1-1), a mutant ethylene receptor first isolated from Arabidopsis, are known to have longer shelf lives. We have generated petunia plants in which the etr1-1 gene was over-expressed under the control of a chemically-inducible promoter, which would allow expression of etr1-1 to be initiated at the desired time and stage of development. Here, we showed that transgenic plants grew and developed normally without a chemical inducer. Semi-quantitative RT-PCR demonstrated that the abundance of transcripts of Arabidopsis etr1-1 gene was substantially induced in flowers with 30 μM dexamethasone (DEX). Consequently, t he life of the flowers was almost doubled and the peak of ethylene production was delayed. We compared gene expression changes of petals with DEX to those without DEX at 24 h and 48 h by microarray. Our results indicated that transcripts of many putative genes encoding transcription factors were down-regulated by etr1-1 induced expression at the early stage. In addition, putative genes involved in gibberellin biosynthesis, response to jasmonic acid/gibberellins stimulus, cell wall modification, ethylene biosynthesis, and cell death were down-regulated associating with etr1-1 induced expression. We investigated time-course gene expression profiles and found two profiles which displayed totally opposite expression patterns under these two treatments. In these profiles, 'the regulation of transcription' was predominant in GO categories. Taking all results together, we concluded those transcription factors down-regulated at early stage might exert a major role in regulating the senescence process which were consequently characterized by cell wall modification and cell death.

MicroRNA-gene expression network in murine liver during Schistosoma japonicum infection

BACKGROUND

Schistosomiasis japonica remains a significant public health problem in China and Southeast Asian countries. The most typical and serious outcome of the chronic oriental schistosomiasis is the progressive granuloma and fibrosis in the host liver, which has been a major medical challenge. However, the molecular mechanism underling the hepatic pathogenesis is still not clear.

METHODOLOGY AND PRINCIPAL FINDINGS

Using microarrays, we quantified the temporal gene expression profiles in the liver of Schistosoma japonicum-infected BALB/c mice at 15, 30, and 45 day post infection (dpi) with that from uninfected mice as controls. Gene expression alternation associated with liver damage was observed in the initial phase of infection (dpi 15), which became more magnificent with the onset of egg-laying. Up-regulated genes were dominantly associated with inflammatory infiltration, whereas down-regulated genes primarily led to the hepatic functional disorders. Simultaneously, microRNA profiles from the same samples were decoded by Solexa sequencing. More than 130 miRNAs were differentially expressed in murine liver during S. japonicum infection. MiRNAs significantly dysregulated in the mid-phase of infection (dpi 30), such as mmu-miR-146b and mmu-miR-155, may relate to the regulation of hepatic inflammatory responses, whereas miRNAs exhibiting a peak expression in the late phase of infection (dpi 45), such as mmu-miR-223, mmu-miR-146a/b, mmu-miR-155, mmu-miR-34c, mmu-miR-199, and mmu-miR-134, may represent a molecular signature of the development of schistosomal hepatopathy. Further, a dynamic miRNA-gene co-expression network in the progression of infection was constructed.

CONCLUSIONS AND SIGNIFICANCE

This study presents a global view of dynamic expression of both mRNA and miRNA transcripts in murine liver during S. japonicum infection, and highlights that miRNAs may play a variety of regulatory roles in balancing the immune responses during the development of hepatic pathology. The data provide robust information for further researches on the pathogenesis and molecular events of hepatopathy induced by schistosome eggs.

Expression changes and bioinformatic analysis of Wallerian degeneration after sciatic nerve injury in rat

Wallerian degeneration (WD) remains an important research topic. Many genes are differentially expressed during the process of WD, but the precise mechanisms responsible for these differentiations are not completely understood. In this study, we used microarrays to analyze the expression changes of the distal nerve stump at 0, 1, 4, 7, 14, 21 and 28 days after sciatic nerve injury in rats. The data revealed 6 076 differentially-expressed genes, with 23 types of expression, specifically enriched in genes associated with nerve development and axonogenesis, cytokine biosynthesis, cell differentiation, cytokine/chemokine production, neuron differentiation, cytokinesis, phosphorylation and axon regeneration. Kyoto Encyclopedia of Genes and Genomes pathway analysis gave findings related mainly to the MAPK signaling pathway, the Jak-STAT signaling pathway, the cell cycle, cytokine-cytokine receptor interaction, the p53 signaling pathway and the Wnt signaling pathway. Some key factors were NGF, MAG, CNTF, CTNNA2, p53, JAK2, PLCB1, STAT3, BDNF, PRKC, collagen II, FGF, THBS4, TNC and c-Src, which were further validated by real-time quantitative PCR, Western blot, and immunohistochemistry. Our findings contribute to a better understanding of the functional analysis of differentially-expressed genes in WD and may shed light on the molecular mechanisms of nerve degeneration and regeneration.

Genome-wide search for the genes accountable for the induced resistance to HIV-1 infection in activated CD4+ T cells: apparent transcriptional signatures, co-expression networks and possible cellular processes

BACKGROUND

Upon co-stimulation with CD3/CD28 antibodies, activated CD4 + T cells were found to lose their susceptibility to HIV-1 infection, exhibiting an induced resistant phenotype. This rather unexpected phenomenon has been repeatedly confirmed but the underlying cell and molecular mechanisms are still unknown.

METHODS

We first replicated the reported system using the specified Dynal beads with PHA/IL-2-stimulated and un-stimulated cells as controls. Genome-wide expression and analysis were then performed by using Agilent whole genome microarrays and established bioinformatics tools.

RESULTS

We showed that following CD3/CD28 co-stimulation, a homogeneous population emerged with uniform expression of activation markers CD25 and CD69 as well as a memory marker CD45RO at high levels. These cells differentially expressed 7,824 genes when compared with the controls on microarrays. Series-Cluster analysis identified 6 distinct expression profiles containing 1,345 genes as the representative signatures in the permissive and resistant cells. Of them, 245 (101 potentially permissive and 144 potentially resistant) were significant in gene ontology categories related to immune response, cell adhesion and metabolism. Co-expression networks analysis identified 137 "key regulatory" genes (84 potentially permissive and 53 potentially resistant), holding hub positions in the gene interactions. By mapping these genes on KEGG pathways, the predominance of actin cytoskeleton functions, proteasomes, and cell cycle arrest in induced resistance emerged. We also revealed an entire set of previously unreported novel genes for further mining and functional validation.

CONCLUSIONS

This initial microarray study will stimulate renewed interest in exploring this system and open new avenues for research into HIV-1 susceptibility and its reversal in target cells, serving as a foundation for the development of novel therapeutic and clinical treatments.

Gene expression profiling of the rat sciatic nerve in early Wallerian degeneration after injury

Wallerian degeneration is an important area of research in modern neuroscience. A large number of genes are differentially regulated in the various stages of Wallerian degeneration, especially during the early response. In this study, we analyzed gene expression in early Wallerian degeneration of the distal nerve stump at 0, 0.5, 1, 6, 12 and 24 hours after rat sciatic nerve injury using gene chip microarrays. We screened for differentially-expressed genes and gene expression patterns. We examined the data for Gene Ontology, and explored the Kyoto Encyclopedia of Genes and Genomes Pathway. This allowed us to identify key regulatory factors and recurrent network motifs. We identified 1 546 differentially-expressed genes and 21 distinct patterns of gene expression in early Wallerian degeneration, and an enrichment of genes associated with the immune response, acute inflammation, apoptosis, cell adhesion, ion transport and the extracellular matrix. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed components involved in the Jak-STAT, ErbB, transforming growth factor-β, T cell receptor and calcium signaling pathways. Key factors included interleukin-6, interleukin-1, integrin, c-sarcoma, carcinoembryonic antigen-related cell adhesion molecules, chemokine (C-C motif) ligand, matrix metalloproteinase, BH3 interacting domain death agonist, baculoviral IAP repeat-containing 3 and Rac. The data were validated with real-time quantitative PCR. This study provides a global view of gene expression profiles in early Wallerian degeneration of the rat sciatic nerve. Our findings provide insight into the molecular mechanisms underlying early Wallerian degeneration, and the regulation of nerve degeneration and regeneration.

Transcriptome analysis of rice mature root tissue and root tips in early development by massive parallel sequencing

Despite the major physiological dissimilarities between mature root regions and their tips, differences in their gene expression profiles remain largely unexplored. In this research, the transcriptome of rice (Oryza sativa L. subsp. japonica) mature root tissue and root tips was monitored using mRNA-Seq at two time points. Almost 50 million 76 bp reads were mapped onto the rice genome sequence, expression patterns for different tissues and time points were investigated, and at least 1106 novel transcriptionally active regions (nTARs) expressed in rice root tissue were detected. More than 30 000 genes were found to be expressed in rice roots, among which were 1761 root-enriched and 306 tip-enriched transcripts. Mature root tissue appears to respond more strongly to external stimuli than tips, showing a higher expression of, for instance, auxin-responsive and abscisic acid-responsive genes, as well as the phenylpropanoid pathway and photosynthesis upon light. The root tip-enriched transcripts are mainly involved in mitochondrial electron transport, organelle development, secondary metabolism, DNA replication and metabolism, translation, and cellular component organization. During root maturation, genes involved in cell wall biosynthesis and modification, response to oxidative stress, and secondary metabolism were activated. For some nTARs, a potential role in root development can be put forward based on homology to genes involved in CLAVATA signalling, cell cycle regulators, and hormone signalling. A subset of differentially expressed genes and novel transcripts was confirmed using (quantitative) reverse transcription-PCR. These results uncover previously unrecognized tissue-specific expression profiles and provide an interesting starting point to study the different regulation of transcribed regions of these tissues.

The cortical and striatal gene expression profile of 100 hz electroacupuncture treatment in 6-hydroxydopamine-induced Parkinson's disease model

Electroacupuncture (EA), especially high-frequency EA, has frequently been used as an alternative therapy for Parkinson disease (PD) and is reportedly effective for alleviating motor symptoms in patients and PD models. However, the molecular mechanism underlying its effectiveness is not completely understood. To implement a full-scale search for the targets of 100 Hz EA, we selected rat models treated with 6-hydroxydopamine into the unilateral MFB, which mimic end-stage PD. High-throughput microarray analysis was then used to uncover the regulated targets in the cortex and striatum after 4-week EA treatment. In the differentially regulated transcripts, the proportion of recovered expression profiles in the genes, the functional categories of targets in different profiles, and the affected pathways were analyzed. Our results suggested that the recovery of homeostasis in the transcript network and many regulated functional clusters in the cortex and striatum after EA treatment may contribute to the behavioral improvement of PD rats.

Gfi1b negatively regulates Rag expression directly and via the repression of FoxO1

Gfi1b negatively regulates Rag expression through direct binding to the Rag locus and through inhibition of Foxo1; mice lacking both Gfi1b and Gfi1 exhibit a block in B cell development.

Precise regulation of Rag (recombination-activating gene) expression is crucial to prevent genomic instability caused by the generation of Rag-mediated DNA breaks. Although mechanisms of Rag activation have been well characterized, the mechanism by which Rag expression is down-regulated in early B cell development has not been fully elucidated. Using a complementary DNA library screen, we identified the transcriptional repressor Gfi1b as negative regulator of the Rag locus. Expression of Gfi1b causes repression of Rag1 and Rag2 in cell lines and primary mouse cells. Conversely, Gfi1b-deficient cell lines exhibit increased Rag expression, double-strand breaks and recombination, and cell cycle defects. In primary cells, transcription of Gfi1b inversely correlates with Rag transcription, and simultaneous inactivation of Gfi1 and Gfi1b leads to an increase in Rag transcription early in B cell development. In addition, deletion of Gfi1 and Gfi1b in vivo results in a severe block in B cell development. Gfi1b orchestrates Rag repression via a dual mechanism. Direct binding of Gfi1b to a site 5′ of the B cell–specific Erag enhancer results in epigenetic changes in the Rag locus, whereas indirect inhibition is achieved through repression of the trans-activator Foxo1. Together, our experiments show that Gfi family members are essential for normal B cell development and play an important role in modulating expression of the V(D)J recombinase.