The transcriptome of equine peripheral blood mononuclear cells.

Complete transcriptomic data at high resolution are available only for a few model organisms with medical importance. The gene structures of non-model organisms are mostly computationally predicted based on comparative genomics with other species. As a result, more than half of the horse gene models are known only by projection. Experimental data supporting these gene models are scarce. Moreover, most of the annotated equine genes are single-transcript genes. Utilizing RNA sequencing (RNA-seq) the experimental validation of predicted transcriptomes has become accessible at reasonable costs. To improve the horse genome annotation we performed RNA-seq on 561 samples of peripheral blood mononuclear cells (PBMCs) derived from 85 Warmblood horses. The mapped sequencing reads were used to build a new transcriptome assembly. The new assembly revealed many alternative isoforms associated to known genes or to those predicted by the Ensembl and/or Gnomon pipelines. We also identified 7,531 transcripts not associated with any horse gene annotated in public databases. Of these, 3,280 transcripts did not have a homologous match to any sequence deposited in the NCBI EST database suggesting horse specificity. The unknown transcripts were categorized as coding and noncoding based on predicted coding potential scores. Among them 230 transcripts had high coding potential score, at least 2 exons, and an open reading frame of at least 300 nt. We experimentally validated 9 new equine coding transcripts using RT-PCR and Sanger sequencing. Our results provide valuable detailed information on many transcripts yet to be annotated in the horse genome.

A transposable element insertion in APOB causes cholesterol deficiency in Holstein cattle

Cholesterol deficiency, a new autosomal recessive inherited genetic defect in Holstein cattle, has been recently reported to have an influence on the rearing success of calves. The affected animals show unresponsive diarrhea accompanied by hypocholesterolemia and usually die within the first weeks or months of life. Here, we show that whole genome sequencing combined with the knowledge about the pedigree and inbreeding status of a livestock population facilitates the identification of the causative mutation. We resequenced the entire genomes of an affected calf and a healthy partially inbred male carrying one copy of the critical 2.24-Mb chromosome 11 segment in its ancestral state and one copy of the same segment with the cholesterol deficiency mutation. We detected a single structural variant, homozygous in the affected case and heterozygous in the non-affected carrier male. The genetic makeup of this key animal provides extremely strong support for the causality of this mutation. The mutation represents a 1.3kb insertion of a transposable LTR element (ERV2-1) in the coding sequence of the APOB gene, which leads to truncated transcripts and aberrant splicing. This finding was further supported by RNA sequencing of the liver transcriptome of an affected calf. The encoded apolipoprotein B is an essential apolipoprotein on chylomicrons and low-density lipoproteins, and therefore, the mutation represents a loss of function mutation similar to autosomal recessive inherited familial hypobetalipoproteinemia-1 (FHBL1) in humans. Our findings provide a direct gene test to improve selection against this deleterious mutation in Holstein cattle.

Gene expression profiling in gills of the great spider crab Hyas araneus in respond to ocean acidification and warming, supplementary data

Hypercapnia and elevated temperatures resulting from climate change may have adverse consequences for many marine organisms. While diverse physiological and ecological effects have been identified, changes in those molecular mechanisms, which shape the physiological phenotype of a species and limit its capacity to compensate, remain poorly understood. Here, we use global gene expression profiling through RNA-Sequencing to study the transcriptional responses to ocean acidification and warming in gills of the boreal spider crab Hyas araneus exposed medium-term (10 weeks) to intermediate (1,120 µatm) and high (1,960 µatm) PCO2 at different temperatures (5°C and 10°C). The analyses reveal shifts in steady state gene expression from control to intermediate and from intermediate to high CO2 exposures. At 5°C acid-base, energy metabolism and stress response related genes were upregulated at intermediate PCO2, whereas high PCO2 induced a relative reduction in expression to levels closer to controls. A similar pattern was found at elevated temperature (10°C). There was a strong coordination between acid-base, metabolic and stress-related processes. Hemolymph parameters at intermediate PCO2 indicate enhanced capacity in acid-base compensation potentially supported by upregulation of a V-ATPase. The likely enhanced energy demand might be met by the upregulation of the electron transport system (ETS), but may lead to increased oxidative stress reflected in upregulated antioxidant defense transcripts. These mechanisms were attenuated by high PCO2, possibly as a result of limited acid-base compensation and metabolic down-regulation. Our findings indicate a PCO2 dependent threshold beyond which compensation by acclimation fails progressively. They also indicate a limited ability of this stenoecious crustacean to compensate for the effects of ocean acidification with and without concomitant warming.

Molecular tools to improve chestnut management: El Bierzo as a case study

The European chestnut (Castanea sativa Mill.) is a multipurpose species that has been widely cultivated around the Mediterranean basin since ancient times. New varieties were brought to the Iberian Peninsula during the Roman Empire, which coexist since then with native populations that survived the last glaciation. The relevance of chestnut cultivation has being steadily growing since the Middle Ages, until the rural decline of the past century put a stop to this trend. Forest fires and diseases were also major factors. Chestnut cultivation is gaining momentum again due to its economic (wood, fruits) and ecologic relevance, and represents currently an important asset in many rural areas of Europe. In this Thesis we apply different molecular tools to help improve current management strategies. For this study we have chosen El Bierzo (Castile and Leon, NW Spain), which has a centenary tradition of chestnut cultivation and management, and also presents several unique features from a genetic perspective (next paragraph). Moreover, its nuts are widely appreciated in Spain and abroad for their organoleptic properties. We have focused our experimental work on two major problems faced by breeders and the industry: the lack of a fine-grained genetic characterization and the need for new strategies to control blight disease. To characterize with sufficient detail the genetic diversity and structure of El Bierzo orchards, we analyzed DNA from 169 trees grafted for nut production covering the entire region. We also analyzed 62 nuts from all traditional varieties. El Bierzo constitutes an outstanding scenario to study chestnut genetics and the influence of human management because: (i) it is located at one extreme of the distribution area; (ii) it is a major glacial refuge for the native species; (iii) it has a long tradition of human management (since Roman times, at least); and (iv) its geographical setting ensures an unusual degree of genetic isolation. Thirteen microsatellite markers provided enough informativeness and discrimination power to genotype at the individual level. Together with an unexpected level of genetic variability, we found evidence of genetic structure, with three major gene pools giving rise to the current population. High levels of genetic differentiation between groups supported this organization. Interestingly, genetic structure does not match with spatial boundaries, suggesting that the exchange of material and cultivation practices have strongly influenced natural gene flow. The microsatellite markers selected for this study were also used to classify a set of 62 samples belonging to all traditional varieties. We identified several cases of synonymies and homonymies, evidencing the need to substitute traditional classification systems with new tools for genetic profiling. Management and conservation strategies should also benefit from these tools. The avenue of high-throughput sequencing technologies, combined with the development of bioinformatics tools, have paved the way to study transcriptomes without the need for a reference genome. We took advantage of RNA sequencing and de novo assembly tools to determine the transcriptional landscape of chestnut in response to blight disease. In addition, we have selected a set of candidate genes with high potential for developing resistant varieties via genetic engineering. Our results evidenced a deep transcriptional reprogramming upon fungal infection. The plant hormones ET and JA appear to orchestrate the defensive response. Interestingly, our results also suggest a role for auxins in modulating such response. Many transcription factors were identified in this work that interact with promoters of genes involved in disease resistance. Among these genes, we have conducted a functional characterization of a two major thaumatin-like proteins (TLP) that belongs to the PR5 family. Two genes encoding chestnut cotyledon TLPs have been previously characterized, termed CsTL1 and CsTL2. We substantiate here their protective role against blight disease for the first time, including in silico, in vitro and in vivo evidence. The synergy between TLPs and other antifungal proteins, particularly endo-p-1,3-glucanases, bolsters their interest for future control strategies based on biotechnological approaches.

Hsf1 and Hsp90 orchestrate temperature-dependent global transcriptional remodelling and chromatin architecture in Candida albicans

We thank Karim Gharbi and Urmi Trivedi for their assistance with RNA sequencing, carried out in the GenePool genomics facility (University of Edinburgh). We also thank Susan Fairley and Eduardo De Paiva Alves (Centre for Genome Enabled Biology and Medicine, University of Aberdeen) for help with the initial bioinformatics analysis. We thank Aaron Mitchell for kindly providing the ALS3 mutant, Julian Naglik for the gift of TR146 cells, and Jon Richardson for technical assistance. We thank the Genomics and Bioinformatics core of the Faculty of Health Sciences for Next Generation Sequencing and Bioinformatics support, the Information and Communication Technology Office at the University of Macau for providing access to a High Performance Computer and Jacky Chan and William Pang for their expert support on the High Performance Computer. Finally, we thank Amanda Veri for generating CaLC2928. M.D.L. is supported by a Sir Henry Wellcome Postdoctoral Fellowship (Wellcome Trust 096072), R.A.F. by a Wellcome Trust-Massachusetts Institute of Technology (MIT) Postdoctoral Fellowship, L.E.C. by a Canada Research Chair in Microbial Genomics and Infectious Disease and by Canadian Institutes of Health Research Grants MOP-119520 and MOP-86452, A.J. P.B. was supported by the UK Biotechnology and Biological Sciences Research Council (BB/F00513X/1) and by the European Research Council (ERC-2009-AdG-249793-STRIFE), KHW is supported by the Science and Technology Development Fund of Macau S.A.R (FDCT) (085/2014/A2) and the Research and Development Administrative Office of the University of Macau (SRG2014-00003-FHS) and R.T.W. by the Burroughs Wellcome fund and NIH R15AO094406. Data availability RNA-sequencing data sets are available at ArrayExpress (www.ebi.ac.uk) under accession code E-MTAB-4075. ChIP-seq data sets are available at the NCBI SRA database (http://www.ncbi.nlm.nih.gov) under accession code SRP071687. The authors declare that all other data supporting the findings of this study are available within the article and its supplementary information files, or from the corresponding author upon request.

Endometrial response toward bovine herpesvirus 4 infection

Metriti ed endometriti sono le patologie maggiormente responsabili delle perdite economiche negli allevamenti bovini da latte, specialmente nel periodo successivo al parto. Mentre le metriti coinvolgono e si sviluppano in tutto l’utero e sono caratterizzate dalla presenza di sintomi sistemici, le endometriti consistono in una infiammazione che riguarda il solo endometrio, con la presenza di perdite purulente, distruzione della superficie epiteliale, congestione vascolare, edema stromale ed accumulo di linfociti e plasmacellule. Queste patologie, inoltre, possono causare, disfunzione ovarica, con conseguente infertilità e riduzione sia dell’efficienza riproduttiva della vacca sia della produzione stessa di latte. Nonostante queste malattie siano, nella maggior parte dei casi, correlate all’instaurarsi di infezioni batteriche, che possono subentrare nell’utero direttamente durante il parto, il ruolo di alcuni virus nello sviluppo di queste patologie è stato recentemente approfondito e la correlazione tra l’ Herpesvirus Bovino 4 e l’insorgere di metriti ed endometriti è stata dimostrata. L’ Herpesvirus Bovino 4 (BoHV-4) è un gamma-herpesvirus ed il suo genoma è costituito da una molecola lineare di DNA a doppio filamento con una struttura genomica di tipo B, caratterizzata dalla presenza di un’unica lunga sequenza centrale (LUR) fiancheggiata da multiple sequenze poli-ripetute (prDNA). BoHV-4 è stato isolato sia da animali sani sia da animali con differenti patologie, tra cui malattie oculari e respiratorie, ma soprattutto da casi di metriti, endometriti, vaginiti o aborti. Generalmente, il ruolo svolto dal virus in questo tipo di patologie è associato alla compresenza di altri tipi di patogeni, che possono essere virus, come nel caso del Virus Della Diarrea Virale Bovina (BVDV), o più frequentemente batteri. Usualmente, l’iniziale difesa dell’endometrio bovino nei confronti dei microbi si fonda sul sistema immunitario innato e l’attivazione di specifici recettori cellulari determina la sintesi e la produzione di citochine e chemochine pro infiammatorie, che possono essere in grado di modulare la replicazione di BoHV-4. Il genoma di BoHV-4 possiede due principali trascritti per i geni Immediate Early (IE), trai quali ORF50/IE2 è il più importante ed il suo prodotto, Rta/ORF50, è fortemente conservato tra tutti gli Herpesvirus. Esso è responsabile della diretta trans-attivazione di numerosi geni virali e cellulari e può essere modulato da differenti stimoli extracellulari. Precedentemente è stato dimostrato come il TNF-, prodotto dalle cellule stromali e dai macrofagi all’interno dell’endometrio, in conseguenza ad infezione batterica, sia in grado di aumentare la replicazione di BoHV-4 attraverso l’attivazione del pathway di NFkB e direttamente agendo sul promotore di IE2. Per queste ragioni, è risultato di forte interesse investigare quali potessero essere, invece, i fattori limitanti la replicazione di BoHV-4. In questo lavoro è stata studiata la relazione tra cellule endometriali stromali bovine infettate con l’Herpesvirus Bovino 4 e l’interferon gamma (IFN-) ed è stata dimostrata la capacità di questa molecola di restringere la replicazione di BoHV-4 in maniera IDO1 indipendente ed IE2 dipendente. Inoltre, la presenza di alcuni elementi in grado di interagire con l’ IFN-γ, all’interno del promotore di IE2 di BoHV-4, ha confermato questa ipotesi. Basandoci su questi dati, abbiamo potuto supporre l’esistenza di uno stretto vincolo tra l’attivazione dell’asse dell’interferon gamma e la ridotta replicazione di BoHV-4, andando a porre le basi per una nuova efficiente cura e prevenzione per le patologie uterine. Poiché il meccanismo corretto attraverso il quale BoHV-4 infetta l’endometrio bovino non è ancora ben compreso, è stato interessante approfondire in maniera più accurata l’interazione presente tra il virus ed il substrato endometriale, analizzando le differenze esistenti tra cellule infettate e non, in termini di espressione genica. Basandoci su dati preliminari ottenuti attraverso analisi con RNA sequencing (RNAseq), abbiamo visto come numerosi geni risultino over-espressi in seguito ad infezione con BoHV-4 e come, tra questi, la Metalloproteasi 1 sia uno dei più interessanti, a causa delle sue possibili implicazioni nello sviluppo delle patologie dell’endometrio uterino bovino. Successive analisi, effettuate tramite westernblotting e real time PCR, sono state in grado di confermare tale dato, sottolineando l’efficacia di un nuovo approccio sperimentale, basato sul RNAseq, per lo studio dell’insorgenza delle patologie.

Role of heme regulated eIF2α kinase in pancreatic beta cell function and viability

The eukaryotic translation initiation factor 2 alpha (eIF2α) is part of the initiation complex that drives the initiator amino acid methionine to the ribosome, a crucial step in protein translation. In stress conditions such as virus infection, endoplasmic reticulum (ER) stress, amino acid or heme deficiency eIF2α can be phosphorylated and thereby inhibit global protein synthesis. This adaptive mechanism prevents protein accumulation and consequent cytotoxic effects. Heme-regulated eIF2α kinase (HRI) is a member of the eIF2α kinase family that regulates protein translation in heme deficiency conditions. Although present in all tissues, HRI is predominantly expressed in erythroid cells where it remains inactive in the presence of normal heme concentrations. In response to heme deficiency, HRI is activated and phosphorylates eIF2α decreasing globin synthesis. This mechanism is important to prevent accumulation of heme-free globin chains which cause ER stress and apoptosis. RNA sequencing data from our group showed that in human islets and in primary rat beta cells HRI is the most expressed eIF2α kinase compared to the other family members. Despite its high expression levels, little is known about HRI function in beta cells. The aim of this project is to identify the role of HRI in pancreatic beta cells. This was investigated taking a loss-of-function approach. HRI knock down (KD) by RNA interference induced beta cell apoptosis in basal condition. HRI KD potentiated the apoptotic effects of palmitate or proinflammatory cytokines, two in vitro models for type 2 and type 1 diabetes, respectively. Increased cytokine-induced apoptosis was also observed in HRI-deficient primary rat beta cells. Unexpectedly, we observed a mild increase in eIF2α phosphorylation in HRI-deficient cells. The levels of mRNA or protein expression of C/EBP homologous protein (CHOP) and activating transcription factor 4 (ATF4) were not modified. HRI KD cells have decreased spliced X-box binding protein 1 (XBP1s), an important branch of the ER stress response. However, overexpression of XBP1s by adenovirus in HRI KD cells did not protect from HRI siRNA-induced apoptosis. HRI deficiency decreased phosphorylation of Akt and its downstream targets glycogen synthase kinase 3 (GSK3), forkhead box protein O1 (FOXO1) and Bcl-2-associated death promoter (BAD). Overexpression of a constitutively active form of Akt by adenovirus in HRI-deficient beta cells partially decreased HRI KD-mediated apoptosis. Interestingly, BAD silencing protected from apoptosis caused by HRI deficiency. HRI silencing in beta cells also induced JNK activation. These results suggest an important role of HRI in beta cell survival through modulation of the Akt/BAD pathway. Thus, HRI may be an interesting target to modulate beta cell fate in diabetic conditions.

PipX, the coactivator of NtcA, is a global regulator in cyanobacteria

To modulate the expression of genes involved in nitrogen assimilation, the cyanobacterial PII-interacting protein X (PipX) interacts with the global transcriptional regulator NtcA and the signal transduction protein PII, a protein found in all three domains of life as an integrator of signals of the nitrogen and carbon balance. PipX can form alternate complexes with NtcA and PII, and these interactions are stimulated and inhibited, respectively, by 2-oxoglutarate, providing a mechanistic link between PII signaling and NtcA-regulated gene expression. Here, we demonstrate that PipX is involved in a much wider interaction network. The effect of pipX alleles on transcript levels was studied by RNA sequencing of S. elongatus strains grown in the presence of either nitrate or ammonium, followed by multivariate analyses of relevant mutant/control comparisons. As a result of this process, 222 genes were classified into six coherent groups of differentially regulated genes, two of which, containing either NtcA-activated or NtcA-repressed genes, provided further insights into the function of NtcA–PipX complexes. The remaining four groups suggest the involvement of PipX in at least three NtcA-independent regulatory pathways. Our results pave the way to uncover new regulatory interactions and mechanisms in the control of gene expression in cyanobacteria.

Analyse du métabolisme du soufre de la bactérie autotrophique acidophile Acidithiobacillus thiooxidans ATCC 19377

Les impacts environnementaux dues à l'extraction minière sont considérables. C'est l'action des microorganismes, en utilisant leur métabolisme du soufre sur les déchets miniers, qui engendre les plus grands défis. Jusqu'à présent, peu de recherches ont été effectués sur les microorganismes environnementaux pour la compréhension globale de l'action du métabolisme du soufre dans une optique de prévention et de rémédiation des impacts environnementaux de l'extraction minière. Dans cette étude, nous avons étudié une bactérie environnementale, Acidithiobacillus thiooxidans, dans le but de comprendre le métabolisme du soufre selon le milieu de culture et le niveau d'acidité du milieu. Nous avons utilisé la transcriptomique à haut débit, RNA-seq, en association avec des techniques de biogéochimie et de microscopie à électrons pour déterminer l'expression des gènes codants les enzymes du métabolisme du soufre. Nous avons trouvé que l'expression des gènes des enzymes du métabolisme du soufre chez ce microorganisme sont dépendantes du milieu, de la phase de croissance et du niveau d'acidité présent dans le milieu. De plus, les analyses biogéochimiques montrent la présence de composés de soufre réduits et d'acide sulfurique dans le milieu. Finalement, une analyse par microscopie électronique révèle que la bactérie emmagasine des réserves de soufre dans son cytoplasme. Ces résultats permettent une meilleure compréhension de son métabolisme et nous rapprochent de la possibilité de développer une technique de prédiction des réactions ayant le potentiel de causer des impacts environnementaux dus à l'extraction minière.

Using the R package Shiny to create web applications that facilitate quantitative gene expression analysis

Thesis (Master's)--University of Washington, 2016-06

Statistical Hurdle Models for Single Cell Gene Expression: Differential Expression and Graphical Modeling

Thesis (Ph.D.)--University of Washington, 2016-06

Understanding Ten-Eleven Translocation-2 in Hematological and Nervous Systems

I proposed the study of two distinct aspects of Ten-Eleven Translocation 2 (TET2) protein for understanding specific functions in different body systems. In Part I, I characterized the molecular mechanisms of Tet2 in the hematological system. As the second member of Ten-Eleven Translocation protein family, TET2 is frequently mutated in leukemic patients. Previous studies have shown that the TET2 mutations frequently occur in 20% myelodysplastic syndrome/myeloproliferative neoplasm (MDS/MPN), 10% T-cell lymphoma leukemia and 2% B-cell lymphoma leukemia. Genetic mouse models also display distinct phenotypes of various types of hematological malignancies. I performed 5-hydroxymethylcytosine (5hmC) chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq) of hematopoietic stem/progenitor cells to determine whether the deletion of Tet2 can affect the abundance of 5hmC at myeloid, T-cell and B-cell specific gene transcription start sites, which ultimately result in various hematological malignancies. Subsequent Exome sequencing (Exome-Seq) showed that disease-specific genes are mutated in different types of tumors, which suggests that TET2 may protect the genome from being mutated. The direct interaction between TET2 and Mutator S Homolog 6 (MSH6) protein suggests TET2 is involved in DNA mismatch repair. Finally, in vivo mismatch repair studies show that the loss of Tet2 causes a mutator phenotype. Taken together, my data indicate that TET2 binds to MSH6 to protect genome integrity. In Part II, I intended to better understand the role of Tet2 in the nervous system. 5-hydroxymethylcytosine regulates epigenetic modification during neurodevelopment and aging. Thus, Tet2 may play a critical role in regulating adult neurogenesis. To examine the physiological significance of Tet2 in the nervous system, I first showed that the deletion of Tet2 reduces the 5hmC levels in neural stem cells. Mice lacking Tet2 show abnormal hippocampal neurogenesis along with 5hmC alternations at different gene promoters and corresponding gene expression downregulation. Through the luciferase reporter assay, two neural factors Neurogenic differentiation 1 (NeuroD1) and Glial fibrillary acidic protein (Gfap) were down-regulated in Tet2 knockout cells. My results suggest that Tet2 regulates neural stem/progenitor cell proliferation and differentiation in adult brain.

Cryptic splicing events in the iron transporter ABCB7 and other key target genes in SF3B1-mutant myelodysplastic syndromes

The splicing factor SF3B1 is the most frequently mutated gene in myelodysplastic syndromes (MDS), and is strongly associated with the presence of ring sideroblasts (RS). We have performed a systematic analysis of cryptic splicing abnormalities from RNA sequencing data on hematopoietic stem cells (HSCs) of SF3B1-mutant MDS cases with RS. Aberrant splicing events in many downstream target genes were identified and cryptic 3' splice site usage was a frequent event in SF3B1-mutant MDS. The iron transporter ABCB7 is a well-recognized candidate gene showing marked downregulation in MDS with RS. Our analysis unveiled aberrant ABCB7 splicing, due to usage of an alternative 3' splice site in MDS patient samples, giving rise to a premature termination codon in the ABCB7 mRNA. Treatment of cultured SF3B1-mutant MDS erythroblasts and a CRISPR/Cas9-generated SF3B1-mutant cell line with the nonsense-mediated decay (NMD) inhibitor cycloheximide showed that the aberrantly spliced ABCB7 transcript is targeted by NMD. We describe cryptic splicing events in the HSCs of SF3B1-mutant MDS, and our data support a model in which NMD-induced downregulation of the iron exporter ABCB7 mRNA transcript resulting from aberrant splicing caused by mutant SF3B1 underlies the increased mitochondrial iron accumulation found in MDS patients with RS.Leukemia advance online publication, 17 June 2016; doi:10.1038/leu.2016.149.

Exploration génétique de l'hypothyroïdie congénitale par dysgénésie thyroïdienne

L'hypothyroïdie congénitale par dysgénésie thyroïdienne (HCDT, ectopie dans plus de 80 %) a une prévalence de 1 cas sur 4000 naissances vivantes. L’HCDT est la conséquence d'une défaillance de la thyroïde embryonnaire à se différencier, à se maintenir ou à migrer vers sa localisation anatomique (partie antérieure du cou), qui aboutit à une absence totale de la thyroïde (athyréose) ou à une ectopie thyroïdienne (linguale ou sublinguale). Les HCDT sont principalement non-syndromiques (soit 98% des cas sont non-familiale), ont un taux de discordance de 92% chez les jumeaux monozygotes, et ont une prédominance féminine et ethnique (i.e., Caucasienne). La majorité des cas d’HCDT n’a pas de cause connue, mais est associée à un déficit sévère en hormones thyroïdiennes (hypothyroïdie). Des mutations germinales dans les facteurs de transcription liés à la thyroïde (NKX2.1, FOXE1, PAX8, NKX2.5) ont été identifiées dans seulement 3% des patients atteints d’HCDT sporadiques et l’analyse de liaisons exclue ces gènes dans les rares familles multiplex avec HCDT. Nous supposons que le manque de transmission familiale claire d’HCDT peut résulter de la nécessité d’au moins deux « hits » génétiques différents dans des gènes importants pour le développement thyroïdien. Pour répondre au mieux nos questions de recherche, nous avons utilisé deux approches différentes: 1) une approche gène candidat, FOXE1, seul gène impliqué dans l’ectopie dans le modèle murin et 2) une approche en utilisant les techniques de séquençage de nouvelle génération (NGS) afin de trouver des variants génétiques pouvant expliquer cette pathologie au sein d’une cohorte de patients avec HCDT. Pour la première approche, une étude cas-contrôles a été réalisée sur le promoteur de FOXE1. Il a récemment été découvert qu’une région du promoteur de FOXE1 est différentiellement méthylée au niveau de deux dinucléotides CpG consécutifs, définissant une zone cruciale de contrôle de l’expression de FOXE1. L’analyse d’association basée sur les haplotypes a révélé qu’un haplotype (Hap1: ACCCCCCdel1C) est associé avec le HCDT chez les Caucasiens (p = 5x10-03). Une réduction significative de l’activité luciférase est observée pour Hap1 (réduction de 68%, p<0.001) comparé au promoteur WT de FOXE1. Une réduction de 50% de l’expression de FOXE1 dans une lignée de cellules thyroïdienne humaine est suffisante pour réduire significativement la migration cellulaire (réduction de 55%, p<0.05). Un autre haplotype (Hap2: ACCCCCCC) est observé moins fréquemment chez les Afro-Américain comparés aux Caucasiens (p = 1.7x10-03) et Hap2 diminue l’activité luciférase (réduction de 26%, p<0.001). Deux haplotypes distincts sont trouvés fréquemment dans les contrôles Africains (Black-African descents). Le premier haplotype (Hap3: GTCCCAAC) est fréquent (30.2%) chez les contrôles Afro-Américains comparés aux contrôles Caucasiens (6.3%; p = 2.59 x 10-9) tandis que le second haplotype (Hap4: GTCCGCAC) est trouvé exclusivement chez les contrôles Afro-Américains (9.4%) et est absent chez les contrôles Caucasiens (P = 2.59 x 10-6). Pour la deuxième approche, le séquençage de l’exome de l’ADN leucocytaire entre les jumeaux MZ discordants n’a révélé aucune différence. D'où l'intérêt du projet de séquençage de l’ADN et l’ARN de thyroïdes ectopiques et orthotopiques dans lesquelles de l'expression monoallélique aléatoire dans a été observée, ce qui pourrait expliquer comment une mutation monoallélique peut avoir des conséquences pathogéniques. Finalement, le séquençage de l’exome d’une cohorte de 36 cas atteints d’HCDT a permis d’identifier de nouveaux variants probablement pathogéniques dans les gènes récurrents RYR3, SSPO, IKBKE et TNXB. Ces quatre gènes sont impliqués dans l’adhésion focale (jouant un rôle dans la migration cellulaire), suggérant un rôle direct dans les défauts de migration de la thyroïde. Les essais de migration montrent une forte diminution (au moins 60% à 5h) de la migration des cellules thyroïdiennes infectées par shRNA comparés au shCtrl dans 2 de ces gènes. Des zebrafish KO (-/- et +/-) pour ces nouveaux gènes seront réalisés afin d’évaluer leur impact sur l’embryologie de la thyroïde.

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

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