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

Genome-wide meta-analysis identifies novel multiple sclerosis susceptibility loci

Objective: To perform a 1-stage meta-analysis of genome-wide association studies (GWAS) of multiple sclerosis (MS) susceptibility and to explore functional consequences of new susceptibility loci. Methods: We synthesized 7 MS GWAS. Each data set was imputed using HapMap phase II, and a per single nucleotide polymorphism (SNP) meta-analysis was performed across the 7 data sets. We explored RNA expression data using a quantitative trait analysis in peripheral blood mononuclear cells (PBMCs) of 228 subjects with demyelinating disease. Results: We meta-analyzed 2,529,394 unique SNPs in 5,545 cases and 12,153 controls. We identified 3 novel susceptibility alleles: rs170934T at 3p24.1 (odds ratio [OR], 1.17; p ¼ 1.6 � 10�8) near EOMES, rs2150702G in the second intron of MLANA on chromosome 9p24.1 (OR, 1.16; p ¼ 3.3 � 10�8), and rs6718520A in an intergenic region on chromosome 2p21, with THADA as the nearest flanking gene (OR, 1.17; p ¼ 3.4 � 10�8). The 3 new loci do not have a strong cis effect on RNA expression in PBMCs. Ten other susceptibility loci had a suggestive p < 1 � 10�6, some of these loci have evidence of association in other inflammatory diseases (ie, IL12B, TAGAP, PLEK, and ZMIZ1). Interpretation: We have performed a meta-analysis of GWAS in MS that more than doubles the size of previous gene discovery efforts and highlights 3 novel MS susceptibility loci. These and additional loci with suggestive evidence of association are excellent candidates for further investigations to refine and validate their role in the genetic architecture of MS.

Sexual selection in true fruit flies (Diptera: Tephritidae) : transcriptome and experimental evidences for phytochemicals increasing male competitive ability

In male tephritid fruit flies of the genus Bactrocera, feeding on secondary plant compounds (sensu lato male lures = methyl eugenol, raspberry ketone and zingerone) increases male mating success. Ingested male lures alter the male pheromonal blend, normally making it more attractive to females and this is considered the primary mechanism for the enhanced mating success. However, the male lures raspberry ketone and zingerone are known, across a diverse range of other organisms, to be involved in increasing energy metabolism. If this also occurs in Bactrocera, then this may represent an additional benefit to males as courtship is metabolically expensive and lure feeding may increase a fly's short-term energy. We tested this hypothesis by performing comparative RNA-seq analysis between zingerone-fed and unfed males of Bactrocera tryoni. We also carried out behavioural assays with zingerone- and cuelure-fed males to test whether they became more active. RNA-seq analysis revealed, in zingerone-fed flies, up-regulation of 3183 genes with homologues transcripts to those known to regulate intermale aggression, pheromone synthesis, mating and accessory gland proteins, along with significant enrichment of several energy metabolic pathways and gene ontology terms. Behavioural assays show significant increases in locomotor activity, weight reduction and successful mating after mounting; all direct/indirect measures of increased activity. These results suggest that feeding on lures leads to complex physiological changes, which result in more competitive males. These results do not negate the pheromone effect, but do strongly suggest that the phytochemical-induced sexual selection is governed by both female preference and male competitive mechanisms.

Combining transcriptome assemblies from multiple de novo assemblers in the allo-tetraploid plant Nicotiana benthamiana

Background Nicotiana benthamiana is an allo-tetraploid plant, which can be challenging for de novo transcriptome assemblies due to homeologous and duplicated gene copies. Transcripts generated from such genes can be distinct yet highly similar in sequence, with markedly differing expression levels. This can lead to unassembled, partially assembled or mis-assembled contigs. Due to the different properties of de novo assemblers, no one assembler with any one given parameter space can re-assemble all possible transcripts from a transcriptome. Results In an effort to maximise the diversity and completeness of de novo assembled transcripts, we utilised four de novo transcriptome assemblers, TransAbyss, Trinity, SOAPdenovo-Trans, and Oases, using a range of k-mer sizes and different input RNA-seq read counts. We complemented the parameter space biologically by using RNA from 10 plant tissues. We then combined the output of all assemblies into a large super-set of sequences. Using a method from the EvidentialGene pipeline, the combined assembly was reduced from 9.9 million de novo assembled transcripts to about 235,000 of which about 50,000 were classified as primary. Metrics such as average bit-scores, feature response curves and the ability to distinguish paralogous or homeologous transcripts, indicated that the EvidentialGene processed assembly was of high quality. Of 35 RNA silencing gene transcripts, 34 were identified as assembled to full length, whereas in a previous assembly using only one assembler, 9 of these were partially assembled. Conclusions To achieve a high quality transcriptome, it is advantageous to implement and combine the output from as many different de novo assemblers as possible. We have in essence taking the ‘best’ output from each assembler while minimising sequence redundancy. We have also shown that simultaneous assessment of a variety of metrics, not just focused on contig length, is necessary to gauge the quality of assemblies.

Identification of a novel fusion transcript between human relaxin-1 (RLN1) and human relaxin-2 (RLN2) in prostate cancer

Simultaneous expression of highly homologous RLN1 and RLN2 genes in prostate impairs their accurate delineation. We used PacBio SMRT sequencing and RNA-Seq in LNCaP cells in order to dissect the expression of RLN1 and RLN2 variants. We identified a novel fusion transcript comprising the RLN1 and RLN2 genes and found evidence of its expression in the normal and prostate cancer tissues. The RLN1-RLN2 fusion putatively encodes RLN2 isoform with the deleted secretory signal peptide. The identification of the fusion transcript provided information to determine unique RLN1-RLN2 fusion and RLN1 regions. The RLN1-RLN2 fusion was co-expressed with RLN1 in LNCaP cells, but the two gene products were inversely regulated by androgens. We showed that RLN1 is underrepresented in common PCa cell lines in comparison to normal and PCa tissue. The current study brings a highly relevant update to the relaxin field, and will encourage further studies of RLN1 and RLN2 in PCa and broader.

Sexual selection in true fruit flies (Diptera: Tephritidae): transcriptome and experimental evidences for phytochemicals increasing male competitive ability

In male tephritid fruit flies of the genus Bactrocera, feeding on secondary plant compounds (sensu lato male lures = methyl eugenol, raspberry ketone and zingerone) increases male mating success. Ingested male lures alter the male pheromonal blend, normally making it more attractive to females and this is considered the primary mechanism for the enhanced mating success. However, the male lures raspberry ketone and zingerone are known, across a diverse range of other organisms, to be involved in increasing energy metabolism. If this also occurs in Bactrocera, then this may represent an additional benefit to males as courtship is metabolically expensive and lure feeding may increase a fly's short-term energy. We tested this hypothesis by performing comparative RNA-seq analysis between zingerone-fed and unfed males of Bactrocera tryoni. We also carried out behavioural assays with zingerone- and cuelure-fed males to test whether they became more active. RNA-seq analysis revealed, in zingerone-fed flies, up-regulation of 3183 genes with homologues transcripts to those known to regulate intermale aggression, pheromone synthesis, mating and accessory gland proteins, along with significant enrichment of several energy metabolic pathways and gene ontology terms. Behavioural assays show significant increases in locomotor activity, weight reduction and successful mating after mounting; all direct/indirect measures of increased activity. These results suggest that feeding on lures leads to complex physiological changes, which result in more competitive males. These results do not negate the pheromone effect, but do strongly suggest that the phytochemical-induced sexual selection is governed by both female preference and male competitive mechanisms.

Identification of I-7 expands the repertoire of genes for resistance to Fusarium wilt in tomato to three resistance gene classes

The tomato I-3 and I-7 genes confer resistance to Fusarium oxysporum f. sp. lycopersici (Fol) race 3 and were introgressed into the cultivated tomato, Solanum lycopersicum, from the wild relative Solanum pennellii. I-3 has been identified previously on chromosome 7 and encodes an S-receptor-like kinase, but little is known about I-7. Molecular markers have been developed for the marker-assisted breeding of I-3, but none are available for I-7. We used an RNA-seq and single nucleotide polymorphism (SNP) analysis approach to map I-7 to a small introgression of S. pennellii DNA (c. 210 kb) on chromosome 8, and identified I-7 as a gene encoding a leucine-rich repeat receptor-like protein (LRR-RLP), thereby expanding the repertoire of resistance protein classes conferring resistance to Fol. Using an eds1 mutant of tomato, we showed that I-7, like many other LRR-RLPs conferring pathogen resistance in tomato, is EDS1 (Enhanced Disease Susceptibility 1) dependent. Using transgenic tomato plants carrying only the I-7 gene for Fol resistance, we found that I-7 also confers resistance to Fol races 1 and 2. Given that Fol race 1 carries Avr1, resistance to Fol race 1 indicates that I-7-mediated resistance, unlike I-2- or I-3-mediated resistance, is not suppressed by Avr1. This suggests that Avr1 is not a general suppressor of Fol resistance in tomato, leading us to hypothesize that Avr1 may be acting against an EDS1-independent pathway for resistance activation. The identification of I-7 has allowed us to develop molecular markers for marker-assisted breeding of both genes currently known to confer Fol race 3 resistance (I-3 and I-7). Given that I-7-mediated resistance is not suppressed by Avr1, I-7 may be a useful addition to I-3 in the tomato breeder's toolbox.

The C. elegans ALA neuron : its transcriptome and role in inducing sleep

A long-standing yet to be accomplished task in understanding behavior is to dissect the function of each gene involved in the development and function of a neuron. The C. elegans ALA neuron was chosen in this study for its known function in sleep, an ancient but less understood animal behavior. Single-cell transcriptome profiling identified 8,133 protein-coding genes in the ALA neuron, of which 57 are neuropeptide-coding genes. The most enriched genes are also neuropeptides. In combination with gain-of-function and loss-of-function assays, here I showed that the ALA-enriched FMRFamide neuropeptides, FLP-7, FLP-13, and FLP-24, are sufficient and necessary for inducing C. elegans sleep. These neuropeptides act as neuromodulators through GPCRs, NPR-7, and NPR-22. Further investigation in zebrafish indicates that FMRFamide neuropeptides are sleep-promoting molecules in animals. To correlate the behavioral outputs with genomic context, I constructed a gene regulatory network of the relevant genes controlling C. elegans sleep behavior through EGFR signaling in the ALA neuron. First, I identified an ALA cell-specific motif to conduct a genome-wide search for possible ALA-expressed genes. I then filtered out non ALA-expressed genes by comparing the motif-search genes with ALA transcriptomes from single-cell profiling. In corroborating with ChIP-seq data from modENCODE, I sorted out direct interaction of ALA-expressed transcription factors and differentiation genes in the EGFR sleep regulation pathway. This approach provides a network reference for the molecular regulation of C. elegans sleep behavior, and serves as an entry point for the understanding of functional genomics in animal behaviors.

Global transcriptional responses of the toxic cyanobacterium, Microcystis aeruginosa, to nitrogen stress, phosphorus stress, and growth on organic matter

Whole transcriptome shotgun sequencing (RNA-seq) was used to assess the transcriptomic response of the toxic cyanobacterium Microcystis aeruginosa during growth with low levels of dissolved inorganic nitrogen (low N), low levels of dissolved inorganic phosphorus (low P), and in the presence of high levels of high molecular weight dissolved organic matter (HMWDOM). Under low N, one third of the genome was differentially expressed, with significant increases in transcripts observed among genes within the nir operon, urea transport genes (urtBCDE), and amino acid transporters while significant decreases in transcripts were observed in genes related to photosynthesis. There was also a significant decrease in the transcription of the microcystin synthetase gene set under low N and a significant decrease in microcystin content per Microcystis cell demonstrating that N supply influences cellular toxicity. Under low P, 27% of the genome was differentially expressed. The Pho regulon was induced leading to large increases in transcript levels of the alkaline phosphatase phoX, the Pst transport system (pstABC), and the sphX gene, and transcripts of multiple sulfate transporter were also significantly more abundant. While the transcriptional response to growth on HMWDOM was smaller (5–22% of genes differentially expressed), transcripts of multiple genes specifically associated with the transport and degradation of organic compounds were significantly more abundant within HMWDOM treatments and thus may be recruited by Microcystis to utilize these substrates. Collectively, these findings provide a comprehensive understanding of the nutritional physiology of this toxic, bloom-forming cyanobacterium and the role of N in controlling microcystin synthesis.

Visualising ribosome profiling and using it for reading frame detection and exploration of eukaryotic translation initiation

Ribosome profiling (ribo-seq) is a recently developed technique that provides genomewide information on protein synthesis (GWIPS) in vivo. The high resolution of ribo-seq is one of the exciting properties of this technique. In Chapter 2, I present a computational method that utilises the sub-codon precision and triplet periodicity of ribosome profiling data to detect transitions in the translated reading frame. Application of this method to ribosome profiling data generated for human HeLa cells allowed us to detect several human genes where the same genomic segment is translated in more than one reading frame. Since the initial publication of the ribosome profiling technique in 2009, there has been a proliferation of studies that have used the technique to explore various questions with respect to translation. A review of the many uses and adaptations of the technique is provided in Chapter 1. Indeed, owing to the increasing popularity of the technique and the growing number of published ribosome profiling datasets, we have developed GWIPS-viz (http://gwips.ucc.ie), a ribo-seq dedicated genome browser. Details on the development of the browser and its usage are provided in Chapter 3. One of the surprising findings of ribosome profiling of initiating ribosomes carried out in 3 independent studies, was the widespread use of non-AUG codons as translation initiation start sites in mammals. Although initiation at non-AUG codons in mammals has been documented for some time, the extent of non-AUG initiation reported by these ribo-seq studies was unexpected. In Chapter 4, I present an approach for estimating the strength of initiating codons based on the leaky scanning model of translation initiation. Application of this approach to ribo-seq data illustrates that initiation at non-AUG codons is inefficient compared to initiation at AUG codons. In addition, our approach provides a probability of initiation score for each start site that allows its strength of initiation to be evaluated.

The relevance of bacterial isoprenoid biosynthetic pathways for host-microbe interactions

This thesis was undertaken to investigate the relevance of two bacterial isoprenoid biosynthetic pathways (Mevalonate (MVAL) and 2-C-methyl-D-erythritol 4-phosphate (MEP)) for host-microbe interactions. We determined a significant reduction in microbial diversity in the murine gut microbiota (by next generation sequencing) following oral administration of a common anti-cholesterol drug Rosuvastatin (RSV) that targets mammalian and bacterial HMG-CoA reductase (HMG-R) for inhibition of MVAL formation. In tandem we identified significant hepatic and intestinal off-target alterations to the murine metabolome indicating alterations in inflammation, bile acid profiles and antimicrobial peptide synthesis with implications on community structure of the gastrointestinal microbiota in statin-treated animals. However we found no effect on local Short Chain Fatty Acid biosynthesis (metabolic health marker in our model). We demonstrated direct inhibition of bacterial growth in-vitro by RSV which correlated with reductions in bacterial MVAL formation. However this was only at high doses of RSV. Our observations demonstrate a significant RSV-associated impact on the gut microbiota prompting similar human analysis. Successful deletion of another MVAL pathway enzyme (HMG-CoA synthase (mvaS)) involved in Listeria monocytogenes EGDe isoprenoid biosynthesis determined that the enzyme is non-essential for normal growth and in-vivo pathogenesis of this pathogen. We highlight potential evidence for alternative means of synthesis of the HMG-CoA substrate that could render mvaS activity redundant under our test conditions. Finally, we showed by global gene expression analysis (Massive Analysis of cDNA Ends (MACE RNA-seq) a significant role for the penultimate MEP pathway metabolite (E)-4-hydroxy-3-methyl-2-but-2-enyl pyrophosphate (HMBPP) in significant up regulation of genes of immunity and antigen presentation in THP-1 cells at nanomolar levels. We infected THP-1 cells with wild type or HMBPP under/over-producing L. monoctyogenes EGDe mutants and determined subtle effects of HMBPP upon overall host responses to Listeria infection. Overall our findings provide greater insights regarding bacterial isoprenoid biosynthetic pathways for host-microbe/microbe-host dialogue.

Hybrid error correction and de novo assembly of single-molecule sequencing reads.

Single-molecule sequencing instruments can generate multikilobase sequences with the potential to greatly improve genome and transcriptome assembly. However, the error rates of single-molecule reads are high, which has limited their use thus far to resequencing bacteria. To address this limitation, we introduce a correction algorithm and assembly strategy that uses short, high-fidelity sequences to correct the error in single-molecule sequences. We demonstrate the utility of this approach on reads generated by a PacBio RS instrument from phage, prokaryotic and eukaryotic whole genomes, including the previously unsequenced genome of the parrot Melopsittacus undulatus, as well as for RNA-Seq reads of the corn (Zea mays) transcriptome. Our long-read correction achieves >99.9% base-call accuracy, leading to substantially better assemblies than current sequencing strategies: in the best example, the median contig size was quintupled relative to high-coverage, second-generation assemblies. Greater gains are predicted if read lengths continue to increase, including the prospect of single-contig bacterial chromosome assembly.