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.

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.

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.

Identification of serum metabolome signatures associated with retinal and renal complications of type 2 diabetes.

BACKGROUND: Type 2 diabetes is a common disease around the world and its major complications are diabetic retinopathy (DR) and diabetic kidney disease (DKD). Persons with type 2 diabetes with complications, especially who have both DR and DKD, have poorer prognoses than those without complications. Therefore, prevention and early identification of the complications of type 2 diabetes are necessary to improve the prognosis of persons with type 2 diabetes. The aim of this study is to identify factors associated with the development of multiple complications of type 2 diabetes. METHODS: We profiled serum metabolites of persons with type 2 diabetes with both DR and DKD (N = 141) and without complications (N = 159) using a comprehensive non-targeted metabolomics approach with mass spectrometry. Based on the serum metabolite profiles, case-control comparisons and metabolite set enrichment analysis (MSEA) were performed. RESULTS: Here we show that five metabolites (cyclohexylamine, P = 4.5 × 10-6; 1,2-distearoyl-glycero-3-phosphocholine, P = 7.3 × 10-6; piperidine, P = 4.8 × 10-4; N-acetylneuraminic acid, P = 5.1 × 10-4; stearoyl ethanolamide, P = 6.8 × 10-4) are significantly increased in those with the complications. MSEA identifies fatty acid biosynthesis as the type 2 diabetes complications-associated biological pathway (P = 0.0020). CONCLUSIONS: Our metabolome analysis identifies the serum metabolite features of the persons with type 2 diabetes with multiple complications, which could potentially be used as biomarkers.

In the management of type 2 diabetes, prevention and early identification of diabetes complications are important. In particular, people with type 2 diabetes with diabetic retinopathy (DR), affecting the eye, and diabetic kidney disease (DKD), have poorer outcomes than those without complications and need early intervention. Here, we comprehensively profiled blood metabolites, or breakdown products of the biological processes occurring in the body, of people with type 2 diabetes with both DR and DKD and those without complications. We found that five metabolites were significantly increased in those with complications, and we identified a specific metabolic pathway associated with having complications. Our analysis identified the blood metabolite features of people with type 2 diabetes with multiple complications, which could potentially be used as markers in the future.

SLE stratification based on BAFF and IFN-I bioactivity for biologics and implications of BAFF produced by glomeruli in lupus nephritis.

OBJECTIVE: B-cell activating factor (BAFF) is implicated in SLE pathogenesis. Blocking BAFF signalling has contributed to reducing glucocorticoid dosage and preventing organ damage. However, clinical characteristics of patients who may benefit from this therapy are not yet fully elucidated. Therefore, we identified patients with high BAFF-bioactivity to investigate their clinical characteristics and BAFF-producing cells. METHODS: We established the reporter cell for BAFF and investigated the clinical characteristics of SLE patients with high BAFF-bioactivity. We identified BAFF-expressing kidney cells using publicly available scRNA-seq data and immunohistological analysis. SLE patients were stratified based on the bioactivity of BAFF and type-I IFN (IFN-I) to identify associated characteristic clinical manifestations. RESULTS: SLE patients, especially patients with LN, had significantly higher serum BAFF-bioactivity than healthy controls (HC) and non-LN patients. Additionally, single-cell-RNA-seq data and immunohistological analysis of kidney samples from LN patients revealed that BAFF is expressed in glomerular macrophages and mesangial cells. Notably, BAFF bioactivity was elevated in the urine of LN patients compared with that of non-LN patients, while no IFN-I bioactivity was detected in the urine. Furthermore, SLE stratification based on bioactivities of serum BAFF and IFN-I revealed the clinical characteristics of patients: high BAFF represented patients with LN and high IFN-I represented patients with blood and skin manifestations. CONCLUSIONS: Monitoring urinary BAFF-bioactivity may be valuable in diagnosing LN. Furthermore, stratification based on serum BAFF and IFN-I bioactivities may allow the identification of appropriate patients for biologics targeting BAFF and IFN-I.

ETS1 governs pathological tissue-remodeling programs in disease-associated fibroblasts.

Fibroblasts, the most abundant structural cells, exert homeostatic functions but also drive disease pathogenesis. Single-cell technologies have illuminated the shared characteristics of pathogenic fibroblasts in multiple diseases including autoimmune arthritis, cancer and inflammatory colitis. However, the molecular mechanisms underlying the disease-associated fibroblast phenotypes remain largely unclear. Here, we identify ETS1 as the key transcription factor governing the pathological tissue-remodeling programs in fibroblasts. In arthritis, ETS1 drives polarization toward tissue-destructive fibroblasts by orchestrating hitherto undescribed regulatory elements of the osteoclast differentiation factor receptor activator of nuclear factor-κB ligand (RANKL) as well as matrix metalloproteinases. Fibroblast-specific ETS1 deletion resulted in ameliorated bone and cartilage damage under arthritic conditions without affecting the inflammation level. Cross-tissue fibroblast single-cell data analyses and genetic loss-of-function experiments lent support to the notion that ETS1 defines the perturbation-specific fibroblasts shared among various disease settings. These findings provide a mechanistic basis for pathogenic fibroblast polarization and have important therapeutic implications.

OMARU: a robust and multifaceted pipeline for metagenome-wide association study.

Microbiome is an essential omics layer to elucidate disease pathophysiology. However, we face a challenge of low reproducibility in microbiome studies, partly due to a lack of standard analytical pipelines. Here, we developed OMARU (Omnibus metagenome-wide association study with robustness), a new end-to-end analysis workflow that covers a wide range of microbiome analysis from phylogenetic and functional profiling to case-control metagenome-wide association studies (MWAS). OMARU rigorously controls the statistical significance of the analysis results, including correction of hidden confounding factors and application of multiple testing comparisons. Furthermore, OMARU can evaluate pathway-level links between the metagenome and the germline genome-wide association study (i.e. MWAS-GWAS pathway interaction), as well as links between taxa and genes in the metagenome. OMARU is publicly available (https://github.com/toshi-kishikawa/OMARU), with a flexible workflow that can be customized by users. We applied OMARU to publicly available type 2 diabetes (T2D) and schizophrenia (SCZ) metagenomic data (n = 171 and 344, respectively), identifying disease biomarkers through comprehensive, multilateral, and unbiased case-control comparisons of metagenome (e.g. increased Streptococcus vestibularis in SCZ and disrupted diversity in T2D). OMARU improves accessibility and reproducibility in the microbiome research community. Robust and multifaceted results of OMARU reflect the dynamics of the microbiome authentically relevant to disease pathophysiology.

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.

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.

Metagenome-wide association study revealed disease-specific landscape of the gut microbiome of systemic lupus erythematosus in Japanese.

OBJECTIVE: Alteration of the gut microbiome has been linked to the pathogenesis of systemic lupus erythematosus (SLE). However, a comprehensive view of the gut microbiome in SLE and its interaction with the host remains to be revealed. This study aimed to reveal SLE-associated changes in the gut microbiome and its interaction with the host by a comprehensive metagenome-wide association study (MWAS) followed by integrative analysis. METHODS: We performed a MWAS of SLE based on shotgun sequencing of the gut microbial DNA from Japanese individuals (Ncase=47, Ncontrol=203). We integrated the result of the MWAS with the genome-wide association study (GWAS) data and plasma metabolite data. RESULTS: Via species level phylogenetic analysis, we identified and validated increases of Streptococcus intermedius and Streptococcus anginosus in the patients with SLE. Microbial gene analysis revealed increases of Streptococcus-derived genes including one involved in redox reaction. Additionally, microbial pathways related to sulfur metabolism and flagella assembly were altered in the patients with SLE. We identified an overlap in the enriched biological pathways between the metagenome and the germline genome by comparing the result of the MWAS and the GWAS of SLE (ie, MWAS-GWAS interaction). α-diversity and ß-diversity analyses provided evidence of dysbiosis in the metagenome of the patients with SLE. Microbiome-metabolome association analysis identified positive dosage correlation of acylcarnitine with Streptococcus intermedius, an SLE-associated taxon. CONCLUSION: Our MWAS followed by integrative analysis revealed SLE-associated changes in the gut microbiome and its interaction with the host, which contribute to our understanding of the relationship between the microbiome and SLE.

Successful treatment of refractory acute lupus haemophagocytic syndrome using rituximab: a case report.

Systemic lupus erythematosus (SLE)-associated haemophagocytic lymphohistiocytosis (HLH) is called acute lupus haemophagocytic syndrome (ALHS), which is relatively rare but life-threatening. We present the case of a 43-year-old woman diagnosed with SLE with panniculitis, pleuritis, and autoimmune hepatitis. She was treated with high-dose glucocorticoids. Although disease activity temporarily improved, she developed fever, elevation of liver enzymes, hyperferritinemia, severe inflammatory response, and thrombocytopenia a month after starting glucocorticoids. Bone marrow biopsy was performed and haemophagocytosis was observed. She was diagnosed with ALHS on day 49. Since she developed ALHS during administration of glucocorticoids, her ALHS was determined to be refractory to glucocorticoid monotherapy; therefore, additional immunosuppressive agents were needed. She was treated with methylprednisolone pulse, plasma exchange and cyclosporine A (CyA). However, CyA was discontinued on day 54 because CyA-induced hypertensive encephalopathy was suspected. Subsequently, rituximab (RTX) was introduced to treat refractory ALHS on day 56; the disease activity subsequently reduced. After four courses of RTX, her ferritin levels and platelet counts were within the normal range and the glucocorticoid dose could be tapered to betamethasone 2.0 mg/day on day 132. No subsequent recurrence of SLE and ALHS was observed until day 132. RTX might therefore be an effective therapeutic option for refractory ALHS.

Fibroblasts as a source of self-antigens for central immune tolerance.

Fibroblasts are one of the most common but also neglected types of stromal cells, the heterogeneity of which underlies the specific function of tissue microenvironments in development and regeneration. In the thymus, autoreactive T cells are thought to be negatively selected by reference to the self-antigens expressed in medullary epithelial cells, but the contribution of other stromal cells to tolerance induction has been poorly examined. In the present study, we report a PDGFR+ gp38+ DPP4- thymic fibroblast subset that is required for T cell tolerance induction. The deletion of the lymphotoxin ß-receptor in thymic fibroblasts caused an autoimmune phenotype with decreased expression of tissue-restricted and fibroblast-specific antigens, offering insight into the long-sought target of lymphotoxin signaling in the context of the regulation of autoimmunity. Thus, thymic medullary fibroblasts play an essential role in the establishment of central tolerance by producing a diverse array of self-antigens.

Chd4 choreographs self-antigen expression for central immune tolerance.

Autoreactive T cells are eliminated in the thymus to prevent autoimmunity by promiscuous expression of tissue-restricted self-antigens in medullary thymic epithelial cells. This expression is dependent on the transcription factor Fezf2, as well as the transcriptional regulator Aire, but the entire picture of the transcriptional program has been obscure. Here, we found that the chromatin remodeler Chd4, also called Mi-2ß, plays a key role in the self-antigen expression in medullary thymic epithelial cells. To maximize the diversity of self-antigen expression, Fezf2 and Aire utilized completely distinct transcriptional mechanisms, both of which were under the control of Chd4. Chd4 organized the promoter regions of Fezf2-dependent genes, while contributing to the Aire-mediated induction of self-antigens via super-enhancers. Mice deficient in Chd4 specifically in thymic epithelial cells exhibited autoimmune phenotypes, including T cell infiltration. Thus, Chd4 plays a critical role in integrating Fezf2- and Aire-mediated gene induction to establish central immune tolerance.

A case of complete atrioventricular block in secondary hemophagocytic syndrome/hemophagocytic lymphohistiocytosis recovered by plasma exchange and cytokine absorbing therapy with AN69ST continuous hemodiafiltration.

We report a case of incipient systemic lupus erythematosus (SLE) that rapidly progressed to complete atrioventricular block (cAVB). A 20-year-old man was admitted with facial erythema, painless oral aphtha, polyarthritis, and myalgia of each extremity. On admission, he developed first-degree atrioventricular block, pericarditis, pleuritis, renal failure, hemophagocytic lymphohistiocytosis, neurophychiatric SLE (left cerebellar infarction), and Staphylococcus aureus bacteremia. He was subsequently diagnosed with SLE based on several positive findings on immunological tests (including positive for antinuclear antibody). Despite immediate glucocorticoid pulse therapy and plasma exchange (PE) along with antibiotic, he developed cAVB that required temporary pacing on day 2. Because it was thought that hypercytokinemia exacerbated pericarditis, which progressed to myocarditis and cAVB, we decided to PE and cytokine-adsorbing therapy with AN69ST-continuous hemodiafiltration (CHDF). Other than renal failure, his organ dysfunctions improved with the multidisciplinary therapy. CAVB improved and temporary pacing was no longer required on day 11. Even a first-degree atrioventricular block can rapidly progress to cAVB; therefore, strict attention to electrocardiogram is necessary in severe SLE cases. When presenting with organ dysfunctions caused by hypercytokinemia such as severe SLE cases or SLE with severe infection cases, use of the combination of PE and AN69ST-CHDF might be beneficial.

Human thymoproteasome variations influence CD8 T cell selection.

The proteasome is a multi-subunit protease complex essential for housekeeping protein degradation and the production of the major histocompatibility complex (MHC) class I-bound antigen peptides that are essential for recognition by CD8 T cells. MHC variations dramatically contribute to T cell selection and autoimmunity, but genetic variations of peptide processing machinery including proteasome genes have been poorly explored in this context. In the computational analysis of human proteasome gene variation, we documented that PSMB11 was highly enriched for nucleotide changes that interfere with protein function. This gene encodes ß5t, a thymus-specific catalytic subunit that regulates positive selection of CD8 T cells by producing a distinct set of MHC class I-bound peptides. The introduction of PSMB11 variations into the mouse genome by genome-editing revealed that these variations impaired the development of CD8 T cells in vivo. One of the PSMB11 polymorphisms altered the CD8 T cell repertoire in mice and was associated with a higher risk of an autoimmune disease in humans. Our findings suggest that, in addition to the MHC haplotype, proteasome variations influence T cell repertoire selection and may contribute to the difference in individual susceptibility to autoimmunity.

Fezf2 Orchestrates a Thymic Program of Self-Antigen Expression for Immune Tolerance.

Self-tolerance to immune reactions is established via promiscuous expression of tissue-restricted antigens (TRAs) in medullary thymic epithelial cells (mTECs), leading to the elimination of T cells that respond to self-antigens. The transcriptional regulator Aire has been thought to be sufficient for the induction of TRAs, despite some indications that other factors may promote TRA expression in the thymus. Here, we show that the transcription factor Fezf2 directly regulates various TRA genes in mTECs independently of Aire. Mice lacking Fezf2 in mTECs displayed severe autoimmune symptoms, including the production of autoantibodies and inflammatory cell infiltration targeted to peripheral organs. These responses differed from those detected in Aire-deficient mice. Furthermore, Fezf2 expression and Aire expression are regulated by distinct signaling pathways and promote the expression of different classes of proteins. Thus, two independent factors, Fezf2 and Aire, permit the expression of TRAs in the thymus to ensure immune tolerance.