Epistatic interaction between KIR3DS1 and HLA-B delays the progression to AIDS

Natural killer (NK) cells provide defense in the early stages of the innate immune response against viral infections by producing cytokines and causing cytotoxicity. The killer immunoglobulin-like receptors (KIRs) on NK cells regulate the inhibition and activation of NK-cell responses through recognition of human leukocyte antigen (HLA) class I molecules on target cells KIR and HLA loci are both highly polymorphic, and some HLA class I products bind and trigger cell-surface receptors specified by KIR genes. Here we report that the activating KIR allele KIR3DS1, in combination with HLA-B alleles that encode molecules with isoleucine at position 80 (HLA-B Bw4-80Ile), is associated with delayed progression to AIDS in individuals infected with human immunodeficiency virus type 1 (HIV-1). In the absence of KIR3DS1, the HLA-B Bw4-80Ile allele was not associated with any of the AIDS outcomes measured. By contrast, in the absence of HLA-B Bw4-80Ile alleles, KIR3DS1 was significantly associated with more rapid progression to AIDS. These observations are strongly suggestive of a model involving an epistatic interaction between the two loci. The strongest synergistic effect of these loci was on progression to depletion of CD4(+) T cells, which suggests that a protective response of NK cells involving KIR3DS1 and its HLA class I ligands begins soon after HIV-1 infection.     

Innate partnership of HLA-B and KIR3DL1 subtypes against HIV-1

Allotypes of the natural killer (NK) cell receptor KIR3DL1 vary in both NK cell expression patterns and inhibitory capacity upon binding to their ligands, HLA-B Bw4 molecules, present on target cells. Using a sample size of over 1,500 human immunodeficiency virus (HIV)+ individuals, we show that various distinct allelic combinations of the KIR3DL1 and HLA-B loci significantly and strongly influence both AIDS progression and plasma HIV RNA abundance in a consistent manner. These genetic data correlate very well with previously defined functional differences that distinguish KIR3DL1 allotypes. The various epistatic effects observed here for common, distinct KIR3DL1 and HLA-B Bw4 combinations are unprecedented with regard to any pair of genetic loci in human disease, and indicate that NK cells may have a critical role in the natural history of HIV infection.     

Global diversity and evidence for coevolution of KIR and HLA

The killer immunoglobulin-like receptor (KIR) gene cluster shows extensive genetic diversity, as do the HLA class I loci, which encode ligands for KIR molecules. We genotyped 1,642 individuals from 30 geographically distinct populations to examine population-level evidence for coevolution of these two functionally related but unlinked gene clusters. We observed strong negative correlations between the presence of activating KIR genes and their corresponding HLA ligand groups across populations, especially KIR3DS1 and its putative HLA-B Bw4-80I ligands (r = -0.66, P = 0.038). In contrast, we observed weak positive relationships between the various inhibitory KIR genes and their ligands. We observed a negative correlation between distance from East Africa and frequency of activating KIR genes and their corresponding ligands, suggesting a balance between selection on HLA and KIR loci. Most KIR-HLA genetic association studies indicate a primary influence of activating KIR-HLA genotypes in disease risk; concomitantly, activating receptor-ligand pairs in this study show the strongest signature of coevolution of these two complex genetic systems as compared with inhibitory receptor-ligand pairs.     

Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility

Ankylosing spondylitis is a common form of inflammatory arthritis predominantly affecting the spine and pelvis that occurs in approximately 5 out of 1,000 adults of European descent. Here we report the identification of three variants in the RUNX3, LTBR-TNFRSF1A and IL12B regions convincingly associated with ankylosing spondylitis (P < 5 × 10(-8) in the combined discovery and replication datasets) and a further four loci at PTGER4, TBKBP1, ANTXR2 and CARD9 that show strong association across all our datasets (P < 5 × 10(-6) overall, with support in each of the three datasets studied). We also show that polymorphisms of ERAP1, which encodes an endoplasmic reticulum aminopeptidase involved in peptide trimming before HLA class I presentation, only affect ankylosing spondylitis risk in HLA-B27-positive individuals. These findings provide strong evidence that HLA-B27 operates in ankylosing spondylitis through a mechanism involving aberrant processing of antigenic peptides.     

Adaptive evolution of the tumour suppressor BRCA1 in humans and chimpanzees. Australian Breast Cancer Family Study

Mutations in BRCA1 (ref. 1) confer an increased risk of female breast cancer. In a genome-wide scan of linkage disequilibrium (LD), a high level of LD was detected among microsatellite markers flanking BRCA1 (ref. 3), raising the prospect that positive natural selection may have acted on this gene. We have used the predictions of evolutionary genetic theory to investigate this further. Using phylogeny-based maximum likelihood analysis of the BRCA1 sequences from primates and other mammals, we found that the ratios of replacement to silent nucleotide substitutions on the human and chimpanzee lineages were not different from one another (P=0.8), were different from those of other primate lineages (P=0.004) and were greater than 1 (P=0.04). This is consistent with the historic occurrence of positive darwinian selection pressure on the BRCA1 protein in the human and chimpanzee lineages. Analysis of genetic variation in a sample of female Australians of Northern European origin showed evidence for Hardy-Weinberg (HW) disequilibrium at polymorphic sites in BRCA1, consistent with the possibility that natural selection is affecting genotype frequencies in modern Europeans. The clustering of between-species variation in the region of the gene encoding the RAD51-interaction domain of BRCA1 suggests the maintenance of genomic integrity as a possible target of selection.     

Epistasis between mouse Klra and major histocompatibility complex class I loci is associated with a new mechanism of natural killer cell-mediated innate resistance to cytomegalovirus infection

Experimental infection with mouse cytomegalovirus (MCMV) has been used to elucidate the intricate host-pathogen mechanisms that determine innate resistance to infection. Linkage analyses in F(2) progeny from MCMV-resistant MA/My (H2 (k)) and MCMV-susceptible BALB/c (H2 (d)) and BALB.K (H2 (k)) mouse strains indicated that only the combination of alleles encoded by a gene in the Klra (also called Ly49) cluster on chromosome 6, and one in the major histocompatibility complex (H2) on chromosome 17, is associated with virus resistance. We found that natural killer cell-activating receptor Ly49P specifically recognized MCMV-infected cells, dependent on the presence of the H2 (k) haplotype. This binding was blocked using antibodies to H-2D(k) but not antibodies to H-2K(k). These results are suggestive of a new natural killer cell mechanism implicated in MCMV resistance, which depends on the functional interaction of the Ly49P receptor and the major histocompatibility complex class I molecule H-2D(k) on MCMV-infected cells.     

A putative RUNX1 binding site variant between SLC9A3R1 and NAT9 is associated with susceptibility to psoriasis

Psoriasis (OMIM 177900) is a chronic inflammatory skin disorder of unknown pathogenesis affecting approximately 2% of the Western population. It occurs more frequently in individuals with human immunodeficiency virus, and 20-30% of individuals with psoriasis have psoriatic arthritis. Psoriasis is associated with HLA class I alleles, and previous linkage analysis by our group identified a second psoriasis locus at 17q24-q25 (PSORS2; ref. 7). Linkage to this locus was confirmed with independent family sets. Additional loci have also been proposed to be associated with psoriasis. Here we describe two peaks of strong association with psoriasis on chromosome 17q25 separated by 6 Mb. Associated single-nucleotide polymorphisms (SNPs) in the proximal peak lie in or near SLC9A3R1 (also called EBP50 and NHERF1) and NAT9, a new member of the N-acetyltransferase family. SLC9A3R1 is a PDZ domain-containing phosphoprotein that associates with members of the ezrin-radixin-moesin family and is implicated in diverse aspects of epithelial membrane biology and immune synapse formation in T cells. The distal peak of association is in RAPTOR (p150 target of rapamycin (TOR)-scaffold protein containing WD-repeats). Expression of SLC9A3R1 is highest in the uppermost stratum Malpighi of psoriatic and normal skin and in inactive versus active T cells. A disease-associated SNP lying between SLC9A3R1 and NAT9 leads to loss of RUNX1 binding. This is the second example of loss of a RUNX1 binding site associated with susceptibility to an autoimmune disease. It also suggests defective regulation of SLC9A3R1 or NAT9 by RUNX1 as a susceptibility factor for psoriasis.     

Linkage disequilibrium of a type 1 diabetes susceptibility locus with a regulatory IL12B allele

Type 1 diabetes (T1D; or insulin-dependent diabetes mellitus, IDDM) is an autoimmune disease with both genetic and environmental components. In addition to the human leukocyte antigen (HLA) complex, the single major genetic contributor of susceptibility, an unknown number of other unidentified genes are required to mediate disease. Although many loci conferring susceptibility to T1D have been mapped, their identification has proven problematic due to the complex nature of this disease. Our strategy for finding T1D susceptibility genes has been to test for human homologues of loci implicated in diabetes-prone NOD (non-obese diabetic) mice, together with application of biologically relevant stratification methods. We report here a new susceptibility locus, IDDM18, located near the interleukin-12 (IL-12)p40 gene, IL12B. Significant bias in transmission of IL12B alleles was observed in affected sibpairs and was confirmed in an independent cohort of simplex families. A single base change in the 3' UTR showed strong linkage disequilibrium with the T1D susceptibility locus. The IL12B 3' UTR alleles showed different levels of expression in cell lines. Variation in IL-12p40 production may influence T-cell responses crucial for either mediating or protecting against this and other autoimmune diseases.     

Rare MTNR1B variants impairing melatonin receptor 1B function contribute to type 2 diabetes

Genome-wide association studies have revealed that common noncoding variants in MTNR1B (encoding melatonin receptor 1B, also known as MT(2)) increase type 2 diabetes (T2D) risk(1,2). Although the strongest association signal was highly significant (P < 1 × 10(-20)), its contribution to T2D risk was modest (odds ratio (OR) of ～1.10-1.15)(1-3). We performed large-scale exon resequencing in 7,632 Europeans, including 2,186 individuals with T2D, and identified 40 nonsynonymous variants, including 36 very rare variants (minor allele frequency (MAF) <0.1%), associated with T2D (OR = 3.31, 95% confidence interval (CI) = 1.78-6.18; P = 1.64 × 10(-4)). A four-tiered functional investigation of all 40 mutants revealed that 14 were non-functional and rare (MAF < 1%), and 4 were very rare with complete loss of melatonin binding and signaling capabilities. Among the very rare variants, the partial- or total-loss-of-function variants but not the neutral ones contributed to T2D (OR = 5.67, CI = 2.17-14.82; P = 4.09 × 10(-4)). Genotyping the four complete loss-of-function variants in 11,854 additional individuals revealed their association with T2D risk (8,153 individuals with T2D and 10,100 controls; OR = 3.88, CI = 1.49-10.07; P = 5.37 × 10(-3)). This study establishes a firm functional link between MTNR1B and T2D risk.     

A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC

The proteins encoded by the classical HLA class I and class II genes in the major histocompatibility complex (MHC) are highly polymorphic and are essential in self versus non-self immune recognition. HLA variation is a crucial determinant of transplant rejection and susceptibility to a large number of infectious and autoimmune diseases. Yet identification of causal variants is problematic owing to linkage disequilibrium that extends across multiple HLA and non-HLA genes in the MHC. We therefore set out to characterize the linkage disequilibrium patterns between the highly polymorphic HLA genes and background variation by typing the classical HLA genes and >7,500 common SNPs and deletion-insertion polymorphisms across four population samples. The analysis provides informative tag SNPs that capture much of the common variation in the MHC region and that could be used in disease association studies, and it provides new insight into the evolutionary dynamics and ancestral origins of the HLA loci and their haplotypes.     

Genome-wide association study identifies susceptibility loci for dengue shock syndrome at MICB and PLCE1

Hypovolemic shock (dengue shock syndrome (DSS)) is the most common life-threatening complication of dengue. We conducted a genome-wide association study of 2,008 pediatric cases treated for DSS and 2,018 controls from Vietnam. Replication of the most significantly associated markers was carried out in an independent Vietnamese sample of 1,737 cases and 2,934 controls. SNPs at two loci showed genome-wide significant association with DSS. We identified a susceptibility locus at MICB (major histocompatibility complex (MHC) class I polypeptide-related sequence B), which was within the broad MHC region on chromosome 6 but outside the class I and class II HLA loci (rs3132468, P(meta) = 4.41 × 10(-11), per-allele odds ratio (OR) = 1.34 (95% confidence interval: 1.23-1.46)). We identified associated variants within PLCE1 (phospholipase C, epsilon 1) on chromosome 10 (rs3765524, P(meta) = 3.08 × 10(-10), per-allele OR = 0.80 (95% confidence interval: 0.75-0.86)). We identify two loci associated with susceptibility to DSS in people with dengue, suggesting possible mechanisms for this severe complication of dengue.     

Loss-of-function mutations in TYROBP (DAP12) result in a presenile dementia with bone cysts

Polycystic lipomembranous osteodysplasia with sclerosing leukoencephalopathy (PLOSL; MIM 221770), also known as Nasu-Hakola disease, is a recessively inherited disease characterized by a combination of psychotic symptoms rapidly progressing to presenile dementia and bone cysts restricted to wrists and ankles. PLOSL has a global distribution, although most of the patients have been diagnosed in Finland and Japan, with an estimated population prevalence of 2x10-6 (ref. 2) in the Finns. We have previously identified a shared 153-kb ancestor haplotype in all Finnish disease alleles between markers D19S1175 and D19S608 on chromosome 19q13.1 (refs 5,6). Here we characterize the molecular defect in PLOSL by identifying one large deletion in all Finnish PLOSL alleles and another mutation in a Japanese patient, both representing loss-of-function mutations, in the gene encoding TYRO protein tyrosine kinase binding protein (TYROBP; formerly DAP12). TYROBP is a transmembrane protein that has been recognized as a key activating signal transduction element in natural killer (NK) cells. On the plasma membrane of NK cells, TYROBP associates with activating receptors recognizing major histocompatibility complex (MHC) class I molecules. No abnormalities in NK cell function were detected in PLOSL patients homozygous for a null allele of TYROBP.     

A predominant role for the HLA class II region in the association of the MHC region with multiple sclerosis

Genetic susceptibility to multiple sclerosis is associated with genes of the major histocompatibility complex (MHC), particularly HLA-DRB1 and HLA-DQB1 (ref. 1). Both locus and allelic heterogeneity have been reported in this genomic region. To clarify whether HLA-DRB1 itself, nearby genes in the region encoding the MHC or combinations of these loci underlie susceptibility to multiple sclerosis, we genotyped 1,185 Canadian and Finnish families with multiple sclerosis (n = 4,203 individuals) with a high-density SNP panel spanning the genes encoding the MHC and flanking genomic regions. Strong associations in Canadian and Finnish samples were observed with blocks in the HLA class II genomic region (P < 4.9 x 10(-13) and P < 2.0 x 10(-16), respectively), but the strongest association was with HLA-DRB1 (P < 4.4 x 10(-17)). Conditioning on either HLA-DRB1 or the most significant HLA class II haplotype block found no additional block or SNP association independent of the HLA class II genomic region. This study therefore indicates that MHC-associated susceptibility to multiple sclerosis is determined by HLA class II alleles, their interactions and closely neighboring variants.