A teca é uma árvore de grande importância comercial pelas características de cor e durabilidade de sua madeira. Devido a sua rusticidade e fácil adaptação ao clima, plantios de teca tornam-se cada vez mais atrativos ao redor do mundo. Contudo, esta espécie apresenta escassez de estudos genéticos moleculares a respeito tanto de sua madeira, quanto de sua tolerância às variações ambientais. Uma vez que o transcriptoma pode apresentar grande quantidade de informação a respeito dos genes expressos por um conjunto celular, neste trabalho foi realizado o primeiro transcriptoma de teca, onde foram sequenciadas flores, folhas, raízes e seedlings pela tecnologia Illumina. A montagem do transcriptoma foi realizada com o programa Trinity acima de 100 milhões de reads e gerou mais de 400 mil contigs, os quais tiveram as anotações funcionais adquiridas com o programa Blas2GO. 51% dos contigs foram anotados, mostrando alta similaridade com as espécies Vitis vinifera e Solanum licopersicum; destes, 78% obtiveram anotações funcionais com o Gene Ontology, totalizando 5.165 termos para Processo Biológico, 2.846 termos para Função Molecular e 742 para Componente Celular. A expressão diferencial foi obtida com o programa edgeR a 5% de probabilidade de erro e mostrou que, para 187.315 contigs montados através da fusão de todas as bibliotecas sequenciadas, 18 mostraram expressão diferencial para flor, 14 para folha, 13 para raiz e 29 para seedling. Após a etapa de caracterização do transcriptoma, foi realizado um experimento de estresse por déficit hídrico em casa-de-vegetação, onde plantas de teca foram submetidas a estresse Moderado (40% de água no substrato por 20 dias), estresse Severo (20 a 40% de água por 30 dias) e tratamento controle (substrato saturado). As medições através de analisador de gases por infravermelho (IRGA) mostraram queda na fotossíntese (até 70% a menos do que o controle), na transpiração (até 77%) e na condutância estomática (até 85%) entre os tratamentos; além disto, o conteúdo relativo de água foliar caiu 13% entre o tratamento severo e o controle, e níveis de prolina livre foram até 3,5 vezes mais altos nos tratamentos de estresse. A temperatura foliar aumentou significativamente com o aumento da irradiância de fótons aplicada. A busca por genes relacionados ao estresse por déficit hídrico na biblioteca de transcritos de Raiz retornou 1.145 sequências, e destas, 4 foram caracterizadas: TgTPS (trealose 6-fosfato sintase), TgPIP (aquaporina, proteína intrínseca de membrana plasmática), TgDREB2 (proteína de ligação a elemento responsivo a desidratação) e TgAREB (proteína de ligação a elemento responsivo a ácido abscísico). Apenas TgTPS, TgPIP e TgDREB2 mostraram alto grau de conservação entre as espécies, podendo ser corretamente amplificadas via PCR e validadas por sequenciamento. Assim, com o banco de dados de transcritos obtido pelo RNA-seq, foi possível identificar genes candidatos ao estudo de características vegetativas e reprodutivas de teca, contribuindo para entender os mecanismos moleculares desta espécie florestal. / Teak is a tree of great commercial importance by the characteristics of color and durability of its wood. Due to its hardiness and easy adaptation to climate, teak plantations become increasingly attractive around the world. However, this species has a lack of molecular genetic studies on both of its wood, as their tolerance to environmental variations. Once the transcriptome can provide lots of information about the genes expressed by a cell group, this work represents the first transcriptome teak, which were sequenced flowers, leaves, roots and seedlings by Illumina technology. The transcriptome assembly was performed with Trinity program above 100 million reads and generated more than 400,000 contigs, which have acquired the functional annotations with Blas2GO program. 51% of the contigs were annotaded, showing high similarity to Vitis vinifera and Solanum licopersicum; of these, 78% had functional annotations with the Gene Ontology, totaling 5,165 terms for Biological Process, 2846 terms for Molecular Function and 742 for Cell Component. The differential expression was obtained with the edgeR program at 5% probability of error and showed that for 187,315 contigs assembled by merging all sequenced libraries, 18 showed differential expression to flower, 14 to leaf, 13 to root and 29 for seedling. After this step of characterization of the transcriptome, we performed a stress experiment by water deficit at greenhouse, where teak plants were subjected to Moderate stress (40% of water in the substrate for 20 days), Severe stress (20 to 40% water for 30 days) and control treatment (saturated substrate). Measurements by infrared gas analyzer (IRGA) showed a decrease in photosynthesis (up to 70% less than the control), transpiration (up 77%) and stomatal conductance (up 85%) between treatments; furthermore, leaf relative water content dropped 13% between the treatment control and severe, and free proline levels were up to 3.5 fold greater in stress treatments. The leaf temperature increased significantly with increasing irradiance of photons applied. The search for genes related to stress by water deficit in the root transcripts library returned 1,145 sequences, and these, 4 were characterized: TgTPS (trehalose 6-phosphate synthase), TgPIP (aquaporin, protein intrinsic of plasma membrane), TgDREB2 (dehydration responsive element binding protein) and TgAREB (abscisic acid responsive element binding protein). Only TgTPS, TgPIP and TgDREB2 showed a high degree of conservation between species, and can be properly amplified by PCR and validated by sequencing. Thus, with the database of transcripts obtained by RNA-seq, candidate genes were identified for the study of vegetative and reproductive characteristics teak, helping to understand the molecular mechanisms of this forest species.

Precision medicine provides targeted treatment for an individual patient based on disease mechanisms, promoting health care. Matched transcriptomes derived from a single subject enable uncovering patient-specific dynamic changes associated with disease status. Yet, conventional statistical methodologies remain largely unavailable for single-subject transcriptome analysis due to the "single-observation" challenge. We hypothesize that, with statistical learning approaches and large-scale inferences, one can learn useful information from single-subject transcriptome data by identifying differentially expressed genes (DEG) / pathways (DEP) between two transcriptomes of an individual. This dissertation is an ensemble of my research work in single-subject transcriptome analytics, including three projects with varying focuses. The first project describes a two-step approach to identify DEPs by employing a parametric Gaussian mixture model followed by Fisher's exact tests. The second project relaxes the parametric assumption and develops a nonparametric algorithm based on k-means, which is more flexible and robust. The third project proposes a novel variance stabilizing framework to transform raw gene counts before identifying DEGs, and the transformation strategically by-passes the challenge of variance estimation in single-subject transcriptome analysis. In this dissertation, I present the main statistical methods and computational algorithms for all the three projects, as well as their real-data applications to personalized treatments.

Doctor of Philosophy / Department of Plant Pathology / Christopher R. Little / Charcoal rot, caused by the necrotrophic fungus, Macrophomina phaseolina (Tassi) Goid., is an important disease in sorghum (Sorghum bicolor (L.) Moench). The molecular interactions between sorghum and M. phaseolina are poorly understood. In this study, a large-scale RNA-Seq experiment and four follow-up functional experiments were conducted to understand the molecular basis of charcoal rot resistance and/or susceptibility in sorghum. In the first experiment, stalk mRNA was extracted from charcoal-rot-resistant (SC599) and susceptible (Tx7000) genotypes and subjected to RNA sequencing. Upon M. phaseolina inoculation, 8560 genes were differentially expressed between the two genotypes, out of which 2053 were components of 200 known metabolic pathways. Many of these pathways were significantly up-regulated in the susceptible genotype and are thought to contribute to enhanced pathogen nutrition and virulence, impeded host basal immunity, and reactive oxygen (ROS) and nitrogen species (RNS)-mediated host cell death. The paradoxical hormonal regulation observed in pathogen-inoculated Tx7000 was characterized by strongly upregulated salicylic acid and down-regulated jasmonic acid pathways. These findings provided useful insights into induced host susceptibility in response to this necrotrophic fungus at the whole-genome scale. The second experiment was conducted to investigate the dynamics of host oxidative stress under pathogen infection. Results showed M. phaseolina’s ability to significantly increase the ROS and RNS content of two charcoal-rot-susceptible genotypes, Tx7000 and BTx3042. Over-accumulation of nitric oxide (NO) in stalk tissues in the pathogen-inoculated susceptible genotypes was confirmed using a NO-specific fluorescent probe and confocal microscopy. Significantly increased malondialdehyde content confirmed the enhanced oxidative stress experienced by the susceptible genotypes after pathogen inoculation. These findings suggested the contribution of oxidative stress-associated induced cell death on charcoal rot susceptibility under infection. In the third functional experiment, the behavior of the sorghum antioxidant system after pathogen inoculation was investigated. M. phaseolina significantly increased the glutathione s-transferase (GST), glutathione peroxidase (GPX), glutathione reductase (GR), and peroxidase activities of the susceptible genotypes (Tx7000, BTx3042) but not in the resistant genotypes (SC599, SC35). Increased activities of these enzymes in susceptible genotypes may contribute to reduced oxidative stress thus lowering charcoal rot susceptibility. The fourth functional experiment was designed to quantify induced host-derived cell wall degrading enzymes (CWDEs) using crude enzyme mixtures from stalks. A gel diffusion assay revealed significantly increased pectinesterase activity in pathogen-inoculated Tx7000 and BTx3042 while significantly increased polygalacturonase activity was determined by absorbance. Fluorimetric determination of cell extracts revealed significantly increased cellulose degrading enzyme activity in M. phaseolina-inoculated Tx7000 and BTx3042. These findings revealed the pathogen’s ability to promote charcoal rot susceptibility in grain sorghum through induced host CWDEs. The last functional study was designed to profile the stalk tissue lipidome of Tx7000 and SC599 after M. phaseolina inoculation using automated direct infusion electrospray ionization-triple quadrupole mass spectrometry (ESI-MS/MS). M. phaseolina significantly decreased the phytosterol, phosphatidylserine, and ox-lipid contents in Tx7000 while significantly increasing stigmasterol:sitosterol ratio. Except for ox-lipid content, none of the above was significantly affected in resistant SC599. Results suggested the lethal impacts of M. phaseolina inoculation on plastid- and cell- membrane integrity and the lipid-based signaling capacity of Tx7000. Findings shed light on the host lipid classes that contribute to induced charcoal rot susceptibility in grain sorghum.

We study the tasks of transcript expression quantification and differential expression analysis based on data from high-throughput sequencing of the transcriptome (RNA-seq). In an RNA-seq experiment subsequences of nucleotides are sampled from a transcriptome specimen, producing millions of short reads. The reads can be mapped to a reference to determine the set of transcripts from which they were sequenced. We can measure the expression of transcripts in the specimen by determining the amount of reads that were sequenced from individual transcripts. In this thesis we propose a new probabilistic method for inferring the expression of transcripts from RNA-seq data. We use a generative model of the data that can account for read errors, fragment length distribution and non-uniform distribution of reads along transcripts. We apply the Bayesian inference approach, using the Gibbs sampling algorithm to sample from the posterior distribution of transcript expression. Producing the full distribution enables assessment of the uncertainty of the estimated expression levels. We also investigate the use of alternative inference techniques for the transcript expression quantification. We apply a collapsed Variational Bayes algorithm which can provide accurate estimates of mean expression faster than the Gibbs sampling algorithm. Building on the results from transcript expression quantification, we present a new method for the differential expression analysis. Our approach utilizes the full posterior distribution of expression from multiple replicates in order to detect significant changes in abundance between different conditions. The method can be applied to differential expression analysis of both genes and transcripts. We use the newly proposed methods to analyse real RNA-seq data and provide evaluation of their accuracy using synthetic datasets. We demonstrate the advantages of our approach in comparisons with existing alternative approaches for expression quantification and differential expression analysis. The methods are implemented in the BitSeq package, which is freely distributed under an open-source license. Our methods can be accessed and used by other researchers for RNA-seq data analysis.

Orientadores: Felipe Rodrigues da Silva, Francisco Pereira Lobo / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Biologia / Made available in DSpace on 2018-08-27T12:42:35Z (GMT). No. of bitstreams: 1 Castro_GiovanniMarquesde_M.pdf: 8706589 bytes, checksum: bd5b5d6428a549bafa82854367f811b9 (MD5) Previous issue date: 2015 / Resumo: A planta Fevillea cordifolia L. possui um grande potencial para produção de biodiesel. Buscando entender o metabolismo foi realizado um experimento exploratório de RNA-seq com sementes inteiras. No entanto, as análises da qualidade na biblioteca indicaram grande quantidade de sequências virais. Após a reconstrução do transcriptoma usando o programa TRINITY, também foi reconstruído o genoma completo do vírus. O vírus reconstruído possui identidade de 96% com o Papaya ringspot virus (PRSV) em um alinhamento global de nucleotídeos dos genomas completos. Avaliando a abundância dos contigs, o PRSV encontrado representa quase 60% das 24,6 milhões de leituras da biblioteca. Para identificar qual a origem do vírus encontrado, este foi comparado com 29 PRSVs existentes no Genbank através de análise filogenética usando do Algerian watermelon mosaic virus (AWMV) para enraizar a árvore. O vírus encontrado agrupa-se com os dois PRSVs do Brasil, dentro do grupo das Américas estando mais próximo do PRSV-W-C (DQ374152). A existência de recombinações entre os PRSVs foi analisada, porém não foi detectada recombinação recente. Devido à profundidade de sequenciamento maior que 10.000x, foi possível analisar as variações existentes do genoma viral reconstruído. Uma análise das variações dos códons virais foi realizada, mostrando uma tendência para ocorrência de indels para a região do genoma no fim do cístron NIb e início do cístron CP. Essas variações ocorrem devido à existência de haplótipos virais na amostra sequenciada. Para se estimar a diversidade de haplótipos virais da amostra, foi realizada uma reconstrução local da região de 500 nt com mais alta entropia. Foram reconstruídos 58 haplótipos, dos quais dois eram predominantes com frequências de 64,84% e 24,34%. Os haplótipos reconstruídos podem ser usados para o desenvolvimento de resistência mediada por RNAi, evitando que uma variante conhecida e preexistente na população possa quebrar a resistência / Abstract: The plant Fevillea cordifolia L. has a great potential for producing biodiesel. In order to understand its metabolic pathways an exploratory RNA-seq experiment was conducted with its whole-seeds. However the quality analysis of the library revealed a substantial amount of virus sequences. After the reconstruction of the transcriptome using the software TRINITY, the complete viral genome was obtained as well. The reconstructed viral genome had an identity of 96% with Papaya ringspot virus (PRSV) in a global nucleotide alignment using whole genomes. When estimating the abundance of the reconstructed sequences, this PRSV had almost 60% of the 24,6 million reads mapping to it. Aiming to elucidate the origin of this virus, its sequence was compared to 29 PRSVs from Genbank using phylogenetic analysis and the Algerian Watermelon Mosaic Virus (AWMV) as an outgroup. This PRSV clustered with the Brazilian isolates, being closer to PRSV-W-C (DQ374152). A recombination analysis was performed within the PRSVs but no recent recombination was detected. Due to the depth of the coverage sequencing being higher than 10.000x, it was possible to analyze the variations existing in the reconstructed genome. An analysis of the codons variations was performed, revealing a tendency for the occurrence of indels in a region at the end of NIb cistron and at the start of the CP cistron. These variations occur due to the existence of viral haplotypes sequenced in this sample. A local reconstruction of the 500nt region with the highest entropy was performed to estimate the diversity of viral haplotypes in this sample. 58 haplotypes were reconstructed, of which 2 were dominant with frequencies of 64,84% and 24,34%. The reconstructed haplotypes may be used for the development of RNAi-mediated resistance, avoiding the breaking of the resistance by variants that are known to exist in the population / Mestrado / Bioinformatica / Mestre em Genética e Biologia Molecular

No description available.

L’analyse des défauts d’épissage est particulièrement complexe. Outre la diversité des transcrits présents à l’état physiologique, les variations nucléotidiques peuvent induire des modifications hétéroclites de l’épissage. Ces variations, appelées variants splicéogéniques, et leur impact au niveau de l’épissage, sont à même de modifier plus ou moins sévèrement le phénotype de l’individu.Au cours de ce travail de thèse, nous nous sommes intéressés à trois grands aspects de l’étude des défauts de l’épissage : (i) la prédiction de ces défauts d’épissage, (ii) l’analyse des données de RNA-seq et (iii) le rôle de l’épissage dans l’interprétation de la pathogénicité d’un variant pour la prédisposition aux cancers du sein et de l’ovaire (syndrome HBOC).Nous avons optimisé les recommandations en vigueur pour identifier les variants splicéogéniques au sein des séquences consensus des sites d’épissage. Ce travail a conduit à la publication d’un nouvel outil SPiCE (Splicing Prediction in Consensus Elements), développé sur 395 variants. SPiCE a le potentiel d’être une aide à la décision pour guider les généticiens vers ces variants splicéogéniques, grâce à une exactitude de 94.4 %. Puis, nous avons comparé les outils de prédiction des points de branchement. Pour cela, une collection sans précédente de 120 variants avec leurs études ARN a été établi dans la région des points de branchements. Nous avons ainsi révélé que ces outils de prédictions sont aptes à prioriser les variants pour des études ARN dans ces régions jusque-là peu étudiées. Pour étendre les prédictions des variants splicéogéniques au-delà d’un motif spécifique, nous avons construit l’outil SPiP (Splicing Prediction Pipeline). SPiP utilise un ensemble d’outils pour prédire un défaut d’épissage quel que soit la position du variant. Ainsi, SPiP peut ainsi s’adresser à la diversité des défauts d’épissage avec une exactitude de 80.21 %, sur une collection de 2 784 variants.Les données issues du RNA-seq sont complexes à analyser, car il existe peu d’outils pour annoter finement les épissages alternatifs. Aussi nous avons publié l’outil SpliceLauncher. Cet outil permet de déterminer une grande diversité de jonctions d’épissage, indépendamment des systèmes RNA-seq utilisés. Cet outil renvoie aussi les résultats sous formes graphiques pour faciliter leur interprétation.Puis nous avons évalué le rôle de l’épissage alternative dans l’interprétation à usage clinique d’un variant. Le gène PALB2, impliqué dans le syndrome HBOC, a été utilisé comme modèle d’étude. Nous avons ainsi démontré que l’épissage alternatif de PALB2 est apte à remettre en cause la pathogénicité de certains variants. La collecte de données fonctionnelles et cliniques sont donc nécessaires pour conclure sur leur pathogénicité.Nos travaux illustrent ainsi l’importance de la caractérisation et de l’interprétation des modifications de l’épissage pour répondre aux défis présents et futurs du diagnostic moléculaire en génétique. / Analysis of splicing defects is particularly complex. In addition to the diversity of physiological transcripts, nucleotidic variations can induce heterogeneous alteration of splicing. These variations, called spliceogenic variants, and their impact on splicing, can involve severe consequences on the individual phenotype.In this thesis work, we focused on three main aspects of the study of splicing defects: (i) the prediction of these splicing defects, (ii) the analysis of RNA-seq data and (iii) the role of splicing in interpreting the pathogenicity of a variant for the hereditary breast and ovarian cancers (HBOC syndrome).We optimized the current recommendations to identify spliceogenic variants within the consensus sequences of splicing sites. This work led to the publication of a new tool, SPiCE (Splicing Prediction in Consensus Elements), developed on 395 variants. SPiCE has the potential to be a decision support tool to guide geneticists towards these spliceogenic variants, with an accuracy of 94.4%. Then, we compared the tools dedicated to branch points prediction. For this purpose, an unprecedented collection of 120 variants with their RNA studies has been established in the branch point region. Thus, we revealed these prediction tools are able to prioritize variants for RNA studies in these hitherto poorly studied regions. To extend the predictions of spliceogenic variants beyond a specific motif, we built SPiP (Splicing Prediction Pipeline) tool. SPiP uses a set of tools to predict a splicing defect regardless of the variant position. Thus, SPiP can address the diversity of splicing defects with an accuracy of 80.21%, on a collection of 2,784 variants.The data from the RNA-seq are complex to analyze, as there are few tools to finely annotate alternative splices. Also we published SpliceLauncher tool. This tool allows to determine a wide variety of splicing junctions, independently of RNA-seq systems used. This tool also returns the results in graphical form to make interpretation user-friendly.Then we evaluated the role of alternative splicing in the clinical interpretation of a variant. The PALB2 gene, involved in HBOC syndrome, was used as a study model. Thus, we demonstrated that the alternative splicing of PALB2 is able of challenging the pathogenicity of certain variants. Collection of functional and clinical data is therefore necessary to conclude on their pathogenicity.Our work thus illustrates the importance of characterizing and interpreting splicing modifications to meet the current and future challenges of molecular diagnosis in human genetics.

Urinary tract infections (UTI) are among the most common infections in humans. Uropathogenic Escherichia coli (UPEC), the main causative agent of UTIs, can invade and replicate within bladder epithelial cells, and recent evidence demonstrated that some UPEC strains also survive within macrophages. To understand the mechanisms of host subversion that enable UPEC to survive within macrophages, and the contribution of macrophages to UPEC-mediated pathology, I performed host-pathogen co-transcriptome analyses using RNA sequencing. I developed an effective computational framework that simultaneously separated, annotated, and quantified the mammalian and bacterial transcriptomes. First, mouse bone morrow-derived macrophages (BMM) were challenged over a 24 h time course with UPEC reference strains, UTI89 (cystitis strain), 83972 and VR50 (asymptomatic bacteriuria strains) that possess contrasting intramacrophage phenotypes. My results showed that BMM responded to the three different UPEC strains with broadly similar gene expression programs. In contrast to the conserved pattern of BMM responses, the transcriptional responses of the different UPEC strains diverged markedly from each other. Hypothesizing that genes upregulated at 24 h post-infection may contribute to intramacrophage survival, I identified UTI89 genes upregulated at this time point, and showed that deletion of one of these genes (pspA) compromised intramacrophage survival of UPEC strain UTI89. Second, human monocyte-derived macrophages (HMDM) and BMM were challenged over a 24 h course with the UPEC strain EC958, a globally disseminated, multi-drug resistant strain. My analysis identified extensive divergence in UPEC-regulated orthologous gene expression between HMDM and BMM, and I validated both known and novel genes in the context of differential regulation. On the contrary, the transcriptional response of EC958 showed a broad conservation across both mammalian intramacrophage environments. My study thus provides both a unique co-culture approach to study infection in vitro and a technological framework for simultaneously capturing global changes in host-pathogen interactions at the transcriptional level in co-cultures. In conclusion, this work has generated new insights into the mechanisms that UPEC strains exploit to persist within the mouse intramacrophage environment, as well as differences in the transcriptional repertoire of HMDM and BMM challenged with the same UPEC strain.

Diplomová práce se zabývá studiem genové regulace v Clostridium beijerinckii NRRL B-598, pro následné odvození genové regulační sítě bakterie C. beijerinckii NRRL B-598. V teoretické části této práce je uvedena obecná nomenklatura problematiky genové regulace se zaměřením na nomenklaturu genových regulačních sítí. Následně jsou zde popsané laboratorní metody, sloužící pro získání vhodných dat popisující expresi genů. Tato data jsou základem pro studium genové regulace a návrhy genových regulačních sítí. Práce se zaměřuje především na technologii RNA-Seq a stručný popis laboratorních dat získaných ze zmíněné bakterie C. beijerinckii NRRL B-598. V praktické části se práce zabývá předzpracováním těchto surových laboratorních dat a následným studiem genové regulace se zaměřením na odvození operonů a vytvoření prvních genových regulačních sítí pomocí různých přístupů pro C. beijerinckii NRRL B-598.

No description available.