Transcriptome profiling of Paracoccidioides brasiliensis yeast-phase cells recovered from infected mice brings new insights into fungal response upon host interaction

Paracoccidioides brasiliensis is a fungal human pathogen with a wide distribution in Latin America. It causes paracoccidioidomycosis, the most widespread systemic mycosis in Latin America. Although gene expression in P. brasiliensis had been studied, little is known about the genome sequences expressed by this species during the infection process. To better understand the infection process, 4934 expressed sequence tags (ESTs) derived from a non-normalized cDNA library from P. brasiliensis (isolate Pb01) yeast-phase cells recovered from the livers of infected mice were annotated and clustered to a UniGene (clusters containing sequences that represent a unique gene) set with 1602 members. A large-scale comparative analysis was performed between the UniGene sequences of P. brasiliensis yeast-phase cells recovered from infected mice and a database constructed with sequences of the yeast-phase and mycelium transcriptome (isolate Pb01) (https://dna.biomol.unb.br/Pb/), as well as with all public ESTs available at GenBank, including sequences of the P. brasiliensis yeast-phase transcriptome (isolate Pb18) (http:// www.ncbi.nlm.nih.gov/). The focus was on the overexpressed and novel genes. From the total, 3184 ESTs (64.53%) were also present in the previously described transcriptome of yeast-form and mycelium cells obtained from in vitro cultures (https://dna.biomol.unb.br/Pb/) and of those, 1172 ESTs (23.75% of the described sequences) represented transcripts overexpressed during the infection process. Comparative analysis identified 1750 ESTs (35.47% of the total), comprising 649 UniGene sequences representing novel transcripts of P. brasiliensis, not previously described for this isolate or for other isolates in public databases. KEGG pathway mapping showed that the novel and overexpressed transcripts represented standard metabolic pathways, including glycolysis, amino acid biosynthesis, lipid and sterol metabolism. The unique and divergent representation of transcripts in the cDNA library of yeast cells recovered from infected mice suggests differential gene expression in response to the host milieu.

Comparative transcriptome profiling of human foreskin fibroblasts infected with the Sylvio and Y strains of Trypanosoma cruzi

Trypanosoma cruzi, the causative agent of Chagas Disease, is phylogenetically distributed into nearly identical genetic strains which show divergent clinical presentations including differences in rates of cardiomyopathy in humans, different vector species and transmission cycles, and differential congenital transmission in a mouse model. The population structure of these strains divides into two groups, which are geographically and clinically distinct. The aim of this study was to compare the transcriptome of two strains of T. cruzi, Sylvio vs. Y to identify differences in expression that could account for clinical and biochemical differences. We collected and sequenced RNA from T. cruzi-infected and control Human Foreskin Fibroblasts at three timepoints. Differential expression analysis identified gene expression profiles at different timepoints in Sylvio infections, and between Sylvio and Y infections in both parasite and host. The Sylvio strain parasite and the host response to Sylvio infection largely mirrored the host-pathogen interaction seen in our previous Y strain work. IL-8 was more highly expressed in Sylvio-infected HFFs than in Y-infected HFFs.

Transcriptome and membrane fatty acid analyses reveal different strategies for responding to permeating and non-permeating solutes in the bacterium Sphingomonas wittichii.

ABSTRACT: BACKGROUND: Sphingomonas wittichii strain RW1 can completely oxidize dibenzo-p-dioxins and dibenzofurans, which are persistent contaminants of soils and sediments. For successful application in soil bioremediation systems, strain RW1 must cope with fluctuations in water availability, or water potential. Thus far, however, little is known about the adaptive strategies used by Sphingomonas bacteria to respond to changes in water potential. To improve our understanding, strain RW1 was perturbed with either the cell-permeating solute sodium chloride or the non-permeating solute polyethylene glycol with a molecular weight of 8000 (PEG8000). These solutes are assumed to simulate the solute and matric components of the total water potential, respectively. The responses to these perturbations were then assessed and compared using a combination of growth assays, transcriptome profiling, and membrane fatty acid analyses. RESULTS: Under conditions producing a similar decrease in water potential but without effect on growth rate, there was only a limited shared response to perturbation with sodium chloride or PEG8000. This shared response included the increased expression of genes involved with trehalose and exopolysaccharide biosynthesis and the reduced expression of genes involved with flagella biosynthesis. Mostly, the responses to perturbation with sodium chloride or PEG8000 were very different. Only sodium chloride triggered the increased expression of two ECF-type RNA polymerase sigma factors and the differential expression of many genes involved with outer membrane and amino acid metabolism. In contrast, only PEG8000 triggered the increased expression of a heat shock-type RNA polymerase sigma factor along with many genes involved with protein turnover and repair. Membrane fatty acid analyses further corroborated these differences. The degree of saturation of membrane fatty acids increased after perturbation with sodium chloride but had the opposite effect and decreased after perturbation with PEG8000. CONCLUSIONS: A combination of growth assays, transcriptome profiling, and membrane fatty acid analyses revealed that permeating and non-permeating solutes trigger different adaptive responses in strain RW1, suggesting these solutes affect cells in fundamentally different ways. Future work is now needed that connects these responses with the responses observed in more realistic scenarios of soil desiccation.

Genome-wide methylation and transcriptome analysis in penile carcinoma: uncovering new molecular markers

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Transcriptome-Wide Mapping of Pea Seed Ageing Reveals a Pivotal Role for Genes Related to Oxidative Stress and Programmed Cell Death

Understanding of seed ageing, which leads to viability loss during storage, is vital for ex situ plant conservation and agriculture alike. Yet the potential for regulation at the transcriptional level has not been fully investigated. Here, we studied the relationship between seed viability, gene expression and glutathione redox status during artificial ageing of pea (Pisum sativum) seeds. Transcriptome-wide analysis using microarrays was complemented with qRT-PCR analysis of selected genes and a multilevel analysis of the antioxidant glutathione. Partial degradation of DNA and RNA occurred from the onset of artificial ageing at 60% RH and 50 degrees C, and transcriptome profiling showed that the expression of genes associated with programmed cell death, oxidative stress and protein ubiquitination were altered prior to any sign of viability loss. After 25 days of ageing viability started to decline in conjunction with progressively oxidising cellular conditions, as indicated by a shift of the glutathione redox state towards more positive values (>-190 mV). The unravelling of the molecular basis of seed ageing revealed that transcriptome reprogramming is a key component of the ageing process, which influences the progression of programmed cell death and decline in antioxidant capacity that ultimately lead to seed viability loss.

Phenotypic consequences of copy number variation: insights from Smith-Magenis and Potocki-Lupski syndrome mouse models.

A large fraction of genome variation between individuals is comprised of submicroscopic copy number variation of genomic DNA segments. We assessed the relative contribution of structural changes and gene dosage alterations on phenotypic outcomes with mouse models of Smith-Magenis and Potocki-Lupski syndromes. We phenotyped mice with 1n (Deletion/+), 2n (+/+), 3n (Duplication/+), and balanced 2n compound heterozygous (Deletion/Duplication) copies of the same region. Parallel to the observations made in humans, such variation in gene copy number was sufficient to generate phenotypic consequences: in a number of cases diametrically opposing phenotypes were associated with gain versus loss of gene content. Surprisingly, some neurobehavioral traits were not rescued by restoration of the normal gene copy number. Transcriptome profiling showed that a highly significant propensity of transcriptional changes map to the engineered interval in the five assessed tissues. A statistically significant overrepresentation of the genes mapping to the entire length of the engineered chromosome was also found in the top-ranked differentially expressed genes in the mice containing rearranged chromosomes, regardless of the nature of the rearrangement, an observation robust across different cell lineages of the central nervous system. Our data indicate that a structural change at a given position of the human genome may affect not only locus and adjacent gene expression but also "genome regulation." Furthermore, structural change can cause the same perturbation in particular pathways regardless of gene dosage. Thus, the presence of a genomic structural change, as well as gene dosage imbalance, contributes to the ultimate phenotype.

Annotating the human genome

2 Abstract2.1 En françaisLe séquençage du génome humain est un pré-requis fondamental à la compréhension de la biologie de l'être humain. Ce projet achevé, les scientifiques ont dû faire face à une tâche aussi importante, comprendre cette suite de 3 milliards de lettres qui compose notre génome. Le consortium ENCODE (ENCyclopedia Of Dna Elements) fût formé comme une suite logique au projet du génome humain. Son rôle est d'identifier tous les éléments fonctionnels de notre génome incluant les régions transcrites, les sites d'attachement des facteurs de transcription, les sites hypersensibles à la DNAse I ainsi que les marqueurs de modification des histones. Dans le cadre de ma thèse doctorale, j'ai participé à 2 sous-projets d'ENCODE. En premier lieu, j'ai eu la tâche de développer et d'optimiser une technique de validation expérimentale à haut rendement de modèles de gènes qui m'a permis d'estimer la qualité de la plus récente annotation manuelle. Ce nouveau processus de validation est bien plus efficace que la technique RNAseq qui est actuellement en train de devenir la norme. Cette technique basée sur la RT-PCR, m'a notamment permis de découvrir de nouveaux exons dans 10% des régions interrogées. En second lieu j'ai participé à une étude ayant pour but d'identifier les extrémités de tous les gènes des chromosomes humains 21 et 22. Cette étude à permis l'identification à large échelle de transcrits chimères comportant des séquences provenant de deux gènes distincts pouvant être à une grande distance l'un de autre.2.2 In EnglishThe completion of the human genome sequence js the prerequisite to fully understand the biology of human beings. This project achieved, scientists had to face another challenging task, understanding the meaning of the 3 billion letters composing this genome. As a logical continuation of the human genome project, the ENCODE (ENCyclopedia Of DNA Elements) consortium was formed with the aim of annotating all its functional elements. These elements include transcribed regions, transcription binding sites, DNAse I hypersensitive sites and histone modification marks. In the frame of my PhD thesis, I was involved in two sub-projects of ENCODE. Firstly I developed and optimized an high throughput method to validate gene models, which allowed me to assess the quality of the most recent manually-curated annotation. This novel experimental validation pipeline is extremely effective, far more so than transcriptome profiling through RNA sequencing, which is becoming the norm. This RT-PCR-seq targeted-approach is likewise particularly efficient in identifying novel exons, as we discovered about 10% of loci with unannotated exons. Secondly, I participated to a study aiming to identify the gene boundaries of all genes in the human chromosome 21 and 22. This study led to the identification of chimeric transcripts that are composed of sequences coming form two distinct genes that can be map far away from each other.

Homer1a is a core brain molecular correlate of sleep loss.

Sleep is regulated by a homeostatic process that determines its need and by a circadian process that determines its timing. By using sleep deprivation and transcriptome profiling in inbred mouse strains, we show that genetic background affects susceptibility to sleep loss at the transcriptional level in a tissue-dependent manner. In the brain, Homer1a expression best reflects the response to sleep loss. Time-course gene expression analysis suggests that 2,032 brain transcripts are under circadian control. However, only 391 remain rhythmic when mice are sleep-deprived at four time points around the clock, suggesting that most diurnal changes in gene transcription are, in fact, sleep-wake-dependent. By generating a transgenic mouse line, we show that in Homer1-expressing cells specifically, apart from Homer1a, three other activity-induced genes (Ptgs2, Jph3, and Nptx2) are overexpressed after sleep loss. All four genes play a role in recovery from glutamate-induced neuronal hyperactivity. The consistent activation of Homer1a suggests a role for sleep in intracellular calcium homeostasis for protecting and recovering from the neuronal activation imposed by wakefulness.

The loss of circadian PAR bZip transcription factors results in epilepsy.

DBP (albumin D-site-binding protein), HLF (hepatic leukemia factor), and TEF (thyrotroph embryonic factor) are the three members of the PAR bZip (proline and acidic amino acid-rich basic leucine zipper) transcription factor family. All three of these transcriptional regulatory proteins accumulate with robust circadian rhythms in tissues with high amplitudes of clock gene expression, such as the suprachiasmatic nucleus (SCN) and the liver. However, they are expressed at nearly invariable levels in most brain regions, in which clock gene expression only cycles with low amplitude. Here we show that mice deficient for all three PAR bZip proteins are highly susceptible to generalized spontaneous and audiogenic epilepsies that frequently are lethal. Transcriptome profiling revealed pyridoxal kinase (Pdxk) as a target gene of PAR bZip proteins in both liver and brain. Pyridoxal kinase converts vitamin B6 derivatives into pyridoxal phosphate (PLP), the coenzyme of many enzymes involved in amino acid and neurotransmitter metabolism. PAR bZip-deficient mice show decreased brain levels of PLP, serotonin, and dopamine, and such changes have previously been reported to cause epilepsies in other systems. Hence, the expression of some clock-controlled genes, such as Pdxk, may have to remain within narrow limits in the brain. This could explain why the circadian oscillator has evolved to generate only low-amplitude cycles in most brain regions.

Host and invader impact of transfer of the clc genomic island into Pseudomonas aeruginosa PAO1

Genomic islands, large potentially mobile regions of bacterial chromosomes, are a major contributor to bacteria evolution. Here, we investigated the fitness cost and phenotypic differences between the bacterium Pseudomonas aeruginosa PAO1 and a derivative carrying one integrated copy of the clc element, a 103-kb genomic island [and integrative and conjugative element (ICE)] originating in Pseudomonas sp. strain B13 and a close relative of genomic islands found in clinical and environmental isolates of P. aeruginosa. By using a combination of whole genome transcriptome profiling, phenotypic arrays, competition experiments, and biofilm formation studies, only few differences became apparent, such as reduced biofilm growth and fourfold stationary phase repression of genes involved in acetoin metabolism in PAO1 containing the clc element. In contrast, PAO1 carrying the clc element acquired the capacity to grow on 3-chlorobenzoate and 2-aminophenol as sole carbon and energy substrates. No fitness loss >1% was detectable in competition experiments between PAO1 and PAO1 carrying the clc element. The genes from the clc element were not silent in PAO1, and excision was observed, although transfer of clc from PAO1 to other recipient bacteria was reduced by two orders of magnitude. Our results indicate that newly acquired mobile DNA not necessarily invoke an important fitness cost on their host. Absence of immediate detriment to the host may have contributed to the wide distribution of genomic islands like clc in bacterial genomes

A multilevel approach to understanding circadian rhythmicity : social and sleep-dependent effects in ants and mice

Résumé françaisLa majorité des organismes vivants sont soumis à l'alternance du jour et de la nuit, conséquence de la rotation de la terre autour de son axe. Ils ont développé un système interne de mesure du temps, appelé horloge circadienne, leur permettant de s'adapter et de synchroniser leur comportement et leur physiologie aux cycles de lumière. Cette dernière est considérée comme étant le signal majeur entraînant l'horloge interne et. par conséquent, les rythmes journaliers d'éveil et de sommeil. Outre sa régulation circadienne, le sommeil est contrôlé par un processus homéostatique qui détermine son besoin. La contribution de ces deux processus dans le fonctionnement cellulaire du cerveau n'a pas encore été investiguée. La mesure de l'amplitude ainsi que de la prévalence des ondes delta de l'EEG (activité delta) constitue un index très fiable du besoin de sommeil. Il a été démontré que cette activité est génétiquement déterminée et associée à un locus de trait quantitatif situé sur le chromosome 13 de la souris.Grâce à des expériences de privation de sommeil et d'analyses de transcriptome du cerveau dans trois souches de souris présentant diverses réponses à la privation de sommeil, nous avons trouvé que Homerla, localisé dans la région d'intérêt du chromosome 13, est le meilleur marqueur du besoin de sommeil. Homerla est impliqué dans la récupération de l'hyperactivité neuronale induite par le glutamate, grâce à son effet tampon sur le calcium intracellulaire. Une fonction fondamentale du sommeil pourrait donc être de protéger le cerveau et de lui permettre de récupérer après une hyperactivité neuronale imposée par une veille prolongée.De plus, nous avons montré que 2032 transcrits sont exprimés rythmiqueraent dans le cerveau de la souris, parmi lesquels seulement 391 le restent après que les animaux aient été privés de sommeil à différents moments au cours des 24 heures. Cette observation montre clairement que la plupart des changements rythmiques au niveau du transcriptome dépendent du sommeil et non de l'horloge circadienne et souligne ainsi l'importance du sommeil dans la physiologie des mammifères.La plupart des expériences concernant les rythmes circadiens ont été réalisées sur des individus isolés en négligeant l'effet du contexte social sur les comportements circadiens. Les espèces sociales, telles que les fourmis, se caractérisent par une division du travail où une répartition des tâches s'effectue entre ses membres. De plus, certaines d'entre elles doivent être pratiquées en continu comme les soins au couvain tandis que d'autres requièrent une activité rythmique comme le fourragement. Ainsi la fourmi est un excellent modèle pour l'étude de 1 influence du contexte social sur les rythmes circadiens.A ces fins, nous avons décidé d'étudier les rythmes circadiens chez une espèce de fourmi Camponotus fellah et de caractériser au niveau moléculaire son horloge circadienne. Nous avons ainsi développé un système vidéo permettant de suivre l'activité locomotrice de tous les individus d'une colonie. Nos résultats montrent que, bien que la plupart des fourmis soient arythmiques à l'intérieur de la colonie, elles développent d'amples rythmes d'activité en isolation. De plus, ces rythmes disparaissent presque aussitôt que la fourmi est réintroduite dans la colonie. Cette rythmicité observée en isolation semble être générée par l'horloge circadienne car elle persiste en condition constante (obscurité totale). Nous avons ensuite regardé si cette apparente arythmie observée dans la colonie résultait d'un effet masquant des interactions sociales sur les rythmes circadiens d'activité. Nos résultats suggèrent que l'horloge interne est fonctionnelle dans la colonie mais que l'expression de ses rythmes au niveau comportemental est inhibée par les interactions sociales. Les analyses moléculaires du statut de l'horloge dans différents contextes sociaux sont actuellement en cours. Le contexte social semble donc un déterminant majeur du comportement circadien chez la fourmi.AbstractAlmost all living organisms on earth are subjected to the alternance of day and night re-sulting from the rotation of the earth around its axis. They have evolved with an internal timing system, termed the circadian clock, enabling them to adapt and synchronize their behavior and physiology to the daily changes in light and related environmental parame¬ters. Light is thought to be the major cue entraining the circadian clock and consequently the rhythms of rest/activity. In addition to its circadian dependent timing, sleep is reg¬ulated by a homeostatic process that determines its need. The contribution of these two processes in the cellular functioning of the brain has not yet been considered. A highly reliable index of the homeostatic process of sleep is the measure of the amplitude and prevalence of the EEG delta waves (delta activity). It has been shown that sleep need, measured by delta activity, is genetically determined and associated with a Quantitative Trait Locus (QTL) located on the mouse chromosome 13. By using sleep deprivation and brain transcriptome profiling in three inbred mouse strains showing different responses to sleep loss, we found that Homerla, localized within this QTL region is the best transcrip¬tional marker of sleep need. Interestingly Homerla is primarily involved in the recovery from glutamate-induced neuronal hyperactivity by its buffering effect on intracellular cal¬cium. A fundamental function of sleep may therefore reside in the protection and recovery of the brain from a neuronal hyperactivity imposed by prolonged wakefulness.Moreover, time course gene expression experiments showed that 2032 brain tran¬scripts present a rhythmic variation, but only 391 of those remain rhythmic when mice are sleep deprived at four time points around the clock. This finding clearly suggests that most changes in gene transcription over the day are sleep-wake dependent rather than clock dependent and underlines the importance of sleep in mammalian physiology.In the second part of this PhD, I was interested in the social influence on circadian behavior. Most experiments done in the circadian field have been performed on isolated individuals and have therefore ignored the effect of the social context on circadian behav-ior. Eusocial insect species such as ants are characterized by a division of labor: colony tasks are distributed among individuals, some of them requiring continuous activity such as nursing or rhythmic ones such as foraging. Thus ants represent a suitable model to study the influence of the social context on the circadian clock and its output rhythms.The aim of this part was to address the effect of social context on circadian rhythms in the ant species Camponotus fellah and to characterize its circadian clock at the molecu¬lar level. We therefore developed a video tracking system to follow the locomotor activity of all individuals in a colony. Our results show that most ants are arrhythmic within the colony, but develop, when subjected to social isolation, strong rhythms of activity that intriguingly disappear when individuals are reintroduced into the colony. The rhythmicity observed in isolated ants seems to be driven by the circadian clock as it persists under constant conditions (complete darkness). We then tested whether the apparent arrhyth- micity in the colony stemmed from a masking effect of social interactions on circadian rhythms. Indeed, we found that circadian clocks of ants in the colony are functional but their expression at the behavioral level is inhibited by social interactions. The molecular assessment of the circadian clock functional state in the different social context is still under investigation. Our results suggest that social context is a major determinant of circadian behavior in ants.

Differential contribution of transcription factors to Arabidopsis thaliana defense against Spodoptera littoralis

In response to insect herbivory, Arabidopsis plants activate the synthesis of the phytohor- mone jasmonate-isoleucine, which binds to a complex consisting of the receptor COI1 and JAZ repressors. Upon proteasome-mediated JAZ degradation, basic helix-loop-helix tran- scription factors (TFs) MYC2, MYC3, and MYC4 become activated and this results in the expression of defense genes. Although the jasmonate (JA) pathway is known to be essen- tial for the massive transcriptional reprogramming that follows herbivory, there is however little information on other TFs that are required for defense against herbivores and whether they contribute significantly to JA-dependent defense gene expression. By transcriptome profiling, we identified 41TFs that were induced in response to herbivory by the generalist Spodoptera littoralis. Among them, nine genes, including WRKY18, WRKY40, ANAC019, ANAC055, ZAT10, ZAT12, AZF2, ERF13, and RRTF1, were found to play a significant role in resistance to S. littoralis herbivory. Compared to the triple mutant myc234 that is as sensitive as coi1-1 to herbivory, knockout lines of these nine TFs were only partially more sensitive to S. littoralis but, however, some displayed distinct gene expression changes at the whole-genome level. Data thus reveal that MYC2, MYC3, and MYC4 are master regu- lators of Arabidopsis resistance to a generalist herbivore and identify new genes involved in insect defense.

Comparative and functional genomics provide insights into the pathogenicity of dermatophytic fungi.

ABSTRACT: BACKGROUND: Millions of humans and animals suffer from superficial infections caused by a group of highly specialized filamentous fungi, the dermatophytes, which exclusively infect keratinized host structures. To provide broad insights into the molecular basis of the pathogenicity-associated traits, we report the first genome sequences of two closely phylogenetically related dermatophytes, Arthroderma benhamiae and Trichophyton verrucosum, both of which induce highly inflammatory infections in humans. RESULTS: 97% of the 22.5 megabase genome sequences of A. benhamiae and T. verrucosum are unambiguously alignable and collinear. To unravel dermatophyte-specific virulence-associated traits, we compared sets of potentially pathogenicity-associated proteins, such as secreted proteases and enzymes involved in secondary metabolite production, with those of closely related onygenales (Coccidioides species) and the mould Aspergillus fumigatus. The comparisons revealed expansion of several gene families in dermatophytes and disclosed the peculiarities of the dermatophyte secondary metabolite gene sets. Secretion of proteases and other hydrolytic enzymes by A. benhamiae was proven experimentally by a global secretome analysis during keratin degradation. Molecular insights into the interaction of A. benhamiae with human keratinocytes were obtained for the first time by global transcriptome profiling. Given that A. benhamiae is able to undergo mating, a detailed comparison of the genomes further unraveled the genetic basis of sexual reproduction in this species. CONCLUSIONS: Our results enlighten the genetic basis of fundamental and putatively virulence-related traits of dermatophytes, advancing future research on these medically important pathogens.

Combining RT-PCR-seq and RNA-seq to catalog all genic elements encoded in the human genome.

Within the ENCODE Consortium, GENCODE aimed to accurately annotate all protein-coding genes, pseudogenes, and noncoding transcribed loci in the human genome through manual curation and computational methods. Annotated transcript structures were assessed, and less well-supported loci were systematically, experimentally validated. Predicted exon-exon junctions were evaluated by RT-PCR amplification followed by highly multiplexed sequencing readout, a method we called RT-PCR-seq. Seventy-nine percent of all assessed junctions are confirmed by this evaluation procedure, demonstrating the high quality of the GENCODE gene set. RT-PCR-seq was also efficient to screen gene models predicted using the Human Body Map (HBM) RNA-seq data. We validated 73% of these predictions, thus confirming 1168 novel genes, mostly noncoding, which will further complement the GENCODE annotation. Our novel experimental validation pipeline is extremely sensitive, far more than unbiased transcriptome profiling through RNA sequencing, which is becoming the norm. For example, exon-exon junctions unique to GENCODE annotated transcripts are five times more likely to be corroborated with our targeted approach than with extensive large human transcriptome profiling. Data sets such as the HBM and ENCODE RNA-seq data fail sampling of low-expressed transcripts. Our RT-PCR-seq targeted approach also has the advantage of identifying novel exons of known genes, as we discovered unannotated exons in ~11% of assessed introns. We thus estimate that at least 18% of known loci have yet-unannotated exons. Our work demonstrates that the cataloging of all of the genic elements encoded in the human genome will necessitate a coordinated effort between unbiased and targeted approaches, like RNA-seq and RT-PCR-seq.

Transcriptional and Epigenetic Regulation of Human CD4+ T Helper Lineage Specification

Activated T helper (Th) cells have ability to differentiate into functionally distinct Th1, Th2 and Th17 subsets through a series of overlapping networks that include signaling and transcriptional control and the epigenetic mechanisms to direct immune responses. However, inappropriate execution in the differentiation process and abnormal function of these Th cells can lead to the development of several immune mediated diseases. Therefore, the thesis aimed at identifying genes and gene regulatory mechanisms responsible for Th17 differentiation and to study epigenetic changes associated with early stage of Th1/Th2 cell differentiation. Genome wide transcriptional profiling during early stages of human Th17 cell differentiation demonstrated differential regulation of several novel and currently known genes associated with Th17 differentiation. Selected candidate genes were further validated at protein level and their specificity for Th17 as compared to other T helper subsets was analyzed. Moreover, combination of RNA interference-mediated downregulation of gene expression, genome-wide transcriptome profiling and chromatin immunoprecipitation followed by massive parallel sequencing (ChIP-seq), combined with computational data integration lead to the identification of direct and indirect target genes of STAT3, which is a pivotal upstream transcription factor for Th17 cell polarization. Results indicated that STAT3 directly regulates the expression of several genes that are known to play a role in activation, differentiation, proliferation, and survival of Th17 cells. These results provide a basis for constructing a network regulating gene expression during early human Th17 differentiation. Th1 and Th2 lineage specific enhancers were identified from genome-wide maps of histone modifications generated from the cells differentiating towards Th1 and Th2 lineages at 72h. Further analysis of lineage-specific enhancers revealed known and novel transcription factors that potentially control lineage-specific gene expression. Finally, we found an overlap of a subset of enhancers with SNPs associated with autoimmune diseases through GWASs suggesting a potential role for enhancer elements in the disease development. In conclusion, the results obtained have extended our knowledge of Th differentiation and provided new mechanistic insights into dysregulation of Th cell differentiation in human immune mediated diseases.