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 sex-biased imbalance between Tfr, Tph, and atypical B cells determines antibody responses in COVID-19 patients.

Sex-biased humoral immune responses to COVID-19 patients have been observed, but the cellular basis for this is not understood. Using single-cell proteomics by mass cytometry, we find disrupted regulation of humoral immunity in COVID-19 patients, with a sex-biased loss of circulating follicular regulatory T cells (cTfr) at a significantly greater rate in male patients. In addition, a male sex-associated cellular network of T-peripheral helper, plasma blasts, proliferating and extrafollicular/atypical CD11c+ memory B cells was strongly positively correlated with neutralizing antibody concentrations and negatively correlated with cTfr frequency. These results suggest that sex-specific differences to the balance of cTfr and a network of extrafollicular antibody production-associated cell types may be a key factor in the altered humoral immune responses between male and female COVID-19 patients.

Galectin-10 in serum extracellular vesicles reflects asthma pathophysiology.

BACKGROUND: Novel biomarkers (BMs) are urgently needed for bronchial asthma (BA) with various phenotypes and endotypes. OBJECTIVE: We sought to identify novel BMs reflecting tissue pathology from serum extracellular vesicles (EVs). METHODS: We performed data-independent acquisition of serum EVs from 4 healthy controls, 4 noneosinophilic asthma (NEA) patients, and 4 eosinophilic asthma (EA) patients to identify novel BMs for BA. We confirmed EA-specific BMs via data-independent acquisition validation in 61 BA patients and 23 controls. To further validate these findings, we performed data-independent acquisition for 6 patients with chronic rhinosinusitis without nasal polyps and 7 patients with chronic rhinosinusitis with nasal polyps. RESULTS: We identified 3032 proteins, 23 of which exhibited differential expression in EA. Ingenuity pathway analysis revealed that protein signatures from each phenotype reflected disease characteristics. Validation revealed 5 EA-specific BMs, including galectin-10 (Gal10), eosinophil peroxidase, major basic protein, eosinophil-derived neurotoxin, and arachidonate 15-lipoxygenase. The potential of Gal10 in EVs was superior to that of eosinophils in terms of diagnostic capability and detection of airway obstruction. In rhinosinusitis patients, 1752 and 8413 proteins were identified from EVs and tissues, respectively. Among 11 BMs identified in EVs and tissues from patients with chronic rhinosinusitis with nasal polyps, 5 (including Gal10 and eosinophil peroxidase) showed significant correlations between EVs and tissues. Gal10 release from EVs was implicated in eosinophil extracellular trapped cell death in vitro and in vivo. CONCLUSION: Novel BMs such as Gal10 from serum EVs reflect disease pathophysiology in BA and may represent a new target for liquid biopsy approaches.

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.

Clinical and molecular significance of flow cytometric analysis for reactive oxygen species production and residual p67phox expression in p67phox-deficient chronic granulomatous disease.

Chronic granulomatous disease (CGD) is a primary immunodeficiency disease caused by molecular defects in nicotinamide adenine dinucleotide phosphate (NADPH) oxidase. p67phox-CGD is an autosomal recessive CGD, which is caused by a defect in the cytosolic components of NADPH oxidase, p67phox, encoded by NCF2. We previously established a flow cytometric analysis for p67phox expression, which allows accurate assessment of residual protein expression in p67phox-CGD. We evaluated the correlation between oxidase function and p67phox expression, and assessed the relevancy to genotypes and clinical phenotypes in 11 patients with p67phox-CGD. Reactive oxygen species (ROS) production by granulocytes was evaluated using dihydrorhodamine-1,2,3 (DHR) assays. p67phox expression was evaluated in the monocyte population. DHR activity and p67phox expression were significantly correlated (r = 0.718, p < 0.0162). Additionally, DHR activity and p67phox expression were significantly higher in patients carrying one missense variant in combination with one nonsense or frameshift variant in the NCF2 gene than in patients with only null variants. The available clinical parameters of our patients (i.e., age at disease onset, number of infectious episodes, and each infection complication) were not linked with DHR activity or p67phox expression levels. In summary, our flow cytometric analysis revealed a significant correlation between residual ROS production and p67phox expression. More deleterious NCF2 genotypes were associated with lower levels of DHR activity and p67phox expression. DHR assays and protein expression analysis by using flow cytometry may be relevant strategies for predicting the genotypes of p67phox-CGD.

Identification of CD14 and lipopolysaccharide-binding protein as novel biomarkers for sarcoidosis using proteomics of serum extracellular vesicles.

Sarcoidosis is a complex, polygenic, inflammatory granulomatous multi-organ disease of unknown cause. The granulomatous inflammation in sarcoidosis is driven by the interplay between T cells and macrophages. Extracellular vesicles (EVs) play important roles in intercellular communication. We subjected serum EVs, isolated by size exclusion chromatography, from seven patients with sarcoidosis and five control subjects to non-targeted proteomics analysis. Non-targeted, label-free proteomics analysis detected 2292 proteins in serum EVs; 42 proteins were up-regulated in patients with sarcoidosis relative to control subjects; and 324 proteins were down-regulated. The protein signature of EVs from patients with sarcoidosis reflected disease characteristics such as antigen presentation and immunological disease. Candidate biomarkers were further verified by targeted proteomics analysis (selected reaction monitoring) in 46 patients and 10 control subjects. Notably, CD14 and lipopolysaccharide-binding protein (LBP) were validated by targeted proteomics analysis. Up-regulation of these proteins was further confirmed by immunoblotting, and their expression was strongly increased in macrophages of lung granulomatous lesions. Consistent with these findings, CD14 levels were increased in lipopolysaccharide-stimulated macrophages during multinucleation, concomitant with increased levels of CD14 and LBP in EVs. The area under the curve values of CD14 and LBP were 0.81 and 0.84, respectively, and further increased to 0.98 in combination with angiotensin-converting enzyme and soluble interleukin-2 receptor. These findings suggest that CD14 and LBP in serum EVs, which are associated with granulomatous pathogenesis, can improve the diagnostic accuracy in patients with sarcoidosis.

Diagnostic significance of secondary bacteremia in patients with COVID-19.

OBJECTIVES: We investigated the occurrence of non-respiratory bacterial and fungal secondary infections, causative organisms, impact on clinical outcomes, and association between the secondary pathogens and mortality in hospitalized patients with coronavirus disease 2019 (COVID-19). METHODS: This was a retrospective cohort study that included data from inpatients with COVID-19 from multiple centers participating in the Japan COVID-19 Taskforce (April 2020 to May 2021). We obtained demographic, epidemiological, and microbiological data throughout the course of hospitalization and analyzed the cases of COVID-19 complicated by non-respiratory bacterial infections. RESULTS: Of the 1914 patients included, non-respiratory bacterial infections with COVID-19 were diagnosed in 81 patients (4.2%). Of these, 59 (3.1%) were secondary infections. Bacteremia was the most frequent bacterial infection, occurring in 33 cases (55.9%), followed by urinary tract infections in 16 cases (27.1%). Staphylococcus epidermidis was the most common causative organism of bacteremia. Patients with COVID-19 with non-respiratory secondary bacterial infections had significantly higher mortality, and a multivariate logistic regression analysis demonstrated that those with bacteremia (aOdds Ratio = 15.3 [5.97-39.1]) were at higher risk of death. Multivariate logistic regression analysis showed that age, male sex, use of steroids to treat COVID-19, and intensive care unit admission increased the risk for nosocomial bacteremia. CONCLUSIONS: Secondary bacteremia is an important complication that may lead to poor prognosis in cases with COVID-19. An appropriate medical management strategy must be established, especially for patients with concomitant predisposing factors.

Impact of respiratory bacterial infections on mortality in Japanese patients with COVID-19: a retrospective cohort study.

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.

Multi-trait and cross-population genome-wide association studies across autoimmune and allergic diseases identify shared and distinct genetic component.

OBJECTIVES: Autoimmune and allergic diseases are outcomes of the dysregulation of the immune system. Our study aimed to elucidate differences or shared components in genetic backgrounds between autoimmune and allergic diseases. METHODS: We estimated genetic correlation and performed multi-trait and cross-population genome-wide association study (GWAS) meta-analysis of six immune-related diseases: rheumatoid arthritis, Graves' disease, type 1 diabetes for autoimmune diseases and asthma, atopic dermatitis and pollinosis for allergic diseases. By integrating large-scale biobank resources (Biobank Japan and UK biobank), our study included 105 721 cases and 433 663 controls. Newly identified variants were evaluated in 21 778 cases and 712 767 controls for two additional autoimmune diseases: psoriasis and systemic lupus erythematosus. We performed enrichment analyses of cell types and biological pathways to highlight shared and distinct perspectives. RESULTS: Autoimmune and allergic diseases were not only mutually classified based on genetic backgrounds but also they had multiple positive genetic correlations beyond the classifications. Multi-trait GWAS meta-analysis newly identified six allergic disease-associated loci. We identified four loci shared between the six autoimmune and allergic diseases (rs10803431 at PRDM2, OR=1.07, p=2.3×10-8, rs2053062 at G3BP1, OR=0.90, p=2.9×10-8, rs2210366 at HBS1L, OR=1.07, p=2.5×10-8 in Japanese and rs4529910 at POU2AF1, OR=0.96, p=1.9×10-10 across ancestries). Associations of rs10803431 and rs4529910 were confirmed at the two additional autoimmune diseases. Enrichment analysis demonstrated link to T cells, natural killer cells and various cytokine signals, including innate immune pathways. CONCLUSION: Our multi-trait and cross-population study should elucidate complex pathogenesis shared components across autoimmune and allergic diseases.

Association of the RPA3-UMAD1 locus with interstitial lung diseases complicated with rheumatoid arthritis in Japanese.

OBJECTIVES: The genetic background of rheumatoid arthritis-interstitial lung disease (RA-ILD) has been evaluated in Europeans, but little knowledge has been obtained in non-Europeans. This study aimed to elucidate genome-wide risk of RA-ILD in non-Europeans. METHODS: We performed an initial genome-wide association study (GWAS) of RA-ILD in the Japanese population. By conducting the meta-analysis of the three GWAS datasets of the RA cohorts and biobank of Japanese, our study included 358 RA-ILD cases and 4550 RA subjects without ILD. We then conducted the stratified analysis of the effect of the GWAS risk allele in each CT image pattern. RESULTS: We identified one novel RA-ILD risk locus at 7p21 that satisfied the genome-wide significance threshold (rs12702634 at RPA3-UMAD1, OR=2.04, 95% CI 1.59 to 2.60, p=1.5×10-8). Subsequent stratified analysis based on the CT image patterns demonstrated that the effect size of the RA-ILD risk allele (rs12702634-C) was large with the UIP pattern (OR=1.86, 95% CI 0.97 to 3.58, p=0.062) and the probable UIP pattern (OR=2.26, 95% CI 1.36 to 3.73, p=0.0015). CONCLUSION: We revealed one novel genetic association with RA-ILD in Japanese. The RA-ILD risk of the identified variant at RPA3-UMAD1 was relatively high in the CT image patterns related to fibrosis. Our study should contribute to elucidation of the complicated aetiology of RA-ILD.

Subtype-specific gout susceptibility loci and enrichment of selection pressure on ABCG2 and ALDH2 identified by subtype genome-wide meta-analyses of clinically defined gout patients.

OBJECTIVES: Genome-wide meta-analyses of clinically defined gout were performed to identify subtype-specific susceptibility loci. Evaluation using selection pressure analysis with these loci was also conducted to investigate genetic risks characteristic of the Japanese population over the last 2000-3000 years. METHODS: Two genome-wide association studies (GWASs) of 3053 clinically defined gout cases and 4554 controls from Japanese males were performed using the Japonica Array and Illumina Array platforms. About 7.2 million single-nucleotide polymorphisms were meta-analysed after imputation. Patients were then divided into four clinical subtypes (the renal underexcretion type, renal overload type, combined type and normal type), and meta-analyses were conducted in the same manner. Selection pressure analyses using singleton density score were also performed on each subtype. RESULTS: In addition to the eight loci we reported previously, two novel loci, PIBF1 and ACSM2B, were identified at a genome-wide significance level (p<5.0×10-8) from a GWAS meta-analysis of all gout patients, and other two novel intergenic loci, CD2-PTGFRN and SLC28A3-NTRK2, from normal type gout patients. Subtype-dependent patterns of Manhattan plots were observed with subtype GWASs of gout patients, indicating that these subtype-specific loci suggest differences in pathophysiology along patients' gout subtypes. Selection pressure analysis revealed significant enrichment of selection pressure on ABCG2 in addition to ALDH2 loci for all subtypes except for normal type gout. CONCLUSIONS: Our findings on subtype GWAS meta-analyses and selection pressure analysis of gout will assist elucidation of the subtype-dependent molecular targets and evolutionary involvement among genotype, phenotype and subtype-specific tailor-made medicine/prevention of gout and hyperuricaemia.

Amphiphilic triphenylamine-benzothiadiazole dyes: preparation, fluorescence and aggregation behavior, and enzyme fluorescence detection.

Fluorescence change systems that can respond to biological objects have attracted attention for use as biological probes and sensors. In this study, we report emission enhancement in a fluorescent aggregate composed of amphiphilic donor-acceptor dye molecules. The emission efficiency of the aggregate was reduced upon introducing a hydrophilic galactopyranose moiety, because of the decrease in the aggregate stability, which in turn was due to disruption of the hydrophilicity-hydrophobicity balance. In contrast, emission enhancement could be achieved by treatment with ß-galactosidase, as a result of the removal of the galactopyranose moiety. The change in aggregate stabilization based on the hydrophilicity-hydrophobicity balance leads to the emission enhancement into detectable ß-galactosidase activity.

Large-scale cross-trait genetic analysis highlights shared genetic backgrounds of autoimmune diseases.

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.

CD98 heavy chain protein is overexpressed in non-small cell lung cancer and is a potential target for CAR T-cell therapy.

Chimeric antigen receptor (CAR) T cells are effective against hematological cancers, but are less effective against solid tumors such as non-small cell lung cancer (NSCLC). One of the reasons is that only a few cell surface targets specific for NSCLC cells have been identified. Here, we report that CD98 heavy chain (hc) protein is overexpressed on the surface of NSCLC cells and is a potential target for CAR T cells against NSCLC. Screening of over 10,000 mAb clones raised against NSCLC cell lines showed that mAb H2A011 bound to NSCLC cells but not normal lung epithelial cells. H2A011 recognized CD98hc. Although CAR T cells derived from H2A011 could not be established presumably due to the high level of H2A011 reactivity in activated T cells, those derived from the anti-CD98hc mAb R8H283, which had been shown to lack reactivity with CD98hc glycoforms expressed on normal hematopoietic cells and some normal tissues, were successfully developed. R8H283 specifically reacted with NSCLC cells in six of 15 patients. R8H283-derived CAR T cells exerted significant anti-tumor effects in a xenograft NSCLC model in vivo. These results suggest that R8H283 CAR T cells may become a new therapeutic tool for NSCLC, although careful testing for off-tumor reactivity should be performed in the future.

SFTPB in serum extracellular vesicles as a biomarker of progressive pulmonary fibrosis.

Progressive pulmonary fibrosis (PPF), defined as the worsening of various interstitial lung diseases (ILDs), currently lacks useful biomarkers. To identify novel biomarkers for early detection of patients at risk of PPF, we performed a proteomic analysis of serum extracellular vesicles (EVs). Notably, the identified candidate biomarkers were enriched for lung-derived proteins participating in fibrosis-related pathways. Among them, pulmonary surfactant-associated protein B (SFTPB) in serum EVs could predict ILD progression better than the known biomarkers, serum KL-6 and SP-D, and it was identified as an independent prognostic factor from ILD-gender-age-physiology index. Subsequently, the utility of SFTPB for predicting ILD progression was evaluated further in 2 cohorts using serum EVs and serum, respectively, suggesting that SFTPB in serum EVs but not in serum was helpful. Among SFTPB forms, pro-SFTPB levels were increased in both serum EVs and lungs of patients with PPF compared with those of the control. Consistently, in a mouse model, the levels of pro-SFTPB, primarily originating from alveolar epithelial type 2 cells, were increased similarly in serum EVs and lungs, reflecting pro-fibrotic changes in the lungs, as supported by single-cell RNA sequencing. SFTPB, especially its pro-form, in serum EVs could serve as a biomarker for predicting ILD progression.

Quantification of escape from X chromosome inactivation with single-cell omics data reveals heterogeneity across cell types and tissues.

Several X-linked genes escape from X chromosome inactivation (XCI), while differences in escape across cell types and tissues are still poorly characterized. Here, we developed scLinaX for directly quantifying relative gene expression from the inactivated X chromosome with droplet-based single-cell RNA sequencing (scRNA-seq) data. The scLinaX and differentially expressed gene analyses with large-scale blood scRNA-seq datasets consistently identified the stronger escape in lymphocytes than in myeloid cells. An extension of scLinaX to a 10x multiome dataset (scLinaX-multi) suggested a stronger escape in lymphocytes than in myeloid cells at the chromatin-accessibility level. The scLinaX analysis of human multiple-organ scRNA-seq datasets also identified the relatively strong degree of escape from XCI in lymphoid tissues and lymphocytes. Finally, effect size comparisons of genome-wide association studies between sexes suggested the underlying impact of escape on the genotype-phenotype association. Overall, scLinaX and the quantified escape catalog identified the heterogeneity of escape across cell types and tissues.

Evaluation of Muscle Oxygen Dynamics in Children's Gait and Its Relationship with the Physiological Cost Index.

The response of muscle oxygen saturation, which is an index for the energy metabolism of muscles during walking in children, and its relationship to the physiological cost index, which indicates walking efficiency, are unknown. This study aimed to evaluate muscle oxygen saturation in lower extremity muscles during walking in children, its changes with age, and the relationship between the physiological cost index. The oxygen saturation was measured by the amount of change during a two-minute walk, and the physiological cost index was calculated from the change in heart rate before and after exercise and walking speed. Results were compared for each muscle, and the correlation between the two was examined. Changes in muscle oxygen saturation were greater in the lower leg muscles, significantly greater in the tibialis anterior at six to seven years, and in the gastrocnemius medial head at eight to ten years. The physiological cost index was significantly correlated with changes in muscle oxygen saturation in the tibialis anterior (r = 0.44, p < 0.001). The lower leg muscles were metabolically active in children's gait, and their response varied with age. Moreover, the muscle oxygenation dynamics of the tibialis anterior may influence walking efficiency.

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.

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.