Estudi de la relació estructura-funció de les proteïnes CPT1

En aquest treball s’ha analitzat la relació estructura-funció dels enzims CPT1, o Carnitina palmitoïltransferasa 1, que catalitza la reacció de transesterificació dels àcids grassos de cadena llarga a acil-carnitines, per tal que puguin accedir a la matriu mitocondrial i ser oxidats. Aquest enzim es troba estrictament regulat per malonil-CoA, primer intermediari de la síntesi d’àcids grassos, establint-se així una regulació coordinada entre la formació i la degradació de grasses. S’han estudiat els tres isotips de CPT1 descrits fins al moment: CPT1A, CPT1B i CPT1C. Mitjançant l’expressió heteròloga de mutants de CPT1A de rata i CPT1B de porc en el llevat P. pastoris, s’ha estudiat l’efecte sobre la inhibició per malonil-CoA de petits canvis en la seva estructura, per tal de trobar una relació entre la seva funció enzimàtica i la disposició conformacional de la proteïna. Segons els resultats obtinguts, el residu Glu590 de CPT1A de rata estaria impedint la unió de l’inhibidor, mentre que el residu Met593 estaria afavorint aquesta unió. Els estudis amb l’enzim CPT1B de porc demostraren l’existència d’un determinant positiu per la sensibilitat al malonil-CoA en els primers 18 residus de la proteïna, i definiren la posició Glu17 com la responsable de l’alta afinitat a la carnitina i la baixa sensibilitat a la inhibició per malonil-CoA (8). Es clonà i caracteritzà la regió promotora del gen de CPT1C humana, amb la intenció d’analitzar la funcionalitat de putatius elements de resposta identificats in silico. Cap dels elements estudiats resultà ser funcional in vivo. A més, es demostrà que la manca d’activitat catalítica de la proteïna no és deguda a l’extensió C-terminal que presenta respecte els isotips A i B, tot i presentar un alt percentatge d’identitat de seqüència. S’ha amplificat una isoforma humana de CPT1C (Pubmed Acc. Num. AK299866), corresponent a la regió carboxiterminal de la proteïna, que es pretén utilitzar per obtenir el primer cristall de la part soluble d’una proteïna CPT1.     

Update of the FANTOM web resource: from mammalian transcriptional landscape to its dynamic regulation.

The international Functional Annotation Of the Mammalian Genomes 4 (FANTOM4) research collaboration set out to better understand the transcriptional network that regulates macrophage differentiation and to uncover novel components of the transcriptome employing a series of high-throughput experiments. The primary and unique technique is cap analysis of gene expression (CAGE), sequencing mRNA 5'-ends with a second-generation sequencer to quantify promoter activities even in the absence of gene annotation. Additional genome-wide experiments complement the setup including short RNA sequencing, microarray gene expression profiling on large-scale perturbation experiments and ChIP-chip for epigenetic marks and transcription factors. All the experiments are performed in a differentiation time course of the THP-1 human leukemic cell line. Furthermore, we performed a large-scale mammalian two-hybrid (M2H) assay between transcription factors and monitored their expression profile across human and mouse tissues with qRT-PCR to address combinatorial effects of regulation by transcription factors. These interdependent data have been analyzed individually and in combination with each other and are published in related but distinct papers. We provide all data together with systematic annotation in an integrated view as resource for the scientific community (http://fantom.gsc.riken.jp/4/). Additionally, we assembled a rich set of derived analysis results including published predicted and validated regulatory interactions. Here we introduce the resource and its update after the initial release.

IL2RA/CD25 gene polymorphisms: uneven association with multiple sclerosis (MS) and type 1 diabetes (T1D).

BACKGROUND: IL-2 receptor (IL2R) alpha is the specific component of the high affinity IL2R system involved in the immune response and in the control of autoimmunity. METHODS AND RESULTS: Here we perform a replication and fine mapping of the IL2RA gene region analyzing 3 SNPs previously associated with multiple sclerosis (MS) and 5 SNPs associated with type 1 diabetes (T1D) in a collection of 798 MS patients and 927 matched Caucasian controls from the south of Spain. We observed association with MS in 6 of 8 SNPs. The rs1570538, at the 3'- UTR extreme of the gene, previously reported to have a weak association with MS, is replicated here (P = 0.032). The most associated T1D SNP (rs41295061) was not associated with MS in the present study. However, the rs35285258, belonging to another independent group of SNPs associated with T1D, showed the maximal association in this study but different risk allele. We replicated the association of only one (rs2104286) of the two IL2RA SNPs identified in the recently performed genome-wide association study of MS. CONCLUSIONS: These findings confirm and extend the association of this gene with MS and reveal a genetic heterogeneity of the associated polymorphisms and risk alleles between MS and T1D suggesting different immunopathological roles of IL2RA in these two diseases.

Epigenetic regulation of human cancer/testis antigen gene, HAGE, in chronic myeloid leukemia.

BACKGROUND AND OBJECTIVES Cancer testis antigens (CTA) provide attractive targets for cancer-specific immunotherapy. Although CTA genes are expressed in some normal tissues, such as the testis, this immunologically protected site lacks MHC I expression and as such, does not present self antigens to T cells. To date, CTA genes have been shown to be expressed in a range of solid tumors via demethylation of their promoter CpG islands, but rarely in chronic myeloid leukemia (CML) or other hematologic malignancies. DESIGN AND METHODS In this study, the methylation status of the HAGE CTA gene promoter was analyzed by quantitative methylation-specific polymerase chain reaction (MSP) and sequencing in four Philadelphia-positive cell lines (TCC-S, K562, KU812 and KYO-1) and in CML samples taken from patients in chronic phase (CP n=215) or blast crisis (BC n=47). HAGE expression was assessed by quantitative reverse transcriptase-polymerase chain reaction. RESULTS The TCC-S cell line showed demethylation of HAGE that was associated with over-expression of this gene. HAGE hypomethylation was significantly more frequent in BC (46%) than in CP (22%) (p=0.01) and was correlated with high expression levels of HAGE transcripts (p<0.0001). Of note, in CP-CML, extensive HAGE hypomethylation was associated with poorer prognosis in terms of cytogenetic response to interferon (p=0.01) or imatinib (p=0.01), molecular response to imatinib (p=0.003) and progression-free survival (p=0.05). INTERPRETATIONS AND CONCLUSION: The methylation status of the HAGE promoter directly correlates with its expression in both CML cell lines and patients and is associated with advanced disease and poor outcome.

The MRC1/CD68 Ratio Is Positively Associated with Adipose Tissue Lipogenesis and with Muscle Mitochondrial Gene Expression in Humans.

BACKGROUND Alternative macrophages (M2) express the cluster differentiation (CD) 206 (MCR1) at high levels. Decreased M2 in adipose tissue is known to be associated with obesity and inflammation-related metabolic disturbances. Here we aimed to investigate MCR1 relative to CD68 (total macrophages) gene expression in association with adipogenic and mitochondrial genes, which were measured in human visceral [VWAT, n = 147] and subcutaneous adipose tissue [SWAT, n = 76] and in rectus abdominis muscle (n = 23). The effects of surgery-induced weight loss were also longitudinally evaluated (n = 6). RESULTS MCR1 and CD68 gene expression levels were similar in VWAT and SWAT. A higher proportion of CD206 relative to total CD68 was present in subjects with less body fat and lower fasting glucose concentrations. The ratio MCR1/CD68was positively associated with IRS1gene expression and with the expression of lipogenic genes such as ACACA, FASN and THRSP, even after adjusting for BMI. The ratio MCR1/CD68 in SWAT increased significantly after the surgery-induced weight loss (+44.7%; p = 0.005) in parallel to the expression of adipogenic genes. In addition, SWAT MCR1/CD68ratio was significantly associated with muscle mitochondrial gene expression (PPARGC1A, TFAM and MT-CO3). AT CD206 was confirmed by immunohistochemistry to be specific of macrophages, especially abundant in crown-like structures. CONCLUSION A decreased ratio MCR1/CD68 is linked to adipose tissue and muscle mitochondrial dysfunction at least at the level of expression of adipogenic and mitochondrial genes.

Regulation of the cardiac voltage-gated sodium channel Nav1.5 : regulation of Nav1.5 by ubiquitin-protein ligases of the Nedd4/Nedd4-like family and characterisation of two novel mutations in the gene encoding Nav1.5 (SCN5A) implicated in the Brugada syndrome triggered by fever

Abstract: The genesis of the cardiac action potential, which accounts for the cardiac contraction, is due to the sodium current INa mediated by the voltage-gated sodium channel Nav1.5. Several cardiac arrhythmias such as the Brugada syndrome are known te be caused by mutations in SCN5A, the gene encoding Nav1.5. Studies of these mutations allowed a better understanding of biophysical and functional properties of Nav1.5. However, only few investigations have been performed in order to understand the regulation of Nav1.5. During my thesis, I investigated different mechanisms of regulation of Nav1.5 using a heterologous expression system, HEK293 cells, coupled with a technique of sodium current recording: the patch clamp in whole cell configuration. In previous studies it has been shown that an enzyme of the Nedd4 family (Nedd4-2) regulates an epithelial sodium channel via the interaction with PY-motifs present in the latter. Interestingly, Nav1.5 contains a similar PY-motif, which motivated us to study the role of Nedd4-2 expressed in heart for the regulation of Nav1.5. In a second study, we investigated the implication of two Nav1.5 mutants, which were either less functional or net functional (Nav1.5 R535X and Nav1.5 L325R respectively) implied in the genesis of the Brugada syndrome by fever. Our results established two mechanisms implied in Nav1.5 regulation. The first one implies that following the interaction between the PY-motif of Nav1.5 and Nedd4- 2 Nav1.5 is ubiquitinated by Nedd4-2. This ubiquitination leads to the internalization of Nav1 .5. The second mechanism is a phenomenon called the "dominant negative" effect of Nav1.5 L325R on Nay1.5 where the decrease of 'Na is potentially due to the retention of Nav1.5 by Nav1.5 L325R in an undefined intracellular compartment. These studies defined two mechanisms of Nav1.5 regulation, which could play an important role for the genesis of cardiac arrhythmias where molecular processes are still poorly understood. Résumé La genèse du potentiel d'action cardiaque, permettant la contraction cardiaque, est due au courant sodique INa issu des canaux sodiques cardiaques dépendants du voltage Nav1.5. Nombreuses arythmies cardiaques telles que le syndrome de Brugada sont connues pour être liées à des mutations du gène SCN5A, codant pour Nav1.5. L'étude de ces mutations a permis une meilleure compréhension des propriétés structurelles et fonctionnelles de Nav1.5 et leurs implications dans la genèse de ces pathologies. Néanmoins peu d'études ont été menées afin de comprendre les mécanismes de régulation de Nav1.5. Mon travail de thèse a consisté à étudier des mécanismes de régulation de Nav1.5 en utilisant un système d'expression hétérologue, les cellules HEK293, couplé à une technique d'enregistrement des courants sodiques, le "patch clamp" en configuration cellule entière. La présence sur Nav1.5 d'un motif-PY similaire à ceux nécessaires pour la régulation d'un canal épithélial sodique par une enzyme de la famille de Nedd4, nous a amenée à étudier le rôle de ces ubiquitine-ligases, en particulier Nedd4-2, dans la régulation de Nav1.5. La seconde étude s'est intéressée aux conséquences de deux mutations de SCN5A codant pour deux mutants peu ou pas fonctionnels (Nav1.5 L325R et Nav1.5 R535X respectivement) retrouvées chez des patients présentant un syndrome de Brugada exacerbé par un état fébrile. Nos résultats ont permis d'établir deux mécanismes de régulation de Nav1.5 L'un par Nedd4-2 qui implique rubiquitination de Nav1.5 par cette ligase suite à l'interaction entre le motif-PY de Nav1.5 et Nedd4-2. Cette modification déclenche l'internalisation du canal impliquée dans la diminution d'INa. Le second mécanisme quant à lui est un effet "dominant négatif" de Nav1.5 L325R sur Nav1.5 aboutissant à une diminution d'INa suite à la séquestration intracellulaire potentielle de Nav1.5 par Nav1.5 L325R. Ces études ont mis en évidence deux mécanismes de régulation de Nav1.5 pouvant jouer un rôle majeur dans la genèse et/ou l'accentuation des arythmies cardiaques dont les processus moléculaires au sein des cardiomyocytes, impliquant des modifications du courant sodiques, sont encore mal compris. Résumé destiné à un large public La dépolarisation électrique de la membrane des cellules cardiaques permet la contraction du coeur. La génèse de cette activité électrique est due au courant sodique issu d'un type de canal à sodium situé dans la membrane des cellules cardiaques. De nombreuses pathologies provoquant des troubles du rythme cardiaque sont issues de mutations du gène qui code pour ce canal à sodium. Ces canaux mutants, entrainant diverses pathologies cardiaques telles que le syndrome de Brugada, ont été largement étudiées. Néanmoins, peu de travaux ont été réalisés sur les mécanismes de régulation de ce canal à sodium non muté. Mon travail de thèse a consisté à étudier certains des mécanismes de régulation de ce canal à sodium en utilisant une technique permettant l'enregistrement des courants sodiques issus de l'expression de ces canaux à sodium à la membrane de cellules mammifères. La présence sur ce canal à sodium d'une structure spécifique, similaire à celle nécessaire pour la régulation d'un canal épithélial à sodium par une enzyme appelée Nedd4-2, nous a amenée à étudier le rôle de cette enzyme dans la régulation de ce canal à sodium. La seconde étude s'est intéressée aux rôles de deux mutations du gène codant pour ce canal à sodium retrouvées chez des patients présentant un syndrome de Brugada exacerbé par la fièvre. Nos résultats nous ont permis d'établir deux mécanismes de régulation de ce canal à sodium diminuant le courant sodique l'un par l'action de l'enzyme Nedd4-2, suite à son interaction avec ce canal, qui modifie ce canal à sodium (ubiquitination) diminuant de ce fait la densité membranaire du canal. L'autre par un mécanisme suggérant un effet négatif de l'un des canaux mutants sur l'expression à la membrane du canal à sodium non muté. Ces études ont mis en évidence deux mécanismes de régulation de ce canal à sodium pouvant jouer un rôle majeur dans la genèse et/ou l'accentuation des troubles du rythme cardiaques dont les mécanismes cellulaires sont encore incompris.

Novel (60%) and recurrent (40%) androgen receptor gene mutations in a series of 59 patients with a 46,XY disorder of sex development.

BACKGROUND Androgen receptor (AR) gene mutations are the most frequent cause of 46,XY disorders of sex development (DSD) and are associated with a variety of phenotypes, ranging from phenotypic women [complete androgen insensitivity syndrome (CAIS)] to milder degrees of undervirilization (partial form or PAIS) or men with only infertility (mild form or MAIS). OBJECTIVE The aim of the study was to characterize the contribution of the AR gene to the molecular cause of 46,XY DSD in a series of Spanish patients. SETTING We studied a series of 133 index patients with 46,XY DSD in whom gonads were differentiated as testes, with phenotypes including varying degrees of undervirilization, and in whom the AR gene was the first candidate for a molecular analysis. METHODS The AR gene was sequenced (exons 1 to 8 with intronic flanking regions) in all patients and in family members of 61% of AR-mutated gene patients. RESULTS AR gene mutations were found in 59 individuals (44.4% of index patients), of whom 46 (78%) were CAIS and 13 (22%) PAIS. Fifty-seven different mutations were found: 21.0% located in exon 1, 15.8% in exons 2 and 3, 57.9% in exons 4-8, and 5.3% intronic. Twenty-three mutations (40.4%) had been previously described and 34 (59.6%) were novel. CONCLUSIONS AR gene mutation is the most frequent cause of 46,XY DSD, with a clearly higher frequency in the complete phenotype. Mutations spread along the whole coding sequence, including exon 1. This series shows that 60% of mutations detected during the period 2002-2009 were novel.

The diamine oxidase gene is associated with hypersensitivity response to non-steroidal anti-inflammatory drugs.

Non-steroidal anti-inflammatory drugs (NSAIDs) are the drugs most frequently involved in hypersensitivity drug reactions. Histamine is released in the allergic response to NSAIDs and is responsible for some of the clinical symptoms. The aim of this study is to analyze clinical association of functional polymorphisms in the genes coding for enzymes involved in histamine homeostasis with hypersensitivity response to NSAIDs. We studied a cohort of 442 unrelated Caucasian patients with hypersensitivity to NSAIDs. Patients who experienced three or more episodes with two or more different NSAIDs were included. If this requirement was not met diagnosis was established by challenge. A total of 414 healthy unrelated controls ethnically matched with patients and from the same geographic area were recruited. Analyses of the SNPs rs17740607, rs2073440, rs1801105, rs2052129, rs10156191, rs1049742 and rs1049793 in the HDC, HNMT and DAO genes were carried out by means of TaqMan assays. The detrimental DAO 16 Met allele (rs10156191), which causes decreased metabolic capacity, is overrepresented among patients with crossed-hypersensitivity to NSAIDs with an OR  = 1.7 (95% CI  = 1.3-2.1; Pc  = 0.0003) with a gene-dose effect (P = 0.0001). The association was replicated in two populations from different geographic areas (Pc  = 0.008 and Pc  = 0.004, respectively). CONCLUSIONS AND IMPLICATIONS: The DAO polymorphism rs10156191 which causes impaired metabolism of circulating histamine is associated with the clinical response in crossed-hypersensitivity to NSAIDs and could be used as a biomarker of response.

A multicenter study confirms CD226 gene association with systemic sclerosis-related pulmonary fibrosis.

INTRODUCTION CD226 genetic variants have been associated with a number of autoimmune diseases and recently with systemic sclerosis (SSc). The aim of this study was to test the influence of CD226 loci in SSc susceptibility, clinical phenotypes and autoantibody status in a large multicenter European population. METHODS A total of seven European populations of Caucasian ancestry were included, comprising 2,131 patients with SSc and 3,966 healthy controls. Three CD226 single nucleotide polymorphisms (SNPs), rs763361, rs3479968 and rs727088, were genotyped using Taqman 5'allelic discrimination assays. RESULTS Pooled analyses showed no evidence of association of the three SNPs, neither with the global disease nor with the analyzed subphenotypes. However, haplotype block analysis revealed a significant association for the TCG haplotype (SNP order: rs763361, rs34794968, rs727088) with lung fibrosis positive patients (PBonf = 3.18E-02 OR 1.27 (1.05 to 1.54)). CONCLUSION Our data suggest that the tested genetic variants do not individually influence SSc susceptibility but a CD226 three-variant haplotype is related with genetic predisposition to SSc-related pulmonary fibrosis.

Gene copy number polymorphisms in an arbuscular mycorrhizal fungal population.

Gene copy number polymorphism was studied in a population of the arbuscular mycorrhizal fungus Glomus intraradices by using a quantitative PCR approach on four different genomic regions. Variation in gene copy number was found for a pseudogene and for three ribosomal genes, providing conclusive evidence for a widespread occurrence of macromutational events in the population.

Response to treatment in Brazilian patients with chronic hepatitis C is associated with a single-nucleotide polymorphism near the interleukin-28B gene

A single-nucleotide polymorphism (SNP) upstream of interleukin (IL)28B was recently identified as an important predictor of the outcome of chronic hepatitis C patients treated with pegylated interferon plus ribavirin (PEG-IFN/RBV). The aim of this study was to investigate the association between the IL28B gene polymorphism (rs12979860) and virological response in chronic hepatitis C patients. Brazilian patients (n = 263) who were infected with hepatitis C virus (HCV) genotype 1 and were receiving PEG-IFN/RBV were genotyped. Early virological response (EVR) (12 weeks), end-of-treatment response (EOTR) (48 weeks), sustained virological response (SVR) (72 weeks) and relapse were evaluated using conventional and quantitative polymerase chain reaction (PCR) assays. The frequency of the C allele in the population was 39%. Overall, 43% of patients experienced SVR. The IL28B CC genotype was significantly associated with higher treatment response rates and a lower relapse rate compared to the other genotypes [84% vs. 58% EVR, 92% vs. 63% EOTR, 76% vs. 38% SVR and 17% vs. 40% relapse rate in CC vs. other genotypes (CT and TT), respectively]. Thus, the IL28B genotype appears to be a strong predictor of SVR following PEG-IFN/RBV therapy in treatment-naïve Brazilian patients infected with HCV genotype 1. This study, together with similar research examining other SNPs, should help to define adequate protocols for the treatment of patients infected with HCV genotype 1, especially those with a poor prognosis.

Equine herpesvirus-2 E10 gene product, but not its cellular homologue, activates NF-kappaB transcription factor and c-Jun N-terminal kinase.

We have previously reported on the death effector domain containing E8 gene product from equine herpesvirus-2, designated FLICE inhibitory protein (v-FLIP), and on its cellular homologue, c-FLIP, which inhibit the activation of caspase-8 by death receptors. Here we report on the structure and function of the E10 gene product of equine herpesvirus-2, designated v-CARMEN, and on its cellular homologue, c-CARMEN, which contain a caspase-recruiting domain (CARD) motif. c-CARMEN is highly homologous to the viral protein in its N-terminal CARD motif but differs in its C-terminal extension. v-CARMEN and c-CARMEN interact directly in a CARD-dependent manner yet reveal different binding specificities toward members of the tumor necrosis factor receptor-associated factor (TRAF) family. v-CARMEN binds to TRAF6 and weakly to TRAF3 and, upon overexpression, potently induces the c-Jun N-terminal kinase (JNK), p38, and nuclear factor (NF)-kappaB transcriptional pathways. c-CARMEN or truncated versions thereof do not appear to induce JNK and NF-kappaB activation by themselves, nor do they affect the JNK and NF-kappaB activating potential of v-CARMEN. Thus, in contrast to the cellular homologue, v-CARMEN may have additional properties in its unique C terminus that allow for an autonomous activator effect on NF-kappaB and JNK. Through activation of NF-kappaB, v-CARMEN may regulate the expression of the cellular and viral genes important for viral replication.

Extended genetic analysis of Brazilian isolates of Bacillus cereus and Bacillus thuringiensis

Multiple locus sequence typing (MLST) was undertaken to extend the genetic characterization of 29 isolates of Bacillus cereus and Bacillus thuringiensis previously characterized in terms of presence/absence of sequences encoding virulence factors and via variable number tandem repeat (VNTR). Additional analysis involved polymerase chain reaction for the presence of sequences (be, cytK, inA, pag, lef, cya and cap), encoding putative virulence factors, not investigated in the earlier study. MLST analysis ascribed novel and unique sequence types to each of the isolates. A phylogenetic tree was constructed from a single sequence of 2,838 bp of concatenated loci sequences. The strains were not monophyletic by analysis of any specific housekeeping gene or virulence characteristic. No clear association in relation to source of isolation or to genotypic profile based on the presence or absence of putative virulence genes could be identified. Comparison of VNTR profiling with MLST data suggested a correlation between these two methods of genetic analysis. In common with the majority of previous studies, MLST was unable to provide clarification of the basis for pathogenicity among members of the B. cereus complex. Nevertheless, our application of MLST served to reinforce the notion that B. cereus and B. thuringiensis should be considered as the same species.

Comparative analyses of classical phenotypic method and ribosomal RNA gene sequencing for identification of medically relevant Candida species

As the distribution of Candida species and their susceptibility to antifungal agents have changed, a new means of accurately and rapidly identifying these species is necessary for the successful early resolution of infection and the subsequent reduction of morbidity and mortality. The current work aimed to evaluate ribosomal RNA gene sequencing for the identification of medically relevant Candida species in comparison with a standard phenotypic method. Eighteen reference strains (RSs), 69 phenotypically identified isolates and 20 inconclusively identified isolates were examined. Internal transcribed spaces (ITSs) and D1/D2 of the 26S ribosomal RNA gene regions were used as targets for sequencing. Additionally, the sequences of the ITS regions were used to establish evolutionary relationships. The sequencing of the ITS regions was successful for 88% (94/107) of the RS and isolates, whereas 100% of the remaining 12% (13/107) of the samples were successfully analysed by sequencing the D1/D2 region. Similarly, genotypic analysis identified all of the RS and isolates, including the 20 isolates that were not phenotypically identified. Phenotypic analysis, however, misidentified 10% (7/69) of the isolates. Phylogenetic analysis allowed the confirmation of the relationships between evolutionarily close species. Currently, the use of genotypic methods is necessary for the correct identification of Candida species.

Fine mapping and functional analysis of the multiple sclerosis risk gene CD6.

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2(nd) SRCR domain with susceptibility to MS (P max(T) permutation = 1×10(-4)). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. - CD4(+) naïve cells, P = 0.0001; CD8(+) naïve cells, P<0.0001; CD4(+) and CD8(+) central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4(+) and CD8(+) T cells.