A comparison of the Thlaspi caerulescens and Thlaspi arvense shoot transcriptomes

Whole-genome transcriptome profiling is revealing how biological systems are regulated at the transcriptional level. This study reports the development of a robust method to profile and compare the transcriptomes of two nonmodel plant species, Thlaspi caerulescens, a zinc (Zn) hyperaccumulator, and Thlaspi arvense, a nonhyperaccumulator, using Affymetrix Arabidopsis thaliana ATH1-121501 GeneChip (R) arrays (Affymetrix, Santa Clara, CA, USA). Transcript abundance was quantified in the shoots of agar- and compost-grown plants of both species. Analyses were optimized using a genomic DNA (gDNA)-based probe-selection strategy based on the hybridization efficiency of Thlaspi gDNA with corresponding A. thaliana probes. In silico alignments of GeneChip (R) probes with Thlaspi gene sequences, and quantitative real-time PCR, confirmed the validity of this approach. Approximately 5000 genes were differentially expressed in the shoots of T. caerulescens compared with T. arvense, including genes involved in Zn transport and compartmentalization. Future functional analyses of genes identified as differentially expressed in the shoots of these closely related species will improve our understanding of the molecular mechanisms of Zn hyperaccumulation.

Comprehensive genome methylation and whole genome expression analysis in penile carcinoma: Uncovering new molecular markers

Background: Penile carcinoma (PeCa) is frequently associated with high morbidity rates. Unlikely of the vast majority of tumors, there is no molecular markers described that are able to assist in diagnosis and prognosis or with potential to be therapeutic targets in PeCa. Patients and methods: DNA methylation status (244K Human DNA Methylation Microarray platform, Agilent Technologies) and large-scale expression analysis (4x44K Whole Human Genome Microarray, Agilent Technologies) were performed in 35 and 37 PeCa, respectively. Quantitative bisulfite pyrosequencing (qBP) and RT-qPCR were used to validate the findings in 93 samples. HPV status was assessed using the Linear Array HPV Genotyping kit (Roche Molecular Diagnostics, CA, USA). Results: Methylome analysis revealed 171 hypermethylated and 449 hypomethylated CpGs sites and the transcriptome profiling showed 2986 down- and 2817 over-expressed genes. HPV positivity was found in 32.7% of the cases, mainly the HPV16. The integrative analysis in 32 PeCa revealed a panel of 96 genes with inverse correlation between methylation and gene expression levels. The CpG hypermetlylation and gene downexpression, was confirmed for TWIST1, RSOP2, SOX3, SOX17, CD133, OTX2, HOXA3 and MEIS. In addition, BIRC5, DNMT1 and DNMT3B presented low levels of methylation and overexpression. The comparison of the results with clinical findings revealed that LIN28A, NKX2.2, NKX2.3, LHX5, BDNF, FOXA1 and CDX2 were associated with poor prognosis features. Conclusion: Putative prognostic markers were detected revealing that DNA methylation modulates the expression of several genes in PeCa. These data may prove instrumental for biomarker discovery in clinics and molecular epidemiology of PeCa.

Overexpression of mitochondrial uncoupling protein 1 (UCP1) induces a hypoxic response in Nicotiana tabacum leaves

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The Role of Acidic Fibroblast Growth Factor and Fibroblast Growth Factor Receptor 4 in High Grade Serous Carcinoma

Background: High grade serous carcinoma whether ovarian, tubal or primary peritoneal, continues to be the most lethal gynecologic malignancy in the USA. Although combination chemotherapy and aggressive surgical resection has improved survival in the past decade the majority of patients still succumb to chemo-resistant disease recurrence. It has recently been reported that amplification of 5q31-5q35.3 is associated with poor prognosis in patients with high grade serous ovarian carcinoma. Although the amplicon contains over 50 genes, it is notable for the presence of several members of the fibroblast growth factor signaling axis. In particular acidic fibroblast growth factor (FGF1) has been demonstrated to be one of the driving genes in mediating the observed prognostic effect of the amplicon in ovarian cancer patients. This study seeks to further validate the prognostic value of fibroblast growth receptor 4 (FGFR4), another candidate gene of the FGF/FGFR axis located in the same amplicon. The emphasis will be delineating the role the FGF1/FGFR4 signaling axis plays in high grade serous ovarian carcinoma; and test the feasibility of targeting the FGF1/FGFR4 axis therapeutically. Materials and Methods: Spearman and Pearson correlation studies on data generated from array CGH and transcriptome profiling analyses on 51 microdissected tumor samples were used to identify genes located on chromosome 5q31-35.3 that showed significant correlation between DNA and mRNA copy numbers. Significant correlation between FGF1 and FGFR4 DNA copy numbers was further validated by qPCR analysis on DNA isolated from 51 microdissected tumor samples. Immunolocalization and quantification of FGFR4 expression were performed on paraffin embedded tissue samples from 183 cases of high-grade serous ovarian carcinoma. The expression was then correlated with clinical data to assess impact on survival. The expression of FGF1 and FGFR4 in vitro was quantified by real-time PCR and western blotting in six high-grade serous ovarian carcinoma cell lines and compared to those in human ovarian surface epithelial cells to identify overexpression. The effect of FGF1 on these cell lines after serum starvation was quantified for in vitro cellular proliferation, migration/invasion, chemoresistance and survival utilizing a combination of commercially available colorimetric, fluorometric and electrical impedance assays. FGFR4 expression was then transiently silenced via siRNA transfection and the effects on response to FGF1, cellular proliferation, and migration were quantified. To identify relevant cellular pathways involved, responsive cell lines were transduced with different transcription response elements using the Cignal-Lenti reporter system and treated with FGF1 with and without transient FGFR4 knock down. This was followed by western blot confirmation for the relevant phosphoproteins. Anti-FGF1 antibodies and FGFR trap proteins were used to attempt inhibition of FGF mediated phenotypic changes and relevant signaling in vitro. Orthotopic intraperitoneal tumors were established in nude mice using serous cell lines that have been previously transfected with luciferase expressing constructs. The mice were then treated with FGFR trap protein. Tumor progression was then followed via bioluminescent imaging. The FGFR4 gene from 52 clinical samples was sequenced to screen for mutations. Results: FGFR4 DNA and mRNA copy numbers were significantly correlated and FGFR4 DNA copy number was significantly correlated with that of FGF1. Survival of patients with high FGFR4 expressing tumors was significantly shorter that those with low expression(median survival 28 vs 55 month p< 0.001) In a multivariate cox regression model FGFR expression significantly increased risk of death (HR 2.1, p<0.001). FGFR4 expression was significantly higher in all cell lines tested compared to HOSE, OVCA432 cell line in particular had very high expression suggesting amplification. FGF1 was also particularly overexpressed in OVCA432. FGF1 significantly increased cell survival after serum deprivation in all cell lines. Transient knock down of FGFR4 caused significant reduction in cell migration and proliferation in vitro and significantly decreased the proliferative effects of FGF1 in vitro. FGFR1, FGFR4 traps and anti-FGF1 antibodies did not show activity in vitro. OVCA432 transfected with the cignal lenti reporter system revealed significant activation of MAPK, NFkB and WNT pathways, western blotting confirmed the results. Reverse phase protein array (RPPA) analysis also showed activation of MAPK, AKT, WNT pathways and down regulation of E Cadherin. FGFR trap protein significantly reduced tumor growth in vivo in an orthotopic mouse model. Conclusions: Overexpression and amplification of several members of the FGF signaling axis present on the amplicon 5q31-35.3 is a negative prognostic indicator in high grade serous ovarian carcinoma and may drive poor survival associated with that amplicon. Activation of The FGF signaling pathway leads to downstream activation of MAPK, AKT, WNT and NFkB pathways leading to a more aggressive cancer phenotype with increased tumor growth, evasion of apoptosis and increased migration and invasion. Inhibition of FGF pathway in vivo via FGFR trap protein leads to significantly decreased tumor growth in an orthotopic mouse model.

Regulation of human (adrenal) androgen biosynthesis-New insights from novel throughput technology studies

Androgens are precursors for sex steroids and are predominantly produced in the human gonads and the adrenal cortex. They are important for intrauterine and postnatal sexual development and human reproduction. Although human androgen biosynthesis has been extensively studied in the past, exact mechanisms underlying the regulation of androgen production in health and disease remain vague. Here, the knowledge on human androgen biosynthesis and regulation is reviewed with a special focus on human adrenal androgen production and the hyperandrogenic disorder of polycystic ovary syndrome (PCOS). Since human androgen regulation is highly specific without a good animal model, most studies are performed on patients harboring inborn errors of androgen biosynthesis, on human biomaterials and human (tumor) cell models. In the past, most studies used a candidate gene approach while newer studies use high throughput technologies to identify novel regulators of androgen biosynthesis. Using genome wide association studies on cohorts of patients, novel PCOS candidate genes have been recently described. Variant 2 of the DENND1A gene was found overexpressed in PCOS theca cells and confirmed to enhance androgen production. Transcriptome profiling of dissected adrenal zones established a role for BMP4 in androgen synthesis. Similarly, transcriptome analysis of human adrenal NCI-H295 cells identified novel regulators of androgen production. Kinase p38α (MAPK14) was found to phosphorylate CYP17 for enhanced 17,20 lyase activity and RARB and ANGPTL1 were detected in novel networks regulating androgens. The discovery of novel players for androgen biosynthesis is of clinical significance as it provides targets for diagnostic and therapeutic use.

Ocean Acidification Affects Redox-Balance and Ion-Homeostasis in the Life-Cycle Stages of Emiliania huxleyi

Ocean Acidification (OA) has been shown to affect photosynthesis and calcification in the coccolithophore Emiliania huxleyi, a cosmopolitan calcifier that significantly contributes to the regulation of the biological carbon pumps. Its non-calcifying, haploid life-cycle stage was found to be relatively unaffected by OA with respect to biomass production. Deeper insights into physiological key processes and their dependence on environmental factors are lacking, but are required to understand and possibly estimate the dynamics of carbon cycling in present and future oceans. Therefore, calcifying diploid and non-calcifying haploid cells were acclimated to present and future CO2 partial pressures (pCO2; 38.5 Pa vs. 101.3 Pa CO2) under low and high light (50 vs. 300 µmol photons/m**2 /s). Comparative microarray-based transcriptome profiling was used to screen for the underlying cellular processes and allowed to follow up interpretations derived from physiological data. In the diplont, the observed increases in biomass production under OA are likely caused by stimulated production of glycoconjugates and lipids. The observed lowered calcification under OA can be attributed to impaired signal-transduction and ion-transport. The haplont utilizes distinct genes and metabolic pathways, reflecting the stage-specific usage of certain portions of the genome. With respect to functionality and energy-dependence, however, the transcriptomic OA-responses resemble those of the diplont. In both life-cycle stages, OA affects the cellular redox-state as a master regulator and thereby causes a metabolic shift from oxidative towards reductive pathways, which involves a reconstellation of carbon flux networks within and across compartments. Whereas signal transduction and ion-homeostasis appear equally OA-sensitive under both light intensities, the effects on carbon metabolism and light physiology are clearly modulated by light availability. These interactive effects can be attributed to the influence of OA and light on the redox equilibria of NAD and NADP, which function as major sensors for energization and stress. This generic mode of action of OA may therefore provoke similar cell-physiological responses in other protists.

Opening Basement Membrane Gaps During C. elegans Uterine-Vulval Attachment

The basement membrane (BM) is a highly conserved form of extracellular matrix that underlies or surrounds and supports most animal tissues. BMs are crossed by cells during various remodeling events in development, immune surveillance, or during cancer metastasis. Because BMs are dense and not easily penetrable, most of these cells must open a gap in order to facilitate their migration. The mechanisms by which cells execute these changes are poorly understood. A developmental event that requires the opening of a BM gap is C. elegans uterine-vulval connection. The anchor cell (AC), a specialized uterine cell, creates a de novo BM gap. Subsequent widening of the BM gap involves the underlying vulval precursor cells (VPCs) and the π cells, uterine neighbors of the AC through non-proteolytic BM sliding. Using forward and reverse genetic screening, transcriptome profiling, and live-cell imaging, I investigated how the cells in these tissues accomplish BM gap formation. In Chapter 2, I identify two potentially novel regulators of BM breaching, isolated through a large-scale forward genetic screen and characterize the invasion defect in these mutants. In Chapter 3, I describe single-cell transcriptome sequencing of the invasive AC. In Chapter 4, I describe the role of the π cells in opening the nascent BM gap. A complete developmental pathway for this process has been elucidated: the AC induces the π fate through Notch signaling, after which the π cells upregulate the Sec14 family protein CTG-1, which in turn restricts the trafficking of DGN-1 (dystroglycan), a laminin receptor, allowing the BM to slide. Chapter 5 outlines the implications of these discoveries.

Genome-wide profiling of the cardiac transcriptome after myocardial infarction identifies novel heart-specific long non-coding RNAs.

AIM: Heart disease is recognized as a consequence of dysregulation of cardiac gene regulatory networks. Previously, unappreciated components of such networks are the long non-coding RNAs (lncRNAs). Their roles in the heart remain to be elucidated. Thus, this study aimed to systematically characterize the cardiac long non-coding transcriptome post-myocardial infarction and to elucidate their potential roles in cardiac homoeostasis. METHODS AND RESULTS: We annotated the mouse transcriptome after myocardial infarction via RNA sequencing and ab initio transcript reconstruction, and integrated genome-wide approaches to associate specific lncRNAs with developmental processes and physiological parameters. Expression of specific lncRNAs strongly correlated with defined parameters of cardiac dimensions and function. Using chromatin maps to infer lncRNA function, we identified many with potential roles in cardiogenesis and pathological remodelling. The vast majority was associated with active cardiac-specific enhancers. Importantly, oligonucleotide-mediated knockdown implicated novel lncRNAs in controlling expression of key regulatory proteins involved in cardiogenesis. Finally, we identified hundreds of human orthologues and demonstrate that particular candidates were differentially modulated in human heart disease. CONCLUSION: These findings reveal hundreds of novel heart-specific lncRNAs with unique regulatory and functional characteristics relevant to maladaptive remodelling, cardiac function and possibly cardiac regeneration. This new class of molecules represents potential therapeutic targets for cardiac disease. Furthermore, their exquisite correlation with cardiac physiology renders them attractive candidate biomarkers to be used in the clinic.

Systematic expression profiling of the mouse transcriptome using RIKEN cDNA mircoarrays

The number of known mRNA transcripts in the mouse has been greatly expanded by the RIKEN Mouse Gene Encyclopedia project. Validation of their reproducible expression in a tissue is an important contribution to the study of functional genomics. In this report, we determine the expression profile of 57,931 clones on 20 mouse tissues using cDNA microarrays. Of these 57,931 clones, 22,928 clones correspond to the FANTOM2 clone set. The set represents 20,234 transcriptional units (TUs) out of 33,409 TUs in the FANTOM2 set. We identified 7206 separate clones that satisfied stringent criteria for tissue-specific expression. Gene Ontology terms were assigned for these 7206 clones, and the proportion of 'molecular function' ontology for each tissue-specific clone was examined. These data will provide insights into the function of each tissue. Tissue-specific gene expression profiles obtained using our cDNA microarrays were also compared with the data extracted from the GNF Expression Atlas based on Affymetrix microarrays. One major outcome of the RIKEN transcriptome analysis is the identification of numerous nonprotein-coding mRNAs. The expression profile was also used to obtain evidence of expression for putative noncoding RNAs. In addition, 1926 clones (70%) of 2768 clones that were categorized as unknown EST, and 1969 (58%) clones of 3388 clones that were categorized as unclassifiable were also shown to be reproducibly expressed.

Retrotransposon expression profiling: Unveiling the hidden SINE transcriptome through Next-Generation Sequencing data analysis

Of the ~1.7 million SINE elements in the human genome, only a tiny number are estimated to be active in transcription by RNA polymerase (Pol) III. Tracing the individual loci from which SINE transcripts originate is complicated by their highly repetitive nature. By exploiting RNA-Seq datasets and unique SINE DNA sequences, we devised a bioinformatic pipeline allowing us to identify Pol III-dependent transcripts of individual SINE elements. When applied to ENCODE transcriptomes of seven human cell lines, this search strategy identified ~1300 Alu loci and ~1100 MIR loci corresponding to detectable transcripts, with ~120 and ~60 respectively Alu and MIR loci expressed in at least three cell lines. In vitro transcription of selected SINEs did not reflect their in vivo expression properties, and required the native 5’-flanking region in addition to internal promoter. We also identified a cluster of expressed AluYa5-derived transcription units, juxtaposed to snaR genes on chromosome 19, formed by a promoter-containing left monomer fused to an Alu-unrelated downstream moiety. Autonomous Pol III transcription was also revealed for SINEs nested within Pol II-transcribed genes raising the possibility of an underlying mechanism for Pol II gene regulation by SINE transcriptional units. Moreover the application of our bioinformatic pipeline to both RNA-seq data of cells subjected to an in vitro pro-oncogenic stimulus and of in vivo matched tumor and non-tumor samples allowed us to detect increased Alu RNA expression as well as the source loci of such deregulation. The ability to investigate SINE transcriptomes at single-locus resolution will facilitate both the identification of novel biologically relevant SINE RNAs and the assessment of SINE expression alteration under pathological conditions.

High-resolution Transcript Profiling Of The Atypical Biotrophic Interaction Between Theobroma Cacao And The Fungal Pathogen Moniliophthora Perniciosa.

Witches' broom disease (WBD), caused by the hemibiotrophic fungus Moniliophthora perniciosa, is one of the most devastating diseases of Theobroma cacao, the chocolate tree. In contrast to other hemibiotrophic interactions, the WBD biotrophic stage lasts for months and is responsible for the most distinctive symptoms of the disease, which comprise drastic morphological changes in the infected shoots. Here, we used the dual RNA-seq approach to simultaneously assess the transcriptomes of cacao and M. perniciosa during their peculiar biotrophic interaction. Infection with M. perniciosa triggers massive metabolic reprogramming in the diseased tissues. Although apparently vigorous, the infected shoots are energetically expensive structures characterized by the induction of ineffective defense responses and by a clear carbon deprivation signature. Remarkably, the infection culminates in the establishment of a senescence process in the host, which signals the end of the WBD biotrophic stage. We analyzed the pathogen's transcriptome in unprecedented detail and thereby characterized the fungal nutritional and infection strategies during WBD and identified putative virulence effectors. Interestingly, M. perniciosa biotrophic mycelia develop as long-term parasites that orchestrate changes in plant metabolism to increase the availability of soluble nutrients before plant death. Collectively, our results provide unique insight into an intriguing tropical disease and advance our understanding of the development of (hemi)biotrophic plant-pathogen interactions.

A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Expression profiling of rectal tumors defines response to neoadjuvant treatment related genes.

To date, no effective method exists that predicts the response to preoperative chemoradiation (CRT) in locally advanced rectal cancer (LARC). Nevertheless, identification of patients who have a higher likelihood of responding to preoperative CRT could be crucial in decreasing treatment morbidity and avoiding expensive and time-consuming treatments. The aim of this study was to identify signatures or molecular markers related to response to pre-operative CRT in LARC. We analyzed the gene expression profiles of 26 pre-treatment biopsies of LARC (10 responders and 16 non-responders) without metastasis using Human WG CodeLink microarray platform. Two hundred and fifty seven genes were differentially over-expressed in the responder patient subgroup. Ingenuity Pathway Analysis revealed a significant ratio of differentially expressed genes related to cancer, cellular growth and proliferation pathways, and c-Myc network. We demonstrated that high Gng4, c-Myc, Pola1, and Rrm1 mRNA expression levels was a significant prognostic factor for response to treatment in LARC patients (p<0.05). Using this gene set, we were able to establish a new model for predicting the response to CRT in rectal cancer with a sensitivity of 60% and 100% specificity. Our results reflect the value of gene expression profiling to gain insight about the molecular pathways involved in the response to treatment of LARC patients. These findings could be clinically relevant and support the use of mRNA levels when aiming to identify patients who respond to CRT therapy.

The transcriptome of the arbuscular mycorrhizal fungus Glomus intraradices (DAOM 197198) reveals functional tradeoffs in an obligate symbiont.

? The arbuscular mycorrhizal symbiosis is arguably the most ecologically important eukaryotic symbiosis, yet it is poorly understood at the molecular level. To provide novel insights into the molecular basis of symbiosis-associated traits, we report the first genome-wide analysis of the transcriptome from Glomus intraradices DAOM 197198. ? We generated a set of 25,906 nonredundant virtual transcripts (NRVTs) transcribed in germinated spores, extraradical mycelium and symbiotic roots using Sanger and 454 sequencing. NRVTs were used to construct an oligoarray for investigating gene expression. ? We identified transcripts coding for the meiotic recombination machinery, as well as meiosis-specific proteins, suggesting that the lack of a known sexual cycle in G. intraradices is not a result of major deletions of genes essential for sexual reproduction and meiosis. Induced expression of genes encoding membrane transporters and small secreted proteins in intraradical mycelium, together with the lack of expression of hydrolytic enzymes acting on plant cell wall polysaccharides, are all features of G. intraradices that are shared with ectomycorrhizal symbionts and obligate biotrophic pathogens. ? Our results illuminate the genetic basis of symbiosis-related traits of the most ancient lineage of plant biotrophs, advancing future research on these agriculturally and ecologically important symbionts.

Hypothalamus transcriptome profile suggests an anorexia-cachexia syndrome in the anx/anx mouse model

The anx/anx mouse displays poor appetite and lean appearance and is considered a good model for the study of anorexia nervosa. To identify new genes involved in feeding behavior and body weight regulation we performed an expression profiling in the hypothalamus of the anx/anx mice. Using commercial microarrays we detected 156 differentially expressed genes and validated 92 of those using TaqMan low-density arrays. The expression of a set of 87 candidate genes selected based on literature evidences was also quantified by TaqMan low-density arrays. Our results showed enrichment in deregulated genes involved in cell death, cell morphology and cancer as well as an alteration of several signaling circuits involved in energy balance including neuropeptide Y and melanocortin signaling. The expression profile along with the phenotype led us to conclude that anx/anx mice resemble the anorexia-cachexia syndrome typically observed in cancer, infection with human immunodeficiency virus or chronic diseases, rather than starvation, and that anx/anx mice could be considered a good model for the treatment and investigation of this condition.