Cell Fate Specification and the Regulation of RNA-dependent DNA Methylation in the Arabidopsis Root Meristem

The Arabidopsis root apical meristem (RAM) is a complex tissue capable of generating all the cell types that ultimately make up the root. The work presented in this thesis takes advantage of the versatility of high-throughput sequencing to address two independent questions about the root meristem. Although a lot of information is known regarding the cell fate decisions that occur at the RAM, cortex specification and differentiation remain poorly understood. In the first part of this thesis, I used an ethylmethanesulfonate (EMS) mutagenized marker line to perform a forward genetics screen. The goal of this screen was to identify novel genes involved in the specification and differentiation of the cortex tissue. Mapping analysis from the results obtained in this screen revealed a new allele of BRASSINOSTEROID4 with abnormal marker expression in the cortex tissue. Although this allele proved to be non-cortex specific, this project highlights new technology that allows mapping of EMS-generated mutations without the need to map-cross or back-cross. In the second part of this thesis, using fluorescence activated cell sorting (FACS) coupled with high throughput sequencing, my collaborators and I generated single-base resolution whole genome DNA methylomes, mRNA transcriptomes, and smallRNA transcriptomes for six different populations of cell types in the Arabidopsis root meristem. We were able to discover that the columella is hypermethylated in the CHH context within transposable elements. This hypermethylation is accompanied by upregulation of the RNA-dependent DNA methylation pathway (RdDM), including higher levels of 24-nt silencing RNAs (siRNAs). In summary, our studies demonstrate the versatility of high-throughput sequencing as a method for identifying single mutations or to perform complex comparative genomic analyses.

Role and modulation of maternal transcripts during the first cleavage divisions in bovine embryos

Ce travail porte sur l’identification, la fonction et la régulation des molécules maternelles d’ARNm qui dirigent la compétence développementale juste après la fécondation chez les bovins. Tout d’abord, en utilisant le modèle du temps écoulé jusqu’au premier clivage zygotique et à travers l’évaluation du transcriptome des embryons à 2-cellules, il fut possible de déterminer la signature moléculaire des niveaux extrêmes de compétence au développement et sélectionner des molécules candidates pour des études postérieures. Les résultats ont montré que les embryons de capacité développementale variable diffèrent dans certaines fonctions comme la réparation de l’ADN, le traitement de l’ARN, la synthèse de protéines et l’expression génique définies par des ARNm synthétisés par l’ovocyte. Pour obtenir une confirmation fonctionnelle, une paire de transcrits maternels (l’un détecté dans notre sondage précédent et l’autre étant une molécule reliée) ont été inhibés par « knock-down » dans des ovocytes. Les effets du knock-down de ces facteurs de transcription sont apparus avant la formation des blastocystes dû à une diminution de la capacité au clivage et celle à progresser après le stage de 8-cellules. L’analyse moléculaire des embryons knock-down survivants suggère qu’un de ces facteurs de transcription est un contrôleur crucial de l’activation du génome embryonnaire, qui représente une fenêtre développementale dans l’embryogenèse précoce. Dans la dernièr étude, nous avons testé si les facteurs de transcription d’intérêt sont modulés au niveau traductionnel. Des ARNm rapporteurs couplés à la GFP (Protéine fluorescente) contenant soit la version courte ou la version longue de la séquence 3’-UTR des deux molécules furent injectées dans des zygotes pour évaluer leur dynamique traductionnelle. Les résultats ont montré que les éléments cis-régulateurs localisés dans les 3’-UTRs contrôlent leur synchronisation traductionnelle et suggèrent une association entre la compétence développementale et la capacité de synthèse de ces protéines. Ceci conduit à l’idée que ces facteurs de transcription cruciaux sont aussi contrôlés au niveau traductionnel chez les embryons précoces. Les connaissances acquises ont joué un rôle essentiel pour définir le contrôle potentiel des molécules maternelles sur les embryons au début de leur développement. Cette étude nous montre aussi une utilisation potentielle de cette information ainsi que les nouveaux défis présents dans le secteur des technologies reproductives.

Phylum-wide transcriptome analysis of oogenesis and early embryogenesis in selected nematode species

Oogenesis is a prerequisite for embryogenesis in Metazoa. During both biological processes important decisions must be made to form the embryo and hence ensure the next generation: (1) Maternal gene products (mRNAs, proteins and nutrients) must be supplied to the embryo. (2) Polarity must be established and axes must be specified. While incorporation of maternal gene products occurs during oogenesis, the time point of polarity establishment and axis specification varies among species, as it is accomplished either prior, during, or after fertilisation. But not only the time point when these events take place varies among species but also the underlying mechanisms by which they are triggered. For the nematode model Caenorhabditis elegans the underlying pathways and gene regulatory networks (GRNs) are well understood. It is known that there the sperm entry point initiates a primary polarity in the 1-celled egg and with it the establishment of the anteroposterior axis. However, studies of other nematodes demonstrated that polarity establishment can be independent of sperm entry (Goldstein et al., 1998; Lahl et al., 2006) and that cleavage patterns, symmetry formation and cell specification also differ from C. elegans. In contrast to the studied Chromadorea (more derived nematodes including C. elegans), embryos of some marine Enoplea (more basal representatives) even show no discernible early polarity and blastomeres can adopt variable cell fates (Voronov and Panchin 1998). The underlying pathways controlling the obviously variant embryonic processes in non-Caenorhabditis nematodes are essentially unknown. In this thesis I addressed this issue by performing a detailed unbiased comparative transcriptome analysis based on microarrays and RNA sequencing of selected developmental stages in a variety of nematodes from different phylogenetic branches with C. elegans as a reference system and a nematomorph as an outgroup representative. In addition, I made use of available genomic data to determine the presence or absence of genes for which no expression had been detected. In particular, I focussed on components of selected pathways or GRNs which are known to play essential roles during C. elegans development and/or other invertebrate or vertebrate model systems. Oogenesis must be regulated differently in non-Caenorhabditis nematodes, as crucial controlling components of Wnt and sex determination signaling are absent in these species. In this respect, I identified female-specific expression of potential polarity associated genes during gonad development and oogenesis in the Enoplean nematode Romanomermis culicivorax. I could show that known downstream components of the polarity complexes PAR-3/-6/PKC-3 and PAR-1/-2 are absent in non-Caenorhabditis species. Even PAR-2 as part of the polarity complex does not exist in these nematodes. Instead, transcriptomes of nematodes (including C. elegans), show expression of other polarity-associated complexes such as the Lgl (Lethal giant larvae) complex. This result could pose an alternative route for nematodes and nematomorphs to initiate polarity during early embryogenesis. I could show that crucial pathways of axis specification, such as Wnt and BMP are very different in C. elegans compared to other nematodes. In the former, Wnt signaling, for instance, is mediated by four paralogous beta-catenins, while other Chromadorea have fewer and Enoplea only one beta-catenin. The transcriptomes of R. culicivorax and the nematomorph show that regulators of BMP (e.g. Chordin), are specifically expressed during early embryogenesis only in Enoplea and the close outgroup of nematomorphs. In conclusion, my results demonstrate that the molecular machinery controlling oogenesis and embryogenesis in nematodes is unexpectedly variable and C. elegans cannot be taken as a general model for nematode development. Under this perspective, Enoplean nematodes show more similarities with outgroups than with C. elegans. It appears that certain pathway components were lost or gained during evolution and others adopted new functions. Based on my findings I can conjecture, which pathway components may be ancestral and which were newly acquired in the course of nematode evolution.

Analyse transcriptomique de l'interaction tripartite "Pseudozyma flocculosa-Blumeria graminis f.sp. hordei-Hordeum vulgare"

Afin d’améliorer nos pratiques agricoles dans le contexte d’une agriculture durable, plusieurs agents de lutte biologique (ALB) ont été développés, testés et sont maintenant utilisés dans le monde pour combattre les pertes de rendements causées par les maladies. Blumeria graminis f. sp. hordei ( Bgh) est l’agent pathogène responsable du blanc de l’orge et peut réduire les rendements de cette culture jusqu’à 40%. Un champignon épiphyte, Pseudozyma flocculosa, a été découvert et identifié en 1987 en association étroite avec le blanc du trèfle. Les chercheurs ont alors remarqué que ce champignon exhibait une forte activité antagoniste contre le blanc en détruisant les structures de l’agent pathogène. Suite à d’autres travaux, il est apparu que ce comportement antagoniste était dirigé contre tous les membres des Erysiphales et semblait lié à la synthèse d’un glycolipide antifongique soit la flocculosine. Toutefois, on n’est toujours pas parvenus à associer l’efficacité de l’ALB avec la production de ce glycolipide. Ces observations suggèrent que d’autres facteurs seraient impliqués lorsque les deux protagonistes, l’ALB et le blanc, sont en contact. L’objectif principal de ce projet était donc de chercher d’autres mécanismes moléculaires pouvant expliquer l’interaction P. flocculosa-blanc et orge, en faisant une analyse transcriptomique complète des trois protagonistes en même temps. L’interaction tripartite a été échantillonnée à différents temps suivant l’inoculation de P. flocculosa sur des feuilles d’orge présentant déjà une intensité de blanc d’environ 50%. Les échantillons de feuilles prélevés ont ensuite été utilisés pour l’extraction de l’ARN qui ont été ensuite transformés en ADNc pour la préparation des librairies. Cinq répliquats ont été effectués pour chaque temps et le tout a été séquencé à l’aide de séquençage par synthèse Illumina HiSeq. Les séquences obtenues (reads) ont ensuite été analysées à l’aide du logiciel CLC Genomics Workbench. Brièvement, les séquences obtenues ont été cartographiées sur les trois génomes de référence. Suite à la cartographie, les analyses d’expression ont été conduites et les gènes exprimés de façon différentielle ont été recherchés. Cette étape a été conduite en portant une attention particulière aux gènes codant pour un groupe de protéines appelées CSEP pour “candidate secreted effector proteins” qui seraient possiblement impliquées dans l’interaction tripartite. Parmi les protéines exprimées de façon différentielle en présence du blanc ou en absence de ce dernier, nous avons pu constater que certaines CSEP étaient fortement exprimées en présence du blanc. Ces résultats sont prometteurs et nous offrent une piste certaine pour l’élucidation des mécanismes impliqués dans cette interaction tripartite.

Doença enxerto-versus hospedeiro: determinação da expressão de genes e alterações de vias associadas ao sistema imunológico usando metodos da biofísica

O transplante de medula óssea (TMO) é um procedimento terapêutico importante em casos relacionados à pacientes com leucemia ou linfoma. Em decorrência desse processo, uma reação conhecida como doença enxerto-versus-hospedeiro (GVHD) pode ocorrer em pacientes susceptíveis como conseqüência da presença de células imunocompetentes do doador. Entretanto, não existe um modelo para descrever completamente as ações relacionadas ao mecanismo imunológico da GVHD desde a fase que inicializa a doença até a fase efetora. O Objetivo geral deste estudo é a investigação da resposta imunológica considerando-se o sistema HLA (antígenos leucocitários humano) em pacientes que desenvolveram a GVHD em decorrência do TMO. O National Cancer Institute (NCI) – Pathway interaction Database e Reactome foram usados como bases de dados com o objetivo de se estudar a expressão de genes e vias relacionados às Classes I e II do sistema HLA (antígenos leucocitários humano). O estudo considerou a mudança de expressão de genes relacionados às 17 vias do sistema imunológico com potencialidade para se expressar em pacientes que desenvolveram a GVHD associada à TMO. Dados referentes aos transcriptomas foram obtidos utilizando-se a plataforma GPL570 Affymetrix Genoma Humano U133 Plus. A atividade relativa foi usada para determinar as alterações das vias em amostras de GVHD em relação ao controle. As análises foram realizadas utilizando-se o software Via Complex e Bioconductor. Observou-se aumento significativo da expressão de genes ralacionados às vias do sistema imune adaptativo, antígenos associados às Classe I e II do HLA, fosforilação de CD3 e CD247, sinalização dos receptores de células T em CD4+ nativas e ativação de NF-kapa β nas células B. Também observou-se alterações significativas na mudança de expressão dos genes associados às vias relacionadas à super família de moléculas B7:CD28\CTLA-4 quando comparadas ao controle. Isso pode indicar a necessidade de geração de um segundo sinal co-estimulador em GVHD, acionado pelas moléculas dessa super família. O aumento da expressão do gene CD69 nas amostras experimentais caracteriza a ativação celular e, portanto, a sinalização de estímulos em GVHD. Os achados obtidos neste estudo contribuem para melhor elucidar o mecanismo imunopatogênico associado à GVHD. P

Transcriptome Analysis of Capsicum Chlorosis Virus-Induced Hypersensitive Resistance Response in Bell Capsicum

Background Capsicum chlorosis virus (CaCV) is an emerging pathogen of capsicum, tomato and peanut crops in Australia and South-East Asia. Commercial capsicum cultivars with CaCV resistance are not yet available, but CaCV resistance identified in Capsicum chinense is being introgressed into commercial Bell capsicum. However, our knowledge of the molecular mechanisms leading to the resistance response to CaCV infection is limited. Therefore, transcriptome and expression profiling data provide an important resource to better understand CaCV resistance mechanisms. Methodology/Principal Findings We assembled capsicum transcriptomes and analysed gene expression using Illumina HiSeq platform combined with a tag-based digital gene expression system. Total RNA extracted from CaCV/mock inoculated CaCV resistant (R) and susceptible (S) capsicum at the time point when R line showed a strong hypersensitive response to CaCV infection was used in transcriptome assembly. Gene expression profiles of R and S capsicum in CaCV- and buffer-inoculated conditions were compared. None of the genes were differentially expressed (DE) between R and S cultivars when mock-inoculated, while 2484 genes were DE when inoculated with CaCV. Functional classification revealed that the most highly up-regulated DE genes in R capsicum included pathogenesis-related genes, cell death-associated genes, genes associated with hormone-mediated signalling pathways and genes encoding enzymes involved in synthesis of defense-related secondary metabolites. We selected 15 genes to confirm DE expression levels by real-time quantitative PCR. Conclusion/Significance DE transcript profiling data provided comprehensive gene expression information to gain an understanding of the underlying CaCV resistance mechanisms. Further, we identified candidate CaCV resistance genes in the CaCV-resistant C. annuum x C. chinense breeding line. This knowledge will be useful in future for fine mapping of the CaCV resistance locus and potential genetic engineering of resistance into CaCV-susceptible crops.

The Kinome of Pacific Oyster Crassostrea gigas, Its Expression during Development and in Response to Environmental Factors

Oysters play an important role in estuarine and coastal marine habitats, where the majority of humans live. In these ecosystems, environmental degradation is substantial, and oysters must cope with highly dynamic and stressful environmental constraints during their lives in the intertidal zone. The availability of the genome sequence of the Pacific oyster Crassostrea gigas represents a unique opportunity for a comprehensive assessment of the signal transduction pathways that the species has developed to deal with this unique habitat. We performed an in silico analysis to identify, annotate and classify protein kinases in C. gigas, according to their kinase domain taxonomy classification, and compared with kinome already described in other animal species. The C. gigas kinome consists of 371 protein kinases, making it closely related to the sea urchin kinome, which has 353 protein kinases. The absence of gene redundancy in some groups of the C. gigas kinome may simplify functional studies of protein kinases. Through data mining of transcriptomes in C. gigas, we identified part of the kinome which may be central during development and may play a role in response to various environmental factors. Overall, this work contributes to a better understanding of key sensing pathways that may be central for adaptation to a highly dynamic marine environment.

Dinoflagellate Genomic Organization and Phylogenetic Marker Discovery Utilizing Deep Sequencing Data

Dinoflagellates possess large genomes in which most genes are present in many copies. This has made studies of their genomic organization and phylogenetics challenging. Recent advances in sequencing technology have made deep sequencing of dinoflagellate transcriptomes feasible. This dissertation investigates the genomic organization of dinoflagellates to better understand the challenges of assembling dinoflagellate transcriptomic and genomic data from short read sequencing methods, and develops new techniques that utilize deep sequencing data to identify orthologous genes across a diverse set of taxa. To better understand the genomic organization of dinoflagellates, a genomic cosmid clone of the tandemly repeated gene Alchohol Dehydrogenase (AHD) was sequenced and analyzed. The organization of this clone was found to be counter to prevailing hypotheses of genomic organization in dinoflagellates. Further, a new non-canonical splicing motif was described that could greatly improve the automated modeling and annotation of genomic data. A custom phylogenetic marker discovery pipeline, incorporating methods that leverage the statistical power of large data sets was written. A case study on Stramenopiles was undertaken to test the utility in resolving relationships between known groups as well as the phylogenetic affinity of seven unknown taxa. The pipeline generated a set of 373 genes useful as phylogenetic markers that successfully resolved relationships among the major groups of Stramenopiles, and placed all unknown taxa on the tree with strong bootstrap support. This pipeline was then used to discover 668 genes useful as phylogenetic markers in dinoflagellates. Phylogenetic analysis of 58 dinoflagellates, using this set of markers, produced a phylogeny with good support of all branches. The Suessiales were found to be sister to the Peridinales. The Prorocentrales formed a monophyletic group with the Dinophysiales that was sister to the Gonyaulacales. The Gymnodinales was found to be paraphyletic, forming three monophyletic groups. While this pipeline was used to find phylogenetic markers, it will likely also be useful for finding orthologs of interest for other purposes, for the discovery of horizontally transferred genes, and for the separation of sequences in metagenomic data sets.

Allocation du carbone et métabolisme azoté chez l'haptophyte Tisochrysis lutea

The aim of this thesis is to improve knowledge on mechanisms involved in the response to nitrogen limitation and in lipid accumulation in the microalgae haptophyte Tisochrysis lutea. The wild type strain and a lipid accumulating mutant strain were grown under different nitrogen limitation and starvation and analyzed by functional genomics. Four genes of high-affinity nitrate/nitrite transporter (Nrt2) were identified and characterized to reveal the mechanisms involved in mineral absorption in this species. Transcriptomes of both strains were sequenced and proteins affected by nitrogen starvation and differentially expressed between the two strains were identified. We so identified the functions regulated by nitrogen deficiency and potentially involved in the accumulation of storage lipids. The responses of both strains to thin variations of nitrogen limitation were studied. The results of high-throughput proteomic analyzes suggest that the lipid-accumulation in the mutant strain is the result of carbon metabolism impacted overall, this spurred on signaling mechanisms. Two proteins have been studied since probably involved in carbon and nitrogen remobilization from amino acids catabolism during nitrogen limitation. This work increases knowledge on haptophytes, and brings assumptions on metabolic key involved in nitrogen limitation and carbon allocation in microalgae.

A Transcriptome Atlas of Physcomitrella patens Provides Insights into the Evolution and Development of Land Plants

Post-print version of the article.

Novel ciliate genetic code variants including the reassignment of all three stop codons to sense codons in Condylostoma magnum

mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum. We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3' end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum. We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3' end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.