Variation in interleukin 7 receptor alpha chain (IL7R) influences risk of multiple sclerosis

Multiple sclerosis is a chronic, often disabling, disease of the central nervous system affecting more than 1 in 1,000 people in most western countries. The inflammatory lesions typical of multiple sclerosis show autoimmune features and depend partly on genetic factors. Of these genetic factors, only the HLA gene complex has been repeatedly confirmed to be associated with multiple sclerosis, despite considerable efforts. Polymorphisms in a number of non-HLA genes have been reported to be associated with multiple sclerosis, but so far confirmation has been difficult. Here, we report compelling evidence that polymorphisms in IL7R, which encodes the interleukin 7 receptor alpha chain (IL7Ralpha), indeed contribute to the non-HLA genetic risk in multiple sclerosis, demonstrating a role for this pathway in the pathophysiology of this disease. In addition, we report altered expression of the genes encoding IL7Ralpha and its ligand, IL7, in the cerebrospinal fluid compartment of individuals with multiple sclerosis.     

Oncogenic IL7R gain-of-function mutations in childhood T-cell acute lymphoblastic leukemia

Interleukin 7 (IL-7) and its receptor, formed by IL-7Rα (encoded by IL7R) and γc, are essential for normal T-cell development and homeostasis. Here we show that IL7R is an oncogene mutated in T-cell acute lymphoblastic leukemia (T-ALL). We find that 9% of individuals with T-ALL have somatic gain-of-function IL7R exon 6 mutations. In most cases, these IL7R mutations introduce an unpaired cysteine in the extracellular juxtamembrane-transmembrane region and promote de novo formation of intermolecular disulfide bonds between mutant IL-7Rα subunits, thereby driving constitutive signaling via JAK1 and independently of IL-7, γc or JAK3. IL7R mutations induce a gene expression profile partially resembling that provoked by IL-7 and are enriched in the T-ALL subgroup comprising TLX3 rearranged and HOXA deregulated cases. Notably, IL7R mutations promote cell transformation and tumor formation. Overall, our findings indicate that IL7R mutational activation is involved in human T-cell leukemogenesis, paving the way for therapeutic targeting of IL-7R-mediated signaling in T-ALL.     

Genetic determinants of ulcerative colitis include the ECM1 locus and five loci implicated in Crohn's disease

We report results of a nonsynonymous SNP scan for ulcerative colitis and identify a previously unknown susceptibility locus at ECM1. We also show that several risk loci are common to ulcerative colitis and Crohn's disease (IL23R, IL12B, HLA, NKX2-3 and MST1), whereas autophagy genes ATG16L1 and IRGM, along with NOD2 (also known as CARD15), are specific for Crohn's disease. These data provide the first detailed illustration of the genetic relationship between these common inflammatory bowel diseases.     

Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47

Genome-wide association studies and candidate gene studies in ulcerative colitis have identified 18 susceptibility loci. We conducted a meta-analysis of six ulcerative colitis genome-wide association study datasets, comprising 6,687 cases and 19,718 controls, and followed up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P < 5 × 10(-8)), increasing the number of ulcerative colitis-associated loci to 47. After annotating associated regions using GRAIL, expression quantitative trait loci data and correlations with non-synonymous SNPs, we identified many candidate genes that provide potentially important insights into disease pathogenesis, including IL1R2, IL8RA-IL8RB, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease risk loci is now 99, including a minimum of 28 shared association signals between Crohn's disease and ulcerative colitis.     

Sequence variants in IL10, ARPC2 and multiple other loci contribute to ulcerative colitis susceptibility

Inflammatory bowel disease (IBD) typically manifests as either ulcerative colitis (UC) or Crohn's disease (CD). Systematic identification of susceptibility genes for IBD has thus far focused mainly on CD, and little is known about the genetic architecture of UC. Here we report a genome-wide association study with 440,794 SNPs genotyped in 1,167 individuals with UC and 777 healthy controls. Twenty of the most significantly associated SNPs were tested for replication in three independent European case-control panels comprising a total of 1,855 individuals with UC and 3,091 controls. Among the four consistently replicated markers, SNP rs3024505 immediately flanking the IL10 (interleukin 10) gene on chromosome 1q32.1 showed the most significant association in the combined verification samples (P = 1.35 x 10(-12); OR = 1.46 (1.31-1.62)). The other markers were located in ARPC2 and in the HLA-BTNL2 region. Association between rs3024505 and CD (1,848 cases, 1,804 controls) was weak (P = 0.013; OR = 1.17 (1.01-1.34)). IL10 is an immunosuppressive cytokine that has long been proposed to influence IBD pathophysiology. Our findings strongly suggest that defective IL10 function is central to the pathogenesis of the UC subtype of IBD.     

Genetic variation in DPP6 is associated with susceptibility to amyotrophic lateral sclerosis

We identified a SNP in the DPP6 gene that is consistently strongly associated with susceptibility to amyotrophic lateral sclerosis (ALS) in different populations of European ancestry, with an overall P value of 5.04 x 10(-8) in 1,767 cases and 1,916 healthy controls and with an odds ratio of 1.30 (95% confidence interval (CI) of 1.18-1.43). Our finding is the first report of a genome-wide significant association with sporadic ALS and may be a target for future functional studies.     

Identification of mutations in CUL7 in 3-M syndrome

Intrauterine growth retardation is caused by maternal, fetal or placental factors that result in impaired endovascular trophoblast invasion and reduced placental perfusion. Although various causes of intrauterine growth retardation have been identified, most cases remain unexplained. Studying 29 families with 3-M syndrome (OMIM 273750), an autosomal recessive condition characterized by severe pre- and postnatal growth retardation, we first mapped the underlying gene to chromosome 6p21.1 and then identified 25 distinct mutations in the gene cullin 7 (CUL7). CUL7 assembles an E3 ubiquitin ligase complex containing Skp1, Fbx29 (also called Fbw8) and ROC1 and promotes ubiquitination. Using deletion analysis, we found that CUL7 uses its central region to interact with the Skp1-Fbx29 heterodimer. Functional studies indicated that the 3-M-associated CUL7 nonsense and missense mutations R1445X and H1464P, respectively, render CUL7 deficient in recruiting ROC1. These results suggest that impaired ubiquitination may have a role in the pathogenesis of intrauterine growth retardation in humans.     

Mutation of the ST7 tumor suppressor gene on 7q31.1 is rare in breast, ovarian and colorectal cancers.

The gene ST7 has recently been implicated as the broad-range tumor suppressor on human chromosome 7q31.1. We did not detect somatic mutations in ST7 in any of 149 primary ovarian, breast or colon carcinomas. These data suggest that epigenetic downregulation or haploinsufficiency, rather than somatic genetic alterations, may be the primary mechanism of abrogation of ST7 function in these tumor types.     

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.     

Genome-wide association identifies a susceptibility locus for coronary artery disease in the Chinese Han population

Coronary artery disease (CAD) causes more than 700,000 deaths each year in China. Previous genome-wide association studies (GWAS) in populations of European ancestry identified several genetic loci for CAD, but no such study has yet been reported in the Chinese population. Here we report a three-stage GWAS in the Chinese Han population. We identified a new association between rs6903956 in a putative gene denoted as C6orf105 on chromosome 6p24.1 and CAD (P = 5.00 × 10?3, stage 2 validation; P = 3.00 × 10?3, P = 1.19 × 10?? and P = 4.00 × 10?3 in three independent stage 3 replication populations; P = 4.87 × 10?12, odds ratio = 1.51 in the combined population). The minor risk allele A of rs6903956 is associated with decreased C6orf105 mRNA expression. We report the first GWAS for CAD in the Chinese Han population and identify a SNP, rs6903956, in C6orf105 associated with susceptibility to CAD in this population.     

Genome-wide association study identifies FCGR2A as a susceptibility locus for Kawasaki disease

Kawasaki disease is a systemic vasculitis of unknown etiology, with clinical observations suggesting a substantial genetic contribution to disease susceptibility. We conducted a genome-wide association study and replication analysis in 2,173 individuals with Kawasaki disease and 9,383 controls from five independent sample collections. Two loci exceeded the formal threshold for genome-wide significance. The first locus is a functional polymorphism in the IgG receptor gene FCGR2A (encoding an H131R substitution) (rs1801274; P = 7.35 × 10(-11), odds ratio (OR) = 1.32), with the A allele (coding for histadine) conferring elevated disease risk. The second locus is at 19q13, (P = 2.51 × 10(-9), OR = 1.42 for the rs2233152 SNP near MIA and RAB4B; P = 1.68 × 10(-12), OR = 1.52 for rs28493229 in ITPKC), which confirms previous findings(1). The involvement of the FCGR2A locus may have implications for understanding immune activation in Kawasaki disease pathogenesis and the mechanism of response to intravenous immunoglobulin, the only proven therapy for this disease.     

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.     

Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants

We have genotyped 14,436 nonsynonymous SNPs (nsSNPs) and 897 major histocompatibility complex (MHC) tag SNPs from 1,000 independent cases of ankylosing spondylitis (AS), autoimmune thyroid disease (AITD), multiple sclerosis (MS) and breast cancer (BC). Comparing these data against a common control dataset derived from 1,500 randomly selected healthy British individuals, we report initial association and independent replication in a North American sample of two new loci related to ankylosing spondylitis, ARTS1 and IL23R, and confirmation of the previously reported association of AITD with TSHR and FCRL3. These findings, enabled in part by increased statistical power resulting from the expansion of the control reference group to include individuals from the other disease groups, highlight notable new possibilities for autoimmune regulation and suggest that IL23R may be a common susceptibility factor for the major 'seronegative' diseases.     

Mutational and functional analyses reveal that ST7 is a highly conserved tumor-suppressor gene on human chromosome 7q31

Loss of heterozygosity (LOH) of markers on human chromosome 7q31 is frequently encountered in a variety of human neoplasias, indicating the presence of a tumor-suppressor gene (TSG). By a combination of microcell-fusion and deletion-mapping studies, we previously established that this TSG resides within a critical region flanked by the genetic markers D7S522 and D7S677. Using a positional cloning strategy and aided by the availability of near-complete sequence of this genomic interval, we have identified a TSG within 7q31, named ST7 (for suppression of tumorigenicity 7; this same gene was recently reported in another context and called RAY1). ST7 is ubiquitously expressed in human tissues. Analysis of a series of cell lines derived from breast tumors and primary colon carcinomas revealed the presence of mutations in ST7. Introduction of the ST7 cDNA into the prostate-cancer-derived cell line PC3 had no effect on the in vitro proliferation of the cells, but abrogated their in vivo tumorigenicity. Our data indicate that ST7 is a TSG within chromosome 7q31 and may have an important role in the development of some types of human cancer.     

Common variants near MC4R are associated with fat mass, weight and risk of obesity

To identify common variants influencing body mass index (BMI), we analyzed genome-wide association data from 16,876 individuals of European descent. After previously reported variants in FTO, the strongest association signal (rs17782313, P = 2.9 x 10(-6)) mapped 188 kb downstream of MC4R (melanocortin-4 receptor), mutations of which are the leading cause of monogenic severe childhood-onset obesity. We confirmed the BMI association in 60,352 adults (per-allele effect = 0.05 Z-score units; P = 2.8 x 10(-15)) and 5,988 children aged 7-11 (0.13 Z-score units; P = 1.5 x 10(-8)). In case-control analyses (n = 10,583), the odds for severe childhood obesity reached 1.30 (P = 8.0 x 10(-11)). Furthermore, we observed overtransmission of the risk allele to obese offspring in 660 families (P (pedigree disequilibrium test average; PDT-avg) = 2.4 x 10(-4)). The SNP location and patterns of phenotypic associations are consistent with effects mediated through altered MC4R function. Our findings establish that common variants near MC4R influence fat mass, weight and obesity risk at the population level and reinforce the need for large-scale data integration to identify variants influencing continuous biomedical traits.     

Mutations in MRAP, encoding a new interacting partner of the ACTH receptor, cause familial glucocorticoid deficiency type 2

Familial glucocorticoid deficiency (FGD), or hereditary unresponsiveness to adrenocorticotropin (ACTH; OMIM 202200), is an autosomal recessive disorder resulting from resistance to the action of ACTH on the adrenal cortex, which stimulates glucocorticoid production. Affected individuals are deficient in cortisol and, if untreated, are likely to succumb to hypoglycemia or overwhelming infection in infancy or childhood. Mutations of the ACTH receptor (melanocortin 2 receptor, MC2R) account for approximately 25% of cases of FGD. FGD without mutations of MC2R is called FGD type 2. Using SNP array genotyping, we mapped a locus involved in FGD type 2 to chromosome 21q22.1. We identified mutations in a gene encoding a 19-kDa single-transmembrane domain protein, now known as melanocortin 2 receptor accessory protein (MRAP). We show that MRAP interacts with MC2R and may have a role in the trafficking of MC2R from the endoplasmic reticulum to the cell surface.     

Genome-wide mapping with biallelic markers in Arabidopsis thaliana.

Single-nucleotide polymorphisms, as well as small insertions and deletions (here referred to collectively as simple nucleotide polymorphisms, or SNPs), comprise the largest set of sequence variants in most organisms. Positional cloning based on SNPs may accelerate the identification of human disease traits and a range of biologically informative mutations. The recent application of high-density oligonucleotide arrays to allele identification has made it feasible to genotype thousands of biallelic SNPs in a single experiment. It has yet to be established, however, whether SNP detection using oligonucleotide arrays can be used to accelerate the mapping of traits in diploid genomes. The cruciferous weed Arabidopsis thaliana is an attractive model system for the construction and use of biallelic SNP maps. Although important biological processes ranging from fertilization and cell fate determination to disease resistance have been modelled in A. thaliana, identifying mutations in this organism has been impeded by the lack of a high-density genetic map consisting of easily genotyped DNA markers. We report here the construction of a biallelic genetic map in A. thaliana with a resolution of 3.5 cM and its use in mapping Eds16, a gene involved in the defence response to the fungal pathogen Erysiphe orontii. Mapping of this trait involved the high-throughput generation of meiotic maps of F2 individuals using high-density oligonucleotide probe array-based genotyping. We developed a software package called InterMap and used it to automatically delimit Eds16 to a 7-cM interval on chromosome 1. These results are the first demonstration of biallelic mapping in diploid genomes and establish means for generalizing SNP-based maps to virtually any genetic organism.