Magdalena Zoledziewska

Variants within the immunoregulatory CBLB gene are associated with Multiple Sclerosis

A genome-wide association scan of ∼6.6 million genotyped or imputed variants in 882 Sardinian individuals with multiple sclerosis (cases) and 872 controls suggested association of CBLB gene variants with disease, which was confirmed in 1,775 cases and 2,005 controls (rs9657904, overall P = 1.60 × 10−10, OR = 1.40). CBLB encodes a negative regulator of adaptive immune responses, and mice lacking the ortholog are prone to experimental autoimmune encephalomyelitis, the animal model of multiple sclerosis.

Low-Pass DNA Sequencing of 1200 Sardinians Reconstructs European Y-Chromosome Phylogeny

Examining Y The evolution of human populations has long been studied with unique sequences from the nonrecombining, male-specific Y chromosome (see the Perspective by Cann). Poznik et al. (p. 562) examined 9.9 Mb of the Y chromosome from 69 men from nine globally divergent populations—identifying population and individual specific sequence variants that elucidate the evolution of the Y chromosome. Sequencing of maternally inherited mitochondrial DNA allowed comparison between the relative rates of evolution, which suggested that the coalescence, or origin, of the human Y chromosome and mitochondria both occurred approximately 120 thousand years ago. Francalacci et al. (p. 565) investigated the sequence divergence of 1204 Y chromosomes that were sampled within the isolated and genetically informative Sardinian population. The sequence analyses, along with archaeological records, were used to calibrate and increase the resolution of the human phylogenetic tree. Local human demographic history is inferred from in-depth DNA sequence analysis of Sardinian mens Y chromosomes. [Also see Perspective by Cann] Genetic variation within the male-specific portion of the Y chromosome (MSY) can clarify the origins of contemporary populations, but previous studies were hampered by partial genetic information. Population sequencing of 1204 Sardinian males identified 11,763 MSY single-nucleotide polymorphisms, 6751 of which have not previously been observed. We constructed a MSY phylogenetic tree containing all main haplogroups found in Europe, along with many Sardinian-specific lineage clusters within each haplogroup. The tree was calibrated with archaeological data from the initial expansion of the Sardinian population ~7700 years ago. The ages of nodes highlight different genetic strata in Sardinia and reveal the presumptive timing of coalescence with other human populations. We calculate a putative age for coalescence of ~180,000 to 200,000 years ago, which is consistent with previous mitochondrial DNA–based estimates.

A genome-wide association scan on the levels of markers of inflammation in sardinians reveals associations that underpin its complex regulation

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals—5 of which were identified only with the custom arrays—and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10−29); for ESR, at the HBB (rs4910472, p = 2.31×10−11) and UCN119B/SPPL3 (rs11829037, p = 8.91×10−10) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10−13) and in CADM3 (rs3026968, p = 7.63×10−13); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10−21), near DARC (rs3845624, p = 1.43×10−10), UNC119B/SPPL3 (rs11829037, p = 1.50×10−14), and ICOSLG/AIRE (rs113459440, p = 1.54×10−08) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.

Variation within the CLEC16A gene shows consistent disease association with both multiple sclerosis and type 1 diabetes in Sardinia

Variation within intron 19 of the CLEC16A (KIAA0350) gene region was recently found to be unequivocally associated with type 1 diabetes (T1D) in genome-wide association (GWA) studies in Northern European populations. A variant in intron 22 that is nearly independent of the intron 19 variant showed suggestive evidence of association with multiple sclerosis (MS). Here, we genotyped the rs725613 polymorphism, representative of the earlier reported associations with T1D within CLEC16A, in 1037 T1D cases, 1498 MS cases and 1706 matched controls, all from the founder, autoimmunity-prone Sardinian population. In these Sardinian samples, allele A of rs725613 is positively associated not only with T1D (odds ratio=1.15, P one-tail=5.1 × 10−3) but also, and with a comparable effect size, with MS (odds ratio=1.21, P one-tail 6.7 × 10−5). Taken together these data provide evidence of joint disease association in T1D and MS within CLEC16A and underline a shared disease pathway.

Genome sequencing elucidates Sardinian genetic architecture and augments association analyses for lipid and blood inflammatory markers

We report ∼17.6 million genetic variants from whole-genome sequencing of 2,120 Sardinians; 22% are absent from previous sequencing-based compilations and are enriched for predicted functional consequences. Furthermore, ∼76,000 variants common in our sample (frequency >5%) are rare elsewhere (<0.5% in the 1000 Genomes Project). We assessed the impact of these variants on circulating lipid levels and five inflammatory biomarkers. We observe 14 signals, including 2 major new loci, for lipid levels and 19 signals, including 2 new loci, for inflammatory markers. The new associations would have been missed in analyses based on 1000 Genomes Project data, underlining the advantages of large-scale sequencing in this founder population.

Overexpression of the Cytokine BAFF and Autoimmunity Risk

BACKGROUND Genomewide association studies of autoimmune diseases have mapped hundreds of susceptibility regions in the genome. However, only for a few association signals has the causal gene been identified, and for even fewer have the causal variant and underlying mechanism been defined. Coincident associations of DNA variants affecting both the risk of autoimmune disease and quantitative immune variables provide an informative route to explore disease mechanisms and drug‐targetable pathways. METHODS Using case–control samples from Sardinia, Italy, we performed a genomewide association study in multiple sclerosis followed by TNFSF13B locus–specific association testing in systemic lupus erythematosus (SLE). Extensive phenotyping of quantitative immune variables, sequence‐based fine mapping, cross‐population and cross‐phenotype analyses, and gene‐expression studies were used to identify the causal variant and elucidate its mechanism of action. Signatures of positive selection were also investigated. RESULTS A variant in TNFSF13B, encoding the cytokine and drug target B‐cell activating factor (BAFF), was associated with multiple sclerosis as well as SLE. The disease‐risk allele was also associated with up‐regulated humoral immunity through increased levels of soluble BAFF, B lymphocytes, and immunoglobulins. The causal variant was identified: an insertion–deletion variant, GCTGT→A (in which A is the risk allele), yielded a shorter transcript that escaped microRNA inhibition and increased production of soluble BAFF, which in turn up‐regulated humoral immunity. Population genetic signatures indicated that this autoimmunity variant has been evolutionarily advantageous, most likely by augmenting resistance to malaria. CONCLUSIONS A TNFSF13B variant was associated with multiple sclerosis and SLE, and its effects were clarified at the population, cellular, and molecular levels. (Funded by the Italian Foundation for Multiple Sclerosis and others.)

Further evidence of a primary, causal association of the PTPN22 620W variant with type 1 diabetes

OBJECTIVE— The minor allele of the nonsynonymous single nucleotide polymorphism (SNP) +1858C>T within the PTPN22 gene is positively associated with type 1 diabetes and other autoimmune diseases. Genetic and functional data underline its causal effect, but some studies suggest that this polymorphism does not entirely explain disease association of the PTPN22 region. The aim of this study was to evaluate type 1 diabetes association within this gene in the Sardinian population. RESEARCH DESIGN AND METHODS— We resequenced the exons and potentially relevant portions of PTPN22 and detected 24 polymorphisms (23 SNPs and 1 deletion insertion polymorphism [DIP]), 8 of which were novel. A representative set of 14 SNPs and the DIP were sequentially genotyped and assessed for disease association in 794 families, 490 sporadic patients, and 721 matched control subjects. RESULTS— The +1858C>T variant, albeit rare in the general Sardinian population (allele frequency 0.014), was positively associated with type 1 diabetes (Pone tail = 3.7 × 10−3). In contrast, the background haplotype in which this mutation occurred was common (haplotype frequency 0.117) and neutrally associated with disease. We did not confirm disease associations reported in other populations for non +1858C>T variants (rs2488457, rs1310182, and rs3811021), although they were present in appreciable frequencies in Sardinia. Additional weak disease associations with rare variants were detected in the Sardinian families but not confirmed in independent case-control sample sets and are most likely spurious. CONCLUSIONS— We provide further evidence that the +1858C>T polymorphism is primarily associated with type 1 diabetes and exclude major contributions from other purportedly relevant variants within this gene.

Rare and functional SIAE variants are not associated with autoimmune disease risk in up to 66,924 individuals of European ancestry

Rare and functional SIAE variants are not associated with autoimmune disease risk in up to 66,924 individuals of European ancestry

Rare variant genotype imputation with thousands of study-specific whole-genome sequences: implications for cost-effective study designs

The utility of genotype imputation in genome-wide association studies is increasing as progressively larger reference panels are improved and expanded through whole-genome sequencing. Developing general guidelines for optimally cost-effective imputation, however, requires evaluation of performance issues that include the relative utility of study-specific compared with general/multipopulation reference panels; genotyping with various array scaffolds; effects of different ethnic backgrounds; and assessment of ranges of allele frequencies. Here we compared the effectiveness of study-specific reference panels to the commonly used 1000 Genomes Project (1000G) reference panels in the isolated Sardinian population and in cohorts of European ancestry including samples from Minnesota (USA). We also examined different combinations of genome-wide and custom arrays for baseline genotypes. In Sardinians, the study-specific reference panel provided better coverage and genotype imputation accuracy than the 1000G panels and other large European panels. In fact, even gene-centered custom arrays (interrogating ~200 000 variants) provided highly informative content across the entire genome. Gain in accuracy was also observed for Minnesotans using the study-specific reference panel, although the increase was smaller than in Sardinians, especially for rare variants. Notably, a combined panel including both study-specific and 1000G reference panels improved imputation accuracy only in the Minnesota sample, and only at rare sites. Finally, we found that when imputation is performed with a study-specific reference panel, cutoffs different from the standard thresholds of MACH-Rsq and IMPUTE-INFO metrics should be used to efficiently filter badly imputed rare variants. This study thus provides general guidelines for researchers planning large-scale genetic studies.

Genome-wide association study of multiple sclerosis confirms a novel locus at 5p13.1.

Multiple Sclerosis (MS) is the most common progressive and disabling neurological condition affecting young adults in the world today. From a genetic point of view, MS is a complex disorder resulting from the combination of genetic and non-genetic factors. We aimed to identify previously unidentified loci conducting a new GWAS of Multiple Sclerosis (MS) in a sample of 296 MS cases and 801 controls from the Spanish population. Meta-analysis of our data in combination with previous GWAS was done. A total of 17 GWAS-significant SNPs, corresponding to three different loci were identified:HLA, IL2RA, and 5p13.1. All three have been previously reported as GWAS-significant. We confirmed our observation in 5p13.1 for rs9292777 using two additional independent Spanish samples to make a total of 4912 MS cases and 7498 controls (ORpooled = 0.84; 95%CI: 0.80–0.89; p = 1.36×10-9). This SNP differs from the one reported within this locus in a recent GWAS. Although it is unclear whether both signals are tapping the same genetic association, it seems clear that this locus plays an important role in the pathogenesis of MS.