Survival impact of sequential chemotherapy following pembrolizumab for recurrent or metastatic head and neck squamous cell carcinoma.

BACKGROUND: Pembrolizumab alone or combined with chemotherapy is the standard of care for first-line treatment of patients with recurrent or metastatic head and neck squamous cell carcinoma (R/M HNSCC) with positive programmed death-ligand 1 combined positive scores. However, data on second-line chemotherapy following pembrolizumab are scarce. METHODS: A single-center, retrospective study was conducted to determine the efficacies of pembrolizumab and pembrolizumab plus chemotherapy as first-line treatments and the efficacy of second-line chemotherapy for patients with R/M HNSCC who were refractory or intolerant to first-line treatment. RESULTS: Fifty-four patients were treated with pembrolizumab, and 29 received second-line therapy, with 27 opting for cetuximab-containing regimens. The median progression-free survival (PFS), overall survival (OS), and PFS on next-line therapy for first-line treatment were 4.7 (95% confidence interval [CI], 2.1-8.7), 22.1 (95% CI, 12.6-not reached), and 15.6 months (95% CI, 9.7-not reached) in the pembrolizumab group and 5.4 (95% CI, 3.3-6.8), 15.8 (95% CI, 8.6-not reached), and 13.7 months (95% CI, 8.1-not reached) in the pembrolizumab plus chemotherapy group, respectively. The overall response rate and median PFS for second-line treatment were 48.3% (95% CI, 30.4-67.0) and 6.1 months (95% CI, 2.30-8.84). The median OS for patients who received second-line treatment was 18.4 months, which was superior to the median OS of 6.0 months for patients who received the best supportive care (log-rank p = 0.10). CONCLUSION: This study indicates that cetuximab-containing second-line chemotherapy can improve outcomes in R/M HNSCC, even after first-line therapy failure or intolerance.

Translational research in head and neck cancer: Molecular and immunological updates.

Head and neck squamous cell carcinoma (HNSCC) has a poor prognosis. Each year, approximately 880,000 patients are newly diagnosed with HNSCC worldwide, and 450,000 patients with HNSCC die. Risk factors for developing HNSCC have been identified, with cigarette smoking, alcohol consumption, and viral infections being the major factors. Owing to the prevalence of human papillomavirus infection, the number of HNSCC cases is increasing considerably. Surgery and chemoradiotherapy are the primary treatments for HNSCC. With advancements in tumor biology, patients are eligible for novel treatment modalities, namely targeted therapies, immunotherapy, and photoimmunotherapy. Because this area of research has rapidly progressed, clinicians should understand the basic biology of HNSCC to choose an appropriate therapy in the upcoming era of personalized medicine. This review summarized recent developments in tumor biology, focusing on epidemiology, genetic/epigenetic factors, the tumor microenvironment, microbiota, immunity, and photoimmunotherapy in HNSCC, as well as how these findings can be translated into clinical settings.

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.

Utilization of ultrasound in photoimmunotherapy for head and neck cancer: a case report.

BACKGROUND: Photoimmunotherapy (PIT) is a new cancer treatment based on a different mechanism from conventional treatments that combines the administration of a photoabsorber and laser illumination. PIT has two characteristics: a high selectivity of target cells and the possibility of effects beyond the illuminated area. It is thus a potentially effective treatment for a wide variety of cancers. CASE PRESENTATION: We herein report a patient with oropharyngeal squamous cell carcinoma with superficial cervical lymph node recurrence. Intraoperative ultrasound confirmed the localization of the lesion and major vessels near the tumor. We punctured the tumor with catheters of laser illumination under ultrasound guidance. Laser illumination was able to be performed safely without causing serious adverse events, and the effect on the illuminated site was fully exhibited. CONCLUSIONS: Ultrasound allows for the intraoperative monitoring of each localized area and the puncture status easily in real-time. Accurate utilization of ultrasound is extremely important in PIT.

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.

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.

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.

S100A8 enhances IL-1ß production from nasal epithelial cells in eosinophilic chronic rhinosinusitis.

BACKGROUND: Chronic rhinosinusitis is classified into eosinophilic chronic rhinosinusitis (ECRS) and non-eosinophilic chronic rhinosinusitis (NECRS). ECRS is a refractory allergic disease involving a variety of immune and epithelial cells. S100A8 is a damage-associated molecular pattern that is closely related to allergic inflammation. However, the pathological implications of S100A8 in ECRS have not been clarified. METHODS: We evaluated the role of S100A8 in the pathogenesis of ECRS. Gene expression profiles of nasal polyps obtained from patients with ECRS or NECRS were evaluated using RNA sequencing. RESULTS: S100A8 was identified as a significantly upregulated gene in nasal polyps associated with ECRS. Immunohistochemistry consistently revealed intense S100A8 staining in nasal polyps from patients with ECRS. Human nasal epithelial cells expressed the receptor for advanced glycation end products and Toll-like receptor 4. Recombinant S100A8 protein induced interleukin-1ß secretion in human nasal epithelial cells. CONCLUSIONS: Our data demonstrate that S100A8 results in production of interleukin-1ß in the nasal epithelium, which may be involved in the pathogenesis of ECRS.

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.

Lysophosphatidylserines derived from microbiota in Crohn's disease elicit pathological Th1 response.

Microbiota alteration and IFN-γ-producing CD4+ T cell overactivation are implicated in Crohn's disease (CD) pathogenesis. However, it remains unclear how dysbiosis enhances Th1 responses, leading to intestinal inflammation. Here, we identified key metabolites derived from dysbiotic microbiota that induce enhanced Th1 responses and exaggerate colitis in mouse models. Patients with CD showed elevated lysophosphatidylserine (LysoPS) concentration in their feces, accompanied by a higher relative abundance of microbiota possessing a gene encoding the phospholipid-hydrolyzing enzyme phospholipase A. LysoPS induced metabolic reprogramming, thereby eliciting aberrant effector responses in both human and mouse IFN-γ-producing CD4+ T cells. Administration of LysoPS into two mouse colitis models promoted large intestinal inflammation. LysoPS-induced aggravation of colitis was impaired in mice lacking P2ry10 and P2ry10b, and their CD4+ T cells were hyporesponsive to LysoPS. Thus, our findings elaborate on the mechanism by which metabolites elevated in patients with CD harboring dysbiotic microbiota promote Th1-mediated intestinal pathology.

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.

Dysbiosis of Gut Microbiome Is Associated With Rupture of Cerebral Aneurysms.

BACKGROUND AND PURPOSE: Environmental factors are important with respect to the rupture of cerebral aneurysms. However, the relationship between the gut microbiome, an environmental factor, and aneurysm rupture is unclear. Therefore, we compared the gut microbiome in patients with unruptured intracranial aneurysms (UIAs) and ruptured aneurysms (RAs) to identify the specific bacteria causing the rupture of cerebral aneurysms. METHODS: A multicenter, prospective case-control study was conducted over one year from 2019 to 2020. The fecal samples of patients with stable UIAs and RAs immediately after onset were collected. Their gut microbiomes were analyzed using 16S rRNA sequencing. Subsequently, a phylogenetic tree was constructed, and polymerase chain reaction was performed to identify the specific species. RESULTS: A total of 28 RAs and 33 UIAs were included in this study. There was no difference in patient characteristics between RAs and UIAs: age, sex, hypertension, dyslipidemia, diabetes status, body mass index, and smoking. No difference was observed in alpha diversity; however, beta diversity was significantly different in the unweighted UniFrac distances. At the phylum level, the relative abundance of Campylobacter in the RA group was larger than that in the UIA group. Furthermore, the gut microbiome in the RA and UIA groups exhibited significantly different taxonomies. However, Campylobacter was focused on because it is widely known as pathogenic among these bacteria. Then, a phylogenetic tree of operational taxonomic units related to Campylobacter was constructed and 4 species were identified. Polymerase chain reaction for these species identified that the abundance of the genus Campylobacter and Campylobacter ureolyticus was significantly higher in the RA group. CONCLUSIONS: The gut microbiome profile of patients with stable UIAs and RAs were significantly different. The genus Campylobacter and Campylobacter ureolyticus may be associated with the rupture of cerebral aneurysms.

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.

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.

Large-scale plasma-metabolome analysis identifies potential biomarkers of psoriasis and its clinical subtypes.

BACKGROUND: Psoriasis is an immune-mediated skin disease for which the crosstalk between genetic and environmental factors is responsible. To date, no definitive diagnostic criteria for psoriasis yet, and specific biomarkers are required. OBJECTIVE: We performed metabolome analysis to identify metabolite biomarkers of psoriasis and its subtypes such as psoriatic arthritis (PsA) and cutaneous psoriasis (PsC). METHODS: We constructed metabolomics profiling of 130 plasma samples (42 PsA patients, 50 PsC patients, and 38 healthy controls) using a nontargeted metabolomics approach. RESULTS: Psoriasis-control association tests showed that one metabolite (ethanolamine phosphate) was significantly increased in psoriasis samples than in the controls, whereas three metabolites decreased (false discovery rate [FDR] < 0.05; XA0019, nicotinic acid, and 20α-hydroxyprogesterone). In the association test between PsA and PsC, tyramine significantly increased in PsA than in PsC, whereas mucic acid decreased (FDR < 0.05). Molecular pathway analysis of the PsA-PsC association test identified enrichment of vitamin digestion and absorption pathway in PsC (P = 1.3 × 10-4). Correlation network analyses elucidated that a subnetwork centered on aspartate was constructed among the psoriasis-associated metabolites; meanwhile, the major subnetwork among metabolites with differences between PsA and PsC was primarily formed from saturated fatty acids. CONCLUSION: Our large-scale metabolome analysis highlights novel characteristics of plasma metabolites in psoriasis and the differences between PsA and PsC, which could be used as potential biomarkers of psoriasis and its clinical subtypes. These findings contribute to our understanding of psoriasis pathophysiology.