Gene expression under thermal stress varies across a geographical range expansion front

Acknowledgements: We would like to thank Hanna Bensch and Hannes Weise for assistance with the collection of samples in the field. This work was supported by the Biodiversity and Ecosystem Services in a Changing Climate (BECC; a joint Lund-Gothenburg University initiative), the Swedish Research Council (EIS, BH), the Crafoord Foundation (EIS, BH), the Swedish Royal Society (EIS), ‘Gyllenstiernska Krapperupstiftelsen (EIS), the Wenner-Gren Foundations (postdoctoral stipend to RYD), EU FP7 (Marie Curie International Incoming Fellowship to RYD), the Kungliga Fysiografiska Sällskapet i Lund (MW) and the Helge Ax:son Johnson Stiftelse (MW). B.H. and E.I.S. conceived of the study. L.L. developed the hypotheses to be tested. L.L. and R.D. collected the field data and samples. All six authors contributed to planning RNA-seq analyses. P.C. and L.L. analysed the data. L.L. wrote the manuscript, which all six authors edited.

Single-Cell Transcriptome Analysis of Olfactory Sensory Neurons

Olfactory sensory neurons (OSNs), which detect a myriad of odorants, are known to express one allele of one olfactory receptor (OR) gene (Olfr) from the largest gene family in the mammalian genome. The OSNs expressing the same OR project their axons to the main olfactory bulb where they converge to form glomeruli. This “One neuron-one receptor rule” makes the olfactory epithelium (OE), which consists of a vast number of OSNs expressing unique ORs, one of the most heterogeneous cell populations. However, the mechanism of how the single OR allele is chosen remains unclear along with the question of whether one OSN only expresses a single OR gene, a hypothesis that has not been rigorously verified while we performed the experiments. Moreover, failure of axonal targeting to single glomerulus was observed in MeCP2 deficient OSNs where delayed development was proposed as an explanation for the phenotype. How Mecp2 mutation caused this aberrant targeting is not entirely understood.

In this dissertation, we explored the transcriptomes of single and mature OSNs by single-cell RNA-Seq to reveal their heterogeneity and further studied the OR gene expression from these isolated OSNs. The singularity of sequenced OSNs was ensured by the observation of monoallelic expression of X-linked genes from the hybrid samples from crosses between mice of different strains where strain-specific polymorphisms could be used to track the allelic origins of SNP-containing reads. The clustering of expression profiles from triplicates that originated from the same cell assured that the transcriptomic identities of OSNs were maintained through the experimental process. The average gene expression profiles of sequenced OSNs correlated well to the conventional transcriptome data of FACS-sorted Omp-positive cells, and the top-ranked expression of OR was conceded in the single-OSN transcriptomes. While exploring cellular diversity, in addition to OR genes, we revealed nearly 200 differentially expressed genes among the sequenced OSNs in this study. Among the 36 sequenced OSNs, eight cells (22.2%) showed multiple OR gene expression and the presences of additional ORs were not restricted to the neighbor loci that shared the transcriptional effect of the primary OR expression, suggesting that the “One neuron-one receptor rule” might not be strictly true at the transcription level. All of the inferable ORs, including additional co-expressed ORs, were shown to be monoallelic. Our sequencing of 21 Mecp2308 mutant OSNs, of which 62% expressed more than one OR genes, and the expression levels of the additional ORs were significantly higher than those in the wild-type, suggested that MeCP2 plays a role in the regulation of singular OR gene expression. Dual label in situ hybridization along with the sequence data revealed that dorsal and ventral ORs were co-expressed in the same Mecp2 mutant OSN, further implying that MeCP2 might be involved in regulation of OR territories in the OE. Our results suggested a new role of MeCP2 in OR gene choice and ratified that this multiple-OR expression caused by Mecp2 mutation did not accompany delayed OSN development that has been observed in the previous studies on the Mecp2 mutants.

Unraveling the function of bacterial-type phosphoenolpyruvate carboxylase in vascular plants

Two distinct phosphoenolpyruvate carboxylase (PEPC) isozymes occur in vascular plants and green algae: plant-type PEPC (PTPC) and bacterial-type PEPC (BTPC). PTPC polypeptides typically form a tightly regulated cytosolic Class-1 PEPC homotetramer. BTPCs, however, appear to be less widely expressed and to exist only as catalytic and regulatory subunits that physically interact with co-expressed PTPC subunits to form hetero-octameric Class-2 PEPC complexes that are highly desensitized to Class-1 PEPC allosteric effectors. Yeast two-hybrid studies indicated that castor plant BTPC (RcPPC4) interacts with all three Arabidopsis thaliana PTPC isozymes, and that it forms stronger interactions with AtPPC2 and AtPPC3, suggesting that specific PTPCs are preferred for Class-2 PEPC formation. In contrast, Arabidopsis BTPC (AtPPC4) appeared to interact very weakly with AtPPC2 and AtPPC3, suggesting that BTPCs from different species may have different physical properties, hypothesized to be due to sequence dissimilarities within their ~10 kDa intrinsically disordered region. Recent RNA-seq and microarray data were analyzed to obtain a better understanding of BTPC expression patterns in different tissues of various monocot and dicot species. High levels of BTPC transcripts, polypeptides and Class-2 PEPC complexes were originally discovered in developing castor seeds, but the analysis revealed a broad range of diverse tissues where abundant BTPC transcripts are also expressed, such as the developing fruits of cucumber, grape, and tomato. Marked BTPC expression correlated well with the presence of ~116 kDa immunoreactive BTPC polypeptides, as well as Class-2 PEPC complexes in the immature fruit of cucumbers and tomatoes. It is therefore hypothesized that in vascular plants BTPC and thus Class-2 PEPC complexes maintain anaplerotic PEP flux in tissues with elevated malate levels that would potently inhibit ‘housekeeping’ Class-1 PEPCs. Elevated levels of malate can be used by biosynthetically active sink tissues such as immature tomatoes and cucumbers for rapid cell expansion, drought or salt stressed roots for osmoregulation, and developing seeds and pollen as a precursor for storage lipid and protein biosynthesis.

REST is a hypoxia-responsive transcriptional repressor

Cellular exposure to hypoxia results in altered gene expression in a range of physiologic and pathophysiologic states. Discrete cohorts of genes can be either up- or down-regulated in response to hypoxia. While the Hypoxia-Inducible Factor (HIF) is the primary driver of hypoxia-induced adaptive gene expression, less is known about the signalling mechanisms regulating hypoxia-dependent gene repression. Using RNA-seq, we demonstrate that equivalent numbers of genes are induced and repressed in human embryonic kidney (HEK293) cells. We demonstrate that nuclear localization of the Repressor Element 1-Silencing Transcription factor (REST) is induced in hypoxia and that REST is responsible for regulating approximately 20% of the hypoxia-repressed genes. Using chromatin immunoprecipitation assays we demonstrate that REST-dependent gene repression is at least in part mediated by direct binding to the promoters of target genes. Based on these data, we propose that REST is a key mediator of gene repression in hypoxia.

Linfócitos T : uma análise do perfil diferencial da expressão gênica na imunorregulação

Tese (doutorado)—Universidade de Brasília, Faculdade de Medicina, Pós-Graduação em Patologia Molecular, 2016.

Pseudomonas Aeruginosa AmpR Transcriptional Regulatory Network

In Enterobacteriaceae, the transcriptional regulator AmpR, a member of the LysR family, regulates the expression of a chromosomal β-lactamase AmpC. The regulatory repertoire of AmpR is broader in Pseudomonas aeruginosa, an opportunistic pathogen responsible for numerous acute and chronic infections including cystic fibrosis. Previous studies showed that in addition to regulating ampC, P. aeruginosa AmpR regulates the sigma factor AlgT/U and production of some quorum sensing (QS)-regulated virulence factors. In order to better understand the ampR regulon, the transcriptional profiles generated using DNA microarrays and RNA-Seq of the prototypic P. aeruginosa PAO1 strain with its isogenic ampR deletion mutant, PAO∆ampR were analyzed. Transcriptome analysis demonstrates that the AmpR regulon is much more extensive than previously thought influencing the differential expression of over 500 genes. In addition to regulating resistance to β-lactam antibiotics via AmpC, AmpR also regulates non-β-lactam antibiotic resistance by modulating the MexEF-OprN efflux pump. Virulence mechanisms including biofilm formation, QS-regulated acute virulence, and diverse physiological processes such as oxidative stress response, heat-shock response and iron uptake are AmpR-regulated. Real-time PCR and phenotypic assays confirmed the transcriptome data. Further, Caenorhabditis elegans model demonstrates that a functional AmpR is required for full pathogenicity of P. aeruginosa. AmpR, a member of the core genome, also regulates genes in the regions of genome plasticity that are acquired by horizontal gene transfer. The extensive AmpR regulon included other transcriptional regulators and sigma factors, accounting for the extensive AmpR regulon. Gene expression studies demonstrate AmpR-dependent expression of the QS master regulator LasR that controls expression of many virulence factors. Using a chromosomally tagged AmpR, ChIP-Seq studies show direct AmpR binding to the lasR promoter. The data demonstrates that AmpR functions as a global regulator in P. aeruginosa and is a positive regulator of acute virulence while negatively regulating chronic infection phenotypes. In summary, my dissertation sheds light on the complex regulatory circuit in P. aeruginosa to provide a better understanding of the bacterial response to antibiotics and how the organism coordinately regulates a myriad of virulence factors.

Understanding ten-eleven translocation-2 in hematological and nervous systems

I proposed the study of two distinct aspects of Ten-Eleven Translocation 2 (TET2) protein for understanding specific functions in different body systems. ^ In Part I, I characterized the molecular mechanisms of Tet2 in the hematological system. As the second member of Ten-Eleven Translocation protein family, TET2 is frequently mutated in leukemic patients. Previous studies have shown that the TET2 mutations frequently occur in 20% myelodysplastic syndrome/myeloproliferative neoplasm (MDS/MPN), 10% T-cell lymphoma leukemia and 2% B-cell lymphoma leukemia. Genetic mouse models also display distinct phenotypes of various types of hematological malignancies. I performed 5-hydroxymethylcytosine (5hmC) chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq) of hematopoietic stem/progenitor cells to determine whether the deletion of Tet2 can affect the abundance of 5hmC at myeloid, T-cell and B-cell specific gene transcription start sites, which ultimately result in various hematological malignancies. Subsequent Exome sequencing (Exome-Seq) showed that disease-specific genes are mutated in different types of tumors, which suggests that TET2 may protect the genome from being mutated. The direct interaction between TET2 and Mutator S Homolog 6 (MSH6) protein suggests TET2 is involved in DNA mismatch repair. Finally, in vivo mismatch repair studies show that the loss of Tet2 causes a mutator phenotype. Taken together, my data indicate that TET2 binds to MSH6 to protect genome integrity. ^ In Part II, I intended to better understand the role of Tet2 in the nervous system. 5-hydroxymethylcytosine regulates epigenetic modification during neurodevelopment and aging. Thus, Tet2 may play a critical role in regulating adult neurogenesis. To examine the physiological significance of Tet2 in the nervous system, I first showed that the deletion of Tet2 reduces the 5hmC levels in neural stem cells. Mice lacking Tet2 show abnormal hippocampal neurogenesis along with 5hmC alternations at different gene promoters and corresponding gene expression downregulation. Through the luciferase reporter assay, two neural factors Neurogenic differentiation 1 (NeuroD1) and Glial fibrillary acidic protein (Gfap) were down-regulated in Tet2 knockout cells. My results suggest that Tet2 regulates neural stem/progenitor cell proliferation and differentiation in adult brain.^

Unraveling the function of bacterial-type phosphoenolpyruvate carboxylase in vascular plants

Thesis (Master, Biology) -- Queen's University, 2016-09-29 20:09:46.997

Mapping new non-genetic dependencies in malignant pleural mesothelioma

Malignant Pleural Mesothelioma (MPM) is a very aggressive cancer whose incidence is growing worldwide. MPM escapes the classical models of carcinogenesis and lacks a distinctive genetic fingerprint, keeping obscure the molecular events that lead to tumorigenesis. This severely impacts on the limited therapeutic options and on the lack of specific biomarkers, concurring to make MPM one of the deadliest cancers. Here we combined a functional genome-wide loss of function CRISPR/Cas9 screening with patients’ transcriptomic and clinical data, to identify genes essential for MPM progression. Besides, we explored the role of non-coding RNAs to MPM progression by analysing gene expression profiles and clinical data from the MESO-TCGA dataset. We identified TRIM28 and the lncRNA LINC00941 as new vulnerabilities of MPM, associated with disease aggressiveness and bad outcome of patients. TRIM28 is a multi-domain protein involved in many processes, including transcription regulation. We showed that TRIM28 silencing impairs MPM cells’ growth and clonogenicity by blocking cells in mitosis. RNA-seq profiling showed that TRIM28 loss abolished the expression of major mitotic players. Our data suggest that TRIM28 is part of the B-MYB/FOXM1-MuvB complex that specifically drives the activation of mitotic genes, keeping the time of mitosis. In parallel, we found LINC00941 as strongly associated with reduced survival probability in MPM patients. LINC00941 KD profoundly reduced MPM cells’ growth, migration and invasion. This is accompanied by changes in morphology, cytoskeleton organization and cell-cell adhesion properties. RNA-seq profiling showed that LINC00941 KD impacts crucial functions of MPM, including HIF1α signalling. Collectively these data provided new insights into MPM biology and demonstrated that the integration of functional screening with patients’ clinical data is a powerful tool to highlight new non-genetic cancer dependencies that associate to a bad outcome in vivo, paving the way to new MPM-oriented targeted strategies and prognostic tools to improve patients risk-based stratification.

Analysis of extracellular vesicles from patients with advanced pancreatic cancer identifies miRNAs with predictive value for treatment with gemcitabine + nab-paclitaxel

Pancreatic cancer (PC) is the seventh leading cause of cancer death. Despite recent therapy advancements, 5-year survival is 11%. Resistance to therapy is common, and no predictive factors, except for BRCA1/2 and PALB2 mutations, can drive treatment selection. Based on the easy isolation of extracellular vesicles (EVs) from blood and the role of EV-borne miRNAs in chemoresistance, we analyzed EVs and their miRNA content in order to identify predictive factors. First, we analyzed samples from 28 PC patients and 7 healthy subjects, in order to establish methods for isolation and analysis of EVs and their miRNA content. We observed a significantly different expression of 28 miRNAs, including oncogenic or tumor suppressor miRNAs, showing the ability of our approach to detect candidate biomarkers. Then, we analyzed samples of 21 advanced PC patients, collected before first-line treatment with gemcitabine + nab-paclitaxel, and compared findings in responders and non-responders. EVs have been analyzed with Nanoparticle tracking analysis, flow cytometry and RNA-Seq; then, laboratory results have been matched with clinical data. Nanoparticle tracking analysis did not show any significant difference. Flow cytometry showed a lower expression of SSE4 and CD81 in responders. Finally, miRNA analysis showed 25 upregulated and 19 downregulated miRNAs in responders. In particular, in responders we observed upregulation of miR-141-3p, miR-141-5p, miR-200a-3p, miR-200b-3p, miR-200c-3p, miR-375-3p, miR-429, miR-545-5p. These miRNAs have targets with a previously reported role in PC. In conclusion, we show the feasibility of the proposed approach to identify EV-derived biomarkers with predictive value for therapy with gemcitabine + nab-paclitaxel in PC. Our findings highlight the possibility to exploit liquid biopsy for personalized treatment in PC, in order to maximize chances of response and patients’ outcome. These findings are worthy of further investigation: in the same setting, with different chemotherapy schedules, and in different disease settings such as preoperative therapy.

Characterisation of RNase Y in H. pylori: a non-essential enzyme involved in the processing of cag-PAI non coding RNA 1 (CncR1) sRNA

The importance of Helicobacter pylori as a human pathogen is underlined by the plethora of diseases it is responsible for. The capacity of H. pylori to adapt to the restricted host-associated environment andto evade the host immune response largely depends on a streamlined signalling network. The peculiar H. pylori small genome size combined with its paucity of transcriptional regulators highlights the relevance of post-transcriptional regulatory mechanisms as small non-coding RNAs (sRNAs). However, among the 8 RNases represented in H. pylori genome, a regulator guiding sRNAs metabolism is still not well studied. We investigated for the first time the physiological role in H. pylori G27 strain of the RNase Y enzyme. In the first line of research we provide a comprehensive characterization of the RNase Y activity by analysing its genomic organization and the factors that orchestrate its expression. Then, based on bioinformatic prediction models, we depict the most relevant determinants of RNase Y function, demonstrating a correlation of both structure and domain organization with orthologues represented in Gram-positive bacteria. To unveil the post-transcriptional regulatory effect exerted by the RNase Y, we compared the transcriptome of an RNase Y knock-out mutant to the parental wild type strain by RNA-seq approach. In the second line of research we characterized the activity of this single strand specific endoribonuclease on cag-PAI non coding RNA 1 (CncR1) sRNA. We found that deletion or inactivation of RNase Y led to the accumulation of a 3’-extended CncR1 (CncR1-L) transcript over time. Moreover, beneath its increased half-life, CncR1-L resembled a CncR1 inactive phenotype. Finally, we focused on the characterization of the in vivo interactome of CncR1. We set up a preliminary MS2-affinity purification coupled with RNA-sequencing (MAPS) approach and we evaluated the enrichment of specific targets, demonstrating the suitability of the technique in the H. pylori G27 strain.

De Novo Assembly And Transcriptome Analysis Of The Rubber Tree (hevea Brasiliensis) And Snp Markers Development For Rubber Biosynthesis Pathways.

Hevea brasiliensis (Willd. Ex Adr. Juss.) Muell.-Arg. is the primary source of natural rubber that is native to the Amazon rainforest. The singular properties of natural rubber make it superior to and competitive with synthetic rubber for use in several applications. Here, we performed RNA sequencing (RNA-seq) of H. brasiliensis bark on the Illumina GAIIx platform, which generated 179,326,804 raw reads on the Illumina GAIIx platform. A total of 50,384 contigs that were over 400 bp in size were obtained and subjected to further analyses. A similarity search against the non-redundant (nr) protein database returned 32,018 (63%) positive BLASTx hits. The transcriptome analysis was annotated using the clusters of orthologous groups (COG), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Pfam databases. A search for putative molecular marker was performed to identify simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs). In total, 17,927 SSRs and 404,114 SNPs were detected. Finally, we selected sequences that were identified as belonging to the mevalonate (MVA) and 2-C-methyl-D-erythritol 4-phosphate (MEP) pathways, which are involved in rubber biosynthesis, to validate the SNP markers. A total of 78 SNPs were validated in 36 genotypes of H. brasiliensis. This new dataset represents a powerful information source for rubber tree bark genes and will be an important tool for the development of microsatellites and SNP markers for use in future genetic analyses such as genetic linkage mapping, quantitative trait loci identification, investigations of linkage disequilibrium and marker-assisted selection.

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.

Análise transcriptômica da interação mamoeiro-Papaya Meleira Virus

O mamoeiro (Carica papaya L.) é uma das fruteiras mais cultivadas nas regiões tropicais e subtropicais do mundo. O Brasil faz parte do grupo dos países que mais produzem e exportam mamão no mundo. O Espírito Santo e a Bahia são responsáveis por mais de 70% da área brasileira produtora deste fruto. Porém, doenças causadas por microrganismos infecciosos afetam de modo considerável sua produção. Entre as principais doenças, destaca-se a meleira do mamoeiro, causada pelo Papaya meleira virus (PMeV), que ainda não possui uma cultivar resistente. Interessantemente os sintomas somente são desencadeados após a frutificação. Os mecanismos moleculares envolvidos no desenvolvimento dos sintomas e na resposta de defesa da planta ao PMeV ainda não foram esclarecidos. Para entender os pontos chaves desta interação, que permitam o desenvolvimento de metodologias de melhoramento genético, um estudo transcriptômico foi abordado. A tecnologia RNA-seq foi usada para o sequenciamento do transcriptoma a partir de plantas com 3, 6 e 8 meses de idade após plantio, inoculadas e não inoculadas com o PMeV. Os genes diferencialmente expressos nos 3 tempos e nas duas condições foram preditos e analisados. Estas análises revelaram um padrão de expressão geral dos genes envolvidos nesta interação. Foram encontrados 21 genes com o perfil de expressão alterado nas plantas inoculadas exclusivamente nos seis meses de idade. Destes, 8 genes envolvidos em processos de respostas de defesa e morte celular, resposta ao estresse e resposta ao estímulo biótico e abiótico foram reprimidos; enquanto os demais (13 genes), envolvidos principalmente em processos metabólicos primários, biogêneses, diferenciação e ciclo celular, comunicação e crescimento celular, bem como processos envolvidos em reprodução, e desenvolvimento da floração, foram superexpressos. Estes resultados sugerem que, aos seis meses de idade, a planta é obrigada a alterar seu programa de expressão gênica, direcionando a resposta para os processos próprios do desenvolvimento, requeridos nesse estádio fisiológico, que primam sob a resposta ao estresse, fato que finalmente leva ao desenvolvimento dos sintomas.

Host-symbiont interactions in the deep-sea vent mussel Bathymodiolus azoricus : a molecular approach

Tese de Doutoramento, Ciências do Mar, especialidade de Biologia Marinha, 19 de Dezembro de 2015, Universidade dos Açores.