Jennifer L. Aponte

An Abundance of Rare Functional Variants in 202 Drug Target Genes Sequenced in 14,002 People

A Deep Look Into Our Genes Recent debates have focused on the degree of genetic variation and its impact upon health at the genomic level in humans (see the Perspective by Casals and Bertranpetit). Tennessen et al. (p. 64, published online 17 May), looking at all of the protein-coding genes in the human genome, and Nelson et al. (p. 100, published online 17 May), looking at genes that encode drug targets, address this question through deep sequencing efforts on samples from multiple individuals. The findings suggest that most human variation is rare, not shared between populations, and that rare variants are likely to play a role in human health. A pharmacogenomics analysis shows how challenging it will be to associate rare variants with phenotypes. 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 (1 every 17 bases) and geographically localized, so 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. We conclude that because of 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.

High-resolution HLA genotyping and severe cutaneous adverse reactions in lamotrigine-treated patients

Background Severe cutaneous adverse reactions (SCARs) are associated with over 200 medicines including lamotrigine, an antiepileptic drug. Previous studies have suggested the involvement of immune mechanisms in the development of drug-induced SCARs. Methods High-resolution HLA genotyping was performed for 65 patients of European ancestry treated with lamotrigine (22 cases with lamotrigine-induced SCARs and 43 controls on lamotrigine without SCAR-related symptoms). Association of HLA genetic variants with SCARs in these patients were evaluated by contrasting allele frequencies between the cases and the controls for each of 112 HLA four-digit alleles. Results Five alleles were observed with higher frequencies in the cases compared with the treated controls with exact P values less than 0.05. These include B*5801 (P = 0.037), previously reported to be associated with allopurinol-induced SCARs. Marginal association evidence was also observed for alleles Cw*0718 and DQB1*0609, both of which were strongly correlated with B*5801. Other alleles identified were A*6801 (P = 0.012) and DRB1*1301 (P = 0.045). In contrast to the study of carbamazepine-induced Stevens–Johnson syndrome in Han Chinese patients, none of the cases carried B*1502. Accounting for the large number of hypothesis tests conducted, none of the associations identified were statistically significant. Conclusion 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.

Performance of genotype imputation for rare variants identified in exons and flanking regions of genes.

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.001<MAF< = 0.005, 54% or 38% with 0.005<MAF< = 0.01, 78% or 57% with 0.01<MAF< = 0.05, and 97% or 86% with MAF>0.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<MAF< = 0.05) was much more dependent on the GWAS panel and the number of reference samples. These results suggest routine use of genotype imputation for extending the assessment of common variants identified in humans via targeted exon resequencing into additional samples with GWAS data, but imputation of very rare variants (MAF< = 0.005) will require reference panels with thousands of subjects.

Deep Resequencing Unveils Genetic Architecture of ADIPOQ and Identifies a Novel Low-Frequency Variant Strongly Associated With Adiponectin Variation

Increased adiponectin levels have been shown to be associated with a lower risk of type 2 diabetes. To understand the relations between genetic variation at the adiponectin-encoding gene, ADIPOQ, and adiponectin levels, and subsequently its role in disease, we conducted a deep resequencing experiment of ADIPOQ in 14,002 subjects, including 12,514 Europeans, 594 African Americans, and 567 Indian Asians. We identified 296 single nucleotide polymorphisms (SNPs), including 30 amino acid changes, and carried out association analyses in a subset of 3,665 subjects from two independent studies. We confirmed multiple genome-wide association study findings and identified a novel association between a low-frequency SNP (rs17366653) and adiponectin levels (P = 2.2E–17). We show that seven SNPs exert independent effects on adiponectin levels. Together, they explained 6% of adiponectin variation in our samples. We subsequently assessed association between these SNPs and type 2 diabetes in the Genetics of Diabetes Audit and Research in Tayside Scotland (GO-DARTS) study, comprised of 5,145 case and 6,374 control subjects. No evidence of association with type 2 diabetes was found, but we were also unable to exclude the possibility of substantial effects (e.g., odds ratio 95% CI for rs7366653 [0.91–1.58]). Further investigation by large-scale and well-powered Mendelian randomization studies is warranted.

Deep sequencing of the LRRK2 gene in 14,002 individuals reveals evidence of purifying selection and independent origin of the p.Arg1628Pro mutation in Europe†

Genetic variation in LRRK2 predisposes to Parkinson disease (PD), which underpins its development as a therapeutic target. Here, we aimed to identify novel genotype–phenotype associations that might support developing LRRK2 therapies for other conditions. We sequenced the 51 exons of LRRK2 in cases comprising 12 common diseases (n = 9,582), and in 4,420 population controls. We identified 739 single‐nucleotide variants, 62% of which were observed in only one person, including 316 novel exonic variants. We found evidence of purifying selection for the LRRK2 gene and a trend suggesting that this is more pronounced in the central (ROC–COR–kinase) core protein domains of LRRK2 than the flanking domains. Population genetic analyses revealed that LRRK2 is not especially polymorphic or differentiated in comparison to 201 other drug target genes. Among Europeans, we identified 17 carriers (0.13%) of pathogenic LRRK2 mutations that were not significantly enriched within any disease or in those reporting a family history of PD. Analysis of pathogenic mutations within Europe reveals that the p.Arg1628Pro (c4883G>C) mutation arose independently in Europe and Asia. Taken together, these findings demonstrate how targeted deep sequencing can help to reveal fundamental characteristics of clinically important loci. Hum Mutat 33:1087–1098, 2012.

A Low-Frequency Variant in MAPK14 Provides Mechanistic Evidence of a Link With Myeloperoxidase: A Prognostic Cardiovascular Risk Marker

BACKGROUND Genetics can be used to predict drug effects and generate hypotheses around alternative indications. To support Losmapimod, a p38 mitogen-activated protein kinase inhibitor in development for acute coronary syndrome, we characterized gene variation in MAPK11/14 genes by exome sequencing and follow-up genotyping or imputation in participants well-phenotyped for cardiovascular and metabolic traits. METHODS AND RESULTS Investigation of genetic variation in MAPK11 and MAPK14 genes using additive genetic models in linear or logistic regression with cardiovascular, metabolic, and biomarker phenotypes highlighted an association of RS2859144 in MAPK14 with myeloperoxidase in a dyslipidemic population (Genetic Epidemiology of Metabolic Syndrome Study), P=2.3×10(-6)). This variant (or proxy) was consistently associated with myeloperoxidase in the Framingham Heart Study and Cardiovascular Health Study studies (replication meta-P=0.003), leading to a meta-P value of 9.96×10(-7) in the 3 dyslipidemic groups. The variant or its proxy was then profiled in additional population-based cohorts (up to a total of 58 930 subjects) including Cohorte Lausannoise, Ely, Fenland, European Prospective Investigation of Cancer, London Life Sciences Prospective Population Study, and the Genetics of Obesity Associations study obesity case-control for up to 40 cardiovascular and metabolic traits. Overall analysis identified the same single nucleotide polymorphisms to be nominally associated consistently with glomerular filtration rate (P=0.002) and risk of obesity (body mass index ≥30 kg/m(2), P=0.004). CONCLUSIONS As myeloperoxidase is a prognostic marker of coronary events, the MAPK14 variant may provide a mechanistic link between p38 map kinase and these events, providing information consistent with current indication of Losmapimod for acute coronary syndrome. If replicated, the association with glomerular filtration rate, along with previous biological findings, also provides support for kidney diseases as alternative indications.

Pharmacogenetic meta-analysis of baseline risk factors, pharmacodynamic, efficacy and tolerability endpoints from two large global cardiovascular outcomes trials for darapladib

Darapladib, a lipoprotein-associated phospholipase A2 (Lp-PLA2) inhibitor, failed to demonstrate efficacy for the primary endpoints in two large phase III cardiovascular outcomes trials, one in stable coronary heart disease patients (STABILITY) and one in acute coronary syndrome (SOLID-TIMI 52). No major safety signals were observed but tolerability issues of diarrhea and odor were common (up to 13%). We hypothesized that genetic variants associated with Lp-PLA2 activity may influence efficacy and tolerability and therefore performed a comprehensive pharmacogenetic analysis of both trials. We genotyped patients within the STABILITY and SOLID-TIMI 52 trials who provided a DNA sample and consent (n = 13,577 and 10,404 respectively, representing 86% and 82% of the trial participants) using genome-wide arrays with exome content and performed imputation using a 1000 Genomes reference panel. We investigated baseline and change from baseline in Lp-PLA2 activity, two efficacy endpoints (major coronary events and myocardial infarction) as well as tolerability parameters at genome-wide and candidate gene level using a meta-analytic approach. We replicated associations of published loci on baseline Lp-PLA2 activity (APOE, CELSR2, LPA, PLA2G7, LDLR and SCARB1) and identified three novel loci (TOMM5, FRMD5 and LPL) using the GWAS-significance threshold P≤5E-08. Review of the PLA2G7 gene (encoding Lp-PLA2) within these datasets identified V279F null allele carriers as well as three other rare exonic null alleles within various ethnic groups, however none of these variants nor any other loci associated with Lp-PLA2 activity at baseline were associated with any of the drug response endpoints. The analysis of darapladib efficacy endpoints, despite low power, identified six low frequency loci with main genotype effect (though with borderline imputation scores) and one common locus (minor allele frequency 0.24) with genotype by treatment interaction effect passing the GWAS-significance threshold. This locus conferred risk in placebo subjects, hazard ratio (HR) 1.22 with 95% confidence interval (CI) 1.11–1.33, but was protective in darapladib subjects, HR 0.79 (95% CI 0.71–0.88). No major loci for tolerability were found. Thus, genetic analysis confirmed and extended the influence of lipoprotein loci on Lp-PLA2 levels, identified some novel null alleles in the PLA2G7 gene, and only identified one potentially efficacious subgroup within these two large clinical trials.

Assessment of rosacea symptom severity by genome-wide association study and expression analysis highlights immuno-inflammatory and skin pigmentation genes

Abstract Rosacea is a common, chronic skin disease of variable severity with limited treatment options. The cause of rosacea is unknown, but it is believed to be due to a combination of hereditary and environmental factors. Little is known about the genetics of the disease. We performed a genome-wide association study (GWAS) of rosacea symptom severity with data from 73 265 research participants of European ancestry from the 23andMe customer base. Seven loci had variants associated with rosacea at the genome-wide significance level (P < 5 × 10−8). Further analyses highlighted likely gene regions or effector genes including IRF4 (P = 1.5 × 10−17), a human leukocyte antigen (HLA) region flanked by PSMB9 and HLA-DMB (P = 2.2 × 10−15), HERC2-OCA2 (P = 4.2 × 10−12), SLC45A2 (P = 1.7 × 10−10), IL13 (P = 2.8 × 10−9), a region flanked by NRXN3 and DIO2 (P = 4.1 × 10−9), and a region flanked by OVOL1and SNX32 (P = 1.2 × 10−8). All associations with rosacea were novel except for the HLA locus. Two of these loci (HERC-OCA2 and SLC45A2) and another precedented variant (rs1805007 in melanocortin 1 receptor) with an association P value just below the significance threshold (P = 1.3 × 10−7) have been previously associated with skin phenotypes and pigmentation, two of these loci are linked to immuno-inflammation phenotypes (IL13 and PSMB9-HLA-DMA) and one has been associated with both categories (IRF4). Genes within three loci (PSMB9-HLA-DMA, HERC-OCA2 and NRX3-DIO2) were differentially expressed in a previously published clinical rosacea transcriptomics study that compared lesional to non-lesional samples. The identified loci provide specificity of inflammatory mechanisms in rosacea, and identify potential pathways for therapeutic intervention.