Digital Gene Expression Data Research Articles

Identification of early jasmonate-responsive genes in Taxus × media cells by analyzing time series digital gene expression data.

Jasmonate, an effective elicitor, can induce the biosynthesis of paclitaxel, a well-known anticancer drug, in Taxus cell culture. The jasmonate signaling pathway has been well studied in Arabidopsis, and many early jasmonate-responsive genes have beenfound to be involved in signaling pathway. In Taxus, only a few late jasmonate-responsive genes related to paclitaxel biosynthesis were identified. So, identification of early responsive genes and knowledge of the jasmonate signaling pathway are essential for understanding the effects of jasmonate on paclitaxel biosynthesis and for improving paclitaxel production in Taxus cells. In this study, total RNA of Taxus × media cells cultured in liquid medium was extracted after 0, 0.5, 3, and 24h of methyl jasmonate treatment. Three biological independent repetitions were performed. The 12 extracted RNA samples were integrated and sequenced on an Illumina HiSeq2500 platform using the paired-end method. A total of 45,583 transcript clusters were obtained by de novo assembly of the sequenced reads. Based on the transcriptome data, the digital gene expressions of each RNA sample were investigated. We found that after 0.5, 3, and 24h of methyl jasmonate treatment; 134, 1008, and 987 unigenes were differentially expressed. For the secondary metabolism pathways, phenylalanine pathway unigenes were responsive to jasmonate after 3h of treatment, while genes related to paclitaxel biosynthesis were induced after 0.5h of treatment. The digital gene expression levels of candidate genes related to paclitaxel biosynthesis were confirmed by qRT-PCR. Transcriptome sequencing and digital gene expression profiling identified early jasmonate-responsive genes in cultured Taxus × media cells. The comprehensive time series jasmonate-responsive gene expression data have provided transcriptome-wide information about the mechanism of paclitaxel biosynthesis regulation by jasmonate signaling.

Building a Competing Endogenous RNA Network to Find Potential Long Non-Coding RNA Biomarkers for Pheochromocytoma

Background/Aims: Accumulating evidence has shown that long non-coding RNAs (lncRNAs) in competing endogenous RNA (ceRNA) networks play crucial roles in tumor survival and patient prognosis; however, studies investigating ceRNA networks in pheochromocytoma (PCC) are lacking. In this study, we investigated the pathogenesis of PCC and whether lncRNAs acting through ceRNAs networks were associated with prognosis. Methods: A total of 183 PCC samples and 3 control samples from The Cancer Genome Atlas database were analyzed. The Empirical Analysis of Digital Gene Expression Data package in R (edgeR) was used to analyze differentially expressed RNAs. Biological processes and pathways functional enrichment analysis were performed based on the Database for Annotation, Visualization, and Integrated Discovery (DAVID) database. LncRNA/mRNA/miRNA ceRNA network was constructed by Cytoscape v3.0 software based on the differentially expressed RNAs Survival package in R was used to perform survival analysis. Results: In total, 554 differentially expressed lncRNAs, 1775 mRNAs and 40 miRNAs were selected for further analysis. Subsequently, 23 lncRNAs, 22 mRNAs, and 6 miRNAs were included in the constructed ceRNA network. Meanwhile, two of the 23 lncRNAs (C9orf147 and BSN-AS2) were identified as independent predictors of overall survival in PCC patients (P< 0.05). Conclusion: This study improves the understanding of lncRNA-related ceRNA networks in PCC and suggests that the lncRNAs C9orf147 and BSN-AS2 could be independent prognostic biomarkers and potential therapeutic targets for PCC.

Assessment of tris (1, 3-dichloro-2-propyl) phosphate toxicology in PC12 cells by using digital gene expression profiling

Tris (1,3-dichloro-2-propyl) phosphate (TDCIPP), one of the most universally used organophosphate flame retardants (OPFRs), is an environmental pollutant. However, limited information is available regarding its toxicity and environmental health risk. In the present study, PC12 cells provided a useful model for the evaluation of the toxic effects of TDCIPP. Exposure to 7.5, 15, 30, or 60 μM TDCIPP for 72 h inhibited cell viability, and enhanced cellular apoptosis and oxidative stress. To further explore the underlying mechanisms, digital gene expression (DGE) technology was used to identify early transcriptional changes following TDCIPP exposure. Expression of the transcripts of 161 genes was significantly altered upon treatment with TDCIPP. Functional and pathway analysis of the transcriptional profile demonstrated that genes showing significant TDCIPP-associated changes in expression were involved in the phosphatidylinositol 3-kinase (PI3K)/protein kinase B (Akt) signaling pathway, extracellular matrix–receptor interactions, protein digestion and absorption, and microRNAs in cancer. Using quantitative real-time PCR, we validated the differential expression of selected genes. These results showed that the expression profiles of cells exposed to 60 μM TDCIPP were consistent with the DGE data. Furthermore, western blotting showed that treatment with TDCIPP reduced the Bcl-2/Bax ratio and attenuated PI3K/Akt/Myc signaling. Taken together, these data suggest that TDCIPP exposure can reduce cell viability and induce apoptosis in PC12 cells by inhibiting activation of the PI3K/Akt/Myc signaling pathway. These observations provide valuable preliminary information regarding the mechanisms of TDCIPP-induced toxicity in PC12 cells and indicate that further study of the toxicity of other environmental OPFRs is warranted.

ITRAQ Mitoproteome Analysis Reveals Mechanisms of Programmed Cell Death in Arabidopsis thaliana Induced by Ochratoxin A.

Ochratoxin A (OTA) is one of the most common and dangerous mycotoxins in the world. Previous work indicated that OTA could elicit spontaneous HR-like lesions formation Arabidopsis thaliana, reactive oxygen species (ROS) play an important role in OTA toxicity, and their major endogenous source is mitochondria. However, there has been no evidence as to whether OTA induces directly PCD in plants until now. In this study, the presence of OTA in Arabidopsis thaliana leaves triggered accelerated respiration, increased production of mitochondrial ROS, the opening of ROS-dependent mitochondrial permeability transition pores and a decrease in mitochondrial membrane potential as well as the release of cytochrome c into the cytosol. There were 42 and 43 significantly differentially expressed proteins identified in response to exposure to OTA for 8 and 24 h, respectively, according to iTRAQ analysis. These proteins were mainly involved in perturbation of the mitochondrial electron transport chain, interfering with ATP synthesis and inducing PCD. Digital gene expression data at transcriptional level was consistent with the cell death induced by OTA being PCD. These results indicated that mitochondrial dysfunction was a prerequisite for OTA-induced PCD and the initiation and execution of PCD via a mitochondrial-mediated pathway.

DGE-seq analysis of MUR3-related Arabidopsis mutants provides insight into how dysfunctional xyloglucan affects cell elongation

Our previous study of the Arabidopsis mur3-3 mutant and mutant plants in which the mur3-3 phenotypes are suppressed (xxt2mur3-3, xxt5mur3-3, xxt1xxt2mur3-3 and 35Spro:XLT2:mur3-3) showed that hypocotyl cell elongation is decreased in plants that synthesize galactose-deficient xyloglucan. To obtain genome-wide insight into the transcriptome changes and regulatory networks that may be involved in this decreased elongation, we performed digital gene expression analyses of the etiolated hypocotyls of wild type (WT), mur3-3 and the four suppressor lines. Numerous differentially expressed genes (DEGs) were detected in comparisons between WT and mur3-3 (1423), xxt2mur3-3 and mur3-3 (675), xxt5mur3-3 and mur3-3 (1272), xxt1xxt2mur3-3 and mur3-3 (1197) and 35Spro:XLT2:mur3-3 vs mur3-3 (121). 550 overlapped DEGs were detected among WT vs mur3-3, xxt2mur3-3 vs mur3-3, xxt5mur3-3 vs mur3-3, and xxt1xxt2mur3-3 vs mur3-3 comparisons. These DEGs include 46 cell wall-related genes, 24 transcription factors, 6 hormone-related genes, 9 protein kinase genes and 9 aquaporin genes. The expression of all of the 550 overlapped genes is restored to near wild-type levels in the four mur3-3 suppressor lines. qRT-PCR of fifteen of these 550 genes showed that their expression levels are consistent with the digital gene expression data. Overexpression of some of these genes (XTH4, XTH30, PME3, EXPA11, MYB88, ROT3, AT5G37790, WAG2 and TIP2;3) that are down-regulated in mur3-3 partially rescued the short hypocotyl phenotype but not the aerial phenotype of mur3-3, indicating that different mechanisms exist between hypocotyl cell elongation and leaf cell elongation.

Comprehensive gene and microRNA expression profiling reveals a role for miRNAs in the oncogenic roles of SphK1 in papillary thyroid cancer.

The oncogenic roles of sphingosine kinase 1 (SphK1) in various cancers, including thyroid cancer, have been well demonstrated. However, the microRNAs (miRNAs) associated with the oncogenic roles of SphK1 remain largely unknown. Global gene and miRNA expression in TPC1-Vector and TPC1-SphK1 cells was analyzed using digital gene expression (DGE) analysis and small RNA-seq, respectively. miRNA-mRNA interactions were explored by microT-CDS, and the predicted networks were visualized using CytoScape®. Cell invasion and migration were assessed by performing Transwell invasion and wound-healing assays. Luciferase reporter and immunoblot assays were used to evaluate the targeting of fibronectin 1 (FN1) by miR-144-3p. In this study, we found that overexpression of SphK1 differentially regulates the expression of 46 miRNAs and 506 mRNAs in papillary thyroid cancer (PTC) TPC1 cells. Combining bioinformatics predictions of mRNA targets with DGE data on mRNA expression allowed us to identify the mRNA targets of deregulated miRNAs. The direct interaction between miR-144-3p and FN1, which mediates the pro-invasive role of SphK1 in PTC cells, was experimentally validated. Our results demonstrated that SphK1 overexpression drives a regulatory network governing miRNA and mRNA expression in PTC cells. We also demonstrated the roles played by miR-144-3p and FN1 in mediating the oncogenic function of SphK1, which enhanced the understanding of the etiology of PTC.

Transcriptional Profiling of Rice Treated with MoHrip1 Reveal the Function of Protein Elicitor in Enhancement of Disease Resistance and Plant Growth.

MoHrip1 is a protein elicitor isolated from Magnaporthe oryzae and was found to induce blast-resistance in rice. To investigate the comprehensive functions of MoHrip1, next-generation sequencing (NGS)-based digital gene expression (DGE) profiling was performed to collect the transcriptional data of differentially expressed genes (DEGs) induced by MoHrip1. A total of 308 genes were identified with differential expression, and 80 genes were predicted to be induced specifically by MoHrip1. Among these 308 genes, a series of genes associated with the salicylic acid (SA) pathway, phytoalexin, transcription factors, and pathogen-related proteins were identified. Both the SA signaling pathway and the gibberellin (GA) pathway were activated, while the jasmonic acid (JA) signaling pathway was repressed. The contents of endogenous SA and GA and the morphological characteristics of the rice after treatment were measured to provide evidence supporting the predictions made based on the DGE data. The 80 genes mentioned above might be candidate genes for studying interactions with MoHrip1. The transcriptional data provided global effect information in rice induced by MoHrip1, and all the results demonstrated that MoHrip1 could induce pathogen resistance and promote plant growth by regulating the contents of SA and GA directly or indirectly.

Targeted myocardial gene expression in failing hearts by RNA sequencing

BackgroundMyocardial recovery with left ventricular assist device (LVAD) therapy is highly variable and difficult to predict. Next generation ribonucleic acid (RNA) sequencing is an innovative, rapid, and quantitative approach to gene expression profiling in small amounts of tissue. Our primary goal was to identify baseline transcriptional profiles in non-ischemic cardiomyopathies that predict myocardial recovery in response to LVAD therapy. We also sought to verify transcriptional differences between failing and non-failing human hearts.MethodsRNA was isolated from failing (n = 16) and non-failing (n = 8) human hearts. RNA from each patient was reverse transcribed and quantitatively sequenced on the personal genome machine (PGM) sequencer (Ion torrent) for 95 heart failure candidate genes. Coverage analysis as well as mapping the reads and alignment was done using the Ion Torrent Browser Suite™. Differential expression analyses were conducted by empirical analysis of digital gene expression data in R (edgeR) to identify differential expressed genes between failing and non-failing groups, and between responder and non-responder groups respectively. Targeted cardiac gene messenger RNA (mRNA) expression was analyzed in proportion to the total number of reads. Gene expression profiles from the PGM sequencer were validated by performing RNA sequencing (RNAseq) with the Illumina Hiseq2500 sequencing system.ResultsThe failing sample population was 75% male with an average age of 50 and a left ventricular ejection fraction (LVEF) of 16%. Myosin light chain kinase (MYLK) and interleukin (IL)-6 genes expression were significantly higher in LVAD responders compared to non-responders. Thirty-six cardiac genes were expressed differentially between failing and non-failing hearts (23 decreased, 13 elevated). MYLK, Beta-1 adrenergic receptor (ADRB1) and myosin heavy chain (MYH)-6 expression were among those significantly decreased in failing hearts compared to non-failing hearts. Natriuretic peptide B (NPPB) and IL-6 were significantly elevated. Targeted gene expression profiles obtained from the Ion torrent PGM sequencer were consistent with those obtained from Illumina HiSeq2500 sequencing system.ConclusionsHeart failure is associated with a network of transcriptional changes involving contractile proteins, metabolism, adrenergic receptors, protein phosphorylation, and signaling factors. Myocardial MYLK and IL-6 expression are positively correlated with ejection fraction (EF) response to LVAD placement. Targeted RNA sequencing of myocardial gene expression can be utilized to predict responders to LVAD therapy and to better characterize transcriptional changes in human heart failure.

Comparative transcriptome profiling of chilling stress responsiveness in grafted watermelon seedlings

Rootstock grafting may improve the resistance of watermelon plants to low temperatures. However, information regarding the molecular responses of rootstock grafted plants to chilling stress is limited. To elucidate the molecular mechanisms of chilling tolerance in grafted plants, the transcriptomic responses of grafted watermelon under chilling stress were analyzed using RNA-seq analysis. Sequencing data were used for digital gene expression (DGE) analysis to characterize the transcriptomic responses in grafted watermelon seedlings. A total of 702 differentially-expressed genes (DEGs) were found in rootstock grafted (RG) watermelon relative to self-grafted (SG) watermelon; among these genes, 522 genes were up-regulated and 180 were down-regulated. Additionally, 164 and 953 genes were found to specifically expressed in RG and SG seedlings under chilling stress, respectively. Functional annotations revealed that up-regulated DEGs are involved in protein processing, plant-pathogen interaction and the spliceosome, whereas down-regulated DEGs are associated with photosynthesis. Moreover, 13 DEGs were randomly selected for quantitative real time PCR (qRT-PCR) analysis. The expression profiles of these 13 DEGs were consistent with those detected by the DGE analysis, supporting the reliability of the DGE data. This work provides additional insight into the molecular basis of grafted watermelon responses to chilling stress.

De novo transcriptome sequencing of black pepper (Piper nigrum L.) and an analysis of genes involved in phenylpropanoid metabolism in response to Phytophthora capsici.

Background Piper nigrum L., or “black pepper”, is an economically important spice crop in tropical regions. Black pepper production is markedly affected by foot rot disease caused by Phytophthora capsici, and genetic improvement of black pepper is essential for combating foot rot diseases. However, little is known about the mechanism of anti- P. capsici in black pepper. The molecular mechanisms underlying foot rot susceptibility were studied by comparing transcriptome analysis between resistant (Piper flaviflorum) and susceptible (Piper nigrum cv. Reyin-1) black pepper species.Results116,432 unigenes were acquired from six libraries (three replicates of resistant and susceptible black pepper samples), which were integrated by applying BLAST similarity searches and noted by adopting Kyoto Encyclopaedia of Genes and Gene Ontology (GO) genome orthology identifiers. The reference transcriptome was mapped using two sets of digital gene expression data. Using GO enrichment analysis for the differentially expressed genes, the majority of the genes associated with the phenylpropanoid biosynthesis pathway were identified in P. flaviflorum. In addition, the expression of genes revealed that after susceptible and resistant species were inoculated with P. capsici, the majority of genes incorporated in the phenylpropanoid metabolism pathway were up-regulated in both species. Among various treatments and organs, all the genes were up-regulated to a relatively high degree in resistant species. Phenylalanine ammonia lyase and peroxidase enzyme activity increased in susceptible and resistant species after inoculation with P. capsici, and the resistant species increased faster. The resistant plants retain their vascular structure in lignin revealed by histochemical analysis.ConclusionsOur data provide critical information regarding target genes and a technological basis for future studies of black pepper genetic improvements, including transgenic breeding.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3155-7) contains supplementary material, which is available to authorized users.

The role of melanin pathways in extremotolerance and virulence of Fonsecaea revealed by de novo assembly transcriptomics using illumina paired-end sequencing.

Melanisation has been considered to be an important virulence factor of Fonsecaea monophora. However, the biosynthetic mechanisms of melanisation remain unknown. We therefore used next generation sequencing technology to investigate the transcriptome and digital gene expression data, which are valuable resources to better understand the molecular and biological mechanisms regulating melanisation in F. monophora. We performed de novo transcriptome assembly and digital gene expression (DGE) profiling analyses of parent (CBS 122845) and albino (CBS 125194) strains using the Illumina RNA-seq system. A total of 17 352 annotated unigenes were found by BLAST search of NR, Swiss-Prot, Gene Ontology, Clusters of Orthologous Groups and Kyoto Encyclopedia of Genes and Genomes (KEGG) (E-value <1e‒5). A total of 2 283 unigenes were judged to be the differentially expressed between the two genotypes. We identified most of the genes coding for key enzymes involved in melanin biosynthesis pathways, including polyketide synthase (pks), multicopper oxidase (mco), laccase, tyrosinase and homogentisate 1,2-dioxygenase (hmgA). DEG analysis showed extensive down-regulation of key genes in the DHN pathway, while up-regulation was noted in the DOPA pathway of the albino mutant. The transcript levels of partial genes were confirmed by real time RT-PCR, while the crucial role of key enzymes was confirmed by either inhibitor or substrate tests in vitro. Meanwhile, numbers of genes involved in light sensing, cell wall synthesis, morphology and environmental stress were identified in the transcriptome of F. monophora. In addition, 3 353 SSRs (Simple Sequence Repeats) markers were identified from 21 600 consensus sequences. Blocking of the DNH pathway is the most likely reason of melanin deficiency in the albino strain, while the production of pheomelanin and pyomelanin were probably regulated by unknown transcription factors on upstream of both pathways. Most of genes involved in environmental tolerance to oxidants, irradiation and extreme temperatures were also assembled and annotated in transcriptomes of F. monophora. In addition, thousands of identified cSSR (combined SSR) markers will favour further genetic linkage studies. In conclusion, these data will contribute to understanding the regulation of melanin biosynthesis and help to improve the studies of pathogenicity of F. monophora.

In silico analysis of Mn transporters (NRAMP1) in various plant species.

Manganese (Mn) is an essential micronutrient in plant life cycle. It may be involved in photosynthesis, carbohydrate and lipid biosynthesis, and oxidative stress protection. Mn deficiency inhibits the plant growth and development, and causes the various plant symptoms such as interveinal chlorosis and tissue necrosis. Despite its importance in plant life cycle, we still have limited knowledge about Mn transporters in many plant species. Therefore, this study aimed to identify and characterize high affinity Arabidopsis Mn root transporter NRAMP1 orthologs in 17 different plant species. Various in silico methods and digital gene expression data were used in identification and characterization of NRAMP1 homologs; physico-chemical properties of sequences were calculated, putative transmembrane domains (TMDs) and conserved motif signatures were determined, phylogenetic tree was constructed, 3D models and interactome map were generated, and gene expression data was analyzed. 49 NRAMP1 homologs were identified from proteome datasets of 17 plant species using AtNRAMP1 as query. Identified sequences were characterized with a NRAMP domain structure, 10-12 putative TMDs with cytosolic N- and C-terminuses, and 10-14 exons encoding a protein of 500-588 amino acids and 53.8-64.3 kDa molecular weight with basic characteristics. Consensus transport residues, GQSSTITGTYAGQY(/F)V(/I)MQGFLD(/E/N) between TMD-8 and 9 were identified in all sequences but putative N-linked glycosylation sites were not highly conserved. In phylogeny, NRAMP1 sequences demonstrated divergence in lower and higher plants as well as in monocots and dicots. Despite divergence of lower plant Physcomitrella patens in phylogeny, it showed similarity in superposed 3D models. Phylogenetic distribution of AtNRAMP1 and 6 homologs inferred a functional relationship to NRAMP6 sequences in Mn transport, while distribution of OsNRAMP1 and 5 homologs implicated an involvement of NRAMP1 sequences in Mn transport or a cross-talk between in Fe-Mn homeostasis. Interactome analysis further confirmed this cross-talk between Mn and Fe pathways. Gene expression profile of AtNRAMP1 under Fe-, K-, P- and S-deficiencies, and cold, drought, heat and salt stresses revealed various proteins involving in transcription regulation, cofactor biosynthesis, diverse developmental roles, carbohydrate metabolism, oxidation-reduction reactions, cellular signaling and protein degradation pathways. Mn deficiency or toxicity could cause serious adverse effects in plants as well as in humans. To reduce these adversities mainly rely on understanding the molecular mechanisms underlying Mn uptake from the soil. However, we still have limited knowledge regarding the structural and functional roles of Mn transporters in many plant species. Therefore, identification and characterization of Mn root uptake transporter, NRAMP1 orthologs in various plant species will provide valuable theoretical knowledge to better understand Mn transporters as well as it may become an insight for future studies aiming to develop genetically engineered and biofortified plants.

Identification of genes regulated by histone acetylation during root development in Populus trichocarpa

BackgroundHistone deacetylases (HDACs) are key enzymes catalyzing the removal of acetyl groups from histones. HDACs act in concert with histone acetyltransferases (HATs) to regulate histone acetylation status, which modifies chromatin structure, affecting gene transcription and thus regulating multiple biological processes such as plant growth and development. Over a decade, certain HDACs in herbaceous plants have been deeply studied. However, functions of HDACs in woody plants are not well understood.ResultsHistone deacetylase specific inhibitor trichostatin A (TSA) was used to investigate the role of HDACs in organogenesis of roots and root development in Populus trochocarpa. The adventitious roots were regenerated and grown on medium supplemented with 0, 1, and 2.5 μM TSA. TSA treatment delayed root regeneration and inhibited primary root growth. To examine the genes modified by TSA in the regenerated roots, tag-based digital gene expression (DGE) analysis was performed using Illumina HiSeqTM 2000. Approximately 4.5 million total clean tags were mapped per library. The distinct clean tags for the three libraries corresponding to 0, 1 and 2.5 μM TSA treatment were 166167, 143103 and 153507, from which 38.45 %, 31.84 % and 38.88 % were mapped unambiguously to the unigene database, respectively. Most of the tags were expressed at similar levels, showing a < 5-fold difference after 1 μM and 2.5 μM TSA treatments and the maximum fold-change of the tag copy number was around 20. The expression levels of many genes in roots were significantly altered by TSA. A total of 36 genes were up-regulated and 1368 genes were down-regulated after 1 μM TSA treatment, while 166 genes were up-regulated and 397 genes were down-regulated after 2.5 μM TSA treatment. Gene ontology (GO) and pathway analyses indicated that the differentially expressed genes were related to many kinds of molecular functions and biological processes. The genes encoding key enzymes catalyzing gibberellin biosynthesis were significantly down-regulated in the roots exposed to 2.5 μM TSA and their expression changes were validated by using real-time PCR.ConclusionsHDACs were required for de novo organogenesis and normal growth of populus roots. DGE data provides the gene profiles in roots probably regulated by histone acetylation during root growth and development, which will lead to a better understanding of the mechanism controlling root development.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2407-x) contains supplementary material, which is available to authorized users.

Abstract B1-40: Biological network analysis using an in vitro model of androgen-resistance in prostate cancer

Abstract The development of androgen resistance is a major limitation to androgen deprivation treatment in prostate cancer (PCa). We have developed an in vitro model of androgen-resistance using the androgen sensitive cell line LNCaP to characterize the phenotypic and transcriptomic changes occurring as androgen resistance develops. Our aim is to understand biological network profiles of transcriptomic changes during the transition to androgen-resistance and to validate these changes between our in vitro model and previously published clinical datasets (paired samples before and after androgen-deprivation therapy of patients with advanced PCa) (1). Methods: PCa cells, LNCaP, expressing mutated AR which are androgen-dependent (2,3), are used in the development of an androgen-resistant subline. Subline cells were established by prolonged cultures in media + 10% CS-FBS to mimic the clinical course of PCa. Cell proliferation, cell motility and invasion, morphology, AR expression were examined. RNA-sequencing was performed using the parental LNCaP cells and an androgen-resistant subline (LNCaP-AI) established by chronic exposure to the androgen-deprivation. Reads from cells and clinical samples (1), pre- vs. post-treatment, were processed through the same standard pipeline and quality control. Data outputs were analysed as differential expression [DEG] (EdgeR) (4) and top scoring protein-protein interaction (PPI) networks [PINA2 (5) and BioNet (6)]. Data outputs from the cell line model and clinical samples were compared. Results: LNCaP cells initially showed poor growth after prolonged exposure to androgen-deprived conditions but later adapted and started to grow well. After 24 weeks of androgen-deprivation, LNCaP-AI's growth was no longer responsive to addition of androgen [0.1 - 10 nM]. AR expression was not different in LNCaP and LNCaP-AI (P&amp;gt;0.05). LNCaP-AI cells had increased proliferation and cell invasion compared to LNCaP. We identified key genes that overlap between our cell line and clinical RNAseq (1) datasets and analyzed the overlapping PPI network that showed the same pattern of behavior in both datasets. The network revealed several potential mechanisms and gene interactions that warrant further investigation, including cooperative behaviors of other nuclear receptors, TP63 mediated signalling pathway and Aryl hydrocarbon receptor transcriptional pathway. Conclusion: Cell line model of androgen-resistance will be used for further longitudinal study of the mechanism of castrate resistant prostate cancer (CRPC). Our approach allows for better characterization of biological processes of CRPC. Knowledge of the genetic profiles during transition to androgen resistance will improve our understanding of this common clinical scenario and may lead to biomarker discovery.

Using digital gene expression profile to detect representational difference of Chlamys farreri genes after laboratory exposure to persistent organic pollutants

Chlamys farreri (C. farreri), as an economic species of marine bivalves, can tolerate manifold anthropogenic stressors including persistent organic pollutants (POPs) stress. A systematic study of transcriptome after POPs exposure may provide insights into the mechanism of acquired pollution tolerance. In this study, we compared C. farreri transcriptome by reanalyzing DGE (digital gene expression) data from previous study. Our results revealed a toxicant-dependent pattern of global transcriptional responses, with 108, 126 and 138 DEGs regulated by BaP, TBBPA and PAHs cocktail (mixed by benzo(a)pyrene, benz(a)anthracene, benzo(b)fluoranthene and chrysene with the ratio of 5:3:1:1) exposure, respectively. Among these DEGs, 219 were commonly regulated by all pollutants, whereas the minority of differences were condition-specific. In addition, we performed the first hierarchical cluster analysis of the common transcripts throughout three types of contaminants stress ulteriorly, leading to new discoveries of genes’ association information with dynamical gene expression data of C. farreri, which goes beyond the focus on individuals. It revealed a cascade of gene expression patterns in the response of C. farreri to different environmental stresses and may be beneficial for further analysis of novel environmental metrics.

Activation of the TOR Signalling Pathway by Glutamine Regulates Insect Fecundity.

The target of rapamycin (TOR) positively controls cell growth in response to nutrients such as amino acids. However, research on the specific nutrients sensed by TOR is limited. Glutamine (Gln), a particularly important amino acid involved in metabolism in organisms, is synthesised and catalysed exclusively by glutamine synthetase (GS), and our previous studies have shown that Gln may regulate fecundity in vivo levels of the brown planthopper (BPH) Nilaparvata lugens. Until now, it has remained unclear whether Gln activates or inhibits the TOR signalling pathway. Here, we performed the combined analyses of iTRAQ (isobaric tags for relative and absolute quantification) and DGE (tag-based digital gene expression) data in N.lugens at the protein and transcript levels after GS RNAi, and we found that 52 pathways overlap, including the TOR pathway. We further experimentally demonstrate that Gln activates the TOR pathway by promoting the serine/threonine protein kinase AKT and inhibiting the 5′AMP-activated protein kinase AMPK phosphorylation activity in the pest. Furthermore, TOR regulates the fecundity of N. lugens probably by mediating vitellogenin (Vg) expression. This work is the first report that Gln activates the TOR pathway in vivo.

Transcriptome-wide identification of salt-responsive members of the WRKY gene family in Gossypium aridum.

WRKY transcription factors are plant-specific, zinc finger-type transcription factors. The WRKY superfamily is involved in abiotic stress responses in many crops including cotton, a major fiber crop that is widely cultivated and consumed throughout the world. Salinity is an important abiotic stress that results in considerable yield losses. In this study, we identified 109 WRKY genes (GarWRKYs) in a salt-tolerant wild cotton species Gossypium aridum from transcriptome sequencing data to elucidate the roles of these factors in cotton salt tolerance. According to their structural features, the predicted members were divided into three groups (Groups I–III), as previously described for Arabidopsis. Furthermore, 28 salt-responsive GarWRKY genes were identified from digital gene expression data and subjected to real-time quantitative RT-PCR analysis. The expression patterns of most GarWRKY genes revealed by this analysis are in good agreement with those revealed by RNA-Seq analysis. RT-PCR analysis revealed that 27 GarWRKY genes were expressed in roots and one was exclusively expressed in roots. Analysis of gene orthology and motif compositions indicated that WRKY members from Arabidopsis, rice and soybean generally shared the similar motifs within the same subgroup, suggesting they have the similar function. Overexpression-GarWRKY17 and –GarWRKY104 in Arabidopsis revealed that they could positively regulate salt tolerance of transgenic Arabidopsis during different development stages. The comprehensive data generated in this study provide a platform for elucidating the functions of WRKY transcription factors in salt tolerance of G. aridum. In addition, GarWRKYs related to salt tolerance identified in this study will be potential candidates for genetic improvement of cultivated cotton salt stress tolerance.

MiRNAome analysis associated with anatomic and transcriptomic investigations reveal the polar exhibition of corky split vein in boron deficient Citrus sinensis.

Corky split vein can develop under long-term boron deficient conditions in Citrus sinensis L. Osbeck cv. Newhall. This symptom only occurs in the upper rather than the lower epidermis of old leaves. Our previous study demonstrated that vascular hypertrophy was involved in the symptoms, and the 3rd developmental stage of corky split vein (BD3) was the critical stage for phenotype formation. Here, we performed an intensive study on the BD3 vein and its control sample (CK3 vein). A lignin test demonstrated that the lignin content in BD3 vein was approximately 1.7 times more than the CK3 vein. Anatomical investigation of the corky split vein indicated that the upper epidermis was destroyed by overgrowing vascular cells, and the increased lignin may contribute to vascular cell differentiation and wounding-induced lignification. In a subsequent small RNA sequencing of the BD3 and CK3 veins, 99 known miRNAs and 22 novel miRNAs were identified. Comparative profiling of these miRNAs demonstrated that the 57 known miRNAs and all novel miRNAs exhibited significant expression differences between the two small RNAs libraries of the BD3 and CK3 veins. Associated with our corresponding digital gene expression data, we propose that the decreased expression of two miRNAs, csi-miR156b and csi-miR164, which leads to the up-regulation of their target genes, SPLs (csi-miR156b-targeted) and CUC2 (csi-miR164-targeted), may promote vascular cell division and orderless stage transition in old leaves.

Expression profiling of cadmium response genes in ramie (Boehmeria nivea L.) root.

Ramie (Boehmeria nivea), a perennial herb belongs to Urticaceae family, is a rapid growth and high biomass crop with highly tolerant and accumulative to heavy metals. However, the gene expression and regulation caused by cadmium (Cd) in ramie has not been well studied. In the present study, a gene expression database of ramie root in the absence (control) or presence of 100 μM Cd was established. Solexa high-throughput sequencing technology showed that 3,654,395 and 3,572,333 tags have been obtained from control and Cd treatment respectively. In total, 3887 genes were detected with significant differential expression levels, in which 2883 genes were up-regulated and 1004 genes were down-regulated. Gene ontology and pathway-based analyses were performed to determine and further to understand the biological functions of those differentially expressed genes. Fifteen genes were selected and their expression levels were confirmed by quantitative RT-PCR, and twelve of them showed consistent expression patterns with the digital gene expression data. Results on these expression profiling of genes lay the basis for biotechnological modification of new transgenic plants with improved phytoremediation capacity.

Genome-wide analysis of the GRAS gene family in physic nut (Jatropha curcas L.).

GRAS proteins play vital roles in plant growth and development. Physic nut (Jatropha curcas L.) was found to have a total of 48 GRAS family members (JcGRAS), 15 more than those found in Arabidopsis. The JcGRAS genes were divided into 12 subfamilies or 15 ancient monophyletic lineages based on the phylogenetic analysis of GRAS proteins from both flowering and lower plants. The functions of GRAS genes in 9 subfamilies have been reported previously for several plants, while the genes in the remaining 3 subfamilies were of unknown function; we named the latter families U1 to U3. No member of U3 subfamily is present in Arabidopsis and Poaceae species according to public genome sequence data. In comparison with the number of GRAS genes in Arabidopsis, more were detected in physic nut, resulting from the retention of many ancient GRAS subfamilies and the formation of tandem repeats during evolution. No evidence of recent duplication among JcGRAS genes was observed in physic nut. Based on digital gene expression data, 21 of the 48 genes exhibited differential expression in four tissues analyzed. Two members of subfamily U3 were expressed only in buds and flowers, implying that they may play specific roles. Our results provide valuable resources for future studies on the functions of GRAS proteins in physic nut.