Polymorphism for a 1.6-Mb deletion of the human Y chromosome persists through balance between recurrent mutation and haploid selection

Many human Y-chromosomal deletions are thought to severely impair reproductive fitness, which precludes their transmission to the next generation and thus ensures their rarity in the population. Here we report a 1.6-Mb deletion that persists over generations and is sufficiently common to be considered a polymorphism. We hypothesized that this deletion might affect spermatogenesis because it removes almost half of the Y chromosome's AZFc region, a gene-rich segment that is critical for sperm production. An association study established that this deletion, called gr/gr, is a significant risk factor for spermatogenic failure. The gr/gr deletion has far lower penetrance with respect to spermatogenic failure than previously characterized Y-chromosomal deletions; it is often transmitted from father to son. By studying the distribution of gr/gr-deleted chromosomes across the branches of the Y chromosome's genealogical tree, we determined that this deletion arose independently at least 14 times in human history. We suggest that the existence of this deletion as a polymorphism reflects a balance between haploid selection, which culls gr/gr-deleted Y chromosomes from the population, and homologous recombination, which continues to generate new gr/gr deletions.     

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.     

A common variant of HMGA2 is associated with adult and childhood height in the general population

Human height is a classic, highly heritable quantitative trait. To begin to identify genetic variants influencing height, we examined genome-wide association data from 4,921 individuals. Common variants in the HMGA2 oncogene, exemplified by rs1042725, were associated with height (P = 4 x 10(-8)). HMGA2 is also a strong biological candidate for height, as rare, severe mutations in this gene alter body size in mice and humans, so we tested rs1042725 in additional samples. We confirmed the association in 19,064 adults from four further studies (P = 3 x 10(-11), overall P = 4 x 10(-16), including the genome-wide association data). We also observed the association in children (P = 1 x 10(-6), N = 6,827) and a tall/short case-control study (P = 4 x 10(-6), N = 3,207). We estimate that rs1042725 explains approximately 0.3% of population variation in height (approximately 0.4 cm increased adult height per C allele). There are few examples of common genetic variants reproducibly associated with human quantitativetraits; these results represent, to our knowledge, the first consistently replicated association with adult and childhood height.     

Identification of ten loci associated with height highlights new biological pathways in human growth

Height is a classic polygenic trait, reflecting the combined influence of multiple as-yet-undiscovered genetic factors. We carried out a meta-analysis of genome-wide association study data of height from 15,821 individuals at 2.2 million SNPs, and followed up the strongest findings in >10,000 subjects. Ten newly identified and two previously reported loci were strongly associated with variation in height (P values from 4 x 10(-7) to 8 x 10(-22)). Together, these 12 loci account for approximately 2% of the population variation in height. Individuals with < or =8 height-increasing alleles and > or =16 height-increasing alleles differ in height by approximately 3.5 cm. The newly identified loci, along with several additional loci with strongly suggestive associations, encompass both strong biological candidates and unexpected genes, and highlight several pathways (let-7 targets, chromatin remodeling proteins and Hedgehog signaling) as important regulators of human stature. These results expand the picture of the biological regulation of human height and of the genetic architecture of this classical complex trait.     

FCGR3B copy number variation is associated with susceptibility to systemic, but not organ-specific, autoimmunity

Naturally occurring variation in gene copy number is increasingly recognized as a heritable source of susceptibility to genetically complex diseases. Here we report strong association between FCGR3B copy number and risk of systemic lupus erythematosus (P = 2.7 x 10(-8)), microscopic polyangiitis (P = 2.9 x 10(-4)) and Wegener's granulomatosis in two independent cohorts from the UK (P = 3 x 10(-3)) and France (P = 1.1 x 10(-4)). We did not observe this association in the organ-specific Graves' disease or Addison's disease. Our findings suggest that low FCGR3B copy number, and in particular complete FCGR3B deficiency, has a key role in the development of systemic autoimmunity.     

Sequence variants in the autophagy gene IRGM and multiple other replicating loci contribute to Crohn's disease susceptibility

A genome-wide association scan in individuals with Crohn's disease by the Wellcome Trust Case Control Consortium detected strong association at four novel loci. We tested 37 SNPs from these and other loci for association in an independent case-control sample. We obtained replication for the autophagy-inducing IRGM gene on chromosome 5q33.1 (replication P = 6.6 x 10(-4), combined P = 2.1 x 10(-10)) and for nine other loci, including NKX2-3, PTPN2 and gene deserts on chromosomes 1q and 5p13.     

Discovery and genotyping of genome structural polymorphism by sequencing on a population scale

Accurate and complete analysis of genome variation in large populations will be required to understand the role of genome variation in complex disease. We present an analytical framework for characterizing genome deletion polymorphism in populations using sequence data that are distributed across hundreds or thousands of genomes. Our approach uses population-level concepts to reinterpret the technical features of sequence data that often reflect structural variation. In the 1000 Genomes Project pilot, this approach identified deletion polymorphism across 168 genomes (sequenced at 4 × average coverage) with sensitivity and specificity unmatched by other algorithms. We also describe a way to determine the allelic state or genotype of each deletion polymorphism in each genome; the 1000 Genomes Project used this approach to type 13,826 deletion polymorphisms (48-995,664 bp) at high accuracy in populations. These methods offer a way to relate genome structural polymorphism to complex disease in populations.     

Loci on 20q13 and 21q22 are associated with pediatric-onset inflammatory bowel disease

Inflammatory bowel disease (IBD) is a common inflammatory disorder with complex etiology that involves both genetic and environmental triggers, including but not limited to defects in bacterial clearance, defective mucosal barrier and persistent dysregulation of the immune response to commensal intestinal bacteria. IBD is characterized by two distinct phenotypes: Crohn's disease (CD) and ulcerative colitis (UC). Previously reported GWA studies have identified genetic variation accounting for a small portion of the overall genetic susceptibility to CD and an even smaller contribution to UC pathogenesis. We hypothesized that stratification of IBD by age of onset might identify additional genes associated with IBD. To that end, we carried out a GWA analysis in a cohort of 1,011 individuals with pediatric-onset IBD and 4,250 matched controls. We identified and replicated significantly associated, previously unreported loci on chromosomes 20q13 (rs2315008[T] and rs4809330[A]; P = 6.30 x 10(-8) and 6.95 x 10(-8), respectively; odds ratio (OR) = 0.74 for both) and 21q22 (rs2836878[A]; P = 6.01 x 10(-8); OR = 0.73), located close to the TNFRSF6B and PSMG1 genes, respectively.     

Evidence that positive selection drives Y-chromosome degeneration in Drosophila miranda

Why does the Y chromosome harbor so few functional loci? Evolutionary theory predicts that Y chromosomes degenerate because they lack genetic recombination. Both positive and negative selection models have been invoked to explain this degeneration, as both can result in the recurrent fixation of linked deleterious mutations on a nonrecombining Y chromosome. To distinguish between these models, I investigated patterns of nucleotide variability along 37 kb of the recently formed neo-Y chromosome in Drosophila miranda. Levels of nucleotide variability on this chromosome are 30 times lower than in highly recombining portions of the genome. Both positive and negative selection models can result in reduced variability levels, but their effects on the frequency spectrum of mutations differ. Using coalescent simulations, I show that the patterns of nucleotide variability on the neo-Y chromosome are unlikely under deleterious mutation models (including background selection and Muller's ratchet) but are expected under recent positive selection. These results implicate positive selection as an important force driving the degeneration of Y chromosomes; adaptation at a few loci, possibly increasing male fitness, occurs at the cost of most other genes on this chromosome.     

Multiple regions within 8q24 independently affect risk for prostate cancer

After the recent discovery that common genetic variation in 8q24 influences inherited risk of prostate cancer, we genotyped 2,973 SNPs in up to 7,518 men with and without prostate cancer from five populations. We identified seven risk variants, five of them previously undescribed, spanning 430 kb and each independently predicting risk for prostate cancer (P = 7.9 x 10(-19) for the strongest association, and P < 1.5 x 10(-4) for five of the variants, after controlling for each of the others). The variants define common genotypes that span a more than fivefold range of susceptibility to cancer in some populations. None of the prostate cancer risk variants aligns to a known gene or alters the coding sequence of an encoded protein.     

Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci

Systemic lupus erythematosus (SLE) is a common systemic autoimmune disease with complex etiology but strong clustering in families (lambda(S) = approximately 30). We performed a genome-wide association scan using 317,501 SNPs in 720 women of European ancestry with SLE and in 2,337 controls, and we genotyped consistently associated SNPs in two additional independent sample sets totaling 1,846 affected women and 1,825 controls. Aside from the expected strong association between SLE and the HLA region on chromosome 6p21 and the previously confirmed non-HLA locus IRF5 on chromosome 7q32, we found evidence of association with replication (1.1 x 10(-7) < P(overall) < 1.6 x 10(-23); odds ratio = 0.82-1.62) in four regions: 16p11.2 (ITGAM), 11p15.5 (KIAA1542), 3p14.3 (PXK) and 1q25.1 (rs10798269). We also found evidence for association (P < 1 x 10(-5)) at FCGR2A, PTPN22 and STAT4, regions previously associated with SLE and other autoimmune diseases, as well as at > or =9 other loci (P < 2 x 10(-7)). Our results show that numerous genes, some with known immune-related functions, predispose to SLE.     

A genome-wide association study identifies six susceptibility loci for chronic lymphocytic leukemia

We conducted a genome-wide association study of 299,983 tagging SNPs for chronic lymphocytic leukemia (CLL) and performed validation in two additional series totaling 1,529 cases and 3,115 controls. We identified six previously unreported CLL risk loci at 2q13 (rs17483466; P = 2.36 x 10(-10)), 2q37.1 (rs13397985, SP140; P = 5.40 x 10(-10)), 6p25.3 (rs872071, IRF4; P = 1.91 x 10(-20)), 11q24.1 (rs735665; P = 3.78 x 10(-12)), 15q23 (rs7176508; P = 4.54 x 10(-12)) and 19q13.32 (rs11083846, PRKD2; P = 3.96 x 10(-9)). These data provide the first evidence for the existence of common, low-penetrance susceptibility to a hematological malignancy and new insights into disease causation in CLL.     

A genome-wide association study for celiac disease identifies risk variants in the region harboring IL2 and IL21

We tested 310,605 SNPs for association in 778 individuals with celiac disease and 1,422 controls. Outside the HLA region, the most significant finding (rs13119723; P = 2.0 x 10(-7)) was in the KIAA1109-TENR-IL2-IL21 linkage disequilibrium block. We independently confirmed association in two further collections (strongest association at rs6822844, 24 kb 5' of IL21; meta-analysis P = 1.3 x 10(-14), odds ratio = 0.63), suggesting that genetic variation in this region predisposes to celiac disease.     

Mutations in the gene encoding the 3'-5' DNA exonuclease TREX1 are associated with systemic lupus erythematosus

TREX1 acts in concert with the SET complex in granzyme A-mediated apoptosis, and mutations in TREX1 cause Aicardi-Goutières syndrome and familial chilblain lupus. Here, we report monoallelic frameshift or missense mutations and one 3' UTR variant of TREX1 present in 9/417 individuals with systemic lupus erythematosus but absent in 1,712 controls (P = 4.1 x 10(-7)). We demonstrate that two mutant TREX1 alleles alter subcellular targeting. Our findings implicate TREX1 in the pathogenesis of SLE.     

PALB2, which encodes a BRCA2-interacting protein, is a breast cancer susceptibility gene

PALB2 interacts with BRCA2, and biallelic mutations in PALB2 (also known as FANCN), similar to biallelic BRCA2 mutations, cause Fanconi anemia. We identified monoallelic truncating PALB2 mutations in 10/923 individuals with familial breast cancer compared with 0/1,084 controls (P = 0.0004) and show that such mutations confer a 2.3-fold higher risk of breast cancer (95% confidence interval (c.i.) = 1.4-3.9, P = 0.0025). The results show that PALB2 is a breast cancer susceptibility gene and further demonstrate the close relationship of the Fanconi anemia-DNA repair pathway and breast cancer predisposition.     

Rheumatoid arthritis susceptibility loci at chromosomes 10p15, 12q13 and 22q13

The WTCCC study identified 49 SNPs putatively associated with rheumatoid arthritis at P = 1 x 10(-4) - 1 x 10(-5) (tier 3). Here we show that three of these SNPs, mapping to chromosome 10p15 (rs4750316), 12q13 (rs1678542) and 22q13 (rs3218253), are also associated (trend P = 4 x 10(-5), P = 4 x 10(-4) and P = 4 x 10(-4), respectively) in a validation study of 4,106 individuals with rheumatoid arthritis and an expanded reference group of 11,238 subjects, confirming them as true susceptibility loci in individuals of European ancestry.     

Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus

We carried out a multistage genome-wide association study of type 2 diabetes mellitus in Japanese individuals, with a total of 1,612 cases and 1,424 controls and 100,000 SNPs. The most significant association was obtained with SNPs in KCNQ1, and dense mapping within the gene revealed that rs2237892 in intron 15 showed the lowest Pvalue (6.7 x 10(-13), odds ratio (OR) = 1.49). The association of KCNQ1 with type 2 diabetes was replicated in populations of Korean, Chinese and European ancestry as well as in two independent Japanese populations, and meta-analysis with a total of 19,930 individuals (9,569 cases and 10,361 controls) yielded a P value of 1.7 x 10(-42) (OR = 1.40; 95% CI = 1.34-1.47) for rs2237892. Among control subjects, the risk allele of this polymorphism was associated with impairment of insulin secretion according to the homeostasis model assessment of beta-cell function or the corrected insulin response. Our data thus implicate KCNQ1 as a diabetes susceptibility gene in groups of different ancestries.     

Identification of loci associated with schizophrenia by genome-wide association and follow-up

We carried out a genome-wide association study of schizophrenia (479 cases, 2,937 controls) and tested loci with P < 10(-5) in up to 16,726 additional subjects. Of 12 loci followed up, 3 had strong independent support (P < 5 x 10(-4)), and the overall pattern of replication was unlikely to occur by chance (P = 9 x 10(-8)). Meta-analysis provided strongest evidence for association around ZNF804A (P = 1.61 x 10(-7)) and this strengthened when the affected phenotype included bipolar disorder (P = 9.96 x 10(-9)).     

Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24

Using a multistage genetic association approach comprising 7,480 affected individuals and 7,779 controls, we identified markers in chromosomal region 8q24 associated with colorectal cancer. In stage 1, we genotyped 99,632 SNPs in 1,257 affected individuals and 1,336 controls from Ontario. In stages 2-4, we performed serial replication studies using 4,024 affected individuals and 4,042 controls from Seattle, Newfoundland and Scotland. We identified one locus on chromosome 8q24 and another on 9p24 having combined odds ratios (OR) for stages 1-4 of 1.18 (trend; P = 1.41 x 10(-8)) and 1.14 (trend; P = 1.32 x 10(-5)), respectively. Additional analyses in 2,199 affected individuals and 2,401 controls from France and Europe supported the association at the 8q24 locus (OR = 1.16, trend; 95% confidence interval (c.i.): 1.07-1.26; P = 5.05 x 10(-4)). A summary across all seven studies at the 8q24 locus was highly significant (OR = 1.17, c.i.: 1.12-1.23; P = 3.16 x 10(-11)). This locus has also been implicated in prostate cancer.