Body mass index stratification optimizes polygenic prediction of type 2 diabetes in cross-biobank analyses.

Type 2 diabetes (T2D) shows heterogeneous body mass index (BMI) sensitivity. Here, we performed stratification based on BMI to optimize predictions for BMI-related diseases. We obtained BMI-stratified datasets using data from more than 195,000 individuals (nT2D = 55,284) from BioBank Japan (BBJ) and UK Biobank. T2D heritability in the low-BMI group was greater than that in the high-BMI group. Polygenic predictions of T2D toward low-BMI targets had pseudo-R2 values that were more than 22% higher than BMI-unstratified targets. Polygenic risk scores (PRSs) from low-BMI discovery outperformed PRSs from high BMI, while PRSs from BMI-unstratified discovery performed best. Pathway-specific PRSs demonstrated the biological contributions of pathogenic pathways. Low-BMI T2D cases showed higher rates of neuropathy and retinopathy. Combining BMI stratification and a method integrating cross-population effects, T2D predictions showed greater than 37% improvements over unstratified-matched-population prediction. We replicated findings in the Tohoku Medical Megabank (n = 26,000) and the second BBJ cohort (n = 33,096). Our findings suggest that target stratification based on existing traits can improve the polygenic prediction of heterogeneous diseases.

Brugada syndrome in Japan and Europe: a genome-wide association study reveals shared genetic architecture and new risk loci.

BACKGROUND AND AIMS: Brugada syndrome (BrS) is an inherited arrhythmia with a higher disease prevalence and more lethal arrhythmic events in Asians than in Europeans. Genome-wide association studies (GWAS) have revealed its polygenic architecture mainly in European populations. The aim of this study was to identify novel BrS-associated loci and to compare allelic effects across ancestries. METHODS: A GWAS was conducted in Japanese participants, involving 940 cases and 1634 controls, followed by a cross-ancestry meta-analysis of Japanese and European GWAS (total of 3760 cases and 11 635 controls). The novel loci were characterized by fine-mapping, gene expression, and splicing quantitative trait associations in the human heart. RESULTS: The Japanese-specific GWAS identified one novel locus near ZSCAN20 (P = 1.0 × 10-8), and the cross-ancestry meta-analysis identified 17 association signals, including six novel loci. The effect directions of the 17 lead variants were consistent (94.1%; P for sign test = 2.7 × 10-4), and their allelic effects were highly correlated across ancestries (Pearson's R = .91; P = 2.9 × 10-7). The genetic risk score derived from the BrS GWAS of European ancestry was significantly associated with the risk of BrS in the Japanese population [odds ratio 2.12 (95% confidence interval 1.94-2.31); P = 1.2 × 10-61], suggesting a shared genetic architecture across ancestries. Functional characterization revealed that a lead variant in CAMK2D promotes alternative splicing, resulting in an isoform switch of calmodulin kinase II-δ, favouring a pro-inflammatory/pro-death pathway. CONCLUSIONS: This study demonstrates novel susceptibility loci implicating potentially novel pathogenesis underlying BrS. Despite differences in clinical expressivity and epidemiology, the polygenic architecture of BrS was substantially shared across ancestries.

Genetic drivers of heterogeneity in type 2 diabetes pathophysiology.

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.

Extracting immunological and clinical heterogeneity across autoimmune rheumatic diseases by cohort-wide immunophenotyping.

OBJECTIVE: Extracting immunological and clinical heterogeneity across autoimmune rheumatic diseases (AIRDs) is essential towards personalised medicine. METHODS: We conducted large-scale and cohort-wide immunophenotyping of 46 peripheral immune cells using Human Immunology Protocol of comprehensive 8-colour flow cytometric analysis. Dataset consisted of >1000 Japanese patients of 11 AIRDs with deep clinical information registered at the FLOW study, including rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE). In-depth clinical and immunological characterisation was conducted for the identified RA patient clusters, including associations of inborn human genetics represented by Polygenic Risk Score (PRS). RESULTS: Multimodal clustering of immunophenotypes deciphered underlying disease-cell type network in immune cell, disease and patient cluster resolutions. This provided immune cell type specificity shared or distinct across AIRDs, such as close immunological network between mixed connective tissue disease and SLE. Individual patient-level clustering dissected patients with AIRD into several clusters with different immunological features. Of these, RA-like or SLE-like clusters were exclusively dominant, showing immunological differentiation between RA and SLE across AIRDs. In-depth clinical analysis of RA revealed that such patient clusters differentially defined clinical heterogeneity in disease activity and treatment responses, such as treatment resistance in patients with RA with SLE-like immunophenotypes. PRS based on RA case-control and within-case stratified genome-wide association studies were associated with clinical and immunological characteristics. This pointed immune cell type implicated in disease biology such as dendritic cells for RA-interstitial lung disease. CONCLUSION: Cohort-wide and cross-disease immunophenotyping elucidate clinically heterogeneous patient subtypes existing within single disease in immune cell type-specific manner.

Analysis of gut microbiome, host genetics, and plasma metabolites reveals gut microbiome-host interactions in the Japanese population.

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.

Estimating gene-level false discovery probability improves eQTL statistical fine-mapping precision.

Statistical fine-mapping prioritizes putative causal variants from a large number of candidate variants, and is widely used in expression quantitative loci (eQTLs) studies. In eQTL fine-mapping, the existence of causal variants for gene expression is not guaranteed, since the genetic heritability of gene expression explained by nearby (cis-) variants is limited. Here we introduce a refined fine-mapping algorithm, named Knockoff-Finemap combination (KFc). KFc estimates the probability that the causal variant(s) exist in the cis-window of a gene through construction of knockoff genotypes (i.e. a set of synthetic genotypes that resembles the original genotypes), and uses it to adjust the posterior inclusion probabilities (PIPs). Utilizing simulated gene expression data, we show that KFc results in calibrated PIP distribution with improved precision. When applied to gene expression data of 465 genotyped samples from the Japan COVID-19 Task Force (JCTF), KFc resulted in significant enrichment of a functional score as well as reporter assay hits in the top PIP bins. When combined with functional priors derived from an external fine-mapping study (GTEx), KFc resulted in a significantly higher proportion of hematopoietic trait putative causal variants in the top PIP bins. Our work presents improvements in the precision of a major fine-mapping algorithm.

A genome-wide association study for allergen component sensitizations identifies allergen component-specific and allergen protein group-specific associations.

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.

Genome-wide association study identifies risk loci within the major histocompatibility complex region for Hunner-type interstitial cystitis.

Hunner-type interstitial cystitis (HIC) is a rare, chronic inflammatory disease of the urinary bladder with unknown etiology and genetic background. Here, we conduct a genome-wide association study of 144 patients with HIC and 41,516 controls of Japanese ancestry. The genetic variant, rs1794275, in the major histocompatibility complex (MHC) region (chromosome 6p21.3) is associated with HIC risk (odds ratio [OR] = 2.32; p = 3.4 × 10-9). The association is confirmed in a replication set of 26 cases and 1,026 controls (p = 0.014). Fine mapping demonstrates the contribution to the disease risk of a completely linked haplotype of three human leukocyte antigen HLA-DQß1 amino acid positions, 71, 74, and 75 (OR = 1.94; p = 5 × 10-8) and of HLA-DPß1 amino acid position 178, which tags HLA-DPB1∗04:02 (OR = 2.35; p = 7.5 × 10-8). The three HLA-DQß1 amino acid positions are located together at the peptide binding groove, suggesting their functional importance in antigen presentation. Our study reveals genetic contributions to HIC risk that may be associated with class II MHC molecule antigen presentation.

Pan-cancer and cross-population genome-wide association studies dissect shared genetic backgrounds underlying carcinogenesis.

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.

Reconstruction of the personal information from human genome reads in gut metagenome sequencing data.

Human DNA present in faecal samples can result in a small number of human reads in gut shotgun metagenomic sequencing data. However, it is presently unclear how much personal information can be reconstructed from such reads, and this has not been quantitatively evaluated. Such a quantitative evaluation is necessary to clarify the ethical concerns related to data sharing and to enable efficient use of human genetic information in stool samples, such as for research and forensics. Here we used genomic approaches to reconstruct personal information from the faecal metagenomes of 343 Japanese individuals with associated human genotype data. Genetic sex could be accurately predicted based on the sequencing depth of sex chromosomes for 97.3% of the samples. Individuals could be re-identified from the matched genotype data based on human reads recovered from the faecal metagenomic data with 93.3% sensitivity using a likelihood score-based method. This method also enabled us to predict the ancestries of 98.3% of the samples. Finally, we performed ultra-deep shotgun metagenomic sequencing of five faecal samples as well as whole-genome sequencing of blood samples. Using genotype-calling approaches, we demonstrated that the genotypes of both common and rare variants could be reconstructed from faecal samples. This included clinically relevant variants. Our approach can be used to quantify personal information contained within gut metagenome data.

Multi-ancestry genome-wide study in >2.5 million individuals reveals heterogeneity in mechanistic pathways of type 2 diabetes and complications.

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes. To characterise the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study (GWAS) data from 2,535,601 individuals (39.7% non-European ancestry), including 428,452 T2D cases. We identify 1,289 independent association signals at genome-wide significance (P<5×10-8) that map to 611 loci, of which 145 loci are previously unreported. We define eight non-overlapping clusters of T2D signals characterised by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial, and enteroendocrine cells. We build cluster-specific partitioned genetic risk scores (GRS) in an additional 137,559 individuals of diverse ancestry, including 10,159 T2D cases, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned GRS are more strongly associated with coronary artery disease and end-stage diabetic nephropathy than an overall T2D GRS across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings demonstrate the value of integrating multi-ancestry GWAS with single-cell epigenomics to disentangle the aetiological heterogeneity driving the development and progression of T2D, which may offer a route to optimise global access to genetically-informed diabetes care.

Single-cell analyses and host genetics highlight the role of innate immune cells in COVID-19 severity.

Mechanisms underpinning the dysfunctional immune response in severe acute respiratory syndrome coronavirus 2 infection are elusive. We analyzed single-cell transcriptomes and T and B cell receptors (BCR) of >895,000 peripheral blood mononuclear cells from 73 coronavirus disease 2019 (COVID-19) patients and 75 healthy controls of Japanese ancestry with host genetic data. COVID-19 patients showed a low fraction of nonclassical monocytes (ncMono). We report downregulated cell transitions from classical monocytes to ncMono in COVID-19 with reduced CXCL10 expression in ncMono in severe disease. Cell-cell communication analysis inferred decreased cellular interactions involving ncMono in severe COVID-19. Clonal expansions of BCR were evident in the plasmablasts of patients. Putative disease genes identified by COVID-19 genome-wide association study showed cell type-specific expressions in monocytes and dendritic cells. A COVID-19-associated risk variant at the IFNAR2 locus (rs13050728) had context-specific and monocyte-specific expression quantitative trait loci effects. Our study highlights biological and host genetic involvement of innate immune cells in COVID-19 severity.

Genetic Risk of Primary Aldosteronism and Its Contribution to Hypertension: A Cross-Ancestry Meta-Analysis of Genome-Wide Association Studies.

BACKGROUND: Hypertension imposes substantial health and economic burden worldwide. Primary aldosteronism (PA) is one of the most common causes of secondary hypertension, causing cardiovascular events at higher risk compared with essential hypertension. However, the germline genetic contribution to the susceptibility of PA has not been well elucidated. METHOD: We conducted a genome-wide association analysis of PA in the Japanese population and a cross-ancestry meta-analysis combined with UK Biobank and FinnGen cohorts (816 PA cases and 425 239 controls) to identify genetic variants that contribute to PA susceptibility. We also performed a comparative analysis for the risk of 42 previously established blood pressure-associated variants between PA and hypertension with the adjustment of blood pressure. RESULTS: In the Japanese genome-wide association study, we identified 10 loci that presented suggestive evidence for the association with the PA risk (P<1.0×10-6). In the meta-analysis, we identified 5 genome-wide significant loci (1p13, 7p15, 11p15, 12q24, and 13q12; P<5.0×10-8), including 3 of the suggested loci in the Japanese genome-wide association study. The strongest association was observed at rs3790604 (1p13), an intronic variant of WNT2B (odds ratio, 1.50 [95% CI, 1.33-1.69]; P=5.2×10-11). We further identified 1 nearly genome-wide significant locus (8q24, CYP11B2), which presented a significant association in the gene-based test (P=7.2×10-7). Of interest, all of these loci were known to be associated with blood pressure in previous studies, presumably because of the prevalence of PA among individuals with hypertension. This assumption was supported by the observation that they had a significantly higher risk effect on PA than on hypertension. We also revealed that 66.7% of the previously established blood pressure-associated variants had a higher risk effect for PA than for hypertension. CONCLUSIONS: This study demonstrates the genome-wide evidence for a genetic predisposition to PA susceptibility in the cross-ancestry cohorts and its significant contribution to the genetic background of hypertension. The strongest association with the WNT2B variants reinforces the implication of the Wnt/ß-catenin pathway in the PA pathogenesis.

Genome-Wide Association Study of Intracranial Artery Stenosis Followed by Phenome-Wide Association Study.

The genetic background of intracranial artery stenosis (ICAS), a major cause of ischemic stroke, remains elusive. We performed the world's first genome-wide association study (GWAS) of ICAS using DNA samples from Japanese subjects, to identify the genetic factors associated with ICAS and their correlation with clinical features. We also conducted a phenome-wide association study (PheWAS) of the top variant identified via GWAS to determine its association with systemic disease. The GWAS involved 408 patients with ICAS and 349 healthy controls and utilized an Asian Screening Array of venous blood samples. The PheWAS was performed using genotypic and phenotypic data of the Biobank Japan Project, which contained information on 46 diseases and 60 quantitative trait data from > 150,000 Japanese individuals. The GWAS revealed that the East Asian-specific functional variant of RNF213, rs112735431 (c.14429G > A, p.Arg4810Lys), was associated with ICAS (odds ratio, 12.3; 95% CI 5.5 to 27.5; P = 7.8 × 10-10). Stratified analysis within ICAS cases demonstrated that clinical features of those with and without the risk allele were different. PheWAS indicated that high blood pressure and angina were significantly associated with RNF213 rs112735431. The first GWAS of ICAS, which stratifies subpopulations within the ICAS cases with distinct clinical features, revealed that RNF213 rs112735431 was the most significant variant associated with ICAS. Thus, RNF213 rs112735431 shows potential as an important clinical biomarker that characterizes pleiotropic risk in various vascular diseases, such as blood pressure and angina, thereby facilitating personalized medicine for systemic vascular diseases in East Asian populations.

Genetic footprints of assortative mating in the Japanese population.

Assortative mating (AM) is a pattern characterized by phenotypic similarities between mating partners. Detecting the evidence of AM has been challenging due to the lack of large-scale datasets that include phenotypic data on both partners, especially in populations of non-European ancestries. Gametic phase disequilibrium between trait-associated alleles is a signature of parental AM on a polygenic trait, which can be detected even without partner data. Here, using polygenic scores for 81 traits in the Japanese population using BioBank Japan Project genome-wide association studies data (n = 172,270), we found evidence of AM on the liability to type 2 diabetes and coronary artery disease, as well as on dietary habits. In cross-population comparison using United Kingdom Biobank data (n = 337,139) we found shared but heterogeneous impacts of AM between populations.

Multidimensional analyses of the pathomechanism caused by the non-catalytic GNE variant, c.620A>T, in patients with GNE myopathy.

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.

Multi-ancestry genome-wide association analyses identify novel genetic mechanisms in rheumatoid arthritis.

Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 × 10-8), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA.

The whole blood transcriptional regulation landscape in 465 COVID-19 infected samples from Japan COVID-19 Task Force.

Coronavirus disease 2019 (COVID-19) is a recently-emerged infectious disease that has caused millions of deaths, where comprehensive understanding of disease mechanisms is still unestablished. In particular, studies of gene expression dynamics and regulation landscape in COVID-19 infected individuals are limited. Here, we report on a thorough analysis of whole blood RNA-seq data from 465 genotyped samples from the Japan COVID-19 Task Force, including 359 severe and 106 non-severe COVID-19 cases. We discover 1169 putative causal expression quantitative trait loci (eQTLs) including 34 possible colocalizations with biobank fine-mapping results of hematopoietic traits in a Japanese population, 1549 putative causal splice QTLs (sQTLs; e.g. two independent sQTLs at TOR1AIP1), as well as biologically interpretable trans-eQTL examples (e.g., REST and STING1), all fine-mapped at single variant resolution. We perform differential gene expression analysis to elucidate 198 genes with increased expression in severe COVID-19 cases and enriched for innate immune-related functions. Finally, we evaluate the limited but non-zero effect of COVID-19 phenotype on eQTL discovery, and highlight the presence of COVID-19 severity-interaction eQTLs (ieQTLs; e.g., CLEC4C and MYBL2). Our study provides a comprehensive catalog of whole blood regulatory variants in Japanese, as well as a reference for transcriptional landscapes in response to COVID-19 infection.

DOCK2 is involved in the host genetics and biology of severe COVID-19.

Identifying the host genetic factors underlying severe COVID-19 is an emerging challenge1-5. Here we conducted a genome-wide association study (GWAS) involving 2,393 cases of COVID-19 in a cohort of Japanese individuals collected during the initial waves of the pandemic, with 3,289 unaffected controls. We identified a variant on chromosome 5 at 5q35 (rs60200309-A), close to the dedicator of cytokinesis 2 gene (DOCK2), which was associated with severe COVID-19 in patients less than 65 years of age. This risk allele was prevalent in East Asian individuals but rare in Europeans, highlighting the value of genome-wide association studies in non-European populations. RNA-sequencing analysis of 473 bulk peripheral blood samples identified decreased expression of DOCK2 associated with the risk allele in these younger patients. DOCK2 expression was suppressed in patients with severe cases of COVID-19. Single-cell RNA-sequencing analysis (n = 61 individuals) identified cell-type-specific downregulation of DOCK2 and a COVID-19-specific decreasing effect of the risk allele on DOCK2 expression in non-classical monocytes. Immunohistochemistry of lung specimens from patients with severe COVID-19 pneumonia showed suppressed DOCK2 expression. Moreover, inhibition of DOCK2 function with CPYPP increased the severity of pneumonia in a Syrian hamster model of SARS-CoV-2 infection, characterized by weight loss, lung oedema, enhanced viral loads, impaired macrophage recruitment and dysregulated type I interferon responses. We conclude that DOCK2 has an important role in the host immune response to SARS-CoV-2 infection and the development of severe COVID-19, and could be further explored as a potential biomarker and/or therapeutic target.

A common deletion at BAK1 reduces enhancer activity and confers risk of intracranial germ cell tumors.

Intracranial germ cell tumors (IGCTs) are rare brain neoplasms that mainly occur in children and adolescents with a particularly high incidence in East Asian populations. Here, we conduct a genome-wide association study (GWAS) of 133 patients with IGCTs and 762 controls of Japanese ancestry. A common 4-bp deletion polymorphism in an enhancer adjacent to BAK1 is significantly associated with the disease risk (rs3831846; P = 2.4 × 10-9, odds ratio = 2.46 [95% CI: 1.83-3.31], minor allele frequency = 0.43). Rs3831846 is in strong linkage disequilibrium with a testicular GCTs susceptibility variant rs210138. In-vitro reporter assays reveal rs3831846 to be a functional variant attenuating the enhancer activity, suggesting its contribution to IGCTs predisposition through altering BAK1 expression. Risk alleles of testicular GCTs derived from the European GWAS show significant positive correlations in the effect sizes with the Japanese IGCTs GWAS (P = 1.3 × 10-4, Spearman's ρ = 0.48). These results suggest the shared genetic susceptibility of GCTs beyond ethnicity and primary sites.