This thesis was performed within the research project “TrophinOak”, which addresses the impact of multitrophic interactions on the pedunculate oak (Quercus robur) clone DF159. In this frame, the present work focuses on the genetic and physiological mechanisms ruling the interaction of the mycorrhiza helper bacterium (MHB) Streptomyces sp. AcH 505 with microcuttings of DF159 either alone or in presence of the ectomycorrhizal fungus Piloderma croceum or the fungal leaf pathogen oak powdery mildew Microsphaera alphitoides. The work consists of 3 chapters. Chapter 1 characterises the growth of AcH 505 and P. croceum in a soil-based culture system used within the TrophinOak project. Besides the establishment and evaluation of quantification methods of these microorganisms by quantitative real-time PCR, the impact of the soil microbial community and the oak on the bacterium-fungus interaction was investigated, and AcH 505 and P. croceum were visualized by scanning electron microscopy. It was observed that the presence of the soil microorganisms and the oak both affect the bacterium-fungus interaction, and that P. croceum enhances the growth of AcH 505. Chapter 2 presents a study with the oak, AcH 505 and the EM fungus P. croceum, enabling to disentangle the direct effect of the MHB on the oak from the indirect one via the EM symbiosis. The used approach was transcriptomic based on RNA sequencing. It was shown that i) differential gene expression occurred between root and the distant leaf tissues (local vs. systemic effects), different developmental stages and treatments, suggesting that oak specifically coordinates its gene expression patterns, and ii) that genes related to plant growth, defence and DNA modification were dominant among the differential expressed genes, suggesting that these processes play essential roles in both symbiotic interactions investigated. Chapter 3 represents a second transcriptome study, addressing how AcH 505 suppresses powdery mildew infection in oak by analysing RNA Sequencing data from singly- and coinoculated oaks. This study combined the systemic impact of the root associated bacterium with local effects of the leaf pathogen, thereby linking belowground and aboveground interactions. Systemic defence response is induced by the bacterium and further enhanced upon pathogen challenge, suggesting that on the leaf level, some bacterial effectors are recognized as harmful for the plant.

Gene expression regulation is dynamic and specific to various factors such as developmental stages, environmental conditions, and stimulation of pathogens. Nowadays, a tremendous amount of transcriptome data sets are available from diverse species. This trend enables us to perform comparative transcriptome analysis that identifies conserved or diverged gene expression responses across species using transcriptome data. The goal of this dissertation is to develop and apply approaches of comparative transcriptomics to transfer knowledge from model species to non-model species with the hope that such an approach can contribute to the improvement of crop yield and human health. First, we presented a comprehensive method to identify cross-species modules between two plant species. We adapted the unsupervised network-based module finding method to identify conserved patterns of co-expression and functional conservation between Arabidopsis, a model species, and soybean, a crop species. Second, we compared drought-responsive genes across Arabidopsis, soybean, rice, corn, and Populus in order to explore the genomic characteristics that are conserved under drought stress across species. We identified hundreds of common gene families and conserved regulatory motifs between monocots and dicots. We also presented a BLS-based clustering method which takes into account evolutionary relationships among species to identify conserved co-expression genes. Last, we analyzed single-cell RNA-seq data from monocytes to attempt to understand regulatory mechanism of innate immune system under low-grade inflammation. We identified novel subpopulations of cells treated with lipopolysaccharide (LPS), that show distinct expression patterns from pro-inflammatory genes. The data revealed that a promising therapeutic reagent, sodium 4-phenylbutyrate, masked the effect of LPS. We inferred the existence of specific cellular transitions under different treatments and prioritized important motifs that modulate the transitions using feature selection by a random forest method. There has been a transition in genomics research from bulk RNA-seq to single-cell RNA-seq, and scRNA-seq has become a widely used approach for transcriptome analysis. With the experience we gained by analyzing scRNA-seq data, we plan to conduct comparative single-cell transcriptome analysis across multiple species. / Doctor of Philosophy / All cells in an organism have the same set of genes, but there are different cell types, tissues, organs with different functions as the organism ages or under different conditions. Gene expression regulation is one mechanism that modulates complex, dynamic, and specific changes in tissues or cell types for any living organisms. Understanding gene regulation is of fundamental importance in biology. With the rapid advancement of sequencing technologies, there is a tremendous amount of gene expression data (transcriptome) from individual species in public repositories. However, major studies have been reported from several model species and research on non-model species have relied on comparison results with a few model species. Comparative transcriptome analysis across species will help us to transform knowledge from model species to non-model species and such knowledge transfer can contribute to the improvement of crop yields and human health. The focus of my dissertation is to develop and apply approaches for comparative transcriptome analysis that can help us better understand what makes each species unique or special, and what kinds of common functions across species have been passed down from ancestors (evolutionarily conserved functions). Three research chapters are presented in this dissertation. First, we developed a method to identify groups of genes that are commonly co-expressed in two species. We chose seed development data from soybean with the hope to contribute to crop improvement. Second, we compared gene expression data across five plant species including soybean, rice, and corn to provide new perspectives about crop plants. We chose drought stress to identify conserved functions and regulatory factors across species since drought stress is one of the major stresses that negatively impact agricultural production. We also proposed a method that groups genes with evolutionary relationships from an unlimited number of species. Third, we analyzed single-cell RNA-seq data from mouse monocytes to understand the regulatory mechanism of the innate immune system under low-grade inflammation. We observed how innate immune cells respond to inflammation that could cause no symptoms but persist for a long period of time. Also, we reported an effect of a promising therapeutic reagent (sodium 4-phenylbutyrate) on chronic inflammatory diseases. The third project will be extended to comparative single-cell transcriptome analysis with multiple species.

No description available.

Le génôme des cellules eucaryotes est condensé au sein d'une structure complexe hiérarchiquement organisée : la chromatine. La chromatine est composée d'ADN, de protéines histone et non-histone. Cette thèse a pour but d'étudier le rôle de la chromatine dans la réponse cellulaire aux dommages de l'ADN (DDR) par les méthodologies de génomique fonctionnelle et de protéomique. Nous avons tout d'abord analysé les modifications post-traductionnelles (PTM) des histones dans le cadre des pontages inter-brins (ou "Interstrand Crosslinks", ICL), type particulier de lésions de l'ADN, en choisissant le modèle de l'Anémie de Fanconi (FA). Ceci a été réalisé grâce aux techniques de protéomique quantitative SILAC (Stable Isotope Labeling of Amino acid during Cell culture) et de spectrométrie de masse (MS). Nous avons ainsi réussi à identifier et à quantifier de nombreuses PTMs dans les histones H3 et H4, et à démontrer que certaines de ces PTM sont dépendantes d'une voie fonctionnelle de la signalisation de FA. Nous avons également approfondi l'étude des DDR dans les cellules de FA par une approche de génomique fonctionnelle. Pour cela, nous avons analysé le profil l'expression d'enzymes associées à l'acétylation et à la méthylation des histones. Nos résultats suggèrent l'existence de corrélations entre le profil d'expression de ces enzymes et les PTMs des histones. Des études complémentaires sont nécessaires en vue de confirmer ces corrélations. Nous avons également comparé le transcriptome de deux lignées cellulaires de FA (mutée en FANCC et corrigée en FANCC) après induction de dommages à l'ADN. Afin de différencier les changements spécifiquement associés à la voie de signalisation de FA en réponse aux ICL de l'ADN des réponses plus générales aux dommages de l'ADN, nous avons inclus des cellules traitées par rayonnement ionisant. En réalisant une analyse d'interactions factorielles, nous avons pu identifier une réponse transcriptionnelle aux dommages de l'ADN nécessitant une voie fonctionnelle de la signalisation de FA. Nous avons également tenté de pallier aux limitations rencontrées dans l'analyse des PTMs des histones. En effet, les PTMs des histones que nous avons identifiées représentent l'ensemble des modifications, c'est-à-dire les PTMs concernant les histones se trouvant immédiatement à proximité du site du dommage et en relation directe avec celui-ci, et les PTMs se trouvant à distance du dommage et pouvant ne pas être en relation directe avec celui-ci. Les approches courantes pour identifier les PTMs se trouvant à des loci particuliers sont basées sur l'immunoprécipitation classique de la chromatine où l'utilisation de formaldéhyde altère les protéines, ce qui en rend impossible l'analyse par MS. Nous avons proposé une nouvelle méthodologie basée sur la biotinylation expérimentale d'histones situées à proximité d'une protéine particulière, suivie de la purification des nucléosomes contenant ces histones biotinylées. Contrairement àl'immunoprécipiatation classique de la chromatine, cette méthode n'induit pas d'altération des protéines, permettant ainsi de purifier les histones à partir d'un locus spécifique et d'analyser à grande échelle leurs PTMs par MS. Cette approche permet aussi de suivre dans le temps les PTMs d'une fraction des histones juste après leur biotinylation. Enfin, elle présente l'avantage de pouvoir étudier le profil des PTMs de différents états fonctionnels de la chromatine grâce à l'utilisation de variants d'histones. / In eukaryotic cells, the genome is packed into chromatin, a hierarchically organized complex composed of DNA and histone and nonhistone proteins. In this thesis we have addressed the role of chromatin in cellular response to DNA damage (DDR) using various methodologies encompassing functional genomics and proteomics. First, we analyzed histone post-translational modifications (PTM) in the context of specific kind of DNA lesions (ICL-Interstrand Crosslinks) in Fanconi anemia using quantitative proteomics methodology, SILAC (Stable Isotope Labeling of Amino acids during Cell Culture). Using mass spectrometry (MS), we have successfully identified and quantified a number of histone PTM marks in histone H3 and H4, mainly acetylations and methylations,which have shown dependence upon functional FA-pathway. As a next step, we applied a functional genomics approach to study DDR in FA cells. In this analysis we first monitored the expression profile of histone modifying enzymes related to histone acetylations and methylations. Our results suggest some correlations between histone PTMs and gene expression of histone modifying enzymes, although conclusive evidence warrants further investigations. Next, we analyzed the total transcriptome after DNA damage induction in FA mutant and wild type cells. We also included in this analysis IR irradiation, in an attempt to dissociate more generic DDR from more specific changes that are associated with the role of FA pathway to the DNA ICLs. By performing a factorial interaction analysis, we were able to isolate the part of transcriptional response to DNA damage that was requiring functional FA pathway, as well as the genes that were sensitized to DNA damage by the inactivation of FA pathway. In the final part of the thesis, we attempted to solve one of the limitations that we encountered in the histone PTM analysis. The current approaches used to study histone PTMs from particular loci involves classical chromatin immunoprecipitation, which due to involvement of formaldehyde crosslinking render the protein part mostly unavailable for MS-based proteomics. We have proposed a novel methodology, which is based upon the biotin tagging of histones proximal to a protein of interest and subsequent purification of nucleosomes carrying the tagged histone. This methodology does not involve any crosslinking, enabling us to purify histones from specific loci, and subject them to large scale MS-based histone PTM analysis. A time dimension can also be added to our approach, as we can follow the modification status of particular fraction of histones once they get biotinylated. Another advantage is the use of alternate variant histones, which allows us to study the PTM profile of different functional states of chromatin. This methodology certainly has an edge on current techniques to study histone PTMs pattern associated with a particular protein of interest or with particular chromatin state.

The first chapter will provide an introduction to the physiology, pathogenesis and biology of mycobacteria. Host-pathogen interactions, different modes of resistance of the bacteria, adaptations for survival under nutrient and oxygen depleted conditions has been discussed. This is followed by a general discussion on gene expression and regulation in the microbe. The physiology of bacteria under stresses from the view of the transcriptional regulation of specific genes has also been discussed. The scope and objective of the present study in M. smegmatis covered in the thesis has been considered at the end. The next chapter discusses the characterization of msdps promoter in vivo with the help of reporter gene assay technologies. With the advent of promoterless E. coli-mycobacterium shuttle vectors, activity assays can be easily performed to characterize unknown upstream putative promoter sequences of genes. Both the 1 kb upstream as well as a 200bp upstream region of msdps gene has been characterized by. Primer extension analysis and subsequent site directed mutagenesis studies reveal +1 transcription start site and the promoter consensus sequence for the msdps gene respectively. Next chapter comprises of the method of constructing heterologous in vitro transcription machinery in mycobacteria. It is followed by characterization of transcription initiation at two dps promoters of M. smegmatis. A novel pull-down assay has been designed which enabled us to identify the sigma factors in the reconstituted RNA polymerases to be associated with the respective dps promoters and to compare the regulation of the two genes at transcription level. Further characterization through single round in vitro transcription at mycobacterial promoters has been attempted. The following two chapters provide some newer insights into the structure-function relationship of the first Dps molecule, MsDps (MsDps1) with respect to its DNA binding activity. The DNA binding activity is associated with the higher oligomeric form only. With the help of time resolved anisotropy and Förster Resonance Energy Transfer (FRET) experiments, we have monitored the nature of Dps dodecamer-DNA complex and mapped the distance between the N and C169 position in the absence and the presence of DNA. A new computational programme, Maximum Entropy Method (MEM) has been applied successfully to analyze data obtained from phase-modulation (Phi-M) lifetime experiments in order to get distribution of lifetime. In the last chapter a new method is adopted to predict amino acids important for stabilizing the interface in a trimeric structure. Subsequently, single and double amino acid mutants of the native MsDps protein has been constructed through site directed mutagenesis and are scored for the ability of the mutants to oligomerize under conditions similar to that of the native protein. This helped us to propose a hypothetical model of the overall mechanism of the protein oligomerization process in solution.

Advanced maternal age has been shown to influence follicular and luteal dynamics in bovine ovary resulting in reduced fertility. The overall objective of the four studies presented in this thesis is to identify the maternal age-associated transcriptional changes in granulosa cells of the dominant follicles during follicle development. In the first study, mRNA expression levels of housekeeping genes were measured by real–time quantitative PCR (RT-qPCR) in granulosa cells of dominant follicles and FSH-stimulated follicles to select and validate suitable reference genes for relative gene expression analyses during maternal and follicular aging. Stability of six reference genes (GAPDH, ACTB, EIF2B2, UBE2D2, SF3A1 and RNF20) was analyzed using GeNorm, DeltaCT and NormFinder programs and comprehensive ranking order was determined based on these programs. Geometric mean of multiple genes (UBE2D2, EIF2B2, GAPDH and SF3A1) was more appropriate reference control than individual genes for the comparison of relative gene expression among dominant and FSH-stimulated follicles during maternal and/or follicular aging studies. In the second study, maternal age-associated changes in the transcriptome of granulosa cells recovered at the time of selection of the dominant follicle from aged (n=3) and young cows (n=3) were determined by EmbryoGENE bovine oligo-microarrays (EMBV3, Agilent Technology). The mRNA expression of five transcripts (CYP19A1, PCNA, GJA1, TPM2, and VNN1) was confirmed in a different set of granulosa cell samples by RT-qPCR to validate microarray data. A total of 169 genes/isoforms were differentially expressed (≥ 2-fold-change; P ≤ 0.05) in aged cows vs. young cows. These transcripts revealed inefficient 1) control of gonadotropins, and gonadotropin-induced changes in the cytoskeleton and extracellular matrix, 2) lipid metabolism and steroidogenesis 3) cell proliferation, cell cycle control and intercellular communication, and 4) higher oxidative stress responses in aged cows vs. young cows. In the third study, changes in the transcriptome of granulosa cells of the preovulatory follicle 24 h after LH treatment from aged (n= 3) and young (n=3) were determined. A total of 1340 genes were expressed differentially (≥ 2-fold change; P ≤ 0.05) in aged cows vs. young cows. The mRNA expression of five transcripts (RGS2, PTGS2, TNFAIP6, VNN1, NR5A2 and GADD45B) was confirmed in a different set of granulosa cell samples to validate microarray data. These transcripts were related to delayed 1) response to LH treatment 2) cellular differentiation and luteinization and 3) progesterone synthesis. Intra-follicle levels of progesterone were lower (P < 0.05) in aged cows compared to young and mid-aged cows. The fourth study compared the aged-associated changes in the transcriptome of granulosa cells during follicle development from the time of dominant follicle selection to preovulatory stage (24 h after LH). In comparison to young cows, aged cows expressed fewer differentially expressed genes/isoforms (1206 vs. 2260, respectively) at ≥ 2-fold-change (P ≤ 0.05) in the granulosa cells of the preovulatory (24 h after LH treatment) vs. the dominant follicle at selection. These transcripts in aged cows were related to late and inefficient 1) organization of cytoskeleton and cytoplasm, 2) differentiation, 3) lipid and cholesterol metabolism, 4) proliferation and 5) higher response to oxidative stress and free radical scavenging in the preovulatory follicles vs. the dominant follicle at selection. In conclusion, maternal age-alters the gene expression of granulosa cells of the dominant follicles during follicle development and results in a compromised follicular environment.

Zelluläre Entscheidungen, wie z. B. die Differenzierung von T-Helferzellen (Th-Zellen) in spezialisierte Effektorlinien, haben großen Einfluss auf die Spezifität von Immunreaktionen. Solche Reaktionen sind das Ergebnis eines komplexen Zusammenspiels einzelner Zellen, die über kleine Signalmoleküle, so genannte Zytokine, kommunizieren. Die hohe Anzahl der Komponenten, sowie deren komplizierte und oft nichtlineare Interaktionen erschweren dabei die Vorhersage, wie bestimmte zelluläre Reaktionen erzeugt werden. Aus diesem Grund sind die globalen Auswirkungen der gezielten Beeinflussung einzelner Zellen oder spezifischer Signalwege nur unzureichend verstanden. So wirken beispielsweise etablierte Behandlungen von Autoimmunkrankheiten oft nur bei einem Teil der Patienten. Durch Einzelzellmethoden wie Live-Cell-Imaging, Massenzytometrie und Einzelzellsequenzierung, können Immunzellen heutzutage quantitativ auf mehreren Ebenen charakterisiert werden. Diese Ansammlung quantitativer Daten erlaubt die Formulierung datengetriebener Modelle zur Vorhersage von zellulären Entscheidungen, allerdings fehlen in vielen Fällen Methoden, um die verschiedenen Daten auf geeignete Weise zu integrieren und zu annotieren. Die vorliegende Arbeit befasst sich mit quantitativen Modellformulierungen für die Entscheidungsfindung von Zellen im Immunsystem mit dem Schwerpunkt auf Lymphozytenproliferation, -differenzierung und -tod. / Cellular decisions, such as the differentiation of T helper (Th) cells into specialized effector lineages, largely impact the direction of immune responses. Such population-level responses are the result of a complex interplay of individual cells which communicate via small signaling molecules called cytokines. The system's complexity, stemming not only from the number of components but also from their intricate and oftentimes non-linear interactions, makes it difficult to develop intuition for how cellular responses are actually generated. Not surprisingly, the global effects of targeting individual cells or specific signaling pathways through perturbations are poorly understood. For instance, common treatments of autoimmune diseases often work for some patients, but not for others. Recently developed methods such as live-cell imaging, mass cytometry and single-cell sequencing now enable quantitative characterization of individual immune cells. This accumulating wealth of quantitative data has laid the basis to derive predictive, data-driven models of immune cell behavior, but in many cases, methods to integrate and annotate the data in a way suitable for model formulation are missing. In this thesis, quantitative workflows and methods are introduced that allow to formulate data-driven models of immune cell decision-making with a particular focus on lymphocyte proliferation, differentiation and death.

Die Funktion und Spezifität von Natürlichen-Killer-T-Zellen (NKT) in angeborener und erworbener Immunität ist nicht vollständig geklärt. Die Mehrheit der NKT-Zellen erkennt alpha-galactosylceramid (alphaGalCer), ein Lipid eines marinen Schwamms mit ungeklärter Relevanz. Verschiedene mykobakterielle Lipide wurden isoliert und auf ihre CD1d-Bindung und NKT-Zell-Aktivierung untersucht. Phospatidylinositol-mannosid (PIM) von Mycobacterium bovis BCG konnte als erstes bakterielles NKT-Zell-Antigen identifiziert werden. PIM aktiviert CD1d-abhängig murine und humane NKT-Zellen zur IFN-gamma aber nicht zur IL-4 Produktion. Mehrere andere Lipid-Fraktionen aktivierten ebenfalls NKT-Zellen. Diese Stimulation war entweder eine direkte, T-Zell-Rezeptor (TZR)-vermittelte und/oder indirekte, Toll-like-receptor 2 (TLR2) vermittelte Aktivierung. Iso-globotrihexosylceramide (iGb3) wurde als das endogene NKT-Zell-Antigen beschrieben. IGb3 ist ubiquitär in Lysosomen vorhanden. Dies wirft die Frage nach der Regulation der Antigen-Verfügkarkeit und der Kontrolle der NKT-Zell-Aktivierung auf. Es konnte gezeigt werden dass die Regulation der Antigen-Verfügbarkeit essentiell für die Regulation der NKT-Zell-Aktivität ist. Unkontrolliertes Auftreten und erhöhte Konzentration von iGb3 führte zu einer substantiellen Reduktion der NKT-Zell-Zahl, vermutlich durch Aktivierungs-induziertem-Zelltod. Mit Hilfe von DNS-Microarray Analysen wurden die Gen-Expressionsprofile von naïven NKT-Zellen und klassischen CD4 T-Zellen, regulatorischen T-Zellen, NK-Zellen und aktivierten NKT-Zellen verglichen. Es konnte sowohl ein NKT-Zell-spezifisches Expressionsmuster etabliert als auch eine gemeinsame Expression von Genen in allen verglichenen Zelltypen identifiziert werden. Naive und aktivierte NKT-Zellen zeigen eine erhöhte Expression von Apoptose-regulierenden Genen welches auf eine starke Selbst-Kontrolle zur präzisen Regulation der eigenen Aktivität hinweist. / The function and specificity of Natural Killer T (NKT) cells in innate and acquired immunity still remains elusive. The vast majority of CD1d restricted NKT cells recognise alpha-galactosylceramid (alphaGalCer), derived from a marine sponge, a lipid of unclear physiological significance. Different mycobycterial glycolipids were isolated and examined for binding to CD1d as well as for their capacity of NKT cell stimulation. Phospatidylinositol-mannoside (PIM) derived from Mycobacterium bovis BCG was identified as the first bacterial lipid antigen presented by CD1d. PIM activated both murine and human NKT cells to secrete IFN-gamma but not IL-4 in a CD1d dependent manner. Additionally, several other lipid fractions with NKT cell activation capacities were identified. This activation was either a direct, T-cell-receptor (TCR) mediated and/or an indirect, toll-like-receptor 2 (TLR2) mediated activation. Iso-globotrihexosylceramide (iGb3) was described as the endogenous NKT cell antigen. iGb3 is a ubiquitously present lysosomal glycolipid which raises the question of regulation of antigen availability and NKT cell activation control. It could be shown that regulation of antigen availability plays a crucial role in regulation of NKT cell activation. Moreover, uncontrolled appearance and increased concentrations of the endogenous antigen iGb3 led to substantial decrease in NKT cell number, presumably by activation induced cell death. Using DNA Microarray analysis, the gene expression profiles of naïve NKT cells and classical CD4 T cells, regulatory T cells, NK cells as well as to activated NKT cells were compared. The profiles revealed a NKT cell specific gene expression pattern as well as expression of genes which NKT cells share with NK cells, conventional CD4+ T cells and Treg cells. Both, naïve and activated NKT cells display elevated expression of apoptosis regulating genes providing NKT cells with high degree of self-control to precisely regulate their own activity.

We have investigated genetic variation in microRNA-mediated regulation as a susceptibility factor for anxiety disorders following two different approaches. We first studied two isoforms of the candidate gene NTRK3 by re-sequencing its different 3'UTRs in patients with Panic (PD) and Obsessive Compulsive disorders (OCD) as well as controls. Two rare variants that altered microRNA-mediated regulation were identified in PD. Conversely, association of a common SNP with OCD hoarding subtype was found. Moreover, we have also studied a possible involvement of microRNAs in anxiety disorders. Consequently, we have analysed the genomic organisation and genetic variation of miRNA-containing regions to construct a panel of SNPs for association analysis. Case-control studies revealed several associations. However, it is worth remarking the associations of miR-22 and miR-488 with PD; two microRNAs for which functional assays and transcriptome analysis after microRNA overexpression showed significant repression of a subset of genes involved in physiological pathways linked to PD development. / Hem investigat la variació genètica a la regulació mediada per microRNAs com a factors de susceptibilitat pels trastorns d'ansietat seguint dues aproximacions diferents. Primer vam estudiar dues isoformes del gen candidat NTRK3 mitjançant la reseqüenciació dels seus diferents 3'UTRs a pacients de pànic (TP), a pacients amb trastorn obsessiu compulsiu (TOC) i a controls. Dues variants rares que alteren la regulació mediada per microRNAs foren identificades per TP. D'altra banda, es trobà associació d'un SNP comú amb el subtipus acumulador de TOC. A més, també hem estudiat la possible implicació dels microRNAs als trastorns d'ansietat. Conseqüentment, hem analitzat l'organització genòmica i la variació genètica a regions que contenen microRNAs per construir un panell d'SNPs per fer anàlisis d'associació. Els estudis cas-control van revelar algunes associacions. Tanmateix, val la pena destacar les associacions del miR-22 i el miR-488 amb TP; dos microRNAs pels quals assajos funcionals i anàlisis de transcriptoma després de la seva sobreexpressió han mostrat una repressió significativa d'un grup de gens implicats en vies fisiològiques lligades al desenvolupament del TP.

Nitrogen (N), an essential plant macronutrient, is among the most limiting factors of crop yield. To sustain modern agriculture, N is often amended in soil in the form of chemical N fertilizer, a major anthropogenic contributor to nutrient pollution that affects climate, biodiversity and human health. To achieve agricultural sustainability, a comprehensive understanding of the regulation of N response in plants is required, in order to engineer crops with higher N use efficiency. Recently, epigenetic mechanisms, such as histone modifications, have gained increasing importance as a new layer of regulation of biological processes. However, our understanding of how epigenetic processes regulate N uptake and assimilation is still in its infancy. To fill this knowledge gap, we first performed a meta-analysis that combined functional genomics and network inference approaches to identify a set of N-responsive epigenetic regulators and predict their effects in regulating epigenome and transcriptome during plant N response. Our analysis suggested that histone modifications could serve as a regulatory mechanism underlying the global transcriptomic reprogramming during plant N response. To test this hypothesis, I applied chromatin immunoprecipitation-sequencing (ChIP-Seq) to monitor the genome-wide changes of four histone marks (H3K27ac, H3K4me3, H3K36me3 and H3K27me3) in response to N supply in tomato plants, followed by RNA-Seq to profile the transcriptomic changes. To investigate the organ specificity of histone modifications, I assayed shoots and roots separately. My results suggest that up to two-thirds of differentially expressed genes (DEGs) are modified in at least one of the four histone marks, supporting an integral role of histone modification in regulating N response. I observed a synergistic modification of active histone marks (H3K27ac, H3K4me3 and H3K36me3) at gene loci functionally relevant to N uptake and assimilation. Surprisingly, I uncovered a non-canonical role of H3K27me3, which is conventionally associated with repressed genes, in modulating active gene expression. Interestingly, such regulatory role of H3K27me3 is specifically associated with highly expressed genes or low expressed genes, depending on the organ context. Overall, I revealed the multi-faceted role of histone marks in mediating the plant N response, which will guide breeding and engineering of better crops with higher N use efficiency