Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution

We show that epigenome- and transcriptome-wide association studies (EWAS and TWAS) are prone to significant inflation and bias of test statistics, an unrecognized phenomenon introducing spurious findings if left unaddressed. Neither GWAS-based methodology nor state-of-the-art confounder adjustment methods completely remove bias and inflation. We propose a Bayesian method to control bias and inflation in EWAS and TWAS based on estimation of the empirical null distribution. Using simulations and real data, we demonstrate that our method maximizes power while properly controlling the false positive rate. We illustrate the utility of our method in large-scale EWAS and TWAS meta-analyses of age and smoking

Thrombus Migration Paradox in Patients With Acute Ischemic Stroke

Background and Purpose- The location of the thrombus as observed on first digital subtraction angiography during endovascular treatment may differ from the initial observation on initial noninvasive imaging. We studied the incidence of thrombus dynamics, its impact on patient outcomes, and its association with intravenous thrombolytics. Methods- We included patients from the MR CLEAN registry (Multicenter Randomized Clinical Trial of Endovascular Treatment of Acute Ischemic Stroke) with an initial target occlusion on computed tomography angiography located in the intracranial internal carotid artery, M1, or M2. The conventional angiography target occlusion was defined during endovascular treatment. Thrombus dynamics were classified as growth, stability, migration, and resolution. The primary outcome was functional outcome at 90 days (modified Rankin Scale). The secondary outcomes were successful and complete reperfusion (extended treatment in cerebral infarction scores of 2b-3 and 3, respectively). Results- The analysis included 1349 patients. Thrombus migration occurred in 302 (22%) patients, thrombus growth in 87 (6%), and thrombus resolution in 39 (3%). Intravenous treatment with alteplase was associated with more thrombus migration (adjusted odds ratio, 2.01; CI, 1.29-3.11) and thrombus resolution (adjusted odds ratio, 1.85; CI, 1.22-2.80). Thrombus migration was associated with a lower chance of complete reperfusion (adjusted odds ratio, 0.57; CI, 0.42-0.78) and successful reperfusion (adjusted odds ratio, 0.74; CI, 0.55-0.99). In the subgroup of patients with M1 initial target occlusion, thrombus migration was associated with better functional outcome (adjusted common odds ratio, 1.49; CI, 1.02-2.17), and there was a trend towards better functional outcome in patients with thrombus resolution (adjusted common odds ratio, 2.23; CI, 0.93-5.37). Conclusions- In patients with acute ischemic stroke, thrombus location regularly changes between computed tomography angiography and digital subtraction angiography. Administration of intravenous alteplase increases the chance of thrombus migration and resolution. Thrombus migration is associated with better functional outcome but reduces the rate of complete reperfusion

Effects of Pharmacogenetic Screening for CYP2D6 Among Elderly Starting Therapy With Nortriptyline or Venlafaxine: A Pragmatic Randomized Controlled Trial (CYSCE Trial)

PURPOSE/BACKGROUND: The duration of untreated depression is a predictor for poor future prognosis, making rapid dose finding essential. Genetic variation of the CYP2D6 isoenzyme can influence the optimal dosage needed for individual patients. The aim of this study was to determine the effectiveness of CYP2D6 pharmacogenetic screening to accelerate drug dosing in older patients with depression initiating nortriptyline or venlafaxine. METHODS/PROCEDURES: In this randomized controlled trial, patients were randomly allocated to one of the study arms. In the intervention arm (DG-I), the specific genotype accompanied by a standardized dosing recommendation based on the patients' genotype and the prescribed drug was directly communicated to the physician of the participant. In both the deviating genotype control arm (DG-C) and the nonrandomized control arm, the physician of the participants was not informed about the genotype and the associated dosing advise. The primary outcome was the time needed to reach adequate drug levels: (1) blood levels within the therapeutic range and (2) no dose adjustments within the previous 3 weeks. FINDINGS/RESULTS: No significant difference was observed in mean time to reach adequate dose or time to adequate dose between DG-I and DG-C. Compared with the nonrandomized control arm group, adequate drug levels were reached significantly faster in the DG-I group (log-rank test; P = 0.004), and there was a similar nonsignificant trend for the DG-C group (log-rank test; P = 0.087). IMPLICATIONS/CONCLUSIONS: The results of this study do not support pharmacogenetic CYP2D6 screening to accelerate dose adjustment for nortriptyline and venlafaxine in older patients with depression

Automated workflow-based exploitation of pathway databases provides new insights into genetic associations of metabolite profiles

Background: Genome-wide association studies (GWAS) have identified many common single nucleotide polymorphisms (SNPs) that associate with clinical phenotypes, but these SNPs usually explain just a small part of the heritability and have relatively modest effect sizes. In contrast, SNPs that associate with metabolite levels generally explain a higher percentage of the genetic variation and demonstrate larger effect sizes. Still, the discovery of SNPs associated with metabolite levels is challenging since testing all metabolites measured in typical metabolomics studies with all SNPs comes with a severe multiple testing penalty. We have developed an automated workflow approach that utilizes prior knowledge of biochemical pathways present in databases like KEGG and BioCyc to generate a smaller SNP set relevant to the metabolite. This paper explores the opportunities and challenges in the analysis of GWAS of metabolomic phenotypes and provides novel insights into the genetic basis of metabolic variation through the re-analysis of published GWAS datasets. Results: Re-analysis of the published GWAS dataset from Illig et al. (Nature Genetics, 2010) using a pathway-based workflow (http://www.myexperiment.org/packs/319.html), confirmed previously identified hits and identified a new locus of human metabolic individuality, associating Aldehyde dehydrogenase family1 L1 (ALDH1L1) with serine/glycine ratios in blood. Replication in an independent GWAS dataset of phospholipids (Demirkan et al., PLoS Genetics, 2012) identified two novel loci supported by additional literature evidence: GPAM (Glycerol-3 phosphate acyltransferase) and CBS (Cystathionine beta-synthase). In addition, the workflow approach provided novel insight into the affected pathways and relevance of some of these gene-metabolite pairs in disease development and progression. Conclusions: We demonstrate the utility of automated exploitation of background knowledge present in pathway databases for the analysis of GWAS datasets of metabolomic phenotypes. We report novel loci and potential biochemical mechanisms that contribute to our understanding of the genetic basis of metabolic variation and its relationship to disease development and progression

Natural History of MYH7-Related Dilated Cardiomyopathy

BACKGROUND Variants in myosin heavy chain 7 (MYH7) are responsible for disease in 1% to 5% of patients with dilated cardiomyopathy (DCM); however, the clinical characteristics and natural history of MYH7-related DCM are poorly described. OBJECTIVES We sought to determine the phenotype and prognosis of MYH7-related DCM. We also evaluated the influence of variant location on phenotypic expression. METHODS We studied clinical data from 147 individuals with DCM-causing MYH7 variants (47.6% female; 35.6 +/- 19.2 years) recruited from 29 international centers. RESULTS At initial evaluation, 106 (72.1%) patients had DCM (left ventricular ejection fraction: 34.5% +/- 11.7%). Median follow-up was 4.5 years (IQR: 1.7-8.0 years), and 23.7% of carriers who were initially phenotype-negative developed DCM. Phenotypic expression by 40 and 60 years was 46% and 88%, respectively, with 18 patients (16%) first diagnosed at <18 years of age. Thirty-six percent of patients with DCM met imaging criteria for LV noncompaction. During follow-up, 28% showed left ventricular reverse remodeling. Incidence of adverse cardiac events among patients with DCM at 5 years was 11.6%, with 5 (4.6%) deaths caused by end-stage heart failure (ESHF) and 5 patients (4.6%) requiring heart transplantation. The major ventricular arrhythmia rate was low (1.0% and 2.1% at 5 years in patients with DCM and in those with LVEF of <= 35%, respectively). ESHF and major ventricular arrhythmia were significantly lower compared with LMNA-related DCM and similar to DCM caused by TTN truncating variants. CONCLUSIONS MYH7-related DCM is characterized by early age of onset, high phenotypic expression, low left ventricular reverse remodeling, and frequent progression to ESHF. Heart failure complications predominate over ventricular arrhythmias, which are rare. (C) 2022 The Authors. Published by Elsevier on behalf of the American College of Cardiology Foundation

The Prospective Dutch Colorectal Cancer (PLCRC) cohort: real-world data facilitating research and clinical care

Real-world data (RWD) sources are important to advance clinical oncology research and evaluate treatments in daily practice. Since 2013, the Prospective Dutch Colorectal Cancer (PLCRC) cohort, linked to the Netherlands Cancer Registry, serves as an infrastructure for scientific research collecting additional patient-reported outcomes (PRO) and biospecimens. Here we report on cohort developments and investigate to what extent PLCRC reflects the “real-world”. Clinical and demographic characteristics of PLCRC participants were compared with the general Dutch CRC population (n = 74,692, Dutch-ref). To study representativeness, standardized differences between PLCRC and Dutch-ref were calculated, and logistic regression models were evaluated on their ability to distinguish cohort participants from the Dutch-ref (AU-ROC 0.5 = preferred, implying participation independent of patient characteristics). Stratified analyses by stage and time-period (2013–2016 and 2017–Aug 2019) were performed to study the evolution towards RWD. In August 2019, 5744 patients were enrolled. Enrollment increased steeply, from 129 participants (1 hospital) in 2013 to 2136 (50 of 75 Dutch hospitals) in 2018. Low AU-ROC (0.65, 95% CI: 0.64–0.65) indicates limited ability to distinguish cohort participants from the Dutch-ref. Characteristics that remained imbalanced in the period 2017–Aug’19 compared with the Dutch-ref were age (65.0 years in PLCRC, 69.3 in the Dutch-ref) and tumor stage (40% stage-III in PLCRC, 30% in the Dutch-ref). PLCRC approaches to represent the Dutch CRC population and will ultimately meet the current demand for high-quality RWD. Efforts are ongoing to improve multidisciplinary recruitment which will further enhance PLCRC’s representativeness and its contribution to a learning healthcare system

Skewed X-inactivation is common in the general female population

X-inactivation is a well-established dosage compensation mechanism ensuring that X-chromosomal genes are expressed at comparable levels in males and females. Skewed X-inactivation is often explained by negative selection of one of the alleles. We demonstrate that imbalanced expression of the paternal and maternal X-chromosomes is common in the general population and that the random nature of the X-inactivation mechanism can be sufficient to explain the imbalance. To this end, we analyzed blood-derived RNA and whole-genome sequencing data from 79 female children and their parents from the Genome of the Netherlands project. We calculated the median ratio of the paternal over total counts at all X-chromosomal heterozygous single-nucleotide variants with coverage ≥10. We identified two individuals where the same X-chromosome was inactivated in all cells. Imbalanced expression of the two X-chromosomes (ratios ≤0.35 or ≥0.65) was observed in nearly 50% of the population. The empirically observed skewing is explained by a theoretical model where X-inactivation takes place in an embryonic stage in which eight cells give rise to the hematopoietic compartment. Genes escaping X-inactivation are expressed from both alleles and therefore demonstrate less skewing than inactivated genes. Using this characteristic, we identified three novel escapee genes (SSR4, REPS2, and SEPT6), but did not find support for many previously reported escapee genes in blood. Our collective data suggest that skewed X-inactivation is common in the general population. This may contribute to manifestation of symptoms in carriers of recessive X-linked disorders. We recommend that X-inactivation results should not be used lightly in the interpretation of X-linked variants

An integrative cross-omics analysis of DNA methylation sites of glucose and insulin homeostasis

Despite existing reports on differential DNA methylation in type 2 diabetes (T2D) and obesity, our understanding of its functional relevance remains limited. Here we show the effect of differential methylation in the early phases of T2D pathology by a blood-based epigenome-wide association study of 4808 non-diabetic Europeans in the discovery phase and 11,750 individuals in the replication. We identify CpGs in LETM1, RBM20, IRS2, MAN2A2 and the 1q25.3 region associated with fasting insulin, and in FCRL6, SLAMF1, APOBEC3H and the 15q26.1 region with fasting glucose. In silico cross-omics analyses highlight the role of differential methylation in the crosstalk between the adaptive immune system and glucose homeostasis. The differential methylation explains at least 16.9% of the association between obesity and insulin. Our study sheds light on the biological interactions between genetic variants driving differential methylation and gene expression in the early pathogenesis of T2D

Plasma lipid profiles discriminate bacterial from viral infection in febrile children

Fever is the most common reason that children present to Emergency Departments. Clinical signs and symptoms suggestive of bacterial infection ar

Integrated analysis of environmental and genetic influences on cord blood DNA methylation in new-borns

Epigenetic processes, including DNA methylation (DNAm), are among the mechanisms allowing integration of genetic and environmental factors to shape cellular function. While many studies have investigated either environmental or genetic contributions to DNAm, few have assessed their integrated effects. Here we examine the relative contributions of prenatal environmental factors and genotype on DNA methylation in neonatal blood at variably methylated regions (VMRs) in 4 independent cohorts (overall n = 2365). We use Akaike’s information criterion to test which factors best explain variability of methylation in the cohort-specific VMRs: several prenatal environmental factors (E), genotypes in cis (G), or their additive (G + E) or interaction (GxE) effects. Genetic and environmental factors in combination best explain DNAm at the majority of VMRs. The CpGs best explained by either G, G + E or GxE are functionally distinct. The enrichment of genetic variants from GxE models in GWAS for complex disorders supports their importance for disease risk