ABSTRACT BACKGROUND Disentangling the genetic constellation underlying Alzheimer’s disease (AD) is important. Doing so allows us to identify biological pathways underlying AD, point towards novel drug targets and use the variants for individualised risk predictions in disease modifying or prevention trials. In the present work we report on the largest genome-wide association study (GWAS) for AD risk to date and show the combined utility of proven AD loci for precision medicine using polygenic risk scores (PRS). METHODS Three sets of summary statistics were included in our meta-GWAS of AD: an Spanish case-control study (GR@ACE/DEGESCO study, n = 12,386), the case-control study of International Genomics of Alzheimer project (IGAP, n = 82,771) and the UK Biobank (UKB) AD-by-proxy case-control study (n=314,278). Using these resources, we performed a fixed-effects inverse-variance-weighted meta-analysis. Detected loci were confirmed in a replication study of 19,089 AD cases and 39,101 controls from 16 European-ancestry cohorts not previously used. We constructed a weighted PRS based on the 39 AD variants. PRS were generated by multiplying the genotype dosage of each risk allele for each variant by its respective weight, and then summing across all variants. We first validated it for AD in independent data (assessing effects of sub-threshold signal, diagnostic certainty, age at onset and sex) and tested its effect on risk (odds for disease) and age at onset in the GR@ACE/DEGESCO study. FINDINGS Using our meta-GWAS approach and follow-up analysis, we identified novel genome-wide significant associations of six genetic variants with AD risk (rs72835061 -CHRNE , rs2154481 -APP , rs876461 -PRKD3/NDUFAF7 , rs3935877 -PLCG2 and two missense variants: rs34173062/rs34674752 in SHARPIN gene) and confirmed a stop codon mutation in the IL34 gene increasing the risk of AD ( IL34-Tyr213Ter ), and two other variants in PLCG2 and HS3ST1 regions. This brings the total number of genetic variants associated with AD to 39 (excluding APOE ). The PRS based on these variants was associated with AD in an independent clinical AD-case control dataset (OR=1.30, per 1-SD increase in the PRS, 95%CI 1.18-1.44, p = 1.1×10 −7 ), a similar effect to that in the GR@ACE/DEGESCO (OR=1.27, 95%CI 1.23-1.32, p = 7.4×10 −39 ). We then explored the combined effects of these 39 variants in a PRS for AD risk and age-at-onset stratification in GR@ACE/DEGESCO. Excluding APOE , we observed a gradual risk increase over the 2% tiles; when comparing the extremes, those with the 2% highest risk had a 2.98-fold (95% CI 2.12–4.18, p = 3.2×10 −10 ) increased risk compared to those with the 2% lowest risk ( p = 5.9×10 −10 ). Using the PRS we identified APOE ε33 carriers with a similar risk as APOE ε 4 heterozygotes carriers, as well as APOE ε4 heterozygote carriers with a similar risk as APOE ε 4 homozygote. Considering age at onset; there was a 9-year difference between median onset of AD the lowest risk group and the highest risk group (82 vs 73 years; p = 1.6×10 −6 ); a 4-year median onset difference (81 vs 77 years; p = 6.9×10 −5 ) within APOE ε4 heterozygotes and a 5.5-year median onset difference (78.5 vs 73 years; p = 4.6×10 −5 ) within APOE ε4 carriers. INTERPRETATION We identified six novel genetic variants associated with AD-risk, among which one common APP variant. A PRS of all genetic loci reported to date could be a robust tool to predict the risk and age at onset of AD, beyond APOE alone. These properties make PRS instrumental in selecting individuals at risk in order to apply preventative strategies and might have potential use in diagnostic work-up.

Long runs of homozygosity (ROH) are contiguous stretches of homozygous genotypes, which are a footprint of inbreeding and recessive inheritance. The presence of recessive loci is suggested for Alzheimer's disease (AD); however, their search has been poorly assessed to date. To investigate homozygosity in AD, here we performed a fine-scale ROH analysis using 10 independent cohorts of European ancestry (11,919 AD cases and 9181 controls.) We detected an increase of homozygosity in AD cases compared to controls [βAVROH (CI 95%) = 0.070 (0.037-0.104); P = 3.91 × 10-5; βFROH (CI95%) = 0.043 (0.009-0.076); P = 0.013]. ROHs increasing the risk of AD (OR > 1) were significantly overrepresented compared to ROHs increasing protection (p < 2.20 × 10-16). A significant ROH association with AD risk was detected upstream the HS3ST1 locus (chr4:11,189,482‒11,305,456), (β (CI 95%) = 1.09 (0.48 ‒ 1.48), p value = 9.03 × 10-4), previously related to AD. Next, to search for recessive candidate variants in ROHs, we constructed a homozygosity map of inbred AD cases extracted from an outbred population and explored ROH regions in whole-exome sequencing data (N = 1449). We detected a candidate marker, rs117458494, mapped in the SPON1 locus, which has been previously associated with amyloid metabolism. Here, we provide a research framework to look for recessive variants in AD using outbred populations. Our results showed that AD cases have enriched homozygosity, suggesting that recessive effects may explain a proportion of AD heritability.

Emerging studies have suggested several chromosomal regions as potential host genetic factors involved in the susceptibility to SARS-CoV-2 infection and disease outcome. We nested a COVID-19 genome-wide association study using the GR@ACE/DEGESCO study, searching for susceptibility factors associated with COVID-19 disease. To this end, we compared 221 COVID-19 confirmed cases with 17,035 individuals in whom the COVID-19 disease status was unknown. Then, we performed a meta-analysis with the publicly available data from the COVID-19 Host Genetics Initiative. Because the APOE locus has been suggested as a potential modifier of COVID-19 disease, we added sensitivity analyses stratifying by dementia status or by disease severity. We confirmed the existence of the 3p21.31 region (LZTFL1, SLC6A20) implicated in the susceptibility to SARS-CoV-2 infection and TYK2 gene might be involved in COVID-19 severity. Nevertheless, no statistically significant association was observed in the COVID-19 fatal outcome or in the stratified analyses (dementia-only and non-dementia strata) for the APOE locus not supporting its involvement in SARS-CoV-2 pathobiology or COVID-19 prognosis.

Liver stiffness (LS) at sustained virological response (SVR) after direct-acting antivirals (DAA)-based therapy is a predictor of liver events in hepatitis C virus (HCV)-infected patients. The study aim was to identify genetic factors associated with LS changes from the moment of starting anti-HCV therapy to SVR. This prospective study included HCV-infected patients from the GEHEP-011 cohort who achieved SVR with DAA-based therapy, with LS pre-treatment ≥ 9.5 kPa and LS measurement available at SVR. Plink and Magma software were used to carry out genome-wide single-nucleotide polymorphism (SNP)-based and gene-based association analyses, respectively. The ShinyGO application was used for exploring enrichment in Gene Ontology (GO) categories for biological processes. Overall, 242 patients were included. Median (quartile 1, quartile 3) LS values at pre-treatment and at SVR were 16.8 (12, 28) kPa and 12.0 (8.5, 19.3) kPa, respectively. Thirty-five SNPs and three genes reached suggestive association with LS changes from the moment of starting anti-HCV therapy to SVR. GO categories related to DNA packaging complex, DNA conformation change, chromosome organization and chromatin organization were significantly enriched. Our study reports possible genetic factors associated with LS changes during HCV-infection cure. In addition, our results suggest that processes related to DNA conformation are also involved in these changes.

ABSTRACT Long runs of homozygosity (ROH) are contiguous stretches of homozygous genotypes, which are a footprint of recent inbreeding and recessive inheritance. The presence of recessive loci is suggested for Alzheimer’s disease (AD). However, the search for recessive variants has been poorly assessed to date. To investigate homozygosity in AD, we performed a fine-scale ROH analysis including 21,100 individuals from 10 cohorts of European ancestry (11,919 AD cases and 9,181 controls). We detected an increase of homozygosity in AD cases compared to controls [β FROH (CI95%) = 0.051 (0.023 – 0.078); P = 3.25 x 10 -4 ]. ROHs increasing the risk of AD (OR > 1) were significantly overrepresented compared to ROHs increasing protection (p < 2.20 x 10 -16 ). The top associated ROH with AD risk (β (CI95%) = 1.09 (0.48 ‒ 1.48), p value = 9.03 x 10 -4 ) was detected upstream the HS3ST1 locus (chr4:11,189,482‒11,305,456), previously related to AD. Next, to construct a homozygosity map of AD cases, we selected ROHs shared by inbred AD cases extracted from an outbred population. We used whole-exome sequencing data from 1,449 individuals from the Knight-ADRC-NIA-LOAD (KANL) cohort to identify potential recessive variants in candidate ROHs. We detected a candidate marker, rs117458494, mapped in the SPON1 locus, which has been previously associated with amyloid metabolism. Here, we provide a research framework to look for recessive variants in AD using outbred populations. Our results showed that AD cases have enriched homozygosity, suggesting that recessive effects may explain a proportion of AD heritability.