Rare genetic variants contribute to complex disease risk; however, the abundance of rare variants in human populations remains unknown. We explored this spectrum of variation by sequencing 202 genes encoding drug targets in 14,002 individuals. We find rare variants are abundant (one every 17 bases) and geographically localized, such that even with large sample sizes, rare variant catalogs will be largely incomplete. We used the observed patterns of variation to estimate population growth parameters, the proportion of variants in a given frequency class that are putatively deleterious, and mutation rates for each gene. Overall we conclude that, due to rapid population growth and weak purifying selection, human populations harbor an abundance of rare variants, many of which are deleterious and have relevance to understanding disease risk.

Fasting glucose and insulin are intermediate traits for type 2 diabetes. Here we explore the role of coding variation on these traits by analysis of variants on the HumanExome BeadChip in 60,564 non-diabetic individuals and in 16,491 T2D cases and 81,877 controls. We identify a novel association of a low-frequency nonsynonymous SNV in GLP1R (A316T; rs10305492; MAF=1.4%) with lower FG (β=-0.09±0.01 mmol L−1, p=3.4×10−12), T2D risk (OR[95%CI]=0.86[0.76-0.96], p=0.010), early insulin secretion (β=-0.07±0.035 pmolinsulin mmolglucose−1, p=0.048), but higher 2-h glucose (β=0.16±0.05 mmol L−1, p=4.3×10−4). We identify a gene-based association with FG at G6PC2 (pSKAT=6.8×10−6) driven by four rare protein-coding SNVs (H177Y, Y207S, R283X and S324P). We identify rs651007 (MAF=20%) in the first intron of ABO at the putative promoter of an antisense lncRNA, associating with higher FG (β=0.02±0.004 mmol L−1, p=1.3×10−8). Our approach identifies novel coding variant associations and extends the allelic spectrum of variation underlying diabetes-related quantitative traits and T2D susceptibility.

These data provide further support for SERPINE2 as a COPD susceptibility gene.

Regulatory authorities have indicated that new drugs to treat type 2 diabetes (T2D) should not be associated with an unacceptable increase in cardiovascular risk. Human genetics may be able to inform development of antidiabetic therapies by predicting cardiovascular and other health endpoints. We therefore investigated the association of variants in 6 genes that encode drug targets for obesity or T2D with a range of metabolic traits in up to 11,806 individuals by targeted exome sequencing, and follow-up in 39,979 individuals by targeted genotyping, with additional in silico follow up in consortia. We used these data to first compare associations of variants in genes encoding drug targets with the effects of pharmacological manipulation of those targets in clinical trials. We then tested the association those variants with disease outcomes, including coronary heart disease, to predict cardiovascular safety of these agents. A low-frequency missense variant (Ala316Thr;rs10305492) in the gene encoding glucagon-like peptide-1 receptor (GLP1R), the target of GLP1R agonists, was associated with lower fasting glucose and lower T2D risk, consistent with GLP1R agonist therapies. The minor allele was also associated with protection against heart disease, thus providing evidence that GLP1R agonists are not likely to be associated with an unacceptable increase in cardiovascular risk. Our results provide an encouraging signal that these agents may be associated with benefit, a question currently being addressed in randomised controlled trials. Genetic variants associated with metabolic traits and multiple disease outcomes can be used to validate therapeutic targets at an early stage in the drug development process.

Genotype imputation has the potential to assess human genetic variation at a lower cost than assaying the variants using laboratory techniques. The performance of imputation for rare variants has not been comprehensively studied. We utilized 8865 human samples with high depth resequencing data for the exons and flanking regions of 202 genes and Genome-Wide Association Study (GWAS) data to characterize the performance of genotype imputation for rare variants. We evaluated reference sets ranging from 100 to 3713 subjects for imputing into samples typed for the Affymetrix (500K and 6.0) and Illumina 550K GWAS panels. The proportion of variants that could be well imputed (true r2>0.7) with a reference panel of 3713 individuals was: 31% (Illumina 550K) or 25% (Affymetrix 500K) with MAF (Minor Allele Frequency) less than or equal 0.001, 48% or 35% with 0.0010.05. The performance for common SNPs (MAF>0.05) within exons and flanking regions is comparable to imputation of more uniformly distributed SNPs. The performance for rare SNPs (0.01

Hereditary tyrosinemia type 1 (HT1) is a severe autosomal recessive metabolic disease associated with point mutations in the human fumarylacetoacetate hydrolase (FAH) gene that disrupt tyrosine catabolism. An acute form of HT1 results in death during the first months of life because of hepatic failure, whereas a chronic form leads to gradual development of liver disease often accompanied by renal dysfunction, childhood rickets, neurological crisis, and hepatocellular carcinoma. Mice homozygous for certain chromosome 7 deletions of the albino Tyr; c locus that also include Fah die perinatally as a result of liver dysfunction and exhibit a complex syndrome characterized by structural abnormalities and alterations in gene expression in the liver and kidney. Here we report that two independent, postnatally lethal mutations induced by N-ethyl-N-nitrosourea and mapped near Tyr are alleles of Fah. The Fah 6287SB allele is a missense mutation in exon 6, and Fah 5961SB is a splice mutation causing loss of exon 7, a subsequent frameshift in the resulting mRNA, and a severe reduction of Fah mRNA levels. Increased levels of the diagnostic metabolite succinylacetone in the urine of the Fah 6287SB and Fah 5961SB mutants indicate that these mutations cause a decrease in Fah enzymatic activity. Thus, the neonatal phenotype present in both mutants is due to a deficiency in Fah caused by a point mutation, and we propose Fah 5961SB and Fah 6287SB as mouse models for acute and chronic forms of human HT1, respectively.

No single major HLA-related genetic risk factor was identified for lamotrigine-induced SCARs in patients of European origin. Only suggestive evidence was obtained for B*5801, A*6801, Cw*0718, DQB1*0609, and DRB1*1301. Confirmation of these results in a larger, independent sample is needed to determine whether any of the HLA alleles identified are truly associated with the development of lamotrigine-induced SCARs.

Elevated plasma levels of lipoprotein-associated phospholipase A2 (Lp-PLA2) activity have been shown to be associated with increased risk of coronary heart disease and an inhibitor of this enzyme is under development for the treatment of that condition. A Val279Phe null allele in this gene, that may influence patient eligibility for treatment, is relatively common in East Asians but has not been observed in Europeans. We investigated the existence and functional effects of low frequency alleles in a Western European population by re-sequencing the exons of PLA2G7 in 2000 samples. In all, 19 non-synonymous single-nucleotide polymorphisms (nsSNPs) were found, 14 in fewer than four subjects (minor allele frequency <0.1%). Lp-PLA2 activity was significantly lower in rare nsSNP carriers compared with non-carriers (167.8±63.2 vs 204.6±41.8, P=0.01) and seven variants had enzyme activities consistent with a null allele. The cumulative frequency of these null alleles was 0.25%, so <1 in 10 000 Europeans would be expected to be homozygous, and thus not potentially benefit from treatment with an Lp-PLA2 inhibitor.