Body positioning-related laryngeal narrowing pattern and expiratory mechanical constraints in advanced chronic obstructive pulmonary disease: A case report.

We previously reported that laryngeal widening led to improved exercise tolerance in COPD. However, it is not clear whether laryngeal narrowing occurs as a compensatory response to tracheal movement or is affected by posture. Here, we report the case of an advanced COPD patient whose more prolonged expiration in a head-forward leaning position compared with that in a neck-extended position occurred with an excessive duration of severe laryngeal narrowing without tracheal obstruction, which led to exercise intolerance with expiratory mechanical constraints. This case provided useful insights into the regulation of the upper airway with body positioning for improving exercise tolerance.

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.

Diagnosis and Management of Drug-Induced Interstitial Lung Disease Associated with Amikacin Liposome Inhalation Suspension in Refractory Mycobacterium Avium Complex Pulmonary Disease: A Case Report.

Amikacin liposome inhalation suspension (ALIS) is a key drug for the treatment of refractory Mycobacterium avium complex pulmonary disease (MAC-PD). Although cases of drug-induced interstitial lung disease (DIILD) by ALIS have been reported, its diagnosis is challenging due to overlapping existing pulmonary shadows, airway bleeding, exacerbation of underlying conditions, and the potential for various concurrent infections. A 72-year-old woman started treatment with ALIS for refractory MAC-PD. Three weeks later, she had a fever, cough, and appetite loss. She was hospitalized because multiple infiltrative opacities were observed on chest X-ray and chest computed tomography. Because the opacities worsened after empiric antibiotic therapy with broad-spectrum antibiotics, we initiated corticosteroid therapy, suspecting DIILD caused by ALIS, although drug lymphocyte stimulation tests for ALIS and amikacin were negative. Three days later, we found signs of improvement and quickly tapered the corticosteroids. After obtaining informed consent, we performed a drug provocation test of ALIS. Seven days later, she exhibited fever, an increased peripheral white blood cell count, and elevated serum C-reactive protein level, all of which returned to baseline 4 days after stopping ALIS, leading to a diagnosis of DIILD caused by ALIS in this patient. DIILD caused by ALIS is rare but should be carefully diagnosed to ensure that patients with refractory MAC-PD do not miss the opportunity to receive ALIS treatment.

Laryngeal widening and adequate ventilation by expiratory pressure load training improve aerobic capacity in COPD: a randomised controlled trial.

RATIONALE: Despite strategies acting on peripheral airway obstruction in chronic obstructive pulmonary disease (COPD), exercise intolerance remains inadequately improved. We hypothesised that laryngeal narrowing is a potential treatment target of expiratory pressure load training (EPT) to improve exercise intolerance in COPD. METHODS: The effect of 3-month EPT was assessed in 47 patients with COPD divided into Global Initiative for Chronic Obstructive Lung Disease (GOLD) mild-to-moderate (I-II) and severe-to-very severe (III-IV), randomly allocating 1:1 to EPT or control groups. The primary outcome was endurance time in the constant work rate exercise test in GOLD III-IV patients. RESULTS: Compared with controls, EPT increased: (1) endurance time, with estimated treatment effect: +703 (95% CI: 379 to 1031) s, p=0.0008 (GOLD I-II); +390 (95% CI: 205 to 574) s, p=0.0006 (GOLD III-IV); (2) peak oxygen uptake (p=0.0086 in GOLD I-II; p=0.0004 in GOLD III-IV); (3) glottic dilatation ratio at maximum collapse on laryngoscopy in the submaximal exercise (p=0.0062 in GOLD I-II; p=0.0001 in GOLD III-IV); and (4) the inflection point of expiratory tidal volume relative to minute ventilation during the incremental exercise (p=0.0015 in GOLD I-II; p=0.0075 in GOLD III-IV). Across GOLD grades, the responses of glottic dilatation ratio at maximum collapse and the expiratory tidal volume at the inflection point were selected as more influential variables correlating with the improvement in peak oxygen uptake and endurance time, respectively. CONCLUSION: These results show that EPT improved aerobic capacity and endurance time with larger laryngeal widening and adequate ventilation despite advanced COPD. TRIAL REGISTRATION NUMBER: UMIN000041250.

A randomized phase 2 study on demeclocycline in patients with mild-to-moderate COVID-19.

Tetracyclines exhibit anti-viral, anti-inflammatory, and immunomodulatory activities via various mechanisms. The present study investigated the efficacy and safety of demeclocycline in patients hospitalized with mild-to-moderate COVID-19 via an open-label, multicenter, parallel-group, randomized controlled phase 2 trial. Primary and secondary outcomes included changes from baseline (day 1, before the study treatment) in lymphocytes, cytokines, and SARS-CoV-2 RNA on day 8. Seven, seven, and six patients in the control, demeclocycline 150 mg daily, and demeclocycline 300 mg daily groups, respectively, were included in the modified intention-to-treat population that was followed until day 29. A significant change of 191.3/µL in the number of CD4+ T cells from day 1 to day 8 was observed in the demeclocycline 150 mg group (95% CI 5.1/µL-377.6/µL) (p = 0.023), whereas that in the control group was 47.8/µL (95% CI - 151.2/µL to 246.8/µL), which was not significant (p = 0.271). The change rates of CD4+ T cells negatively correlated with those of IL-6 in the demeclocycline-treated groups (R = - 0.807, p = 0.009). All treatment-emergent adverse events were of mild-to-moderate severity. The present results indicate that the treatment of mild-to-moderate COVID-19 patients with demeclocycline elicits immune responses conducive to recovery from COVID-19 with good tolerability.Trial registration: This study was registered with the Japan Registry of Clinical Trials (Trial registration number: jRCTs051200049; Date of the first registration: 26/08/2020).

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.

MGIT-seq for the Identification of Nontuberculous Mycobacteria and Drug Resistance: a Prospective Study.

Because nontuberculous mycobacterial pulmonary disease is a considerable health burden, a simple and clinically applicable analytical protocol enabling the identification of subspecies and drug-resistant disease is required to determine the treatment strategy. We aimed to develop a simplified workflow consisting only of direct sequencing of mycobacterial growth indicator tube cultures (MGIT-seq). In total, 138 patients were prospectively enrolled between April 2021 and May 2022, and culture-positive MGIT broths were subjected to sequencing using MinION, a portable next-generation sequencer. Sequence analysis was conducted to identify species using core genome multilocus sequence typing and to predict macrolide and amikacin (AMK) resistance based on previously reported mutations in rrl, rrs, and erm(41). The results were compared to clinical tests for species identification and drug susceptibility. A total of 116 patients with positive MGIT cultures were included in the analysis. MGIT-seq yielded 99.1% accuracy in species-level identification and identified 98 isolates (84.5%) at the subspecies level. Macrolide and AMK resistance were detected in 19.4% and 1.9% of Mycobacterium avium complex (MAC) and Mycobacterium abscessus isolates. The predicted macrolide and AMK resistance was consistent with the results of conventional drug susceptibility tests, with specificities of 97.6% and 100.0%, respectively. Direct MGIT-seq has achieved comprehensive identification and drug resistance detection of nontuberculous mycobacteria, which could be applicable to determine the treatment strategy by a single test in clinical practice.

Real-world impact of antifibrotics on prognosis in patients with progressive fibrosing interstitial lung disease.

OBJECTIVE: No studies have demonstrated the real-world efficacy of antifibrotics for progressive fibrosing interstitial lung disease (PF-ILD). Therefore, we evaluated the efficacy of antifibrotics in patients with PF-ILD. METHODS: We retrospectively reviewed the medical records of patients with ILD from January 2012 to July 2021. Patients were diagnosed with PF-ILD if they had ≥10% fibrosis on high-resolution CT (HRCT) and a relative forced vital capacity (FVC) decline of either ≥10% or >5% to <10% with clinical deterioration or progression of fibrosis on HRCT during overlapping windows of 2 years and with a %FVC of ≥45%. We compared FVC changes and overall survival (OS) between patients with and without antifibrotics. FVC changes were analysed using generalised estimating equations. We used inverse probability weighting (IPW) and statistical matching to adjust for covariates. RESULTS: Of the 574 patients, 167 were diagnosed with PF-ILD (idiopathic pulmonary fibrosis (IPF), n=64; non-IPF, n=103). Antifibrotics improved the FVC decline in both IPF (p=0.002) and non-IPF (p=0.05) (IPW: IPF, p=0.015; non-IPF, p=0.031). Among patients with IPF, OS was longer in the antifibrotic group (log-rank p=0.001). However, among patients with non-IPF, OS was not longer in the antifibrotic group (p=0.3263) (IPW and statistical matching: IPF, p=0.0534 and p=0.0018; non-IPF, p=0.5663 and p=0.5618). CONCLUSION: This is the first real-world study to show that antifibrotics improve the FVC decline in PF-ILD. However, among patients with non-IPF, we found no significant difference in mortality between those with and without antifibrotics. Future studies must clarify whether antifibrotics improve the prognosis of non-IPF.

Genomic repertoires linked with pathogenic potency of arthritogenic Prevotella copri isolated from the gut of patients with rheumatoid arthritis.

OBJECTIVES: Prevotella copri is considered to be a contributing factor in rheumatoid arthritis (RA). However, in some non-Westernised countries, healthy individuals also harbour an abundance of P. copri in the intestine. This study investigated the pathogenicity of RA patient-derived P. copri (P. copri RA) compared with healthy control-derived P. copri (P. copri HC). METHODS: We obtained 13 P. copri strains from the faeces of patients with RA and healthy controls. Following whole genome sequencing, the sequences of P. copri RA and P. copri HC were compared. To analyse the arthritis-inducing ability of P. copri, we examined two arthritis models (1) a collagen-induced arthritis model harbouring P. copri under specific-pathogen-free conditions and (2) an SKG mouse arthritis model under P. copri-monocolonised conditions. Finally, to evaluate the ability of P. copri to activate innate immune cells, we performed in vitro stimulation of bone marrow-derived dendritic cells (BMDCs) by P. copri RA and P. copri HC. RESULTS: Comparative genomic analysis revealed no apparent differences in the core gene contents between P. copri RA and P. copri HC, but pangenome analysis revealed the high genome plasticity of P. copri. We identified a P. copri RA-specific genomic region as a conjugative transposon. In both arthritis models, P. copri RA-induced more severe arthritis than P. copri HC. In vitro BMDC stimulation experiments revealed the upregulation of IL-17 and Th17-related cytokines (IL-6, IL-23) by P. copri RA. CONCLUSION: Our findings reveal the genetic diversity of P. copri, and the genomic signatures associated with strong arthritis-inducing ability of P. copri RA. Our study contributes towards elucidation of the complex pathogenesis of RA.

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.

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.

Bronchoalveolar lavage fluid reveals factors contributing to the efficacy of PD-1 blockade in lung cancer.

Bronchoalveolar lavage is commonly performed to assess inflammation and identify responsible pathogens in lung diseases. Findings from bronchoalveolar lavage might be used to evaluate the immune profile of the lung tumor microenvironment (TME). To investigate whether bronchoalveolar lavage fluid (BALF) analysis can help identify patients with non-small cell lung cancer (NSCLC) who respond to immune checkpoint inhibitors (ICIs), BALF and blood were prospectively collected before initiating nivolumab. The secreted molecules, microbiome, and cellular profiles based on BALF and blood analysis of 12 patients were compared with regard to therapeutic effect. Compared with ICI nonresponders, responders showed significantly higher CXCL9 levels and a greater diversity of the lung microbiome profile in BALF, along with a greater frequency of the CD56+ subset in blood T cells, whereas no significant difference in PD-L1 expression was found in tumor cells. Antibiotic treatment in a preclinical lung cancer model significantly decreased CXCL9 in the lung TME, resulting in reduced sensitivity to anti-PD-1 antibody, which was reversed by CXCL9 induction in tumor cells. Thus, CXCL9 might be associated with the lung TME microbiome, and the balance of CXCL9 and lung TME microbiome could contribute to nivolumab sensitivity in patients with NSCLC. BALF analysis can help predict the efficacy of ICIs when performed along with currently approved examinations.

Increased Oxygen Extraction by Pulmonary Rehabilitation Improves Exercise Tolerance and Ventilatory Efficiency in Advanced Chronic Obstructive Pulmonary Disease.

BACKGROUND: In cardiopulmonary exercise testing (CPET), oxygen uptake (V'O2) is calculated using the product of minute ventilation (V'E) and the difference between inspiratory and expiratory O2 concentrations (ΔFO2). However, little is known about the response of ΔFO2 to pulmonary rehabilitation (PR). The aim of the present study was (1) to investigate whether PR increases peak V'O2, based on whether ΔFO2 or V'E at peak exercise increase after PR, and (2) to investigate whether an improvement in ΔFO2 correlates with an improvement in ventilatory efficiency. METHODS: A total of 38 patients with severe and very severe COPD, whose PR responses were evaluated by CPET, were retrospectively analyzed. RESULTS: After PR, peak V'O2 was increased in 14 patients. The difference in ΔFO2 at peak exercise following PR correlated with the difference in peak V'O2 (r = 0.4884, p = 0.0019), the difference in V'E/V'CO2-nadir (r = -0.7057, p < 0.0001), and the difference in V'E-V'CO2 slope (r = -0.4578, p = 0.0039), but it did not correlate with the difference in peak V'E. CONCLUSIONS: The increased O2 extraction following PR correlated with improved exercise tolerance and ventilatory efficiency. In advanced COPD patients, a new strategy for improving O2 extraction ability might be effective in those in whom ventilatory ability can be only minimally increased.

The Diagnosis of Nontuberculous Mycobacterial Pulmonary Disease by Single Bacterial Isolation Plus Anti-GPL-Core IgA Antibody.

Although serum anti-glycopeptidolipid (GPL)-core IgA antibody is a highly specific test for infection with Mycobacterium avium complex (MAC), Mycobacterium abscessus, and its subspecies abscessus, subsp. massiliense, and subsp. bolletii (MAB), its use for the definitive diagnosis of MAC pulmonary disease (PD) and MAB-PD are unknown. To clarify the diagnostic accuracy of the anti-GPL-core IgA antibody test among patients with radiologically suspected MAC-PD or MAB-PD who already have a single positive sputum culture test. The first isolations of MAC and MAB from patients with radiologically suspected MAC-PD or MAB-PD at the Osaka Toneyama Medical Center between January 2006 and December 2020 were collected. Patients were enrolled when their serum anti-GPL-core IgA antibody was measured during the 3 months before and after the first isolation. We retrospectively compared the results of anti-GPL-core IgA antibody testing with the final diagnoses based on the current guidelines. We included 976 patients for analysis. The serum anti-GPL-core IgA antibody was positive in 699 patients (71.6%). The positive predictive value of anti-GPL-core IgA antibody for the diagnosis of MAC-PD or MAB-PD was 97.4%. The median time required for the second positive culture after the first isolation was 51 days (interquartile range 12 to 196 days). The positive serum anti-GPL-core IgA antibody test allowed an early and definitive diagnosis of MAC-PD or MAB-PD in those who already had a single positive sputum culture test. IMPORTANCE To satisfy the microbiologic criteria of the current diagnostic guideline for nontuberculous mycobacterial pulmonary disease (PD), at least two positive sputum cultures of the same species of mycobacteria from sputum are required to avoid the casual isolation of mycobacteria. This study showed that the positivity of a serum anti-glycopeptidolipid (GPL)-core IgA antibody test has an excellent diagnostic ability among patients with radiologically suspected Mycobacterium avium complex (MAC)-PD or Mycobacterium abscessus (MAB)-PD who already had a single positive sputum culture test. The usage of single culture isolation plus anti-GPL-core IgA antibody as another diagnostic criterion has a time, cost, and effort-saving effect. Furthermore, it will facilitate the diagnosis of MAC-PD or MAB-PD in the early stage of disease because serum anti-GPL-core IgA antibody becomes high in these patients. Therefore, we proposed adding single culture isolation plus anti-GPL-core IgA antibody as "combined microbiological and serological criteria" to the diagnostic guidelines for MAC-PD and MAB-PD.

Whole gut virome analysis of 476 Japanese revealed a link between phage and autoimmune disease.

OBJECTIVE: The relationship between autoimmune diseases and the gut microbiome has been intensively studied, and several autoimmunity-associated bacterial taxa have been identified. However, much less is known about the roles of the gut virome in autoimmune diseases. METHODS: Here, we performed a whole gut virome analysis based on the shotgun sequencing of 476 Japanese which included patients with rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), multiple sclerosis and healthy control subjects. RESULTS: Our case-control comparison of the viral abundance revealed that crAss-like phages, which are one of the main components of a healthy gut virome, significantly decreased in the gut of the patients with autoimmune disease, specifically the patients with RA and SLE. In addition, Podoviridae significantly decreased in the gut of the patients with SLE. To understand how these viruses affected the bacteriome, we performed a quantitative virus-bacterium association analysis and clustered regularly interspaced short palindromic repeat-based virus-bacterium interaction analysis. We identified a symbiosis between Podoviridae and Faecalibacterium. In addition, multiple bacterial targets of crAss-like phages were identified (eg, Ruminococcus spp). CONCLUSION: Our data suggest that the gut virome can affect our body either directly or via bacteria. Our analyses have elucidated a previously missing part of the autoimmunity-associated gut microbiome and presented new candidates that contribute to the development of autoimmune diseases.

Genetic architecture of microRNA expression and its link to complex diseases in the Japanese population.

Understanding the genetic effects on non-coding RNA (ncRNA) expression facilitates functional characterization of disease-associated genetic loci. Among several classes of ncRNAs, microRNAs (miRNAs) are key post-transcriptional gene regulators. Despite its biological importance, previous studies on the genetic architecture of miRNA expression focused mostly on the European individuals, underrepresented in other populations. Here, we mapped miRNA expression quantitative trait loci (miRNA-eQTL) for 343 miRNAs in 141 Japanese using small RNA sequencing and whole-genome sequencing, identifying 1275 cis-miRNA-eQTL variants for 40 miRNAs (false discovery rate < 0.2). Of these, 25 miRNAs having eQTL were unreported in the European studies, including 5 miRNAs with their lead variant monomorphic in the European populations, which demonstrates the value of miRNA-eQTL analysis in diverse ancestral populations. MiRNAs with eQTL effect showed allele-specific expression (ASE; e.g. miR-146a-3p), and ASE analysis further detected cis-regulatory variants not captured by the conventional miRNA-eQTL mapping (e.g. miR-933). We identified a copy number variation associated with miRNA expression (e.g. miR-570-3p, P = 7.2 × 10-6), which contributes to a more comprehensive landscape of miRNA-eQTLs. To elucidate a post-transcriptional modification in miRNAs, we created a catalog of miRNA-editing sites, including 10 canonical and 6 non-canonical sites. Finally, by integrating the miRNA-eQTLs and Japanese genome-wide association studies of 25 complex traits (mean n = 192 833), we conducted a transcriptome-wide association study, identifying miR-1908-5p as a potential mediator for adult height, colorectal cancer and type 2 diabetes (P < 9.1 × 10-5). Our study broadens the population diversity in ncRNA-eQTL studies and contributes to functional annotation of disease-associated loci found in non-European populations.

Prokaryotic and viral genomes recovered from 787 Japanese gut metagenomes revealed microbial features linked to diets, populations, and diseases.

We reconstructed 19,084 prokaryotic and 31,395 viral genomes from 787 Japanese gut metagenomes as Japanese metagenome-assembled genomes (JMAG) and Japanese Virus Database (JVD), which are large microbial genome datasets for a single population. Population-specific enrichment of the Bacillus subtilis and ß-porphyranase among the JMAG could derive from the Japanese traditional food natto (fermented soybeans) and nori (laver), respectively. Dairy-related Enterococcus_B lactis and Streptococcus thermophilus were nominally associated with the East Asian-specific missense variant rs671:G>A in ALDH2, which was associated with dairy consumption. Of the species-level viral genome clusters in the JVD, 62.9% were novel. The ß crAss-like phage composition was low among the Japanese but relatively high among African and Oceanian peoples. Evaluations of the association between crAss-like phages and diseases showed significant disease-specific associations. Our large catalog of virus-host pairs identified the positive correlation between the abundance of the viruses and their hosts.