Assessing Digital Phenotyping to Enhance Genetic Studies of Human Diseases

Ancestry-independent risk of venous thromboembolism in individuals with sickle cell trait vs factor V Leiden

ABSTRACT

Sickle cell trait (SCT) is a risk factor for venous thromboembolism (VTE). Prior studies investigating the association between SCT and VTE have been performed nearly exclusively in Black populations. However, race-based research can contribute to systemic racism in medicine. We leveraged data from the 23andMe research cohort (4 184 082 participants) to calculate the ancestry-independent risk of VTE associated with SCT as well as comparative risk estimates for heterozygous factor V Leiden (FVL). Odds ratios (ORs) were calculated using a meta-analysis of 3 genetic ancestry groups (European [n = 3 183 142], Latine [n = 597 539], and African [n = 202 281]) and a secondary full-cohort analysis including 2 additional groups (East Asian [n = 159 863] and South Asian [n = 41 257]). Among the full cohort, 94 323 participants (2.25%) reported a history of VTE. On meta-analysis, individuals with SCT had a 1.45-fold (confidence interval [CI], 1.32-1.60) increased risk of VTE compared with SCT noncarriers, which was similar to the full-cohort estimate. The risk of pulmonary embolism (PE) in SCT (OR, 1.95; CI, 1.72-2.20) was higher than that of isolated deep venous thrombosis (DVT; OR, 1.04; CI, 0.90-1.21). FVL carriers had 3.30-fold (CI, 3.24-3.37) increased risk of VTE compared with FVL noncarriers, with a higher risk of isolated DVT (OR, 3.59; CI, 3.51-3.68) than PE (OR, 2.72; CI, 2.64-2.81). In this large, diverse cohort, the risk of VTE was increased among individuals with SCT compared with those without, independent of race or genetic ancestry. The risk of VTE with SCT was lower than that observed in FVL; however, the pattern of VTE in SCT was PE predominant, which is the opposite to that observed in FVL.

A genome-wide Association study of the Count of Codeine prescriptions

Opioid prescription records in existing electronic health record (EHR) databases are a potentially useful, high-fidelity data source for opioid use-related risk phenotyping in genetic analyses. Prescriptions for codeine derived from EHR records were used as targeting traits by screening 16 million patient-level medication records. Genome-wide association analyses were then conducted to identify genomic loci and candidate genes associated with different count patterns of codeine prescriptions. Both low- and high-prescription counts were captured by developing 8 types of phenotypes with selected ranges of prescription numbers to reflect potentially different levels of opioid risk severity. We identified one significant locus associated with low-count codeine prescriptions (1, 2 or 3 prescriptions), while up to 7 loci were identified for higher counts (≥ 4, ≥ 5, ≥6, or ≥ 7 prescriptions), with a strong overlap across different thresholds. We identified 9 significant genomic loci with all-count phenotype. Further, using the polygenic risk approach, we identified a significant correlation (Tau = 0.67, p = 0.01) between an externally derived polygenic risk score for opioid use disorder and numbers of codeine prescriptions. As a proof-of-concept study, our research provides a novel and generalizable phenotyping pipeline for the genomic study of opioid-related risk traits.

Guidance for estimating penetrance of monogenic disease-causing variants in population cohorts

Penetrance is the probability that an individual with a pathogenic genetic variant develops a specific disease. Knowing the penetrance of variants for monogenic disorders is important for counseling of individuals. Until recently, estimates of penetrance have largely relied on affected individuals and their at-risk family members being clinically referred for genetic testing, a 'phenotype-first' approach. This approach substantially overestimates the penetrance of variants because of ascertainment bias. The recent availability of whole-genome sequencing data in individuals from very-large-scale population-based cohorts now allows 'genotype-first' estimates of penetrance for many conditions. Although this type of population-based study can underestimate penetrance owing to recruitment biases, it provides more accurate estimates of penetrance for secondary or incidental findings. Here, we provide guidance for the conduct of penetrance studies to ensure that robust genotypes and phenotypes are used to accurately estimate penetrance of variants and groups of similarly annotated variants from population-based studies.

Multi-ancestry meta-analysis of tobacco use disorder identifies 461 potential risk genes and reveals associations with multiple health outcomes

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviours and although strides have been made using genome-wide association studies to identify risk variants, most variants identified have been for nicotine consumption, rather than TUD. Here we leveraged four US biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records) in 653,790 individuals (495,005 European, 114,420 African American and 44,365 Latin American) and data from UK Biobank (ncombined = 898,680). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviours in children and hundreds of medical outcomes, including HIV infection, heart disease and pain. This work furthers our biological understanding of TUD and establishes electronic health records as a source of phenotypic information for studying the genetics of TUD.

Lifetime risk and genetic predisposition to post-traumatic OA of the knee in the UK Biobank

OBJECTIVE

Acute knee injury is associated with post-traumatic OA (PTOA). Very little is known about the genome-wide associations of PTOA when compared with idiopathic OA (iOA). Our objective was to describe the development of knee OA after a knee injury and its genetic associations in UK Biobank (UKB).

DESIGN

Clinically significant structural knee injuries in those ≤50 years were identified from electronic health records and self-reported data in 502,409 UKB participants. Time-to-first knee osteoarthritis (OA) code was compared in injured cases and age-/sex-matched non-injured controls using Cox Proportional Hazards models. A time-to-OA genome-wide association study (GWAS) sought evidence for PTOA risk variants 6 months to 20 years following injury. Evidence for associations of two iOA polygenic risk scores (PRS) was sought.

RESULTS

Of 4233 knee injury cases, 1896 (44.8%) were female (mean age at injury 34.1 years [SD 10.4]). Over a median of 30.2 (IQR 19.5-45.4) years, 1096 (25.9%) of injured cases developed knee OA. The overall hazards ratio (HR) for knee OA after injury was 1.81 (1.70,1.93), P = 8.9 × 10-74. Female sex and increasing age at injury were associated with knee OA following injury (HR 1.15 [1.02,1.30];1.07 [1,07,1.07] respectively). OA risk was highest in the first 5 years after injury (HR 3.26 [2.67,3.98]), persisting for 40 years. In 3074 knee injury cases included in the time-to-OA GWAS, no variants reached genome-wide significance. iOA PRS was not associated with time-to-OA (HR 0.43 [0.02,8.41]).

CONCLUSIONS

Increasing age at injury and female sex appear to be associated with future development of PTOA in UKB, the risk of which was greatest in the 5 years after injury. Further international efforts towards a better-powered meta-analysis will definitively elucidate genetic similarities and differences of PTOA and iOA.

Multi-ancestry meta-analysis of tobacco use disorder prioritizes novel candidate risk genes and reveals associations with numerous health outcomes

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviors, and although strides have been made using genome-wide association studies (GWAS) to identify risk variants, the majority of variants identified have been for nicotine consumption, rather than TUD. We leveraged five biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records, EHR) in 898,680 individuals (739,895 European, 114,420 African American, 44,365 Latin American). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviors in children, and hundreds of medical outcomes, including HIV infection, heart disease, and pain. This work furthers our biological understanding of TUD and establishes EHR as a source of phenotypic information for studying the genetics of TUD.

Comparative neurogenetics of dog behavior complements efforts towards human neuropsychiatric genetics

Domestic dogs display a wide array of heritable behaviors that have intermediate genetic complexity thanks to a long history of human-influenced selection. Comparative genetics in dogs could address the scarcity of non-human neurogenetic systems relevant to human neuropsychiatric disorders, which are characterized by mental, emotional, and behavioral symptoms and involve vastly complex genetic and non-genetic risk factors. Our review describes the diverse behavioral "phenome" of domestic dogs, past and ongoing sources of behavioral selection, and the state of canine behavioral genetics. We highlight two naturally disordered behavioral domains that illustrate how dogs may prove useful as a comparative, forward neurogenetic system: canine age-related cognitive dysfunction, which can be examined more rapidly given the attenuated lifespan of dogs, and compulsive disorders, which may have genetic roots in purpose-bred behaviors. Growing community science initiatives aimed at the companion dog population will be well suited to investigating such complex behavioral phenotypes and offer a comparative resource that parallels human genomic initiatives in scale and dimensionality.

Explainable multi-task learning improves the parallel estimation of polygenic risk scores for many diseases through shared genetic basis

Many complex diseases share common genetic determinants and are comorbid in a population. We hypothesized that the co-occurrences of diseases and their overlapping genetic etiology can be exploited to simultaneously improve multiple diseases' polygenic risk scores (PRS). This hypothesis was tested using a multi-task learning (MTL) approach based on an explainable neural network architecture. We found that parallel estimations of the PRS for 17 prevalent cancers in a pan-cancer MTL model were generally more accurate than independent estimations for individual cancers in comparable single-task learning (STL) models. Such performance improvement conferred by positive transfer learning was also observed consistently for 60 prevalent non-cancer diseases in a pan-disease MTL model. Interpretation of the MTL models revealed significant genetic correlations between the important sets of single nucleotide polymorphisms used by the neural network for PRS estimation. This suggested a well-connected network of diseases with shared genetic basis.

Leveraging family history data to disentangle time-varying effects on disease risk using lifecourse mendelian randomization

Lifecourse Mendelian randomization is a causal inference technique which harnesses genetic variants with time-varying effects to develop insight into the influence of age-dependent lifestyle factors on disease risk. Here, we apply this approach to evaluate whether childhood body size has a direct consequence on 8 major disease endpoints by analysing parental history data from the UK Biobank study.Our findings suggest that, whilst childhood body size increases later risk of outcomes such as heart disease (odds ratio (OR) = 1.15, 95% CI = 1.07 to 1.23, P = 7.8 × 10- 5) and diabetes (OR = 1.43, 95% CI = 1.31 to 1.56, P = 9.4 × 10- 15) based on parental history data, these findings are likely attributed to a sustained influence of being overweight for many years over the lifecourse. Likewise, we found evidence that remaining overweight throughout the lifecourse increases risk of lung cancer, which was partially mediated by lifetime smoking index. In contrast, using parental history data provided evidence that being overweight in childhood may have a protective effect on risk of breast cancer (OR = 0.87, 95% CI = 0.78 to 0.97, P = 0.01), corroborating findings from observational studies and large-scale genetic consortia.Large-scale family disease history data can provide a complementary source of evidence for epidemiological studies to exploit, particularly given that they are likely more robust to sources of selection bias (e.g. survival bias) compared to conventional case control studies. Leveraging these data using approaches such as lifecourse Mendelian randomization can help elucidate additional layers of evidence to dissect age-dependent effects on disease risk.

Disentangling the aetiological pathways between body mass index and site-specific cancer risk using tissue-partitioned Mendelian randomisation

BACKGROUND

Body mass index (BMI) is known to influence the risk of various site-specific cancers, however, dissecting which subcomponents of this heterogenous risk factor are predominantly responsible for driving disease effects has proven difficult to establish. We have leveraged tissue-specific gene expression to separate the effects of distinct phenotypes underlying BMI on the risk of seven site-specific cancers.

METHODS

SNP-exposure estimates were weighted in a multivariable Mendelian randomisation analysis by their evidence for colocalization with subcutaneous adipose- and brain-tissue-derived gene expression using a recently developed methodology.

RESULTS

Our results provide evidence that brain-tissue-derived BMI variants are predominantly responsible for driving the genetically predicted effect of BMI on lung cancer (OR: 1.17; 95% CI: 1.01-1.36; P = 0.03). Similar findings were identified when analysing cigarettes per day as an outcome (Beta = 0.44; 95% CI: 0.26-0.61; P = 1.62 × 10-6), highlighting a possible shared aetiology or mediator effect between brain-tissue BMI, smoking and lung cancer. Our results additionally suggest that adipose-tissue-derived BMI variants may predominantly drive the effect of BMI and increased risk for endometrial cancer (OR: 1.71; 95% CI: 1.07-2.74; P = 0.02), highlighting a putatively important role in the aetiology of endometrial cancer.

CONCLUSIONS

The study provides valuable insight into the divergent underlying pathways between BMI and the risk of site-specific cancers.

Exploring the Potential of Apple SensorKit and Digital Phenotyping Data as New Digital Biomarkers for Mental Health Research

The use of digital phenotyping continues to expand across all fields of health. By collecting quantitative data in real-time using devices such as smartphones or smartwatches, researchers and clinicians can develop a profile of a wide range of conditions. Smartphones contain sensors that collect data, such as GPS or accelerometer data, which can inform secondary metrics such as time spent at home, location entropy, or even sleep duration. These metrics, when used as digital biomarkers, are not only used to investigate the relationship between behavior and health symptoms but can also be used to support personalized and preventative care. Successful phenotyping requires consistent long-term collection of relevant and high-quality data. In this paper, we present the potential of newly available, for approved research, opt-in SensorKit sensors on iOS devices in improving the accuracy of digital phenotyping. We collected opt-in sensor data over 1 week from a single person with depression using the open-source mindLAMP app developed by the Division of Digital Psychiatry at Beth Israel Deaconess Medical Center. Five sensors from SensorKit were included. The names of the sensors, as listed in official documentation, include the following: phone usage, messages usage, visits, device usage, and ambient light. We compared data from these five new sensors from SensorKit to our current digital phenotyping data collection sensors to assess similarity and differences in both raw and processed data. We present sample data from all five of these new sensors. We also present sample data from current digital phenotyping sources and compare these data to SensorKit sensors when applicable. SensorKit offers great potential for health research. Many SensorKit sensors improve upon previously accessible features and produce data that appears clinically relevant. SensorKit sensors will likely play a substantial role in digital phenotyping. However, using these data requires advanced health app infrastructure and the ability to securely store high-frequency data.

Inverse Association of APOE ε4 and Glaucoma Modified by Systemic Hypertension: The Canadian Longitudinal Study on Aging

Purpose

Studies examining the apolipoprotein E (APOE) ε4 allele and glaucoma are inconsistent, which could be due to interactions with other factors. We examined the relationship between the APOE ε4 allele and glaucoma and intraocular pressure in a large, population-based random sample and explored whether the APOE ε4 allele interacted with systemic hypertension.

Methods

Data came from the Canadian Longitudinal Study on Aging, a population-based study that included 24,655 adults ages 45 to 85 years old in the European ancestry cohort. APOE genotypes were derived from single-nucleotide polymorphisms rs429358 and rs7412. Participants were asked about a prior diagnosis of glaucoma from a doctor. Corneal compensated intraocular pressure (IOP) was measured using the Reichart Ocular Response Analyzer.

Results

Having an APOE ε4 allele was associated with a lower odds of glaucoma after adjusting for age, sex, IOP, and the top 10 population structure principal components (odds ratio [OR] = 0.83; 95% confidence interval [CI], 0.69-0.98; P = 0.033). A novel statistically significant interaction was found in that having an APOE ε4 allele was only associated with glaucoma in those without systemic hypertension (OR = 0.62; 95% CI, 0.46-0.85) although it was not associated in those with it (OR = 0.97; 95% CI, 0.79-1.21) (interaction term P value = 0.017). APOE ε4 was not associated with IOP (β = -0.01; 95% CI, -0.13 to 0.10).

Conclusions

Evidence increasingly points to the APOE ε4 allele having protective benefits against glaucoma, but this association was limited to those without systemic hypertension. Further research is needed to understand the biological mechanisms for these findings and the treatment potential they hold.

Genetics of irritable bowel syndrome: shifting gear via biobank-scale studies

The pathophysiology of irritable bowel syndrome (IBS) is multifactorial and probably involves genetic predisposition and the effect of environmental factors. Unlike other gastrointestinal diseases with a heritable component, genetic research in IBS has been scarce and mostly characterized by small underpowered studies, leading to inconclusive results. The availability of genomic and health-related data from large international cohorts and population-based biobanks offers unprecedented opportunities for long-awaited, well-powered genetic studies in IBS. This Review focuses on the latest advances that provide compelling evidence for the importance of genes involved in the digestion of carbohydrates, ion channel function, neurotransmitters and their receptors, neuronal pathways and the control of gut motility. These discoveries have generated novel information that might be further refined for the identification of predisposed individuals and selection of management strategies for patients. This Review presents a conceptual framework, the advantages and potential limitations of modern genetic research in IBS, and a summary of available evidence.

Obesity-related biomarkers underlie a shared genetic architecture between childhood body mass index and childhood asthma

Obesity and asthma are both common diseases with high population burden worldwide. Recent genetic association studies have shown that obesity is associated with asthma in adults. The relationship between childhood obesity and childhood asthma, and the underlying mechanisms linking obesity to asthma remain to be clarified. In the present study, leveraging large-scale genetic data from UK biobank and several other data sources, we investigated the shared genetic components between body mass index (BMI, n = 39620) in children and childhood asthma (n_case = 10524, n_control = 373393). We included GWAS summary statistics for nine obesity-related biomarkers to evaluate potential biological mediators underlying obesity and asthma. We found a genetic correlation (Rg = 0.10, P = 0.02) between childhood BMI and childhood asthma, whereas the genetic correlation between adult BMI (n = 371541) and childhood asthma was null (Rg = -0.03, P = 0.21). Genomic structural equation modeling analysis further provided evidence that the genetic effect of childhood BMI on childhood asthma (standardized effect size 0.17, P = 0.009) was not driven by the genetic component of adult BMI. Bayesian colocalization analysis identified a shared causal variant rs12436181 that was mapped to gene AMN using gene expression data in lung tissue. Mendelian randomization showed that the odds ratio of childhood asthma for one standard deviation higher of childhood BMI was 1.13 (95% confidence interval: 0.96-1.34). A systematic survey of obesity-related biomarkers showed that IL-6 and adiponectin are potential biological mediators linking obesity and asthma in children. This large-scale genetic study provides evidence that unique childhood obesity pathways could lead to childhood asthma. The findings shed light on childhood asthma pathogenic mechanisms and prevention.

Novel digital approaches to the assessment of problematic opioid use

The opioid epidemic continues to contribute to loss of life through overdose and significant social and economic burdens. Many individuals who develop problematic opioid use (POU) do so after being exposed to prescribed opioid analgesics. Therefore, it is important to accurately identify and classify risk factors for POU. In this review, we discuss the etiology of POU and highlight novel approaches to identifying its risk factors. These approaches include the application of polygenic risk scores (PRS) and diverse machine learning (ML) algorithms used in tandem with data from electronic health records (EHR), clinical notes, patient demographics, and digital footprints. The implementation and synergy of these types of data and approaches can greatly assist in reducing the incidence of POU and opioid-related mortality by increasing the knowledge base of patient-related risk factors, which can help to improve prescribing practices for opioid analgesics.

Digital tools for the assessment of pharmacological treatment for depressive disorder: State of the art

Depression is an invalidating disorder, marked by phenotypic heterogeneity. Clinical assessments for treatment adjustments and data-collection for pharmacological research often rely on subjective representations of functioning. Better phenotyping through digital applications may add unseen information and facilitate disentangling the clinical characteristics and impact of depression and its pharmacological treatment in everyday life. Researchers, physicians, and patients benefit from well-understood digital phenotyping approaches to assess the treatment efficacy and side-effects. This review discusses the current possibilities and pitfalls of wearables and technology for the assessment of the pharmacological treatment of depression. Their applications in the whole spectrum of treatment for depression, including diagnosis, treatment of an episode, and monitoring of relapse risk and prevention are discussed. Multiple aspects are to be considered, including concerns that come with collecting sensitive data and health recordings. Also, privacy and trust are addressed. Available applications range from questionnaire-like apps to objective assessment of behavioural patterns and promises in handling suicidality. Nonetheless, interpretation and integration of this high-resolution information with other phenotyping levels, remains challenging. This review provides a state-of-the-art description of wearables and technology in digital phenotyping for monitoring pharmacological treatment in depression, focusing on the challenges and opportunities of its application in clinical trials and research.

Environmental variables and genome-environment interactions predicting IBD diagnosis in large UK cohort

A combination of genetic susceptibility and environmental exposure is thought to cause inflammatory bowel disease (IBD), but the non-genetic component remains poorly characterized. We therefore undertook a search for environmental variables and gene-environment interactions associated with future IBD diagnosis in a large UK cohort. Using self-report and electronic health records, we identified 1946 Crohn's disease (CD) and 3715 ulcerative colitis (UC) patients after quality control in the UK Biobank. Based on prior literature and biological plausibility , we tested 38 candidate environmental variables for association with CD, UC, and overall IBD using Cox proportional hazard regressions. We also tested whether these variables interacted with polygenic risk in predicting disease, following up significant (FDR < 0.05) results with tests for SNP-environment associations. We performed robustness analyses on all significant results. As in previous reports, appendectomy protected against UC, smoking (both current and previous) elevated risk for CD, current smoking protected against UC, and previous smoking imparted a risk for UC. Childhood antibiotic use associated with IBD, as did sun exposure during the winter. Socioeconomic deprivation was conferred a risk for IBD, CD, and UC. We uncovered negative interactions between polygenic risk and previous oral contraceptive use for IBD and UC. Polygenic risk also interacted negatively with previous smoking in predicting UC. There were no individually significant SNP-environment interactions. Thus, for a limited set of environmental variables, there was strong evidence of association with IBD diagnosis in the UK Biobank, and interaction with polygenic risk was minimal.

Ancestry-inclusive dog genomics challenges popular breed stereotypes

Behavioral genetics in dogs has focused on modern breeds, which are isolated subgroups with distinctive physical and, purportedly, behavioral characteristics. We interrogated breed stereotypes by surveying owners of 18,385 purebred and mixed-breed dogs and genotyping 2155 dogs. Most behavioral traits are heritable [heritability (h²) > 25%], and admixture patterns in mixed-breed dogs reveal breed propensities. Breed explains just 9% of behavioral variation in individuals. Genome-wide association analyses identify 11 loci that are significantly associated with behavior, and characteristic breed behaviors exhibit genetic complexity. Behavioral loci are not unusually differentiated in breeds, but breed propensities align, albeit weakly, with ancestral function. We propose that behaviors perceived as characteristic of modern breeds derive from thousands of years of polygenic adaptation that predates breed formation, with modern breeds distinguished primarily by aesthetic traits.

Assessing the causal role of epigenetic clocks in the development of multiple cancers: a Mendelian randomization study

Background

Epigenetic clocks have been associated with cancer risk in several observational studies. Nevertheless, it is unclear whether they play a causal role in cancer risk or if they act as a non-causal biomarker.

Methods

We conducted a two-sample Mendelian randomization (MR) study to examine the genetically predicted effects of epigenetic age acceleration as measured by HannumAge (nine single-nucleotide polymorphisms (SNPs)), Horvath Intrinsic Age (24 SNPs), PhenoAge (11 SNPs), and GrimAge (4 SNPs) on multiple cancers (i.e. breast, prostate, colorectal, ovarian and lung cancer). We obtained genome-wide association data for biological ageing from a meta-analysis (N = 34,710), and for cancer from the UK Biobank (N cases = 2671-13,879; N controls = 173,493-372,016), FinnGen (N cases = 719-8401; N controls = 74,685-174,006) and several international cancer genetic consortia (N cases = 11,348-122,977; N controls = 15,861-105,974). Main analyses were performed using multiplicative random effects inverse variance weighted (IVW) MR. Individual study estimates were pooled using fixed effect meta-analysis. Sensitivity analyses included MR-Egger, weighted median, weighted mode and Causal Analysis using Summary Effect Estimates (CAUSE) methods, which are robust to some of the assumptions of the IVW approach.

Results

Meta-analysed IVW MR findings suggested that higher GrimAge acceleration increased the risk of colorectal cancer (OR = 1.12 per year increase in GrimAge acceleration, 95% CI 1.04-1.20, p = 0.002). The direction of the genetically predicted effects was consistent across main and sensitivity MR analyses. Among subtypes, the genetically predicted effect of GrimAge acceleration was greater for colon cancer (IVW OR = 1.15, 95% CI 1.09-1.21, p = 0.006), than rectal cancer (IVW OR = 1.05, 95% CI 0.97-1.13, p = 0.24). Results were less consistent for associations between other epigenetic clocks and cancers.

Conclusions

GrimAge acceleration may increase the risk of colorectal cancer. Findings for other clocks and cancers were inconsistent. Further work is required to investigate the potential mechanisms underlying the results.

Funding

FMB was supported by a Wellcome Trust PhD studentship in Molecular, Genetic and Lifecourse Epidemiology (224982/Z/22/Z which is part of grant 218495/Z/19/Z). KKT was supported by a Cancer Research UK (C18281/A29019) programme grant (the Integrative Cancer Epidemiology Programme) and by the Hellenic Republic's Operational Programme 'Competitiveness, Entrepreneurship & Innovation' (OΠΣ 5047228). PH was supported by Cancer Research UK (C18281/A29019). RMM was supported by the NIHR Biomedical Research Centre at University Hospitals Bristol and Weston NHS Foundation Trust and the University of Bristol and by a Cancer Research UK (C18281/A29019) programme grant (the Integrative Cancer Epidemiology Programme). RMM is a National Institute for Health Research Senior Investigator (NIHR202411). The views expressed are those of the author(s) and not necessarily those of the NIHR or the Department of Health and Social Care. GDS and CLR were supported by the Medical Research Council (MC_UU_00011/1 and MC_UU_00011/5, respectively) and by a Cancer Research UK (C18281/A29019) programme grant (the Integrative Cancer Epidemiology Programme). REM was supported by an Alzheimer's Society project grant (AS-PG-19b-010) and NIH grant (U01 AG-18-018, PI: Steve Horvath). RCR is a de Pass Vice Chancellor's Research Fellow at the University of Bristol.

Significant sparse polygenic risk scores across 813 traits in UK Biobank

We present a systematic assessment of polygenic risk score (PRS) prediction across more than 1,500 traits using genetic and phenotype data in the UK Biobank. We report 813 sparse PRS models with significant (p < 2.5 x 10-5) incremental predictive performance when compared against the covariate-only model that considers age, sex, types of genotyping arrays, and the principal component loadings of genotypes. We report a significant correlation between the number of genetic variants selected in the sparse PRS model and the incremental predictive performance (Spearman's ⍴ = 0.61, p = 2.2 x 10-59 for quantitative traits, ⍴ = 0.21, p = 9.6 x 10-4 for binary traits). The sparse PRS model trained on European individuals showed limited transferability when evaluated on non-European individuals in the UK Biobank. We provide the PRS model weights on the Global Biobank Engine (https://biobankengine.stanford.edu/prs).

Physician-Confirmed and Administrative Definitions of Stroke in UK Biobank Reflect the Same Underlying Genetic Trait

Background

Stroke in UK Biobank (UKB) is ascertained via linkages to coded administrative datasets and self-report. We studied the accuracy of these codes using genetic validation.

Methods

We compiled stroke-specific and broad cerebrovascular disease (CVD) code lists (Read V2/V3, ICD-9/-10) for medical settings (hospital, death record, primary care) and self-report. Among 408,210 UKB participants, we identified all with a relevant code, creating 12 stroke definitions based on the code type and source. We performed genome-wide association studies (GWASs) for each definition, comparing summary results against the largest published stroke GWAS (MEGASTROKE), assessing genetic correlations, and replicating 32 stroke-associated loci.

Results

The stroke case numbers identified varied widely from 3,976 (primary care stroke-specific codes) to 19,449 (all codes, all sources). All 12 UKB stroke definitions were significantly correlated with the MEGASTROKE summary GWAS results (rg.81-1) and each other (rg.4-1). However, Bonferroni-corrected confidence intervals were wide, suggesting limited precision of some results. Six previously reported stroke-associated loci were replicated using ≥1 UKB stroke definition.

Conclusions

Stroke case numbers in UKB depend on the code source and type used, with a 5-fold difference in the maximum case-sample size. All stroke definitions are significantly genetically correlated with the largest stroke GWAS to date.

Genome-wide analysis of 53,400 people with irritable bowel syndrome highlights shared genetic pathways with mood and anxiety disorders

Irritable bowel syndrome (IBS) results from disordered brain-gut interactions. Identifying susceptibility genes could highlight the underlying pathophysiological mechanisms. We designed a digestive health questionnaire for UK Biobank and combined identified cases with IBS with independent cohorts. We conducted a genome-wide association study with 53,400 cases and 433,201 controls and replicated significant associations in a 23andMe panel (205,252 cases and 1,384,055 controls). Our study identified and confirmed six genetic susceptibility loci for IBS. Implicated genes included NCAM1, CADM2, PHF2/FAM120A, DOCK9, CKAP2/TPTE2P3 and BAG6. The first four are associated with mood and anxiety disorders, expressed in the nervous system, or both. Mirroring this, we also found strong genome-wide correlation between the risk of IBS and anxiety, neuroticism and depression (rg > 0.5). Additional analyses suggested this arises due to shared pathogenic pathways rather than, for example, anxiety causing abdominal symptoms. Implicated mechanisms require further exploration to help understand the altered brain-gut interactions underlying IBS.

APOC3 genetic variation, serum triglycerides, and risk of coronary artery disease in Asian Indians, Europeans, and other ethnic groups

BACKGROUND

Hypertriglyceridemia has emerged as a critical coronary artery disease (CAD) risk factor. Rare loss-of-function (LoF) variants in apolipoprotein C-III have been reported to reduce triglycerides (TG) and are cardioprotective in American Indians and Europeans. However, there is a lack of data in other Europeans and non-Europeans. Also, whether genetically increased plasma TG due to ApoC-III is causally associated with increased CAD risk is still unclear and inconsistent. The objectives of this study were to verify the cardioprotective role of earlier reported six LoF variants of APOC3 in South Asians and other multi-ethnic cohorts and to evaluate the causal association of TG raising common variants for increasing CAD risk.

METHODS

We performed gene-centric and Mendelian randomization analyses and evaluated the role of genetic variation encompassing APOC3 for affecting circulating TG and the risk for developing CAD.

RESULTS

One rare LoF variant (rs138326449) with a 37% reduction in TG was associated with lowered risk for CAD in Europeans (p = 0.007), but we could not confirm this association in Asian Indians (p = 0.641). Our data could not validate the cardioprotective role of other five LoF variants analysed. A common variant rs5128 in the APOC3 was strongly associated with elevated TG levels showing a p-value 2.8 × 10- 424. Measures of plasma ApoC-III in a small subset of Sikhs revealed a 37% increase in ApoC-III concentrations among homozygous mutant carriers than the wild-type carriers of rs5128. A genetically instrumented per 1SD increment of plasma TG level of 15 mg/dL would cause a mild increase (3%) in the risk for CAD (p = 0.042).

CONCLUSIONS

Our results highlight the challenges of inclusion of rare variant information in clinical risk assessment and the generalizability of implementation of ApoC-III inhibition for treating atherosclerotic disease. More studies would be needed to confirm whether genetically raised TG and ApoC-III concentrations would increase CAD risk.

Genomewide Association Studies in Pharmacogenomics

The increasing availability of genotype data linked with information about drug-response phenotypes has enabled genomewide association studies (GWAS) that uncover genetic determinants of drug response. GWAS have discovered associations between genetic variants and both drug efficacy and adverse drug reactions. Despite these successes, the design of GWAS in pharmacogenomics (PGx) faces unique challenges. In this review, we analyze the last decade of GWAS in PGx. We review trends in publications over time, including the drugs and drug classes studied and the clinical phenotypes used. Several data sharing consortia have contributed substantially to the PGx GWAS literature. We anticipate increased focus on biobanks and highlight phenotypes that would best enable future PGx discoveries.

Digital Phenotyping in Livestock Farming

Currently, large volumes of data are being collected on farms using multimodal sensor technologies. These sensors measure the activity, housing conditions, feed intake, and health of farm animals. With traditional methods, the data from farm animals and their environment can be collected intermittently. However, with the advancement of wearable and non-invasive sensing tools, these measurements can be made in real-time for continuous quantitation relating to clinical biomarkers, resilience indicators, and behavioral predictors. The digital phenotyping of humans has drawn enormous attention recently due to its medical significance, but much research is still needed for the digital phenotyping of farm animals. Implications from human studies show great promise for the application of digital phenotyping technology in modern livestock farming, but these technologies must be directly applied to animals to understand their true capacities. Due to species-specific traits, certain technologies required to assess phenotypes need to be tailored efficiently and accurately. Such devices allow for the collection of information that can better inform farmers on aspects of animal welfare and production that need improvement. By explicitly addressing farm animals' individual physiological and mental (affective states) needs, sensor-based digital phenotyping has the potential to serve as an effective intervention platform. Future research is warranted for the design and development of digital phenotyping technology platforms that create shared data standards, metrics, and repositories.

Advancing drug discovery using the power of the human genome

Human genetics plays an increasingly important role in drug development and population health. Here we review the history of human genetics in the context of accelerating the discovery of therapies, present examples of how human genetics evidence supports successful drug targets, and discuss how polygenic risk scores could be beneficial in various clinical settings. We highlight the value of direct-to-consumer platforms in the era of fast-paced big data biotechnology, and how diverse genetic and health data can benefit society. © 2021 23andMe, Inc. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Large-scale machine-learning-based phenotyping significantly improves genomic discovery for optic nerve head morphology

Genome-wide association studies (GWASs) require accurate cohort phenotyping, but expert labeling can be costly, time intensive, and variable. Here, we develop a machine learning (ML) model to predict glaucomatous optic nerve head features from color fundus photographs. We used the model to predict vertical cup-to-disc ratio (VCDR), a diagnostic parameter and cardinal endophenotype for glaucoma, in 65,680 Europeans in the UK Biobank (UKB). A GWAS of ML-based VCDR identified 299 independent genome-wide significant (GWS; p ≤ 5 × 10-8) hits in 156 loci. The ML-based GWAS replicated 62 of 65 GWS loci from a recent VCDR GWAS in the UKB for which two ophthalmologists manually labeled images for 67,040 Europeans. The ML-based GWAS also identified 93 novel loci, significantly expanding our understanding of the genetic etiologies of glaucoma and VCDR. Pathway analyses support the biological significance of the novel hits to VCDR: select loci near genes involved in neuronal and synaptic biology or harboring variants are known to cause severe Mendelian ophthalmic disease. Finally, the ML-based GWAS results significantly improve polygenic prediction of VCDR and primary open-angle glaucoma in the independent EPIC-Norfolk cohort.

Associations of sleep apnoea with glaucoma and age-related macular degeneration: an analysis in the United Kingdom Biobank and the Canadian Longitudinal Study on Aging

BACKGROUND

Sleep apnoea, a common sleep-disordered breathing condition, is characterised by upper airway collapse during sleep resulting in transient hypoxia, hypoperfusion of the optic nerve, and spike in intracranial pressure. Previous studies have reported conflicting findings on the association of sleep apnoea with glaucoma, and there are limited reports on the link between sleep apnoea and age-related macular degeneration (AMD).

METHODS

Middle-aged and older participants from the longitudinal United Kingdom (UK) Biobank (n = 502,505) and the Canadian Longitudinal Study on Aging (CLSA; n = 24,073) were included in this analysis. Participants in the UK Biobank and the CLSA were followed for 8 and 3 years, respectively. Participants with diagnosed glaucoma or AMD at baseline were excluded from the analysis. In the UK Biobank, sleep apnoea and incident cases of glaucoma and AMD were identified through hospital inpatient admission, primary care records, and self-reported data. Multivariable Cox proportional hazards models were used to explore associations of sleep apnoea with incidence of glaucoma or AMD.

RESULTS

During the 8-year follow-up in the UK Biobank, glaucoma incidence rates per 1000 person-years were 2.46 and 1.59 for participants with and without sleep apnoea, and the AMD incidence rates per 1000 person-years were 2.27 and 1.42 for participants with and without sleep apnoea, respectively. Multivariable adjusted hazard ratios of glaucoma and AMD risk for sleep apnoea were 1.33 (95% confidence interval [CI] 1.10-1.60, P = 0.003) and 1.39 (95% CI 1.15-1.68, P <  0.001) relative to participants without sleep apnoea. In the CLSA cohort, disease information was collected through in-person interview questionnaires. During the 3-year follow-up, glaucoma incidence rates per 1000 person-years for those with and without sleep apnoea were 9.31 and 6.97, and the AMD incidence rates per 1000 person-years were 8.44 and 6.67, respectively. In the CLSA, similar associations were identified, with glaucoma and AMD odds ratios of 1.43 (95% CI 1.13-1.79) and 1.39 (95% CI 1.08-1.77), respectively, in participants with sleep apnoea compared to those without sleep apnoea (both P <  0.001).

CONCLUSIONS

In two large-scale prospective cohort studies, sleep apnoea is associated with a higher risk of both glaucoma and AMD. These findings indicate that patients with sleep apnoea might benefit from regular ophthalmologic examinations.

Genome-wide analysis of 944 133 individuals provides insights into the etiology of haemorrhoidal disease

OBJECTIVE

Haemorrhoidal disease (HEM) affects a large and silently suffering fraction of the population but its aetiology, including suspected genetic predisposition, is poorly understood. We report the first genome-wide association study (GWAS) meta-analysis to identify genetic risk factors for HEM to date.

DESIGN

We conducted a GWAS meta-analysis of 218 920 patients with HEM and 725 213 controls of European ancestry. Using GWAS summary statistics, we performed multiple genetic correlation analyses between HEM and other traits as well as calculated HEM polygenic risk scores (PRS) and evaluated their translational potential in independent datasets. Using functional annotation of GWAS results, we identified HEM candidate genes, which differential expression and coexpression in HEM tissues were evaluated employing RNA-seq analyses. The localisation of expressed proteins at selected loci was investigated by immunohistochemistry.

RESULTS

We demonstrate modest heritability and genetic correlation of HEM with several other diseases from the GI, neuroaffective and cardiovascular domains. HEM PRS validated in 180 435 individuals from independent datasets allowed the identification of those at risk and correlated with younger age of onset and recurrent surgery. We identified 102 independent HEM risk loci harbouring genes whose expression is enriched in blood vessels and GI tissues, and in pathways associated with smooth muscles, epithelial and endothelial development and morphogenesis. Network transcriptomic analyses highlighted HEM gene coexpression modules that are relevant to the development and integrity of the musculoskeletal and epidermal systems, and the organisation of the extracellular matrix.

CONCLUSION

HEM has a genetic component that predisposes to smooth muscle, epithelial and connective tissue dysfunction.

Polygenic risk modeling with latent trait-related genetic components

Polygenic risk models have led to significant advances in understanding complex diseases and their clinical presentation. While polygenic risk scores (PRS) can effectively predict outcomes, they do not generally account for disease subtypes or pathways which underlie within-trait diversity. Here, we introduce a latent factor model of genetic risk based on components from Decomposition of Genetic Associations (DeGAs), which we call the DeGAs polygenic risk score (dPRS). We compute DeGAs using genetic associations for 977 traits and find that dPRS performs comparably to standard PRS while offering greater interpretability. We show how to decompose an individual's genetic risk for a trait across DeGAs components, with examples for body mass index (BMI) and myocardial infarction (heart attack) in 337,151 white British individuals in the UK Biobank, with replication in a further set of 25,486 non-British white individuals. We find that BMI polygenic risk factorizes into components related to fat-free mass, fat mass, and overall health indicators like physical activity. Most individuals with high dPRS for BMI have strong contributions from both a fat-mass component and a fat-free mass component, whereas a few "outlier" individuals have strong contributions from only one of the two components. Overall, our method enables fine-scale interpretation of the drivers of genetic risk for complex traits.

Genetics of 35 blood and urine biomarkers in the UK Biobank

Clinical laboratory tests are a critical component of the continuum of care. We evaluate the genetic basis of 35 blood and urine laboratory measurements in the UK Biobank (n = 363,228 individuals). We identify 1,857 loci associated with at least one trait, containing 3,374 fine-mapped associations and additional sets of large-effect (>0.1 s.d.) protein-altering, human leukocyte antigen (HLA) and copy number variant (CNV) associations. Through Mendelian randomization (MR) analysis, we discover 51 causal relationships, including previously known agonistic effects of urate on gout and cystatin C on stroke. Finally, we develop polygenic risk scores (PRSs) for each biomarker and build 'multi-PRS' models for diseases using 35 PRSs simultaneously, which improved chronic kidney disease, type 2 diabetes, gout and alcoholic cirrhosis genetic risk stratification in an independent dataset (FinnGen; n = 135,500) relative to single-disease PRSs. Together, our results delineate the genetic basis of biomarkers and their causal influences on diseases and improve genetic risk stratification for common diseases.

Unique roles of rare variants in the genetics of complex diseases in humans

Genome-wide association studies have identified >10,000 genetic variants associated with various phenotypes and diseases. Although the majority are common variants, rare variants with >0.1% of minor allele frequency have been investigated by imputation and using disease-specific custom SNP arrays. Rare variants sequencing analysis mainly revealed have played unique roles in the genetics of complex diseases in humans due to their distinctive features, in contrast to common variants. Unique roles are hypothesis-free evidence for gene causality, a precise target of functional analysis for understanding disease mechanisms, a new favorable target for drug development, and a genetic marker with high disease risk for personalized medicine. As whole-genome sequencing continues to identify more rare variants, the roles associated with rare variants will also increase. However, a better estimation of the functional impact of rare variants across whole genome is needed to enhance their contribution to improvements in human health.

Rare protein-altering variants in ANGPTL7 lower intraocular pressure and protect against glaucoma

Protein-altering variants that are protective against human disease provide in vivo validation of therapeutic targets. Here we use genotyping data from UK Biobank (n = 337,151 unrelated White British individuals) and FinnGen (n = 176,899) to conduct a search for protein-altering variants conferring lower intraocular pressure (IOP) and protection against glaucoma. Through rare protein-altering variant association analysis, we find a missense variant in ANGPTL7 in UK Biobank (rs28991009, p.Gln175His, MAF = 0.8%, genotyped in 82,253 individuals with measured IOP and an independent set of 4,238 glaucoma patients and 250,660 controls) that significantly lowers IOP (β = -0.53 and -0.67 mmHg for heterozygotes, -3.40 and -2.37 mmHg for homozygotes, P = 5.96 x 10-9 and 1.07 x 10-13 for corneal compensated and Goldman-correlated IOP, respectively) and is associated with 34% reduced risk of glaucoma (P = 0.0062). In FinnGen, we identify an ANGPTL7 missense variant at a greater than 50-fold increased frequency in Finland compared with other populations (rs147660927, p.Arg220Cys, MAF Finland = 4.3%), which was genotyped in 6,537 glaucoma patients and 170,362 controls and is associated with a 29% lower glaucoma risk (P = 1.9 x 10-12 for all glaucoma types and also protection against its subtypes including exfoliation, primary open-angle, and primary angle-closure). We further find three rarer variants in UK Biobank, including a protein-truncating variant, which confer a strong composite lowering of IOP (P = 0.0012 and 0.24 for Goldman-correlated and corneal compensated IOP, respectively), suggesting the protective mechanism likely resides in the loss of interaction or function. Our results support inhibition or down-regulation of ANGPTL7 as a therapeutic strategy for glaucoma.