Quantification de la relation séquence-activité de l’ARN par prédiction de structure tridimensionnelle

Dans un premier temps, nous avons modélisé la structure d’une famille d’ARN avec une grammaire de graphes afin d’identifier les séquences qui en font partie. Plusieurs autres méthodes de modélisation ont été développées, telles que des grammaires stochastiques hors-contexte, des modèles de covariance, des profils de structures secondaires et des réseaux de contraintes. Ces méthodes de modélisation se basent sur la structure secondaire classique comparativement à nos grammaires de graphes qui se basent sur les motifs cycliques de nucléotides. Pour exemplifier notre modèle, nous avons utilisé la boucle E du ribosome qui contient le motif Sarcin-Ricin qui a été largement étudié depuis sa découverte par cristallographie aux rayons X au début des années 90. Nous avons construit une grammaire de graphes pour la structure du motif Sarcin-Ricin et avons dérivé toutes les séquences qui peuvent s’y replier. La pertinence biologique de ces séquences a été confirmée par une comparaison des séquences d’un alignement de plus de 800 séquences ribosomiques bactériennes. Cette comparaison a soulevée des alignements alternatifs pour quelques unes des séquences que nous avons supportés par des prédictions de structures secondaires et tertiaires. Les motifs cycliques de nucléotides ont été observés par les membres de notre laboratoire dans l'ARN dont la structure tertiaire a été résolue expérimentalement. Une étude des séquences et des structures tertiaires de chaque cycle composant la structure du Sarcin-Ricin a révélé que l'espace des séquences dépend grandement des interactions entre tous les nucléotides à proximité dans l’espace tridimensionnel, c’est-à-dire pas uniquement entre deux paires de bases adjacentes. Le nombre de séquences générées par la grammaire de graphes est plus petit que ceux des méthodes basées sur la structure secondaire classique. Cela suggère l’importance du contexte pour la relation entre la séquence et la structure, d’où l’utilisation d’une grammaire de graphes contextuelle plus expressive que les grammaires hors-contexte. Les grammaires de graphes que nous avons développées ne tiennent compte que de la structure tertiaire et négligent les interactions de groupes chimiques spécifiques avec des éléments extra-moléculaires, comme d’autres macromolécules ou ligands. Dans un deuxième temps et pour tenir compte de ces interactions, nous avons développé un modèle qui tient compte de la position des groupes chimiques à la surface des structures tertiaires. L’hypothèse étant que les groupes chimiques à des positions conservées dans des séquences prédéterminées actives, qui sont déplacés dans des séquences inactives pour une fonction précise, ont de plus grandes chances d’être impliqués dans des interactions avec des facteurs. En poursuivant avec l’exemple de la boucle E, nous avons cherché les groupes de cette boucle qui pourraient être impliqués dans des interactions avec des facteurs d'élongation. Une fois les groupes identifiés, on peut prédire par modélisation tridimensionnelle les séquences qui positionnent correctement ces groupes dans leurs structures tertiaires. Il existe quelques modèles pour adresser ce problème, telles que des descripteurs de molécules, des matrices d’adjacences de nucléotides et ceux basé sur la thermodynamique. Cependant, tous ces modèles utilisent une représentation trop simplifiée de la structure d’ARN, ce qui limite leur applicabilité. Nous avons appliqué notre modèle sur les structures tertiaires d’un ensemble de variants d’une séquence d’une instance du Sarcin-Ricin d’un ribosome bactérien. L’équipe de Wool à l’université de Chicago a déjà étudié cette instance expérimentalement en testant la viabilité de 12 variants. Ils ont déterminé 4 variants viables et 8 létaux. Nous avons utilisé cet ensemble de 12 séquences pour l’entraînement de notre modèle et nous avons déterminé un ensemble de propriétés essentielles à leur fonction biologique. Pour chaque variant de l’ensemble d’entraînement nous avons construit des modèles de structures tertiaires. Nous avons ensuite mesuré les charges partielles des atomes exposés sur la surface et encodé cette information dans des vecteurs. Nous avons utilisé l’analyse des composantes principales pour transformer les vecteurs en un ensemble de variables non corrélées, qu’on appelle les composantes principales. En utilisant la distance Euclidienne pondérée et l’algorithme du plus proche voisin, nous avons appliqué la technique du « Leave-One-Out Cross-Validation » pour choisir les meilleurs paramètres pour prédire l’activité d’une nouvelle séquence en la faisant correspondre à ces composantes principales. Finalement, nous avons confirmé le pouvoir prédictif du modèle à l’aide d’un nouvel ensemble de 8 variants dont la viabilité à été vérifiée expérimentalement dans notre laboratoire. En conclusion, les grammaires de graphes permettent de modéliser la relation entre la séquence et la structure d’un élément structural d’ARN, comme la boucle E contenant le motif Sarcin-Ricin du ribosome. Les applications vont de la correction à l’aide à l'alignement de séquences jusqu’au design de séquences ayant une structure prédéterminée. Nous avons également développé un modèle pour tenir compte des interactions spécifiques liées à une fonction biologique donnée, soit avec des facteurs environnants. Notre modèle est basé sur la conservation de l'exposition des groupes chimiques qui sont impliqués dans ces interactions. Ce modèle nous a permis de prédire l’activité biologique d’un ensemble de variants de la boucle E du ribosome qui se lie à des facteurs d'élongation.

Prévalence d’excrétion de Giardia et Cryptosporidium chez les humains, les animaux domestiques et les lémuriens de l’écosystème du Parc National de Ranomafana, Madagascar

L’augmentation des interactions entre humains et animaux sauvages en lisière des habitats naturels pourrait faciliter la transmission d’agents pathogènes entre les humains et les différentes espèces animales d’un écosystème et ainsi favoriser l’émergence de maladies. Nous avons effectué une étude transversale portant sur l’infection par Giardia et Cryptosporidium chez les humains, les animaux domestiques, les rongeurs et les lémuriens au sein de l’écosystème de Ranomafana, Madagascar. Des échantillons de fèces ont étés collectés de manière non invasive chez des personnes volontaires, des mammifères domestiques et des rongeurs introduits habitant trois villages situés en lisière du Parc National de Ranomafana (PNR) ainsi que quatre espèces de lémuriens (Propithecus edwardsii, Prolemur simus, Eulemur rubriventer et Microcebus rufus) du PNR. Des analyses coproscopiques par la technique d’immunofluorescence directe ont été réalisées afin de détecter la présence de Cryptosporidium et Giardia. Leur prévalence a été estimée et certaines variables reliées à l’infection par les parasites ont été identifiées. Cryptosporidium et Giardia ont été détectés avec une prévalence estimée à 22,9 % et 13,6 % respectivement chez les humains. La prévalence de ces deux parasites variait de 0 % à 60 % chez les animaux domestiques et les rongeurs au sein des villages. L’espèce hôte, l’âge ainsi que la co-infection par un autre protozoaire sont les seules variables associées à l’infection par Cryptosporidium et Giardia dans cet écosystème tandis qu’aucune association avec une coinfection par un ordre de nématode n’a été détecté. De plus, Cryptosporidium a été détecté chez 10,5 % des lémuriens du PNR. Cette étude documente pour la première fois la présence de Cryptosporidium chez deux espèces de lémuriens du PNR. Par contre, Giardia n’a pas été détecté dans les échantillons issus de lémuriens du PNR.

Rythme veille-sommeil et dimensions cliniques dans le trouble de personnalité limite à l’adolescence

Cette thèse examine le rythme veille-sommeil et son association avec l’instabilité émotionnelle, l’agressivité et l’impulsivité dans le trouble de personnalité limite (TPL) à l’adolescence. Dans un premier temps, la revue de la littérature sur les perturbations objectives du sommeil dans le TPL a mis en lumière plusieurs difficultés similaires, évaluées par polysomnographie, à celles observées dans la dépression adulte. De 1980 à 2010, aucune recherche n’a examiné le rythme veille-sommeil, aucune n’a étudié les adolescents TPL et plusieurs n’ont pas contrôlé l’état dépressif comme facteur de confusion. De ce constat, il s’avérait pertinent de mener une étude sur le rythme veille-sommeil dans le TPL à l’adolescence en l’absence de dépression co-occurrente. L’adolescence comportant plusieurs caractéristiques physiologiques, psychologiques et sociales, tenir compte des aspects développementaux était essentiel. Dans un second temps, un protocole de recherche fût mis en place à la Clinique des troubles de l’humeur et le recrutement a été réalisé auprès d’adolescents souffrant d’un TPL et sans état dépressif actuel. Ils devaient porter pendant plus de neuf jours (période comprenant deux fins de semaine) un actigraphe, appareil non invasif évaluant l’alternance veille-sommeil dans l’environnement naturel. L’abandon précoce au traitement étant prévalent chez les patients TPL, la fiabilité de l’étude a été examinée afin de déterminer les raisons favorisant et celles nuisant au recrutement et à la collecte des données. La réflexion sur les aspects méthodologiques de l’étude actigraphique a permis d’expliquer les limites de ce type de protocole. Dans un troisième temps, le rythme veille-sommeil des adolescents TPL (n=18) a été caractérisé et comparé à celui des jeunes ayant un trouble bipolaire (n=6), trouble psychiatrique partageant plusieurs manifestations communes avec le TPL, et à celui des adolescents sans trouble de santé mentale (n=20). Les résultats suggèrent que l’adolescent TPL passe plus de temps en éveil durant la période de repos que les jeunes appartenant aux deux autres groupes. De plus, les adolescents TPL présentent une plus grande variabilité inter journalière des heures de lever et du temps total de sommeil que les autres adolescents. Ils se réveillent une heure de plus, et dorment donc une heure supplémentaire, que les adolescents sans trouble mental lors des journées sans routine. Dans un quatrième temps, les analyses corrélationnelles entre les données actigraphiques et les scores aux questionnaires auto-rapportés évaluant l’instabilité émotionnelle, l’agressivité et l’impulsivité suggèrent que plus l’adolescent TPL passe du temps éveillé alors qu’il est au lit, plus il déclare présenter des comportements agressifs, surtout physiques, durant le jour. En résumé, cette thèse contribue à la littérature scientifique en explorant pour la première fois le rythme veille-sommeil et son lien avec les manifestations symptomatiques dans le TPL à l’adolescence. Les résultats suggèrent fortement l’importance d’évaluer et de traiter les problèmes du rythme veille-sommeil que présentent ces jeunes lors de la prise en charge.

Analyse des processus de dérive lors de la gravure profonde du silicium dans des plasmas SF6 et C4F8

L’objectif de ce mémoire de maîtrise est de développer des outils de diagnostics non-invasifs et de caractériser in-situ les dérives de procédé dans un réacteur industriel utilisé en production pour la gravure profonde du silicium par le procédé Bosch. Ce dernier repose sur l’alternance d’un plasma de SF6 pour la gravure isotrope du Si et d’un plasma de C4F8 pour la passivation des parois dans l’optique d’obtenir des tranchées profondes et étroites. Dans un premier temps, nous avons installé une sonde courant-tension sur la ligne de transmission du signal rf au porte-substrat pour l’étude de son impédance caractéristique et un spectromètre optique pour l’étude de l’émission optique du plasma. Nos travaux ont montré que l’évolution temporelle de l’impédance constitue un excellent moyen pour identifier des changements dans la dynamique du procédé, notamment une gravure complète de la photorésine. De plus, à partir des spectres d’émission, nous avons pu montrer que des produits carbonés sont libérés du substrat et des parois lors de l’alternance passivation/gravure et que ceux-ci modifient considérablement la concentration de fluor atomique dans le plasma. Dans un second temps, nous avons développé un réacteur à « substrat-tournant » pour l’analyse in-situ des interactions plasma-parois dans le procédé Bosch. Nos travaux sur ce réacteur visaient à caractériser par spectrométrie de masse l’évolution temporelle des populations de neutres réactifs et d’ions positifs. Dans les conditions opératoires étudiées, le SF6 se dissocie à près de 45% alors que le degré de dissociation du C4F8 atteint 70%. Le SF6 est avant tout dissocié en F et SF3 et l’ion dominant est le SF3+ alors que le C4F8 est fragmenté en CF, CF3 et CF4 et nous mesurons plusieurs ions significatifs. Dans les deux cas, la chaîne de dissociation demeure loin d’être complète. Nous avons noté une désorption importante des parois de CF4 lors du passage du cycle de passivation au cycle de gravure. Un modèle d’interactions plasmas-parois est proposé pour expliquer cette observation.

Evaluating perceptual maps of asymmetries for gait symmetry quantification and pathology detection

Le mouvement de la marche est un processus essentiel de l'activité humaine et aussi le résultat de nombreuses interactions collaboratives entre les systèmes neurologiques, articulaires et musculo-squelettiques fonctionnant ensemble efficacement. Ceci explique pourquoi une analyse de la marche est aujourd'hui de plus en plus utilisée pour le diagnostic (et aussi la prévention) de différents types de maladies (neurologiques, musculaires, orthopédique, etc.). Ce rapport présente une nouvelle méthode pour visualiser rapidement les différentes parties du corps humain liées à une possible asymétrie (temporellement invariante par translation) existant dans la démarche d'un patient pour une possible utilisation clinique quotidienne. L'objectif est de fournir une méthode à la fois facile et peu dispendieuse permettant la mesure et l'affichage visuel, d'une manière intuitive et perceptive, des différentes parties asymétriques d'une démarche. La méthode proposée repose sur l'utilisation d'un capteur de profondeur peu dispendieux (la Kinect) qui est très bien adaptée pour un diagnostique rapide effectué dans de petites salles médicales car ce capteur est d'une part facile à installer et ne nécessitant aucun marqueur. L'algorithme que nous allons présenter est basé sur le fait que la marche saine possède des propriétés de symétrie (relativement à une invariance temporelle) dans le plan coronal.

Dynamic magnetic resonance spectroscopy of phosphate energetics during muscle exercise and recovery

Entailing of phosphorus exchanges in most bio-chemicals as a key factor in disease, increases researcher’s interest to develop the technologies capable of detecting this metabolite. Phosphorus magnetic resonance spectroscopy is able to detect key metabolites in a non-invasive manner. Particularly, it offers the ability to measure the dynamic rate of phosphocreatine(PCr) degeneration through the exercise and recovery. This metric as a valid indication of mitochondrial oxidative metabolism in muscle, differentiate between normal and pathological state. To do magnetic resonance imaging and spectroscopy, clinical research tools provide a wide variety of anatomical and functional contrasts, however they are typically restricted to the tissues containing water or hydrogen atoms and they are still blind to the biochemicals of other atoms of interests. Through this project we intended to obtain the phosphorus spectrum in human body – specificadenerativelly in muscle – using 31P spectroscopy. To do so a double loop RF surface coil, tuned to phosphorus frequency, is designed and fabricated using bench work facilities and then validated through in vitro spectroscopy using 3 Tesla Siemens scanner. We acquired in vitro as well as in vivo phosphorus spectrum in a 100 mM potassium phosphate phantom and human calf muscle in rest-exercise-recovery phase in a 3T MR scanner. The spectrum demonstrates the main constituent in high-energy phosphate metabolism. We also observed the dynamic variation of PCr for five young healthy subjects who performed planter flexions using resistance band during exercise and recovery. The took steps in this project pave the way for future application of spectroscopic quantification of phosphate metabolism in patients affected by carotid artery disease as well as in age-matched control subjects.

Influence of plant species and richness on pollutant removal in treatment wetlands

Les marais filtrants artificiels sont des écosystèmes recréés par l’homme dans le but d’optimiser l’épuration des eaux usées. Lors de la sélection d’espèces végétales pour la mise en place de ces marais filtrants, l’utilisation d’une polyculture ainsi que d’espèces indigènes non invasives est de plus en plus recommandée. Néanmoins, la plupart des marais filtrants existants sont des monocultures utilisant des plantes envahissantes, probablement à cause du manque d’évidences scientifiques sur les avantages de la diversité végétale et de la performance des espèces locales. Ainsi, les questions de recherche autour desquelles s’oriente ma thèse sont: Les polycultures présentent-elles un potentiel épuratoire aussi ou plus grand que les monocultures, et une espèce indigène est-elle aussi efficace et performante qu’une espèce exotique envahissante dans des marais filtrants ? Trois expériences ont été conduites afin de répondre à ces questions. J’ai d’abord testé l’influence de la richesse végétale sur l’élimination des polluants en deux dispositifs expérimentaux: 1) comparant deux espèces de plantes émergentes en monoculture ou combinées séquentiellement, et 2) évaluant la performance de quatre espèces flottantes plantées en monoculture par rapport à des associations de deux (avec toutes les combinaisons possibles) et de quatre espèces. Une troisième expérience a été réalisée afin de comparer l’efficacité épuratoire de l’haplotype européen envahissant du roseau commun (Phragmites australis) et de la sous-espèce locale non-invasive (P. australis subsp. americanus). La composition en espèces végétales a produit un effet notable sur la performance des marais filtrants. La comparaison des performances en mono- et en polyculture n’a pas permis de démontrer clairement les avantages de la diversité végétale pour l’élimination des polluants dans les marais filtrants. Toutefois, les marais filtrants plantés avec une combinaison d’espèces étaient aussi efficaces que les monocultures des espèces les plus performantes. La comparaison entre les deux sous-espèces de P. australis indiquent que la sous-espèce indigène pourrait remplacer le roseau exotique envahissant, évitant ainsi les potentiels risques environnementaux sans toutefois compromettre l’efficacité du traitement. Les résultats prometteurs de la sous-espèce indigène de P. australis doivent encore être testés dans des expériences à grande échelle avant d’utiliser largement cette espèce dans les marais filtrants. Nos résultats suggèrent que, dans des conditions où la performance des macrophytes disponibles est inconnue ou ne peut être déterminée, l’utilisation d’une combinaison d’espèces présente les meilleures chances d’accomplir le plus haut niveau possible d’élimination de polluants. De plus, même si la diversité végétale ne présente pas un avantage mesurable en termes d’efficacité épuratoire, celle-ci améliore la résilience des marais filtrants et leur résistance aux stress et aux maladies.

Hemodynamic changes during posterior epilepsies: an eeg-fnirs study

Posterior epilepsies are relatively rare, mainly suspected clinically by the presence of visual auras. Functional near-infrared spectroscopy (fNIRS) is an emerging non-invasive imaging technique that has the potential to monitor hemodynamic changes during epileptic activity. Combined with electroencephalography (EEG), 9 patients with posterior epilepsies were recorded using EEG-fNIRS with large sampling (19 EEG electrodes and over 100 fNIRS channels). Spikes and seizures were carefully marked on EEG traces, and convolved with a standard hemodynamic response function for general linear model (GLM) analysis. GLM results for seizures (in 3 patients) and spikes (7 patients) were broadly sensitive to the epileptic focus in 7/9 patients, and specific in 5/9 patients with fNIRS deoxyhemoglobin responses lateralized to the correct lobe, and to plausible locations within the occipital or parietal lobes. This work provides evidence that EEG-fNIRS is a sensitive technique for monitoring posterior epileptic activity.

fNIRS-EEG study of focal interictal epileptiform discharges

Functional near-infrared spectroscopy (fNIRS) acquired with electroencephalography (EEG) is a relatively new non-invasive neuroimaging technique with potential for long term monitoring of the epileptic brain. Simultaneous EEG-fNIRS recording allows the spatio-temporal reconstruction of the hemodynamic response in terms of the concentration changes in oxy-hemoglobin (HbO) and deoxy-hemoglobin (HbR) associated with recorded epileptic events such as interictal epileptic discharges (IEDs) or seizures. While most previous studies investigating fNIRS in epilepsy had limitations due to restricted spatial coverage and small sample sizes, this work includes a sufficiently large number of channels to provide an extensive bilateral coverage of the surface of the brain for a sample size of 40 patients with focal epilepsies. Topographic maps of significant activations due to each IED type were generated in four different views (dorsal, frontal, left and right) and were compared with the epileptic focus previously identified by an epileptologist. After excluding 5 patients due to the absence of IEDs and 6 more with mesial temporal foci too deep for fNIRS, we report that significant HbR (respectively HbO) concentration changes corresponding to IEDs were observed in 62% (resp. 38%) of patients with neocortical epilepsies. This HbR/HbO response was most significant in the epileptic focus region among all the activations in 28%/21% of patients.

Les effets de la stimulation électrique transcrânienne à courant direct appliquée au cortex somatosensoriel primaire sur la perception vibrotactile

La stimulation électrique transcrânienne à courant direct (tDCS) est une technique non invasive de neuromodulation qui modifie l’excitabilité corticale via deux grosses électrodes de surface. Les effets dépendent de la polarité du courant, anodique = augmentation de l’excitabilité corticale et cathodique = diminution. Chez l’humain, il n’existe pas de consensus sur des effets de la tDCS appliquée au cortex somatosensoriel primaire (S1) sur la perception somesthésique. Nous avons étudié la perception vibrotactile (20 Hz, amplitudes variées) sur le majeur avant, pendant et après la tDCS appliquée au S1 controlatéral (anodale, a; cathodale, c; sham, s). Notre hypothèse « shift-gain » a prédit une diminution des seuils de détection et de discrimination pour la tDCS-a (déplacement vers la gauche de la courbe stimulus-réponse et une augmentation de sa pente). On attendait les effets opposés avec la tDCS-c, soit une augmentation des seuils (déplacement à droite et diminution de la pente). Chez la majorité des participants, des diminutions des seuils ont été observées pendant et immédiatement suivant la tDCS-a (1 mA, 20 min) en comparaison à la stimulation sham. Les effets n’étaient plus présents 30 min plus tard. Une diminution du seuil de discrimination a également été observée pendant, mais non après la tDCS-c (aucun effet pour détection). Nos résultats supportent notre hypothèse, uniquement pour la tDCS-a. Une suite logique serait d’étudier si des séances répétées de tDCS-a mènent à des améliorations durables sur la perception tactile. Ceci serait bénéfique pour la réadaptation sensorielle (ex. suite à un accident vasculaire cérébral).

Modélisation de l’activité électrique des oreillettes avant et après ablation par cathéter

- Réalisé au centre de recherche de l'hospital du Sacré-Coeur de Montréal. - Programme conjoint entre Université de Montréal et École Polytechnique de Montréal.

Scoliosis curve type classification from 3D trunk image

Adolescent idiopathic scoliosis (AIS) is a deformity of the spine manifested by asymmetry and deformities of the external surface of the trunk. Classification of scoliosis deformities according to curve type is used to plan management of scoliosis patients. Currently, scoliosis curve type is determined based on X-ray exam. However, cumulative exposure to X-rays radiation significantly increases the risk for certain cancer. In this paper, we propose a robust system that can classify the scoliosis curve type from non invasive acquisition of 3D trunk surface of the patients. The 3D image of the trunk is divided into patches and local geometric descriptors characterizing the surface of the back are computed from each patch and forming the features. We perform the reduction of the dimensionality by using Principal Component Analysis and 53 components were retained. In this work a multi-class classifier is built with Least-squares support vector machine (LS-SVM) which is a kernel classifier. For this study, a new kernel was designed in order to achieve a robust classifier in comparison with polynomial and Gaussian kernel. The proposed system was validated using data of 103 patients with different scoliosis curve types diagnosed and classified by an orthopedic surgeon from the X-ray images. The average rate of successful classification was 93.3% with a better rate of prediction for the major thoracic and lumbar/thoracolumbar types.

Evaluation of Reducibility of Trunk Asymmetry in Lateral Bending

The value of the lateral bending test is important in the assessment of spinal curve mobility and prediction of surgical outcome in the treatment of adolescent idiopathic scoliosis (AIS). However, radiographic bending tests are unable to assess the reducibility of trunk asymmetry. This study aims to exploit surface topography measurement in order to evaluate the changes in shape of the trunk (a) between bending and neutral standing positions, and (b) between standing pre- and post-operative visits, in a cohort of adolescents with AIS having undergone surgical correction; and to correlate the differences measured in cases (a) and (b). Our cohort includes 13 patients with right thoracic AIS. Each patient had their 3D trunk surface digitized with a multi-head InSpeck system in standing posture (at the pre-op and post-op visits) and in maximum voluntary right and left bending (at the pre-op visit). We developed a novel trunk shape analysis method which produces a set of inclined trunk cross-sections allowing comparison between different postures. Two asymmetry indices, trunk rotation (TR) and back surface rotation (BSR), were computed in all cases and a statistical analysis was performed. Our correlation study (Pearson test) showed fair correlations in most cases between the changes in side-bending and those post-surgery, with the strongest relationship (p-value < 0.01) when combining the TR measurements from both bendings. These results provide evidence that the bending test can be used to assess trunk asymmetry reducibility. The proposed approach could provide a non-invasive trunk asymmetry reducibility test for routine clinical use in AIS surgery planning.

Noninvasive Clinical Assessment of Trunk Deformities Associated With Scoliosis

Besides the spinal deformity, scoliosis modifies notably the general appearance of the trunk resulting in trunk rotation, imbalance, and asymmetries that constitutes patients' major concern. Existing classifications of scoliosis, based on the type of spinal curve as depicted on radiographs, are currently used to guide treatment strategies. Unfortunately, even though a perfect correction of the spinal curve is achieved, some trunk deformities remain, making patients dissatisfied with the treatment received. The purpose of this study is to identify possible shape patterns of trunk surface deformity associated with scoliosis. First, trunk surface is represented by a multivariate functional trunk shape descriptor based on 3-D clinical measurements computed on cross sections of the trunk. Then, the classical formulation of hierarchical clustering is adapted to the case of multivariate functional data and applied to a set of 236 trunk surface 3-D reconstructions. The highest internal validity is obtained when considering 11 clusters that explain up to 65% of the variance in our dataset. Our clustering result shows a concordance with the radiographic classification of spinal curves in 68% of the cases. As opposed to radiographic evaluation, the trunk descriptor is 3-D and its functional nature offers a compact and elegant description of not only the type, but also the severity and extent of the trunk surface deformity along the trunk length. In future work, new management strategies based on the resulting trunk shape patterns could be thought of in order to improve the esthetic outcome after treatment, and thus patients satisfaction.

Prediction of scoliosis curve type based on the analysis of trunk surface topography

Scoliosis treatment strategy is generally chosen according to the severity and type of the spinal curve. Currently, the curve type is determined from X-rays whose acquisition can be harmful for the patient. We propose in this paper a system that can predict the scoliosis curve type based on the analysis of the surface of the trunk. The latter is acquired and reconstructed in 3D using a non invasive multi-head digitizing system. The deformity is described by the back surface rotation, measured on several cross-sections of the trunk. A classifier composed of three support vector machines was trained and tested using the data of 97 patients with scoliosis. A prediction rate of 72.2% was obtained, showing that the use of the trunk surface for a high-level scoliosis classification is feasible and promising.