Automating sequence-based detection and genotyping of SNPs from diploid samples

The detection of sequence variation, for which DNA sequencing has emerged as the most sensitive and automated approach, forms the basis of all genetic analysis. Here we describe and illustrate an algorithm that accurately detects and genotypes SNPs from fluorescence-based sequence data. Because the algorithm focuses particularly on detecting SNPs through the identification of heterozygous individuals, it is especially well suited to the detection of SNPs in diploid samples obtained after DNA amplification. It is substantially more accurate than existing approaches and, notably, provides a useful quantitative measure of its confidence in each potential SNP detected and in each genotype called. Calls assigned the highest confidence are sufficiently reliable to remove the need for manual review in several contexts. For example, for sequence data from 47-90 individuals sequenced on both the forward and reverse strands, the highest-confidence calls from our algorithm detected 93% of all SNPs and 100% of high-frequency SNPs, with no false positive SNPs identified and 99.9% genotyping accuracy. This algorithm is implemented in a software package, PolyPhred version 5.0, which is freely available for academic use.     

Genome-wide genetic association of complex traits in heterogeneous stock mice

Difficulties in fine-mapping quantitative trait loci (QTLs) are a major impediment to progress in the molecular dissection of complex traits in mice. Here we show that genome-wide high-resolution mapping of multiple phenotypes can be achieved using a stock of genetically heterogeneous mice. We developed a conservative and robust bootstrap analysis to map 843 QTLs with an average 95% confidence interval of 2.8 Mb. The QTLs contribute to variation in 97 traits, including models of human disease (asthma, type 2 diabetes mellitus, obesity and anxiety) as well as immunological, biochemical and hematological phenotypes. The genetic architecture of almost all phenotypes was complex, with many loci each contributing a small proportion to the total variance. Our data set, freely available at http://gscan.well.ox.ac.uk, provides an entry point to the functional characterization of genes involved in many complex traits.     

Irf6 is a key determinant of the keratinocyte proliferation-differentiation switch

The epidermis is a highly organized structure, the integrity of which is central to the protection of an organism. Development and subsequent maintenance of this tissue depends critically on the intricate balance between proliferation and differentiation of a resident stem cell population; however, the signals controlling the proliferation-differentiation switch in vivo remain elusive. Here, we show that mice carrying a homozygous missense mutation in interferon regulatory factor 6 (Irf6), the homolog of the gene mutated in the human congenital disorders Van der Woude syndrome and popliteal pterygium syndrome, have a hyperproliferative epidermis that fails to undergo terminal differentiation, resulting in soft tissue fusions. We further demonstrate that mice that are compound heterozygotes for mutations in Irf6 and the gene encoding the cell cycle regulator protein stratifin (Sfn; also known as 14-3-3sigma) show similar defects of keratinizing epithelia. Our results indicate that Irf6 is a key determinant of the keratinocyte proliferation-differentiation switch and that Irf6 and Sfn interact genetically in this process.     

Conserved requirement for a plant host cell protein in powdery mildew pathogenesis

In the fungal phylum Ascomycota, the ability to cause disease in plants and animals has been gained and lost repeatedly during phylogenesis. In monocotyledonous barley, loss-of-function mlo alleles result in effective immunity against the Ascomycete Blumeria graminis f. sp. hordei, the causal agent of powdery mildew disease. However, mlo-based disease resistance has been considered a barley-specific phenomenon to date. Here, we demonstrate a conserved requirement for MLO proteins in powdery mildew pathogenesis in the dicotyledonous plant species Arabidopsis thaliana. Epistasis analysis showed that mlo resistance in A. thaliana does not involve the signaling molecules ethylene, jasmonic acid or salicylic acid, but requires a syntaxin, glycosyl hydrolase and ABC transporter. These findings imply that a common host cell entry mechanism of powdery mildew fungi evolved once and at least 200 million years ago, suggesting that within the Erysiphales (powdery mildews) the ability to cause disease has been a stable trait throughout phylogenesis.     

Constitutional 11p15 abnormalities, including heritable imprinting center mutations, cause nonsyndromic Wilms tumor

Constitutional abnormalities at the imprinted 11p15 growth regulatory region cause syndromes characterized by disordered growth, some of which include a risk of Wilms tumor. We explored their possible contribution to nonsyndromic Wilms tumor and identified constitutional 11p15 abnormalities in genomic lymphocyte DNA from 13 of 437 individuals (3%) with sporadic Wilms tumor without features of growth disorders, including 12% of bilateral cases (P = 0.001) and in one familial Wilms tumor pedigree. No abnormality was detected in 220 controls (P = 0.006). Abnormalities identified included H19 DMR epimutations, uniparental disomy 11p15 and H19 DMR imprinting center mutations (one microinsertion and one microdeletion), thus identifying microinsertion as a new class of imprinting center mutation. Our data identify constitutional 11p15 defects as one of the most common known causes of Wilms tumor, provide mechanistic insights into imprinting disruption and reveal clinically important epigenotype-phenotype associations. The impact on clinical management dictates that constitutional 11p15 analysis should be considered in all individuals with Wilms tumor.     

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.     

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.     

SLC2A9 is a newly identified urate transporter influencing serum urate concentration, urate excretion and gout

Uric acid is the end product of purine metabolism in humans and great apes, which have lost hepatic uricase activity, leading to uniquely high serum uric acid concentrations (200-500 microM) compared with other mammals (3-120 microM). About 70% of daily urate disposal occurs via the kidneys, and in 5-25% of the human population, impaired renal excretion leads to hyperuricemia. About 10% of people with hyperuricemia develop gout, an inflammatory arthritis that results from deposition of monosodium urate crystals in the joint. We have identified genetic variants within a transporter gene, SLC2A9, that explain 1.7-5.3% of the variance in serum uric acid concentrations, following a genome-wide association scan in a Croatian population sample. SLC2A9 variants were also associated with low fractional excretion of uric acid and/or gout in UK, Croatian and German population samples. SLC2A9 is a known fructose transporter, and we now show that it has strong uric acid transport activity in Xenopus laevis oocytes.