Faculty Bibliography

BACKGROUND:Statins increase the risk of new-onset type 2 diabetes mellitus. We aimed to assess whether this increase in risk is a consequence of inhibition of 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), the intended drug target. METHODS:We used single nucleotide polymorphisms in the HMGCR gene, rs17238484 (for the main analysis) and rs12916 (for a subsidiary analysis) as proxies for HMGCR inhibition by statins. We examined associations of these variants with plasma lipid, glucose, and insulin concentrations; bodyweight; waist circumference; and prevalent and incident type 2 diabetes. Study-specific effect estimates per copy of each LDL-lowering allele were pooled by meta-analysis. These findings were compared with a meta-analysis of new-onset type 2 diabetes and bodyweight change data from randomised trials of statin drugs. The effects of statins in each randomised trial were assessed using meta-analysis. FINDINGS/RESULTS:Data were available for up to 223 463 individuals from 43 genetic studies. Each additional rs17238484-G allele was associated with a mean 0·06 mmol/L (95% CI 0·05-0·07) lower LDL cholesterol and higher body weight (0·30 kg, 0·18-0·43), waist circumference (0·32 cm, 0·16-0·47), plasma insulin concentration (1·62%, 0·53-2·72), and plasma glucose concentration (0·23%, 0·02-0·44). The rs12916 SNP had similar effects on LDL cholesterol, bodyweight, and waist circumference. The rs17238484-G allele seemed to be associated with higher risk of type 2 diabetes (odds ratio [OR] per allele 1·02, 95% CI 1·00-1·05); the rs12916-T allele association was consistent (1·06, 1·03-1·09). In 129 170 individuals in randomised trials, statins lowered LDL cholesterol by 0·92 mmol/L (95% CI 0·18-1·67) at 1-year of follow-up, increased bodyweight by 0·24 kg (95% CI 0·10-0·38 in all trials; 0·33 kg, 95% CI 0·24-0·42 in placebo or standard care controlled trials and -0·15 kg, 95% CI -0·39 to 0·08 in intensive-dose vs moderate-dose trials) at a mean of 4·2 years (range 1·9-6·7) of follow-up, and increased the odds of new-onset type 2 diabetes (OR 1·12, 95% CI 1·06-1·18 in all trials; 1·11, 95% CI 1·03-1·20 in placebo or standard care controlled trials and 1·12, 95% CI 1·04-1·22 in intensive-dose vs moderate dose trials). INTERPRETATION/CONCLUSIONS:The increased risk of type 2 diabetes noted with statins is at least partially explained by HMGCR inhibition. FUNDING/BACKGROUND:The funding sources are cited at the end of the paper.

Genome-wide association studies (GWASs) have identified many genetic variants underlying complex traits. Many detected genetic loci harbor variants that associate with multiple-even distinct-traits. Most current analysis approaches focus on single traits, even though the final results from multiple traits are evaluated together. Such approaches miss the opportunity to systemically integrate the phenome-wide data available for genetic association analysis. In this study, we propose a general approach that can integrate association evidence from summary statistics of multiple traits, either correlated, independent, continuous, or binary traits, which might come from the same or different studies. We allow for trait heterogeneity effects. Population structure and cryptic relatedness can also be controlled. Our simulations suggest that the proposed method has improved statistical power over single-trait analysis in most of the cases we studied. We applied our method to the Continental Origins and Genetic Epidemiology Network (COGENT) African ancestry samples for three blood pressure traits and identified four loci (CHIC2, HOXA-EVX1, IGFBP1/IGFBP3, and CDH17; p < 5.0 × 10(-8)) associated with hypertension-related traits that were missed by a single-trait analysis in the original report. Six additional loci with suggestive association evidence (p < 5.0 × 10(-7)) were also observed, including CACNA1D and WNT3. Our study strongly suggests that analyzing multiple phenotypes can improve statistical power and that such analysis can be executed with the summary statistics from GWASs. Our method also provides a way to study a cross phenotype (CP) association by using summary statistics from GWASs of multiple phenotypes.

BACKGROUND:Leber congenital amaurosis (LCA) is a severe form of retinal dystrophy with marked underlying genetic heterogeneity. Until recently, allele-specific assays and Sanger sequencing of targeted segments were the only available approaches for attempted genetic diagnosis in this condition. A broader next-generation sequencing (NGS) strategy, such as whole exome sequencing, provides an improved molecular genetic diagnostic capacity for patients with these conditions. MATERIALS AND METHODS/METHODS:In a child with LCA, an allele-specific assay analyzing 135 known LCA-causing variations, followed by targeted segment sequencing of 61 regions in 14 causative genes was performed. Subsequently, exome sequencing was undertaken in the proband, unaffected consanguineous parents and two unaffected siblings. Bioinformatic analysis used two independent pipelines, BWA-GATK and SOAP, followed by Annovar and SnpEff to annotate the variants. RESULTS:No disease-causing variants were found using the allele-specific or targeted segment Sanger sequencing assays. Analysis of variants in the exome sequence data revealed a novel homozygous nonsense mutation (c.1081C > T, p.Arg361*) in TULP1, a gene with roles in photoreceptor function where mutations were previously shown to cause LCA and retinitis pigmentosa. The identified homozygous variant was the top candidate using both bioinformatic pipelines. CONCLUSIONS:This study highlights the value of the broad sequencing strategy of exome sequencing for disease gene identification in LCA, over other existing methods. NGS is particularly beneficial in LCA where there are a large number of causative disease genes, few distinguishing clinical features for precise candidate disease gene selection, and few mutation hotspots in any of the known disease genes.

This study aimed to identify the genetics underlying dominant forms of inherited retinal dystrophies using whole exome sequencing (WES) in six families extensively screened for known mutations or genes. Thirty-eight individuals were subjected to WES. Causative variants were searched among single nucleotide variants (SNVs) and insertion/deletion variants (indels) and whenever no potential candidate emerged, copy number variant (CNV) analysis was performed. Variants or regions harboring a candidate variant were prioritized and segregation of the variant with the disease was further assessed using Sanger sequencing in case of SNVs and indels, and quantitative PCR (qPCR) for CNVs. SNV and indel analysis led to the identification of a previously reported mutation in PRPH2. Two additional mutations linked to different forms of retinal dystrophies were identified in two families: a known frameshift deletion in RPGR, a gene responsible for X-linked retinitis pigmentosa and p.Ser163Arg in C1QTNF5 associated with Late-Onset Retinal Degeneration. A novel heterozygous deletion spanning the entire region of PRPF31 was also identified in the affected members of a fourth family, which was confirmed with qPCR. This study allowed the identification of the genetic cause of the retinal dystrophy and the establishment of a correct diagnosis in four families, including a large heterozygous deletion in PRPF31, typically considered one of the pitfalls of this method. Since all findings in this study are restricted to known genes, we propose that targeted sequencing using gene-panel is an optimal first approach for the genetic screening and that once known genetic causes are ruled out, WES might be used to uncover new genes involved in inherited retinal dystrophies.

BACKGROUND:Despite heritability estimates of 40-70 % for obesity, less than 2 % of its variation is explained by Body Mass Index (BMI) associated loci that have been identified so far. Epistasis, or gene-gene interactions are a plausible source to explain portions of the missing heritability of BMI. METHODS:Using genotypic data from 18,686 individuals across five study cohorts - ARIC, CARDIA, FHS, CHS, MESA - we filtered SNPs (Single Nucleotide Polymorphisms) using two parallel approaches. SNPs were filtered either on the strength of their main effects of association with BMI, or on the number of knowledge sources supporting a specific SNP-SNP interaction in the context of BMI. Filtered SNPs were specifically analyzed for interactions that are highly associated with BMI using QMDR (Quantitative Multifactor Dimensionality Reduction). QMDR is a nonparametric, genetic model-free method that detects non-linear interactions associated with a quantitative trait. RESULTS:We identified seven novel, epistatic models with a Bonferroni corrected p-value of association < 0.1. Prior experimental evidence helps explain the plausible biological interactions highlighted within our results and their relationship with obesity. We identified interactions between genes involved in mitochondrial dysfunction (POLG2), cholesterol metabolism (SOAT2), lipid metabolism (CYP11B2), cell adhesion (EZR), cell proliferation (MAP2K5), and insulin resistance (IGF1R). Moreover, we found an 8.8 % increase in the variance in BMI explained by these seven SNP-SNP interactions, beyond what is explained by the main effects of an index FTO SNP and the SNPs within these interactions. We also replicated one of these interactions and 58 proxy SNP-SNP models representing it in an independent dataset from the eMERGE study. CONCLUSION/CONCLUSIONS:This study highlights a novel approach for discovering gene-gene interactions by combining methods such as QMDR with traditional statistics.

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22,981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10(-8)). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine.

White blood cell (WBC) count is a common clinical measure used as a predictor of certain aspects of human health, including immunity and infection status. WBC count is also a complex trait that varies among individuals and ancestry groups. Differences in linkage disequilibrium structure and heterogeneity in allelic effects are expected to play a role in the associations observed between populations. Prior genome-wide association study (GWAS) meta-analyses have identified genomic loci associated with WBC and its subtypes, but much of the heritability of these phenotypes remains unexplained. Using GWAS summary statistics for over 50 000 individuals from three diverse populations (Japanese, African-American and European ancestry), a Bayesian model methodology was employed to account for heterogeneity between ancestry groups. This approach was used to perform a trans-ethnic meta-analysis of total WBC, neutrophil and monocyte counts. Ten previously known associations were replicated and six new loci were identified, including several regions harboring genes related to inflammation and immune cell function. Ninety-five percent credible interval regions were calculated to narrow the association signals and fine-map the putatively causal variants within loci. Finally, a conditional analysis was performed on the most significant SNPs identified by the trans-ethnic meta-analysis (MA), and nine secondary signals within loci previously associated with WBC or its subtypes were identified. This work illustrates the potential of trans-ethnic analysis and ascribes a critical role to multi-ethnic cohorts and consortia in exploring complex phenotypes with respect to variants that lie outside the European-biased GWAS pool.

BACKGROUND:Phosphoribosyl pyrophosphate synthetase (PRS) I deficiency is a rare medical condition caused by missense mutations in PRPS1 that lead to three different phenotypes: Arts Syndrome (MIM 301835), X-linked Charcot-Marie-Tooth (CMTX5, MIM 311070) or X-linked non-syndromic sensorineural deafness (DFN2, MIM 304500). All three are X-linked recessively inherited and males affected display variable degree of central and peripheral neuropathy. We applied whole exome sequencing to a three-generation family with optic atrophy followed by retinitis pigmentosa (RP) in all three cases, and ataxia, progressive peripheral neuropathy and hearing loss with variable presentation. METHODS:Whole exome sequencing was performed in two affecteds and one unaffected member of the family. Sanger sequencing was used to validate and segregate the 12 candidate mutations in the family and to confirm the absence of the novel variant in PRPS1 in 191 controls. The pathogenic role of the novel mutation in PRPS1 was assessed in silico and confirmed by enzymatic determination of PRS activity, mRNA expression and sequencing, and X-chromosome inactivation. RESULTS:A novel missense mutation was identified in PRPS1 in the affected females. Age of onset, presentation and severity of the phenotype are highly variable in the family: both the proband and her mother have neurological and ophthalmological symptoms, whereas the phenotype of the affected sister is milder and currently confined to the eye. Moreover, only the proband displayed a complete lack of expression of the wild type allele in leukocytes that seems to correlate with the degree of PRS deficiency and the severity of the phenotype. Interestingly, optic atrophy and RP are the only common manifestations to all three females and the only phenotype correlating with the degree of enzyme deficiency. CONCLUSIONS:These results are in line with recent evidence of the existence of intermediate phenotypes in PRS-I deficiency syndromes and demonstrate that females can exhibit a disease phenotype as severe and complex as their male counterparts.