High-density SNP mapping of the HLA region identifies multiple independent susceptibility loci associated with selective IgA deficiency

Selective IgA deficiency (IgAD; serum IgA<0.07 g/l) is the most common form of human primary immune deficiency, affecting approximately 1∶600 individuals in populations of Northern European ancestry. The polygenic nature of IgAD is underscored by the recent identification of several new risk genes in a genome-wide association study. Among the characterized susceptibility loci, the association with specific HLA haplotypes represents the major genetic risk factor for IgAD. Despite the robust association, the nature and location of the causal variants in the HLA region remains unknown. To better characterize the association signal in this region, we performed a high-density SNP mapping of the HLA locus and imputed the genotypes of common HLA-B, -DRB1, and -DQB1 alleles in a combined sample of 772 IgAD patients and 1,976 matched controls from 3 independent European populations. We confirmed the complex nature of the association with the HLA locus, which is the result of multiple effects spanning the entire HLA region. The primary association signal mapped to the HLA-DQB1*02 allele in the HLA Class II region (combined P = 7.69×10⁻⁵⁷; OR = 2.80) resulting from the combined independent effects of the HLA-B*0801-DRB1*0301-DQB1*02 and -DRB1*0701-DQB1*02 haplotypes, while additional secondary signals were associated with the DRB1*0102 (combined P = 5.86×10⁻¹⁷; OR = 4.28) and the DRB1*1501 (combined P = 2.24×10⁻³⁵; OR = 0.13) alleles. Despite the strong population-specific frequencies of HLA alleles, we found a remarkable conservation of these effects regardless of the ethnic background, which supports the use of large multi-ethnic populations to characterize shared genetic association signals in the HLA region. We also provide evidence for the location of association signals within the specific extended haplotypes, which will guide future sequencing studies aimed at characterizing the precise functional variants contributing to disease pathogenesis.     

Epigenomic priming of immune genes implicates oligodendroglia in multiple sclerosis susceptibility

1. Download : Download high-res image (178KB)
2. Download : Download full-size image

     

Heterogeneity for multiple disease loci in linkage analysis.

Many mendelian traits show heterogeneity; that is, the disease phenotype in different families may be caused by genes at different locations. In linkage analysis, this admixture type of heterogeneity (locus heterogeneity) has often been accommodated with one of the HOMOG programs, which thus far have been restricted to at most two disease gene locations. Here, an extension to an arbitrary number of disease locations is described. It has been implemented in a computer program, HOMOGM. This approach is also suitable as an approximation to the situation of complex traits, in which multiple disease genes may occur in the same family.     

Linkage analysis of "necessary" disease loci versus "susceptibility" loci.

The association of some diseases with specific alleles of certain genetic markers has been difficult to explain. Several explanations have been proposed for the phenomenon of association, e.g. the existence of multiple, interacting genes (epistasis) or a disease locus in linkage disequilibrium with the marker locus. One might suppose that when marker data from families with associated diseases are analyzed for linkage, the existence of the association would assure that linkage will be found, and found at a tight recombination fraction. In fact, however, linkage analyses of some diseases associated with HLA, as well as diseases associated with alleles at other loci located throughout the genome, show significant evidence against linkage, and others show loose linkage, to the puzzlement of many researchers. In part, the puzzlement arises because linkage analysis is ideal for looking for loci that are necessary, even if not sufficient, for disease expression but may be much less useful for finding loci that are neither necessary nor sufficient for disease expression (so-called susceptibility loci). This work explores what happens when one looks for linkage to susceptibility loci. A susceptibility locus in this case means that the allele increases risk but is neither necessary nor sufficient for disease expression. It might be either an allele at the marker locus itself that is increasing susceptibility or an allele at a locus in linkage disequilibrium with the marker. This work uses computer simulation to examine how linkage analyses behave when confronted with data from such a model.(ABSTRACT TRUNCATED AT 250 WORDS)     

Prevalence of celiac disease in multiple sclerosis

Background  

Celiac disease (CD) is a common systemic disease related to a permanent intolerance to gluten and is often associated with different autoimmune and neurological diseases. Its mean prevalence in the general population is 1-2% worldwide. Our aim was to study the prevalence of celiac disease in a prospective series of Multiple Sclerosis (MS) patients and their first-degree relatives.     

Combined analysis of genome-wide-linked susceptibility loci to Kawasaki disease in Han Chinese

Kawasaki disease (KD) is a dominant cause of acquired heart disease in children due to frequent complicating coronary artery lesions (CALs). Genome-wide association study and linkage analysis have recently identified 6 susceptibility loci at genome-wide significance of P < 5.0 × 10^?8 in subjects of Japanese, Taiwanese and European. In present study, we analysed the variants of 6 single nucleotide polymorphisms (SNPs) in the genetic loci to investigate their potential effect on KD susceptibility and outcomes in Han Chinese population. As a result, the risk alleles of rs1801274 and rs2254546 were observed significant effect on KD with higher frequencies in 358 patients than those in 815 controls. The significant role of rs1801274, rs2857151 and rs2254546 in KD was found in the multi-variable logistic regression analysis of the SNPs. Two 2-locus and one 3-locus combinations of the SNPs showed significant effect on KD with stronger association with KD relative to comparable single SNP or 2-locus combinations. Significant susceptibility to CALs was found in KD patients with high-risk genotypes at both rs1801274 and rs2857151. The meta-analyses first revealed significant risk for CALs in KD patients carrying risk allele of rs11340705, and the association of rs28493229 with KD was not observed in the Han Chinese. In conclusion, the findings demonstrated that 5 of the 6 genetic loci influence the risk for KD and 3 of them may be involved in secondary CALs formation in Han Chinese. The additive effects of 3 multi-locus combinations on KD/CALs imply that some loci may participate together in certain unknown gene networks related to KD/CALs. Further function studies of the genetic loci are helpful for better understanding the pathophysiology of KD.     

Functional annotation signatures of disease susceptibility loci improve SNP association analysis

Background

Genetic association studies are conducted to discover genetic loci that contribute to an inherited trait, identify the variants behind these associations and ascertain their functional role in determining the phenotype. To date, functional annotations of the genetic variants have rarely played more than an indirect role in assessing evidence for association. Here, we demonstrate how these data can be systematically integrated into an association study’s analysis plan.

Results

We developed a Bayesian statistical model for the prior probability of phenotype–genotype association that incorporates data from past association studies and publicly available functional annotation data regarding the susceptibility variants under study. The model takes the form of a binary regression of association status on a set of annotation variables whose coefficients were estimated through an analysis of associated SNPs in the GWAS Catalog (GC). The functional predictors examined included measures that have been demonstrated to correlate with the association status of SNPs in the GC and some whose utility in this regard is speculative: summaries of the UCSC Human Genome Browser ENCODE super–track data, dbSNP function class, sequence conservation summaries, proximity to genomic variants in the Database of Genomic Variants and known regulatory elements in the Open Regulatory Annotation database, PolyPhen–2 probabilities and RegulomeDB categories. Because we expected that only a fraction of the annotations would contribute to predicting association, we employed a penalized likelihood method to reduce the impact of non–informative predictors and evaluated the model’s ability to predict GC SNPs not used to construct the model. We show that the functional data alone are predictive of a SNP’s presence in the GC. Further, using data from a genome–wide study of ovarian cancer, we demonstrate that their use as prior data when testing for association is practical at the genome–wide scale and improves power to detect associations.

Conclusions

We show how diverse functional annotations can be efficiently combined to create ‘functional signatures’ that predict the a priori odds of a variant’s association to a trait and how these signatures can be integrated into a standard genome–wide–scale association analysis, resulting in improved power to detect truly associated variants.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-398) contains supplementary material, which is available to authorized users.     

Association of polymorphisms in the apolipoprotein E region with susceptibility to and progression of multiple sclerosis

Multiple sclerosis (MS) is a chronic inflammatory disorder of the central nervous system, with a complex etiology that includes a strong genetic component. The contribution of the major histocompatibility complex (MHC) has been established in numerous genetic linkage and association studies. In addition to the MHC, the chromosome 19q13 region surrounding the apolipoprotein E (APOE) gene has shown consistent evidence of involvement in MS when family-based analyses were conducted. Furthermore, several clinical reports have suggested that the APOE-4 allele may be associated with more-severe disease and faster progression of disability. To thoroughly examine the role of APOE in MS, we genotyped its functional alleles, as well as seven single-nucleotide polymorphisms (SNPs) located primarily within 13 kb of APOE, in a data set of 398 families. Using family-based association analysis, we found statistically significant evidence that an SNP haplotype near APOE is associated with MS susceptibility (P=.005). An analysis of disease progression in 614 patients with MS from 379 families indicated that APOE-4 carriers are more likely to be affected with severe disease (P=.03), whereas a higher proportion of APOE-2 carriers exhibit a mild disease course (P=.02).     

Investigation of multiple susceptibility loci for inflammatory bowel disease in an Italian cohort of patients

Background

Recent GWAs and meta-analyses have outlined about 100 susceptibility genes/loci for inflammatory bowel diseases (IBD). In this study we aimed to investigate the influence of SNPs tagging the genes/loci PTGER4, TNFSF15, NKX2-3, ZNF365, IFNG, PTPN2, PSMG1, and HLA in a large pediatric- and adult-onset IBD Italian cohort.

Methods

Eight SNPs were assessed in 1,070 Crohn''s disease (CD), 1,213 ulcerative colitis (UC), 557 of whom being diagnosed at the age of ≤16 years, and 789 healthy controls. Correlations with sub-phenotypes and major variants of NOD2 gene were investigated.

Results

The SNPs tagging the TNFSF15, NKX2-3, ZNF365, and PTPN2 genes were associated with CD (P values ranging from 0.037 to 7×10⁻⁶). The SNPs tagging the PTGER4, NKX2-3, ZNF365, IFNG, PSMG1, and HLA area were associated with UC (P values 0.047 to 4×10⁻⁵). In the pediatric cohort the associations of TNFSF15, NKX2-3 with CD, and PTGER4, NKX2-3, ZNF365, IFNG, PSMG1 with UC, were confirmed. Association with TNFSF15 and pediatric UC was also reported. A correlation with NKX2-3 and need for surgery (P  =  0.038), and with HLA and steroid-responsiveness (P  =  0.024) in UC patients was observed. Moreover, significant association in our CD cohort with TNFSF15 SNP and colonic involvement (P  =  0.021), and with ZNF365 and ileal location (P  =  0.024) was demonstrated.

Conclusions

We confirmed in a large Italian cohort the associations with CD and UC of newly identified genes, both in adult and pediatric cohort of patients, with some influence on sub-phenotypes.     

Association of polymorphisms in the HLA-B region with extended haplotypes

Genomic probes from the HLA-B region of the major histocompatibility complex (MHC) were used to study the association of restriction fragment length polymorphisms (RFLPs) with various MHC alleles, complotypes, and extended haplotypes. The two DNA probes, M20A and R5A, were derived from previously cloned cosmids and are located 38 and 110 kilobases (kb) centromeric to HLa-B, respectively. Five different RFLP variants occuring in five different haplotypic combinations were detected within a panel of 40 homozygous-typing cells and cells from 21 families using Bst EII. In two informative families with HLA-B/DR recombinations the inheritance of the RFLP variants was consistent with their mapping between HLA-B and complotypes. The R5A/M20A haplotypic pattern 6.5 kb/3.0 kb (A) had a normal Caucasian frequency of approximately 0.43 and was found in all independent examples of the extended haplotypes [HLA-B8,SC01,DR3], [HLA-B18,F1C30, DR3], [HLa-Bw62,SC33,Dr4], [HLa-B44,SC30,Dr4], and [HLA-Bw47,FC91,0DR7]. The patterns of 6.9 kb/ 3.0 kb (B), 6.5 kb/4.7 kb (C), 1.45 kb/3.0 kb (D), and 6.9 kb/4.7 kb (E) had normal Caucasian frequencies of approximately 0.23, 0.15, 0.15, and 0.04 and were found on all independent examples of [HLA-B38,SC21, DR4], [HLA-Bw57,SC61,DR7], [HLA-B7,SC31,DR2], and [HLA-B44,FC31,DR7], respectively. Individual complotypes or HLA-B alleles which were markers of extended haplotypes showed variable associations. For example, HLA-B7 and the complotype SC31 were associated with all R5A/M20A RFLP haplotypes except haplotype E, although [HLA-B7,SC31,DR7] was associated exclusively with haplotype D. HLA-B27, not known to be part of an extended haplotype, was suprisingly exclusively associated with the 6.5 kb/4.7 kb or C haplotypic pattern in all five instances tested. These findings support the concept of regional conservation of DNA on independent examples of extended haplotypes. The results also further characterize these haplotypes.     

复杂疾病的遗传易感基因区域的精细定位

以单核苷酸多态性(Single-nucleotide polymorphism, SNP)为遗传标记, 采用全基因组关联研究(Genome-wide association studies, GWAS)的策略, 已经在660多种疾病(或性状)中发现了3800多个遗传易感基因区域。但是, 其中最显著关联的遗传变异或致病性的遗传变异位点及其生物学功能并不完全清楚。这些位点的鉴定有助于阐明复杂疾病的生物学机制, 以及发现新的疾病标记物。后GWAS时代的主要任务之一就是通过精细定位研究找到复杂疾病易感基因区域内最显著关联的易感位点或致病性的易感位点并阐明其生物学功能。针对常见变异, 可通过推断或重测序增加SNP密度, 寻找最显著关联的SNP位点, 并通过功能元件分析、表达数量性状位点(Expression quantitative trait locus, eQTL)分析和单体型分析等方法寻找功能性的SNP位点和易感基因。针对罕见变异, 则可采用重测序、罕见单体型分析、家系分析和负荷检验等方法进行精细定位。文章对这些策略和所面临的问题进行了综述。     

Genome-wide association analysis reveals putative Alzheimer's disease susceptibility loci in addition to APOE

Alzheimer's disease (AD) is a genetically complex and heterogeneous disorder. To date four genes have been established to either cause early-onset autosomal-dominant AD (APP, PSEN1, and PSEN2(1-4)) or to increase susceptibility for late-onset AD (APOE5). However, the heritability of late-onset AD is as high as 80%, (6) and much of the phenotypic variance remains unexplained to date. We performed a genome-wide association (GWA) analysis using 484,522 single-nucleotide polymorphisms (SNPs) on a large (1,376 samples from 410 families) sample of AD families of self-reported European descent. We identified five SNPs showing either significant or marginally significant genome-wide association with a multivariate phenotype combining affection status and onset age. One of these signals (p = 5.7 x 10(-14)) was elicited by SNP rs4420638 and probably reflects APOE-epsilon4, which maps 11 kb proximal (r2 = 0.78). The other four signals were tested in three additional independent AD family samples composed of nearly 2700 individuals from almost 900 families. Two of these SNPs showed significant association in the replication samples (combined p values 0.007 and 0.00002). The SNP (rs11159647, on chromosome 14q31) with the strongest association signal also showed evidence of association with the same allele in GWA data generated in an independent sample of approximately 1,400 AD cases and controls (p = 0.04). Although the precise identity of the underlying locus(i) remains elusive, our study provides compelling evidence for the existence of at least one previously undescribed AD gene that, like APOE-epsilon4, primarily acts as a modifier of onset age.     

Association between polymorphisms in the promoter region of interleukin-10 and susceptibility to inflammatory bowel disease

The aim of this study was to assess the association of polymorphisms in the promoter region of the IL-10 gene with the risk of inflammatory bowel disease (IBD), including Crohn’s disease (CD) and ulcerative colitis (UC). Fifteen studies (3,693 cases and 4,574 controls) were included in a meta-analysis of association between IL-10 ?1082G/A, ?819C/T and ?592C/A polymorphisms, and IBD, CD and UC using allele contrast and the recessive, dominant, and additive models. Hardy–Weinberg equilibrium was confirmed for each study. Heterogeneity and study quality were investigated using stratification analyses and sensitivity analyses. Polymorphism ?1082G/A showed significant association with CD, with odds ratios (ORs) for the GG + GA genotype and GG versus AA genotype of 1.278 (1.004–1.627) and 1.238 (1.027–1.492) in all subjects. Significant associations were found in the Caucasian subgroup using the allele contrast, dominant, and additive models. C-allele carriers of the ?819C/T polymorphism were at increased risk of IBD (OR 1.093, 95 % CI 1.004–1.190). Association with the ?819C/T polymorphism was also found in Caucasians with CD (C vs. T: OR 1.104, 95 % CI 1.010–1.206; CC + CT vs. TT: OR 1.328, 95 % CI 1.006–1.754; CC vs. TT: OR 1.339, 95 % CI 1.008–1.778), and with UC (CC vs. CT + TT: OR 1.188, 95 % CI 1.019–1.385). No significant association was found between the ?592C/A polymorphism and IBD, CD or UC. In conclusion, the meta-analysis demonstrated clear association between the IL-10 polymorphisms ?1082G/A and ?819C/T and the risk of IBD.     

Large scale association analysis identifies three susceptibility loci for coronary artery disease

Genome wide association studies (GWAS) and their replications that have associated DNA variants with myocardial infarction (MI) and/or coronary artery disease (CAD) are predominantly based on populations of European or Eastern Asian descent. Replication of the most significantly associated polymorphisms in multiple populations with distinctive genetic backgrounds and lifestyles is crucial to the understanding of the pathophysiology of a multifactorial disease like CAD. We have used our Lebanese cohort to perform a replication study of nine previously identified CAD/MI susceptibility loci (LTA, CDKN2A-CDKN2B, CELSR2-PSRC1-SORT1, CXCL12, MTHFD1L, WDR12, PCSK9, SH2B3, and SLC22A3), and 88 genes in related phenotypes. The study was conducted on 2,002 patients with detailed demographic, clinical characteristics, and cardiac catheterization results. One marker, rs6922269, in MTHFD1L was significantly protective against MI (OR=0.68, p=0.0035), while the variant rs4977574 in CDKN2A-CDKN2B was significantly associated with MI (OR=1.33, p=0.0086). Associations were detected after adjustment for family history of CAD, gender, hypertension, hyperlipidemia, diabetes, and smoking. The parallel study of 88 previously published genes in related phenotypes encompassed 20,225 markers, three quarters of which with imputed genotypes The study was based on our genome-wide genotype data set, with imputation across the whole genome to HapMap II release 22 using HapMap CEU population as a reference. Analysis was conducted on both the genotyped and imputed variants in the 88 regions covering selected genes. This approach replicated HNRNPA3P1-CXCL12 association with CAD and identified new significant associations of CDKAL1, ST6GAL1, and PTPRD with CAD. Our study provides evidence for the importance of the multifactorial aspect of CAD/MI and describes genes predisposing to their etiology.     

Mapping of a novel susceptibility gene for rheumatoid arthritis in the telomeric MHC region

Rheumatoid arthritis (RA) is a complex heterogeneous disease with an estimated genetic contribution to of 30-50%. Approximately one third arises from the major histocompatibility complex (MHC) at 6p21.3. The contribution of specific DRB1 alleles encoding the shared epitope has been well described, however, several recent studies have suggested that additional telomeric genetic influences may exist. This region is difficult to study as a result of the presence of strong linkage disequilibrium (LD) within the MHC and high gene density particularly in the central class III region. In this article we review the current data supporting the existence of a non-DRB1 susceptibility gene for rheumatoid arthritis, in particular within the class III region.     

Analysis of a marsupial MHC region containing two recently duplicated class I loci

A 37-kb cosmid containing two complete major histocompatibility complex (MHC) class I chain loci from the opossum Monodelphis domestica was isolated, fully sequenced, and characterized. This sequence represents the largest contiguous genomic sequence reported for the MHC region of a nonplacental mammal. Based on particular conserved amino acid residues, and limited expression analyses, the two MHC-I loci, designated ModoUB and ModoUC, appear to encode functional MHC-I molecules. The two coding regions are 98% identical at the nucleotide level; however, their promoter regions differ significantly. Two CpG islands present in the cosmid sequence correspond to the two coding regions. Twelve microsatellites and six retroelements were also present in the cosmid. The retroelements share highest sequence homology to the CORE–SINE family of retroelements. Due to high sequence identity, it is very likely that ModoUB and ModoUC loci are products of recent gene duplication that occurred less than 4 million years ago.     

Detection of susceptibility loci by genome-wide linkage analysis

The objective of this study is to evaluate the efficacy of a model-free linkage statistics for finding evidence of linkage using two different maps and to illustrate how the comparison of results from several populations might provide insight into the underlying genetic etiology of the disease of interest. The results obtained in terms of detection of the risk loci and threshold for declaring linkage and power are very similar for a dense SNP map and a sparser microsatellite map. The populations differed in terms of family ascertainment and diagnosis criteria, leading to different power to detect the individual underlying disease loci. Our results for the individual replicates are consistent with the disease model used in the simulation.     

Presence of multiple independent effects in risk loci of common complex human diseases

Many genetic loci and SNPs associated with many common complex human diseases and traits are now identified. The total genetic variance explained by these loci for a trait or disease, however, has often been very small. Much of the "missing heritability" has been revealed to be hidden in the genome among the large number of variants with small effects. Several recent studies have reported the presence of multiple independent SNPs and genetic heterogeneity in trait-associated loci. It is therefore reasonable to speculate that such a phenomenon could be common among loci known to be associated with a complex trait or disease. For testing this hypothesis, a total of 117 loci known to be associated with rheumatoid arthritis (RA), Crohn disease (CD), type 1 diabetes (T1D), or type 2 diabetes (T2D) were selected. The presence of multiple independent effects was assessed in the case-control samples genotyped by the Wellcome Trust Case Control Consortium study and imputed with SNP genotype information from the HapMap Project and the 1000 Genomes Project. Eleven loci with evidence of multiple independent effects were identified in the study, and the number was expected to increase at larger sample sizes and improved statistical power. The variance explained by the multiple effects in a locus was much higher than the variance explained by the single reported SNP effect. The results thus significantly improve our understanding of the allelic structure of these individual disease-associated loci, as well as our knowledge of the general genetic mechanisms of common complex traits and diseases.     

Genetic identification of multiple loci that control breast cancer susceptibility in the rat.

We have used a rat model of induced mammary carcinomas in an effort to identify breast cancer susceptibility genes. Using genetic crosses between the carcinoma-resistant Copenhagen (COP) and carcinoma-sensitive Wistar-Furth rats, we have confirmed the identification of the Mcs1 locus that modulates tumor number. We have now also identified two additional loci, Mcs2 and Mcs3. These three loci map to chromosomes 2, 7, and 1, respectively, and interact additively to suppress mammary carcinoma development in the COP strain. They are responsible for a major portion of the tumor-resistant phenotype of the COP rat. No loss of heterozygosity was observed surrounding the three loci. A fourth COP locus, Mcs4, has also been identified on chromosome 8 and acts in contrast to increase the number of carcinomas. These results show that mammary carcinoma susceptibility in the COP rat is a polygenic trait. Interestingly, a polymorphism in the human genomic region homologous to the rat Mcs4 region is associated with an increased breast cancer risk in African-American women. The isolation of the Mcs genes may help elucidate novel mechanisms of carcinogenesis, provide information important for human breast cancer risk estimation, and also provide unique drug discovery targets for breast cancer prevention.