Search

Abstract The tripartite ancestral structure is a recently proposed model for the genetic origin of modern Japanese, comprising indigenous Jomon hunter-gatherers and two additional continental ancestors from Northeast Asia and East Asia. To investigate the impact of the tripartite structure on genetic and phenotypic variation today, we conducted biobank-scale analyses by merging Biobank Japan (BBJ; n = 171,287) with ancient Japanese and Eurasian genomes (n = 22). We demonstrate the applicability of the tripartite model to Japanese populations throughout the archipelago, with an extremely strong correlation between Jomon ancestry and genomic variation among individuals. We also find that the genetic legacy of Jomon ancestry underlies an elevated body mass index (BMI). Genome-wide association analysis with rigorous adjustments for geographical and ancestral substructures identifies 132 variants that are informative for predicting individual Jomon ancestry. This prediction model is validated using independent Japanese cohorts (Nagahama cohort, n = 2993; the second cohort of BBJ, n = 72,695). We further confirm the phenotypic association between Jomon ancestry and BMI using East Asian individuals from UK Biobank (n = 2286). Our extensive analysis of ancient and modern genomes, involving over 250,000 participants, provides valuable insights into the genetic legacy of ancient hunter-gatherers in contemporary populations.

Abstract Background Although polygenic risk scores (PRSs) are expected to be helpful in precision medicine, it remains unclear whether high-PRS groups are more likely to benefit from preventive interventions for diseases. Recent methodological advancements enable us to predict treatment effects at the individual level. Methods We employed causal forest to explore the relationship between PRSs and individual risk of diseases associated with certain environmental factors. Following simulations illustrating its performance, we applied our approach to investigate the individual risk of cardiometabolic diseases, including coronary artery diseases (CAD) and type 2 diabetes (T2D), associated with obesity and smoking among individuals from UK Biobank (UKB; n = 369,942) and BioBank Japan (BBJ; n = 149,421). Results Here we find the heterogeneous association of obesity and smoking with diseases across PRS values, complicated by the multi-dimensional combination of individual characteristics such as age and sex. The highest positive correlations of PRSs and the exposure-related disease risks are observed between obesity and T2D in UKB and between smoking and CAD in BBJ (Spearman’s ρ = 0.61 and 0.32, respectively). However, most relationships are weak or negative, suggesting that high-PRS groups will not necessarily benefit most from environmental factor prevention. Conclusions Our study highlights the importance of individual-level prediction of disease risks associated with target exposure in precision medicine.

Disorders associated with the immune system burden multiple organs, although the shared biology exists across the diseases. Preceding family-based studies reveal that immune diseases are heritable to varying degrees, providing the basis for immunogenomics. The recent cost reduction in genetic analysis intensively promotes biobank-scale studies and the development of frameworks for statistical genetics. The accumulating multi-layer omics data, including genome-wide association studies (GWAS) and RNA-sequencing at single-cell resolution, enable us to dissect the genetic backgrounds of immune-related disorders. Although autoimmune and allergic diseases are generally categorized into different disease categories, epidemiological studies reveal the high incidence of autoimmune and allergic disease complications, suggesting the shared genetics and biology between the disease categories. Biobank resources and consortia cover multiple immune-related disorders to accumulate phenome-wide associations of genetic variants and enhance researchers to analyze the shared and heterogeneous genetic backgrounds. The emerging post-GWAS and integrative multi-omics analyses provide genetic and biological insights into the multicategorical disease associations.

Abstract Elevated blood pressure (BP) is major risk factor for cardiovascular diseases (CVD). Genome-wide association studies (GWAS) conducted predominantly in populations of European ancestry have identified >2,000 BP-associated loci, but other ancestries have been less well-studied. We conducted GWAS of systolic, diastolic, pulse, and mean arterial BP in 100,453 Chinese adults. We identified 128 non-overlapping loci associated with one or more BP traits, including 74 newly-reported associations. Despite strong genetic correlations between populations, we identified appreciably higher heritability and larger variant effect sizes in Chinese compared with European or Japanese ancestry populations. Using instruments derived from these GWAS, multivariable Mendelian randomisation demonstrated that BP traits contribute differently to the causal associations of BP with CVD. In particular, only pulse pressure was independently causally associated with carotid plaque. These findings reinforce the need for studies in diverse populations to understand the genetic determinants of BP traits and their roles in disease risk.

Abstract Recurrent pregnancy loss (RPL) is a major reproductive health issue with multifactorial causes, affecting 2.6% of all pregnancies worldwide. Nearly half of the RPL cases lack clinically identifiable causes (e.g., antiphospholipid syndrome, uterine anomalies, and parental chromosomal abnormalities), referred to as unexplained RPL (uRPL). Here, we perform a genome-wide association study focusing on uRPL in 1,728 cases and 24,315 female controls of Japanese ancestry. We detect significant associations in the major histocompatibility complex (MHC) region at 6p21 (lead variant=rs9263738; P = 1.4 × 10−10; odds ratio [OR] = 1.51 [95% CI: 1.33–1.72]; risk allele frequency = 0.871). The MHC associations are fine-mapped to the classical HLA alleles, HLA-C*12:02, HLA-B*52:01, and HLA-DRB1*15:02 (P = 1.1 × 10−10, 1.5 × 10−10, and 1.2 × 10−9, respectively), which constitute a population-specific common long-range haplotype with a protective effect (P = 2.8 × 10−10; OR = 0.65 [95% CI: 0.57–0.75]; haplotype frequency=0.108). Genome-wide copy-number variation (CNV) calling demonstrates rare predicted loss-of-function (pLoF) variants of the cadherin-11 gene (CDH11) conferring the risk of uRPL (P = 1.3 × 10−4; OR = 3.29 [95% CI: 1.78–5.76]). Our study highlights the importance of reproductive immunology and rare variants in the uRPL etiology.

Abstract Polygenic scores (PGSs) offer the ability to predict genetic risk for complex diseases across the life course; a key benefit over short-term prediction models. To produce risk estimates relevant to clinical and public health decision-making, it is important to account for varying effects due to age and sex. Here, we develop a novel framework to estimate country-, age-, and sex-specific estimates of cumulative incidence stratified by PGS for 18 high-burden diseases. We integrate PGS associations from seven studies in four countries (N = 1,197,129) with disease incidences from the Global Burden of Disease. PGS has a significant sex-specific effect for asthma, hip osteoarthritis, gout, coronary heart disease and type 2 diabetes (T2D), with all but T2D exhibiting a larger effect in men. PGS has a larger effect in younger individuals for 13 diseases, with effects decreasing linearly with age. We show for breast cancer that, relative to individuals in the bottom 20% of polygenic risk, the top 5% attain an absolute risk for screening eligibility 16.3 years earlier. Our framework increases the generalizability of results from biobank studies and the accuracy of absolute risk estimates by appropriately accounting for age- and sex-specific PGS effects. Our results highlight the potential of PGS as a screening tool which may assist in the early prevention of common diseases.

Abstract Background Renal impairment is a predictor of coronavirus disease (COVID-19) severity. No studies have compared COVID-19 outcomes in patients with chronic kidney disease (CKD) and patients with impaired renal function without a prior diagnosis of CKD. This study aimed to identify the impact of pre-existing impaired renal function without CKD on COVID-19 outcomes. Methods This retrospective study included 3,637 patients with COVID-19 classified into three groups by CKD history and estimated glomerular filtration rate (eGFR) on referral: Group 1 (n = 2,460), normal renal function without a CKD history; Group 2 (n = 905), impaired renal function without a CKD history; and Group 3 (n = 272), history of CKD. We compared the clinical characteristics of these groups and assessed the effect of CKD and impaired renal function on critical outcomes (requirement for respiratory support with high-flow oxygen devices, invasive mechanical ventilation, or extracorporeal membrane oxygen, and death during hospitalization) using multivariable logistic regression. Results The prevalence of comorbidities (hypertension, diabetes, and cardiovascular disease) and incidence of inflammatory responses (white blood counts, and C-reactive protein, procalcitonin, and D-dimer levels) and complications (bacterial infection and heart failure) were higher in Groups 2 and 3 than that in Group 1. The incidence of critical outcomes was 10.8%, 17.7%, and 26.8% in Groups 1, 2, and 3, respectively. The mortality rate and the rate of requiring IMV support was lowest in Group 1 and highest in Group 3. Compared with Group 1, the risk of critical outcomes was higher in Group 2 (adjusted odds ratio [aOR]: 1.32, 95% confidence interval [CI]: 1.03–1.70, P = 0.030) and Group 3 (aOR: 1.94, 95% CI: 1.36–2.78, P

Background Observational studies have reported associations between primary aldosteronism (PA) and cardiovascular outcomes, including coronary artery diseases (CAD), congestive heart failure (CHF), and stroke. However, establishing causality remains a challenge due to the lack of randomized controlled trial data on this topic. We thus aimed to investigate the causal relationship between PA and the risk of developing CAD, CHF, and stroke. Methods and Results Cross‐ancestry meta‐analysis of genome‐wide association studies combining East Asian and European ancestry (1560 PA cases and 742 139 controls) was conducted to identify single‐nucleotide variants that are associated with PA. Then, using the identified genetic variants as instrumental variables, we conducted the 2‐sample Mendelian randomization analysis to investigate the causal relationship between PA and incident CAD, CHF, and stroke among both East Asian and European ancestry. Summary association results were extracted from large genome‐wide association studies consortia. Our cross‐ancestry meta‐analysis of East Asian and European populations identified 7 genetic loci significantly associated with the risk of PA, for which the genes nearest to the lead variants were CASZ1, WNT2B, HOTTIP, LSP1, TBX3, RXFP2, and NDP. Among the East Asian population, the pooled odds ratio estimates using these 7 genetic instruments of PA were 1.07 (95% CI, 1.03–1.11) for CAD, 1.10 (95% CI, 1.01–1.20) for CHF, and 1.13 (95% CI, 1.09–1.18) for stroke. The results were consistent among the European population. Conclusions Our 2‐sample Mendelian randomization study revealed that PA had increased risks of CAD, CHF, and stroke. These findings highlight that early and active screening of PA is critical to prevent future cardiovascular events.

Objective This study aimed to investigate the utility of CT quantification of lung volume for predicting critical outcomes in COVID-19 patients.Methods This retrospective cohort study included 1200 hospitalised patients with COVID-19 from 4 hospitals. Lung fields were extracted using artificial intelligence-based segmentation, and the percentage of the predicted (%pred) total lung volume (TLC (%pred)) was calculated. The incidence of critical outcomes and posthospitalisation complications was compared between patients with low and high CT lung volumes classified based on the median percentage of predicted TLCct (n=600 for each). Prognostic factors for residual lung volume loss were investigated in 208 patients with COVID-19 via a follow-up CT after 3 months.Results The incidence of critical outcomes was higher in the low TLCct (%pred) group than in the high TLCct (%pred) group (14.2% vs 3.3%, p

Abstract To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Abstract Here we report the largest Asian genome-wide association study (GWAS) for systemic sclerosis performed to date, based on data from Japanese subjects and comprising of 1428 cases and 112,599 controls. The lead SNP is in the FCGR/FCRL region, which shows a penetrating association in the Asian population, while a complete linkage disequilibrium SNP, rs10917688, is found in a cis-regulatory element for IRF8. IRF8 is also a significant locus in European GWAS for systemic sclerosis, but rs10917688 only shows an association in the presence of the risk allele of IRF8 in the Japanese population. Further analysis shows that rs10917688 is marked with H3K4me1 in primary B cells. A meta-analysis with a European GWAS detects 30 additional significant loci. Polygenic risk scores constructed with the effect sizes of the meta-analysis suggest the potential portability of genetic associations beyond populations. Prioritizing the top 5% of SNPs of IRF8 binding sites in B cells improves the fitting of the polygenic risk scores, underscoring the roles of B cells and IRF8 in the development of systemic sclerosis. The results also suggest that systemic sclerosis shares a common genetic architecture across populations.

Abstract X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Abstract The severity of chest X-ray (CXR) findings is a prognostic factor in patients with coronavirus disease 2019 (COVID-19). We investigated the clinical and genetic characteristics and prognosis of patients with worsening CXR findings during early hospitalization. We retrospectively included 1656 consecutive Japanese patients with COVID-19 recruited through the Japan COVID-19 Task Force. Rapid deterioration of CXR findings was defined as increased pulmonary infiltrates in ≥ 50% of the lung fields within 48 h of admission. Rapid deterioration of CXR findings was an independent risk factor for death, most severe illness, tracheal intubation, and intensive care unit admission. The presence of consolidation on CXR, comorbid cardiovascular and chronic obstructive pulmonary diseases, high body temperature, and increased serum aspartate aminotransferase, potassium, and C-reactive protein levels were independent risk factors for rapid deterioration of CXR findings. Risk variant at the ABO locus (rs529565-C) was associated with rapid deterioration of CXR findings in all patients. This study revealed the clinical features, genetic features, and risk factors associated with rapid deterioration of CXR findings, a poor prognostic factor in patients with COVID-19.

Abstract The immunological basis of the clinical heterogeneity in autoimmune vasculitis remains poorly understood. In this study, we conduct single-cell transcriptome analyses on peripheral blood mononuclear cells (PBMCs) from newly-onset patients with microscopic polyangiitis (MPA). Increased proportions of activated CD14+ monocytes and CD14+ monocytes expressing interferon signature genes (ISGs) are distinctive features of MPA. Patient-specific analysis further classifies MPA into two groups. The MPA-MONO group is characterized by a high proportion of activated CD14+ monocytes, which persist before and after immunosuppressive therapy. These patients are clinically defined by increased monocyte ratio in the total PBMC count and have a high relapse rate. The MPA-IFN group is characterized by a high proportion of ISG+ CD14+ monocytes. These patients are clinically defined by high serum interferon-alpha concentrations and show good response to immunosuppressive therapy. Our findings identify the immunological phenotypes of MPA and provide clinical insights for personalized treatment and accurate prognostic prediction.

Abstract Atopic dermatitis (AD) is a skin disease that is heterogeneous both in terms of clinical manifestations and molecular profiles. It is increasingly recognized that AD is a systemic rather than a local disease and should be assessed in the context of whole-body pathophysiology. Here we show, via integrated RNA-sequencing of skin tissue and peripheral blood mononuclear cell (PBMC) samples along with clinical data from 115 AD patients and 14 matched healthy controls, that specific clinical presentations associate with matching differential molecular signatures. We establish a regression model based on transcriptome modules identified in weighted gene co-expression network analysis to extract molecular features associated with detailed clinical phenotypes of AD. The two main, qualitatively differential skin manifestations of AD, erythema and papulation are distinguished by differential immunological signatures. We further apply the regression model to a longitudinal dataset of 30 AD patients for personalized monitoring, highlighting patient heterogeneity in disease trajectories. The longitudinal features of blood tests and PBMC transcriptome modules identify three patient clusters which are aligned with clinical severity and reflect treatment history. Our approach thus serves as a framework for effective clinical investigation to gain a holistic view on the pathophysiology of complex human diseases.

Abstract Background Computed tomography (CT) imaging and artificial intelligence (AI)-based analyses have aided in the diagnosis and prediction of the severity of COVID-19. However, the potential of AI-based CT quantification of pneumonia in assessing patients with COVID-19 has not yet been fully explored. This study aimed to investigate the potential of AI-based CT quantification of COVID-19 pneumonia to predict the critical outcomes and clinical characteristics of patients with residual lung lesions. Methods This retrospective cohort study included 1,200 hospitalized patients with COVID-19 from four hospitals. The incidence of critical outcomes (requiring the support of high-flow oxygen or invasive mechanical ventilation or death) and complications during hospitalization (bacterial infection, renal failure, heart failure, thromboembolism, and liver dysfunction) was compared between the groups of pneumonia with high/low-percentage lung lesions, based on AI-based CT quantification. Additionally, 198 patients underwent CT scans 3 months after admission to analyze prognostic factors for residual lung lesions. Results The pneumonia group with a high percentage of lung lesions (N = 400) had a higher incidence of critical outcomes and complications during hospitalization than the low percentage group (N = 800). Multivariable analysis demonstrated that AI-based CT quantification of pneumonia was independently associated with critical outcomes (adjusted odds ratio [aOR] 10.5, 95% confidence interval [CI] 5.59–19.7), as well as with oxygen requirement (aOR 6.35, 95% CI 4.60–8.76), IMV requirement (aOR 7.73, 95% CI 2.52–23.7), and mortality rate (aOR 6.46, 95% CI 1.87–22.3). Among patients with follow-up CT scans (N = 198), the multivariable analysis revealed that the pneumonia group with a high percentage of lung lesions on admission (aOR 4.74, 95% CI 2.36–9.52), older age (aOR 2.53, 95% CI 1.16–5.51), female sex (aOR 2.41, 95% CI 1.13–5.11), and medical history of hypertension (aOR 2.22, 95% CI 1.09–4.50) independently predicted persistent residual lung lesions. Conclusions AI-based CT quantification of pneumonia provides valuable information beyond qualitative evaluation by physicians, enabling the prediction of critical outcomes and residual lung lesions in patients with COVID-19.

Objectives: Although SARS-CoV-2 RNAemia has been reported to strongly impact patients with severe COVID-19, the clinical characteristics of patients with COVID-19 harboring detectable intracellular SARS-CoV-2 RNA remain unknown. Methods: We included adult patients who had developed COVID-19 between February and September 2020. Total white blood cells derived from the buffy coat of peripheral whole blood were used to detect SARS-CoV-2 RNA using the Illumina COVIDSeq test. We compared the clinical characteristics between patients with and without detected viral RNA (detected and undetected groups). Results: Among the 390 patients included, 17 harbored SARS-CoV-2 RNA in peripheral white blood cells. All 17 patients required oxygen support during the disease course and had higher intensive care unit admission (52.9% vs 28.9%, P = 0.035), mortality (17.7% vs 3.5%, P = 0.004), kidney dysfunction (severe, 23.5% vs 6.4%, P = 0.029), and corticosteroid treatment rates (76.5% vs 46.5%, P = 0.016) than those of patients in the undetected group. Conclusion: We propose that patients with circulating intracellular SARS-CoV-2 RNA in the peripheral blood exhibited the most severe disease course.

Objectives: COVID-19 was severe in the Delta variant-dominated epidemic wave (fifth wave) in Japan. The clinical characteristics and effectiveness of COVID-19 vaccination are not fully understood in the Omicron variant-dominated wave (sixth and seventh waves), especially in hospitalized patients. We investigated the relationship between vaccination and disease severity in the Omicron-dominated wave and compared these variant-dominated waves. Methods: The nationwide COVID-19 database (Japan COVID-19 Task Force) was used to compare clinical characteristics and critical outcomes in patients hospitalized with Delta (fifth, N = 735) vs Omicron-dominated waves (sixth, N = 495; seventh, N = 128). Results: Patients in the sixth and seventh waves had a lower incidence of critical outcomes and respiratory outcomes, and a higher incidence of bacterial infection, although the mortality rate did not differ significantly between waves. In the sixth and seventh waves, 138 (27.9%) and 29 (22.7%) patients with COVID-19 were unvaccinated, respectively. Multivariable analysis adjusted with previously reported factors revealed that the proportion of (1) critical outcomes and (2) respiratory outcomes decreased in a frequency-dependent manner. Thus, (1) (the number of vaccinations): 1-2 times: adjusted odds ratio (aOR) 0.37 (95% confidence interval [CI]; 0.20-0.69); 3-4 times: aOR 0.25 (95% CI; 0.11-0.58); and (2) 1-2 times: aOR 0.43 (95% CI; 0.27-0.66); 3-4 times: aOR 0.36 (95% CI; 0.21-0.60). Conclusions: Patients hospitalized with COVID-19 with Omicron infections showed a lower incidence of critical outcomes than those with Delta infections, and COVID-19 vaccination may contribute to preventing respiratory failure.

Abstract Integrating genomic data of multiple cancers allows de novo cancer grouping and elucidating the shared genetic basis across cancers. Here, we conduct the pan-cancer and cross-population genome-wide association study (GWAS) meta-analysis and replication studies on 13 cancers including 250,015 East Asians (Biobank Japan) and 377,441 Europeans (UK Biobank). We identify ten cancer risk variants including five pleiotropic associations (e.g., rs2076295 at DSP on 6p24 associated with lung cancer and rs2525548 at TRIM4 on 7q22 nominally associated with six cancers). Quantifying shared heritability among the cancers detects positive genetic correlations between breast and prostate cancer across populations. Common genetic components increase the statistical power, and the large-scale meta-analysis of 277,896 breast/prostate cancer cases and 901,858 controls identifies 91 newly genome-wide significant loci. Enrichment analysis of pathways and cell types reveals shared genetic backgrounds across said cancers. Focusing on genetically correlated cancers can contribute to enhancing our insights into carcinogenesis.

Background Krebs von den Lungen-6 (KL-6) is a known biomarker for diagnosis and monitoring of interstitial lung diseases. However, the role of serum KL-6 and the mucin 1 (MUC1) variant (rs4072037) in COVID-19 outcomes remains to be elucidated. We aimed to evaluate the relationships among serum KL-6 levels, critical outcomes and the MUC1 variant in Japanese patients with COVID-19.Methods This is a secondary analysis of a multicentre retrospective study using data from the Japan COVID-19 Task Force collected from February 2020 to November 2021, including 2226 patients with COVID-19 whose serum KL-6 levels were measured. An optimal serum KL-6 level cut-off to predict critical outcomes was determined and used for multivariable logistic regression analysis. Furthermore, the relationship among the allele dosage of the MUC1 variant, calculated from single nucleotide polymorphism typing data of genome-wide association studies using the imputation method, serum KL-6 levels and COVID-19 critical outcomes was evaluated.Results Serum KL-6 levels were significantly higher in patients with COVID-19 with critical outcomes (511±442 U/mL) than those without (279±204 U/mL) (p

Abstract Background Although cases of respiratory bacterial infections associated with coronavirus disease 2019 (COVID-19) have often been reported, their impact on the clinical course remains unclear. Herein, we evaluated and analyzed the complication rates of bacterial infections, causative organisms, patient backgrounds, and clinical outcome in Japanese patients with COVID-19. Methods We performed a retrospective cohort study that included inpatients with COVID-19 from multiple centers participating in the Japan COVID-19 Taskforce (April 2020 to May 2021) and obtained demographic, epidemiological, and microbiological results and the clinical course and analyzed the cases of COVID-19 complicated by respiratory bacterial infections. Results Of the 1,863 patients with COVID-19 included in the analysis, 140 (7.5%) had respiratory bacterial infections. Community-acquired co-infection at COVID-19 diagnosis was uncommon (55/1,863, 3.0%) and was mainly caused by Staphylococcus aureus, Klebsiella pneumoniae and Streptococcus pneumoniae. Hospital-acquired bacterial secondary infections, mostly caused by Staphylococcus aureus, Pseudomonas aeruginosa, and Stenotrophomonas maltophilia, were diagnosed in 86 patients (4.6%). Severity-associated comorbidities were frequently observed in hospital-acquired secondary infection cases, including hypertension, diabetes, and chronic kidney disease. The study results suggest that the neutrophil–lymphocyte ratio (> 5.28) may be useful in diagnosing complications of respiratory bacterial infections. COVID-19 patients with community-acquired or hospital-acquired secondary infections had significantly increased mortality. Conclusions Respiratory bacterial co-infections and secondary infections are uncommon in patients with COVID-19 but may worsen outcomes. Assessment of bacterial complications is important in hospitalized patients with COVID-19, and the study findings are meaningful for the appropriate use of antimicrobial agents and management strategies.

Objectives: Smoking and chronic obstructive pulmonary disease (COPD) are risk factors for severe COVID-19. However, limited literature exists on the effect of COPD and smoking on COVID-19 outcomes. This study examined the impact of smoking exposure in pack-years (PY) and COPD on COVID-19 outcomes among smokers in Japan. Methods: The study included 1266 smokers enrolled by the Japan COVID-19 task force between February 2020 and December 2021. PY and COPD status was self-reported by patients. Patients were classified into the non-COPD (n = 1151) and COPD (n = 115) groups; the non-COPD group was further classified into 30 PY (n = 361). The study outcome was the need for invasive mechanical ventilation (IMV). Results: The incidence of IMV increased with increasing PY and was highest in the COPD group (30 PY = 15.2%, COPD = 26.1%; P 30 PY and COPD groups through univariate (odds ratio [OR]: >30 PY = 2.11, COPD = 4.14) and multivariate (OR: >30 PY = 2.38; COPD = 7.94) analyses. Increasing PY number was also associated with increased IMV requirement in patients aged

Summary: Interaction between the gut microbiome and host plays a key role in human health. Here, we perform a metagenome shotgun-sequencing-based analysis of Japanese participants to reveal associations between the gut microbiome, host genetics, and plasma metabolome. A genome-wide association study (GWAS) for microbial species (n = 524) identifies associations between the PDE1C gene locus and Bacteroides intestinalis and between TGIF2 and TGIF2-RAB5IF gene loci and Bacteroides acidifiaciens. In a microbial gene ortholog GWAS, agaE and agaS, which are related to the metabolism of carbohydrates forming the blood group A antigen, are associated with blood group A in a manner depending on the secretor status determined by the East Asian-specific FUT2 variant. A microbiome-metabolome association analysis (n = 261) identifies associations between bile acids and microbial features such as bile acid metabolism gene orthologs including bai and 7β-hydroxysteroid dehydrogenase. Our publicly available data will be a useful resource for understanding gut microbiome-host interactions in an underrepresented population.

Summary: Polygenic risk scores (PRSs) developed from multi-ancestry genome-wide association studies (GWASs), PRSmulti, hold promise for improving PRS accuracy and generalizability across populations. To establish best practices for leveraging the increasing diversity of genomic studies, we investigated how various factors affect the performance of PRSmulti compared with PRSs constructed from single-ancestry GWASs (PRSsingle). Through extensive simulations and empirical analyses, we showed that PRSmulti overall outperformed PRSsingle in understudied populations, except when the understudied population represented a small proportion of the multi-ancestry GWAS. Furthermore, integrating PRSs based on local ancestry-informed GWASs and large-scale, European-based PRSs improved predictive performance in understudied African populations, especially for less polygenic traits with large-effect ancestry-enriched variants. Our work highlights the importance of diversifying genomic studies to achieve equitable PRS performance across ancestral populations and provides guidance for developing PRSs from multiple studies.

A genome-wide association study in the BioBank Japan cohort reveals new locus, BET1L, which is associated with survival time in a Japanese population.

Tomofuji et al. perform non-targeted metabolomics on serum samples from people with type 2 diabetes with or without complications. The authors identify metabolites associated with retinal and renal complications.

Abstract Background Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery. Results To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3–5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism. Conclusions Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Abstract Background We aimed to elucidate differences in the characteristics of patients with coronavirus disease 2019 (COVID-19) requiring hospitalization in Japan, by COVID-19 waves, from conventional strains to the Delta variant. Methods We used secondary data from a database and performed a retrospective cohort study that included 3261 patients aged ≥ 18 years enrolled from 78 hospitals that participated in the Japan COVID-19 Task Force between February 2020 and September 2021. Results Patients hospitalized during the second (mean age, 53.2 years [standard deviation {SD}, ± 18.9]) and fifth (mean age, 50.7 years [SD ± 13.9]) COVID-19 waves had a lower mean age than those hospitalized during the other COVID-19 waves. Patients hospitalized during the first COVID-19 wave had a longer hospital stay (mean, 30.3 days [SD ± 21.5], p

Abstract GNE myopathy is a distal myopathy caused by biallelic variants in GNE, which encodes a protein involved in sialic acid biosynthesis. Compound heterozygosity of the second most frequent variant among Japanese GNE myopathy patients, GNE c.620A>T encoding p.D207V, occurs in the expected number of patients; however, homozygotes for this variant are rare; three patients identified while 238 homozygotes are estimated to exist in Japan. The aim of this study was to elucidate the pathomechanism caused by c.620A>T. Identity-by-descent mapping indicated two distinct c.620A>T haplotypes, which were not correlated with age onset or development of myopathy. Patients homozygous for c.620A>T had mildly decreased sialylation, and no additional pathogenic variants in GNE or abnormalities in transcript structure or expression of other genes related to sialic acid biosynthesis in skeletal muscle. Structural modeling of full-length GNE dimers revealed that the variant amino acid localized close to the monomer interface, but far from catalytic sites, suggesting functions in enzymatic product transfer between the epimerase and kinase domains on GNE oligomerization. In conclusion, homozygotes for c.620A>T rarely develop myopathy, while symptoms occur in compound heterozygotes, probably because of mildly decreased sialylation, due to partial defects in oligomerization and product trafficking by the mutated GNE protein.

Abstract Background The coronavirus disease 2019 (COVID-19) pandemic is widespread; however, accurate predictors of refractory cases have not yet been established. Circulating extracellular vesicles, involved in many pathological processes, are ideal resources for biomarker exploration. Methods To identify potential serum biomarkers and examine the proteins associated with the pathogenesis of refractory COVID-19, we conducted high-coverage proteomics on serum extracellular vesicles collected from 12 patients with COVID-19 at different disease severity levels and 4 healthy controls. Furthermore, single-cell RNA sequencing of peripheral blood mononuclear cells collected from 10 patients with COVID-19 and 5 healthy controls was performed. Results Among the 3046 extracellular vesicle proteins that were identified, expression of MACROH2A1 was significantly elevated in refractory cases compared to non-refractory cases; moreover, its expression was increased according to disease severity. In single-cell RNA sequencing of peripheral blood mononuclear cells, the expression of MACROH2A1 was localized to monocytes and elevated in critical cases. Consistently, single-nucleus RNA sequencing of lung tissues revealed that MACROH2A1 was highly expressed in monocytes and macrophages and was significantly elevated in fatal COVID-19. Moreover, molecular network analysis showed that pathways such as “estrogen signaling pathway,” “p160 steroid receptor coactivator (SRC) signaling pathway,” and “transcriptional regulation by STAT” were enriched in the transcriptome of monocytes in the peripheral blood mononuclear cells and lungs, and they were also commonly enriched in extracellular vesicle proteomics. Conclusions Our findings highlight that MACROH2A1 in extracellular vesicles is a potential biomarker of refractory COVID-19 and may reflect the pathogenesis of COVID-19 in monocytes.

Abstract Background Respiratory symptoms are associated with coronavirus disease 2019 (COVID-19) outcomes. However, the impacts of upper and lower respiratory symptoms on COVID-19 outcomes in the same population have not been compared. The objective of this study was to characterize upper and lower respiratory symptoms and compare their impacts on outcomes of hospitalized COVID-19 patients. Methods This was a multicenter, retrospective cohort study; the database from the Japan COVID-19 Task Force was used. A total of 3314 COVID-19 patients were included in the study, and the data on respiratory symptoms were collected. The participants were classified according to their respiratory symptoms (Group 1: no respiratory symptoms, Group 2: only upper respiratory symptoms, Group 3: only lower respiratory symptoms, and Group 4: both upper and lower respiratory symptoms). The impacts of upper and lower respiratory symptoms on the clinical outcomes were compared. The primary outcome was the percentage of patients with poor clinical outcomes, including the need for oxygen supplementation via high-flow oxygen therapy, mechanical ventilation, and extracorporeal membrane oxygenation or death. Results Of the 3314 COVID-19 patients, 605, 1331, 1229, and 1149 were classified as Group 1, Group 2, Group 3, and Group 4, respectively. In univariate analysis, patients in Group 2 had the best clinical outcomes among all groups (odds ratio [OR]: 0.21, 95% confidence interval [CI]: 0.11–0.39), while patients in Group 3 had the worst outcomes (OR: 3.27, 95% CI: 2.43–4.40). Group 3 patients had the highest incidence of pneumonia, other complications due to secondary infections, and thrombosis during the clinical course. Conclusions Upper and lower respiratory tract symptoms had vastly different impacts on the clinical outcomes of COVID-19.

Abstract Background The clinical course of coronavirus disease (COVID-19) is diverse, and the usefulness of phenotyping in predicting the severity or prognosis of the disease has been demonstrated overseas. This study aimed to investigate clinically meaningful phenotypes in Japanese COVID-19 patients using cluster analysis. Methods From April 2020 to May 2021, data from inpatients aged ≥ 18 years diagnosed with COVID-19 and who agreed to participate in the study were collected. A total of 1322 Japanese patients were included. Hierarchical cluster analysis was performed using variables reported to be associated with COVID-19 severity or prognosis, namely, age, sex, obesity, smoking history, hypertension, diabetes mellitus, malignancy, chronic obstructive pulmonary disease, hyperuricemia, cardiovascular disease, chronic liver disease, and chronic kidney disease. Results Participants were divided into four clusters: Cluster 1, young healthy (n = 266, 20.1%); Cluster 2, middle-aged (n = 245, 18.5%); Cluster 3, middle-aged obese (n = 435, 32.9%); and Cluster 4, elderly (n = 376, 28.4%). In Clusters 3 and 4, sore throat, dysosmia, and dysgeusia tended to be less frequent, while shortness of breath was more frequent. Serum lactate dehydrogenase, ferritin, KL-6, d-dimer, and C-reactive protein levels tended to be higher in Clusters 3 and 4. Although Cluster 3 had a similar age as Cluster 2, it tended to have poorer outcomes. Both Clusters 3 and 4 tended to exhibit higher rates of oxygen supplementation, intensive care unit admission, and mechanical ventilation, but the mortality rate tended to be lower in Cluster 3. Conclusions We have successfully performed the first phenotyping of COVID-19 patients in Japan, which is clinically useful in predicting important outcomes, despite the simplicity of the cluster analysis method that does not use complex variables.

Objectives: This study aimed to identify the relationship between abnormal serum uric acid levels or a history of hyperuricemia and COVID-19 severity in the Japanese population. Methods: We included 1523 patients enrolled in the Japan COVID-19 Task Force cohort between February 2020 and May 2021. We compared the clinical characteristics, including co-morbidities, laboratory findings, and outcomes, particularly invasive mechanical ventilation (IMV), among patients with and without abnormal uric acid levels or a history of hyperuricemia. Results: Patients with high serum uric acid levels were older and had higher body weight and body mass index than those without. In addition, the multiple logistic regression analysis revealed a significant association between high serum uric acid levels or a history of hyperuricemia and an increased risk of IMV (odds ratio [OR] = 1.77; P = 0.03/OR = 1.56; P = 0.04). Moreover, patients with low uric acid levels on admission were also associated significantly with the requirement of IMV (OR = 5.09; P

Purpose: The aim of this study was to understand the income and employment status of patients at the start of and during follow-up after palliative radiation therapy for bone metastasis. Methods and Materials: From December 2020 to March 2021, a prospective multi-institutional observational study was conducted to investigate income and employment of patients at the start of administration of radiation therapy for bone metastasis and at 2 and 6 months after treatment. Of 333 patients referred to radiation therapy for bone metastasis, 101 were not registered, mainly because of their poor general condition, and another 8 were excluded from the follow-up analysis owing to ineligibility. Results: In 224 patients analyzed, 108 had retired for reasons unrelated to cancer, 43 had retired for reasons related to cancer, 31 were taking leave, and 2 had lost their jobs at the time of registration. The number of patients who were in the working group was 40 (30 with no change in income and 10 with decreased income) at registration, 35 at 2 months, and 24 at 6 months. Younger patients (P = 0), patients with better performance status (P = 0), patients who were ambulatory (P = .008), and patients with lower scores on a numerical rating scale of pain (P = 0) were significantly more likely to be in the working group at registration. There were 9 patients who experienced improvements in their working status or income at least once in the follow-up after radiation therapy. Conclusions: The majority of patients with bone metastasis were not working at the start of or after radiation therapy, but the number of patients who were working was not negligible. Radiation oncologists should be aware of the working status of patients and provide appropriate support for each patient. The benefit of radiation therapy to support patients continuing their work and returning to work should be investigated further in prospective studies.

Intracranial germ cell tumors (IGCTs) are rare brain tumors mainly diagnosed in children and young adults. Here, the authors conduct a genome-wide association study for IGCTs, identify a risk locus at BAK1, and characterize its functional consequences.

A genome-wide association study identifies MFGE8 as protective against coronary atherosclerosis in European and East Asian populations.

Abstract Background Obesity is reported to be a risk factor for severe disease in patients with coronavirus disease 2019 (COVID-19). However, there are no specific reports on the risk of severe disease according to body mass index (BMI) in Japan. Thus, this study aimed to investigate the effect of obesity stratified by BMI on the severity of COVID-19 in the general Japanese population. Methods From February 2020 to May 2021, 1 837 patients aged ≥18 years were enrolled in the Japan COVID-19 Task Force. Patients with known BMI and disease severity were analyzed. Severity was defined as critical if the patient was treated in the intensive care unit, required invasive mechanical ventilation, or died. Results Class 1 obesity (25.0 ≤ BMI

Genetic mechanisms influencing COVID-19 susceptibility are not well understood. Here, the authors analyzed whole blood RNA-seq data of 465 Japanese individuals with COVID-19, highlighting thousands of fine-mapped variants affecting expression and splicing of genes, as well as the presence of COVID-19 severity-interaction eQTLs.

A large-scale GWAS provides insight on diabetes-dependent genetic effects on the glomerular filtration rate, a common metric to monitor kidney health in disease.

Objective The effects of coffee consumption on serum uric acid (SUA) levels and gout risk are controversial. There have hitherto been no reports based on Mendelian randomization (MR) analysis of its effects that consider pleiotropy. Here, we evaluated the effects of coffee consumption across ancestry populations, taking pleiotropy into account. Methods We performed the first MR analyses for coffee consumption on SUA levels and gout, considering pleiotropy. We used the following summary statistics of genome‐wide association studies from a Japanese population: habitual coffee consumption (152,634 subjects), gout (3053 cases and 4554 controls), and SUA levels (121,745 subjects). In addition to fixed‐effect inverse variance weighted (IVW) meta‐analysis, we performed a robust evaluation of heterogeneity and removed several instruments for reasons of possible pleiotropy. Previous European datasets were also reevaluated while heterogeneity was considered. Results Habitual coffee consumption was significantly and inversely associated with gout (odds ratio [OR] = 0.29, 95% confidence interval [95% CI]: 0.16‐0.51, P = 1.9 × 10−5) in random‐effect IVW (Phet = 5.5 × 10−19). Excluding pleiotropic instruments, the protective effect on gout was confirmed in fixed‐effect IVW analysis (OR = 0.75, 95% CI: 0.58‐0.97, P = 0.026) without heterogeneity (Phet = 0.39). However, we observed no significance in the previous European datasets when heterogeneity was considered. Associations were not observed between coffee consumption and SUA levels in either ancestry in MR analyses that considered pleiotropy. Multivariable MR analysis showed that increased coffee consumption significantly reduced gout risk, even after adjusting for SUA levels (OR = 0.50, 95% CI: 0.31‐0.81, P = 0.0046). Conclusion With pleiotropy taken into account, our MR analyses revealed that coffee consumption can causally reduce gout risk, and that it may reduce gout risk independently of SUA levels.

Background: Allergic diseases are some of the most common diseases worldwide. Genome-wide association studies (GWASs) have been conducted to elucidate the genetic factors of allergic diseases. However, no GWASs for allergen component sensitization have been performed. Objective: We sought to detect genetic variants associated with differences in immune responsiveness against allergen components. Methods: The participants of the present study were recruited from the Tokyo Children’s Health, Illness, and Development study, and allergen component–specific IgE level at age 9 years was measured by means of allergen microarray immunoassays. We performed GWASs for allergen component sensitization against each allergen (single allergen component sensitization, number of allergen components analyzed, n = 31), as well as against allergen protein families (allergen protein group sensitization, number of protein groups analyzed, n = 16). Results: We performed GWAS on 564 participants of the Tokyo Children’s Health, Illness, and Development study and found associations between Amb a 1 sensitization and the immunoglobulin heavy-chain variable gene on chromosome 14 and between Phl p 1 sensitization and the HLA class II region on chromosome 6 (P < 5.0 × 10−8). A GWAS-significant association was also observed between the HLA class II region and profilin sensitization (P < 5.0 × 10−8). Conclusions: Our data provide the first demonstration of genetic risk for allergen component sensitization and show that this genetic risk is related to immune response genes including immunoglobulin heavy-chain variable gene and HLA.

Summary: Polygenic risk scores (PRSs) have been widely explored in precision medicine. However, few studies have thoroughly investigated their best practices in global populations across different diseases. We here utilized data from Global Biobank Meta-analysis Initiative (GBMI) to explore methodological considerations and PRS performance in 9 different biobanks for 14 disease endpoints. Specifically, we constructed PRSs using pruning and thresholding (P + T) and PRS-continuous shrinkage (CS). For both methods, using a European-based linkage disequilibrium (LD) reference panel resulted in comparable or higher prediction accuracy compared with several other non-European-based panels. PRS-CS overall outperformed the classic P + T method, especially for endpoints with higher SNP-based heritability. Notably, prediction accuracy is heterogeneous across endpoints, biobanks, and ancestries, especially for asthma, which has known variation in disease prevalence across populations. Overall, we provide lessons for PRS construction, evaluation, and interpretation using GBMI resources and highlight the importance of best practices for PRS in the biobank-scale genomics era.

Background and design: The coronavirus disease (COVID-19) pandemic is having a devastating effect worldwide. Host genome differences between populations may influence the severity of COVID-19.The Japan COVID-19 Task Force is conducting host genome analysis of hospitalized patients with COVID-19 from more than 70 institutions nationwide in Japan. This report describes the clinical characteristics of patients enrolled to date. Results: The median (interquartile range) age of the 1674 patients included in the analysis was 59 (45–71) years, and more than half of the patients (66.2%) were male. Less than half of the patients (41.2%) had severe disease. The case fatality rate was 3.2%. Conclusions: Since this is a hospital-based study, the number of severe cases was relatively high, but the case fatality rate was relatively low, when compared to that of other countries. In the future, we will continue to enroll patients and conduct genome analyses of patients with COVID-19.

Abstract Detailed analysis of the cells that infiltrate lesional skin cannot be performed in skin biopsy specimens using immunohistochemistry or cell separation techniques because enzyme treatments applied during the isolation step can destroy small amounts of protein and minor cell populations in the biopsy specimen. Here, we describe a method for isolating T cells from drops of whole blood obtained from lesions during skin biopsy in patients with cutaneous T-cell lymphoma. Lesional blood is assumed to contain lesional resident cells, cells from capillary vessels, and blood overflowing from capillary vessels into the lesion area. The lesional blood showed substantial increases in distinct cell populations, chemokines, and the expression of various genes. The proportion of CD8+CD45RO+ T cells in the lesional blood negatively correlated with the modified severity-weighted assessment tool scores. CD4+CD45RO+ T cells in the lesional blood expressed genes associated with the development of cancer and progression of cutaneous T-cell lymphoma. In addition, CD8+CD45RO+ T cells in lesional blood had unique T-cell receptor repertoires in lesions of each stage. Assessment of lesional blood drops might provide new insight into the pathogenesis of mycosis fungoides and facilitate evaluation of the treatment efficacy for mycosis fungoides as well as other skin inflammatory diseases.

Consecutive mRNA vaccinations against SARS-CoV-2 reinforced both innate and adaptive immune responses. However, it remains unclear whether the enhanced innate immune responses are mediated by epigenetic regulation and, if so, whether these effects persist. Using mass cytometry, RNA-Seq, and ATAC-Seq, we show that BNT162b2 mRNA vaccination upregulated antiviral and IFN-stimulated gene expression in monocytes with greater effects after the second vaccination than those after the first vaccination. Transcription factor–binding motif analysis also revealed enriched IFN regulatory factors and PU.1 motifs in accessible chromatin regions. Importantly, although consecutive BNT162b2 mRNA vaccinations boosted innate immune responses and caused epigenetic changes in isolated monocytes, we show that these effects occurred only transiently and disappeared 4 weeks after the second vaccination. Furthermore, single-cell RNA-Seq analysis revealed that a similar gene signature was impaired in the monocytes of unvaccinated patients with COVID-19 with acute respiratory distress syndrome. These results reinforce the importance of the innate immune response in the determination of COVID-19 severity but indicate that, unlike adaptive immunity, innate immunity is not unexpectedly sustained even after consecutive vaccination. This study, which focuses on innate immune memory, may provide novel insights into the vaccine development against infectious diseases.

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

Abstract The prediction of disease risks is an essential part of personalized medicine, which includes early disease detection, prevention, and intervention. The polygenic risk score (PRS) has become the standard for quantifying genetic liability in predicting disease risks. PRS utilizes single-nucleotide polymorphisms (SNPs) with genetic risks elucidated by genome-wide association studies (GWASs) and is calculated as weighted sum scores of these SNPs with genetic risks using their effect sizes from GWASs as their weights. The utilities of PRS have been explored in many common diseases, such as cancer, coronary artery disease, obesity, and diabetes, and in various non-disease traits, such as clinical biomarkers. These applications demonstrated that PRS could identify a high-risk subgroup of these diseases as a predictive biomarker and provide information on modifiable risk factors driving health outcomes. On the other hand, there are several limitations to implementing PRSs in clinical practice, such as biased sensitivity for the ethnic background of PRS calculation and geographical differences even in the same population groups. Also, it remains unclear which method is the most suitable for the prediction with high accuracy among numerous PRS methods developed so far. Although further improvements of its comprehensiveness and generalizability will be needed for its clinical implementation in the future, PRS will be a powerful tool for therapeutic interventions and lifestyle recommendations in a wide range of diseases. Thus, it may ultimately improve the health of an entire population in the future.

Finding causal variants and genes from GWAS loci results remains a challenge. Here, the authors train a model to predict if a variant affects nearby gene expression, and apply the method to identify new possible causal eQTLs and mechanisms of GWAS loci.