Genetic Epidemiology最新文献

Persistent opioid use after surgery is a common morbidity outcome associated with subsequent opioid use disorder, overdose, and death. While phenotypic associations have been described, genetic associations remain unidentified. Here, we conducted the largest genetic study of persistent opioid use after surgery, comprising ~40,000 non-Hispanic, European-ancestry Michigan Genomics Initiative participants (3198 cases and 36,321 surgically exposed controls). Our study primarily focused on the reproducibility and reliability of 72 genetic studies of opioid use disorder phenotypes. Nominal associations (p < 0.05) occurred at 12 of 80 unique (r² < 0.8) signals from these studies. Six occurred in OPRM1 (most significant: rs79704991-T, OR = 1.17, p = 8.7 × 10^-5), with two surviving multiple testing correction. Other associations were rs640561-LRRIQ3 (p = 0.015), rs4680-COMT (p = 0.016), rs9478495 (p = 0.017, intergenic), rs10886472-GRK5 (p = 0.028), rs9291211-SLC30A9/BEND4 (p = 0.043), and rs112068658-KCNN1 (p = 0.048). Two highly referenced genes, OPRD1 and DRD2/ANKK1, had no signals in MGI. Associations at previously identified OPRM1 variants suggest common biology between persistent opioid use and opioid use disorder, further demonstrating connections between opioid dependence and addiction phenotypes. Lack of significant associations at other variants challenges previous studies' reliability.

Biological datasets often consist of thousands or millions of variables, e.g. genetic variants or biomarkers, and when sample sizes are large it is common to find many associated with an outcome of interest, for example, disease risk in a GWAS, at high levels of statistical significance, but with very small effects. The False Discovery Rate (FDR) is used to identify effects of interest based on ranking variables according to their statistical significance. Here, we develop a complementary measure to the FDR, the priorityFDR, that ranks variables by a combination of effect size and significance, allowing further prioritisation among a set of variables that pass a significance or FDR threshold. Applying to the largest GWAS of type 1 diabetes to date (15,573 cases and 158,408 controls), we identified 26 independent genetic associations, including two newly-reported loci, with qualitatively lower priorityFDRs than the remaining 175 signals. We detected putatively causal type 1 diabetes risk genes using Mendelian Randomisation, and found that these were located disproportionately close to low priorityFDR signals (p = 0.005), as were genes in the IL-2 pathway (p = 0.003). Selecting variables on both effect size and significance can lead to improved prioritisation for mechanistic follow-up studies from genetic and other large biological datasets.

In large cohort studies the number of unaffected individuals outnumbers the number of affected individuals, and the power can be low to detect associations for outcomes with low prevalence. We consider how including recorded family history in regression models increases the power to detect associations between genetic variants and disease risk. We show theoretically and using Monte-Carlo simulations that including a family history of the disease, with a weighting of 0.5 compared with true cases, increases the power to detect associations. This is a powerful approach for detecting variants with moderate effects, but for larger effect sizes a weighting of > 0.5 can be more powerful. We illustrate this both for common variants and for exome sequencing data for over 400,000 individuals in UK Biobank to evaluate the association between the burden of protein-truncating variants in genes and risk for four cancer types.

Gorfine, M., Qu, C.,Peters, U., & Hsu, L. (2024). Unveiling challenges in Mendelian randomization for gene-environment interaction. Genetic Epidemiology, 48, 164–189. https://doi.org/10.1002/gepi.22552

Genome-wide association studies (GWAS) have provided large numbers of genetic markers that can be used as instrumental variables in a Mendelian Randomisation (MR) analysis to assess the causal effect of a risk factor on an outcome. An extension of MR analysis, multivariable MR, has been proposed to handle multiple risk factors. However, adjusting or stratifying the outcome on a variable that is associated with it may induce collider bias. For an outcome that represents progression of a disease, conditioning by selecting only the cases may cause a biased MR estimation of the causal effect of the risk factor of interest on the progression outcome. Recently, we developed instrument effect regression and corrected weighted least squares (CWLS) to adjust for collider bias in observational associations. In this paper, we highlight the importance of adjusting for collider bias in MR with a risk factor of interest and disease progression as the outcome. A generalised version of the instrument effect regression and CWLS adjustment is proposed based on a multivariable MR model. We highlight the assumptions required for this approach and demonstrate its utility for bias reduction. We give an illustrative application to the effect of smoking initiation and smoking cessation on Crohn's disease prognosis, finding no evidence to support a causal effect.

(2024), The 2024 Annual Meeting of the International Genetic Epidemiology Society. Genetic Epidemiology, 48: 344-398. https://doi.org/10.1002/gepi.22598

In the originally-published article, several abstracts were inadvertently left out. They appear on the following pages.

We apologize for this error.

The main goal of fine-mapping is the identification of relevant genetic variants that have a causal effect on some trait of interest, such as the presence of a disease. From a statistical point of view, fine mapping can be seen as a variable selection problem. Fine-mapping methods are often challenging to apply because of the presence of linkage disequilibrium (LD), that is, regions of the genome where the variants interrogated have high correlation. Several methods have been proposed to address this issue. Here we explore the ‘Sum of Single Effects’ (SuSiE) method, applied to real data (summary statistics) from a genome-wide meta-analysis of the autoimmune liver disease primary biliary cholangitis (PBC). Fine-mapping in this data set was previously performed using the FINEMAP program; we compare these previous results with those obtained from SuSiE, which provides an arguably more convenient and principled way of generating ‘credible sets’, that is set of predictors that are correlated with the response variable. This allows us to appropriately acknowledge the uncertainty when selecting the causal effects for the trait. We focus on the results from SuSiE-RSS, which fits the SuSiE model to summary statistics, such as z-scores, along with a correlation matrix. We also compare the SuSiE results to those obtained using a more recently developed method, h2-D2, which uses the same inputs. Overall, we find the results from SuSiE-RSS and, to a lesser extent, h2-D2, to be quite concordant with those previously obtained using FINEMAP. The resulting genes and biological pathways implicated are therefore also similar to those previously obtained, providing valuable confirmation of these previously reported results. Detailed examination of the credible sets identified suggests that, although for the majority of the loci (33 out of 56) the results from SuSiE-RSS seem most plausible, there are some loci (5 out of 56 loci) where the results from h2-D2 seem more compelling. Computer simulations suggest that, overall, SuSiE-RSS generally has slightly higher power, better precision, and better ability to identify the true number of causal variants in a region than h2-D2, although there are some scenarios where the power of h2-D2 is higher. Thus, in real data analysis, the use of complementary approaches such as both SuSiE and h2-D2 is potentially warranted.

Many statistical genetics analysis methods make use of GWAS summary statistics. Best statistical practice requires evaluating these methods in realistic simulation experiments. However, simulating summary statistics by first simulating individual genotype and phenotype data is extremely computationally demanding. This high cost may force researchers to conduct overly simplistic simulations that fail to accurately measure method performance. Alternatively, summary statistics can be simulated directly from their theoretical distribution. Although this is a common need among statistical genetics researchers, no software packages exist for comprehensive GWAS summary statistic simulation. We present GWASBrewer, an open source R package for direct simulation of GWAS summary statistics. We show that statistics simulated by GWASBrewer have the same distribution as statistics generated from individual level data, and can be produced at a fraction of the computational expense. Additionally, GWASBrewer can simulate standard error estimates, something that is typically not done when sampling summary statistics directly. GWASBrewer is highly flexible, allowing the user to simulate data for multiple traits connected by causal effects and with complex distributions of effect sizes. We demonstrate example uses of GWASBrewer for evaluating Mendelian randomization, polygenic risk score, and heritability estimation methods.