Genes da família homeobox têm sido alvo de intensas pesquisas científicas relacionadas ao câncer. Recentemente, mostramos que o gene homeobox TGIF1 está expresso no carcinoma epidermóide de boca na sua forma genérica. Porém, não foi feita uma discriminação entre quais variantes transcricionais, ou subtipos, do TGIF1 estariam expressas. Neste estudo, propusemo-nos a verificar diferenças na freqüência de expressão das variantes transcricionais do TGIF1 em carcinomas epidermóides de boca (CEB) em relação a tecidos morfologicamente não tumorais (TN), avaliar possíveis associações entre essas variantes nos pacientes portadores de CEB e relacionar o grau de expressão dessas variantes com aspectos clínicos, histológicos e com a sobrevida dos pacientes. Adicionalmente, procuramos analisar a expressão da proteína do TGIF1. Foram analisadas 48 amostras congeladas de CEB e 12 de TN. O RNA total de cada amostra foi extraído utilizando-se solução de TRizol®. Os transcritos do TGIF1 foram amplificados por RT-PCR para cada caso de CEB e TN utilizando-se inicialmente um par de iniciadores genéricos para todas as suas variantes. Após essa triagem, os casos positivos foram amplificados utilizando-se pares de iniciadores específicos para cada variante. Não houve diferença estatística entre a freqüência de expressão das variantes do TGIF1 nos grupos CEB e TN, porém as variantes 4 e 8 foram as que apresentaram p-valores menores. Dentro do grupo de CEB, houve associação significativa entre algumas variantes entre si (sete dos vinte e um cruzamentos), de forma que as variantes 1 e 4 foram as que mais tiverem associação com outras variantes. Dessa forma, optamos por prosseguir o estudo utilizando somente as variantes 1, 4 e 8. Não houve correlação entre a expressão das variantes selecionadas e variáveis clinicas e histológicas propostas. Em relação à sobrevida, a expressão das variantes 1 e 8 mostraram correlação univariada com o desfecho óbito, de maneira que o grupo de indivíduos que mais expressou ambas as variantes foi o que apresentou menor risco de óbito. Através da análise multivariada das variantes 1 e 8 e aspectos clínicos e histológicos, somente o tamanho da lesão e a invasão vascular sanguínea tiveram significado estatístico. A expressão protéica do TGIF1 foi analisada em 46 amostras parafinadas considerando-se a diferenciação celular, a graduação imunoistoquímica e o compartimento celular. Os resultados obtidos mostraram que as variantes do TGIF1 estão expressas diferentemente no grupo dos pacientes com CEB e sugerese que a expressão das variantes 1 e 8 estejam relacionadas, de maneira semelhante, a um menor risco de óbito para pacientes portadores de CEB. Adicionalmente, sugere-se que o tamanho da lesão e a invasão vascular sanguinea representem fatores de risco relevantes para o óbito de pacientes portadores de CEB. Sugere-se também que a expressão simultânea da proteína do TGIF1 tanto no núcleo quanto no citoplasma da célula esteja correlacionada a lesões pobremente diferenciadas. / Genes of the homeobox family have been the subject of intense scientific research related to cancer. Recently, we show that this gene is expressed in oral squamous cell carcinoma in its generic form. However, it was not made a discrimination between which transcript variants, or subtypes\", of TGIF1 were expressed. In this study, we aim to verify differences in the expression of the eight transcript variants described for the homeobox gene TGIF1 in oral squamous cell carcinomas (OSCC) related with morphologically non-tumoral tissues (NT), evaluate possible associations between these variants in patients with OSCC and relate the degree of expression of these variants with clinical, histological and survival aspects of patients. Additionally, we analyzed the expression of TGIF1 protein. It was analyzed 48 frozen samples of OSCC and 12 of NT. The total RNA from each sample was extracted using TRizol solution. TGIF1 transcripts were first amplified by RT-PCR for each case of OSCC and NT using a generic pair of primers. After these screening, the positive cases were amplified using specific pair of primers for each variant. There was no statistical difference between the frequency of expression of TGIF1 transcript variants in OSSC and NT groups, but the variants 4 and 8 had the lowest p-values. Within the OSCC group, there was a significant association between some variants among themselves (seven of the twenty-one associations), in a way that the variants 1 and 4 were the ones that had most association with other variants. Thus, we chose to continue the study using only variants 1, 4 and 8. There was no correlation between the expression of the selected variants with the clinical and histological aspects. Regarding survival, the expression of variants 1 and 8 showed statistical correlation with the outcome death, in a way that the group of patients that most expressed both variants was that one with the lower risk of death. By multivariate analysis of variants 1 and 8 and clinical and histological aspects, only the size of the lesion and vascular invasion blood were statistically significant. The protein expression of TGIF1 was analyzed in 46 paraffin-embedded tissues considering the cell differentiation, the immunohistochemistry graduation and the cell compartment. The results showed that the variants of TGIF1 are differently expressed in the OSCC group of patients and it is suggested that the expression of variants 1 and 8 may be related, in a similar way, to a lower risk of death for patients with OSCC. Additionally, it is suggested that the size of the lesion and the vascular invasion of blood represent relevant risk factors for death in patients with OSCC. It is also suggested that the simultaneous expression of TGIF1 protein not only in the nucleus but also in the cytoplasm of the cell is correlated with the poorly differentiated lesions of OSCC.

L'objectif principal de cette thèse est le développement, l'amélioration et l'évaluation de méthodes de traitement de données massives de séquençage, principalement des lectures de séquençage d'ARN courtes et longues, pour éventuellement aider la communauté à répondre à certaines questions biologiques, en particulier dans les contextes de transcriptomique et d'épissage alternatif. Notre objectif initial était de développer des méthodes pour traiter les données d'ARN-seq de deuxième génération à l'aide de graphes de De Bruijn afin de contribuer à la littérature sur l'épissage alternatif, qui a été exploré dans les trois premiers travaux. Le premier article (Chapitre 3, article [77]) a exploré le problème que les répétitions apportent aux assembleurs de transcriptome si elles ne sont pas correctement traitées. Nous avons montré que la sensibilité et la précision de notre assembleur local d'épissage alternatif augmentaient considérablement lorsque les répétitions étaient formellement modélisées. Le second (Chapitre 4, article [11]) montre que l'annotation d'événements d'épissage alternatifs avec une seule approche conduit à rater un grand nombre de candidats, dont beaucoup sont importants. Ainsi, afin d'explorer de manière exhaustive les événements d'épissage alternatifs dans un échantillon, nous préconisons l'utilisation combinée des approches mapping-first et assembly-first. Étant donné que nous avons une énorme quantité de bulles dans les graphes de De Bruijn construits à partir de données réelles d'ARN-seq, qui est impossible à analyser dans la pratique, dans le troisième travail (Chapitre 5, articles [1, 2]), nous avons exploré théoriquement la manière de représenter efficacement et de manière compacte l'espace des bulles via un générateur des bulles. L'exploration et l'analyse des bulles dans le générateur sont réalisables dans la pratique et peuvent être complémentaires aux algorithmes de l'état de l'art qui analysent un sous-ensemble de l'espace des bulles. Les collaborations et les avancées sur la technologie de séquençage nous ont incités à travailler dans d'autres sous-domaines de la bioinformatique, tels que: études d'association à l'échelle des génomes, correction d'erreur et assemblage hybride. Notre quatrième travail (Chapitre 6, article [48]) décrit une méthode efficace pour trouver et interpréter des unitigs fortement associées à un phénotype, en particulier la résistance aux antibiotiques, ce qui rend les études d'association à l'échelle des génomes plus accessibles aux panels bactériens, surtout ceux qui contiennent des bactéries plastiques. Dans notre cinquième travail (Chapitre 7, article [76]), nous évaluons dans quelle mesure les méthodes existantes de correction d'erreur ADN à lecture longue sont capables de corriger les lectures longues d'ARN-seq à taux d'erreur élevé. Nous concluons qu'aucun outil ne surpasse tous les autres pour tous les indicateurs et est le mieux adapté à toutes les situations, et que le choix devrait être guidé par l'analyse en aval. Les lectures longues d'ARN-seq fournissent une nouvelle perspective sur la manière d'analyser les données transcriptomiques, puisqu'elles sont capables de décrire les séquences complètes des ARN messagers, ce qui n'était pas possible avec des lectures courtes dans plusieurs cas, même en utilisant des assembleurs de transcriptome de l'état de l'art. En tant que tel, dans notre dernier travail (Chapitre 8, article [75]), nous explorons une méthode hybride d'assemblage d'épissages alternatifs qui utilise des lectures à la fois courtes et longues afin de répertorier les événements d'épissage alternatifs de manière complète, grâce aux lectures courtes, guidé par le contexte intégral fourni par les lectures longues / The main goal of this thesis is the development, improvement and evaluation of methods to process massively sequenced data, mainly short and long RNA-sequencing reads, to eventually help the community to answer some biological questions, especially in the transcriptomic and alternative splicing contexts. Our initial objective was to develop methods to process second-generation RNA-seq data through de Bruijn graphs to contribute to the literature of alternative splicing, which was explored in the first three works. The first paper (Chapter 3, paper [77]) explored the issue that repeats bring to transcriptome assemblers if not addressed properly. We showed that the sensitivity and the precision of our local alternative splicing assembler increased significantly when repeats were formally modeled. The second (Chapter 4, paper [11]), shows that annotating alternative splicing events with a single approach leads to missing out a large number of candidates, many of which are significant. Thus, to comprehensively explore the alternative splicing events in a sample, we advocate for the combined use of both mapping-first and assembly-first approaches. Given that we have a huge amount of bubbles in de Bruijn graphs built from real RNA-seq data, which are unfeasible to be analysed in practice, in the third work (Chapter 5, papers [1, 2]), we explored theoretically how to efficiently and compactly represent the bubble space through a bubble generator. Exploring and analysing the bubbles in the generator is feasible in practice and can be complementary to state-of-the-art algorithms that analyse a subset of the bubble space. Collaborations and advances on the sequencing technology encouraged us to work in other subareas of bioinformatics, such as: genome-wide association studies, error correction, and hybrid assembly. Our fourth work (Chapter 6, paper [48]) describes an efficient method to find and interpret unitigs highly associated to a phenotype, especially antibiotic resistance, making genome-wide association studies more amenable to bacterial panels, especially plastic ones. In our fifth work (Chapter 7, paper [76]), we evaluate the extent to which existing long-read DNA error correction methods are capable of correcting high-error-rate RNA-seq long reads. We conclude that no tool outperforms all the others across all metrics and is the most suited in all situations, and that the choice should be guided by the downstream analysis. RNA-seq long reads provide a new perspective on how to analyse transcriptomic data, since they are able to describe the full-length sequences of mRNAs, which was not possible with short reads in several cases, even by using state-of-the-art transcriptome assemblers. As such, in our last work (Chapter 8, paper [75]) we explore a hybrid alternative splicing assembly method, which makes use of both short and long reads, in order to list alternative splicing events in a comprehensive manner, thanks to short reads, guided by the full-length context provided by the long reads

Alternativ splicing av pre-mRNA ger upphov till proteindiversitet. Histonmodifieringar kopplas till den alternativa splicingens reglering genom adaptorsystem som overfor den epigenetiska informationen direkt till splicingfaktorerna. De cis- agerande RNA- elementen pa exoner och introner med tillhorande trans- reglerande splicingfaktorer paverkas darfor direkt av specifika histonmodifieringar. En sammankopplande integrerad modell over en rad DNA- baserade processer foreslas. Denna komplexa modell ger en bild av interaktioner och paverkan mellan dessa delar. Kromatin remodellering kravs for bildandet av eukromatin. Nukleosomers placering vid exonrika regioner med specifika modifieringsmonster pekar ut exonerna samt mojliggor inbindning av RNA polymeras II som med sin CTD doman rekryterar bade splicing- och modifieringsfaktorer. Transkriptionshastigheten paverkas av nukleosomplaceringen vilket i sin tur paverkar rekrytering av spliceosomens komponenter, andra trans- agerande regulatorer och aven pre-mRNA sekvensens sekundarstruktur. Kromatin- adaptorkomplex laser av specifika histonmodifieringar och overfor informationen till splicingapparaten. Detta skapar mojlighet till den viktiga cell- och vavnadsspecifika alternativa splicingens reglering. I den integrerade modellen blir komplexiteten tydligare dar alla dessa processer interagerar med varandra och de cis- regulatoriska sekvenserna pa premRNA transkriptet. / Alternative splicing of pre-mRNA generates protein diversity. Histone modifications are connected to the regulation of alternative splicing through adaptor systems that transfers the epigenetic information directly to the splicing factors. The cis- acting RNA elements on the exons and introns together with the trans- regulating splicing factors are therefore directly affected of specific histone modifications. An integrated model over several DNA process mechanisms is suggested. This complex model explains the interactions of the different parts and how they affect each other. Chromatin remodelers are required to obtain euchromatin. Nucleosome positioning at exon rich regions with a specific modification pattern point out where the exons are, and this enable the RNA polymerase II to find and bind to the DNA. It’s CTD domain recruits both splicing- and modifications factors. The transcription rate is also affected of the nucleosome positioning and that in turn affects the recruitment of the components of the spliceosomen, other trans- acting regulators and even the formation of the secondary structure of the pre-mRNA transcript. Chromatin- adaptor complex reads specific histone modifications and transfers this information to the splicing apparatus. All this creates the possibility to regulate important cell- and tissue specific alternative splicing patterns. The integrated model makes the complex processes more clearer when all these integrates with each other and the cis- acting regulating elements on the pre-mRNA transcript.

Background: Genome-wide association studies (GWAs) have identified almost one hundred genetic loci associated with variances in human blood lipid phenotypes including very low-density lipoprotein cholesterol, low-density lipoprotein cholesterol, high-density lipoprotein cholesterol, total cholesterol and triglycerides. Nevertheless the revealed loci only explain a small fraction of heritability and therefore a subtile phenotype of cholesterol homoestasis was examined in our study for the very first time. Methods and Results: Using a multi-stage approach of a GWA, firstly, a genome-wide analysis (Affymetrix 500K GeneChip) for serum lanosterol and serum total cholesterol using LC-MS/MS was conducted in 1495 participants of the KORA-S3/F3 cohort with subsequent replication in two additional independent samples of the the KORA-S3/F3 cohort (n = 1157) and CARLA cohort (n = 1760). Two genetic variants, SNP rs7703051 and rs17562686, in the HMGCR locus were significantly associated with serum lanosterol and showed similar effects of elevated serum lanosterol for each minor allele (combined n = 4412: p = 1,4 x 10-10, +7,1% and p = 4,3 x 10-6, +7,8%). Furthermore, rs7703051 showed a nominal statistical significance to serum cholesterol (p = 0,04). A combined analysis of both SNPs demonstrated that observed associations of rs17562686 can be partly explained by LD with rs7703051 being the primary polymorphism in that study. Nevertheless, rs17562686 shows consistent independent effects on serum lanosterol, thus being associated to a lipid phenotype for the very first time. The following SNP-fine mapping of the HMGCR locus was carried out in the CARLA cohort with subsequent validation in the LE-Heart cohort (n = 1895). The recently published SNP rs3846662 being in tight LD with rs7703051 could be associated with variances of serum lanosterol in both cohorts and functional in vivo studies of gen expression using qRT-PCR assays demonstrated a highly significant association of higher expression of alternatively spliced HMGCR mRNA lacking exon 13 with homozygosity for the rs3846662 major allele in 51 human liver samples (p < 0,01) and 958 human PBMCs (p = 2,1 x 10-7). The overall HMGCR gen expression was not affected. Further investigation of in vivo HMG-CoA reductase enzyme activity in both human samples (n = 48 and n = 55) using anionic exchange column chromatography and scintillation counting of [3-14C]-HMG-CoA and [5-3H]-mevalonolacton did not show any significant results. In addition there was not any association in the LE-Heart cohort between these SNPs and the development of CAD. Finally, rs7703051 could be replicated for already published total cholesterol (combined n = 4412) and rs3846662 for LDL-cholesterol (LE-Heart n = 1895). Since fine mapping in CARLA showed several SNPs throughout the HMGCR locus being in LD with rs17562686 we performed a DNA sequencing of the extended 5´-HMGCR promotor region in six human liver samples. A unknown SNP was discovered in the promotor but could not be associated with any of the examined phenotypes mentioned above. The minor allele of SNP rs5909 situated next to the stop codon and being in high LD with rs17562686 was associated with elevated serum lanosterol and slightly reduced HMGCR gen expression, but further studies including the above mentioned as well as measurement of 3’-UTR transcript lengths using qRT-PCR assays did not produce significant results. Conclusion: The phenotype serum lanosterol could be associated with genetic polymorphisms (e.g. rs7703051) in the HMGCR locus. Therefore already published associations of HMGCR with total cholesterol and LDL-cholesterol can be explained by variances of cholesterol homeostasis. The SNP rs17562686 could be associated with a phenotype of human blood lipids for the very first time. Subsequent gen expression analyses demonstrated a highly significant association of rs3846662 with variant patterns of HMGCR alternative splicing. A significant effect of alternatively spliced protein on enzyme activity and a association of these SNPs with CAD could not be shown.

No description available.

Le cœur des vertébrés est un organe modulaire qui requiert le " patterning " complexe des champs morphogénétiques cardiogènes et la convergence coordonnée des diverses sous-populations de progéniteurs cardiogéniques. Au moins 7 facteurs de transcription de la famille T-box coopèrent au sein de ces nombreuses sous-populations de progéniteurs cardiogéniques afin de réguler la morphogenèse et l’agencement de multiples structures le long de l’ébauche cardiaque, ce qui explique que les mutations humaines de ces gènes engendrent diverses malformations congénitales cardiaques (MCCs). L’un de ces gènes T-box, Tbx5, dont l’haploinsuffisance génère le syndrome de Holt-Oram (SHO), intervient dans une grande variété de réseaux de régulation géniques (RRGs) qui orchestrent la morphogenèse des oreillettes, du ventricule gauche, de la valve mitrale, des septums inter-auriculaire et inter-ventriculaire, ainsi que du système de conduction cardiaque. La diversité des RRGs impliqués dans la formation de ces structures cardiaques suggère que Tbx5 détient une profusion de fonctions qui ne seront identifiables qu’en répertoriant ses activités moléculaires dans chaque lignée cardiaque examinée isolément. Afin d’aborder cette problématique, une ablation génétique de Tbx5 dans l’endocarde a été réalisée. Cette expérience a démontré le rôle crucial de Tbx5 dans la survie des cellules endocardiques bordant le septum primum et des cardiomyocytes au sein de cette structure embryonnaire qui contribuera à la morphogenèse du septum inter-auriculaire. En outre, cette étude a révélé l’existence d’une communication croisée entre la sous-population de cellules endocardiques Tbx5+ et le myocarde au niveau du septum primum, afin d’assurer la survie des cardiomyocytes, et ultimement de garantir la maturation du septum inter-auriculaire. Nos résultats confirment aussi l’importance de l’interdépendance génétique (Tbx5 et Gata4 ainsi que Tbx5 et Nos3) entre différents loci dans la morphogenèse de la cloison inter-auriculaire, et particulièrement de l’influence que peut avoir l’environnement sur la pénétrance et l’expressivité des communications inter-auriculaires (CIAs) dans le SHO. En outre, puisque les fonctions d’un gène dépendent ordinairement des différents isoformes qu’il peut générer, une deuxième étude a focalisé davantage sur l’aspect transcriptionnel de Tbx5. Cette approche a mené à la découverte de 6 transcrits alternatifs exhibant des fonctions à la fois communes et divergentes. La caractérisation de 2 de ces isoformes a révélé le rôle de l’isoforme long (Tbx5_v1) dans la régulation de la croissance des cardiomyocytes durant la cardiogénèse, tandis que l’isoforme court (Tbx5_v2), préférentiellement exprimé dans le cœur mature, réprime la croissance cellulaire. Il est donc entièrement concevable que les mutations de TBX5 entraînant une troncation de la région C-terminale accroissent la concentration d’une protéine mutée qui, à l’instar de Tbx5_v2, interfère avec la croissance de certaines structures cardiaques. En revanche, la divergence de fonctions de ces isoformes, caractérisée par les disparités de localisation subcellulaire et de d’interaction avec d’autres cofacteurs cardiaques, suggère que les mutations affectant davantage un isoforme favoriseraient l’émergence d’un type particulier de MCC. Finalement, un dernier objectif était d’identifier le ou les mécanisme(s) moléculaire(s) par le(s)quel(s) Tbx5 régule son principal gène cible, Nppa, et d’en extraire les indices qui éclairciraient sa fonction transcriptionnelle. Cet objectif nécessitait dans un premier lieu d’identifier les différents modules cis-régulateurs (MCRs) coordonnant la régulation transcriptionnelle de Nppa et Nppb, deux gènes natriurétiques dont l’organisation en tandem et le profil d’expression durant la cardiogénèse sont conservés dans la majorité des vertébrés. L’approche d’empreinte phylogénétique employée pour scanner le locus Nppb/Nppa a permis d’identifier trois MCRs conservés entre diverses espèces de mammifères, dont un (US3) est spécifique aux euthériens. Cette étude a corroboré que la régulation de l’expression du tandem génique Nppb/Nppa requérait l’activité transcriptionnelle d’enhancers en complément aux promoteurs de Nppa et Nppb. La concordance quasiment parfaite entre les profils d’expression de Tbx5 et de ces deux gènes natriurétiques chez les mammifères, suggère que le gradient d’expression ventriculaire de Tbx5 est interprété par le recrutement de ce facteur au niveau des différents enhancers identifiés. En somme, les études présentées dans cette thèse ont permis de clarifier la profusion de fonctions cardiaques que possède Tbx5. Certaines de ces fonctions émanent de l’épissage alternatif de Tbx5, qui favorise la synthèse d’isoformes dotés de propriétés spécifiques. Les diverses interactions combinatoires entre ces isoformes et d’autres facteurs cardiaques au sein des diverses sous-populations de progéniteurs cardiogènes contribuent à l’émergence de RRGs cardiaques divergents. / The vertebrate heart is a modular organ, which requires the complex patterning of the morphogenetic heart fields and the coordinated convergence of the diverse subpopulations of cardiogenic progenitors. At least 7 transcription factors of the T-box family cooperate within these numerous subpopulations of cardiogenic progenitors to regulate the morphogenesis and the layout of multiple structures along the primordial heart tube, which explains that the human mutations of these genes induce various congenital heart defects (CHDs). One of these T-box genes, Tbx5, whose haploinsufficiency generates the Holt-Oram syndrome (HOS), intervenes in a wide variety of gene regulatory networks (GRNs) that orchestrate the morphogenesis of the atria, the left ventricle, the mitral valve, the inter-atrial and inter-ventricular septa, as well as the cardiac conduction system. The diversity of GRNs involved in the formation of these cardiac structures suggests that Tbx5 holds a profusion of functions which will be identifiable only by indexing its molecular activities in each separately examined cardiac lineage. To address this problem, a conditional knockout of Tbx5 in the endocardium was generated. This experiment demonstrated a crucial role of Tbx5 in the survival of the endocardial cells lining the septum primum and the cardiomyocytes within this embryonic structure, which will contribute to the morphogenesis of the inter-atrial septum. Moreover, this study revealed a crosstalk between the Tbx5-positive endocardial cells subpopulation and the myocardium at the level of the septum primum to ensure the survival of cardiomyocytes, and ultimately to guarantee the maturation of the inter-atrial septum. Our results also confirmed the importance of genetic interdependence (Tbx5 and Gata4 as well as Tbx5 and Nos3) between different loci in the morphogenesis of the inter-atrial septum, and particularly the influence that the environment can have on the penetrance and the expressivity of atrial septal defects (ASDs) in the HOS. Besides, since the functions of a gene usually depend on the different isoforms it can generate, a second study focused more on the transcriptional aspect of Tbx5. This approach led to the discovery of 6 alternative transcripts exhibiting both common and specific functions. The characterization of 2 of these isoforms revealed the role of the long isoform (Tbx5_v1) in the regulation of cardiomyocytes growth during cardiogenesis, whereas the short isoform (Tbx5_v2), preferentially expressed in the mature heart, represses cell growth. It is thus entirely conceivable that TBX5 mutations leading to a C-terminal truncation increase the concentration of a mutated protein, which, like Tbx5_v2, interferes with the growth of certain cardiac structures. On the other hand, the divergence of functions of these isoforms, characterized by the disparities of subcellular localization and interaction with other cardiac cofactors, suggests that mutations affecting more one isoform would favor the emergence of a particular type of CHD. Finally, a last objective was to identify one or several molecular mechanism(s) by which Tbx5 regulates its main target gene, Nppa, and to extract clues that might clarify its transcriptional function. This objective required in a first place to identify the various cis-regulatory modules (CRMs) coordinating the transcriptional regulation of Nppa and Nppb, two natriuretic genes whose tandem organization and expression pattern during cardiogenesis are preserved in most vertebrates. The phylogenetic footprint approach employed to scan the Nppb/Nppa locus allowed the identification of three CRMs evolutionary conserved between different mammals species, one of which (US3) is specific to eutherians. This study confirmed that the regulation of the tandem genes Nppb/Nppa required the transcriptional activity of enhancers in complement to Nppa and Nppb promoters. The almost perfect concordance between the expression profiles of Tbx5 and these two natriuretic genes in mammals, suggests that the ventricular expression gradient of Tbx5 is interpreted by the recruitment of this factor to the identified enhancers. Altogether, the studies presented in this thesis allowed clarifying the profusion of Tbx5 cardiac functions. Some of these functions emanate from the alternative splicing of Tbx5, which favors the synthesis of isoforms endowed with specific properties. The diverse combinatorial interactions between these isoforms and other cardiac factors within the various cardiogenic progenitor subpopulations contribute to the emergence of distinct cardiac RRGs.

In this thesis, we address the problem of identifying and quantifying variants (alternative splicing and genomic polymorphism) in RNA-seq data when no reference genome is available, without assembling the full transcripts. Based on the idea that each variant corresponds to a recognizable pattern, a bubble, in a de Bruijn graph constructed from the RNA-seq reads, we propose a general model for all variants in such graphs. We then introduce an exact method, called KisSplice, to extract alternative splicing events and show that it outperforms general purpose transcriptome assemblers. We put an extra effort to make KisSplice as scalable as possible. In order to improve the running time, we propose a new polynomial delay algorithm to enumerate bubbles. We show that it is several orders of magnitude faster than previous approaches. In order to reduce its memory consumption, we propose a new compact way to build and represent a de Bruijn graph. We show that our approach uses 30% to 40% less memory than the state of the art, with an insignificant impact on the construction time. Additionally, we apply the techniques developed to list bubbles in two classical problems: cycle enumeration and the K-shortest paths problem. We give the first optimal algorithm to list cycles in undirected graphs, improving over Johnson's algorithm. This is the first improvement to this problem in almost 40 years. We then consider a different parameterization of the K-shortest (simple) paths problem: instead of bounding the number of st-paths, we bound the weight of the st-paths. We present new algorithms using exponentially less memory than previous approaches

Οι μεταγραφικοί παράγοντες COUP-TF (Chicken Ovalbumin Upstream Promoter Transcription Factor) ανήκουν στην υπεροικογένεια των υποδοχέων στεροειδών/ θυρεοειδών ορμονών και θεωρούνται “ορφανοί”, αφού μέχρι στιγμής δεν έχει βρεθεί το υπεύθυνο πρόσδεμα για την ενεργοποίηση τους. Πειραματικές διαδικασίες στο εργαστήριο μας έχουν δείξει ότι στο πρωτογενές μετάγραφο των COUP-TFs συμβαίνει εναλλακτικό μάτισμα που έχει ως αποτέλεσμα την παραγωγή δύο mRNAs που κωδικοποιούν δύο πρωτεΐνες οι οποίες διαφέρουν ως προς το μέγεθος λόγω της εισαγωγής 21 επίπρόσθετων αμινοξέων στην καρβοξυτελική περιοχή (Carboxy Terminal Extension) της περιοχής πρόσδεσης στο DNA (DBD). Πειράματα EMSA με τη χρήση in vitro μεταφρασμένων πρωτεϊνών, αποκάλυψαν ότι η μεγάλη πρωτεΐνη δεν μπορεί να προσδεθεί σε κανένα COUP-TF στοιχείο απόκρισης. Επίσης παρατηρήθηκε ότι παρουσία της μεγάλης πρωτεΐνης, η ικανότητα της μικρής πρωτεΐνης να προσδένεται στο DNA μειώνεται με ανταγωνιστικό τρόπο, οδηγώντας στο συμπέρασμα ότι το ετεροδιμερές πρωτεϊνικό σύμπλοκο δεν μπορεί να προσδεθεί στο DNA. Σκοπός μας είναι να ερευνήσουμε το ρόλο της ένθεσης των 21 αμινοξέων στην μεγάλη πρωτεΐνη, ως προς την ικανότητα πρόσδεσης της στο DNA. Στον αχινό Paracentrotus lividus, η αμινοξική ένθεση στην καρβοξυτελική περιοχή (CTE) της μεγάλης πρωτεΐνης περιέχει δύο προλίνες. Η υπόθεση μας είναι ότι αυτές οι δύο προλίνες παίζουν σημαντικό ρόλο στην στερεοδιαμόρφωση της πρωτεΐνης, επηρεάζοντας την ικανότητα πρόσδεσης στο DNA. Για να ελέγξουμε την υπόθεση αυτή, προκαλέσαμε σημειακές μεταλλάξεις, μεταλλάσσοντας συγχρόνως τις δύο προλίνες σε αλανίνες αλλά κάθε μία προλίνη σε αλανίνη μεμονωμένα, καθώς επίσης μελετήσαμε και μια σειρά εσωτερικών αμινοξικών ελλειμάτων μέσα στην ένθεση των 21 αμινοξέων στην καρβοξυτελική περιοχή. Το σύνολο των μεταλλάξεων έδειξε ότι οι μεταλλαγμένες μεγάλες πρωτεΐνες δεν προσδένονται στο DNA καθώς επίσης ότι οι μεταλλαγμένες μεγάλες πρωτεΐνες ετεροδιμερίζονται πιο αποτελεσματικά με την μικρή ισομορφή πιθανότατα λόγω επαγόμενης αλλαγής της στερεοδιαμόρφωσης της μεταλλαγμένης πρωτεΐνης. Επίσης μελετήσαμε την εξάπλωση του εναλλακτικού ματίσματός στα Δευτεροστόμια, χρησιμοποιώντας ειδικούς εκφυλισμένους εκκινητές σε πειράματα PCR. Εναλλακτικό μάτισμα των μεταγράφων COUP-TF παρατηρήθηκε επίσης στους οργανισμούς Sphaerechinus granularis (Εχινόδερμο) και Saccoglossus kowalevskii (Ημιχορδωτό). / COUP-TFs (Chicken Ovalbumin Upstream Promoter- Transcription Factors) belong to the superfamily of steroid/thyroid hormone receptors and they are consider orphans since the proper ligand that activates them is not yet found. Experimental procedures in our laboratory have shown that in Echinoderms the alternative splicing of the COUP-TF primary transcript results in two mRNAs which encode two protein variants that differ by a 21 amino acid insertion in the Carboxy Terminal Extension (CTE) of the DNA Binding Domain (DBD). EMSA experiments with the use of in vitro translated proteins revealed that the large protein variant is incapable of binding any COUP-TF response elements. Furthermore, in the presence of the large variant, the small COUP-TF protein ability to bind DNA is diminished in an antagonistic way, suggesting that the heterodimeric protein is also incapable of DNA binding. Our aim is to investigate the role of the 21aa insertion in the large variant regarding the DNA binding affinity. In sea urchins the CTE insertion in the large variant contains two prolines. Our hypothesis is that these two prolines play an important role in the protein‟s conformation which in turn is responsible for the loss of DNA binding. To check this, we created point mutations by mutating both prolines to alanines simultaneously and then each proline to alanine separately. We also analyzed a series of internal amino acid deletions within the 21aa insertion of the CTE. All the mutations proved that the large mutated proteins are incapable of binding DNA and that they heterodimerize more effectively with the small protein variant possibly because of the changed conformation of the large protein variant. We also studied the alternative splicing among Deuterostomes, by using degenerate primers in PCR experiments. We observed that alternative splicing of COUP-TF transcripts occurs in the sea urchin Sphaerechinus granularis and in the hemichordate Saccoglossus kowalevskii.

Les rétrotransposons sont des éléments génétiques mobiles qui constituent presque la moitié de notre génome. Seule la sous-famille L1HS appartenant à la classe des Long Interspersed Element-1(LINE-1 ou L1) a gardé une capacité de mobilité autonome chez l’Homme. Leur mobilisation dans la lignée germinale, mais Aussi dans certains tissus somatiques, contribue à la diversité du génome humain ainsi qu’à certaines maladies comme le cancer. Ainsi, de nouvelles copies de L1 peuvent directement s'intégrer dans des séquences codantes ou régulatrices, et altérer leur fonction. De plus, les séquences L1 contiennent elles-mêmes plusieurs éléments cis-régulateurs et leur insertion à proximité ou dans un gène peut produire des altérations génétiques plus subtiles. Afin d'explorer l'ensemble de ces altérations à l'échelle du génome, nous avons développé un logiciel dédié à l’analyse des données de séquençage d'ARN qui permet d'identifier des transcrits chimériques ou antisens impliquant les L1 et d'annoter ces isoformes en fonction des différents événements d’épissage alternatif subits. Au cours de ce travail, il est apparu que la compréhension du lien entre polymorphisme des insertions et phénotype nécessite une vue complète des différentes copies L1HS présentes chez un individu donné. Afin de disposer d'un catalogue aussi complet que possible de ces polymorphismes identifiés dans des échantillons humains sains ou pathologiques et publiés dans des journaux scientifiques, nous avons développé euL1db, la base de données des insertions de rétrotransposon L1HS chez l’Homme. En conclusion, ce travail aidera à comprendre l’impact des L1 sur l’expression des gènes, à l'échelle du génome. / Retrotransposons are mobile genetics elements, which form almost half of our genome. Only the L1HS subfamily of the Long Interspersed Element-1 class (LINE-1 or L1) has retained the ability to jump autonomously in humans. Their mobilization in the germline – but also in some somatic tissues – contributes to human genetic diversity and to diseases, such as cancer. L1 reactivation can be directly mutagenic by disrupting genes or regulatory sequences. In addition, L1 sequences themselves contain many regulatory cis-elements. Thus, L1 insertions near a gene or within intronic sequences can also produce more subtle genic alterations. To explore L1-mediated genic alterations in a genome-wide manner, we have developed a dedicated RNA-seq analysis software able to identify L1 chimeric or antisense transcripts and to annotate these novel isoforms with their associated alternative splicing events. During the course of this work, it appeared that understanding the link between L1HS insertion polymorphisms and phenotype or disease requires a comprehensive view of the different L1HS copies present in a given individual or sample. To provide a comprehensive summary of L1HS insertion polymorphisms identified in healthy or pathological human samples and published in peer-reviewed journals, we developed euL1db, the European database of L1HS retrotransposon insertions in humans. This work will help understanding the overall impact of L1 insertions on gene expression, at a genome-wide scale.

The extracellular matrix (ECM) is a complex network of molecules secreted by cells, which is essential for providing structural support and facilitating cell processes including adhesion, migration and survival. Tenascin-C is an immunomodulatory ECM protein that exhibits limited expression in healthy tissues, but is transiently elevated at sites of tissue injury, and is persistently expressed in chronic inflammatory diseases and tumours. Alternative splicing of 9 of tenascin-C's fibronectin type III-like domains (FnIII- A1, A2, A3, A4, B, AD2, AD1, C and D) generates enormous diversity in form; yielding 511 possible isoforms. Post-transcriptional modification of tenascin-C has been studied in cancer and during development where disease and tissue specific isoforms exhibit distinct adhesive, migratory and proliferative effects. However, little is known of how tenascin-C is expressed or alternatively spliced during inflammation. This study characterises inflammation and disease specific tenascin-C isoforms made by immune cells and fibroblasts, and investigates their functional relevance. Biosynthesis and alternative splicing of tenascin-C was examined using standard curve qPCR, ELISA, Western blot and confocal immunocytochemistry in resting and activated primary human immune cells, dermal fibroblasts, and in synovial fibroblasts isolated from healthy controls and from osteoarthritis (OA) and rheumatoid arthritis (RA) patients. Based on these data, three recombinant proteins comprising FnIII domains AD2-AD1, B-C-D and B-AD2-AD1-C-D were cloned, expressed and purified, and their impact on cell behaviour including adhesion, morphology and migration was assessed. Basal tenascin-C expression was lower in myeloid and lymphoid cells than fibroblasts, and was induced in all following inflammatory stimulation. Tenascin-C expression was elevated in disease with RA and OA synovial fibroblasts containing higher levels than healthy controls. Alternative splicing following cell activation was cell-type specific: all FnIII except AD2 and AD1 were upregulated in dendritic cells and macrophages, in T-cells all FnIII remained unchanged with FnIII A1 absent; and no change in splicing was observed in activated dermal fibroblasts. Normal and OA synovial fibroblasts exhibited similar tenascin-C splicing patterns, but FnIII B and D were specifically elevated in RA. Functional analysis revealed differences in the adhesion, morphology and migration of myeloid cells and dermal fibroblasts cultured on FnIII AD2-AD1, B-C-D, B-AD2-AD1-C-D and full length tenascin-C substrates; FnIII B-C-D promoted MDDC migration while B-AD2-AD1-C-D promoted fibroblast adhesion, compared to full length tenascin-C. For the first time, this study reveals differences in tenascin-C biosynthesis and alternative splicing by immune cells and fibroblasts following activation with inflammatory stimuli; and starts to reveal how alternative splicing of tenascin-C may influence the behaviours of both stromal and immune cells types during inflammation and in inflammatory diseases.