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.     

Osteoclast-poor human osteopetrosis due to mutations in the gene encoding RANKL

Autosomal recessive osteopetrosis is usually associated with normal or elevated numbers of nonfunctional osteoclasts. Here we report mutations in the gene encoding RANKL (receptor activator of nuclear factor-KB ligand) in six individuals with autosomal recessive osteopetrosis whose bone biopsy specimens lacked osteoclasts. These individuals did not show any obvious defects in immunological parameters and could not be cured by hematopoietic stem cell transplantation; however, exogenous RANKL induced formation of functional osteoclasts from their monocytes, suggesting that they could, theoretically, benefit from exogenous RANKL administration.     

Assessing the impact of population stratification on genetic association studies

Population stratification refers to differences in allele frequencies between cases and controls due to systematic differences in ancestry rather than association of genes with disease. It has been proposed that false positive associations due to stratification can be controlled by genotyping a few dozen unlinked genetic markers. To assess stratification empirically, we analyzed data from 11 case-control and case-cohort association studies. We did not detect statistically significant evidence for stratification but did observe that assessments based on a few dozen markers lack power to rule out moderate levels of stratification that could cause false positive associations in studies designed to detect modest genetic risk factors. After increasing the number of markers and samples in a case-cohort study (the design most immune to stratification), we found that stratification was in fact present. Our results suggest that modest amounts of stratification can exist even in well designed studies.     

Inherited variants of MYH associated with somatic G:C-->T:A mutations in colorectal tumors

Inherited defects of base excision repair have not been associated with any human genetic disorder, although mutations of the genes mutM and mutY, which function in Escherichia coli base excision repair, lead to increased transversions of G:C to T:A. We have studied family N, which is affected with multiple colorectal adenomas and carcinoma but lacks an inherited mutation of the adenomatous polyposis coli gene (APC) that is associated with familial adenomatous polyposis. Here we show that 11 tumors from 3 affected siblings contain 18 somatic inactivating mutations of APC and that 15 of these mutations are G:C-->A transversions--a significantly greater proportion than is found in sporadic tumors or in tumors associated with familial adenomatous polyposis. Analysis of the human homolog of mutY, MYH, showed that the siblings were compound heterozygotes for the nonconservative missense variants Tyr165Cys and Gly382Asp. These mutations affect residues that are conserved in mutY of E. coli (Tyr82 and Gly253). Tyrosine 82 is located in the pseudo-helix-hairpin-helix (HhH) motif and is predicted to function in mismatch specificity. Assays of adenine glycosylase activity of the Tyr82Cys and Gly253Asp mutant proteins with 8-oxoG:A and G:A substrates show that their activity is reduced significantly. Our findings link the inherited variants in MYH to the pattern of somatic APC mutation in family N and implicate defective base excision repair in predisposition to tumors in humans.     

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.     

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.     

The creation of 'consistent areas through time' (CATTs) in Scotland, 1981-2001

Small area analysis of sociodemographic change over time is often impeded by the changing configuration of census zones for each decennial UK census. Approaches for matching zone configurations exist, but these typically require population estimation techniques, which inevitably involve some error. Because of the postcode-based geographical zoning system used in Scottish censuses, it is possible to create a local-area geography which is consistent for recent censuses. Here we present a methodology used to create three sets of consistent areas through time (CATTs) that can be used for analysing comparable small area data output from the 1981, 1991 and 2001 Censuses.     

Detection in vivo of protein-DNA interactions within the lac operon of Escherichia coli

Studies of the sequence-specific binding of proteins to DNA have so far relied on in vitro experiments using cloned restriction fragments containing the relevant DNA sequences. We have applied the genomic sequencing technique of Church and Gilbert to show that the interactions observed in vitro occur in vivo. We use this approach to study the binding of regulatory proteins to the lac operon in vivo and detect changes in the reactivity (inhibition or enhancement) of guanines to methylation by dimethyl sulphate caused by the proximity of proteins to the N-7 atom of these guanines. We can detect the simultaneous binding of the catobolite gene activator protein (CAP) and the Lac repressor to their specific recognition sequences, and following induction of the lac operon we observe effects that are related to RNA polymerase binding or RNA elongation. We have successfully used oligonucleotide probes as short as 17 bases to display genomic sequence.