Identification by denaturing high-performance liquid chromatography of numerous polymorphisms in a candidate region for multiple sclerosis susceptibility

Genetic association analysis of candidate regions where evidence of linkage has accumulated is becoming a key issue in the study of complex diseases. A high density of markers, at least one per centimorgan, is required to improve the chances of observing linkage disequilibrium with disease alleles. A recently available single nucleotide polymorphism (SNP) map designed to cover the whole genome provides an average density of one marker per 2 cM. In the present study we show that the number of markers can be approximately doubled in a selected region, thus reaching a density suitable for association studies, by applying a completely automated technique for polymorphism detection, denaturing high-performance liquid chromatography (DHPLC). A systematic search for SNPs was performed in the region 5ptel-q13, where weak but convergent evidence for linkage with multiple sclerosis has accumulated. Screening for polymorphisms was performed on 124 sequence tagged sites (STSs) in the 3'UTR ends of expressed sequence tags totaling about 30,000 bp. Thirty SNPs in 28 STSs were found with less than 10% overlap with the markers already detected in the same region. The data confirm the validity of the approach using DHPLC on expressed gene sequences tagged by a set of standard commercially available primers.     

A functional variant in ERAP1 predisposes to multiple sclerosis

The ERAP1 gene encodes an aminopeptidase involved in antigen processing. A functional polymorphism in the gene (rs30187, Arg528Lys) associates with susceptibility to ankylosying spondylitis (AS), whereas a SNP in the interacting ERAP2 gene increases susceptibility to another inflammatory autoimmune disorder, Crohn''s disease (CD). We analysed rs30187 in 572 Italian patients with CD and in 517 subjects suffering from multiple sclerosis (MS); for each cohort, an independent sex- and age-matched control group was genotyped. The frequency of the 528Arg allele was significantly higher in both disease cohorts compared to the respective control population (for CD, OR = 1.20 95%CI: 1.01–1.43, p = 0.036; for RRMS, OR = 1.26; 95%CI: 1.04–1.51, p = 0.01). Meta-analysis with the Wellcome Trust Cases Control Consortium GWAS data confirmed the association with MS (p_meta = 0.005), but not with CD. In AS, the rs30187 variant has a predisposing effect only in an HLA-B27 allelic background. It remains to be evaluated whether interaction between ERAP1 and distinct HLA class I alleles also affects the predisposition to MS, and explains the failure to provide definitive evidence for a role of rs30187 in CD. Results herein support the emerging concept that a subset of master-regulatory genes underlay the pathogenesis of autoimmunity.     

Marburg variant multiple sclerosis - a case report

In this case report we describe the case of a 24 year-old female with a fulminant demyelinating disease of white matter. Disease progression was most probably consistent with the Marburg variant (malignant form) of multiple sclerosis with rapid deterioration of the patient's clinical condition, including bulbar symptoms and epileptic paroxysms and ending with persistent coma and tetraparesis, over the course of 6 months from first symptoms. Repeated Magnetic Resonance Imaging (MRI) examination showed progression of multiple demyelinating lesions culminating in a contiguous focal disorder of the white matter extending both supratentorially and infratentorially. The serial MRI changes closely mapped the deterioration in the patients clinical status. Our patient showed no response to repeated pulse corticotherapy, administration of intravenous immunoglobulins, serial plasmapheresis, and combined high-dose pulse immunosuppression (specify what was used here) and mitoxantrone.     

A search for susceptibility genes in multiple sclerosis

Multiple sclerosis is a chronic inflammatory, putatively autoimmune disease characterized by multifocal demyelination in the central nervous system. Two main strategies are used to identify genes influencing the susceptibility to multiple sclerosis: (i) elucidation of the role of a candidate gene chosen on the basis of the possible function of the encoded protein in etiology and/or pathogenesis of the disease, and (ii) complete genomic screen using a panel of anonymous genetic markers for identification of the chromosome regions involved in the disease development. The complete genomic search revealed multiple loci for multiple sclerosis on thirteen chromosomes, and analysis of the candidate genes added three more chromosomes to this list. The combined data prove the polygenic nature of this complex disease. Detection of individual genes responsible for susceptibility to multiple sclerosis is complicated by the genetic heterogeneity of analyzed populations and families, which is determined both by the ethnic heterogeneity and the peculiarity of clinical forms of the disease. However, it seems highly probable that HLA and non-HLA genes of the major histocompatibility complex, as well as some unidentified genes on chromosomes 5p and 17q, are involved in the disease development. In addition to HLA, some authors have also shown that a role in the disease development is played by the genes encoding other components of the trimolecular complex involved in antigen presentation: those of the T-cell receptor and the best studied autoantigen, the myelin basic protein. The most promising for further studies of the genetic susceptibility to multiple sclerosis are approaches that combine the candidate-gene strategy with the complete genomic search as well as distinguish the genetically differentiated forms of the disease.     

TCRBV20S1 polymorphism does not influence the susceptibility to type 1 diabetes and multiple sclerosis in the Sardinian population

Among the different T-cell receptor (TCR) BV20S1 polymorphisms, nucleotide substitution at position 524 results in the introduction of a stop codon, whose potential functional relevance is still unknown. We have recently showed in Sardinian subjects the most elevated allele frequency ever reported worldwide for this “null allele” (0.44). As this variant generates a gap in the TCR repertoire, this preliminary finding prompted us to further analyze the role of this polymorphism in the susceptibility to type 1 diabetes (T1D) and multiple sclerosis (MS), which are extremely common in this population. With this aim, we evaluated the influence of the TCRBV20S1 polymorphism by assessing it with the transmission disequilibirum test (TDT) in 652 T1D and 616 MS families, without detecting any significant difference. We conclude that the high frequency of this null allele in Sardinia is not directly related to the high incidence of these autoimmune diseases observed in this founder population.     

The genetic basis of multiple sclerosis: a model for MS susceptibility

Background  

MS-pathogenesis is known to involve both multiple environmental events, and several independent genetic risk-factors.     

HLA-DRB1-DQB1 haplotypes confer susceptibility and resistance to multiple sclerosis in Sardinia

Introduction

Genetic predisposition to multiple sclerosis (MS) in Sardinia (Italy) has been associated with five DRB1*-DQB1* haplotypes of the human leukocyte antigen (HLA). Given the complexity of these associations, an in-depth re-analysis was performed with the specific aims of confirming the haplotype associations; establishing the independence of the associated haplotypes; and assessing patients'' genotypic risk of developing MS.

Methods and Results

A transmission disequilibrium test (TDT) of the DRB1*-DQB1* haplotypes in 943 trio families, confirmed a higher than expected transmission rate (over-transmission) of the *13:03-*03:01 (OR = 2.9, P = 7.6×10⁻³), *04:05-*03:01 (OR = 2.4, P = 4.4×10⁻⁶) and *03:01-*02:01 (OR = 2.1, P = 1.0×10⁻¹⁵) haplotype. In contrast, the *16:01-*05:02 (OR = 0.5, P = 5.4×10⁻¹¹) and the *15:02-*06:01 (OR = 0.3, P = 1.5×10⁻³) haplotypes exhibited a lower than expected transmission rate (under-transmission). The independence of the transmission of each positively and negatively associated haplotype was confirmed relative to all positively associated haplotypes, and to the negatively associated *16:01-*05:02 haplotype. In patients, carriage of two predisposing haplotypes, or of protective haplotypes, respectively increased or decreased the patient''s risk of developing MS. The risk of MS followed a multiplicative model of genotypes, which was, in order of decreasing ORs: *04:05-*0301/*03:01-*02:01 (OR = 4.5); *03:01-*02:01/*03:01-*02:01 (OR = 4.1); and the *16:01-*05:02/*16:01-*0502 (OR = 0.2) genotypes. Analysis of DRB1 and DQB1 protein chain residues showed that the Val/Gly residue at position 86 of the DRB1 chain was the only difference between the protective *16:01- *15:02 alleles and the predisposing *15:01 one. Similarly, the Ala/Val residue at position 38 of the DQB1 chain differentiated the positively associated *06:02 allele and the negatively associated *05:02, *06:01 alleles.

Conclusions

These findings show that the association of specific, independent DRB1*-DQB1* haplotypes confers susceptibility or resistance to MS in the MS-prone Sardinian population. The data also supports a functional role for specific residues of the DRB1 and DQB1 proteins in predisposing patients to MS.     

Identification of novel susceptibility Loci for kawasaki disease in a Han chinese population by a genome-wide association study

Kawasaki disease (KD) is an acute systemic vasculitis syndrome that primarily affects infants and young children. Its etiology is unknown; however, epidemiological findings suggest that genetic predisposition underlies disease susceptibility. Taiwan has the third-highest incidence of KD in the world, after Japan and Korea. To investigate novel mechanisms that might predispose individuals to KD, we conducted a genome-wide association study (GWAS) in 250 KD patients and 446 controls in a Han Chinese population residing in Taiwan, and further validated our findings in an independent Han Chinese cohort of 208 cases and 366 controls. The most strongly associated single-nucleotide polymorphisms (SNPs) detected in the joint analysis corresponded to three novel loci. Among these KD-associated SNPs three were close to the COPB2 (coatomer protein complex beta-2 subunit) gene: rs1873668 (p = 9.52×10⁻⁵), rs4243399 (p = 9.93×10⁻⁵), and rs16849083 (p = 9.93×10⁻⁵). We also identified a SNP in the intronic region of the ERAP1 (endoplasmic reticulum amino peptidase 1) gene (rs149481, p_best = 4.61×10⁻⁵). Six SNPs (rs17113284, rs8005468, rs10129255, rs2007467, rs10150241, and rs12590667) clustered in an area containing immunoglobulin heavy chain variable regions genes, with p_best-values between 2.08×10⁻⁵ and 8.93×10⁻⁶, were also identified. This is the first KD GWAS performed in a Han Chinese population. The novel KD candidates we identified have been implicated in T cell receptor signaling, regulation of proinflammatory cytokines, as well as antibody-mediated immune responses. These findings may lead to a better understanding of the underlying molecular pathogenesis of KD.     

Genome-wide interaction-based association analysis identified multiple new susceptibility Loci for common diseases

Genome-wide interaction-based association (GWIBA) analysis has the potential to identify novel susceptibility loci. These interaction effects could be missed with the prevailing approaches in genome-wide association studies (GWAS). However, no convincing loci have been discovered exclusively from GWIBA methods, and the intensive computation involved is a major barrier for application. Here, we developed a fast, multi-thread/parallel program named “pair-wise interaction-based association mapping” (PIAM) for exhaustive two-locus searches. With this program, we performed a complete GWIBA analysis on seven diseases with stringent control for false positives, and we validated the results for three of these diseases. We identified one pair-wise interaction between a previously identified locus, C1orf106, and one new locus, TEC, that was specific for Crohn''s disease, with a Bonferroni corrected P<0.05 (P = 0.039). This interaction was replicated with a pair of proxy linked loci (P = 0.013) on an independent dataset. Five other interactions had corrected P<0.5. We identified the allelic effect of a locus close to SLC7A13 for coronary artery disease. This was replicated with a linked locus on an independent dataset (P = 1.09×10⁻⁷). Through a local validation analysis that evaluated association signals, rather than locus-based associations, we found that several other regions showed association/interaction signals with nominal P<0.05. In conclusion, this study demonstrated that the GWIBA approach was successful for identifying novel loci, and the results provide new insights into the genetic architecture of common diseases. In addition, our PIAM program was capable of handling very large GWAS datasets that are likely to be produced in the future.     

Computer programs for population genetics data analysis: a survival guide

The analysis of genetic diversity within species is vital for understanding evolutionary processes at the population level and at the genomic level. A large quantity of data can now be produced at an unprecedented rate, requiring the use of dedicated computer programs to extract all embedded information. Several statistical packages have been recently developed, which offer a panel of standard and more sophisticated analyses. We describe here the functionalities, special features and assumptions of more than 20 such programs, indicate how they can interoperate, and discuss new directions that could lead to improved software and analyses.     

A polymorphism in the HLA-DPB1 gene is associated with susceptibility to multiple sclerosis

We conducted an association study across the human leukocyte antigen (HLA) complex to identify loci associated with multiple sclerosis (MS). Comparing 1927 SNPs in 1618 MS cases and 3413 controls of European ancestry, we identified seven SNPs that were independently associated with MS conditional on the others (each P ≤ 4 x 10(-6)). All associations were significant in an independent replication cohort of 2212 cases and 2251 controls (P ≤ 0.001) and were highly significant in the combined dataset (P ≤ 6 x 10(-8)). The associated SNPs included proxies for HLA-DRB1*15:01 and HLA-DRB1*03:01, and SNPs in moderate linkage disequilibrium (LD) with HLA-A*02:01, HLA-DRB1*04:01 and HLA-DRB1*13:03. We also found a strong association with rs9277535 in the class II gene HLA-DPB1 (discovery set P = 9 x 10(-9), replication set P = 7 x 10(-4), combined P = 2 x 10(-10)). HLA-DPB1 is located centromeric of the more commonly typed class II genes HLA-DRB1, -DQA1 and -DQB1. It is separated from these genes by a recombination hotspot, and the association is not affected by conditioning on genotypes at DRB1, DQA1 and DQB1. Hence rs9277535 represents an independent MS-susceptibility locus of genome-wide significance. It is correlated with the HLA-DPB1*03:01 allele, which has been implicated previously in MS in smaller studies. Further genotyping in large datasets is required to confirm and resolve this association.     

Identification of coronin‐1a as a novel antibody target for clinically isolated syndrome and multiple sclerosis

Recently, we identified the mimotope UH‐CIS6 as a novel candidate antibody target for clinically isolated syndrome (CIS) and relapsing‐remitting (RR) multiple sclerosis (MS). The purpose of this study was to further validate UH‐CIS6 as an antibody target for CIS and MS and to identify the in vivo antibody target of UH‐CIS6. First, a UH‐CIS6 peptide ELISA was optimized. Next, we investigated the antibody response toward UH‐CIS6 in cerebrospinal fluid (CSF) from patients with CIS (n = 20), MS (n = 43) and other neurological diseases (n = 42). Immunoprecipitation of anti‐UH‐CIS6 antibodies on a normal human brain lysate was performed to identify the in vivo antibody target of UH‐CIS6. The cellular expression of an in vivo candidate target was investigated by immunohistochemistry using MS brain tissue sections. Antibody reactivity toward UH‐CIS6 was detected in a significantly increased proportion of CSF samples from CIS and RR‐MS patients as compared with neurological controls (p = 0.046). We identified and confirmed coronin‐1a as the in vivo antibody target for UH‐CIS6. Furthermore, coronin‐1a was expressed by T cells and macrophages in an active MS lesion. Together, these results demonstrate that coronin‐1a is a novel antibody target for CIS and MS.     

A new era of human population genetics

A growing resource of methicillin-resistant Staphylococcus aureus (MRSA) genomes uncovers intriguing phylogeographic and recombination patterns and highlights challenges in identifying the source of these phenomena.     

Measles antigen in multiple sclerosis: Identification in the jejunum by immunofluorescence

Measles virus antigen has been detected by the fluorescent antibody method in jejunal mucosa obtained from 24 patients with clinical multiple sclerosis. Deposits of antigen, various components of the complement system and on occasion immunoglobulins were identified in the epithelial basement membrane. Antigen and complement were usually found in the lamina propria, and sometimes could be seen within epithelial cells and in capillary walls. These findings support the concept that multiple sclerosis is caused by persistent measles infection, and indicate that the virus is harbored in the wall of the small intestine.     

Identification of spatial genetic boundaries using a multifractal model in human population genetics

There are two purposes in displaying spatial genetic structure. One is that a visual representation of the variation of the genetic variable should be provided in the contour map. The other is that spatial genetic structure should be reflected by the patterns or the gradients with genetic boundaries in the map. Nevertheless, most conventional interpolation methods, such as Cavalli-Sforza's method in genography, inverse distance-weighted methods, and the Kriging technique, focus only on the first primary purpose because of their arbitrary thresholds marked on the maps. In this paper we present an application of the contour area multifractal model (CAMM) to human population genetics. The method enables the analysis of the geographic distribution of a genetic marker and provides an insight into the spatial and geometric properties of obtained patterns. Furthermore, the CAMM may overcome some of the limitations of other interpolation techniques because no arbitrary thresholds are necessary in the computation of genetic boundaries. The CAMM is built by establishing power law relationships between the area A (> or =rho) in the contour map and the value p itself after plotting these values on a log-log graph. A series of straight-line segments can be fitted to the points on the log-log graph, each representing a power law relationship between the area A (> or =rho) and the cutoff genetic variable value for rho in a particular range. These straight-line segments can yield a group of cutoff values, which can be identified as the genetic boundaries that can classify the map of genetic variable into discrete genetic zones. These genetic zones usually correspond to spatial genetic structure on the landscape. To provide a better understanding of the interest in the CAMM approach, we analyze the spatial genetic structures of three loci (ABO, HLA-A, and TPOX) in China using the CAMM. Each synthetic principal component (SPC) contour map of the three loci is created by using both Han and minority groups data together. These contour maps all present an obvious geographic diversity, which gradually increases from north to south, and show that the genetic differences among populations in different districts of the same nationality are greater than those among different nationalities of the same district. It is surprising to find that both the value of p and the fractal dimension alpha have a clear north to south gradient for each locus, and the same clear boundary between southern and northern Asians in each contour map is still seen in the zone of the Yangtze River, although substantial population migrations have occurred because of war or famine in the last 2,000 or 3,000 years. A clear genetic boundary between Europeans and Asians in each contour map is still seen in northwestern China with a small value of alpha, although the genetic gradient caused by gene flow between Europeans and Asians has tended to show expansion from northwestern China. From the three contour maps another interesting result can be found: The values of alpha north of the Yangtze River are generally less than those south of the Yangtze River. This indicates that the genetic differences among the populations north of the Yangtze River are generally smaller than those in populations south of the Yangtze River.     

Identification of a novel prostate cancer susceptibility variant in the KLK3 gene transcript

Genome-wide association studies (GWAS) have identified more than 30 prostate cancer (PrCa) susceptibility loci. One of these (rs2735839) is located close to a plausible candidate susceptibility gene, KLK3, which encodes prostate-specific antigen (PSA). PSA is widely used as a biomarker for PrCa detection and disease monitoring. To refine the association between PrCa and variants in this region, we used genotyping data from a two-stage GWAS using samples from the UK and Australia, and the Cancer Genetic Markers of Susceptibility (CGEMS) study. Genotypes were imputed for 197 and 312 single nucleotide polymorphisms (SNPs) from HapMap2 and the 1000 Genome Project, respectively. The most significant association with PrCa was with a previously unidentified SNP, rs17632542 (combined P?=?3.9?×?10(-22)). This association was confirmed by direct genotyping in three stages of the UK/Australian GWAS, involving 10,405 cases and 10,681 controls (combined P?=?1.9?×?10(-34)). rs17632542 is also shown to be associated with PSA levels and it is a non-synonymous coding SNP (Ile179Thr) in KLK3. Using molecular dynamic simulation, we showed evidence that this variant has the potential to introduce alterations in the protein or affect RNA splicing. We propose that rs17632542 may directly influence PrCa risk.     

Investigation of the role of mitochondrial DNA in multiple sclerosis susceptibility

Several lines of evidence suggest that mitochondrial genetic factors may influence susceptibility to multiple sclerosis. To explore this hypothesis further, we re-sequenced the mitochondrial genome (mtDNA) from 159 patients with multiple sclerosis and completed a haplogroup analysis including a further 835 patients and 1,506 controls. A trend towards over-representation of super-haplogroup U was the only evidence for association with mtDNA that we identified in these samples. In a parallel analysis of nuclear encoded mitochondrial genes, we also found a trend towards association with the complex I gene, NDUFS2. These results add to the evidence suggesting that variation in mtDNA and nuclear encoded mitochondrial genes may contribute to disease susceptibility in multiple sclerosis.