Differential confounding of rare and common variants in spatially structured populations

Well-powered genome-wide association studies, now made possible through advances in technology and large-scale collaborative projects, promise to characterize the contribution of rare variants to complex traits and disease. However, while population structure is a known confounder of association studies, it remains unknown whether methods developed to control stratification are equally effective for rare variants. Here, we demonstrate that rare variants can show a stratification that is systematically different from, and typically stronger than, common variants, and this is not necessarily corrected by existing methods. We show that the same process leads to inflation for load-based tests and can obscure signals at truly associated variants. Furthermore, we show that populations can display spatial structure in rare variants, even when Wright's fixation index F(ST) is low, but that allele frequency-dependent metrics of allele sharing can reveal localized stratification. These results underscore the importance of collecting and integrating spatial information in the genetic analysis of complex traits.     

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.     

Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs

Accurate and complete measurement of single nucleotide (SNP) and copy number (CNV) variants, both common and rare, will be required to understand the role of genetic variation in disease. We present Birdsuite, a four-stage analytical framework instantiated in software for deriving integrated and mutually consistent copy number and SNP genotypes. The method sequentially assigns copy number across regions of common copy number polymorphisms (CNPs), calls genotypes of SNPs, identifies rare CNVs via a hidden Markov model (HMM), and generates an integrated sequence and copy number genotype at every locus (for example, including genotypes such as A-null, AAB and BBB in addition to AA, AB and BB calls). Such genotypes more accurately depict the underlying sequence of each individual, reducing the rate of apparent mendelian inconsistencies. The Birdsuite software is applied here to data from the Affymetrix SNP 6.0 array. Additionally, we describe a method, implemented in PLINK, to utilize these combined SNP and CNV genotypes for association testing with a phenotype.     

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.     

Common and rare variants in multifactorial susceptibility to common diseases

Here, we give a historical overview of the search for genetic variants that influence the susceptibility of an individual to a chronic disease, from RA Fisher's seminal work to the current excitement of whole-genome association studies (WGAS). We then discuss the concepts behind the identification of common variants as disease causal factors and contrast them to the basic ideas that underlie the rare variant hypothesis. The identification of rare variants involves the careful selection of candidate genes to examine, the availability of highly efficient resequencing techniques and the appropriate assessment of the functional consequences of the implicated variant. We believe that this strategy can be successfully applied at present in order to unravel the contribution of rare variants to the multifactorial inheritance of common diseases, which could lead to the implementation of much needed preventative screening schemes.     

Comprehensive analysis of the gene encoding filaggrin uncovers prevalent and rare mutations in ichthyosis vulgaris and atopic eczema

We recently reported two common filaggrin (FLG) null mutations that cause ichthyosis vulgaris and predispose to eczema and secondary allergic diseases. We show here that these common European mutations are ancestral variants carried on conserved haplotypes. To facilitate comprehensive analysis of other populations, we report a strategy for full sequencing of this large, highly repetitive gene, and we describe 15 variants, including seven that are prevalent. All the variants are either nonsense or frameshift mutations that, in representative cases, resulted in loss of filaggrin production in the epidermis. In an Irish case-control study, the five most common European mutations showed a strong association with moderate-to-severe childhood eczema (chi2 test: P = 2.12 x 10(-51); Fisher's exact test: heterozygote odds ratio (OR) = 7.44 (95% confidence interval (c.i.) = 4.9-11.3), and homozygote OR = 151 (95% c.i. = 20-1,136)). We found three additional rare null mutations in this case series, suggesting that the genetic architecture of filaggrin-related atopic dermatitis consists of both prevalent and rare risk alleles.     

Truncating mutations in the last exon of NOTCH2 cause a rare skeletal disorder with osteoporosis

Hajdu-Cheney syndrome is a rare autosomal dominant skeletal disorder with facial anomalies, osteoporosis and acro-osteolysis. We sequenced the exomes of six unrelated individuals with this syndrome and identified heterozygous nonsense and frameshift mutations in NOTCH2 in five of them. All mutations cluster to the last coding exon of the gene, suggesting that the mutant mRNA products escape nonsense-mediated decay and that the resulting truncated NOTCH2 proteins act in a gain-of-function manner.     

A rare variant in MYH6 is associated with high risk of sick sinus syndrome

Through complementary application of SNP genotyping, whole-genome sequencing and imputation in 38,384 Icelanders, we have discovered a previously unidentified sick sinus syndrome susceptibility gene, MYH6, encoding the alpha heavy chain subunit of cardiac myosin. A missense variant in this gene, c.2161C>T, results in the conceptual amino acid substitution p.Arg721Trp, has an allelic frequency of 0.38% in Icelanders and associates with sick sinus syndrome with an odds ratio = 12.53 and P = 1.5 × 10?2?. We show that the lifetime risk of being diagnosed with sick sinus syndrome is around 6% for non-carriers of c.2161C>T but is approximately 50% for carriers of the c.2161C>T variant.     

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.     

MYH9 is associated with nondiabetic end-stage renal disease in African Americans

As end-stage renal disease (ESRD) has a four times higher incidence in African Americans compared to European Americans, we hypothesized that susceptibility alleles for ESRD have a higher frequency in the West African than the European gene pool. We carried out a genome-wide admixture scan in 1,372 ESRD cases and 806 controls and found a highly significant association between excess African ancestry and nondiabetic ESRD (lod score = 5.70) but not diabetic ESRD (lod = 0.47) on chromosome 22q12. Each copy of the European ancestral allele conferred a relative risk of 0.50 (95% CI = 0.39-0.63) compared to African ancestry. Multiple common SNPs (allele frequencies ranging from 0.2 to 0.6) in the gene encoding nonmuscle myosin heavy chain type II isoform A (MYH9) were associated with two to four times greater risk of nondiabetic ESRD and accounted for a large proportion of the excess risk of ESRD observed in African compared to European Americans.     

Gene polymorphism in Netherton and common atopic disease.

Atopic dermatitis (AD) and asthma are characterized by IgE-mediated atopic (allergic) responses to common proteins (allergens), many of which are proteinases. Loci influencing atopy have been localized to a number of chromosomal regions, including the chromosome 5q31 cytokine cluster. Netherton disease is a rare recessive skin disorder in which atopy is a universal accompaniment. The gene underlying Netherton disease (SPINK5) encodes a 15-domain serine proteinase inhibitor (LEKTI) which is expressed in epithelial and mucosal surfaces and in the thymus. We have identified six coding polymorphisms in SPINK5 (Table 1) and found that a Glu420-->Lys variant shows significant association with atopy and AD in two independent panels of families. Our results implicate a previously unrecognized pathway for the development of common allergic illnesses.     

Mutations in KERA, encoding keratocan, cause cornea plana

Specialized collagens and small leucine-rich proteoglycans (SLRPs) interact to produce the transparent corneal structure. In cornea plana, the forward convex curvature is flattened, leading to a decrease in refraction. A more severe, recessively inherited form (CNA2; MIM 217300) and a milder, dominantly inherited form (CNA1; MIM 121400) exist. CNA2 is a rare disorder with a worldwide distribution, but a high prevalence in the Finnish population. The gene mutated in CNA2 was assigned by linkage analysis to 12q (refs 4, 5), where there is a cluster of several SLRP genes. We cloned two additional SLRP genes highly expressed in cornea: KERA (encoding keratocan) in 12q and OGN (encoding osteoglycin) in 9q. Here we report mutations in KERA in 47 CNA2 patients: 46 Finnish patients are homozygous for a founder missense mutation, leading to the substitution of a highly conserved amino acid; and one American patient is homozygous for a mutation leading to a premature stop codon that truncates the KERA protein. Our data establish that mutations in KERA cause CNA2. CNA1 patients had no mutations in these proteoglycan genes.     

A degradation-sensitive anionic trypsinogen (PRSS2) variant protects against chronic pancreatitis

Chronic pancreatitis is a common inflammatory disease of the pancreas. Mutations in the genes encoding cationic trypsinogen (PRSS1) and the pancreatic secretory trypsin inhibitor (SPINK1) are associated with chronic pancreatitis. Because increased proteolytic activity owing to mutated PRSS1 enhances the risk for chronic pancreatitis, mutations in the gene encoding anionic trypsinogen (PRSS2) may also predispose to disease. Here we analyzed PRSS2 in individuals with chronic pancreatitis and controls and found, to our surprise, that a variant of codon 191 (G191R) is overrepresented in control subjects: G191R was present in 220/6,459 (3.4%) controls but in only 32/2,466 (1.3%) affected individuals (odds ratio 0.37; P = 1.1 x 10(-8)). Upon activation by enterokinase or trypsin, purified recombinant G191R protein showed a complete loss of trypsin activity owing to the introduction of a new tryptic cleavage site that renders the enzyme hypersensitive to autocatalytic proteolysis. In conclusion, the G191R variant of PRSS2 mitigates intrapancreatic trypsin activity and thereby protects against chronic pancreatitis.     

A single EFEMP1 mutation associated with both Malattia Leventinese and Doyne honeycomb retinal dystrophy.

Malattia Leventinese (ML) and Doyne honeycomb retinal dystrophy (DHRD) refer to two autosomal dominant diseases characterized by yellow-white deposits known as drusen that accumulate beneath the retinal pigment epithelium (RPE). Both loci were mapped to chromosome 2p16-21 (refs 5,6) and this genetic interval has been subsequently narrowed. The importance of these diseases is due in large part to their close phenotypic similarity to age-related macular degeneration (AMD), a disorder with a strong genetic component that accounts for approximately 50% of registered blindness in the Western world. Just as in ML and DHRD, the early hallmark of AMD is the presence of drusen. Here we use a combination of positional and candidate gene methods to identify a single non-conservative mutation (Arg345Trp) in the gene EFEMP1 (for EGF-containing fibrillin-like extracellular matrix protein 1) in all families studied. This change was not present in 477 control individuals or in 494 patients with age-related macular degeneration. Identification of this mutation may aid in the development of an animal model for drusen, as well as in the identification of other genes involved in human macular degeneration.     

Extremely low-coverage sequencing and imputation increases power for genome-wide association studies

Genome-wide association studies (GWAS) have proven to be a powerful method to identify common genetic variants contributing to susceptibility to common diseases. Here, we show that extremely low-coverage sequencing (0.1-0.5×) captures almost as much of the common (>5%) and low-frequency (1-5%) variation across the genome as SNP arrays. As an empirical demonstration, we show that genome-wide SNP genotypes can be inferred at a mean r(2) of 0.71 using off-target data (0.24× average coverage) in a whole-exome study of 909 samples. Using both simulated and real exome-sequencing data sets, we show that association statistics obtained using extremely low-coverage sequencing data attain similar P values at known associated variants as data from genotyping arrays, without an excess of false positives. Within the context of reductions in sample preparation and sequencing costs, funds invested in extremely low-coverage sequencing can yield several times the effective sample size of GWAS based on SNP array data and a commensurate increase in statistical power.