BR-squared: a practical solution to the winner’s curse in genome-wide scans

The detrimental effects of the winner’s curse, including overestimation of the genetic effects of associated variants and underestimation of sufficient sample sizes for replication studies are well-recognized in genome-wide association studies (GWAS). These effects can be expected to worsen as the field moves from GWAS into whole genome sequencing. To date, few studies have reported statistical adjustments to the naive estimates, due to the lack of suitable statistical methods and computational tools. We have developed an efficient genome-wide non-parametric method that explicitly accounts for the threshold, ranking, and allele frequency effects in whole genome scans. Here, we implement the method to provide bias-reduced estimates via bootstrap re-sampling (BR-squared) for association studies of both disease status and quantitative traits, and we report the results of applying BR-squared to GWAS of psoriasis and HbA1c. We observed over 50% reduction in the genetic effect size estimation for many associated SNPs. This translates into a greater than fourfold increase in sample size requirements for successful replication studies, which in part explains some of the apparent failures in replicating the original signals. Our analysis suggests that adjusting for the winner’s curse is critical for interpreting findings from whole genome scans and planning replication and meta-GWAS studies, as well as in attempts to translate findings into the clinical setting

Common Genetic Variation Near the Phospholamban Gene Is Associated with Cardiac Repolarisation: Meta-Analysis of Three Genome-Wide Association Studies

To identify loci affecting the electrocardiographic QT interval, a measure of cardiac repolarisation associated with risk of ventricular arrhythmias and sudden cardiac death, we conducted a meta-analysis of three genome-wide association studies (GWAS) including 3,558 subjects from the TwinsUK and BRIGHT cohorts in the UK and the DCCT/EDIC cohort from North America. Five loci were significantly associated with QT interval at P<1×10$^{−6}$. To validate these findings we performed an in silico comparison with data from two QT consortia: QTSCD (n = 15,842) and QTGEN (n = 13,685). Analysis confirmed the association between common variants near NOS1AP (P = 1.4×10$^{−83}$) and the phospholamban (PLN) gene (P = 1.9×10$^{−29}$). The most associated SNP near NOS1AP (rs12143842) explains 0.82% variance; the SNP near PLN (rs11153730) explains 0.74% variance of QT interval duration. We found no evidence for interaction between these two SNPs (P = 0.99). PLN is a key regulator of cardiac diastolic function and is involved in regulating intracellular calcium cycling, it has only recently been identified as a susceptibility locus for QT interval. These data offer further mechanistic insights into genetic influence on the QT interval which may predispose to life threatening arrhythmias and sudden cardiac death

Genetic association study of QT interval highlights role for calcium signaling pathways in myocardial repolarization.

The QT interval, an electrocardiographic measure reflecting myocardial repolarization, is a heritable trait. QT prolongation is a risk factor for ventricular arrhythmias and sudden cardiac death (SCD) and could indicate the presence of the potentially lethal mendelian long-QT syndrome (LQTS). Using a genome-wide association and replication study in up to 100,000 individuals, we identified 35 common variant loci associated with QT interval that collectively explain ∼8-10% of QT-interval variation and highlight the importance of calcium regulation in myocardial repolarization. Rare variant analysis of 6 new QT interval-associated loci in 298 unrelated probands with LQTS identified coding variants not found in controls but of uncertain causality and therefore requiring validation. Several newly identified loci encode proteins that physically interact with other recognized repolarization proteins. Our integration of common variant association, expression and orthogonal protein-protein interaction screens provides new insights into cardiac electrophysiology and identifies new candidate genes for ventricular arrhythmias, LQTS and SCD

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms

New susceptibility loci associated with kidney disease in type 1 diabetes

WOS:000309817900008Diabetic kidney disease, or diabetic nephropathy (DN), is a major complication of diabetes and the leading cause of end-stage renal disease (ESRD) that requires dialysis treatment or kidney transplantation. In addition to the decrease in the quality of life, DN accounts for a large proportion of the excess mortality associated with type 1 diabetes (T1D). Whereas the degree of glycemia plays a pivotal role in DN, a subset of individuals with poorly controlled T1D do not develop DN. Furthermore, strong familial aggregation supports genetic susceptibility to DN. However, the genes and the molecular mechanisms behind the disease remain poorly understood, and current therapeutic strategies rarely result in reversal of DN. In the GEnetics of Nephropathy: an International Effort (GENIE) consortium, we have undertaken a meta-analysis of genome-wide association studies (GWAS) of T1D DN comprising ∼2.4 million single nucleotide polymorphisms (SNPs) imputed in 6,691 individuals. After additional genotyping of 41 top ranked SNPs representing 24 independent signals in 5,873 individuals, combined meta-analysis revealed association of two SNPs with ESRD: rs7583877 in the AFF3 gene (P = 1.2×10(-8)) and an intergenic SNP on chromosome 15q26 between the genes RGMA and MCTP2, rs12437854 (P = 2.0×10(-9)). Functional data suggest that AFF3 influences renal tubule fibrosis via the transforming growth factor-beta (TGF-β1) pathway. The strongest association with DN as a primary phenotype was seen for an intronic SNP in the ERBB4 gene (rs7588550, P = 2.1×10(-7)), a gene with type 2 diabetes DN differential expression and in the same intron as a variant with cis-eQTL expression of ERBB4. All these detected associations represent new signals in the pathogenesis of DN.Peer reviewe

Systematic evaluation of medium-throughput mRNA abundance platforms

Profiling of mRNA abundances with high-throughput platforms such as microarrays and RNA-seq has become an important tool in both basic and biomedical research. However, these platforms remain prone to systematic errors and have challenges in clinical and industrial applications. As a result, it is standard practice to validate a subset of key results using alternate technologies. Similarly, clinical and industrial applications typically involve transitions from a high-throughput discovery platform to medium-throughput validation ones. These medium-throughput validation platforms have high technical reproducibility and reduced sample input needs, and low sensitivity to sample quality (e.g., for processing FFPE specimens). Unfortunately, while medium-throughput platforms have proliferated, there are no comprehensive comparisons of them. Here we fill that gap by comparing two key medium-throughput platforms—NanoString's nCounter Analysis System and ABI's OpenArray System—to gold-standard quantitative real-time RT-PCR. We quantified 38 genes and positive and negative controls in 165 samples. Signal:noise ratios, correlations, dynamic range, and detection accuracy were compared across platforms. All three measurement technologies showed good concordance, but with divergent price/time/sensitivity trade-offs. This study provides the first detailed comparison of medium-throughput RNA quantification platforms and provides a template and a standard data set for the evaluation of additional technologies

Pilot genome-wide association search identifies potential loci for risk of erectile dysfunction in type 1 diabetes using the DCCT/EDIC study cohort.

PURPOSE: To identify genetic predictors of diabetes-associated ED using genome wide and candidate gene approaches in a cohort of men with type I diabetes. METHODS: We examined 528 white men with T1D (125 with ED) from the DCCT and its observational follow up EDIC Study. ED was defined from a single item of the IIEF. An Illumina Human1M BeadChip was used for genotyping. 867,125 single nucleotide polymorphisms (SNPs) were subjected to analysis. Whole genome and candidate gene approaches tested the hypothesis that genetic polymorphisms may predispose men with T1D to ED. Univariate and multivariate models were used controlling for age, HbA1c, diabetes duration, and prior randomization to intensive or conventional insulin therapy during DCCT. A stratified false discovery rate was used to perform the candidate gene approach. RESULTS: Two SNPs located on chromosome 3 in one genomic loci were associated with ED with p < 1×10(−6). rs9810233 had a p-value of 7 × 10(−7) and rs1920201 had a p-value of 9×10(−7) The nearest gene to these two SNPs is ALCAM. The genetic association results at these loci were similar in univariate and multivariate analysis. No candidate genes met criteria for statistical significance. CONCLUSIONS: Two SNPs, rs9810233 and rs1920101, which are 25 kb apart, are both associated with ED, albeit not meeting the standard GWAS significance criteria of p < 5 × 10(−8). Other studies with larger sample sizes will be required to determine whether ALCAM represents a novel gene in the pathogenesis of diabetes associated ED

Identification of a microRNA signature associated with risk of distant metastasis in nasopharyngeal carcinoma.

PurposeDespite significant improvement in locoregional control in the contemporary era of nasopharyngeal carcinoma (NPC) treatment, patients still suffer from a significant risk of distant metastasis (DM). Identifying those patients at risk of DM would aid in personalized treatment in the future. MicroRNAs (miRNAs) play many important roles in human cancers; hence, we proceeded to address the primary hypothesis that there is a miRNA expression signature capable of predicting DM for NPC patients.Methods and resultsThe expression of 734 miRNAs was measured in 125 (Training) and 121 (Validation) clinically annotated NPC diagnostic biopsy samples. A 4-miRNA expression signature associated with risk of developing DM was identified by fitting a penalized Cox Proportion Hazard regression model to the Training data set (HR 8.25; p &lt; 0.001), and subsequently validated in an independent Validation set (HR 3.2; p = 0.01). Pathway enrichment analysis indicated that the targets of miRNAs associated with DM appear to be converging on cell-cycle pathways.ConclusionsThis 4-miRNA signature adds to the prognostic value of the current "gold standard" of TNM staging. In-depth interrogation of these 4-miRNAs will provide important biological insights that could facilitate the discovery and development of novel molecularly targeted therapies to improve outcome for future NPC patients