Genetic predisposition and antipsychotic treatment effect on metabolic syndrome in schizophrenia: a ten-year follow-up study using the Estonian Biobank.

Background: Schizophrenia (SCZ) patients exhibit 30% higher prevalence of metabolic syndrome (MetS) compared to the general population with its suboptimal management contributing to increased mortality. Large-scale studies providing real-world evidence of the underlying causes remain limited. Methods: To address this gap, we used real-world health data from the Estonian Biobank, spanning a median follow-up of ten years, to investigate the impact of genetic predisposition and antipsychotic treatment on the development of MetS in SCZ patients. Specifically, we set out to characterize antipsychotic treatment patterns, genetic predisposition of MetS traits, MetS prognosis, and body mass index (BMI) trajectories, comparing SCZ cases (n = 677) to age- and sex-matched controls (n = 2708). Findings: SCZ cases exhibited higher genetic predisposition to SCZ (OR = 1.75, 95% CI 1.58-1.94), but lower polygenic burden for increased BMI (OR = 0.88, 95% CI 0.88-0.96) and C-reactive protein (OR = 0.88, 95% CI 0.81-0.97) compared to controls. While SCZ cases showed worse prognosis of MetS (HR 1.95, 95% CI 1.54-2.46), higher antipsychotic adherence within the first treatment year was associated with reduced long-term MetS incidence. Linear mixed modelling, incorporating multiple BMI timepoints, underscored the significant contribution of both, antipsychotic medication, and genetic predisposition to higher BMI, driving the substantially upward trajectory of BMI in SCZ cases. Interpretation: These findings contribute to refining clinical risk prediction and prevention strategies for MetS among SCZ patients and emphasize the significance of incorporating genetic information, long-term patient tracking, and employing diverse perspectives when analyzing real-world health data. Funding: EU Horizon 2020, Swedish Research Council, Estonian Research Council, Estonian Ministry of Education and Research, University of Tartu.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects.

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

Evaluation of noninvasive biospecimens for transcriptome studies.

Transcriptome studies disentangle functional mechanisms of gene expression regulation and may elucidate the underlying biology of disease processes. However, the types of tissues currently collected typically assay a single post-mortem timepoint or are limited to investigating cell types found in blood. Noninvasive tissues may improve disease-relevant discovery by enabling more complex longitudinal study designs, by capturing different and potentially more applicable cell types, and by increasing sample sizes due to reduced collection costs and possible higher enrollment from vulnerable populations. Here, we develop methods for sampling noninvasive biospecimens, investigate their performance across commercial and in-house library preparations, characterize their biology, and assess the feasibility of using noninvasive tissues in a multitude of transcriptomic applications. We collected buccal swabs, hair follicles, saliva, and urine cell pellets from 19 individuals over three to four timepoints, for a total of 300 unique biological samples, which we then prepared with replicates across three library preparations, for a final tally of 472 transcriptomes. Of the four tissues we studied, we found hair follicles and urine cell pellets to be most promising due to the consistency of sample quality, the cell types and expression profiles we observed, and their performance in disease-relevant applications. This is the first study to thoroughly delineate biological and technical features of noninvasive samples and demonstrate their use in a wide array of transcriptomic and clinical analyses. We anticipate future use of these biospecimens will facilitate discovery and development of clinical applications.

The role of depression and antidepressant treatment in antihypertensive medication adherence and persistence: Utilising electronic health record data.

Higher blood pressure levels in patients with depression may be associated with lower adherence to antihypertensive medications (AHMs). Here, we use electronic health record (EHR) data from the Estonian Biobank (EstBB) to investigate the role of lifetime depression in AHM adherence and persistence. We also explore the relationship between antidepressant initiation and intraindividual change in AHM adherence among hypertension (HTN) patients with newly diagnosed depression. Diagnosis and pharmacy refill data were obtained from the National Health Insurance database. Adherence and persistence to AHMs were determined for hypertension (HTN) patients initiating treatment between 2009 and 2017 with a three-year follow-up period. Multivariable regression was used to explore the associations between depression and AHM adherence or persistence, adjusting for sociodemographic, genetic, and health-related factors. A linear mixed-effects model was used to estimate the effect of antidepressant treatment initiation on antihypertensive medication adherence, adjusting for age and sex. We identified 20,724 individuals with newly diagnosed HTN (6294 depression cases and 14,430 controls). Depression was associated with 6% lower probability of AHM adherence (OR = 0.943, 95%CI = 0.909-0.979) and 12% lower odds of AHM persistence (OR = 0.876, 95%CI = 0.821-0.936). Adjusting for sociodemographic, genetic, and health-related factors did not significantly influence these associations. AHM adherence increased 8% six months after initiating antidepressant therapy (N = 132; ß = 0.078; 95%CI = 0.025-0.131). Based on the EHR data on EstBB participants, depression is associated with lower AHM adherence and persistence. Additionally, antidepressant therapy may help improve AHM adherence in patients with depression.

The functional impact of rare variation across the regulatory cascade.

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Multiset correlation and factor analysis enables exploration of multi-omics data.

Multi-omics datasets are becoming more common, necessitating better integration methods to realize their revolutionary potential. Here, we introduce multi-set correlation and factor analysis (MCFA), an unsupervised integration method tailored to the unique challenges of high-dimensional genomics data that enables fast inference of shared and private factors. We used MCFA to integrate methylation markers, protein expression, RNA expression, and metabolite levels in 614 diverse samples from the Trans-Omics for Precision Medicine/Multi-Ethnic Study of Atherosclerosis multi-omics pilot. Samples cluster strongly by ancestry in the shared space, even in the absence of genetic information, while private spaces frequently capture dataset-specific technical variation. Finally, we integrated genetic data by conducting a genome-wide association study (GWAS) of our inferred factors, observing that several factors are enriched for GWAS hits and trans-expression quantitative trait loci. Two of these factors appear to be related to metabolic disease. Our study provides a foundation and framework for further integrative analysis of ever larger multi-modal genomic datasets.

Psychiatric disorders and subsequent risk of cardiovascular disease: a longitudinal matched cohort study across three countries.

Background: Several psychiatric disorders have been associated with increased risk of cardiovascular disease (CVD), however, the role of familial factors and the main disease trajectories remain unknown. Methods: In this longitudinal cohort study, we identified a cohort of 900,240 patients newly diagnosed with psychiatric disorders during January 1, 1987 and December 31, 2016, their 1,002,888 unaffected full siblings, and 1:10 age- and sex-matched reference population from nationwide medical records in Sweden, who had no prior diagnosis of CVD at enrolment. We used flexible parametric models to determine the time-varying association between first-onset psychiatric disorders and incident CVD and CVD death, comparing rates of CVD among patients with psychiatric disorders to the rates of unaffected siblings and matched reference population. We also used disease trajectory analysis to identify main disease trajectories linking psychiatric disorders to CVD. Identified associations and disease trajectories of the Swedish cohort were validated in a similar cohort from nationwide medical records in Denmark (N = 875,634 patients, same criteria during January 1, 1969 and December 31, 2016) and in Estonian cohorts from the Estonian Biobank (N = 30,656 patients, same criteria during January 1, 2006 and December 31, 2020), respectively. Findings: During up to 30 years of follow-up of the Swedish cohort, the crude incidence rate of CVD was 9.7, 7.4 and 7.0 per 1000 person-years among patients with psychiatric disorders, their unaffected siblings, and the matched reference population. Compared with their siblings, patients with psychiatric disorders experienced higher rates of CVD during the first year after diagnosis (hazard ratio [HR], 1.88; 95% confidence interval [CI], 1.79-1.98) and thereafter (1.37; 95% CI, 1.34-1.39). Similar rate increases were noted when comparing with the matched reference population. These results were replicated in the Danish cohort. We identified several disease trajectories linking psychiatric disorders to CVD in the Swedish cohort, with or without mediating medical conditions, including a direct link between psychiatric disorders and hypertensive disorder, ischemic heart disease, venous thromboembolism, angina pectoris, and stroke. These trajectories were validated in the Estonian Biobank cohort. Interpretation: Independent of familial factors, patients with psychiatric disorders are at an elevated risk of subsequent CVD, particularly during first year after diagnosis. Increased surveillance and treatment of CVDs and CVD risk factors should be considered as an integral part of clinical management, in order to reduce risk of CVD among patients with psychiatric disorders. Funding: This research was supported by EU Horizon 2020 Research and Innovation Action Grant, European Research Council Consolidator grant, Icelandic Research fund, Swedish Research Council, US NIMH, the Outstanding Clinical Discipline Project of Shanghai Pudong, the Fundamental Research Funds for the Central Universities, and the European Union through the European Regional Development Fund; the Research Council of Norway; the South-East Regional Health Authority, the Stiftelsen Kristian Gerhard Jebsen, and the EEA-RO-NO-2018-0535.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects.

Bulk tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, while context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell type proportions, we demonstrate that cell type iQTLs could be considered as proxies for cell type-specific QTL effects. The interpretation of age iQTLs, however, warrants caution as the moderation effect of age on the genotype and molecular phenotype association may be mediated by changes in cell type composition. Finally, we show that cell type iQTLs contribute to cell type-specific enrichment of diseases that, in combination with additional functional data, may guide future functional studies. Overall, this study highlights iQTLs to gain insights into the context-specificity of regulatory effects.

Pulmonary emphysema subtypes defined by unsupervised machine learning on CT scans.

BACKGROUND: Treatment and preventative advances for chronic obstructive pulmonary disease (COPD) have been slow due, in part, to limited subphenotypes. We tested if unsupervised machine learning on CT images would discover CT emphysema subtypes with distinct characteristics, prognoses and genetic associations. METHODS: New CT emphysema subtypes were identified by unsupervised machine learning on only the texture and location of emphysematous regions on CT scans from 2853 participants in the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS), a COPD case-control study, followed by data reduction. Subtypes were compared with symptoms and physiology among 2949 participants in the population-based Multi-Ethnic Study of Atherosclerosis (MESA) Lung Study and with prognosis among 6658 MESA participants. Associations with genome-wide single-nucleotide-polymorphisms were examined. RESULTS: The algorithm discovered six reproducible (interlearner intraclass correlation coefficient, 0.91-1.00) CT emphysema subtypes. The most common subtype in SPIROMICS, the combined bronchitis-apical subtype, was associated with chronic bronchitis, accelerated lung function decline, hospitalisations, deaths, incident airflow limitation and a gene variant near DRD1, which is implicated in mucin hypersecretion (p=1.1 ×10-8). The second, the diffuse subtype was associated with lower weight, respiratory hospitalisations and deaths, and incident airflow limitation. The third was associated with age only. The fourth and fifth visually resembled combined pulmonary fibrosis emphysema and had distinct symptoms, physiology, prognosis and genetic associations. The sixth visually resembled vanishing lung syndrome. CONCLUSION: Large-scale unsupervised machine learning on CT scans defined six reproducible, familiar CT emphysema subtypes that suggest paths to specific diagnosis and personalised therapies in COPD and pre-COPD.

Genetic regulators of sputum mucin concentration and their associations with COPD phenotypes.

Hyper-secretion and/or hyper-concentration of mucus is a defining feature of multiple obstructive lung diseases, including chronic obstructive pulmonary disease (COPD). Mucus itself is composed of a mixture of water, ions, salt and proteins, of which the gel-forming mucins, MUC5AC and MUC5B, are the most abundant. Recent studies have linked the concentrations of these proteins in sputum to COPD phenotypes, including chronic bronchitis (CB) and acute exacerbations (AE). We sought to determine whether common genetic variants influence sputum mucin concentrations and whether these variants are also associated with COPD phenotypes, specifically CB and AE. We performed a GWAS to identify quantitative trait loci for sputum mucin protein concentration (pQTL) in the Sub-Populations and InteRmediate Outcome Measures in COPD Study (SPIROMICS, n = 708 for total mucin, n = 215 for MUC5AC, MUC5B). Subsequently, we tested for associations of mucin pQTL with CB and AE using regression modeling (n = 822-1300). Replication analysis was conducted using data from COPDGene (n = 5740) and by examining results from the UK Biobank. We identified one genome-wide significant pQTL for MUC5AC (rs75401036) and two for MUC5B (rs140324259, rs10001928). The strongest association for MUC5B, with rs140324259 on chromosome 11, explained 14% of variation in sputum MUC5B. Despite being associated with lower MUC5B, the C allele of rs140324259 conferred increased risk of CB (odds ratio (OR) = 1.42; 95% confidence interval (CI): 1.10-1.80) as well as AE ascertained over three years of follow up (OR = 1.41; 95% CI: 1.02-1.94). Associations between rs140324259 and CB or AE did not replicate in COPDGene. However, in the UK Biobank, rs140324259 was associated with phenotypes that define CB, namely chronic mucus production and cough, again with the C allele conferring increased risk. We conclude that sputum MUC5AC and MUC5B concentrations are associated with common genetic variants, and the top locus for MUC5B may influence COPD phenotypes, in particular CB.

Canonical correlation analysis for multi-omics: Application to cross-cohort analysis.

Integrative approaches that simultaneously model multi-omics data have gained increasing popularity because they provide holistic system biology views of multiple or all components in a biological system of interest. Canonical correlation analysis (CCA) is a correlation-based integrative method designed to extract latent features shared between multiple assays by finding the linear combinations of features-referred to as canonical variables (CVs)-within each assay that achieve maximal across-assay correlation. Although widely acknowledged as a powerful approach for multi-omics data, CCA has not been systematically applied to multi-omics data in large cohort studies, which has only recently become available. Here, we adapted sparse multiple CCA (SMCCA), a widely-used derivative of CCA, to proteomics and methylomics data from the Multi-Ethnic Study of Atherosclerosis (MESA) and Jackson Heart Study (JHS). To tackle challenges encountered when applying SMCCA to MESA and JHS, our adaptations include the incorporation of the Gram-Schmidt (GS) algorithm with SMCCA to improve orthogonality among CVs, and the development of Sparse Supervised Multiple CCA (SSMCCA) to allow supervised integration analysis for more than two assays. Effective application of SMCCA to the two real datasets reveals important findings. Applying our SMCCA-GS to MESA and JHS, we identified strong associations between blood cell counts and protein abundance, suggesting that adjustment of blood cell composition should be considered in protein-based association studies. Importantly, CVs obtained from two independent cohorts also demonstrate transferability across the cohorts. For example, proteomic CVs learned from JHS, when transferred to MESA, explain similar amounts of blood cell count phenotypic variance in MESA, explaining 39.0% ~ 50.0% variation in JHS and 38.9% ~ 49.1% in MESA. Similar transferability was observed for other omics-CV-trait pairs. This suggests that biologically meaningful and cohort-agnostic variation is captured by CVs. We anticipate that applying our SMCCA-GS and SSMCCA on various cohorts would help identify cohort-agnostic biologically meaningful relationships between multi-omics data and phenotypic traits.

Distinct COPD subtypes in former smokers revealed by gene network perturbation analysis.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) varies significantly in symptomatic and physiologic presentation. Identifying disease subtypes from molecular data, collected from easily accessible blood samples, can help stratify patients and guide disease management and treatment. METHODS: Blood gene expression measured by RNA-sequencing in the COPDGene Study was analyzed using a network perturbation analysis method. Each COPD sample was compared against a learned reference gene network to determine the part that is deregulated. Gene deregulation values were used to cluster the disease samples. RESULTS: The discovery set included 617 former smokers from COPDGene. Four distinct gene network subtypes are identified with significant differences in symptoms, exercise capacity and mortality. These clusters do not necessarily correspond with the levels of lung function impairment and are independently validated in two external cohorts: 769 former smokers from COPDGene and 431 former smokers in the Multi-Ethnic Study of Atherosclerosis (MESA). Additionally, we identify several genes that are significantly deregulated across these subtypes, including DSP and GSTM1, which have been previously associated with COPD through genome-wide association study (GWAS). CONCLUSIONS: The identified subtypes differ in mortality and in their clinical and functional characteristics, underlining the need for multi-dimensional assessment potentially supplemented by selected markers of gene expression. The subtypes were consistent across cohorts and could be used for new patient stratification and disease prognosis.

Arsenic Exposure, Blood DNA Methylation, and Cardiovascular Disease.

BACKGROUND: Epigenetic dysregulation has been proposed as a key mechanism for arsenic-related cardiovascular disease (CVD). We evaluated differentially methylated positions (DMPs) as potential mediators on the association between arsenic and CVD. METHODS: Blood DNA methylation was measured in 2321 participants (mean age 56.2, 58.6% women) of the Strong Heart Study, a prospective cohort of American Indians. Urinary arsenic species were measured using high-performance liquid chromatography coupled to inductively coupled plasma mass spectrometry. We identified DMPs that are potential mediators between arsenic and CVD. In a cross-species analysis, we compared those DMPs with differential liver DNA methylation following early-life arsenic exposure in the apoE knockout (apoE-/-) mouse model of atherosclerosis. RESULTS: A total of 20 and 13 DMPs were potential mediators for CVD incidence and mortality, respectively, several of them annotated to genes related to diabetes. Eleven of these DMPs were similarly associated with incident CVD in 3 diverse prospective cohorts (Framingham Heart Study, Women's Health Initiative, and Multi-Ethnic Study of Atherosclerosis). In the mouse model, differentially methylated regions in 20 of those genes and DMPs in 10 genes were associated with arsenic. CONCLUSIONS: Differential DNA methylation might be part of the biological link between arsenic and CVD. The gene functions suggest that diabetes might represent a relevant mechanism for arsenic-related cardiovascular risk in populations with a high burden of diabetes.

Transcription factor regulation of eQTL activity across individuals and tissues.

Tens of thousands of genetic variants associated with gene expression (cis-eQTLs) have been discovered in the human population. These eQTLs are active in various tissues and contexts, but the molecular mechanisms of eQTL variability are poorly understood, hindering our understanding of genetic regulation across biological contexts. Since many eQTLs are believed to act by altering transcription factor (TF) binding affinity, we hypothesized that analyzing eQTL effect size as a function of TF level may allow discovery of mechanisms of eQTL variability. Using GTEx Consortium eQTL data from 49 tissues, we analyzed the interaction between eQTL effect size and TF level across tissues and across individuals within specific tissues and generated a list of 10,098 TF-eQTL interactions across 2,136 genes that are supported by at least two lines of evidence. These TF-eQTLs were enriched for various TF binding measures, supporting with orthogonal evidence that these eQTLs are regulated by the implicated TFs. We also found that our TF-eQTLs tend to overlap genes with gene-by-environment regulatory effects and to colocalize with GWAS loci, implying that our approach can help to elucidate mechanisms of context-specificity and trait associations. Finally, we highlight an interesting example of IKZF1 TF regulation of an APBB1IP gene eQTL that colocalizes with a GWAS signal for blood cell traits. Together, our findings provide candidate TF mechanisms for a large number of eQTLs and offer a generalizable approach for researchers to discover TF regulators of genetic variant effects in additional QTL datasets.

Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression.

Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation.

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

Genetic regulation of spermine oxidase activity and cancer risk: a Mendelian randomization study.

Spermine oxidase (SMOX) catalyzes the oxidation of spermine to spermidine. Observational studies have reported SMOX as a source of reactive oxygen species associated with cancer, implying that inhibition of SMOX could be a target for chemoprevention. Here we test causality of SMOX levels with cancer risk using a Mendelian randomization analysis. We performed a GWAS of spermidine/spermine ratio to identify genetic variants associated with regulation of SMOX activity. Replication analysis was performed in two datasets of SMOX gene expression. We then did a Mendelian randomization analysis by testing the association between the SMOX genetic instrument and neuroblastoma, gastric, lung, breast, prostate, and colorectal cancers using GWAS summary statistics. GWAS of spermidine/spermine ratio identified SMOX locus (P = 1.34 × 10-49) explaining 32% of the variance. The lead SNP rs1741315 was also associated with SMOX gene expression in newborns (P = 8.48 × 10-28) and adults (P = 2.748 × 10-8) explaining 37% and 6% of the variance, respectively. Genetically determined SMOX activity was not associated with neuroblastoma, gastric, lung, breast, prostate nor colorectal cancer (P > 0.05). A PheWAS of rs1741315 did not reveal any relevant associations. Common genetic variation in the SMOX gene was strongly associated with SMOX activity in newborns, and less strongly in adults. Genetic down-regulation of SMOX was not significantly associated with lower odds of neuroblastoma, gastric, lung, breast, prostate and colorectal cancer. These results may inform studies of SMOX inhibition as a target for chemoprevention.

Epigenome-wide association study of mitochondrial genome copy number.

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD.

Integrative approach identifies SLC6A20 and CXCR6 as putative causal genes for the COVID-19 GWAS signal in the 3p21.31 locus.

To date, the locus with the most robust human genetic association to COVID-19 severity is 3p21.31. Here, we integrate genome-scale CRISPR loss-of-function screens and eQTLs in diverse cell types and tissues to pinpoint genes underlying COVID-19 risk. Our findings identify SLC6A20 and CXCR6 as putative causal genes that modulate COVID-19 risk and highlight the usefulness of this integrative approach to bridge the divide between correlational and causal studies of human biology.