Genetic mechanisms of critical illness in COVID-19.

Host-mediated lung inflammation is present1, and drives mortality2, in the critical illness caused by coronavirus disease 2019 (COVID-19). Host genetic variants associated with critical illness may identify mechanistic targets for therapeutic development3. Here we report the results of the GenOMICC (Genetics Of Mortality In Critical Care) genome-wide association study in 2,244 critically ill patients with COVID-19 from 208 UK intensive care units. We have identified and replicated the following new genome-wide significant associations: on chromosome 12q24.13 (rs10735079, P = 1.65 × 10-8) in a gene cluster that encodes antiviral restriction enzyme activators (OAS1, OAS2 and OAS3); on chromosome 19p13.2 (rs74956615, P = 2.3 × 10-8) near the gene that encodes tyrosine kinase 2 (TYK2); on chromosome 19p13.3 (rs2109069, P = 3.98 ×  10-12) within the gene that encodes dipeptidyl peptidase 9 (DPP9); and on chromosome 21q22.1 (rs2236757, P = 4.99 × 10-8) in the interferon receptor gene IFNAR2. We identified potential targets for repurposing of licensed medications: using Mendelian randomization, we found evidence that low expression of IFNAR2, or high expression of TYK2, are associated with life-threatening disease; and transcriptome-wide association in lung tissue revealed that high expression of the monocyte-macrophage chemotactic receptor CCR2 is associated with severe COVID-19. Our results identify robust genetic signals relating to key host antiviral defence mechanisms and mediators of inflammatory organ damage in COVID-19. Both mechanisms may be amenable to targeted treatment with existing drugs. However, large-scale randomized clinical trials will be essential before any change to clinical practice

Mapping the human genetic architecture of COVID-19

The genetic make-up of an individual contributes to the susceptibility and response to viral infection. Although environmental, clinical and social factors have a role in the chance of exposure to SARS-CoV-2 and the severity of COVID-191,2, host genetics may also be important. Identifying host-specific genetic factors may reveal biological mechanisms of therapeutic relevance and clarify causal relationships of modifiable environmental risk factors for SARS-CoV-2 infection and outcomes. We formed a global network of researchers to investigate the role of human genetics in SARS-CoV-2 infection and COVID-19 severity. Here we describe the results of three genome-wide association meta-analyses that consist of up to 49,562 patients with COVID-19 from 46 studies across 19 countries. We report 13 genome-wide significant loci that are associated with SARS-CoV-2 infection or severe manifestations of COVID-19. Several of these loci correspond to previously documented associations to lung or autoimmune and inflammatory diseases3–7. They also represent potentially actionable mechanisms in response to infection. Mendelian randomization analyses support a causal role for smoking and body-mass index for severe COVID-19 although not for type II diabetes. The identification of novel host genetic factors associated with COVID-19 was made possible by the community of human genetics researchers coming together to prioritize the sharing of data, results, resources and analytical frameworks. This working model of international collaboration underscores what is possible for future genetic discoveries in emerging pandemics, or indeed for any complex human disease

Pretreatment prediction of response to ursodeoxycholic acid in primary biliary cholangitis: development and validation of the UDCA Response Score.

BACKGROUND Treatment guidelines recommend a stepwise approach to primary biliary cholangitis: all patients begin treatment with ursodeoxycholic acid (UDCA) monotherapy and those with an inadequate biochemical response after 12 months are subsequently considered for second-line therapies. However, as a result, patients at the highest risk can wait the longest for effective treatment. We determined whether UDCA response can be accurately predicted using pretreatment clinical parameters. METHODS We did logistic regression analysis of pretreatment variables in a discovery cohort of patients in the UK with primary biliary cholangitis to derive the best-fitting model of UDCA response, defined as alkaline phosphatase less than 1·67 times the upper limit of normal (ULN), measured after 12 months of treatment with UDCA. We validated the model in an external cohort of patients with primary biliary cholangitis and treated with UDCA in Italy. Additionally, we assessed correlations between model predictions and key histological features, such as biliary injury and fibrosis, on liver biopsy samples. FINDINGS 2703 participants diagnosed with primary biliary cholangitis between Jan 1, 1998, and May 31, 2015, were included in the UK-PBC cohort for derivation of the model. The following pretreatment parameters were associated with lower probability of UDCA response: higher alkaline phosphatase concentration (p<0·0001), higher total bilirubin concentration (p=0·0003), lower aminotransferase concentration (p=0·0012), younger age (p<0·0001), longer interval from diagnosis to the start of UDCA treatment (treatment time lag, p<0·0001), and worsening of alkaline phosphatase concentration from diagnosis (p<0·0001). Based on these variables, we derived a predictive score of UDCA response. In the external validation cohort, 460 patients diagnosed with primary biliary cholangitis were treated with UDCA, with follow-up data until May 31, 2016. In this validation cohort, the area under the receiver operating characteristic curve for the score was 0·83 (95% CI 0·79-0·87). In 20 liver biopsy samples from patients with primary biliary cholangitis, the UDCA response score was associated with ductular reaction (r=-0·556, p=0·0130) and intermediate hepatocytes (probability of response was 0·90 if intermediate hepatocytes were absent vs 0·51 if present). INTERPRETATION We have derived and externally validated a model based on pretreatment variables that accurately predicts UDCA response. Association with histological features provides face validity. This model provides a basis to explore alternative approaches to treatment stratification in patients with primary biliary cholangitis. FUNDING UK Medical Research Council and University of Milan-Bicocca

Pretreatment prediction of response to ursodeoxycholic acid in primary biliary cholangitis: development and validation of the UDCA Response Score.

BACKGROUND: Treatment guidelines recommend a stepwise approach to primary biliary cholangitis: all patients begin treatment with ursodeoxycholic acid (UDCA) monotherapy and those with an inadequate biochemical response after 12 months are subsequently considered for second-line therapies. However, as a result, patients at the highest risk can wait the longest for effective treatment. We determined whether UDCA response can be accurately predicted using pretreatment clinical parameters. METHODS: We did logistic regression analysis of pretreatment variables in a discovery cohort of patients in the UK with primary biliary cholangitis to derive the best-fitting model of UDCA response, defined as alkaline phosphatase less than 1·67 times the upper limit of normal (ULN), measured after 12 months of treatment with UDCA. We validated the model in an external cohort of patients with primary biliary cholangitis and treated with UDCA in Italy. Additionally, we assessed correlations between model predictions and key histological features, such as biliary injury and fibrosis, on liver biopsy samples. FINDINGS: 2703 participants diagnosed with primary biliary cholangitis between Jan 1, 1998, and May 31, 2015, were included in the UK-PBC cohort for derivation of the model. The following pretreatment parameters were associated with lower probability of UDCA response: higher alkaline phosphatase concentration (p<0·0001), higher total bilirubin concentration (p=0·0003), lower aminotransferase concentration (p=0·0012), younger age (p<0·0001), longer interval from diagnosis to the start of UDCA treatment (treatment time lag, p<0·0001), and worsening of alkaline phosphatase concentration from diagnosis (p<0·0001). Based on these variables, we derived a predictive score of UDCA response. In the external validation cohort, 460 patients diagnosed with primary biliary cholangitis were treated with UDCA, with follow-up data until May 31, 2016. In this validation cohort, the area under the receiver operating characteristic curve for the score was 0·83 (95% CI 0·79-0·87). In 20 liver biopsy samples from patients with primary biliary cholangitis, the UDCA response score was associated with ductular reaction (r=-0·556, p=0·0130) and intermediate hepatocytes (probability of response was 0·90 if intermediate hepatocytes were absent vs 0·51 if present). INTERPRETATION: We have derived and externally validated a model based on pretreatment variables that accurately predicts UDCA response. Association with histological features provides face validity. This model provides a basis to explore alternative approaches to treatment stratification in patients with primary biliary cholangitis. FUNDING: UK Medical Research Council and University of Milan-Bicocca

X Chromosome Contribution to the Genetic Architecture of Primary Biliary Cholangitis.

BACKGROUND & AIMS: Genome-wide association studies in primary biliary cholangitis (PBC) have failed to find X chromosome (chrX) variants associated with the disease. Here, we specifically explore the chrX contribution to PBC, a sexually dimorphic complex autoimmune disease. METHODS: We performed a chrX-wide association study, including genotype data from 5 genome-wide association studies (from Italy, United Kingdom, Canada, China, and Japan; 5244 case patients and 11,875 control individuals). RESULTS: Single-marker association analyses found approximately 100 loci displaying P < 5 × 10(-4), with the most significant being a signal within the OTUD5 gene (rs3027490; P = 4.80 × 10(-6); odds ratio [OR], 1.39; 95% confidence interval [CI], 1.028-1.88; Japanese cohort). Although the transethnic meta-analysis evidenced only a suggestive signal (rs2239452, mapping within the PIM2 gene; OR, 1.17; 95% CI, 1.09-1.26; P = 9.93 × 10(-8)), the population-specific meta-analysis showed a genome-wide significant locus in East Asian individuals pointing to the same region (rs7059064, mapping within the GRIPAP1 gene; P = 6.2 × 10(-9); OR, 1.33; 95% CI, 1.21-1.46). Indeed, rs7059064 tags a unique linkage disequilibrium block including 7 genes: TIMM17B, PQBP1, PIM2, SLC35A2, OTUD5, KCND1, and GRIPAP1, as well as a superenhancer (GH0XJ048933 within OTUD5) targeting all these genes. GH0XJ048933 is also predicted to target FOXP3, the main T-regulatory cell lineage specification factor. Consistently, OTUD5 and FOXP3 RNA levels were up-regulated in PBC case patients (1.75- and 1.64-fold, respectively). CONCLUSIONS: This work represents the first comprehensive study, to our knowledge, of the chrX contribution to the genetics of an autoimmune liver disease and shows a novel PBC-related genome-wide significant locus.The article is available via Open Access. Click on the 'Additional link' above to access the full-text.Published version, accepted versio

GWAS and meta-analysis identifies 49 genetic variants underlying critical COVID-19

Data availability: Downloadable summary data are available through the GenOMICC data site (https://genomicc.org/data). Summary statistics are available, but without the 23andMe summary statistics, except for the 10,000 most significant hits, for which full summary statistics are available. The full GWAS summary statistics for the 23andMe discovery dataset will be made available through 23andMe to qualified researchers under an agreement with 23andMe that protects the privacy of the 23andMe participants. For further information and to apply for access to the data, see the 23andMe website (https://research.23andMe.com/dataset-access/). All individual-level genotype and whole-genome sequencing data (for both academic and commercial uses) can be accessed through the UKRI/HDR UK Outbreak Data Analysis Platform (https://odap.ac.uk). A restricted dataset for a subset of GenOMICC participants is also available through the Genomics England data service. Monocyte RNA-seq data are available under the title ‘Monocyte gene expression data’ within the Oxford University Research Archives (https://doi.org/10.5287/ora-ko7q2nq66). Sequencing data will be made freely available to organizations and researchers to conduct research in accordance with the UK Policy Framework for Health and Social Care Research through a data access agreement. Sequencing data have been deposited at the European Genome–Phenome Archive (EGA), which is hosted by the EBI and the CRG, under accession number EGAS00001007111.Extended data figures and tables are available online at https://www.nature.com/articles/s41586-023-06034-3#Sec21 .Supplementary information is available online at https://www.nature.com/articles/s41586-023-06034-3#Sec22 .Code availability: Code to calculate the imputation of P values on the basis of SNPs in linkage disequilibrium is available at GitHub (https://github.com/baillielab/GenOMICC_GWAS).Acknowledgements: We thank the members of the Banco Nacional de ADN and the GRA@CE cohort group; and the research participants and employees of 23andMe for making this work possible. A full list of contributors who have provided data that were collated in the HGI project, including previous iterations, is available online (https://www.covid19hg.org/acknowledgements).Change history: 11 July 2023: A Correction to this paper has been published at: https://doi.org/10.1038/s41586-023-06383-z. -- In the version of this article initially published, the name of Ana Margarita Baldión-Elorza, of the SCOURGE Consortium, appeared incorrectly (as Ana María Baldion) and has now been amended in the HTML and PDF versions of the article.Copyright © The Author(s) 2023, Critical illness in COVID-19 is an extreme and clinically homogeneous disease phenotype that we have previously shown1 to be highly efficient for discovery of genetic associations2. Despite the advanced stage of illness at presentation, we have shown that host genetics in patients who are critically ill with COVID-19 can identify immunomodulatory therapies with strong beneficial effects in this group3. Here we analyse 24,202 cases of COVID-19 with critical illness comprising a combination of microarray genotype and whole-genome sequencing data from cases of critical illness in the international GenOMICC (11,440 cases) study, combined with other studies recruiting hospitalized patients with a strong focus on severe and critical disease: ISARIC4C (676 cases) and the SCOURGE consortium (5,934 cases). To put these results in the context of existing work, we conduct a meta-analysis of the new GenOMICC genome-wide association study (GWAS) results with previously published data. We find 49 genome-wide significant associations, of which 16 have not been reported previously. To investigate the therapeutic implications of these findings, we infer the structural consequences of protein-coding variants, and combine our GWAS results with gene expression data using a monocyte transcriptome-wide association study (TWAS) model, as well as gene and protein expression using Mendelian randomization. We identify potentially druggable targets in multiple systems, including inflammatory signalling (JAK1), monocyte–macrophage activation and endothelial permeability (PDE4A), immunometabolism (SLC2A5 and AK5), and host factors required for viral entry and replication (TMPRSS2 and RAB2A).GenOMICC was funded by Sepsis Research (the Fiona Elizabeth Agnew Trust), the Intensive Care Society, a Wellcome Trust Senior Research Fellowship (to J.K.B., 223164/Z/21/Z), the Department of Health and Social Care (DHSC), Illumina, LifeArc, the Medical Research Council, UKRI, a BBSRC Institute Program Support Grant to the Roslin Institute (BBS/E/D/20002172, BBS/E/D/10002070 and BBS/E/D/30002275) and UKRI grants MC_PC_20004, MC_PC_19025, MC_PC_1905 and MRNO2995X/1. A.D.B. acknowledges funding from the Wellcome PhD training fellowship for clinicians (204979/Z/16/Z), the Edinburgh Clinical Academic Track (ECAT) programme. This research is supported in part by the Data and Connectivity National Core Study, led by Health Data Research UK in partnership with the Office for National Statistics and funded by UK Research and Innovation (grant MC_PC_20029). Laboratory work was funded by a Wellcome Intermediate Clinical Fellowship to B.F. (201488/Z/16/Z). We acknowledge the staff at NHS Digital, Public Health England and the Intensive Care National Audit and Research Centre who provided clinical data on the participants; and the National Institute for Healthcare Research Clinical Research Network (NIHR CRN) and the Chief Scientist’s Office (Scotland), who facilitate recruitment into research studies in NHS hospitals, and to the global ISARIC and InFACT consortia. GenOMICC genotype controls were obtained using UK Biobank Resource under project 788 funded by Roslin Institute Strategic Programme Grants from the BBSRC (BBS/E/D/10002070 and BBS/E/D/30002275) and Health Data Research UK (HDR-9004 and HDR-9003). UK Biobank data were used in the GSMR analyses presented here under project 66982. The UK Biobank was established by the Wellcome Trust medical charity, Medical Research Council, Department of Health, Scottish Government and the Northwest Regional Development Agency. It has also had funding from the Welsh Assembly Government, British Heart Foundation and Diabetes UK. The work of L.K. was supported by an RCUK Innovation Fellowship from the National Productivity Investment Fund (MR/R026408/1). J.Y. is supported by the Westlake Education Foundation. SCOURGE is funded by the Instituto de Salud Carlos III (COV20_00622 to A.C., PI20/00876 to C.F.), European Union (ERDF) ‘A way of making Europe’, Fundación Amancio Ortega, Banco de Santander (to A.C.), Cabildo Insular de Tenerife (CGIEU0000219140 ‘Apuestas científicas del ITER para colaborar en la lucha contra la COVID-19’ to C.F.) and Fundación Canaria Instituto de Investigación Sanitaria de Canarias (PIFIISC20/57 to C.F.). We also acknowledge the contribution of the Centro National de Genotipado (CEGEN) and Centro de Supercomputación de Galicia (CESGA) for funding this project by providing supercomputing infrastructures. A.D.L. is a recipient of fellowships from the National Council for Scientific and Technological Development (CNPq)-Brazil (309173/2019-1 and 201527/2020-0)

Mapping the human genetic architecture of COVID-19

The genetic make-up of an individual contributes to the susceptibility and response to viral infection. Although environmental, clinical and social factors have a role in the chance of exposure to SARS-CoV-2 and the severity of COVID-191,2, host genetics may also be important. Identifying host-specific genetic factors may reveal biological mechanisms of therapeutic relevance and clarify causal relationships of modifiable environmental risk factors for SARS-CoV-2 infection and outcomes. We formed a global network of researchers to investigate the role of human genetics in SARS-CoV-2 infection and COVID-19 severity. Here we describe the results of three genome-wide association meta-analyses that consist of up to 49,562 patients with COVID-19 from 46 studies across 19 countries. We report 13 genome-wide significant loci that are associated with SARS-CoV-2 infection or severe manifestations of COVID-19. Several of these loci correspond to previously documented associations to lung or autoimmune and inflammatory diseases3,4,5,6,7. They also represent potentially actionable mechanisms in response to infection. Mendelian randomization analyses support a causal role for smoking and body-mass index for severe COVID-19 although not for type II diabetes. The identification of novel host genetic factors associated with COVID-19 was made possible by the community of human genetics researchers coming together to prioritize the sharing of data, results, resources and analytical frameworks. This working model of international collaboration underscores what is possible for future genetic discoveries in emerging pandemics, or indeed for any complex human disease