Unusual selection on the KIR3DL1/S1 natural killer cell receptor in Africans

Interactions of killer cell immunoglobulin-like receptors (KIRs) with major histocompatibility complex (MHC) class I ligands diversify natural killer cell responses to infection. By analyzing sequence variation in diverse human populations, we show that the KIR3DL1/S1 locus encodes two lineages of polymorphic inhibitory KIR3DL1 allotypes that recognize Bw4 epitopes of protein">HLA-A and HLA-B and one lineage of conserved activating KIR3DS1 allotypes, also implicated in Bw4 recognition. Balancing selection has maintained these three lineages for over 3 million years. Variation was selected at D1 and D2 domain residues that contact HLA class I and at two sites on D0, the domain that enhances the binding of KIR3D to HLA class I. HLA-B variants that gained Bw4 through interallelic microconversion are also products of selection. A worldwide comparison uncovers unusual KIR3DL1/S1 evolution in modern sub-Saharan Africans. Balancing selection is weak and confined to D0, KIR3DS1 is rare and KIR3DL1 allotypes with similar binding sites predominate. Natural killer cells express the dominant KIR3DL1 at a high frequency and with high surface density, providing strong responses to cells perturbed in Bw4 expression.     

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.     

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.     

Role for DNA methylation in the control of cell type specific maspin expression

The major histocompatibility complex (MHC) is a source of unique individual odors that influence individual recognition, mating preferences, nesting behavior and selective block of pregnancy in animals. Such phenomena have been difficult to study in humans, because the human leukocyte antigen (HLA, human MHC) loci are the most polymorphic loci in the human genome, with the potential to generate millions of unique combinations of genotypes. In addition, high variability in background odors, encoded by the rest of the genome and influenced by cultural practices, contribute to a low signal-to-noise ratio that could mask HLA-based olfactory cues. Here we show that women can detect differences of one HLA allele among male odor donors with different MHC genotypes. Notably, the mechanism for a woman's ability to discriminate and choose odors is based on HLA alleles inherited from her father but not her mother. The parents' HLA alleles that she does not inherit show no relationship with odor choice, despite exposure to these HLA-encoded odors throughout her life. Our data indicate that paternally inherited HLA-associated odors influence odor preference and may serve as social cues.     

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.     

Paternal transmission of the very common class I INS VNTR alleles predisposes to childhood obesity

To identify some of the genetic factors that contribute to obesity in children of Central European and North African descent, we studied the parental transmission of alleles at the insulin locus to offspring with early-onset obesity. A variable nucleotide tandem repeat (VNTR) polymorphism upstream of the insulin gene (INS) is associated with variations in the expression of INS and the nearby gene encoding insulin-like growth factor 2 (IGF2). We found an excess of paternal transmission of class I VNTR alleles to obese children: children who inherited a class I allele from their father (but not those inheriting it from their mother) had a relative risk of early-onset obesity of 1.8. Due to the frequency of class I alleles in this population, this risk concerns 65-70% of all infants. These results suggest that increased in utero expression of paternal INS or IGF2 due to the class I INS VNTR allele may predispose offspring to postnatal fat deposition.     

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.     

Common variants in P2RY11 are associated with narcolepsy

Growing evidence supports the hypothesis that narcolepsy with cataplexy is an autoimmune disease. We here report genome-wide association analyses for narcolepsy with replication and fine mapping across three ethnic groups (3,406 individuals of European ancestry, 2,414 Asians and 302 African Americans). We identify a SNP in the 3' untranslated region of P2RY11, the purinergic receptor subtype P2Y?? gene, which is associated with narcolepsy (rs2305795, combined P = 6.1 × 10?1?, odds ratio = 1.28, 95% CI 1.19-1.39, n = 5689). The disease-associated allele is correlated with reduced expression of P2RY11 in CD8(+) T lymphocytes (339% reduced, P = 0.003) and natural killer (NK) cells (P = 0.031), but not in other peripheral blood mononuclear cell types. The low expression variant is also associated with reduced P2RY11-mediated resistance to ATP-induced cell death in T lymphocytes (P = 0.0007) and natural killer cells (P = 0.001). These results identify P2RY11 as an important regulator of immune-cell survival, with possible implications in narcolepsy and other autoimmune diseases.     

Genome-wide association study identifies 12 new susceptibility loci for primary biliary cirrhosis

In addition to the HLA locus, six genetic risk factors for primary biliary cirrhosis (PBC) have been identified in recent genome-wide association studies (GWAS). To identify additional loci, we carried out a GWAS using 1,840 cases from the UK PBC Consortium and 5,163 UK population controls as part of the Wellcome Trust Case Control Consortium 3 (WTCCC3). We followed up 28 loci in an additional UK cohort of 620 PBC cases and 2,514 population controls. We identified 12 new susceptibility loci (at a genome-wide significance level of P < 5 × 10??) and replicated all previously associated loci. We identified three further new loci in a meta-analysis of data from our study and previously published GWAS results. New candidate genes include STAT4, DENND1B, CD80, IL7R, CXCR5, TNFRSF1A, CLEC16A and NFKB1. This study has considerably expanded our knowledge of the genetic architecture of PBC.     

Identification of Stk6/STK15 as a candidate low-penetrance tumor-susceptibility gene in mouse and human

Linkage analysis and haplotype mapping in interspecific mouse crosses (Mus musculus x Mus spretus) identified the gene encoding Aurora2 (Stk6 in mouse and STK15 in human) as a candidate skin tumor susceptibility gene. The Stk6 allele inherited from the susceptible M. musculus parent was overexpressed in normal cells and preferentially amplified in tumor cells from F(1) hybrid mice. We identified a common genetic variant in STK15 (resulting in the amino acid substitution F31I) that is preferentially amplified and associated with the degree of aneuploidy in human colon tumors. The Ile31 variant transforms rat1 cells more potently than the more common Phe31 variant. The E2 ubiquitin-conjugating enzyme UBE2N was a preferential binding partner of the 'weak' STK15 Phe31 variant form in yeast two-hybrid screens and in human cells. This interaction results in colocalization of UBE2N with STK15 at the centrosomes during mitosis. These results are consistent with an important role for the Ile31 variant of STK15 in human cancer susceptibility.     

Heterogeneous mutation processes in human microsatellite DNA sequences

Although microsatellite polymorphisms are one of the most commonly used tools in genetic analyses, it remains to be understood how microsatellite DNA has evolved as a ubiquitous and highly abundant class of repetitive sequences in eukaryotic genomes. On the basis of analyses of spontaneous human microsatellite mutations of germline origin, I show here that different mutation biases underlie the evolution of microsatellite repeats. The within-locus mutation rate increases with allele length, but is not affected by the size difference between an individual's two alleles (allele span). Within loci, long alleles tend to mutate to shorter lengths, thereby acting to prevent infinite growth. Expansions are more common than contractions among dinucleotide repeats, whereas no such trend is evident among tetranucleotide repeats. This observation is consistent with the longer repeat lengths and higher frequency of di- compared with tetranucleotide repeats. An excess of paternally transmitted mutations (male-to-female ratio of 4.9) supports a male-biased mutation rate in the human genome.     

TLR4 mutations are associated with endotoxin hyporesponsiveness in humans

There is much variability between individuals in the response to inhaled toxins, but it is not known why certain people develop disease when challenged with environmental agents and others remain healthy. To address this, we investigated whether TLR4 (encoding the toll-like receptor-4), which has been shown to affect lipopolysaccharide (LPS) responsiveness in mice, underlies the variability in airway responsiveness to inhaled LPS in humans. Here we show that common, co-segregating missense mutations (Asp299Gly and Thr399Ile) affecting the extracellular domain of the TLR4 receptor are associated with a blunted response to inhaled LPS in humans. Transfection of THP-1 cells demonstrates that the Asp299Gly mutation (but not the Thr399Ile mutation) interrupts TLR4-mediated LPS signalling. Moreover, the wild-type allele of TLR4 rescues the LPS hyporesponsive phenotype in either primary airway epithelial cells or alveolar macrophages obtained from individuals with the TLR4 mutations. Our findings provide the first genetic evidence that common mutations in TLR4 are associated with differences in LPS responsiveness in humans, and demonstrate that gene-sequence changes can alter the ability of the host to respond to environmental stress.     

A viable allele of Mcm4 causes chromosome instability and mammary adenocarcinomas in mice

Mcm4 (minichromosome maintenance-deficient 4 homolog) encodes a subunit of the MCM2-7 complex (also known as MCM2-MCM7), the replication licensing factor and presumptive replicative helicase. Here, we report that the mouse chromosome instability mutation Chaos3 (chromosome aberrations occurring spontaneously 3), isolated in a forward genetic screen, is a viable allele of Mcm4. Mcm4(Chaos3) encodes a change in an evolutionarily invariant amino acid (F345I), producing an apparently destabilized MCM4. Saccharomyces cerevisiae strains that we engineered to contain a corresponding allele (resulting in an F391I change) showed a classical minichromosome loss phenotype. Whereas homozygosity for a disrupted Mcm4 allele (Mcm4(-)) caused preimplantation lethality, Mcm(Chaos3/-) embryos died late in gestation, indicating that Mcm4(Chaos3) is hypomorphic. Mutant embryonic fibroblasts were highly susceptible to chromosome breaks induced by the DNA replication inhibitor aphidicolin. Most notably, >80% of Mcm4(Chaos3/Chaos3) females succumbed to mammary adenocarcinomas with a mean latency of 12 months. These findings suggest that hypomorphic alleles of the genes encoding the subunits of the MCM2-7 complex may increase breast cancer risk.     

A genome-wide association study identifies new psoriasis susceptibility loci and an interaction between HLA-C and ERAP1

To identify new susceptibility loci for psoriasis, we undertook a genome-wide association study of 594,224 SNPs in 2,622 individuals with psoriasis and 5,667 controls. We identified associations at eight previously unreported genomic loci. Seven loci harbored genes with recognized immune functions (IL28RA, REL, IFIH1, ERAP1, TRAF3IP2, NFKBIA and TYK2). These associations were replicated in 9,079 European samples (six loci with a combined P < 5 × 10?? and two loci with a combined P < 5 × 10??). We also report compelling evidence for an interaction between the HLA-C and ERAP1 loci (combined P = 6.95 × 10??). ERAP1 plays an important role in MHC class I peptide processing. ERAP1 variants only influenced psoriasis susceptibility in individuals carrying the HLA-C risk allele. Our findings implicate pathways that integrate epidermal barrier dysfunction with innate and adaptive immune dysregulation in psoriasis pathogenesis.     

Genome-wide association study identifies susceptibility loci for IgA nephropathy

We carried out a genome-wide association study of IgA nephropathy, a major cause of kidney failure worldwide. We studied 1,194 cases and 902 controls of Chinese Han ancestry, with targeted follow up in Chinese and European cohorts comprising 1,950 cases and 1,920 controls. We identified three independent loci in the major histocompatibility complex, as well as a common deletion of CFHR1 and CFHR3 at chromosome 1q32 and a locus at chromosome 22q12 that each surpassed genome-wide significance (P values for association between 1.59 × 10?2? and 4.84 × 10?? and minor allele odds ratios of 0.63-0.80). These five loci explain 4-7% of the disease variance and up to a tenfold variation in interindividual risk. Many of the alleles that protect against IgA nephropathy impart increased risk for other autoimmune or infectious diseases, and IgA nephropathy risk allele frequencies closely parallel the variation in disease prevalence among Asian, European and African populations, suggesting complex selective pressures.     

Failure to confirm NOTCH4 association with schizophrenia in a large population-based sample from Scotland

The NOTCH4 gene was recently reported to be associated with schizophrenia based on TDT analysis of 80 British trios. The strongest evidence for association derived from two microsatellites. We genotyped both loci in a large sample of unrelated Scottish schizophrenics and controls, but failed to replicate the reported association, finding instead that each putative schizophrenia-associated allele had a somewhat lower frequency in schizophrenics than in controls.