Molecular genetics of narcolepsy

1 Abstract Sleep is a vital necessity, yet its basic physiological function is still unknown, despite numerous studies both in healthy humans and animal models. The study of patients with sleep disorders may help uncover major biological pathways in sleep regulation and thus shed light on the actual function of sleep. Narcolepsy is a well defined but rare sleep disorder characterized by excessive daytime sleepiness and cataplexy, thought to be caused by a combination of genetic and environmental factors. The aim of this work was to identify genes or genetic variants, which contribute to the pathogenesis of sporadic and familial narcolepsy. Sporadic narcolepsy is the disorder with the strongest human leukocyte antigen (HLA) association ever reported. Since the associated HLA-DRB1 *1501-DQB1 *0602 haplotype is common in the general population (15-25%), it has been suggested that it is necessary but not sufficient for developing narcolepsy. To further define the genetic basis of narcolepsy risk, we performed a genome-wide association study (GWAS) in 562 European individuals with narcolepsy (cases) and 702 ethnically matched controls, with independent replication in 370 cases and 495 controls, all heterozygous for DRB1*1501-DQB1*0602. We found association with a protective variant near HLA-DQA2. Further analysis revealed that the identified SNP is strongly linked to DRB1*03-DQB1*02 and DRBΠ 301-DQB1*0603. Cases almost never carried a trans DRB1*1301-DQB1*0603 haplotype. This unexpected protective HLA haplotype suggests a causal involvement of the HLA region in narcolepsy susceptibility. Familial cases of narcolepsy account for 10% of all narcolepsy cases. However, due to low number of affected family members, narcolepsy families are usually not eligible for genetic linkage studies. We identified and characterized a large Spanish family with 11 affected family members representing the largest ever reported narcolepsy family. We ran a genetic linkage analysis using DNA of 11 affected and 15 unaffected family members and hereby identified a chromosomal candidate region on chromosome 6 encompassing 163 kb with a maximum multipoint LOD score of 5.02. The coding sequences of 4 genes within this haplotype block as well as 2 neighboring genes were screened for pathogenetic mutations in 2 affected and 1 healthy family members. So far no pathogenic mutation could be identified. Further in-depth sequencing of our candidate region as well as whole genome exome sequencing are underway to identify the pathogenic mutation(s) in this family and will further improve our understanding of the genetic basis of narcolepsy. 2 Résumé Le sommeil est un processus vital, dont la fonction physiologique est encore inconnue, malgré de nombreuses études chez des sujets humains sains ainsi que dans des modèles animaux. L'étude de patients souffrant de troubles du sommeil peut permettre la découverte de voies biologiques jouant un rôle majeur dans la régulation du sommeil. L'un de ces troubles, la narcolepsie, est une maladie rare mais néanmoins bien définie, caractérisée par une somnolence diurne excessive accompagnée de cataplexies. Les connaissances actuelles suggèrent qu'une combinaison de facteurs génétiques et environnementaux en est à l'origine. Le but du présent travail était d'identifier !e(s) gène(s) ou les polymorphismes constituant des facteurs de risque dans les formes sporadique et familiale de narcolepsie. La narcolepsie sporadique est la maladie possédant la plus forte association avec le complexe majeur d'histocompatibilité humain (HLA) jamais reportée. La fréquence au sein de la population générale de l'haplotype associé HLA-DRB1*1501- DQB1*0602 (15-25%) suggère que ce dernier est nécessaire, mais pas suffisant, pour (e développement de la maladie. Nous avons voulu approfondir la recherche de facteurs génétiques augmentant le risque de la narcolepsie. A cette fin, nous avons entrepris une étude d'association à l'échelle du génome (genome-wide association study, GWAS) parmi 562 sujets narcoleptiques européens (cas) et 702 individus contrôle de même origine ethnique et nous avons trouvé une association avec un variant protecteur près du gène HLA- DQA2. Ce résultat a été répliqué indépendamment dans 370 cas et 495 contrôles, tous hétérozygotes au locus DRB1*1501-DQB1*0602. Une analyse plus fine montre que le polymorphisme identifié est fortement lié aux allèles DRB1*03-DQB1*02 et DRB1*1301-DQB1*0603. Nous notons que seul un cas était porteur d'un haplotype en trans DRB1*1301-DQBr0603. La découverte de cet allele HLA protecteur suggère que la région HLA joue un rôle causal dans la susceptibilité à la narcolepsie. Dix pourcents des cas de narcolepsie sont familiaux. Cependant, le faible nombre de membres affectés rend ces familles inéligibles pour des études de liaison génétique. Nous avons identifié et caractérisé une grande famille espagnole, dont 11 membres sont atteints par la maladie, ce qui représente la plus grande famille narcoleptique rapportée jusqu'à ce jour. A partir de l'ADN de 11 membres atteints et 15 non- atteints, nous avons identifié par étude de liaison une région candidate de 163 kîlobases (kb) sur le chromosome 6, correspondant à un LOD score multipoints de 5.02. Nous avons cherché, sans succès, des mutations pathogéniques dans la séquence codante de deux gènes situés à l'intérieur de ce segment, ainsi que 4 gènes adjacents. Un séquençage plus approfondi de la région ainsi que le séquençage des exons de tout le génome est en cours et doit s'avérer plus fructueux et révéler la ou tes mutation(s) pathogénique(s) dans cette famille, ce qui contribuerait à une meilleure compréhension des causes génétiques de la narcolepsie. 3 Résumé pour un large public Le sommeil est une nécessité vitale, dont le rôle physiologique exact reste inconnu malgré de nombreuses études sur des sujets humains sains ainsi que sur des modèles animaux. C'est pourquoi les troubles du sommeil intéressent les chercheurs, car l'élucidation des mécanismes responsables peut permettre de mieux comprendre le fonctionnement du sommeil normal. La narcolepsie est une maladie du sommeil caractérisée par une somnolence diurne excessive. Les personnes atteintes peuvent s'endormir involontairement à tout moment de la journée, et souffrent également de pertes du tonus musculaire (cataplexie) lors de fortes émotions, par exemple un fou rire. La narcolepsie est une maladie rare, apparaissant dans 1 personne sur 2000. Les connaissances actuelles suggèrent qu'une combinaison de facteurs génétiques et environnementaux en est à l'origine. Nous avons voulu identifier les facteurs génétiques influençant le déclenchement de la maladie, d'abord dans sa forme sporadique, puis dans une famille comptant de nombreux membres atteints. En comparant les variations génétiques de près de 1000 sujets narcoleptiques européens avec ceux de 1200 individus sains, nous avons trouvé chez 30% de ces derniers un variant protecteur, qui diminue de 50 fois le risque de développer la maladie, ce qui constitue le plus puissant facteur génétique protecteur décrit à ce jour. Nous avons ensuite étudié une grande famille espagnole comptant une trentaine de membres, dont 11 sont atteints de narcolepsie. De nouveau, nous avons comparé les variations génétiques des membres atteints avec ceux des membres sains. Nous avons ainsi pu identifier une région dans le génome où se trouverait le(s) gène(s) impliqué(s) dans la maladie dans cette famille, mais n'avons pas encore trouvé le(s) variant(s) exact(s). Une étude plus approfondie devrait permettre de P(les) identifier et ainsi contribuer à l'élucidation des mécanismes menant au développement de la narcolepsie.

The complex SNP and CNV genetic architecture of the increased risk of congenital heart defects in Down syndrome.

Congenital heart defect (CHD) occurs in 40% of Down syndrome (DS) cases. While carrying three copies of chromosome 21 increases the risk for CHD, trisomy 21 itself is not sufficient to cause CHD. Thus, additional genetic variation and/or environmental factors could contribute to the CHD risk. Here we report genomic variations that in concert with trisomy 21, determine the risk for CHD in DS. This case-control GWAS includes 187 DS with CHD (AVSD = 69, ASD = 53, VSD = 65) as cases, and 151 DS without CHD as controls. Chromosome 21-specific association studies revealed rs2832616 and rs1943950 as CHD risk alleles (adjusted genotypic P-values <0.05). These signals were confirmed in a replication cohort of 92 DS-CHD cases and 80 DS-without CHD (nominal P-value 0.0022). Furthermore, CNV analyses using a customized chromosome 21 aCGH of 135K probes in 55 DS-AVSD and 53 DS-without CHD revealed three CNV regions associated with AVSD risk (FDR ≤ 0.05). Two of these regions that are located within the previously identified CHD region on chromosome 21 were further confirmed in a replication study of 49 DS-AVSD and 45 DS- without CHD (FDR ≤ 0.05). One of these CNVs maps near the RIPK4 gene, and the second includes the ZBTB21 (previously ZNF295) gene, highlighting the potential role of these genes in the pathogenesis of CHD in DS. We propose that the genetic architecture of the CHD risk of DS is complex and includes trisomy 21, and SNP and CNV variations in chromosome 21. In addition, a yet-unidentified genetic variation in the rest of the genome may contribute to this complex genetic architecture.

Short communication. SNP analyses of the 5HT1R and STAT genes in Pacific white shrimp, Litopenaeus vannamei

La industria de la producción de camarón es una de las industrias acuícolas que se encuentra en más crecimiento en la actualidad. Los estudios para encontrar marcadores genéticos son muy efectivos para la mejora de sus propiedades y de gran interés para los productores de camarón. En este trabajo se utilizaron seis individuos de una población de Litopenaeus vannamei, donde se encontraron cuatro polimorfismos de nucleótido único (SNPs) en el gen 5HT1R (5-hidroxitriptamina receptor1) y un SNP en el gen STAT (transductor de señal y activador de la transcripción). Sin embargo, el polimorfismo en el gen STAT resultó ser homocigoto en una población diferente utilizada para análisis de asociación. Los presentes análisis revelaron que el alelo C, en dos polimorfismos SNP (C109T y C395G) del gen 5HT1R, tiende a estar asociado con el aumento del peso corporal. Consideramos que hay necesidad de hacer nuevos estudios utilizando una muestra más amplia y diversa de la población en cuestión.

Copy number variation in the porcine genome inferred from a 60 k SNP BeadChip

Background: Recent studies in pigs have detected copy number variants (CNVs) using the Comparative Genomic Hybridization technique in arrays designed to cover specific porcine chromosomes. The goal of this study was to identify CNV regions (CNVRs) in swine species based on whole genome SNP genotyping chips. Results: We used predictions from three different programs (cnvPartition, PennCNV and GADA) to analyze data from the Porcine SNP60 BeadChip. A total of 49 CNVRs were identified in 55 animals from an Iberian x Landrace cross (IBMAP) according to three criteria: detected in at least two animals, contained three or more consecutive SNPs and recalled by at least two programs. Mendelian inheritance of CNVRs was confirmed in animals belonging to several generations of the IBMAP cross. Subsequently, a segregation analysis of these CNVRs was performed in 372 additional animals from the IBMAP cross and its distribution was studied in 133 unrelated pig samples from different geographical origins. Five out of seven analyzed CNVRs were validated by real time quantitative PCR, some of which coincide with well known examples of CNVs conserved across mammalian species. Conclusions: Our results illustrate the usefulness of Porcine SNP60 BeadChip to detect CNVRs and show that structural variants can not be neglected when studying the genetic variability in this species.

Fast and Rigorous Computation of Gene and Pathway Scores from SNP-Based Summary Statistics.

Integrating single nucleotide polymorphism (SNP) p-values from genome-wide association studies (GWAS) across genes and pathways is a strategy to improve statistical power and gain biological insight. Here, we present Pascal (Pathway scoring algorithm), a powerful tool for computing gene and pathway scores from SNP-phenotype association summary statistics. For gene score computation, we implemented analytic and efficient numerical solutions to calculate test statistics. We examined in particular the sum and the maximum of chi-squared statistics, which measure the strongest and the average association signals per gene, respectively. For pathway scoring, we use a modified Fisher method, which offers not only significant power improvement over more traditional enrichment strategies, but also eliminates the problem of arbitrary threshold selection inherent in any binary membership based pathway enrichment approach. We demonstrate the marked increase in power by analyzing summary statistics from dozens of large meta-studies for various traits. Our extensive testing indicates that our method not only excels in rigorous type I error control, but also results in more biologically meaningful discoveries.

Impact of a cis-associated gene expression SNP on chromosome 20q11.22 on bipolar disorder susceptibility, hippocampal structure and cognitive performance.

BackgroundBipolar disorder is a highly heritable polygenic disorder. Recent enrichment analyses suggest that there may be true risk variants for bipolar disorder in the expression quantitative trait loci (eQTL) in the brain.AimsWe sought to assess the impact of eQTL variants on bipolar disorder risk by combining data from both bipolar disorder genome-wide association studies (GWAS) and brain eQTL.MethodTo detect single nucleotide polymorphisms (SNPs) that influence expression levels of genes associated with bipolar disorder, we jointly analysed data from a bipolar disorder GWAS (7481 cases and 9250 controls) and a genome-wide brain (cortical) eQTL (193 healthy controls) using a Bayesian statistical method, with independent follow-up replications. The identified risk SNP was then further tested for association with hippocampal volume (n = 5775) and cognitive performance (n = 342) among healthy individuals.ResultsIntegrative analysis revealed a significant association between a brain eQTL rs6088662 on chromosome 20q11.22 and bipolar disorder (log Bayes factor = 5.48; bipolar disorder P = 5.85×10(-5)). Follow-up studies across multiple independent samples confirmed the association of the risk SNP (rs6088662) with gene expression and bipolar disorder susceptibility (P = 3.54×10(-8)). Further exploratory analysis revealed that rs6088662 is also associated with hippocampal volume and cognitive performance in healthy individuals.ConclusionsOur findings suggest that 20q11.22 is likely a risk region for bipolar disorder; they also highlight the informative value of integrating functional annotation of genetic variants for gene expression in advancing our understanding of the biological basis underlying complex disorders, such as bipolar disorder.

Ocorrência do SNP c.767G>T no gene DNM1responsável pelo colapso induzido pelo exercício em cães da raça Labrador Retriever no Estado de São Paulo

Resumo:O colapso induzido pelo exercício (EIC) é considerado uma síndrome autossômica recessiva que afeta principalmente cães da raça Labrador Retriever. A doença é caracterizada por fraqueza muscular e colapso após exercício intenso. Usualmente, ocorre recuperação clínica após o episódio, mas alguns animais podem vir a óbito. Os sinais clínicos são decorrentes do polimorfismo de base única (SNP) c.767G>T no gene Dynamin 1 (DNM1). O objetivo deste trabalho foi determinar a ocorrência deste SNP em 321 cães da raça Labrador Retriever do Estado de São Paulo. Primers específicos para a amplificação de todo o exon 6 do gene DNM1 foram usados nas PCRs utilizando DNA a partir de amostras de sangue ou swab bucal, a avaliação final foi realizada com sequenciamento direto dos produtos da PCR. Dentre os 321 animais estudados, 3,4 % (11/321) eram homozigotos para o SNP c.767G>T no gene DNM1 e 24,6% (79/321) eram heterozigotos. Somente um dos 11 animais homozigotos apresentavam sinais clínicos compatíveis com a EIC. Este é o primeiro estudo sobre a ocorrência deste SNP no Brasil e considerando que quase 25% dos animais estudados eram heterozigotos, a genotipagem dos animais para este SNP pode ser importante antes dos acasalamentos para cães desta raça. A EIC deve ser considerada nos diagnósticos diferenciais de enfermidades neuromusculares em cães da raça Labrador Retriever.

Prevalence study of SNP c.421G>T in the ADAMTS2 gene responsible for dermatosparaxis in White Dorper sheep in Brazil

Abstract: Dermatosparaxis is an autosomal recessive disorder of connective tissue; the disorder is clinically characterized by skin fragility and hyperextensibility. Dermatosparaxis in White Dorper sheep is caused by a single nucleotide polymorphism (SNP) (c.421G>T) in the ADAM metalloproteinase with thrombospondin type 1 motif, 2 (ADAMTS2) gene. The aim of this study was to investigate the prevalence of this SNP in a White Dorper herd in São Paulo state, Brazil. In this study, we collected blood DNA samples from 303 White Dorper sheep and performed polymerase chain reaction to amplify the SNP region. The samples were sequenced to determine the presence of the SNP in the ADAMTS2 gene. The SNP prevalence in the studied population was 15.5%; this finding indicates that more effective control measures should be used to prevent the inheritance of SNP c.421G>T in the ADAMTS2 gene in Brazilian White Dorper herds.

Étude du polymorphisme intra- et inter-spécifique du gène β-tubuline chez des espèces de champignons mycorhiziens à arbuscules en vue de développer des marqueurs moléculaires

Les champignons mycorhiziens à arbuscules (CMA), classés dans le phylum Glomeromycota, ne peuvent pas être facilement identifiés par la morphologie de leurs spores et leurs mycélia à l'intérieur ou à l'extérieur des racines de leurs hôtes. Ce problème fondamental d'identification rend l'étude de leur diversité, en particulier dans leur habitat naturel (sol et racine) extrêmement difficile. Les gènes ribosomaux ont été largement utilisés pour développer des amorces spécifiques et en inférer des arbres phylogénétiques. Cependant, ces gènes sont très polymorphes et existent en plusieurs copies dans le génome des CMA, ce qui complique l’interprétation des résultats. Dans notre étude, nous avons étudié le polymorphisme intra- et inter-spécifique du gène β-tubuline, présent en faible nombre de copies dans le génome des CMA, afin d’obtenir de nouvelles séquences nucléotidiques pour développer des marqueurs moléculaires. Les gènes β-tubuline amplifiés à partir de l'ADN génomique de cinq espèces du genre Glomus ont été clonés et séquencés. L’analyse des séquences indique un polymorphisme intraspécifique chez trois espèces de CMA. Deux séquences paralogues très variables ont été nouvellement identifiées chez les G. aggregatum, G. fasciculatum et G. cerebriforme. Aucun polymorphisme n’a été détecté chez les G. clarum et G. etunicatum. Toutes les séquences montrent la présence de deux introns hautement variables. La majorité des substitutions ont été localisées dans les exons et sont synonymes à 90%. La conservation des acides aminés suggère un niveau élevé de sélection négative sur le gène β-tubuline et nous permet de confirmer que les CMA représentent un ancien groupe fongique (400 million d’années). L’analyse phylogénétique, réalisée avec vingt et une séquences nucléotidiques du gène β-tubuline, a révélé que les séquences des Glomaceae forment un groupe monophylétique bien supporté, avec les Acaulosporaceae et Gigasporaceae comme groupe frère. Les séquences paralogues nouvellement identifiées chez les G. aggregatum et G. fasciculatum n'ont pas été monophylétiques au sein de chaque espèce. Les oligonucléotides ont été choisis sur la base des régions variables et conservées du gène β-tubuline. Le test PCR des amorces β-Tub.cerb.F/ β-Tub.cerb.R a révélé des bandes spécifiques de 401 pb pour les séquences paralogues du G. cerebriforme. Deux paires d’amorces ont été développées afin d’identifier les séquences du groupe nommé Tub.1. Les tests PCR nous ont permis d’identifier certaines séquences du groupe Tub.1. Une paire d’amorce β-Tub.2.F/ β-Tub.2.R nous a permis d’identifier certaines séquences paralogues du groupe nommé Tub.2. L’analyse d’autres gènes combinée à celle du gène β-tubuline permettra le développement de marqueurs moléculaires plus spécifiques pour l’identification de CMA.

Les polymorphismes des gènes encodant les protéines apoptotiques Bim et Bax : leur rôle dans la réponse thérapeutique chez les enfants ayant la leucémie lymphoblastique aiguë

INTRODUCTION Des réponses thérapeutiques variables aux glucocorticoïdes (GCs) sont observées parmi les patients atteints de la leucémie lymphoblastique aiguë (LLA). Les protéines Bax et Bim ont déjà montré un rôle important dans l’apoptose des cellules leucémiques. L’expression de Bax était plus basse chez les patients leucémiques résistants au médicament, de même une sensibilité diminuée aux GCs a été associée avec une expression réduite de Bim. La différence dans l’expression pourrait être due à des polymorphismes présents dans ces gènes et donc être associés avec la résistance aux GCs. MÉTHODE Dix-huit polymorphismes en régions régulatrices, 2 polymorphismes exoniques et 7 polymorphismes en région 3’UTR de ces gènes ont été analysés chez les témoins (n=50) et ont permis de déterminer un nombre minimal de polymorphismes suffisants pour définir les haplotypes (tagSNPs). Ces 8 polymorphismes ont ensuite été génotypés chez 286 enfants atteints de la LLA et ont été testés pour l’issue de la maladie par l’analyse de survie. RÉSULTATS Une survie sans évènement et une survie sans rechute diminuées ont été observées pour l’haplotype 3 (p=0,03 et p=0,02). Une survie globale diminuée a été associée avec l’homozygotie pour l’allèle exonique T298C>T (p=0,03), de même que pour les haplotypes 1 et 4 (p=0,04 et p=0,02) du gène Bim. CONCLUSION Les polymorphismes ont été associés avec une survie diminuée chez des enfants atteints de LLA. Il reste à tester d’autres polymorphismes présents dans ces deux gènes ainsi qu’à définir leurs fonctions afin de comprendre leurs rôles dans la réponse aux GCs.

L’impact du polymorphisme du gène E2 sur la quantification de la charge virale du VPH-16 dans les maladies précancéreuses du col utérin

Le VPH-16 de même que certains VPH, dont le VPH-18, causent le cancer du col utérin. Son intégration dans le génome humain pourrait être un marqueur de progression de l’infection. Les charges virales totale et intégrée sont présentement mesurées en quantifiant par PCR en temps réel les gènes E6 (RT-E6) et E2 (RT-E2-1) du VPH-16. Nous avons évalué l’impact du polymorphisme du gène E2 sur la quantification de l’ADN du VPH-16 dans des spécimens cliniques. Dans un premier temps, le gène E2 de 135 isolats de VPH-16 (123 appartenaient au clade Européen et 12 à des clades non- Européens) fut séquencé. Ensuite, un test de PCR en temps réel ciblant les séquences conservées dans E2 (RT-E2-2) fut développé et optimisé. Cent trente-neuf spécimens (lavages cervicaux et vaginaux) provenant de 74 participantes (58 séropositives pour le VIH, 16 séronégatives pour le VIH) ont été étudiés avec les trois tests E2 (RT-E2-2), E6 (RT-E6) et E2 (RT-E2-1). Les ratios de la quantité d’ADN de VPH-16 mesuré avec RT-E2-2 et RT-E2-1 dans les isolats Européens (médiane, 1.02; intervalle, 0.64-1.80) et Africains 1 (médiane, 0.80; intervalle, 0.53-1.09) sont similaires (P=0.08). Par contre, les ratios mesurés avec les isolats Africains 2 (médiane, 3.23; intervalle, 1.92-3.49) ou Asiatique- Américains (médiane, 3.78; intervalle, 1.47-37) sont nettement supérieurs à ceux obtenus avec les isolats Européens (P<0.02 pour chaque comparaison). Les distributions des quantités de E2 contenues dans les 139 échantillons mesurées avec RT-E2-2 (médiane, 6150) et RT-E2-1 (médiane, 8960) étaient statistiquement différentes (P<0.0001). Nous avons observé que les charges virales totale (odds ratio (OR) OR, 2.16 95% intervalle de confiance (IC) 1.11-4.19), et épisomale du VPH-16 (OR, 2.14 95% IC 1.09-4.19), mais pas la présence de formes intégrées (OR, 3.72 95% IC 1.03-13.4), sont associées aux néoplasies intraepitheliales cervicales de haut grade (CIN-2,3), et ce, en contrôlant pour des facteurs confondants tels que l’âge, le taux de CD4 sanguin, l’infection au VIH, et le polymorphisme de VPH-16. La proportion des échantillons ayant un ratio E6/E2 > 2 pour les femmes sans lésion intraépithéliale (7 de 35) est similaire à celle des femmes avec CIN-2,3 (5 de 11, p=0.24) ou avec CIN- 1 (4 de14, P=0.65). Le polymorphisme du gène E2 est un facteur qui influence la quantification des charges intégrées de VPH-16.

Étude du polymorphisme du gène majeur d’histocompatibilité de classe IIb (MHIIb) chez l’omble de fontaine (Salvelinus fontinalis)

Les molécules classiques du complexe majeur d’histocompatibilité de classe II (CMHII) sont des glycoprotéines de surface spécialisées dans la présentation de peptides, principalement dérivés de pathogènes extracellulaires, aux récepteurs des lymphocytes T CD4+ afin d’initier la réponse immunitaire adaptative. Elles sont encodées, avec celles du CMH de classe I, par les gènes les plus polymorphiques identifiés jusqu’à maintenant, avec plusieurs loci et une grande diversité allélique à chacun d’eux. De plus, le polymorphisme des gènes du CMHII n’est pas limité qu’aux séquences codantes. Il est également observé dans les promoteurs où on a démontré ses effets sur le niveau d’expression des gènes. La variation de la régulation d’un gène est considérée comme un facteur important et pour laquelle des modifications morphologiques, physiologiques et comportementales sont observées chez tous les organismes. Des séquences d’ADN répétées impliquées dans cette régulation ont été identifiées dans les régions non-codantes des génomes. D’un autre côté, la sélection par les pathogènes permettrait l’évolution et le maintien du polymorphisme des gènes du CMH chez les vertébrés. À ce sujet, plusieurs études ont montré l’implication de différents allèles du CMH dans la résistance ou la susceptibilité aux maladies. Cette étude avait pour objectifs de caractériser le polymorphisme du gène MHIIb chez l’omble de fontaine (Salvelinus fontinalis) et de documenter ses effets au niveau de la survie conférée par des allèles et/ou génotypes particuliers lors d’une infection, ainsi que sur la variation du niveau d’expression du gène dans différentes conditions. Dans une première partie, nous avons identifié un total de 6 allèles du gène MHIIb, désignés Safo-DAB*0101 à Safo-DAB*0601, qui montrent une grande similarité avec les séquences codantes provenant de poissons téléostéens et de l’humain. L’analyse des séquences du domaine b1 a permis de détecter l’effet d’une pression sélective positive pour maintenir le polymorphisme dans cette région de la molécule. Quatre de ces allèles ont été testés lors d’une expérience d’infection avec le pathogène Aeromonas salmonicida afin d’évaluer l’effet qu’ils pouvaient avoir sur la survie des poissons. Nous avons trouvé que l’allèle DAB*0101 était significativement associé à la résistance à la furonculose. En plus d’avoir été identifié chez les individus homozygotes pour cet allèle, l’effet a également été remarqué au niveau de la survie les poissons de génotype DAB*0101/*0201. À l’opposé, les facteurs de risque élevé obtenus pour les génotypes DAB*0201/*0301 et DAB*0301/*0401 suggèrent plutôt une association à la susceptibilité. Étant donné la faible fréquence à laquelle l’allèle DAB*0101 a été retrouvé dans la population, le modèle de la sélection dépendante de la fréquence pourrait expliquer l’avantage conféré par ce dernier et souligne l’importance de ce mécanisme pour le maintien du polymorphisme du gène MHIIb chez l’omble de fontaine. Dans une seconde partie, nous avons rapporté la présence d’un minisatellite polymorphique formé d’un motif de 32 nucléotides dans le second intron du gène MHIIb, et pour lequel un nombre exclusif de répétitions du motif a été associé à chaque allèle (69, 27, 20, 40, 19 et 25 répétitions pour les allèles DAB*0101 à DAB*0601 respectivement). L’expression relative de quatre allèles a été évaluée dans des poissons hétérozygotes aux températures de 6 ºC et 18 ºC. Les résultats indiquent que les allèles possédant un long minisatellite montrent une réduction de l’expression du gène d’un facteur 1,67 à 2,56 par rapport aux allèles qui en contiennent un court. De même, des allèles qui incluent des minisatellites de tailles similaires n’affichent pas de différence significative au niveau de l’abondance du transcrit aux deux températures. De plus, l’effet répressif associé aux longs minisatellites est amplifié à la température de 18 ºC dans des poissons de trois génotypes différents. Nous avons finalement observé une augmentation significative par un facteur 2,08 de l’expression totale du gène MHIIb à la température de 6 ºC. Ces résultats appuient l’implication des séquences d’ADN répétées dans la régulation de l’activité transcriptionnelle d’un gène et suggèrent qu’un minisatellite sensible aux différences de températures pourrait être soumis aux forces sélectives et jouer un rôle important dans l’expression de gènes et l’évolution des organismes poïkilothermes.

Polymorphisme de la capside des papillomavirus appartenant à l’espèce 9 des alphapapillomavirus

Le gène L1 encode pour la protéine majeure de la capside des papillomavirus humains (VPH). L’information relative au polymorphisme de L1 pour les types autres que VPH- 16 est jusqu’ici limitée. Cet ouvrage explore le polymorphisme de L1 en comparant les séquences des types phylogénétiquement apparentés VPH-31, -33, -35 à VPH-16. Des spécimens génitaux recueillis de 732 femmes VIH-séropositives et 323 VIHséronégatives ont été criblés à le recherche d’ADN de VPH par PCR consensus au niveau du gène L1. Les échantillons positifs pour VPH-16 (n=74), -31 (n=74), -33 (n=37) et -35 (n=58) étaient analysés par PCR-séquençage pour la totalité du gène L1. Le nombre de nucléotides substitués pour L1 variait de 19 pour VPH-33 à 52 pour VPH-31. Le rapport du nombre de variantes sur le nombre d’isolats testés était plus élevé pour VPH-31 (56.4%, p=0.05) et VPH-35 (60.3%, p=0.04) comparativement à VPH-16 (40.5%), alors que ce ratio était inférieur pour VPH-33 mais sans différence statistiquement significative (24.3%, p=0.14). La distance entre les variantes était plus grande à l’intérieur des cinq boucles présumément exposées à la surface de la protéine L1 que dans la séquence à l’extérieur (p<0.01) Des variations synonymes étaient observées chez 1.7% (95% CI 1.1- 2.3) des nucléotides intra-boucles et 2.4% (95% CI 1.2-3.7) de ceux extra-boucles. Les variations non-synonymes étaient rencontrées pour 1.8% (95% CI 1.1-2.5) des nucléotides intra-boucles et 0.2% (95% CI 0-0.4) pour les nucléotides extra-boucles. Les ratios dN/dS étaient inférieurs à 1.0 pour les régions extra-boucles et encore davantage pour les régions intra-boucles. Ces résultats suggèrent que les séquences des régions hypervariables de L1 ont été sélectionnées positivement.

Alternative strategies for deciphering the genetic architecture of childhood Pre-B acute lymphoblastic leukemia

La leucémie lymphoblastique aigüe (LLA) est une maladie génétique complexe. Malgré que cette maladie hématologique soit le cancer pédiatrique le plus fréquent, ses causes demeurent inconnues. Des études antérieures ont démontrées que le risque à la LLA chez l’enfant pourrait être influencé par des gènes agissant dans le métabolisme des xénobiotiques, dans le maintient de l’intégrité génomique et dans la réponse au stress oxydatif, ainsi que par des facteurs environnementaux. Au cours de mes études doctorales, j’ai tenté de disséquer davantage les bases génétiques de la LLA de l’enfant en postulant que la susceptibilité à cette maladie serait modulée, au moins en partie, par des variants génétiques agissant dans deux voies biologiques fondamentales : le point de contrôle G1/S du cycle cellulaire et la réparation des cassures double-brin de l’ADN. En utilisant une approche unique reposant sur l’analyse d’une cohorte cas-contrôles jumelée à une cohorte de trios enfants-parents, j’ai effectué une étude d’association de type gènes/voies biologiques candidats. Ainsi, j’ai évaluer le rôle de variants provenant de la séquence promotrice de 12 gènes du cycle cellulaire et de 7 gènes de la voie de réparation de l’ADN, dans la susceptibilité à la LLA. De tels polymorphismes dans la région promotrice (pSNPs) pourraient perturber la liaison de facteurs de transcription et mener à des différences dans les niveaux d’expression des gènes pouvant influencer le risque à la maladie. En combinant différentes méthodes analytiques, j’ai évalué le rôle de différents mécanismes génétiques dans le développement de la LLA chez l’enfant. J’ai tout d’abord étudié les associations avec gènes/variants indépendants, et des essaies fonctionnels ont été effectués afin d’évaluer l’impact des pSNPs sur la liaison de facteurs de transcription et l’activité promotrice allèle-spécifique. Ces analyses ont mené à quatre publications. Il est peu probable que ces gènes de susceptibilité agissent seuls; j’ai donc utilisé une approche intégrative afin d’explorer la possibilité que plusieurs variants d’une même voie biologique ou de voies connexes puissent moduler le risque de la maladie; ces travaux ont été soumis pour publication. En outre, le développement précoce de la LLA, voir même in utero, suggère que les parents, et plus particulièrement la mère, pourraient jouer un rôle important dans le développement de cette maladie chez l’enfant. Dans une étude par simulations, j’ai évalué la performance des méthodes d’analyse existantes de détecter des effets fœto-maternels sous un design hybride trios/cas-contrôles. J’ai également investigué l’impact des effets génétiques agissant via la mère sur la susceptibilité à la LLA. Cette étude, récemment publiée, fût la première à démontrer que le risque de la leucémie chez l’enfant peut être modulé par le génotype de sa mère. En conclusions, mes études doctorales ont permis d’identifier des nouveaux gènes de susceptibilité pour la LLA pédiatrique et de mettre en évidence le rôle du cycle cellulaire et de la voie de la réparation de l’ADN dans la leucémogenèse. À terme, ces travaux permettront de mieux comprendre les bases génétiques de la LLA, et conduiront au développement d’outils cliniques qui amélioreront la détection, le diagnostique et le traitement de la leucémie chez l’enfant.