On Mendelian randomization analysis of case-control study

Mendelian randomization (MR) analysis uses genotypes as instruments to estimate the causal effect of an exposure in the presence of unobserved confounders. The existing MR methods focus on the data generated from prospective cohort studies. We develop a procedure for studying binary outcomes under a case-control design. The proposed procedure is built upon two working models commonly used for MR analyses and adopts a quasi-empirical likelihood framework to address the ascertainment bias from case-control sampling. We derive various approaches for estimating the causal effect and hypothesis testing under the empirical likelihood framework. We conduct extensive simulation studies to evaluate the proposed methods. We find that the proposed empirical likelihood estimate is less biased than the existing estimates. Among all the approaches considered, the Lagrange multiplier (LM) test has the highest power, and the confidence intervals derived from the LM test have the most accurate coverage. We illustrate the use of our method in MR analysis of prostate cancer case-control data with vitamin D level as exposure and three single nucleotide polymorphisms as instruments.     

Assortative mating and within-spouse pair comparisons

Spousal comparisons have been proposed as a design that can both reduce confounding and estimate effects of the shared adulthood environment. However, assortative mating, the process by which individuals select phenotypically (dis)similar mates, could distort associations when comparing spouses. We evaluated the use of spousal comparisons, as in the within-spouse pair (WSP) model, for aetiological research such as genetic association studies. We demonstrated that the WSP model can reduce confounding but may be susceptible to collider bias arising from conditioning on assorted spouse pairs. Analyses using UK Biobank spouse pairs found that WSP genetic association estimates were smaller than estimates from random pairs for height, educational attainment, and BMI variants. Within-sibling pair estimates, robust to demographic and parental effects, were also smaller than random pair estimates for height and educational attainment, but not for BMI. WSP models, like other within-family models, may reduce confounding from demographic factors in genetic association estimates, and so could be useful for triangulating evidence across study designs to assess the robustness of findings. However, WSP estimates should be interpreted with caution due to potential collider bias.     

Cardiometabolic risk factors for COVID-19 susceptibility and severity: A Mendelian randomization analysis

BackgroundEpidemiological studies report associations of diverse cardiometabolic conditions including obesity with COVID-19 illness, but causality has not been established. We sought to evaluate the associations of 17 cardiometabolic traits with COVID-19 susceptibility and severity using 2-sample Mendelian randomization (MR) analyses.Methods and findingsWe selected genetic variants associated with each exposure, including body mass index (BMI), at p < 5 × 10⁻⁸ from genome-wide association studies (GWASs). We then calculated inverse-variance-weighted averages of variant-specific estimates using summary statistics for susceptibility and severity from the COVID-19 Host Genetics Initiative GWAS meta-analyses of population-based cohorts and hospital registries comprising individuals with self-reported or genetically inferred European ancestry. Susceptibility was defined as testing positive for COVID-19 and severity was defined as hospitalization with COVID-19 versus population controls (anyone not a case in contributing cohorts). We repeated the analysis for BMI with effect estimates from the UK Biobank and performed pairwise multivariable MR to estimate the direct effects and indirect effects of BMI through obesity-related cardiometabolic diseases. Using p < 0.05/34 tests = 0.0015 to declare statistical significance, we found a nonsignificant association of genetically higher BMI with testing positive for COVID-19 (14,134 COVID-19 cases/1,284,876 controls, p = 0.002; UK Biobank: odds ratio 1.06 [95% CI 1.02, 1.10] per kg/m²; p = 0.004]) and a statistically significant association with higher risk of COVID-19 hospitalization (6,406 hospitalized COVID-19 cases/902,088 controls, p = 4.3 × 10⁻⁵; UK Biobank: odds ratio 1.14 [95% CI 1.07, 1.21] per kg/m², p = 2.1 × 10⁻⁵). The implied direct effect of BMI was abolished upon conditioning on the effect on type 2 diabetes, coronary artery disease, stroke, and chronic kidney disease. No other cardiometabolic exposures tested were associated with a higher risk of poorer COVID-19 outcomes. Small study samples and weak genetic instruments could have limited the detection of modest associations, and pleiotropy may have biased effect estimates away from the null.ConclusionsIn this study, we found genetic evidence to support higher BMI as a causal risk factor for COVID-19 susceptibility and severity. These results raise the possibility that obesity could amplify COVID-19 disease burden independently or through its cardiometabolic consequences and suggest that targeting obesity may be a strategy to reduce the risk of severe COVID-19 outcomes.

Aaron Leong and co-workers investigate causal risk factors for COVID-10 illness and severity.     

Mendelian randomization analysis of the causal impact of body mass index and waist-hip ratio on rates of hospital admission

We analyze how measures of adiposity – body mass index (BMI) and waist hip ratio (WHR) – causally influence rates of hospital admission. Conventional analyses of this relationship are susceptible to omitted variable bias from variables that jointly influence both hospital admission and adipose status. We implement a novel quasi-Poisson instrumental variable model in a Mendelian randomization framework, identifying causal effects from random perturbations to germline genetic variation. We estimate the individual and joint effects of BMI, WHR, and WHR adjusted for BMI. We also implement multivariable instrumental variable methods in which the causal effect of one exposure is estimated conditionally on the causal effect of another exposure. Data on 310,471 participants and over 550,000 inpatient admissions in the UK Biobank were used to perform one-sample and two-sample Mendelian randomization analyses. The results supported a causal role of adiposity on hospital admissions, with consistency across all estimates and sensitivity analyses. Point estimates were generally larger than estimates from comparable observational specifications. We observed an attenuation of the BMI effect when adjusting for WHR in the multivariable Mendelian randomization analyses, suggesting that an adverse fat distribution, rather than a higher BMI itself, may drive the relationship between adiposity and risk of hospital admission.     

Associations between multimorbidity and adverse health outcomes in UK Biobank and the SAIL Databank: A comparison of longitudinal cohort studies

BackgroundCohorts such as UK Biobank are increasingly used to study multimorbidity; however, there are concerns that lack of representativeness may lead to biased results. This study aims to compare associations between multimorbidity and adverse health outcomes in UK Biobank and a nationally representative sample.Methods and findingsThese are observational analyses of cohorts identified from linked routine healthcare data from UK Biobank participants (n = 211,597 from England, Scotland, and Wales with linked primary care data, age 40 to 70, mean age 56.5 years, 54.6% women, baseline assessment 2006 to 2010) and from the Secure Anonymised Information Linkage (SAIL) databank (n = 852,055 from Wales, age 40 to 70, mean age 54.2, 50.0% women, baseline January 2011). Multimorbidity (n = 40 long-term conditions [LTCs]) was identified from primary care Read codes and quantified using a simple count and a weighted score. Individual LTCs and LTC combinations were also assessed. Associations with all-cause mortality, unscheduled hospitalisation, and major adverse cardiovascular events (MACEs) were assessed using Weibull or negative binomial models adjusted for age, sex, and socioeconomic status, over 7.5 years follow-up for both datasets.Multimorbidity was less common in UK Biobank than SAIL (26.9% and 33.0% with ≥2 LTCs in UK Biobank and SAIL, respectively). This difference was attenuated, but persisted, after standardising by age, sex, and socioeconomic status. The association between increasing multimorbidity count and mortality, hospitalisation, and MACE was similar between both datasets at LTC counts of ≤3; however, above this level, UK Biobank underestimated the risk associated with multimorbidity (e.g., mortality hazard ratio for 2 LTCs 1.62 (95% confidence interval 1.57 to 1.68) in SAIL and 1.51 (1.43 to 1.59) in UK Biobank, hazard ratio for 5 LTCs was 3.46 (3.31 to 3.61) in SAIL and 2.88 (2.63 to 3.15) in UK Biobank). Absolute risk of mortality, hospitalisation, and MACE, at all levels of multimorbidity, was lower in UK Biobank than SAIL (adjusting for age, sex, and socioeconomic status). Both cohorts produced similar hazard ratios for some LTCs (e.g., hypertension and coronary heart disease), but UK Biobank underestimated the risk for others (e.g., alcohol-related disorders or mental health conditions). Hazard ratios for some LTC combinations were similar between the cohorts (e.g., cardiovascular conditions); however, UK Biobank underestimated the risk for combinations including other conditions (e.g., mental health conditions). The main limitations are that SAIL databank represents only part of the UK (Wales only) and that in both cohorts we lacked data on severity of the LTCs included.ConclusionsIn this study, we observed that UK Biobank accurately estimates relative risk of mortality, unscheduled hospitalisation, and MACE associated with LTC counts ≤3. However, for counts ≥4, and for some LTC combinations, estimates of magnitude of association from UK Biobank are likely to be conservative. Researchers should be mindful of these limitations of UK Biobank when conducting and interpreting analyses of multimorbidity. Nonetheless, the richness of data available in UK Biobank does offers opportunities to better understand multimorbidity, particularly where complementary data sources less susceptible to selection bias can be used to inform and qualify analyses of UK Biobank.

Peter Hanlon and colleagues compare the associations between multimorbidity and adverse health outcomes in UK Biobank and the SAIL Databank.     

Integrating comprehensive functional annotations to boost power and accuracy in gene-based association analysis

Gene-based association tests aggregate genotypes across multiple variants for each gene, providing an interpretable gene-level analysis framework for genome-wide association studies (GWAS). Early gene-based test applications often focused on rare coding variants; a more recent wave of gene-based methods, e.g. TWAS, use eQTLs to interrogate regulatory associations. Regulatory variants are expected to be particularly valuable for gene-based analysis, since most GWAS associations to date are non-coding. However, identifying causal genes from regulatory associations remains challenging and contentious. Here, we present a statistical framework and computational tool to integrate heterogeneous annotations with GWAS summary statistics for gene-based analysis, applied with comprehensive coding and tissue-specific regulatory annotations. We compare power and accuracy identifying causal genes across single-annotation, omnibus, and annotation-agnostic gene-based tests in simulation studies and an analysis of 128 traits from the UK Biobank, and find that incorporating heterogeneous annotations in gene-based association analysis increases power and performance identifying causal genes.     

Causal inference for heritable phenotypic risk factors using heterogeneous genetic instruments

Over a decade of genome-wide association studies (GWAS) have led to the finding of extreme polygenicity of complex traits. The phenomenon that “all genes affect every complex trait” complicates Mendelian Randomization (MR) studies, where natural genetic variations are used as instruments to infer the causal effect of heritable risk factors. We reexamine the assumptions of existing MR methods and show how they need to be clarified to allow for pervasive horizontal pleiotropy and heterogeneous effect sizes. We propose a comprehensive framework GRAPPLE to analyze the causal effect of target risk factors with heterogeneous genetic instruments and identify possible pleiotropic patterns from data. By using GWAS summary statistics, GRAPPLE can efficiently use both strong and weak genetic instruments, detect the existence of multiple pleiotropic pathways, determine the causal direction and perform multivariable MR to adjust for confounding risk factors. With GRAPPLE, we analyze the effect of blood lipids, body mass index, and systolic blood pressure on 25 disease outcomes, gaining new information on their causal relationships and potential pleiotropic pathways involved.     

Using genetic variants to evaluate the causal effect of cholesterol lowering on head and neck cancer risk: A Mendelian randomization study

Head and neck squamous cell carcinoma (HNSCC), which includes cancers of the oral cavity and oropharynx, is a cause of substantial global morbidity and mortality. Strategies to reduce disease burden include discovery of novel therapies and repurposing of existing drugs. Statins are commonly prescribed for lowering circulating cholesterol by inhibiting HMG-CoA reductase (HMGCR). Results from some observational studies suggest that statin use may reduce HNSCC risk. We appraised the relationship of genetically-proxied cholesterol-lowering drug targets and other circulating lipid traits with oral (OC) and oropharyngeal (OPC) cancer risk using two-sample Mendelian randomization (MR). For the primary analysis, germline genetic variants in HMGCR, NPC1L1, CETP, PCSK9 and LDLR were used to proxy the effect of low-density lipoprotein cholesterol (LDL-C) lowering therapies. In secondary analyses, variants were used to proxy circulating levels of other lipid traits in a genome-wide association study (GWAS) meta-analysis of 188,578 individuals. Both primary and secondary analyses aimed to estimate the downstream causal effect of cholesterol lowering therapies on OC and OPC risk. The second sample for MR was taken from a GWAS of 6,034 OC and OPC cases and 6,585 controls (GAME-ON). Analyses were replicated in UK Biobank, using 839 OC and OPC cases and 372,016 controls and the results of the GAME-ON and UK Biobank analyses combined in a fixed-effects meta-analysis. We found limited evidence of a causal effect of genetically-proxied LDL-C lowering using HMGCR, NPC1L1, CETP or other circulating lipid traits on either OC or OPC risk. Genetically-proxied PCSK9 inhibition equivalent to a 1 mmol/L (38.7 mg/dL) reduction in LDL-C was associated with an increased risk of OC and OPC combined (OR 1.8 95%CI 1.2, 2.8, p = 9.31 x10^-05), with good concordance between GAME-ON and UK Biobank (I² = 22%). Effects for PCSK9 appeared stronger in relation to OPC (OR 2.6 95%CI 1.4, 4.9) than OC (OR 1.4 95%CI 0.8, 2.4). LDLR variants, resulting in genetically-proxied reduction in LDL-C equivalent to a 1 mmol/L (38.7 mg/dL), reduced the risk of OC and OPC combined (OR 0.7, 95%CI 0.5, 1.0, p = 0.006). A series of pleiotropy-robust and outlier detection methods showed that pleiotropy did not bias our findings. We found limited evidence for a role of cholesterol-lowering in OC and OPC risk, suggesting previous observational results may have been confounded. There was some evidence that genetically-proxied inhibition of PCSK9 increased risk, while lipid-lowering variants in LDLR, reduced risk of combined OC and OPC. This result suggests that the mechanisms of action of PCSK9 on OC and OPC risk may be independent of its cholesterol lowering effects; however, this was not supported uniformly across all sensitivity analyses and further replication of this finding is required.     

Impact of nonrandom selection mechanisms on the causal effect estimation for two-sample Mendelian randomization methods

Nonrandom selection in one-sample Mendelian Randomization (MR) results in biased estimates and inflated type I error rates only when the selection effects are sufficiently large. In two-sample MR, the different selection mechanisms in two samples may more seriously affect the causal effect estimation. Firstly, we propose sufficient conditions for causal effect invariance under different selection mechanisms using two-sample MR methods. In the simulation study, we consider 49 possible selection mechanisms in two-sample MR, which depend on genetic variants (G), exposures (X), outcomes (Y) and their combination. We further compare eight pleiotropy-robust methods under different selection mechanisms. Results of simulation reveal that nonrandom selection in sample II has a larger influence on biases and type I error rates than those in sample I. Furthermore, selections depending on X+Y, G+Y, or G+X+Y in sample II lead to larger biases than other selection mechanisms. Notably, when selection depends on Y, bias of causal estimation for non-zero causal effect is larger than that for null causal effect. Especially, the mode based estimate has the largest standard errors among the eight methods. In the absence of pleiotropy, selections depending on Y or G in sample II show nearly unbiased causal effect estimations when the casual effect is null. In the scenarios of balanced pleiotropy, all eight MR methods, especially MR-Egger, demonstrate large biases because the nonrandom selections result in the violation of the Instrument Strength Independent of Direct Effect (InSIDE) assumption. When directional pleiotropy exists, nonrandom selections have a severe impact on the eight MR methods. Application demonstrates that the nonrandom selection in sample II (coronary heart disease patients) can magnify the causal effect estimation of obesity on HbA1c levels. In conclusion, nonrandom selection in two-sample MR exacerbates the bias of causal effect estimation for pleiotropy-robust MR methods.     

Mendelian randomization mixed-scale treatment effect robust identification and estimation for causal inference

Standard Mendelian randomization (MR) analysis can produce biased results if the genetic variant defining an instrumental variable (IV) is confounded and/or has a horizontal pleiotropic effect on the outcome of interest not mediated by the treatment variable. We provide novel identification conditions for the causal effect of a treatment in the presence of unmeasured confounding by leveraging a possibly invalid IV for which both the IV independence and exclusion restriction assumptions may be violated. The proposed Mendelian randomization mixed-scale treatment effect robust identification (MR MiSTERI) approach relies on (i) an assumption that the treatment effect does not vary with the possibly invalid IV on the additive scale; (ii) that the confounding bias does not vary with the possibly invalid IV on the odds ratio scale; and (iii) that the residual variance for the outcome is heteroskedastic with respect to the possibly invalid IV. Although assumptions (i) and (ii) have, respectively, appeared in the IV literature, assumption (iii) has not; we formally establish that their conjunction can identify a causal effect even with an invalid IV. MR MiSTERI is shown to be particularly advantageous in the presence of pervasive heterogeneity of pleiotropic effects on the additive scale. We propose a simple and consistent three-stage estimator that can be used as a preliminary estimator to a carefully constructed efficient one-step-update estimator. In order to incorporate multiple, possibly correlated, and weak invalid IVs, a common challenge in MR studies, we develop a MAny Weak Invalid Instruments (MR MaWII MiSTERI) approach for strengthened identification and improved estimation accuracy. Both simulation studies and UK Biobank data analysis results demonstrate the robustness of the proposed methods.     

Robust inference of bi-directional causal relationships in presence of correlated pleiotropy with GWAS summary data

To infer a causal relationship between two traits, several correlation-based causal direction (CD) methods have been proposed with the use of SNPs as instrumental variables (IVs) based on GWAS summary data for the two traits; however, none of the existing CD methods can deal with SNPs with correlated pleiotropy. Alternatively, reciprocal Mendelian randomization (MR) can be applied, which however may perform poorly in the presence of (unknown) invalid IVs, especially for bi-directional causal relationships. In this paper, first, we propose a CD method that performs better than existing CD methods regardless of the presence of correlated pleiotropy. Second, along with a simple but yet effective IV screening rule, we propose applying a closely related and state-of-the-art MR method in reciprocal MR, showing its almost identical performance to that of the new CD method when their model assumptions hold; however, if the modeling assumptions are violated, the new CD method is expected to better control type I errors. Notably bi-directional causal relationships impose some unique challenges beyond those for uni-directional ones, and thus requiring special treatments. For example, we point out for the first time several scenarios where a bi-directional relationship, but not a uni-directional one, can unexpectedly cause the violation of some weak modeling assumptions commonly required by many robust MR methods. We also offer some numerical support and a modeling justification for the application of our new methods (and more generally MR) to binary traits. Finally we applied the proposed methods to 12 risk factors and 4 common diseases, confirming mostly well-known uni-directional causal relationships, while identifying some novel and plausible bi-directional ones such as between body mass index and type 2 diabetes (T2D), and between diastolic blood pressure and stroke.     

Genetic associations of adult height with risk of cardioembolic and other subtypes of ischemic stroke: A mendelian randomization study in multiple ancestries

BackgroundTaller adult height is associated with lower risks of ischemic heart disease in mendelian randomization (MR) studies, but little is known about the causal relevance of height for different subtypes of ischemic stroke. The present study examined the causal relevance of height for different subtypes of ischemic stroke.Methods and findingsHeight-associated genetic variants (up to 2,337) from previous genome-wide association studies (GWASs) were used to construct genetic instruments in different ancestral populations. Two-sample MR approaches were used to examine the associations of genetically determined height with ischemic stroke and its subtypes (cardioembolic stroke, large-artery stroke, and small-vessel stroke) in multiple ancestries (the MEGASTROKE consortium, which included genome-wide studies of stroke and stroke subtypes: 60,341 ischemic stroke cases) supported by additional cases in individuals of white British ancestry (UK Biobank [UKB]: 4,055 cases) and Chinese ancestry (China Kadoorie Biobank [CKB]: 10,297 cases). The associations of genetically determined height with established cardiovascular and other risk factors were examined in 336,750 participants from UKB and 58,277 participants from CKB. In MEGASTROKE, genetically determined height was associated with a 4% lower risk (odds ratio [OR] 0.96; 95% confidence interval [CI] 0.94, 0.99; p = 0.007) of ischemic stroke per 1 standard deviation (SD) taller height, but this masked a much stronger positive association of height with cardioembolic stroke (13% higher risk, OR 1.13 [95% CI 1.07, 1.19], p < 0.001) and stronger inverse associations with large-artery stroke (11% lower risk, OR 0.89 [0.84, 0.95], p < 0.001) and small-vessel stroke (13% lower risk, OR 0.87 [0.83, 0.92], p < 0.001). The findings in both UKB and CKB were directionally concordant with those observed in MEGASTROKE, but did not reach statistical significance: For presumed cardioembolic stroke, the ORs were 1.08 (95% CI 0.86, 1.35; p = 0.53) in UKB and 1.20 (0.77, 1.85; p = 0.43) in CKB; for other subtypes of ischemic stroke in UKB, the OR was 0.97 (95% CI 0.90, 1.05; p = 0.49); and for other nonlacunar stroke and lacunar stroke in CKB, the ORs were 0.89 (0.80, 1.00; p = 0.06) and 0.99 (0.88, 1.12; p = 0.85), respectively. In addition, genetically determined height was also positively associated with atrial fibrillation (available only in UKB), and with lean body mass and lung function, and inversely associated with low-density lipoprotein (LDL) cholesterol in both British and Chinese ancestries. Limitations of this study include potential bias from assortative mating or pleiotropic effects of genetic variants and incomplete generalizability of genetic instruments to different populations.ConclusionsThe findings provide support for a causal association of taller adult height with higher risk of cardioembolic stroke and lower risk of other ischemic stroke subtypes in diverse ancestries. Further research is needed to understand the shared biological and physical pathways underlying the associations between height and stroke risks, which could identify potential targets for treatments to prevent stroke.

In a Mendelian randomization study, Andrew B. Linden and colleagues study the relationship between height and risk of stroke subtypes among individuals from the MEGASTROKE consortium, China Kadoorie Biobank, and UK Biobank.     

A Bayesian nonparametric model for zero‐inflated outcomes: Prediction,clustering, and causal estimation

Researchers are often interested in predicting outcomes, detecting distinct subgroups of their data, or estimating causal treatment effects. Pathological data distributions that exhibit skewness and zero‐inflation complicate these tasks—requiring highly flexible, data‐adaptive modeling. In this paper, we present a multipurpose Bayesian nonparametric model for continuous, zero‐inflated outcomes that simultaneously predicts structural zeros, captures skewness, and clusters patients with similar joint data distributions. The flexibility of our approach yields predictions that capture the joint data distribution better than commonly used zero‐inflated methods. Moreover, we demonstrate that our model can be coherently incorporated into a standardization procedure for computing causal effect estimates that are robust to such data pathologies. Uncertainty at all levels of this model flow through to the causal effect estimates of interest—allowing easy point estimation, interval estimation, and posterior predictive checks verifying positivity, a required causal identification assumption. Our simulation results show point estimates to have low bias and interval estimates to have close to nominal coverage under complicated data settings. Under simpler settings, these results hold while incurring lower efficiency loss than comparator methods. We use our proposed method to analyze zero‐inflated inpatient medical costs among endometrial cancer patients receiving either chemotherapy or radiation therapy in the SEER‐Medicare database.     

An empirical comparison of several clustered data approaches under confounding due to cluster effects in the analysis of complications of coronary angioplasty

In the analysis of binary response data from many types of large studies, the data are likely to have arisen from multiple centers, resulting in a within-center correlation for the response. Such correlation, or clustering, occurs when outcomes within centers tend to be more similar to each other than to outcomes in other centers. In studies where there is also variability among centers with respect to the exposure of interest, analysis of the exposure-outcome association may be confounded, even after accounting for within-center correlations. We apply several analytic methods to compare the risk of major complications associated with two strategies, staged and combined procedures, for performing percutaneous transluminal coronary angioplasty (PTCA), a mechanical means of relieving blockage of blood vessels due to atherosclerosis. Combined procedures are used in some centers as a cost-cutting strategy. We performed a number of population-averaged and cluster-specific (conditional) analyses, which (a) make no adjustments for center effects of any kind; (b) make adjustments for the effect of center on only the response; or (c) make adjustments for both the effect of center on the response and the relationship between center and exposure. The method used for this third approach decomposes the procedure type variable into within-center and among-center components, resulting in two odds ratio estimates. The naive analysis, ignoring clusters, gave a highly significant effect of procedure type (OR = 1.6). Population average models gave marginally to very nonsignificant estimates of the OR for treatment type ranging from 1.6 to 1.2 with adjustment only for the effect of centers on response. These results depended on the assumed correlation structure. Conditional (cluster-specific) models and other methods that decomposed the treatment type variable into among- and within-center components all found no within-center effect of procedure type (OR = 1.02, consistently) and a considerable among-center effect. This among-center variability in outcomes was related to the proportion of patients who receive combined procedures and was found even when conditioned on procedure type (within-center) and other patient- and center-level covariates. This example illustrates the importance of addressing the potential for center effects to confound an outcome-exposure association when average exposure varies across clusters. While conditional approaches provide estimates of the within-cluster effect, they do not provide information about among-center effects. We recommend using the decomposition approach, as it provides both types of estimates.     

On cross-ancestry cancer polygenic risk scores

Polygenic risk scores (PRS) can provide useful information for personalized risk stratification and disease risk assessment, especially when combined with non-genetic risk factors. However, their construction depends on the availability of summary statistics from genome-wide association studies (GWAS) independent from the target sample. For best compatibility, it was reported that GWAS and the target sample should match in terms of ancestries. Yet, GWAS, especially in the field of cancer, often lack diversity and are predominated by European ancestry. This bias is a limiting factor in PRS research. By using electronic health records and genetic data from the UK Biobank, we contrast the utility of breast and prostate cancer PRS derived from external European-ancestry-based GWAS across African, East Asian, European, and South Asian ancestry groups. We highlight differences in the PRS distributions of these groups that are amplified when PRS methods condense hundreds of thousands of variants into a single score. While European-GWAS-derived PRS were not directly transferrable across ancestries on an absolute scale, we establish their predictive potential when considering them separately within each group. For example, the top 10% of the breast cancer PRS distributions within each ancestry group each revealed significant enrichments of breast cancer cases compared to the bottom 90% (odds ratio of 2.81 [95%CI: 2.69,2.93] in European, 2.88 [1.85, 4.48] in African, 2.60 [1.25, 5.40] in East Asian, and 2.33 [1.55, 3.51] in South Asian individuals). Our findings highlight a compromise solution for PRS research to compensate for the lack of diversity in well-powered European GWAS efforts while recruitment of diverse participants in the field catches up.     

Assessment of the causal relevance of ECG parameters for risk of atrial fibrillation: A mendelian randomisation study

BackgroundAtrial electrical and structural remodelling in older individuals with cardiovascular risk factors has been associated with changes in surface electrocardiographic (ECG) parameters (e.g., prolongation of the PR interval) and higher risks of atrial fibrillation (AF). However, it has been difficult to establish whether altered ECG parameters are the cause or a consequence of the myocardial substrate leading to AF. This study aimed to examine the potential causal relevance of ECG parameters on risk of AF using mendelian randomisation (MR).Methods and findingsWeighted genetic scores explaining lifelong differences in P-wave duration, PR interval, and QT interval were constructed, and associations between these ECG scores and risk of AF were estimated among 278,792 UK Biobank participants (mean age: 57 years at recruitment; 19,132 AF cases). The independent genetic variants contributing to each of the separate ECG scores, and their corresponding weights, were based on published genome-wide association studies. In UK Biobank, genetic scores representing a 5 ms longer P-wave duration or PR interval were significantly associated with lower risks of AF (odds ratio [OR] 0.91; 95% confidence interval [CI]: 0.87–0.96, P = 2 × 10⁻⁴ and OR 0.94; 95% CI: 0.93–0.96, P = 2 × 10⁻¹⁹, respectively), while longer QT interval was not significantly associated with AF. These effects were independently replicated among a further 17,931 AF cases from the AFGen Consortium. Investigation of potential mechanistic pathways showed that differences in ECG parameters associated with specific ion channel genes had effects on risk of AF consistent with the overall scores, while the overall scores were not associated with changes in left atrial size. Limitations of the study included the inherent assumptions of MR, restriction to individuals of European ancestry, and possible restriction of results to the normal ECG ranges represented in UK Biobank.ConclusionsIn UK Biobank, we observed evidence suggesting a causal relationship between lifelong differences in ECG parameters (particularly PR interval) that reflect longer atrial conduction times and a lower risk of AF. These findings, which appear to be independent of atrial size and concomitant cardiovascular comorbidity, support the relevance of varying mechanisms underpinning AF and indicate that more individualised treatment strategies warrant consideration.

In a Mendelian randomization analysis, Parag Gajendragadkar and colleagues investigate associations between genetically-predicted EEG parameters and risk of atrial fibrillation among UK Biobank participants.     

Vitamin D and COVID-19 susceptibility and severity in the COVID-19 Host Genetics Initiative: A Mendelian randomization study

BackgroundIncreased vitamin D levels, as reflected by 25-hydroxy vitamin D (25OHD) measurements, have been proposed to protect against COVID-19 based on in vitro, observational, and ecological studies. However, vitamin D levels are associated with many confounding variables, and thus associations described to date may not be causal. Vitamin D Mendelian randomization (MR) studies have provided results that are concordant with large-scale vitamin D randomized trials. Here, we used 2-sample MR to assess evidence supporting a causal effect of circulating 25OHD levels on COVID-19 susceptibility and severity.Methods and findingsGenetic variants strongly associated with 25OHD levels in a genome-wide association study (GWAS) of 443,734 participants of European ancestry (including 401,460 from the UK Biobank) were used as instrumental variables. GWASs of COVID-19 susceptibility, hospitalization, and severe disease from the COVID-19 Host Genetics Initiative were used as outcome GWASs. These included up to 14,134 individuals with COVID-19, and up to 1,284,876 without COVID-19, from up to 11 countries. SARS-CoV-2 positivity was determined by laboratory testing or medical chart review. Population controls without COVID-19 were also included in the control groups for all outcomes, including hospitalization and severe disease. Analyses were restricted to individuals of European descent when possible. Using inverse-weighted MR, genetically increased 25OHD levels by 1 standard deviation on the logarithmic scale had no significant association with COVID-19 susceptibility (odds ratio [OR] = 0.95; 95% CI 0.84, 1.08; p = 0.44), hospitalization (OR = 1.09; 95% CI: 0.89, 1.33; p = 0.41), and severe disease (OR = 0.97; 95% CI: 0.77, 1.22; p = 0.77). We used an additional 6 meta-analytic methods, as well as conducting sensitivity analyses after removal of variants at risk of horizontal pleiotropy, and obtained similar results. These results may be limited by weak instrument bias in some analyses. Further, our results do not apply to individuals with vitamin D deficiency.ConclusionsIn this 2-sample MR study, we did not observe evidence to support an association between 25OHD levels and COVID-19 susceptibility, severity, or hospitalization. Hence, vitamin D supplementation as a means of protecting against worsened COVID-19 outcomes is not supported by genetic evidence. Other therapeutic or preventative avenues should be given higher priority for COVID-19 randomized controlled trials.

In a Mendelian randomization analysis, Guillaume Butler-Laporte, Tomoki Nakanishi, and colleages study genetic evidence for a relationship between vitamin D and COVID-19 outcomes.     

Operating characteristics of the rank-based inverse normal transformation for quantitative trait analysis in genome-wide association studies

Quantitative traits analyzed in Genome-Wide Association Studies (GWAS) are often nonnormally distributed. For such traits, association tests based on standard linear regression are subject to reduced power and inflated type I error in finite samples. Applying the rank-based inverse normal transformation (INT) to nonnormally distributed traits has become common practice in GWAS. However, the different variations on INT-based association testing have not been formally defined, and guidance is lacking on when to use which approach. In this paper, we formally define and systematically compare the direct (D-INT) and indirect (I-INT) INT-based association tests. We discuss their assumptions, underlying generative models, and connections. We demonstrate that the relative powers of D-INT and I-INT depend on the underlying data generating process. Since neither approach is uniformly most powerful, we combine them into an adaptive omnibus test (O-INT). O-INT is robust to model misspecification, protects the type I error, and is well powered against a wide range of nonnormally distributed traits. Extensive simulations were conducted to examine the finite sample operating characteristics of these tests. Our results demonstrate that, for nonnormally distributed traits, INT-based tests outperform the standard untransformed association test, both in terms of power and type I error rate control. We apply the proposed methods to GWAS of spirometry traits in the UK Biobank. O-INT has been implemented in the R package RNOmni , which is available on CRAN.     

Overcoming the winner's curse: estimating penetrance parameters from case-control data

Genomewide association studies are now a widely used approach in the search for loci that affect complex traits. After detection of significant association, estimates of penetrance and allele-frequency parameters for the associated variant indicate the importance of that variant and facilitate the planning of replication studies. However, when these estimates are based on the original data used to detect the variant, the results are affected by an ascertainment bias known as the "winner's curse." The actual genetic effect is typically smaller than its estimate. This overestimation of the genetic effect may cause replication studies to fail because the necessary sample size is underestimated. Here, we present an approach that corrects for the ascertainment bias and generates an estimate of the frequency of a variant and its penetrance parameters. The method produces a point estimate and confidence region for the parameter estimates. We study the performance of this method using simulated data sets and show that it is possible to greatly reduce the bias in the parameter estimates, even when the original association study had low power. The uncertainty of the estimate decreases with increasing sample size, independent of the power of the original test for association. Finally, we show that application of the method to case-control data can improve the design of replication studies considerably.     

Causal association between heart failure and bone mineral density: Insights from a two-sample bidirectional Mendelian randomization study

In recent times, the association between HF and BMD has attracted enormous interest in the scientific community. However, published epidemiological observational studies on the relationship between HF and BMD remain inconclusive. Herein, we evaluated from the analytical perspective a two-sample bidirectional MR study to analyze the causal association between HF and BMD using a summary-level GWAS Catalog. To select instrumental SNPs strongly associated with exposure, we took a series of rigorous quality control steps at the time of analysis. The causal MR assessment of HF on the risk of BMD was performed first and then in the opposite direction. To make the conclusions more reliable and robust, the fixed-effects IVW, weighted median-based method, MR–Egger, simple mode and weighted mode were utilized. A maximum likelihood model was also used if necessary. MR–Egger regression, IVW “leave-one-out” sensitivity analysis, MR-PRESSO, MR–Egger intercept test and Cochran's Q statistic methods were used to assess heterogeneity and pleiotropy. Our MR studies supported a meaningful causal association between HF and TB-BMD (IVW: OR = 0.78, 95% CI: 0.68–0.87, p = 0.00588). At the same time, we did not find a significant causal relationship between HF and FA-BMD, FN-BMD or LS-BMD. No significant causal relationships between BMD and HF were observed. This bidirectional MR analysis suggested a causal association of HF with only low TB-BMD, while the reverse causality hypothesis was not found. Studies of the prevention and treatment of total bone mineral density decline in patients with heart failure need to be performed.