Analyse neurogénétique du trouble désintégratif de l’enfant

Sujets de l’étude

Les caractéristiques cliniques de la cohorte CDD et le nombre de sujets examinés par chaque modalité d’étude sont indiqués dans le tableau 1. Le sex-ratio de 3,25 hommes pour 1 femme est similaire à celui rapporté pour les TSA. L’âge moyen et médian au moment de l’apparition des symptômes était de 46 et 40 mois, respectivement, avec une fourchette de 28 à 84 mois. Soixante-dix pour cent des sujets ont connu un prodrome d’anxiété et de terreur. Trente pour cent des sujets ont connu plusieurs épisodes de régression. La durée du premier épisode de régression allait de 2 mois à presque 7 ans chez un sujet. La plupart des sujets ont une DI sévère à profonde, le QI moyen et médian étant respectivement de 30 et 26, avec une fourchette de 8 à 74. Tous présentaient une perte du langage, une perte des aptitudes sociales ou du comportement adaptatif, et une perte des aptitudes au jeu. Soixante-cinq pour cent d’entre eux avaient perdu le contrôle de leurs intestins ou de leur vessie, et la même proportion avait perdu leurs capacités motrices. Bien que l’on ait rapporté que le CDD était presque toujours sporadique, quelques-uns de nos sujets avaient des membres de leur famille immédiate présentant des TSA ou des caractéristiques autistiques, notamment deux paires de jumeaux monozygotes. Les deux membres d’une paire (CDD13-03/04) sont atteints de CDD ; dans l’autre paire (CDD20-03/04), l’un est atteint de CDD et l’autre de TSA.

Génétique

Compte tenu de la rareté, de la sévérité et de la transmission apparemment sporadique observée dans la plupart des cas de CDD, nous avons émis l’hypothèse que des variantes rares à effet important contribuent à l’étiologie. En effet, il existe de nombreuses preuves de la contribution des variants rares aux TSA. Comme le montre le tableau 2, nous avons trouvé un ou plusieurs variants rares pour tous les proband sauf un, qui n’étaient pas partagés par les frères et sœurs témoins non affectés (fichier supplémentaire 2 : tableau S2). Nous avons également recherché des variantes hétérozygotes composées chez les sujets en recherchant des variantes supplémentaires dans les gènes affectés par les variantes de novo, mais nous n’en avons trouvé aucune. Les taux de tous les variants de novo à forte probabilité (score de qualité bayésien ≥ 50) étaient de 0,80/exome de bande et de 0,92/exome de fratrie (fichier supplémentaire 2 : tableau S3), ce qui est similaire aux taux globaux calculés à partir de 11 études WES récentes sur les troubles du développement neurologique : 1,00/exome de bande, n = 2358 ; 0,82/exome de contrôle, n = 731 . Il n’y avait pas de différences significatives dans les taux de variants non synonymes de novo, homozygotes et hémizygotes (fichier supplémentaire 2 : tableau S3) ; le taux de gènes exprimés dans le cerveau affectés ; les scores de conservation de la valeur P phylogénétique (PhyloP) aux positions des variants ; les scores d’intolérance à la variation résiduelle (RVIS) ; et les scores de phénotypage du polymorphisme v2 (PolyPhen-2) (fichier supplémentaire 2 : tableau S2) entre les probands et les frères et sœurs.

Tableau 2 Variants rares uniques aux probands CDD

Nous avons trouvé un CNV génique de novo chez un proband (0,07/proband, tableau 2), ce qui est similaire aux taux précédemment rapportés pour les TSA , et aucun chez les frères et sœurs. Le CNV du proband est une délétion hétérozygote de 2 kb de la 3′UTR de OGDHL, qui code pour un composant d’un complexe protéique mitochondrial impliqué dans la neurodégénération. Un gène, SUPT20HL2, et deux familles de gènes, USP et BBS, sont affectés chez plus d’un proband CDD. Deux variants hémizygotes faux-sens ont été identifiés dans SUPT20HL2, qui code pour un facteur de transcription putatif mais pourrait être un pseudogène selon la base de données UniProtKB (http://www.uniprot.org/). Trois membres de la famille des gènes USP (ubiquitin-specific peptidase) sont affectés chez les probands CDD : USP9X (faux sens hémizygote), USP9Y (non-sens d’origine paternelle) et USP26 (faux sens hémizygote). Ils codent pour des enzymes de déubiquitination qui empêchent la dégradation des protéines. Deux membres de la famille de gènes du syndrome de Bardet-Biedel (BBS), qui est impliquée dans la ciliogenèse, présentent des variantes faux-sens de novo chez les probands CDD : BBS5 et BBS9. Bien que les variantes spécifiques modifiant les protéines identifiées chez les sujets CDD soient rares et n’aient pas été précédemment associées à la maladie, nous avons examiné la littérature et trouvé un certain chevauchement entre les gènes candidats CDD et les gènes potentiellement associés à d’autres troubles neurologiques (tableau 2).

Il n’y avait pas de variantes clairement délétères chez les probands CDD. Pour identifier les variants potentiellement pathogènes, nous avons considéré une combinaison de facteurs : (1) expression cérébrale positive, (2) score PhyloP ≥ 1,30 (P = 0,05 pour la conservation), (3) RVIS négatif (gène intolérant à la variation), et (4) classification PolyPhen-2 de missense probablement dommageable (ou n/a dû à une variante autre que missense). Sur les 47 gènes candidats CDD, 14 remplissaient tous ces critères : NRK, TBC1D8B, TRRAP, NAV2, OGDHL, ZNF236, PRKCSH, MTMR8, BCOR, SRPK3, USP9Y, KIAA2018, CXorf57 et ALG13 (tableau 2). Pour affiner cette liste, l’inspection des données de séquençage de l’Exome Aggregation Consortium (http://exac.broadinstitute.org/) a révélé que : (1) les variants de tous les gènes, à l’exception de NAV2, MTMR8 et ALG13, sont nouveaux ou trouvés au plus une fois dans l’ensemble de données et (2) parmi les 11 gènes restants, 4 font partie des 5 % les plus intolérants : TRRAP, ZNF236, BCOR et KIAA2018.

TRRAP (transformation/transcription domain-associated protein)affectée par un variant faux-sens de novo chez un proband CDD masculin ; elle code pour un composant des complexes histone acétyltransférase et est impliquée dans la transcription et la réparation de l’ADN. Il n’est pas associé à un trouble OMIM, mais des variants de novo ont été identifiés dans d’autres troubles neurologiques (tableau 2). ZNF236 (Zinc Finger Protein 236) est également affecté par un variant faux-sens de novo chez un proband mâle ; il peut être impliqué dans la régulation transcriptionnelle (UniProtKB) mais n’est pas associé à un trouble connu. BCOR (BCL6 Corepressor) est affecté par un variant faux-sens hémizygote chez un proband mâle CDD ; il code pour un corepresseur transcriptionnel. Il est associé à une microphtalmie syndromique, qui peut avoir la caractéristique d’une ID mais ne caractérise pas autrement le proband CDD et est généralement causé par des mutations tronquantes chez les femelles. KIAA2018 est affecté par une délétion homozygote d’un acide aminé chez un proband mâle CDD ; il est également connu sous le nom de USF3 (upstream transcription factor 3). Il n’est pas associé à un trouble OMIM, mais des variants de novo ont été identifiés dans d’autres troubles neurologiques (tableau 2). Il convient de noter que tous ces gènes candidats jouent un rôle ou peuvent être impliqués dans la transcription, ce qui caractérise également de nombreux gènes associés aux TSA.

En utilisant l’ensemble de données transcriptomiques BrainSpan exon-array humain, nous avons tracé le niveau d’expression médian des gènes candidats CDD en tant que groupe pour toutes les régions du cerveau disponibles des stades embryonnaires à la fin de l’âge adulte (n = 40 gènes représentés une fois dans l’ensemble de sondes de base, fichier supplémentaire 2 : tableau S4). Comme le montre le profil d’expression de la figure 1, les gènes candidats CDD sont plus fortement exprimés dans les régions non néocorticales que dans les régions néocorticales tout au long de la vie (fichier complémentaire 2 : tableau S5). De plus, on observe des niveaux d’expression croissants dans l’AMY, le STR et le HIP au cours des périodes 10 (1-6 ans) et 11 (6-12 ans), la fourchette qui englobe l’âge d’apparition des symptômes dans notre cohorte CDD.

Fig. 1

Niveaux d’expression médians des gènes candidats CDD (n = 40) par région cérébrale et période de temps (fichier supplémentaire 2 : tableau S5) en utilisant l’ensemble de données transcriptomiques BrainSpan exon-array humain . La ligne verticale foncée indique la naissance. Une intensité de signal transformée en log2 ≥ 6 dans au moins un échantillon est considérée comme une expression positive . AMY amygdale, CBC cortex cérébelleux, HIP hippocampe MD noyau médio-dorsal du thalamus, NCX néocortex, STR striatum

Compte tenu de cette observation, nous avons comparé la différence des niveaux d’expression médians entre les régions non néocorticales et néocorticales pour les gènes affectés par des variants non synonymes et synonymes chez les probands CDD, leurs frères et sœurs non affectés, et les probands TSA de la Simons Simplex Collection (SSC) avec et sans régression appariés par le sexe, l’âge au moment de l’évaluation, le QI et la gravité des symptômes de l’autisme (voir le fichier supplémentaire 1 : Informations complémentaires pour les détails de la sélection de la cohorte). Le profil d’expression des gènes candidats CDD est qualitativement distinct des autres ensembles de gènes, sauf qu’il est similaire au profil des gènes affectés par des variants non synonymes chez les probands SSC avec régression, même s’ils n’ont qu’un seul gène, NAV2, en commun (Fig. 2, Additional file 2 : Tables S4 et S6). La différence d’expression, non-néocorticale moins néocorticale, atteint une valeur positive maximale au stade mi-fœtal. Pour les gènes candidats CDD, cela se produit à la période six ; un test de permutation avec 100 000 itérations de 40 gènes choisis au hasard dans l’ensemble de données BrainSpan a confirmé la signification de cette expression différentielle (P = 0,0022). Nous avons étendu l’analyse à plusieurs autres ensembles de gènes, tels que ceux identifiés chez les probands du SSC et les frères et sœurs non affectés par des variants non synonymes, synonymes et LGD ; les gènes les plus significativement associés aux TSA par trois grandes études WES et CNV récentes ; et tous les gènes de l’ensemble de données BrainSpan. Le profil d’expression des gènes affectés par des variants non synonymes chez les probands CDD et les probands SSC avec régression est qualitativement distinct de ces autres ensembles également (fichier supplémentaire 1 : figure S1, fichier supplémentaire 2 : tableaux S4 et S6).

Fig. 2

Niveaux d’expression différentiels de divers ensembles de gènes. La différence des niveaux d’expression médians (régions cérébrales non néocorticales moins néocorticales) est indiquée pour les gènes affectés par des variants non synonymes ou synonymes chez les probands CDD, leurs frères et sœurs non affectés, les probands SSC avec régression et les probands SSC sans régression. Le nombre entre parenthèses indique le nombre de sujets ou de variants, et la ligne verticale sombre dans chaque panneau indique la naissance. Pour les gènes candidats potentiels du CDD, la différence atteint une valeur positive maximale à la période six (stades mi-fœtaux) ; la signification a été confirmée par un test de permutation avec 100 000 itérations de 40 gènes sélectionnés au hasard (P = 0,0022). CDD childhood disintegrative disorder, SSC Simons Simplex Collection

Nous avons également cherché à savoir si les gènes candidats CDD sont coexprimés entre eux. Sur les 40 gènes candidats, 11 sont coexprimés avec au moins un autre gène candidat dans toutes les régions du cerveau et toutes les périodes de temps, avec un coefficient de corrélation de Pearson r ≥ 0,7 (figure 3, fichier additionnel 2 : tableau S7). Il existe 23 connexions de ce type, pour une moyenne de 2,09 corrélations/gène et un coefficient moyen de 0,779. Un test de permutation avec 100 000 itérations de 40 gènes choisis au hasard dans l’ensemble de données BrainSpan a révélé que l’observation de 11 gènes avec au moins 2,09 corrélations/gène est significative (P = 0,036), tout comme l’observation de 11 gènes avec un coefficient de corrélation moyen d’au moins 0,779 (P = 0,019). Le fait d’atteindre les deux seuils est également significatif (P = 0,0059). Étant donné que les 11 gènes candidats CDD qui sont coexprimés les uns avec les autres ont une expression cérébrale positive selon l’ensemble de données BrainSpan, un test de permutation avec 100 000 itérations a également été effectué avec 40 gènes exprimés dans le cerveau sélectionnés au hasard dans BrainSpan. Si l’observation de 11 gènes présentant au moins 2,09 corrélations/gène n’est pas significative (P = 0,066), l’observation de 11 gènes présentant un coefficient de corrélation moyen d’au moins 0,779 est significative (P = 0,022), tout comme le fait d’atteindre les deux seuils (P = 0,011). La comparaison de l’ensemble de 11 gènes candidats CDD coexprimés avec l’ensemble restant de 29 qui ne sont pas coexprimés n’a révélé aucune différence significative entre le taux de gènes exprimés dans le cerveau, les scores PhyloP ou les scores PolyPhen-2 ; cependant, les gènes coexprimés sont significativement plus intolérants à la variation (RVIS moyen -1,42 contre -0,15, t(35) = -2,91, P = 0,0062, test t indépendant, bilatéral). L’analyse d’enrichissement de l’ontologie des gènes utilisant la base de données pour l’annotation, la visualisation et la découverte intégrée v6.8 (https://david.ncifcrf.gov/) pour l’ensemble des gènes candidats CDD, et le sous-ensemble de 11 gènes coexprimés n’a pas identifié d’enrichissement significatif des termes GO après correction de Benjamini-Hochberg des valeurs P.

Fig. 3

Analyse du réseau de coexpression des gènes. Onze des 40 gènes candidats CDD sont coexprimés avec au moins un autre gène candidat dans toutes les régions du cerveau et toutes les périodes de temps avec un coefficient de corrélation de Pearson r ≥ 0,7 (fichier supplémentaire 2 : tableau S7), une moyenne de 2,09 corrélations/gène (P = 0,036) et un coefficient moyen de 0,779 (P = 0,019, test de permutation avec 100 000 itérations de 40 gènes sélectionnés au hasard). Les corrélations positives sont indiquées en bleu, et les corrélations négatives en rouge. Plus l’amplitude du coefficient est grande, plus les bords sont larges et sombres. La taille d’un nœud est proportionnelle au nombre d’arêtes qu’il possède

Systèmes neuronaux

Vu l’universalité des déficits sociaux dans les TSA, le dysfonctionnement des systèmes cérébraux subsistant dans la perception sociale, y compris la perception des visages, est un point clé de la recherche sur les TSA. Les stimuli visuels des visages et des maisons activent et dissocient de manière fiable les systèmes impliqués dans le traitement de l’information socio-émotionnelle (visages craintifs) et non-socio-émotionnelle (maisons). Nous avons étudié quatre cohortes : CDD (n = 7), LFASD (n = 7), HFASD (n = 14), et TD (n = 19). Même si les personnes atteintes de TALF sont plus nombreuses que celles atteintes de TDC, la taille de notre échantillon a été limitée par la difficulté d’obtenir des données de neuroimagerie (et d’oculométrie) de haute qualité chez les sujets à faible niveau de fonctionnement. Cela dit, à notre connaissance, il s’agit de la toute première présentation de données d’IRMf non sédatées provenant de personnes atteintes de TSA et de DI marquée.

Il n’y avait pas de différences significatives entre les quatre cohortes en ce qui concerne le sexe, l’âge, le volume intracrânien et le mouvement de la tête dans le scanner. Les groupes CDD et LFASD n’étaient pas non plus significativement différents en fonction du QI et de la sévérité de l’autisme, et les groupes HFASD et TD n’étaient pas significativement différents en fonction du QI (fichier supplémentaire 2 : tableaux S8 et S9). Tout d’abord, nous avons utilisé un échantillon de découverte de 12 de nos 19 sujets atteints de TD dans une analyse du cerveau entier pour une localisation indépendante des régions d’intérêt impliquées dans le traitement des visages par rapport aux maisons. La figure 4a illustre les régions du cortex occipitotemporal ventrolatéral où les sujets atteints de TD présentaient une activation significative des visages > des maisons (fichier supplémentaire 2 : tableau S10). Ces régions comprenaient les emplacements prévus des nœuds bien connus du réseau occipitotemporal sensible aux visages, notamment l’aire fusiforme du visage et l’aire occipitale du visage . Comme le montrent la figure 4b et le fichier supplémentaire 2 : tableau S11, l’extraction du pourcentage moyen de changement de signal (visages > maisons) pour chacun des quatre groupes a indiqué une absence de différences entre les groupes dans la réponse aux visages par rapport aux maisons dans ces régions d’intérêt définies de manière indépendante lors de la comparaison des groupes TD:validation et HFASD et lors de la comparaison des groupes LFASD et CDD. La réponse aux visages > maisons dans le groupe CDD n’était pas significativement supérieure à zéro, ce qui suggère un manque global de sensibilité aux visages dans le réseau occipito-temporal sensible aux visages dans son ensemble. Il existe un constat bien établi d’hypoactivation aux visages (par rapport aux maisons) dans le gyrus fusiforme moyen droit dans le HFASD par rapport au TD . Nous avons pu reproduire ce résultat dans nos cohortes. Cependant, la comparaison de l’activité des visages > maisons dans le CDD par rapport au TD n’a révélé aucune différence significative, tout comme la comparaison des groupes LFASD et TD (fichier supplémentaire 1 : Figure S2 et fichier supplémentaire 2 : Tableau S12).

Fig. 4

Régions cérébrales d’intérêt (ROI) impliquées dans le traitement des stimuli visuels socio-émotionnels (visage craintif) par rapport aux stimuli visuels non socio-émotionnels (maison). a La carte cérébrale de couleur verte indique les régions d’activation significative des visages >maisons dans un échantillon de découverte de 12 sujets TD (Z >3,09, cerveau entier corrigé au niveau du cluster P <0,05). b Ces ROIs définies indépendamment ont ensuite été utilisées pour des comparaisons entre les quatre cohortes restantes, un échantillon TD:validation (n = 7), HFASD (n = 14), LFASD (n = 7), et CDD (n = 7). Le graphique à barres indique le % moyen de changement de signal (faces > maisons) pour chaque cohorte. Les différences entre les groupes n’étaient pas significatives lorsqu’on comparait les groupes TD:validation et HFASD et lorsqu’on comparait les groupes LFASD et CDD . La réponse des visages > maisons dans le groupe CDD n’était pas significativement supérieure à zéro. Les barres d’erreur indiquent l’erreur standard de la moyenne. Toutes les valeurs P ont été calculées par un test t indépendant et sont bilatérales. FFG gyrus fusiforme, L gauche, LOC cortex occipital latéral, MTG gyrus temporal moyen, R droit

Compte tenu du manque possible de sensibilité aux visages dans le cortex occipitotemporal ventrolatéral dans le CDD, nous avons ensuite effectué une évaluation du cerveau entier des sujets CDD pour localiser les substrats neuroanatomiques de la perception des visages chez ces personnes. Comme l’illustre la figure 5a, les sujets atteints de TDC présentent une activité des maisons des visages > dans le gyrus frontal moyen, le gyrus précentral, le caudate (striatum), le thalamus, l’hippocampe et le cervelet (fichier supplémentaire 2 : tableau S13). Ces régions se chevauchent avec les régions du cerveau déterminées comme ayant les niveaux les plus élevés d’expression des gènes candidats CDD (Fig. 1). Comme le montrent la figure 5b et le fichier supplémentaire 2 : tableau S14, la comparaison du pourcentage moyen de changement de signal (faces > maisons) de ces régions d’intérêt a révélé une différence significative entre le groupe CDD et le groupe HFASD , mais aucune différence significative entre le groupe CDD et le groupe LFASD . Le groupe LFASD a montré un phénotype intermédiaire à celui des groupes HFASD et CDD (Fig. 5b).

Fig. 5

Analyse IRMF du cerveau entier par la DDC. a La carte cérébrale de couleur rouge indique les régions d’activation significative des faces > maisons chez les sujets CDD (Z > 3.09, cerveau entier corrigé au niveau du groupe P < 0,05). b Le graphique à barres indique le % moyen de changement de signal (visages > maisons) dans ces régions pour chaque cohorte : TD:découverte (n = 12), TD:validation (n = 7), HFASD (n = 14), LFASD (n = 7) et CDD (n = 7). La cohorte CDD différait significativement de la cohorte HFASD mais pas de la cohorte LFASD. Les barres d’erreur indiquent l’erreur standard de la moyenne. Toutes les valeurs P ont été calculées par un test t indépendant et sont bilatérales. MFG gyrus frontal moyen, PG gyrus précentral

Comportement oculomoteur

Nous avons recueilli des données oculomotrices pour quantifier le phénotype social de nos quatre cohortes lorsqu’elles regardaient des visages émotionnels . Comme le montre le fichier supplémentaire 2 : tableaux S15 et S16, les groupes n’étaient pas significativement différents selon le sexe, l’âge et la durée totale de fixation sur l’image. Les groupes CDD et LFASD n’étaient pas non plus significativement différents selon le QI et la sévérité de l’autisme, et les groupes HFASD et TD n’étaient pas significativement différents selon le QI. Comme le montre la figure 6, nous reproduisons les résultats antérieurs d’une diminution de la fixation sur les yeux et d’une augmentation de la fixation sur la bouche dans le groupe HFASD par rapport au groupe TD. Cependant, alors que le pourcentage de temps passé par les sujets atteints de TALF à regarder les yeux ne diffère pas significativement du groupe TALF, les sujets CDD fixent les yeux significativement plus que le groupe TALF. Comparés les uns aux autres, les sujets CDD et LFASD ne diffèrent pas significativement dans le temps passé à regarder les yeux. Comme pour les résultats de l’IRMf (Fig. 5b), le groupe LFASD a montré un phénotype intermédiaire à celui des groupes HFASD et CDD (rapport œil-bouche, fichier additionnel 2 : tableau S16).

Fig. 6

Analyse comportementale par le suivi des yeux. Les barres jaunes et vertes du graphique représentent le % moyen de temps passé à fixer (axe y) les yeux et la bouche des visages, respectivement, par cohorte (axe x) : TD (n = 14), HFASD (n = 32), LFASD (n = 7), CDD (n = 5). Les cartes thermiques du regard illustrent les données du regard au niveau du groupe superposées à l’une des images que les sujets ont regardées. Par rapport aux sujets TD, les sujets HFASD montrent une diminution de la fixation sur les yeux et une augmentation de la fixation sur la bouche. Le pourcentage de temps passé par les sujets atteints de TALF à regarder les yeux ne diffère pas de celui des sujets atteints de TSAH, mais les sujets atteints de TDC fixent les yeux beaucoup plus que les sujets atteints de TSAH. Les barres d’erreur indiquent l’erreur standard de la moyenne. Toutes les valeurs P ont été calculées par un test t indépendant et sont bilatérales

.

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée.