<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2024.1310444</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Exploring the causal relationship between gut microbiota and multiple myeloma risk based on Mendelian randomization and biological annotation</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Feng</surname> <given-names>Zuxi</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1798977/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Liao</surname> <given-names>Minjing</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Bai</surname> <given-names>Jun</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2343839/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Li</surname> <given-names>Yanhong</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Chen</surname> <given-names>Yue</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhang</surname> <given-names>Li</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2158485/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Guo</surname> <given-names>Xuege</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes" equal-contrib="yes">
<name><surname>Li</surname> <given-names>Lijuan</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<xref ref-type="author-notes" rid="fn002"><sup>&#x2020;</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes" equal-contrib="yes">
<name><surname>Zhang</surname> <given-names>Liansheng</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c002"><sup>&#x002A;</sup></xref>
<xref ref-type="author-notes" rid="fn002"><sup>&#x2020;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2403844/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Hematology, Lanzhou University Second Hospital</institution>, <addr-line>Lanzhou</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>Second Clinical Medical College, Lanzhou University</institution>, <addr-line>Lanzhou</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>Key Laboratory of the Hematology of Gansu Province, Lanzhou University Second Hospital, Lanzhou University</institution>, <addr-line>Lanzhou</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Mamoru Yamada, Yamaguchi University, Japan</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Shintaro Maeno, Yamaguchi University, Japan</p><p>Ying Liao, Sun Yat-sen University Cancer Center (SYSUCC), China</p><p>Ruizheng Sun, Central South University, China</p></fn>
<corresp id="c001">&#x002A;Correspondence: Lijuan Li, <email>doctorjuan@sina.com</email></corresp>
<corresp id="c002">Liansheng Zhang, <email>doctorzhanglsh@sina.com</email></corresp>
<fn fn-type="equal" id="fn002"><p><sup>&#x2020;</sup>These authors have contributed equally to this work</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>12</day>
<month>02</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>15</volume>
<elocation-id>1310444</elocation-id>
<history>
<date date-type="received">
<day>09</day>
<month>10</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>15</day>
<month>01</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2024 Feng, Liao, Bai, Li, Chen, Zhang, Guo, Li and Zhang.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Feng, Liao, Bai, Li, Chen, Zhang, Guo, Li and Zhang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>The microbial genome-wide association studies (mbGWAS) have highlighted significant host-microbiome interactions based on microbiome heritability. However, establishing causal relationships between particular microbiota and multiple myeloma (MM) remains challenging due to limited sample sizes.</p>
</sec>
<sec>
<title>Methods</title>
<p>Gut microbiota data from a GWAS with 18,340 participants and MM summary statistics from 456,348 individuals. The inverse variance-weighted (IVW) method was used as the main bidirectional Mendelian randomization (MR) analysis. To assess the robustness of our results, we further performed supplementary analyses, including MR pleiotropy residual sum and outlier (MR-PRESSO) test, MR-Egger, Weighted median, Simple mode, and Weighted mode. Moreover, a backward MR analysis was conducted to investigate the potential for reverse causation. Finally, gene and gene-set-based analyses were then conducted to explore the common biological factors connecting gut microbiota and MM.</p>
</sec>
<sec>
<title>Results</title>
<p>We discovered that 10 gut microbial taxa were causally related to MM risk. Among them, family <italic>Acidaminococcaceae</italic>, <italic>Bacteroidales</italic> family S24-7, family <italic>Porphyromonadaceae</italic>, genus <italic>Eubacterium ruminantium group</italic>, genus <italic>Parabacteroides</italic>, and genus <italic>Turicibacter</italic> were positively correlated with MM. Conversely, class <italic>Verrucomicrobia</italic>, family <italic>Verrucomicrobiaceae</italic>, genus <italic>Akkermansia</italic>, and order <italic>Verrucomicrobiales</italic> were negatively correlated with MM. The heterogeneity test revealed no Heterogeneity. MR-Egger and MR-PRESSO tests showed no significant horizontal pleiotropy. Importantly, leave-one-out analysis confirmed the robustness of MR results. In the backward MR analysis, no statistically significant associations were discovered between MM and 10 gut microbiota taxa. Lastly, we identified novel host-microbiome shared genes (AUTS2, CDK2, ERBB3, IKZF4, PMEL, SUOX, and RAB5B) that are associated with immunoregulation and prognosis in MM through biological annotation.</p>
</sec>
<sec>
<title>Discussion</title>
<p>Overall, this study provides evidence supporting a potential causal relationship between gut microbiota and MM risk, while also revealing novel host-microbiome shared genes relevant to MM immunoregulation and clinical prognosis.</p>
</sec>
</abstract>
<kwd-group>
<kwd>multiple myeloma</kwd>
<kwd>gut microbiota</kwd>
<kwd>Mendelian randomization</kwd>
<kwd>biological annotation</kwd>
<kwd>immunoregulation</kwd>
</kwd-group>
<counts>
<fig-count count="7"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="58"/>
<page-count count="14"/>
<word-count count="8132"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Evolutionary and Genomic Microbiology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="S1" sec-type="intro">
<title>1 Introduction</title>
<p>Multiple myeloma (MM) is a clonal plasma cell malignant tumor originating in the bone marrow, characterized by the excessive production of monoclonal immunoglobulins (<xref ref-type="bibr" rid="B19">Joshua et al., 2019</xref>). Its global incidence has steadily increased, making it the second most prevalent blood cancer after leukemia. Importantly, MM is not a single entity but rather a heterogeneous disease (<xref ref-type="bibr" rid="B31">Owens, 2020</xref>). Clinical presentations vary from asymptomatic cases to life-threatening conditions, with common manifestations including fatigue, bone pain, hypercalcemia with renal insufficiency, anemia, and hyperproteinemia (<xref ref-type="bibr" rid="B2">Angtuaco et al., 2004</xref>). Diagnosis relies on multidisciplinary approaches encompassing orthopedics, radiology, nuclear medicine, hematology, and oncology (<xref ref-type="bibr" rid="B7">Caers et al., 2018</xref>). In recent decades, patient survival rates have significantly increased due to the implementation of proteasome inhibitors and immunomodulatory drugs. However, these treatments are not curative, and MM often recurs, necessitating further intervention (<xref ref-type="bibr" rid="B21">Laubach et al., 2015</xref>). Consequently, the quest for optimized treatment strategies for MM remains a significant long-term challenge.</p>
<p>The human intestine hosts a diverse and intricate microbial community known as gut microbiota, which derives from various sources, including dietary intake and probiotic supplements. These microbiotas are essential for maintaining mucosal barrier integrity and offering numerous other health benefits (<xref ref-type="bibr" rid="B3">Badgeley et al., 2021</xref>). Currently, gut microbiota is recognized as both a contributor to and a safeguard against various human diseases, including its intricate involvement in the pathogenesis of MM (<xref ref-type="bibr" rid="B33">Park et al., 2022</xref>). Conversely, cancer can influence the host&#x2019;s gut microbiota, potentially leading to microbiota disruption and tumor progression (<xref ref-type="bibr" rid="B51">Yonekura et al., 2022</xref>). The gut microbiota within a human host represents a partially heritable phenotype, and correlations exist between the host&#x2019;s genotype and gut microbiota variations (<xref ref-type="bibr" rid="B12">Goodrich et al., 2014</xref>). The interactions between hosts and microorganisms are intricate, with distinct mechanisms governing gut microbiota&#x2019;s impact on host health and disease. Current scientific research is rapidly advancing toward enhancing cancer treatment strategies by assessing human microbial composition and function and subsequently implementing targeted regulation (<xref ref-type="bibr" rid="B33">Park et al., 2022</xref>). Although limited sample sizes have hindered research in this area, resulting in limited published mbGWAS studies, mbGWAS have still revealed significant insights (<xref ref-type="bibr" rid="B26">Lopera-Maya et al., 2022</xref>). Over the past 5 years, research into the relationship between gut microbiota and MM has steadily progressed. The latest research reveals that in MM patients, the gut microbiome&#x2019;s enrichment of <italic>Citrobacter freundii</italic> induces drug resistance by increasing ammonium levels through the transmembrane channel protein SLC12A2 (<xref ref-type="bibr" rid="B58">Zhu et al., 2024</xref>). Simultaneously, a potential therapeutic strategy is proposed, utilizing furosemide sodium to inhibit ammonium uptake and enhance treatment efficacy. Gut microbiota diversity in autologous hematopoietic cell transplant recipients is initially reduced, resembling allogeneic transplants (<xref ref-type="bibr" rid="B20">Khan et al., 2021</xref>). Higher diversity during transplantation is associated with a lower risk of progression or death, highlighting a potential connection between gut microbiota and patient outcomes. Gut microbiota alterations and changes in short-chain fatty acid levels are associated with the progression of multiple myeloma, indicating potential therapeutic targets and predictors of treatment response (<xref ref-type="bibr" rid="B37">Rodriguez-Garcia et al., 2023</xref>). Thus, understanding the role of gut microbiota in MM is crucial, as investigating the causal relationship between gut microbiota and MM holds significant implications for advancing our knowledge of disease mechanisms and potential therapeutic interventions.</p>
<p>In recent years, the role of immunology in tumor pathogenesis has gained widespread attention, introducing a new era in cancer treatment through immunotherapy. Medical research has elucidated the vital role of microbiota in regulating tumor immune surveillance and participating in cancer pathogenesis and progression (<xref ref-type="bibr" rid="B23">Liu et al., 2021</xref>). The latest findings demonstrate that fiber-rich diets and fermented foods can modulate microbial composition, activity, and host immune status (<xref ref-type="bibr" rid="B45">Wastyk et al., 2021</xref>). Consequently, gut microbiota-targeted immunotherapy has recently surfaced as an innovative approach in treating cancer (<xref ref-type="bibr" rid="B57">Zhou et al., 2021</xref>). <xref ref-type="bibr" rid="B8">Calcinotto et al. (2018)</xref> discovered that the gut microbiota plays a role in the progression of MM by affecting the differentiation and migration of Th17 cells in the bone marrow. Furthermore, an immunological connection between the gut and the transition from asymptomatic to symptomatic MM has been confirmed (<xref ref-type="bibr" rid="B8">Calcinotto et al., 2018</xref>). In MM patients, the gut microbiome exerts influence not only on disease progression but also on treatment response and treatment-related toxicity (<xref ref-type="bibr" rid="B1">Ahmed et al., 2020</xref>). Recent research highlights microbiota&#x2019;s role in promoting inflammation and influencing MM&#x2019;s disease development, including transplantation (<xref ref-type="bibr" rid="B35">Pianko and Golob, 2022</xref>).</p>
<p>Despite the challenges of obtaining real-time microbiome data in clinical practice, retrospective studies have started illuminating the relationship between changes in particular microbial characteristic and MM (<xref ref-type="bibr" rid="B56">Zhang et al., 2022</xref>). Due to confounding factors such as lifestyle, environment, host gene mutations, and potential reverse causality, the precise association between genes in MM patients and gut microbiota remains underexplored. Mendelian Randomization (MR) using two-sample analysis is a robust approach that leverages inherited genetic variations as instrumental variables (IVs) to investigate causal relationships between exposure and outcome traits. MR leverages large GWAS data to study the associations between risk factors and diseases while effectively controlling non-genetic confounding variables, leading to more precise and reliable results. Currently, MR is extensively utilized for investigating the causal link between gut microbiota and diseases (<xref ref-type="bibr" rid="B40">Shi et al., 2023</xref>; <xref ref-type="bibr" rid="B54">Zeng Y. et al., 2023</xref>). Although current research results suggest a potential interaction between gut microbiota and MM, there is still considerable research potential. Thus, we aim to delve deeper and comprehensively explore the causal relationship between gut microbiota and MM through MR analysis. Identifying therapeutic targets, promoting precision medicine, and uncovering prognostic indicators are crucial aspects in the diagnosis and treatment of MM.</p>
<p>In our study, we utilized bidirectional two-sample MR analyses to investigate potential causal links between microbial characteristics and MM. Furthermore, we offered a biological annotation for the noteworthy connections between gut microbiota and MM.</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>2 Materials and methods</title>
<sec id="S2.SS1">
<title>2.1 Study design</title>
<p><xref ref-type="fig" rid="F1">Figure 1</xref> depicts the comprehensive study design. In summary, our study involved extracting genetic variants linked to the respective exposure by analyzing GWAS summary statistics. We meticulously selected instrumental variables based on stringent criteria. We then conducted a bidirectional two-sample MR analysis using five distinct MR methods in a sequential manner. Next, a comprehensive series of sensitivity analyses, comprising the examination of heterogeneity, pleiotropy, and the leave-one-out test, were conducted to scrutinize the robustness of statistically significant associations. Lastly, we explored positional mapped genes and gene-set-based analyses using biological annotation.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>Visual depiction of the research study. GWAS, genome-wide association study; IV, instrumental variables; MR, Mendelian randomization; IVW, inverse-variance weighted; MR-PRESSO, MR pleiotropy residual sum and outlier; GEO, gene expression omnibus.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1310444-g001.tif"/>
</fig>
</sec>
<sec id="S2.SS2">
<title>2.2 Data sources</title>
<p>In this study, we acquired summary statistics for human MM from a previously published GWAS.<sup><xref ref-type="fn" rid="footnote1">1</xref></sup> This GWAS incorporated data from the UK Biobank and encompassed 456,348 individuals with European ancestry (<xref ref-type="bibr" rid="B17">Jiang et al., 2021</xref>). This dataset includes 564 cases of European ancestry and 455,784 controls of European ancestry. Additionally, summary statistics concerning the human gut microbiome in this study were derived from a GWAS involving 18,340 participants.<sup><xref ref-type="fn" rid="footnote2">2</xref></sup> Our study encompasses 24 cohorts, samples from single-ancestry cohorts include European (16 cohorts, <italic>N</italic> = 13,266), Middle-Eastern (1 cohort, <italic>N</italic> = 481), East Asian (1 cohort, <italic>N</italic> = 811), American Hispanic/Latin (1 cohort, <italic>N</italic> = 1,097), and African American (1 cohort, <italic>N</italic> = 114). Additionally, four cohorts were multi-ancestry (<italic>N</italic> = 2,571). The gene expression data for 138 samples were downloaded from the Gene Expression Omnibus (GEO)<sup><xref ref-type="fn" rid="footnote3">3</xref></sup> database with accession number <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="GSE13591">GSE13591</ext-link>. To ascertain the relationship between the expression of hub genes and the overall survival of patients with MM, we utilized the Kaplan-Meier plotter database.<sup><xref ref-type="fn" rid="footnote4">4</xref></sup> Note: The acquisition dates for databases used in this study are as follows: MM GWAS data (August 20, 2023), gut microbiome data (September 7, 2023), GEO data (September 10, 2023), and Kaplan-Meier plotter database (September 12, 2023).</p>
</sec>
<sec id="S2.SS3">
<title>2.3 Instrumental variable (IV)</title>
<p>At first, we omitted 15 microbial characteristics that did not have specific names, resulting in 196 microbial features. Later, we chose IVs using a less strict significance threshold of <italic>p</italic> &#x003C; 1.0 &#x00D7; 10<sup>&#x2013;5</sup> (<xref ref-type="bibr" rid="B24">Liu et al., 2023</xref>; <xref ref-type="bibr" rid="B53">Zeng C. et al., 2023</xref>), which was determined to be the best <italic>p</italic>-value threshold for choosing genetic predictors that account for a greater amount of variance in the outcome. Additionally, we employed several alternative thresholds 5.0 &#x00D7; 10<sup>&#x2013;6</sup> (<xref ref-type="bibr" rid="B14">He et al., 2022</xref>; <xref ref-type="bibr" rid="B30">Nolde et al., 2022</xref>) and 5.0 &#x00D7; 10<sup>&#x2013;8</sup> (<xref ref-type="bibr" rid="B32">Papiol et al., 2021</xref>; <xref ref-type="bibr" rid="B52">Yuan and Larsson, 2022</xref>) to choose IVs for gut microbial (<xref ref-type="supplementary-material" rid="TS1">Supplementary Table 1</xref>). To prevent biased results stemming from linkage disequilibrium (LD), we performed clumping analysis on European samples from the 1000 Genomes project using the TwoSampleMR package, with stringent criteria (<italic>R</italic><sup>2</sup> &#x003C; 0.001 and a window size of 10,000 kb). With an <italic>R</italic><sup>2</sup> value of less than 0.001, we retained only the single nucleotide polymorphism (SNP) possessing the lowest <italic>p</italic>-value. Additionally, we excluded SNPs not found in the LD reference panel. To gauge the strength of the IVs, we computed F-statistics using the provided equation <italic>R</italic><sup>2</sup> = 2 &#x00D7; MAF &#x00D7; (1-MAF) &#x00D7; &#x03B2;<sup>2</sup>, <italic>F</italic> = <italic>R</italic><sup>2</sup> &#x00D7; (N&#x2013;1&#x2013;k)/ k (1&#x2013;<italic>R</italic><sup>2</sup>), where <italic>R</italic><sup>2</sup> represents the variance explained by the IVs, &#x201C;N&#x201D; denotes the data sample size, and &#x201C;k&#x201D; represents the count of SNPs included in the instrument (<xref ref-type="bibr" rid="B6">Bowden et al., 2016</xref>). Ultimately, we opted for independent SNPs linked to MM as IVs, employing a threshold of 1.0 &#x00D7; 10<sup>&#x2013;5</sup> (<xref ref-type="supplementary-material" rid="TS1">Supplementary Table 2</xref>), which we utilized in subsequent MR analyses. In the backward MR analysis, we chose IVs linked to MM at various cutoff levels (1.0 &#x00D7; 10<sup>&#x2013;5</sup>, 5.0 &#x00D7; 10<sup>&#x2013;6</sup>, and 5.0 &#x00D7; 10<sup>&#x2013;8</sup>) (<xref ref-type="supplementary-material" rid="TS1">Supplementary Table 3</xref>). In contrast to the forward analyses, the results at a threshold of 5.0 &#x00D7; 10<sup>&#x2013;8</sup> as the discovery set, while treating the other thresholds as validation criteria.</p>
</sec>
<sec id="S2.SS4">
<title>2.4 Statistical analysis</title>
<p>In this study, we employed various approaches to investigate potential causal connections between gut microbiota and MM, such as MR-PRESSO, IVW, weighted median, MR-Egger regression, simple mode, and weighted model. Our primary MR analysis relied on the IVW method (<xref ref-type="bibr" rid="B50">Yavorska and Burgess, 2017</xref>). The IVW method integrates the Delta method and Wald estimates for every SNP, producing a holistic assessment of the gut microbiota&#x2019;s influence on MM (<xref ref-type="bibr" rid="B41">Song et al., 2023</xref>). Without horizontal pleiotropy, IVW results stay unbiased. We assessed the presence of overall horizontal pleiotropy by measuring genetic variant heterogeneity using the Q statistic and examining the intercept in the MR-Egger test to validate MR assumptions (<xref ref-type="bibr" rid="B5">Bowden et al., 2015</xref>). To address horizontal pleiotropy, we conducted sensitivity assessments employing simple mode, weighted median, and weighted mode. We also used the MR-PRESSO approach to conduct global and outlier tests in order to detect potential outliers. Subsequently, we obtained corrected association results by excluding these potential outliers (<xref ref-type="bibr" rid="B43">Verbanck et al., 2018</xref>). To assess the potential prevalence of particular variants in the process of estimating causality, we conducted leave-one-out examinations. This involved systematically omitting each SNP from the evaluation and subsequently reassessing the causative connection (<xref ref-type="bibr" rid="B9">Corbin et al., 2016</xref>).</p>
<p>Additionally, we explored the potential impact of directional pleiotropy by investigating the secondary phenotypes associated with each of the SNPs utilized as instrumental variables. This analysis was conducted using the PhenoScanner database, and the data was accessed on August 30, 2023, via the following link: <ext-link ext-link-type="uri" xlink:href="http://www.phenoscanner.medschl.cam.ac.uk/">http://www.phenoscanner.medschl.cam.ac.uk/</ext-link>. After confirming the connection between microbiota features and MM through MR methods, we proceeded with backward MR analyses. In this analysis, we examined MM as the exposure variable and assessed microbial features as the outcomes. &#x201C;Mendelian Randomization,&#x201D; &#x201C;TwoSampleMR,&#x201D; and &#x201C;MRPRESSO&#x201D; packages in the open-source statistical software R (version 4.2.3) were employed for our MR analyses. The connections between the human gut microbiota and MM risk were quantified using odds ratios (ORs) accompanied by 95% confidence intervals (CIs). Statistical significance was confirmed when two-sided <italic>p</italic>-values were below 0.05.</p>
</sec>
<sec id="S2.SS5">
<title>2.5 Biological annotation</title>
<p>Initially, we annotated the genes mapped to instrumental variants&#x2019; positions between gut microbiota and MM utilizing the GWAS catalog provided with FUMA (<xref ref-type="bibr" rid="B48">Welter et al., 2014</xref>). We then employed integrated hypergeometric tests within the FUMA platform, further enhancing our understanding of the previously reported associations between mapped genes and various phenotypes (<xref ref-type="bibr" rid="B47">Watanabe et al., 2017</xref>). Subsequently, we constructed protein&#x2013;protein interaction (PPI) networks using genes mapped to instrumental variants for both MM and our MR analysis results for mbGWAS. We utilized STRING version 12.0<sup><xref ref-type="fn" rid="footnote5">5</xref></sup> and Cytoscape (version 3.8.2) to obtain and visualize the connected PPI network. We have determined the top ten hub genes through the Maximal Clique Centrality (MCC) approach, given that genes with numerous interconnecting links are pivotal for maintaining network stability (<xref ref-type="bibr" rid="B49">Xiao et al., 2023</xref>). Furthermore, we examined the gene expression profiles of central genes within bulk tissue through the FUMA portal. We assessed the impact of these hub genes on a wide range of characteristics by conducting phenome-wide association studies (PheWAS). We explored their pleiotropic effects in summary data for 4,756 complex traits and illnesses across 28 different domains, utilizing the GWAS ATLAS (<xref ref-type="bibr" rid="B46">Watanabe et al., 2019</xref>). To perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis for the hub genes, we used the &#x201C;ClusterProfiler&#x201D; and &#x201C;org.Hs.eg.db&#x201D; package in R software (version 4.2.3). In this study, we employed both <italic>p</italic>-value = 0.05 and <italic>q</italic>-value = 0.05 as filtering criteria for the analysis of GO terms and KEGG pathways. Additionally, we utilized the Kaplan-Meier plotter database to establish the association between hub gene expression and the prognosis of MM patients. Finally, we retrieved the expression data of hub genes from the GEO database to investigate their immune regulatory functions in MM.</p>
</sec>
</sec>
<sec id="S3" sec-type="results">
<title>3 Results</title>
<sec id="S3.SS1">
<title>3.1 Forward MR analyses</title>
<p>In our forward MR analyses, we employed the IVW method to successfully identified potential causal links between 10 specific bacterial traits and MM risk (<xref ref-type="fig" rid="F2">Figure 2</xref> and <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 4</xref>). Notably, we discovered positive relationships between MM risk and the following: family <italic>Acidaminococcaceae</italic> (OR, 2.3947; 95% CI, 1.2937&#x2013;4.4326; <italic>P</italic> = 0.0054), <italic>Bacteroidales</italic> family S24-7 (OR, 1.7964; 95% CI, 1.0797&#x2013;2.9890; <italic>P</italic> = 0.0241), family <italic>Porphyromonadaceae</italic> (OR, 2.7848; 95% CI, 1.1246&#x2013;6.8959; <italic>P</italic> = 0.0268), genus <italic>Eubacterium ruminantium group</italic> (OR, 1.4440; 95% CI, 1.0016&#x2013;2.0817; <italic>P</italic> = 0.0490), genus <italic>Parabacteroides</italic> (OR, 2.6896; 95% CI, 1.1425&#x2013;6.3316; <italic>P</italic> = 0.0235), and genus <italic>Turicibacter</italic> (OR, 1.6478; 95% CI, 1.0061&#x2013;2.6986; <italic>P</italic> = 0.0472) (<xref ref-type="fig" rid="F2">Figure 2</xref>). On the contrary, we found a negative association between MM risk and the following: class <italic>Verrucomicrobia</italic> (OR, 0.5476; 95% CI, 0.3216&#x2013;0.9322; <italic>P</italic> = 0.0265), family <italic>Verrucomicrobiaceae</italic> (OR, 0.5476; 95% CI, 0.3216&#x2013;0.9324; <italic>P</italic> = 0.0266), genus <italic>Akkermansia</italic> (OR, 0.5476; 95% CI, 0.3216&#x2013;0.9324; <italic>P</italic> = 0.0266), and order <italic>Verrucomicrobiales</italic> (OR, 0.5476; 95% CI, 0.3216&#x2013;0.9322; <italic>P</italic> = 0.0265) (<xref ref-type="fig" rid="F2">Figure 2</xref>). Furthermore, the association between <italic>Bacteroidales</italic> family S24-7 and MM remained consistent when applying the weighted median approach (OR, 2.0346; 95% CI, 1.0714&#x2013;3.8635; <italic>P</italic> = 0.0300). The scatter plot and funnel plot were displayed <xref ref-type="supplementary-material" rid="FS1">Supplementary Figures 1</xref>, <xref ref-type="supplementary-material" rid="FS1">2</xref>. Comprehensive genetic instruments for evaluating causal effects can be found in <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 2</xref>. Notably, The lowest F-statistic observed among the instruments stood at 19.5, signifying robust associations between all IVs and the microbiome features. Results from the MR-PRESSO test corroborated our findings (<xref ref-type="table" rid="T1">Table 1</xref>). Our findings did not suggest the existence of horizontal pleiotropy in the intercept of the MR-Egger regression analysis (<xref ref-type="table" rid="T1">Table 1</xref>). Additionally, heterogeneity tests showed no noteworthy heterogeneity among these independent variables (<xref ref-type="table" rid="T1">Table 1</xref>). Importantly, leave-one-out analyses failed to identify any highly influential variants among the instrumental SNPs (<xref ref-type="fig" rid="F3">Figure 3</xref>). In an effort to address potential pleiotropy, we scanned all SNPs used as IVs in our study using the PhenoScanner database. This led to the identification of 8 SNPs associated with secondary traits, and relevant literature was reviewed to explore their connections to MM (<xref ref-type="supplementary-material" rid="TS1">Supplementary Table 5</xref>). While smoking is associated with various cancers, its link to MM remains unclear (<xref ref-type="bibr" rid="B36">Psaltopoulou et al., 2013</xref>). For the other traits in the table, there is currently no substantial evidence linking them to MM. Nevertheless, we conducted MR analyses using the IVW method after excluding these pleiotropic SNPs. Our results continued to support associations between family <italic>Acidaminococcaceae</italic> (OR, 2.2100; 95% CI, 1.1498&#x2013;4.2479; <italic>P</italic> = 0.0174), <italic>Bacteroidales</italic> family S24-7 (OR, 1.7948; 95% CI, 1.0779&#x2013;2.9884; <italic>P</italic> = 0.0245), family <italic>Porphyromonadaceae</italic> (OR, 3.71309; 95% CI, 1.2368&#x2013;11.1471; <italic>P</italic> = 0.01934), and genus <italic>Turicibacter.id.2162</italic> (OR, 1.6475; 95% CI, 1.0041&#x2013;2.7032; <italic>P</italic> = 0.0481) with MM risk did not change significantly in the IVW method. However, the relationships between family <italic>Verrucomicrobiaceae</italic> (OR, 0.5891; 95% CI, 0.3383&#x2013;1.0257; <italic>P</italic> = 0.0614), genus <italic>Akkermansia.id.4037</italic> (OR, 0.5892; 95% CI, 0.3384&#x2013;1.0259; <italic>P</italic> = 0.0615), and order <italic>Verrucomicrobiales</italic> (OR, 0.5891; 95% CI, 0.3383&#x2013;1.0256; <italic>P</italic> = 0.0614) with MM appeared less stable. These three unstable associations were influenced by a common SNP (rs4936098), but direct evidence linking this SNP&#x2019;s traits to MM is lacking. Detailed information regarding the results of MR analyses with updated SNPs are provided in <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 6</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p>Forest plot depicting the forward MR findings. This forest plot illustrates potential causal connections between microbial characteristics and MM. The information is displayed in the form of odds ratios (OR) accompanied by their respective 95% confidence intervals (CI). nSNP, number of single nucleotide polymorphism.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1310444-g002.tif"/>
</fig>
<table-wrap position="float" id="T1">
<label>TABLE 1</label>
<caption><p>Effect estimates of the relationships between MM risk and gut microbiotas in the MR analyses.</p></caption>
<table cellspacing="5" cellpadding="5" frame="box" rules="all">
<thead>
<tr>
<td valign="top" align="left" style="color:#ffffff;background-color: #7f8080;">Gut microbiota</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;">IVW</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;">MR-PRESSO<break/>test (IVW)</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;">Pleiotropy<break/>Test (IVW)</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;">Heterogeneity<break/>Test (MR-Egger)</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Class <italic>Verrucomicrobia</italic></td>
<td valign="top" align="center">0.0265</td>
<td valign="top" align="center">0.8560</td>
<td valign="top" align="center">0.6543718</td>
<td valign="top" align="center">0.8330287</td>
</tr>
<tr>
<td valign="top" align="left">Order <italic>Verrucomicrobiales</italic></td>
<td valign="top" align="center">0.0265</td>
<td valign="top" align="center">0.8536</td>
<td valign="top" align="center">0.6543718</td>
<td valign="top" align="center">0.8330287</td>
</tr>
<tr>
<td valign="top" align="left">Family <italic>Acidaminococcaceae</italic></td>
<td valign="top" align="center">0.0054</td>
<td valign="top" align="center">0.7668</td>
<td valign="top" align="center">0.9276273</td>
<td valign="top" align="center">0.7405918</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Bacteroidales</italic> family S24-7</td>
<td valign="top" align="center">0.0241</td>
<td valign="top" align="center">0.9683</td>
<td valign="top" align="center">0.2707392</td>
<td valign="top" align="center">0.9649069</td>
</tr>
<tr>
<td valign="top" align="left">Family <italic>Porphyromonadaceae</italic></td>
<td valign="top" align="center">0.0268</td>
<td valign="top" align="center">0.5729</td>
<td valign="top" align="center">0.8650975</td>
<td valign="top" align="center">0.5399784</td>
</tr>
<tr>
<td valign="top" align="left">Family <italic>Verrucomicrobiaceae</italic></td>
<td valign="top" align="center">0.0266</td>
<td valign="top" align="center">0.8528</td>
<td valign="top" align="center">0.6532106</td>
<td valign="top" align="center">0.8328674</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Eubacterium ruminantium group</italic></td>
<td valign="top" align="center">0.0490</td>
<td valign="top" align="center">0.9377</td>
<td valign="top" align="center">0.6252739</td>
<td valign="top" align="center">0.9310561</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Akkermansia</italic></td>
<td valign="top" align="center">0.0266</td>
<td valign="top" align="center">0.8515</td>
<td valign="top" align="center">0.6523758</td>
<td valign="top" align="center">0.8327713</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Parabacteroides</italic></td>
<td valign="top" align="center">0.0235</td>
<td valign="top" align="center">0.5564</td>
<td valign="top" align="center">0.5362979</td>
<td valign="top" align="center">0.4810427</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Turicibacter</italic></td>
<td valign="top" align="center">0.0472</td>
<td valign="top" align="center">0.3613</td>
<td valign="top" align="center">0.9869001</td>
<td valign="top" align="center">0.3164009</td>
</tr>
</tbody>
</table></table-wrap>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p>Leave-one-out analysis for MR results. The leave-one-out analysis was conducted by iteratively recalculating MR estimates through the IVW method, systematically omitting one SNP at each iteration. This process aimed to investigate whether any individual SNP with a substantial horizontal pleiotropic effect could significantly impact the MR estimates.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1310444-g003.tif"/>
</fig>
</sec>
<sec id="S3.SS2">
<title>3.2 Backward MR analyses</title>
<p>We further conducted backward MR analyses to investigate potential inverse associations between ten microbial characteristics and MM. Utilizing the IVW method, our analysis failed to reveal any substantiating evidence regarding the similar cause-and-effect relationships between MM and the ten microbial characteristics. We employed three distinct thresholds (SNP <italic>p</italic>-values less than 1.0 &#x00D7; 10<sup>&#x2013;5</sup>, 5.0 &#x00D7; 10<sup>&#x2013;6</sup>, or 5.0 &#x00D7; 10<sup>&#x2013;8</sup>) to choose IVs linked to MM. These three thresholds are commonly used in current MR analysis for SNP selection. Here, we opted for a lower threshold of 1.0 &#x00D7; 10<sup>&#x2013;5</sup> and a higher threshold of 5.0 &#x00D7; 10<sup>&#x2013;8</sup> for analysis and presentation of results. With a threshold of 5.0 &#x00D7; 10<sup>&#x2013;8</sup>, only 1 SNP remained, with no statistically meaningful correlations were observed between MM and any of the ten bacterial traits (<xref ref-type="table" rid="T2">Table 2</xref>). Our results indicate a lack of substantial evidence supporting the notion that MM leads to changes in microbial characteristics using thresholds of 5.0 &#x00D7; 10<sup>&#x2013;8</sup> and 1.0 &#x00D7; 10<sup>&#x2013;5</sup> (<xref ref-type="table" rid="T2">Table 2</xref> and <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 7</xref>). The findings remained consistent across sensitivity analyses, as outlined in <xref ref-type="table" rid="T2">Table 2</xref> and <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 7</xref>.</p>
<table-wrap position="float" id="T2">
<label>TABLE 2</label>
<caption><p>Effect estimates of the relationships between MM and gut microbiota were examined in the backward MR analyses.</p></caption>
<table cellspacing="5" cellpadding="5" frame="box" rules="all">
<thead>
<tr>
<td valign="top" align="left" style="color:#ffffff;background-color: #7f8080;">Gut microbiota</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;">Threshold</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;">Methods</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;">nSNP</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;">OR</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;">95% CI</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;"><italic>p</italic>-value</td>
<td valign="top" align="center" style="color:#ffffff;background-color: #7f8080;"><italic>p</italic>-value (Intercept)</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Class <italic>Verrucomicrobia</italic></td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">1.0301</td>
<td valign="top" align="center">0.9568&#x2013;1.1091</td>
<td valign="top" align="center">0.4303</td>
<td valign="top" align="center">/</td>
</tr>
<tr>
<td valign="top" align="left">Order <italic>Verrucomicrobiales</italic></td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">1.0301</td>
<td valign="top" align="center">0.9568&#x2013;1.1091</td>
<td valign="top" align="center">0.4303</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">Family <italic>Acidaminococcaceae</italic></td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0.9911</td>
<td valign="top" align="center">0.9222&#x2013;1.0652</td>
<td valign="top" align="center">0.808</td>
<td valign="top" align="center">/</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Bacteroidales</italic> family S24-7</td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0.9667</td>
<td valign="top" align="center">0.8831&#x2013;1.0581</td>
<td valign="top" align="center">0.4623</td>
<td valign="top" align="center">/</td>
</tr>
<tr>
<td valign="top" align="left">Family <italic>Porphyromonadaceae</italic></td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0.9748</td>
<td valign="top" align="center">0.9180&#x2013;1.0351</td>
<td valign="top" align="center">0.4048</td>
<td valign="top" align="center">/</td>
</tr>
<tr>
<td valign="top" align="left">Family <italic>Verrucomicrobiaceae</italic></td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">1.0301</td>
<td valign="top" align="center">0.9568&#x2013;1.1091</td>
<td valign="top" align="center">0.4303</td>
<td valign="top" align="center">/</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Eubacterium ruminantium group</italic></td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">1.0034</td>
<td valign="top" align="center">0.9156&#x2013;1.0995</td>
<td valign="top" align="center">0.9427</td>
<td valign="top" align="center">/</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Akkermansia</italic></td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">1.0308</td>
<td valign="top" align="center">0.9574&#x2013;1.1098</td>
<td valign="top" align="center">0.4207</td>
<td valign="top" align="center">/</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Parabacteroides</italic></td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0.9859</td>
<td valign="top" align="center">0.9273&#x2013;1.0482</td>
<td valign="top" align="center">0.6494</td>
<td valign="top" align="center">/</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Turicibacter</italic></td>
<td valign="top" align="center">5 &#x00D7; 10<sup>&#x2212;8</sup></td>
<td valign="top" align="center">Wald ratio</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">1.0162</td>
<td valign="top" align="center">0.9335&#x2013;1.1062</td>
<td valign="top" align="center">0.7104</td>
<td valign="top" align="center">/</td>
</tr>
<tr>
<td valign="top" align="left">Class <italic>Verrucomicrobia</italic></td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">1.0001</td>
<td valign="top" align="center">0.9672&#x2013;1.0341</td>
<td valign="top" align="center">0.9959</td>
<td valign="top" align="center">0.0741</td>
</tr>
<tr>
<td valign="top" align="left">Order <italic>Verrucomicrobiales</italic></td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">1.0001</td>
<td valign="top" align="center">0.9672&#x2013;1.0341</td>
<td valign="top" align="center">0.9959</td>
<td valign="top" align="center">0.0741</td>
</tr>
<tr>
<td valign="top" align="left">Family <italic>Acidaminococcaceae</italic></td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">1.0121</td>
<td valign="top" align="center">0.9827&#x2013;1.0424</td>
<td valign="top" align="center">0.4249</td>
<td valign="top" align="center">0.3991</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Bacteroidales</italic> family S24-7</td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">0.9867</td>
<td valign="top" align="center">0.9510&#x2013;1.0237</td>
<td valign="top" align="center">0.4745</td>
<td valign="top" align="center">0.694</td>
</tr>
<tr>
<td valign="top" align="left">Family <italic>Porphyromonadaceae</italic></td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">0.9992</td>
<td valign="top" align="center">0.9754&#x2013;1.0235</td>
<td valign="top" align="center">0.9462</td>
<td valign="top" align="center">0.3557</td>
</tr>
<tr>
<td valign="top" align="left">Family <italic>Verrucomicrobiaceae</italic></td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">1.0001</td>
<td valign="top" align="center">0.9672&#x2013;1.0342</td>
<td valign="top" align="center">0.9932</td>
<td valign="top" align="center">0.0736</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Eubacterium ruminantium group</italic></td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">0.9688</td>
<td valign="top" align="center">0.9314&#x2013;1.0078</td>
<td valign="top" align="center">0.1159</td>
<td valign="top" align="center">0.4394</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Akkermansia</italic></td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">1.0004</td>
<td valign="top" align="center">0.9675&#x2013;1.0344</td>
<td valign="top" align="center">0.9829</td>
<td valign="top" align="center">0.0742</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Parabacteroides</italic></td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">0.9927</td>
<td valign="top" align="center">0.9686&#x2013;1.0174</td>
<td valign="top" align="center">0.5584</td>
<td valign="top" align="center">0.9469</td>
</tr>
<tr>
<td valign="top" align="left">Genus <italic>Turicibacter</italic></td>
<td valign="top" align="center">1 &#x00D7; 10<sup>&#x2212;5</sup></td>
<td valign="top" align="center">IVW</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">1.0194</td>
<td valign="top" align="center">0.9846&#x2013;1.0555</td>
<td valign="top" align="center">0.2778</td>
<td valign="top" align="center">0.3992</td>
</tr>
</tbody>
</table></table-wrap>
</sec>
<sec id="S3.SS3">
<title>3.3 Biological annotation</title>
<p>To elucidate the biological mechanisms linking gut microbiota to MM, we annotated the instrumental variables mapping onto 74 genes (<xref ref-type="supplementary-material" rid="TS1">Supplementary Table 8</xref>). A tightly connected network of 74 shared proteins was revealed through PPI network analysis (<xref ref-type="fig" rid="F4">Figure 4A</xref>). To gain further insights, we used the MCC method to rank the top 10 nodes as hub genes in the PPI network (<xref ref-type="fig" rid="F4">Figure 4B</xref> and <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 9</xref>). To uncover the possible molecular mechanisms connecting gut microbiota and MM, we performed gene-set enrichment analysis on these hub genes. KEGG pathway enrichment analysis revealed that hub genes were enriched in pathways related to Th17 cell differentiation and sulfur metabolism (<xref ref-type="fig" rid="F4">Figure 4C</xref> and <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 10</xref>). Additionally, we observed substantial enrichment in processes associated with regulating T cell differentiation and activation through GO analysis (<xref ref-type="fig" rid="F4">Figure 4D</xref> and <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 11</xref>). Significantly, every hub gene displayed connections with numerous characteristics, even when subjected to the stringent significance threshold of 1.05 &#x00D7; 10<sup>&#x2212;5</sup>. Gene-based PheWAS revealed that 9 out of 10 hub genes displayed heightened genetic signals associated with immunological domains (<xref ref-type="fig" rid="F5">Figure 5</xref>, <xref ref-type="supplementary-material" rid="FS1">Supplementary Figure 3</xref> and <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 12</xref>).</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption><p>Biological annotation of 74 mapped genes between MM and gut microbial. <bold>(A)</bold> PPI networks of 74 common genes from the STRING database were visualized in Cytoscape, with disconnected nodes omitted for a concise overview. <bold>(B)</bold> The Maximal Clique Centrality (MCC) method implemented in Cytoscape was utilized to ascertain the top 10 hub genes. <bold>(C)</bold> Enrichment analysis of hub genes using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. <bold>(D)</bold> Gene Ontology (GO) pathways analysis of hub genes. BP, biological process; CC, cellular component; MF, molecular function.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1310444-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption><p>Gene-based phenome-wide association studies. Hub genes displayed heightened genetic signals associated with various domains.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1310444-g005.tif"/>
</fig>
<p>According to the above results, we hypothesize that hub genes exert a substantial influence on the progression of MM. Thus, we utilized the Kaplan-Meier plotter to investigate the association between the expression profiles of hub genes and the overall survival (OS) rates in individuals with MM. We found that higher expression levels of six genes were associated with improved OS in MM patients, whereas high expression of RAB5B was linked to poorer OS (<xref ref-type="fig" rid="F6">Figure 6</xref>). Our GO, KEGG, and PheWAS analyses indicate that hub genes may play a pivotal role in immunological. We leveraged the GSE13591 dataset, which included 133 MM samples and 5 normal donor samples (ND), as a validation set. Gene expression levels underwent normalization through the utilization of the &#x201C;limma&#x201D; package within the R software environment. Immune cell infiltration in MM samples vs. normal donors was assessed using the CIBERSORT method, and the results are detailed in <xref ref-type="supplementary-material" rid="TS1">Supplementary Table 13</xref>. We investigated the correlation between seven hub genes and the abundance of distinct immune cell types within MM samples. The findings were visualized using the &#x201C;vioplot&#x201D; packages. Our analysis revealed a notable positive relationship between AUTS2 and activated Mast cells as well as naive B cells, whereas it demonstrated an inverse correlation with resting Mast cells (<xref ref-type="supplementary-material" rid="FS1">Supplementary Figure 4A</xref>). CDK2 exhibited a positive association with M2 Macrophages and Plasma cells while demonstrating a negative association with T follicular helper cells (<xref ref-type="fig" rid="F7">Figure 7A</xref> and <xref ref-type="supplementary-material" rid="FS1">Supplementary Figure 4B</xref>). The expression of IKZF4 is positively associated with CD8 T cells and inversely associated with M0 Macrophages and CD4 memory resting T cells (<xref ref-type="fig" rid="F7">Figures 7B, C</xref> and <xref ref-type="supplementary-material" rid="FS1">Supplementary Figure 4C</xref>). The occurrence of plasma cells exhibited a favorable correlation with PMEL, while it displayed an adverse link with neutrophils, CD4 memory quiescent T cells, memory B cells, and gamma delta T cells (<xref ref-type="fig" rid="F7">Figures 7D, E</xref> and <xref ref-type="supplementary-material" rid="FS1">Supplementary Figure 4D</xref>). Memory B cells and eosinophils showed a positive correlation with RAB5B, whereas activated mast cells exhibited a negative correlation with RAB5B (<xref ref-type="supplementary-material" rid="FS1">Supplementary Figure 4E</xref>). SUOX exhibited a positive correlation with memory resting T cells CD4 and Neutrophils, as well as memory B cells, while displaying a negative association with naive B cells and M0 Macrophages in our findings (<xref ref-type="fig" rid="F7">Figure 7F</xref> and <xref ref-type="supplementary-material" rid="FS1">Supplementary Figure 4F</xref>). These results suggest that these newly identified hub genes are linked to various immune cells in MM, with a particular emphasis on T cells.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption><p>Kaplan Meier survival curves were utilized to evaluate overall survival (OS) among patients diagnosed with multiple myeloma, while also examining the expression levels of hub genes <bold>(A&#x2013;J)</bold>. Statistically significant findings were determined when all of the following criteria were satisfied: FDR &#x003C; 5%, HR &#x2260; 1, <italic>P</italic> &#x003C; 0.01. FDR, false discovery rate; HR, hazard ratio.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1310444-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption><p>Analysis of immune infiltration in MM using seven hub genes. <bold>(A&#x2013;F)</bold> Exploring the relationships between hub genes and immune cells. Statistical significance is achieved when <italic>p</italic> &#x003C; 0.05.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1310444-g007.tif"/>
</fig>
</sec>
</sec>
<sec id="S4" sec-type="discussion">
<title>4 Discussion</title>
<p>Our study conducted two-sample Mendelian randomization analysis and biological annotation by utilizing the existing cost and statistical power of mbGWAS. We obtained the gut microbiota data from MibioGen and the GWAS summary statistics of MM from GWAS Catalog. Next, we conducted a comprehensive investigation into potential causal associations between microbiota features and MM. Finally, we identified novel hub genes that are associated with immunoregulation and clinical prognosis in MM through biological annotation. The results of our forward MR analyses unveiled compelling evidence of both positive and negative causal effects of ten distinct microbial features on MM risk. Specifically, we observed that family <italic>Acidaminococcaceae</italic>, <italic>Bacteroidales</italic> family S24-7, family <italic>Porphyromonadaceae</italic>, genus <italic>Eubacterium ruminantium group</italic>, genus <italic>Parabacteroides</italic>, and genus <italic>Turicibacter</italic> exhibited positive correlations with MM risk. Conversely, class <italic>Verrucomicrobia</italic>, family <italic>Verrucomicrobiaceae</italic>, genus <italic>Akkermansia</italic> and order <italic>Verrucomicrobiales</italic> demonstrated negative correlations with MM risk. While previous studies have highlighted associations between intestinal microbiota and MM, such as <italic>Klebsiella pneumoniae</italic> (<xref ref-type="bibr" rid="B16">Jian et al., 2020</xref>), <italic>Clostridium butyricum</italic> (<xref ref-type="bibr" rid="B16">Jian et al., 2020</xref>), <italic>Faecalibacterium prausnitzii</italic> (<xref ref-type="bibr" rid="B55">Zhang et al., 2019</xref>), <italic>Pseudomonas aeruginosa</italic> (<xref ref-type="bibr" rid="B55">Zhang et al., 2019</xref>), and <italic>Prevotella heparinolytica</italic> (<xref ref-type="bibr" rid="B8">Calcinotto et al., 2018</xref>). The microbiota revealed in our investigation has not been previously documented in the milieu of MM. <italic>Acidaminococcaceae</italic>, a family in the <italic>Firmicutes</italic> phylum, includes the beer-associated genus <italic>Pectinatus</italic>. Despite being classified as low GC Gram-positive bacteria, <italic>Pectinatus</italic> exhibits unique traits, such as a Gram-negative-like outer membrane and distinctive lipopolysaccharides with remarkable structural heterogeneity (<xref ref-type="bibr" rid="B15">Helander et al., 2004</xref>). IL-22 neutralization was found to lead to a reduction in the abundance of the family <italic>Acidaminococcaceae</italic> (<xref ref-type="bibr" rid="B27">Nagao-Kitamoto et al., 2020</xref>). Without innate immune defense mechanisms, <italic>Porphyromonadaceae</italic> was discovered to thrive, and it exhibited a higher presence within the male microbiota (<xref ref-type="bibr" rid="B11">Fransen et al., 2017</xref>). <italic>Parabacteroides</italic>, a commensal bacterium identified in the gut microbiota, plays a crucial role in mitigating acute pancreatitis, particularly in the context of heparanase-induced exacerbation (<xref ref-type="bibr" rid="B22">Lei et al., 2021</xref>). The administration of <italic>Parabacteroides</italic> has been shown to alleviate acute pancreatitis by producing acetate and reducing neutrophil infiltration in both wild-type and heparanase-transgenic mice. <italic>Turicibacter</italic> has shown the ability to predict both immune-related adverse events and the efficacy of immune checkpoint inhibitors (<xref ref-type="bibr" rid="B13">Hamada et al., 2023</xref>). <italic>Verrucomicrobia</italic>, <italic>Verrucomicrobiales</italic>, and <italic>Verrucomicrobiaceae</italic> are also related to inflammation (<xref ref-type="bibr" rid="B4">Borton et al., 2017</xref>). <xref ref-type="bibr" rid="B25">Liu et al. (2022)</xref> found that TLR4 has a substantial influence on ROR&#x03B3;t++ regulatory T cell reactions. This interaction with <italic>Akkermansia muciniphila</italic> was shown to exert a noteworthy influence on susceptibility to colon inflammation (<xref ref-type="bibr" rid="B25">Liu et al., 2022</xref>). Although our research indicates positive/negative correlations between the aforementioned significant bacteria and MM risk, experimental studies on significant gut microbiota in relation to MM are currently lacking. Therefore, further empirical research is required to substantiate the precise roles and mechanisms of these gut microbiota in the development of MM.</p>
<p>Biological annotation analyses suggest that 74 shared genes are mapped between 10 gut microbiomes and MM. Using the STRING database for PPI network analysis and visualizing it with Cytoscape, we identified top 10 shared genes as hub genes connecting MM and gut microbiota. Enrichment analysis of genes has uncovered connections between hub genes and processes such as T cell differentiation and activation, Th17 cell differentiation, and sulfur metabolism. Additionally, PheWAS analysis revealed a correlation between hub genes and several characteristics, such as metabolic and immunological characteristics, consistent with prior studies (<xref ref-type="bibr" rid="B10">Fan and Pedersen, 2021</xref>). Among these hub genes, seven genes (AUTS2, CDK2, ERBB3, IKZF4, PMEL, RAB5B, SUOX) were identified as having prognostic significance in MM patients. AUTS2 demonstrated significantly higher mutation frequencies among African American cases compared to Caucasians. A different investigative study indicates that Hispanic individuals with AUTS2 mutations among MM patients exhibited poorer overall survival (<xref ref-type="bibr" rid="B34">Peres et al., 2022</xref>). CDK2 plays a carcinogenic role in MM by participating in multiple signaling pathways, and inhibitors targeting CDK2 have the potential to inhibit the growth of MM (<xref ref-type="bibr" rid="B42">Tu et al., 2011</xref>). The ErbB3 atypical expression in MM may have potential effects on tumor cell growth or survival (<xref ref-type="bibr" rid="B44">Walters and Jelinek, 2004</xref>). However, the roles of IKZF4, PMEL, RAB5B, and SUOX in the pathogenesis of MM have not been reported yet. T reg cell effector IKZF4 mediates biological processes such as Nr4a factor maintaining T reg cell lineage stability and inhibitory activity (<xref ref-type="bibr" rid="B39">Sekiya et al., 2015</xref>). PMEL, also known as gp100, participates in cancer immune infiltration and tumor-specific immunotherapy, especially in targeting T cell therapy (<xref ref-type="bibr" rid="B38">Scheffel et al., 2016</xref>). RAB5B, a member of the RAB GTPase family, is pivotal in multiple immune processes, including endocytosis, antigen presentation, and immune signal transduction. SUOX has shown great promise as a diagnostic and prognostic biomarker for various cancers, such as oral squamous cell carcinoma (<xref ref-type="bibr" rid="B28">Nakamura et al., 2018</xref>), and hepatocellular carcinoma (<xref ref-type="bibr" rid="B18">Jin et al., 2013</xref>). Further bioinformatics analysis demonstrated that six hub genes have a strong connection with different immune cells, particularly T cells, in MM samples. For example, CDK2 exhibited a noteworthy inverse correlation with T follicular helper cells, consistent with previous studies (<xref ref-type="bibr" rid="B29">Niu et al., 2023</xref>). IKZF4 exhibited a robust positive correlation with CD8 T cells but displayed a notable inverse relationship with CD4 memory resting T cells. PMEL displayed a significant negative correlation with CD4 memory resting T cells and gamma delta T cells. SUOX exhibited a noteworthy positive association with resting memory CD4 T cells. Therefore, a potential mechanism for the causal link between microbial features and MM may involving the mediation of immune system changes by hub genes, particularly in the context of T cell differentiation and activation.</p>
<p>Despite these valuable insights, it is crucial to recognize several constraints inherent to the current study. To begin with, the examination of bacterial taxa was restricted to the genus category, with no investigation conducted at finer levels like species or strain. Next, although the majority of participants in this GWAS hail from European backgrounds, it&#x2019;s worth noting that the inclusion of a limited number of individuals from diverse ethnicities could potentially impact the outcomes. Consequently, it is important to recognize that our findings may have certain constraints when applied to different racial demographics. Individuals of European ancestry constitute approximately 78% of the total population, while other ancestries make up the remaining 22%. Thus, a noteworthy limitation arises from the disparity in ancestral backgrounds between the data in MiBioGen, which may pose a potential threat to the independence assumption due to population stratification. We chose IVs associated with gut microbiota at a threshold of <italic>p</italic> &#x003C; 1.0 &#x00D7; 10<sup>&#x2212;5</sup>, which exceeds the conventional genome-wide significance level (<italic>p</italic> &#x003C; 5 &#x00D7; 10<sup>&#x2212;8</sup>). This higher threshold was necessary to ensure we had a sufficient number of IVs for our MR analyses, as only one IV remained when selecting IVs at the genome-wide significance level (<italic>p</italic> &#x003C; 5 &#x00D7; 10<sup>&#x2212;8</sup>). In the future, we anticipate achieving more dependable outcomes with MR by significantly expanding the sample size while adhering to stricter criteria. While we have pinpointed several possible causal connections in our MR analyses, it is crucial to emphasize that these findings demand additional confirmation. In our pursuit of a more robust understanding, not only should we contemplate the expansion of our sample size, but we must also direct our attention toward the advancement of novel techniques aimed at bolstering the statistical potency of mbGWAS. Furthermore, it is imperative that we embark on further investigations into the biological mechanisms at play, with a specific focus on elucidating the central genes that act as pivotal links in the intricate interplay between the host genome of patients with MM and their gut microbiome.</p>
<p>In summary, our study has identified several gut microbiota that may possess the potential to impact the risk of MM. We also identified novel host-microbiome shared genes linked to immune regulation and clinical prognosis in MM. These results could hold considerable clinical implications for MM prevention and therapy. Nevertheless, further investigations are needed to elucidate the specific mechanisms that could serve as potential intervention targets for MM.</p>
</sec>
<sec id="S5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The data presented in this study are available at public database, accession link: human multiple myeloma GWAS (data available at: <ext-link ext-link-type="uri" xlink:href="http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90042001-GCST90043000/GCST90042760/">http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST 90042001-GCST90043000/GCST90042760/</ext-link>); human gut microbiome GWAS (data available at: <ext-link ext-link-type="uri" xlink:href="https://mibiogen.gcc.rug.nl">https://mibiogen.gcc.rug.nl</ext-link>); Gene Expression Omnibus (data available at: <ext-link ext-link-type="uri" xlink:href="https://www.ncbi.nlm.nih.gov/geo/">https://www.ncbi. nlm.nih.gov/geo/</ext-link>); Kaplan-Meier plotter database (data available at: <ext-link ext-link-type="uri" xlink:href="http://kmplot.com/analysis/index.php?p=service&#x0026;cancer">http://kmplot.com/analysis/index.php?p=service&#x0026;cancer</ext-link>); STRING database (data available at: <ext-link ext-link-type="uri" xlink:href="https://string-db.org/">https://string-db.org/</ext-link>).</p>
</sec>
<sec id="S6" sec-type="author-contributions">
<title>Author contributions</title>
<p>ZF: Conceptualization, Methodology, Supervision, Writing &#x2013; original draft, Writing &#x2013; review and editing. ML: Conceptualization, Writing &#x2013; review and editing. JB: Data curation, Methodology, Writing &#x2013; review and editing. YL: Data curation, Methodology, Writing &#x2013; review and editing. YC: Investigation, Supervision, Writing &#x2013; review and editing. LiZ: Methodology, Supervision, Validation, Writing &#x2013; review and editing. XG: Investigation, Visualization, Writing &#x2013; review and editing. LL: Conceptualization, Funding acquisition, Supervision, Validation, Writing &#x2013; original draft, Writing &#x2013; review and editing. LSZ: Conceptualization, Funding acquisition, Supervision, Validation, Writing &#x2013; original draft, Writing &#x2013; review and editing.</p>
</sec>
</body>
<back>
<sec id="S7" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (82360029), the National Blood System Disease Clinical Medical Research Center Commissioned Project (2021WWA01), the Construction of Clinical Medical Research Center for Gansu Province Science and Technology Plan Project (21JR7RA435), the Gansu Province Science and Technology Planning Project Natural Science Fund (21JR7RA394), and the Gansu Provincial Science and Technology Plan Project Natural Science Foundation (21JR11RA104).</p>
</sec>
<ack><p>We would like to express our sincere gratitude to the GWAS consortium, Kaplan-Meier Plotter, and GEO database initiative for their dedicated efforts in developing and meticulously curating their valuable data resources for multiple myeloma (MM). Additionally, our heartfelt thanks go out to the MiBioGen consortium for generously providing us with gut microbiota GWAS summary statistics data, which was instrumental in conducting our analyses.</p>
</ack>
<sec id="S8" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="S9" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="S10" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2024.1310444/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmicb.2024.1310444/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="FS1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_1.XLSX" id="TS1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<fn-group>
<fn id="footnote1">
<label>1</label>
<p><ext-link ext-link-type="uri" xlink:href="http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90042001-GCST90043000/GCST90042760/">http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90042001-GCST90043000/GCST90042760/</ext-link></p></fn>
<fn id="footnote2">
<label>2</label>
<p><ext-link ext-link-type="uri" xlink:href="https://mibiogen.gcc.rug.nl">https://mibiogen.gcc.rug.nl</ext-link></p></fn>
<fn id="footnote3">
<label>3</label>
<p><ext-link ext-link-type="uri" xlink:href="https://www.ncbi.nlm.nih.gov/geo/">https://www.ncbi.nlm.nih.gov/geo/</ext-link></p></fn>
<fn id="footnote4">
<label>4</label>
<p><ext-link ext-link-type="uri" xlink:href="http://kmplot.com/analysis/index.php?p=service&#x0026;cancer">http://kmplot.com/analysis/index.php?p=service&#x0026;cancer</ext-link></p></fn>
<fn id="footnote5">
<label>5</label>
<p><ext-link ext-link-type="uri" xlink:href="https://string-db.org/">https://string-db.org/</ext-link></p></fn>
</fn-group>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ahmed</surname> <given-names>N.</given-names></name> <name><surname>Ghannoum</surname> <given-names>M.</given-names></name> <name><surname>Gallogly</surname> <given-names>M.</given-names></name> <name><surname>De Lima</surname> <given-names>M.</given-names></name> <name><surname>Malek</surname> <given-names>E.</given-names></name></person-group> (<year>2020</year>). <article-title>Influence of gut microbiome on multiple myeloma: friend or foe?</article-title> <source><italic>J. Immunother. Cancer</italic></source> <volume>8</volume>:<issue>e000576</issue>.</citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Angtuaco</surname> <given-names>E. J.</given-names></name> <name><surname>Fassas</surname> <given-names>A. B.</given-names></name> <name><surname>Walker</surname> <given-names>R.</given-names></name> <name><surname>Sethi</surname> <given-names>R.</given-names></name> <name><surname>Barlogie</surname> <given-names>B.</given-names></name></person-group> (<year>2004</year>). <article-title>Multiple myeloma: clinical review and diagnostic imaging.</article-title> <source><italic>Radiology</italic></source> <volume>231</volume> <fpage>11</fpage>&#x2013;<lpage>23</lpage>.</citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Badgeley</surname> <given-names>A.</given-names></name> <name><surname>Anwar</surname> <given-names>H.</given-names></name> <name><surname>Modi</surname> <given-names>K.</given-names></name> <name><surname>Murphy</surname> <given-names>P.</given-names></name> <name><surname>Lakshmikuttyamma</surname> <given-names>A.</given-names></name></person-group> (<year>2021</year>). <article-title>Effect of probiotics and gut microbiota on anti-cancer drugs: mechanistic perspectives.</article-title> <source><italic>Biochim. Biophys. Acta Rev. Cancer</italic></source> <volume>1875</volume>:<issue>188494</issue>. <pub-id pub-id-type="doi">10.1016/j.bbcan.2020.188494</pub-id> <pub-id pub-id-type="pmid">33346129</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Borton</surname> <given-names>M. A.</given-names></name> <name><surname>Sabag-Daigle</surname> <given-names>A.</given-names></name> <name><surname>Wu</surname> <given-names>J.</given-names></name> <name><surname>Solden</surname> <given-names>L. M.</given-names></name> <name><surname>O&#x2019;banion</surname> <given-names>B. S.</given-names></name> <name><surname>Daly</surname> <given-names>R. A.</given-names></name></person-group> (<year>2017</year>). <article-title>Chemical and pathogen-induced inflammation disrupt the murine intestinal microbiome.</article-title> <source><italic>Microbiome</italic></source> <volume>5</volume>:<issue>47</issue>. <pub-id pub-id-type="doi">10.1186/s40168-017-0264-8</pub-id> <pub-id pub-id-type="pmid">28449706</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bowden</surname> <given-names>J.</given-names></name> <name><surname>Davey Smith</surname> <given-names>G.</given-names></name> <name><surname>Burgess</surname> <given-names>S.</given-names></name></person-group> (<year>2015</year>). <article-title>Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression.</article-title> <source><italic>Int. J. Epidemiol.</italic></source> <volume>44</volume> <fpage>512</fpage>&#x2013;<lpage>525</lpage>.</citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bowden</surname> <given-names>J.</given-names></name> <name><surname>Davey Smith</surname> <given-names>G.</given-names></name> <name><surname>Haycock</surname> <given-names>P. C.</given-names></name> <name><surname>Burgess</surname> <given-names>S.</given-names></name></person-group> (<year>2016</year>). <article-title>Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator.</article-title> <source><italic>Genet. Epidemiol.</italic></source> <volume>40</volume> <fpage>304</fpage>&#x2013;<lpage>314</lpage>.</citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Caers</surname> <given-names>J.</given-names></name> <name><surname>Garderet</surname> <given-names>L.</given-names></name> <name><surname>Kortum</surname> <given-names>K. M.</given-names></name> <name><surname>O&#x2019;dwyer</surname> <given-names>M. E.</given-names></name> <name><surname>van de Donk</surname> <given-names>N.</given-names></name> <name><surname>Binder</surname> <given-names>M.</given-names></name></person-group> (<year>2018</year>). <article-title>European myeloma network recommendations on tools for the diagnosis and monitoring of multiple myeloma: what to use and when.</article-title> <source><italic>Haematologica</italic></source> <volume>103</volume> <fpage>1772</fpage>&#x2013;<lpage>1784</lpage>.</citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Calcinotto</surname> <given-names>A.</given-names></name> <name><surname>Brevi</surname> <given-names>A.</given-names></name> <name><surname>Chesi</surname> <given-names>M.</given-names></name> <name><surname>Ferrarese</surname> <given-names>R.</given-names></name> <name><surname>Garcia Perez</surname> <given-names>L.</given-names></name> <name><surname>Grioni</surname> <given-names>M.</given-names></name></person-group> (<year>2018</year>). <article-title>Microbiota-driven interleukin-17-producing cells and eosinophils synergize to accelerate multiple myeloma progression.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>9</volume>:<issue>4832</issue>. <pub-id pub-id-type="doi">10.1038/s41467-018-07305-8</pub-id> <pub-id pub-id-type="pmid">30510245</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Corbin</surname> <given-names>L. J.</given-names></name> <name><surname>Richmond</surname> <given-names>R. C.</given-names></name> <name><surname>Wade</surname> <given-names>K. H.</given-names></name> <name><surname>Burgess</surname> <given-names>S.</given-names></name> <name><surname>Bowden</surname> <given-names>J.</given-names></name> <name><surname>Smith</surname> <given-names>G. D.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>BMI as a modifiable risk factor for type 2 diabetes: refining and understanding causal estimates using mendelian randomization.</article-title> <source><italic>Diabetes</italic></source> <volume>65</volume> <fpage>3002</fpage>&#x2013;<lpage>3007</lpage>. <pub-id pub-id-type="doi">10.2337/db16-0418</pub-id> <pub-id pub-id-type="pmid">27402723</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fan</surname> <given-names>Y.</given-names></name> <name><surname>Pedersen</surname> <given-names>O.</given-names></name></person-group> (<year>2021</year>). <article-title>Gut microbiota in human metabolic health and disease.</article-title> <source><italic>Nat. Rev. Microbiol.</italic></source> <volume>19</volume> <fpage>55</fpage>&#x2013;<lpage>71</lpage>.</citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fransen</surname> <given-names>F.</given-names></name> <name><surname>Van Beek</surname> <given-names>A. A.</given-names></name> <name><surname>Borghuis</surname> <given-names>T.</given-names></name> <name><surname>Meijer</surname> <given-names>B.</given-names></name> <name><surname>Hugenholtz</surname> <given-names>F.</given-names></name> <name><surname>van der Gaast-de</surname></name><etal/></person-group> (<year>2017</year>). <article-title>The impact of gut microbiota on gender-specific differences in immunity.</article-title> <source><italic>Front. Immunol.</italic></source> <volume>8</volume>:<issue>754</issue>. <pub-id pub-id-type="doi">10.3389/fimmu.2017.00754</pub-id> <pub-id pub-id-type="pmid">28713378</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Goodrich</surname> <given-names>J. K.</given-names></name> <name><surname>Waters</surname> <given-names>J. L.</given-names></name> <name><surname>Poole</surname> <given-names>A. C.</given-names></name> <name><surname>Sutter</surname> <given-names>J. L.</given-names></name> <name><surname>Koren</surname> <given-names>O.</given-names></name> <name><surname>Blekhman</surname> <given-names>R.</given-names></name></person-group> (<year>2014</year>). <article-title>Human genetics shape the gut microbiome.</article-title> <source><italic>Cell</italic></source> <volume>159</volume> <fpage>789</fpage>&#x2013;<lpage>799</lpage>.</citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hamada</surname> <given-names>K.</given-names></name> <name><surname>Isobe</surname> <given-names>J.</given-names></name> <name><surname>Hattori</surname> <given-names>K.</given-names></name> <name><surname>Hosonuma</surname> <given-names>M.</given-names></name> <name><surname>Baba</surname> <given-names>Y.</given-names></name> <name><surname>Murayama</surname> <given-names>M.</given-names></name></person-group> (<year>2023</year>). <article-title>Turicibacter and Acidaminococcus predict immune-related adverse events and efficacy of immune checkpoint inhibitor.</article-title> <source><italic>Front. Immunol.</italic></source> <volume>14</volume>:<issue>1164724</issue>. <pub-id pub-id-type="doi">10.3389/fimmu.2023.1164724</pub-id> <pub-id pub-id-type="pmid">37207204</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>He</surname> <given-names>H.</given-names></name> <name><surname>Liao</surname> <given-names>S.</given-names></name> <name><surname>Zeng</surname> <given-names>Y.</given-names></name> <name><surname>Liang</surname> <given-names>L.</given-names></name> <name><surname>Chen</surname> <given-names>J.</given-names></name> <name><surname>Tao</surname> <given-names>C.</given-names></name></person-group> (<year>2022</year>). <article-title>Causal relationships between metabolic-associated fatty liver disease and iron status: two-sample Mendelian randomization.</article-title> <source><italic>Liver Int.</italic></source> <volume>42</volume> <fpage>2759</fpage>&#x2013;<lpage>2768</lpage>. <pub-id pub-id-type="doi">10.1111/liv.15455</pub-id> <pub-id pub-id-type="pmid">36226474</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Helander</surname> <given-names>I. M.</given-names></name> <name><surname>Haikara</surname> <given-names>A.</given-names></name> <name><surname>Sadovskaya</surname> <given-names>I.</given-names></name> <name><surname>Vinogradov</surname> <given-names>E.</given-names></name> <name><surname>Salkinoja-Salonen</surname> <given-names>M. S.</given-names></name></person-group> (<year>2004</year>). <article-title>Lipopolysaccharides of anaerobic beer spoilage bacteria of the genus <italic>Pectinatus&#x2013;lipopolysaccharides</italic> of a Gram-positive genus.</article-title> <source><italic>FEMS Microbiol. Rev.</italic></source> <volume>28</volume> <fpage>543</fpage>&#x2013;<lpage>552</lpage>. <pub-id pub-id-type="doi">10.1016/j.femsre.2004.05.001</pub-id> <pub-id pub-id-type="pmid">15539073</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jian</surname> <given-names>X.</given-names></name> <name><surname>Zhu</surname> <given-names>Y.</given-names></name> <name><surname>Ouyang</surname> <given-names>J.</given-names></name> <name><surname>Wang</surname> <given-names>Y.</given-names></name> <name><surname>Lei</surname> <given-names>Q.</given-names></name> <name><surname>Xia</surname> <given-names>J.</given-names></name></person-group> (<year>2020</year>). <article-title>Alterations of gut microbiome accelerate multiple myeloma progression by increasing the relative abundances of nitrogen-recycling bacteria.</article-title> <source><italic>Microbiome</italic></source> <volume>8</volume>:<issue>74</issue>. <pub-id pub-id-type="doi">10.1186/s40168-020-00854-5</pub-id> <pub-id pub-id-type="pmid">32466801</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jiang</surname> <given-names>L.</given-names></name> <name><surname>Zheng</surname> <given-names>Z.</given-names></name> <name><surname>Fang</surname> <given-names>H.</given-names></name> <name><surname>Yang</surname> <given-names>J.</given-names></name></person-group> (<year>2021</year>). <article-title>A generalized linear mixed model association tool for biobank-scale data.</article-title> <source><italic>Nat. Genet.</italic></source> <volume>53</volume> <fpage>1616</fpage>&#x2013;<lpage>1621</lpage>.</citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jin</surname> <given-names>G.-Z.</given-names></name> <name><surname>Yu</surname> <given-names>W.-L.</given-names></name> <name><surname>Dong</surname> <given-names>H.</given-names></name> <name><surname>Zhou</surname> <given-names>W.-P.</given-names></name> <name><surname>Gu</surname> <given-names>Y.-J.</given-names></name> <name><surname>Yu</surname> <given-names>H.</given-names></name></person-group> (<year>2013</year>). <article-title>SUOX is a promising diagnostic and prognostic biomarker for hepatocellular carcinoma.</article-title> <source><italic>J. Hepatol.</italic></source> <volume>59</volume> <fpage>510</fpage>&#x2013;<lpage>517</lpage>. <pub-id pub-id-type="doi">10.1016/j.jhep.2013.04.028</pub-id> <pub-id pub-id-type="pmid">23665285</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Joshua</surname> <given-names>D. E.</given-names></name> <name><surname>Bryant</surname> <given-names>C.</given-names></name> <name><surname>Dix</surname> <given-names>C.</given-names></name> <name><surname>Gibson</surname> <given-names>J.</given-names></name> <name><surname>Ho</surname> <given-names>J.</given-names></name></person-group> (<year>2019</year>). <article-title>Biology and therapy of multiple myeloma.</article-title> <source><italic>Med. J. Australia</italic></source> <volume>210</volume> <fpage>375</fpage>&#x2013;<lpage>380</lpage>.</citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Khan</surname> <given-names>N.</given-names></name> <name><surname>Lindner</surname> <given-names>S.</given-names></name> <name><surname>Gomes</surname> <given-names>A. L. C.</given-names></name> <name><surname>Devlin</surname> <given-names>S. M.</given-names></name> <name><surname>Shah</surname> <given-names>G. L.</given-names></name> <name><surname>Sung</surname> <given-names>A. D.</given-names></name></person-group> (<year>2021</year>). <article-title>Fecal microbiota diversity disruption and clinical outcomes after auto-HCT: a multicenter observational study.</article-title> <source><italic>Blood</italic></source> <volume>137</volume> <fpage>1527</fpage>&#x2013;<lpage>1537</lpage>. <pub-id pub-id-type="doi">10.1182/blood.2020006923</pub-id> <pub-id pub-id-type="pmid">33512409</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Laubach</surname> <given-names>J. P.</given-names></name> <name><surname>Moreau</surname> <given-names>P.</given-names></name> <name><surname>San-Miguel</surname> <given-names>J. F.</given-names></name> <name><surname>Richardson</surname> <given-names>P. G.</given-names></name></person-group> (<year>2015</year>). <article-title>Panobinostat for the treatment of multiple Myeloma.</article-title> <source><italic>Clin. Cancer Res.</italic></source> <volume>21</volume> <fpage>4767</fpage>&#x2013;<lpage>4773</lpage>.</citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lei</surname> <given-names>Y.</given-names></name> <name><surname>Tang</surname> <given-names>L.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name> <name><surname>Hu</surname> <given-names>S.</given-names></name> <name><surname>Wu</surname> <given-names>L.</given-names></name> <name><surname>Liu</surname> <given-names>Y.</given-names></name></person-group> (<year>2021</year>). <article-title>Parabacteroides produces acetate to alleviate heparanase-exacerbated acute pancreatitis through reducing neutrophil infiltration.</article-title> <source><italic>Microbiome</italic></source> <volume>9</volume>:<issue>115</issue>. <pub-id pub-id-type="doi">10.1186/s40168-021-01065-2</pub-id> <pub-id pub-id-type="pmid">34016163</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>X.</given-names></name> <name><surname>Chen</surname> <given-names>Y.</given-names></name> <name><surname>Zhang</surname> <given-names>S.</given-names></name> <name><surname>Dong</surname> <given-names>L.</given-names></name></person-group> (<year>2021</year>). <article-title>Gut microbiota-mediated immunomodulation in tumor.</article-title> <source><italic>J. Exp. Clin. Cancer Res.</italic></source> <volume>40</volume>:<issue>221</issue>.</citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>X.</given-names></name> <name><surname>Qi</surname> <given-names>X.</given-names></name> <name><surname>Han</surname> <given-names>R.</given-names></name> <name><surname>Mao</surname> <given-names>T.</given-names></name> <name><surname>Tian</surname> <given-names>Z.</given-names></name></person-group> (<year>2023</year>). <article-title>Gut microbiota causally affects cholelithiasis: a two-sample Mendelian randomization study.</article-title> <source><italic>Front. Cell Infect Microbiol.</italic></source> <volume>13</volume>:<issue>1253447</issue>. <pub-id pub-id-type="doi">10.3389/fcimb.2023.1253447</pub-id> <pub-id pub-id-type="pmid">37876873</pub-id></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>Y.</given-names></name> <name><surname>Yang</surname> <given-names>M.</given-names></name> <name><surname>Tang</surname> <given-names>L.</given-names></name> <name><surname>Wang</surname> <given-names>F.</given-names></name> <name><surname>Huang</surname> <given-names>S.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name></person-group> (<year>2022</year>). <article-title>TLR4 regulates ROR&#x03B3;t+ regulatory T-cell responses and susceptibility to colon inflammation through interaction with Akkermansia muciniphila.</article-title> <source><italic>Microbiome</italic></source> <volume>10</volume>:<issue>98</issue>.</citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lopera-Maya</surname> <given-names>E. A.</given-names></name> <name><surname>Kurilshikov</surname> <given-names>A.</given-names></name> <name><surname>van der Graaf</surname> <given-names>A.</given-names></name> <name><surname>Hu</surname> <given-names>S.</given-names></name> <name><surname>Andreu-Sanchez</surname> <given-names>S.</given-names></name> <name><surname>Chen</surname> <given-names>L.</given-names></name></person-group> (<year>2022</year>). <article-title>Effect of host genetics on the gut microbiome in 7,738 participants of the Dutch Microbiome project.</article-title> <source><italic>Nat. Genet.</italic></source> <volume>54</volume> <fpage>143</fpage>&#x2013;<lpage>151</lpage>.</citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nagao-Kitamoto</surname> <given-names>H.</given-names></name> <name><surname>Leslie</surname> <given-names>J. L.</given-names></name> <name><surname>Kitamoto</surname> <given-names>S.</given-names></name> <name><surname>Jin</surname> <given-names>C.</given-names></name> <name><surname>Thomsson</surname> <given-names>K. A.</given-names></name> <name><surname>Gillilland</surname> <given-names>M. G.</given-names></name></person-group> (<year>2020</year>). <article-title>Interleukin-22-mediated host glycosylation prevents Clostridioides difficile infection by modulating the metabolic activity of the gut microbiota.</article-title> <source><italic>Nat. Med.</italic></source> <volume>26</volume> <fpage>608</fpage>&#x2013;<lpage>617</lpage>. <pub-id pub-id-type="doi">10.1038/s41591-020-0764-0</pub-id> <pub-id pub-id-type="pmid">32066975</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nakamura</surname> <given-names>K.</given-names></name> <name><surname>Akiba</surname> <given-names>J.</given-names></name> <name><surname>Ogasawara</surname> <given-names>S.</given-names></name> <name><surname>Naito</surname> <given-names>Y.</given-names></name> <name><surname>Nakayama</surname> <given-names>M.</given-names></name> <name><surname>Abe</surname> <given-names>Y.</given-names></name></person-group> (<year>2018</year>). <article-title>SUOX is negatively associated with multistep carcinogenesis and proliferation in oral squamous cell carcinoma.</article-title> <source><italic>Med. Mol. Morphol.</italic></source> <volume>51</volume> <fpage>102</fpage>&#x2013;<lpage>110</lpage>. <pub-id pub-id-type="doi">10.1007/s00795-017-0177-4</pub-id> <pub-id pub-id-type="pmid">29280012</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Niu</surname> <given-names>X.</given-names></name> <name><surname>Liu</surname> <given-names>Z.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Wu</surname> <given-names>D.</given-names></name></person-group> (<year>2023</year>). <article-title>Green tea EGCG inhibits na&#x00EF;ve CD4+ T cell division and progression in mice: an integration of network pharmacology, molecular docking and experimental validation.</article-title> <source><italic>Curr. Res. Food Sci.</italic></source> <volume>7</volume>:<issue>100537</issue>. <pub-id pub-id-type="doi">10.1016/j.crfs.2023.100537</pub-id> <pub-id pub-id-type="pmid">37441168</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nolde</surname> <given-names>M.</given-names></name> <name><surname>Holtfreter</surname> <given-names>B.</given-names></name> <name><surname>Kocher</surname> <given-names>T.</given-names></name> <name><surname>Alayash</surname> <given-names>Z.</given-names></name> <name><surname>Reckelkamm</surname> <given-names>S. L.</given-names></name> <name><surname>Ehmke</surname> <given-names>B.</given-names></name></person-group> (<year>2022</year>). <article-title>No bidirectional relationship between depression and periodontitis: a genetic correlation and Mendelian randomization study.</article-title> <source><italic>Front. Immunol.</italic></source> <volume>13</volume>:<issue>918404</issue>.</citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Owens</surname> <given-names>B.</given-names></name></person-group> (<year>2020</year>). <article-title>Outlook: multiple myeloma.</article-title> <source><italic>Nature</italic></source> <volume>587</volume>:<issue>S55</issue>.</citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Papiol</surname> <given-names>S.</given-names></name> <name><surname>Schmitt</surname> <given-names>A.</given-names></name> <name><surname>Maurus</surname> <given-names>I.</given-names></name> <name><surname>Rossner</surname> <given-names>M. J.</given-names></name> <name><surname>Schulze</surname> <given-names>T. G.</given-names></name> <name><surname>Falkai</surname> <given-names>P.</given-names></name></person-group> (<year>2021</year>). <article-title>Association between physical activity and schizophrenia: results of a 2-sample mendelian randomization analysis.</article-title> <source><italic>JAMA Psychiatry</italic></source> <volume>78</volume> <fpage>441</fpage>&#x2013;<lpage>444</lpage>.</citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Park</surname> <given-names>E. M.</given-names></name> <name><surname>Chelvanambi</surname> <given-names>M.</given-names></name> <name><surname>Bhutiani</surname> <given-names>N.</given-names></name> <name><surname>Kroemer</surname> <given-names>G.</given-names></name> <name><surname>Zitvogel</surname> <given-names>L.</given-names></name> <name><surname>Wargo</surname> <given-names>J. A.</given-names></name></person-group> (<year>2022</year>). <article-title>Targeting the gut and tumor microbiota in cancer.</article-title> <source><italic>Nat. Med.</italic></source> <volume>28</volume> <fpage>690</fpage>&#x2013;<lpage>703</lpage>.</citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peres</surname> <given-names>L. C.</given-names></name> <name><surname>Colin-Leitzinger</surname> <given-names>C. M.</given-names></name> <name><surname>Teng</surname> <given-names>M.</given-names></name> <name><surname>Dutil</surname> <given-names>J.</given-names></name> <name><surname>Alugubelli</surname> <given-names>R. R.</given-names></name> <name><surname>Deavila</surname> <given-names>G.</given-names></name></person-group> (<year>2022</year>). <article-title>Racial and ethnic differences in clonal hematopoiesis, tumor markers, and outcomes of patients with multiple myeloma.</article-title> <source><italic>Blood Adv.</italic></source> <volume>6</volume> <fpage>3767</fpage>&#x2013;<lpage>3778</lpage>. <pub-id pub-id-type="doi">10.1182/bloodadvances.2021006652</pub-id> <pub-id pub-id-type="pmid">35500227</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pianko</surname> <given-names>M. J.</given-names></name> <name><surname>Golob</surname> <given-names>J. L.</given-names></name></person-group> (<year>2022</year>). <article-title>Host-microbe interactions and outcomes in multiple myeloma and hematopoietic stem cell transplantation.</article-title> <source><italic>Cancer Metastasis Rev.</italic></source> <volume>41</volume> <fpage>367</fpage>&#x2013;<lpage>382</lpage>. <pub-id pub-id-type="doi">10.1007/s10555-022-10033-7</pub-id> <pub-id pub-id-type="pmid">35488106</pub-id></citation></ref>
<ref id="B36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Psaltopoulou</surname> <given-names>T.</given-names></name> <name><surname>Sergentanis</surname> <given-names>T. N.</given-names></name> <name><surname>Kanellias</surname> <given-names>N.</given-names></name> <name><surname>Kanavidis</surname> <given-names>P.</given-names></name> <name><surname>Terpos</surname> <given-names>E.</given-names></name> <name><surname>Dimopoulos</surname> <given-names>M. A.</given-names></name></person-group> (<year>2013</year>). <article-title>Tobacco smoking and risk of multiple myeloma: a meta-analysis of 40 observational studies.</article-title> <source><italic>Int. J. Cancer</italic></source> <volume>132</volume> <fpage>2413</fpage>&#x2013;<lpage>2431</lpage>. <pub-id pub-id-type="doi">10.1002/ijc.27898</pub-id> <pub-id pub-id-type="pmid">23047159</pub-id></citation></ref>
<ref id="B37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rodriguez-Garcia</surname> <given-names>A.</given-names></name> <name><surname>Arroyo</surname> <given-names>A.</given-names></name> <name><surname>Garcia-Vicente</surname> <given-names>R.</given-names></name> <name><surname>Morales</surname> <given-names>M. L.</given-names></name> <name><surname>Gomez-Gordo</surname> <given-names>R.</given-names></name> <name><surname>Justo</surname> <given-names>P.</given-names></name></person-group> (<year>2023</year>). <article-title>Short-chain fatty acid production by gut microbiota predicts treatment response in multiple myeloma.</article-title> <source><italic>Clin. Cancer Res.</italic></source> Online ahead of print. <pub-id pub-id-type="doi">10.1158/1078-0432.CCR-23-0195</pub-id> <pub-id pub-id-type="pmid">38109212</pub-id></citation></ref>
<ref id="B38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Scheffel</surname> <given-names>M. J.</given-names></name> <name><surname>Scurti</surname> <given-names>G.</given-names></name> <name><surname>Simms</surname> <given-names>P.</given-names></name> <name><surname>Garrett-Mayer</surname> <given-names>E.</given-names></name> <name><surname>Mehrotra</surname> <given-names>S.</given-names></name> <name><surname>Nishimura</surname> <given-names>M. I.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>Efficacy of adoptive T-cell therapy is improved by treatment with the antioxidant N-Acetyl cysteine, which limits activation-induced T-cell death.</article-title> <source><italic>Cancer Res.</italic></source> <volume>76</volume> <fpage>6006</fpage>&#x2013;<lpage>6016</lpage>. <pub-id pub-id-type="doi">10.1158/0008-5472.CAN-16-0587</pub-id> <pub-id pub-id-type="pmid">27742673</pub-id></citation></ref>
<ref id="B39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sekiya</surname> <given-names>T.</given-names></name> <name><surname>Kondo</surname> <given-names>T.</given-names></name> <name><surname>Shichita</surname> <given-names>T.</given-names></name> <name><surname>Morita</surname> <given-names>R.</given-names></name> <name><surname>Ichinose</surname> <given-names>H.</given-names></name> <name><surname>Yoshimura</surname> <given-names>A.</given-names></name></person-group> (<year>2015</year>). <article-title>Suppression of Th2 and Tfh immune reactions by Nr4a receptors in mature T reg cells.</article-title> <source><italic>J. Exp. Med.</italic></source> <volume>212</volume> <fpage>1623</fpage>&#x2013;<lpage>1640</lpage>. <pub-id pub-id-type="doi">10.1084/jem.20142088</pub-id> <pub-id pub-id-type="pmid">26304965</pub-id></citation></ref>
<ref id="B40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shi</surname> <given-names>H.</given-names></name> <name><surname>Zhao</surname> <given-names>T.</given-names></name> <name><surname>Geng</surname> <given-names>R.</given-names></name> <name><surname>Sun</surname> <given-names>L.</given-names></name> <name><surname>Fan</surname> <given-names>H.</given-names></name></person-group> (<year>2023</year>). <article-title>The associations between gut microbiota and chronic respiratory diseases: a Mendelian randomization study.</article-title> <source><italic>Front. Microbiol.</italic></source> <volume>14</volume>:<issue>1200937</issue>. <pub-id pub-id-type="doi">10.3389/fmicb.2023.1200937</pub-id> <pub-id pub-id-type="pmid">37333634</pub-id></citation></ref>
<ref id="B41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>J.</given-names></name> <name><surname>Wu</surname> <given-names>Y.</given-names></name> <name><surname>Yin</surname> <given-names>X.</given-names></name> <name><surname>Ma</surname> <given-names>H.</given-names></name> <name><surname>Zhang</surname> <given-names>J.</given-names></name></person-group> (<year>2023</year>). <article-title>The causal links between gut microbiota and COVID-19: a Mendelian randomization study.</article-title> <source><italic>J. Med. Virol.</italic></source> <volume>95</volume>:<issue>e28784</issue>. <pub-id pub-id-type="doi">10.1002/jmv.28784</pub-id> <pub-id pub-id-type="pmid">37219044</pub-id></citation></ref>
<ref id="B42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tu</surname> <given-names>Y.-S.</given-names></name> <name><surname>Kang</surname> <given-names>X.-L.</given-names></name> <name><surname>Zhou</surname> <given-names>J.-G.</given-names></name> <name><surname>Lv</surname> <given-names>X.-F.</given-names></name> <name><surname>Tang</surname> <given-names>Y.-B.</given-names></name> <name><surname>Guan</surname> <given-names>Y.-Y.</given-names></name></person-group> (<year>2011</year>). <article-title>Involvement of Chk1-Cdc25A-cyclin A/CDK2 pathway in simvastatin induced S-phase cell cycle arrest and apoptosis in multiple myeloma cells.</article-title> <source><italic>Eur. J. Pharmacol.</italic></source> <volume>670</volume> <fpage>356</fpage>&#x2013;<lpage>364</lpage>. <pub-id pub-id-type="doi">10.1016/j.ejphar.2011.09.031</pub-id> <pub-id pub-id-type="pmid">21958871</pub-id></citation></ref>
<ref id="B43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Verbanck</surname> <given-names>M.</given-names></name> <name><surname>Chen</surname> <given-names>C.-Y.</given-names></name> <name><surname>Neale</surname> <given-names>B.</given-names></name> <name><surname>Do</surname> <given-names>R.</given-names></name></person-group> (<year>2018</year>). <article-title>Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases.</article-title> <source><italic>Nat. Genet.</italic></source> <volume>50</volume> <fpage>693</fpage>&#x2013;<lpage>698</lpage>.</citation></ref>
<ref id="B44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Walters</surname> <given-names>D. K.</given-names></name> <name><surname>Jelinek</surname> <given-names>D. F.</given-names></name></person-group> (<year>2004</year>). <article-title>A role for Janus kinases in crosstalk between ErbB3 and the interferon-alpha signaling complex in myeloma cells.</article-title> <source><italic>Oncogene</italic></source> <volume>23</volume> <fpage>1197</fpage>&#x2013;<lpage>1205</lpage>. <pub-id pub-id-type="doi">10.1038/sj.onc.1207203</pub-id> <pub-id pub-id-type="pmid">14647450</pub-id></citation></ref>
<ref id="B45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wastyk</surname> <given-names>H. C.</given-names></name> <name><surname>Fragiadakis</surname> <given-names>G. K.</given-names></name> <name><surname>Perelman</surname> <given-names>D.</given-names></name> <name><surname>Dahan</surname> <given-names>D.</given-names></name> <name><surname>Merrill</surname> <given-names>B. D.</given-names></name> <name><surname>Yu</surname> <given-names>F. B.</given-names></name></person-group> (<year>2021</year>). <article-title>Gut-microbiota-targeted diets modulate human immune status.</article-title> <source><italic>Cell</italic></source> <volume>184</volume> <fpage>4137</fpage>&#x2013;<lpage>4153e14</lpage>.</citation></ref>
<ref id="B46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Watanabe</surname> <given-names>K.</given-names></name> <name><surname>Stringer</surname> <given-names>S.</given-names></name> <name><surname>Frei</surname> <given-names>O.</given-names></name> <name><surname>Umi&#x0107;evi&#x0107; Mirkov</surname> <given-names>M.</given-names></name> <name><surname>de Leeuw</surname> <given-names>C.</given-names></name> <name><surname>Polderman</surname> <given-names>T. J. C.</given-names></name></person-group> (<year>2019</year>). <article-title>A global overview of pleiotropy and genetic architecture in complex traits.</article-title> <source><italic>Nat. Genet.</italic></source> <volume>51</volume> <fpage>1339</fpage>&#x2013;<lpage>1348</lpage>.</citation></ref>
<ref id="B47"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Watanabe</surname> <given-names>K.</given-names></name> <name><surname>Taskesen</surname> <given-names>E.</given-names></name> <name><surname>van Bochoven</surname> <given-names>A.</given-names></name> <name><surname>Posthuma</surname> <given-names>D.</given-names></name></person-group> (<year>2017</year>). <article-title>Functional mapping and annotation of genetic associations with FUMA.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>8</volume>:<issue>1826</issue>.</citation></ref>
<ref id="B48"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Welter</surname> <given-names>D.</given-names></name> <name><surname>Macarthur</surname> <given-names>J.</given-names></name> <name><surname>Morales</surname> <given-names>J.</given-names></name> <name><surname>Burdett</surname> <given-names>T.</given-names></name> <name><surname>Hall</surname> <given-names>P.</given-names></name> <name><surname>Junkins</surname> <given-names>H.</given-names></name></person-group> (<year>2014</year>). <article-title>The NHGRI GWAS Catalog, a curated resource of SNP-trait associations.</article-title> <source><italic>Nucleic Acids Res.</italic></source> <volume>42</volume> <fpage>D1001</fpage>&#x2013;<lpage>D1006</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gkt1229</pub-id> <pub-id pub-id-type="pmid">24316577</pub-id></citation></ref>
<ref id="B49"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xiao</surname> <given-names>L.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name> <name><surname>Wu</surname> <given-names>Y.</given-names></name> <name><surname>Huang</surname> <given-names>Y.</given-names></name> <name><surname>Tao</surname> <given-names>S.</given-names></name> <name><surname>Liu</surname> <given-names>Y.</given-names></name></person-group> (<year>2023</year>). <article-title>The interactions between host genome and gut microbiome increase the risk of psychiatric disorders: mendelian randomization and biological annotation.</article-title> <source><italic>Brain Behav. Immun.</italic></source> <volume>113</volume> <fpage>389</fpage>&#x2013;<lpage>400</lpage>.</citation></ref>
<ref id="B50"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yavorska</surname> <given-names>O. O.</given-names></name> <name><surname>Burgess</surname> <given-names>S.</given-names></name></person-group> (<year>2017</year>). <article-title>MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data.</article-title> <source><italic>Int. J. Epidemiol.</italic></source> <volume>46</volume> <fpage>1734</fpage>&#x2013;<lpage>1739</lpage>.</citation></ref>
<ref id="B51"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yonekura</surname> <given-names>S.</given-names></name> <name><surname>Terrisse</surname> <given-names>S.</given-names></name> <name><surname>Alves Costa, Silva</surname> <given-names>C.</given-names></name> <name><surname>Lafarge</surname> <given-names>A.</given-names></name> <name><surname>Iebba</surname> <given-names>V.</given-names></name><etal/></person-group> (<year>2022</year>). <article-title>Cancer induces a stress ileopathy depending on beta-adrenergic receptors and promoting dysbiosis that contributes to carcinogenesis.</article-title> <source><italic>Cancer Discov.</italic></source> <volume>12</volume> <fpage>1128</fpage>&#x2013;<lpage>1151</lpage>. <pub-id pub-id-type="doi">10.1158/2159-8290.CD-21-0999</pub-id> <pub-id pub-id-type="pmid">34930787</pub-id></citation></ref>
<ref id="B52"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yuan</surname> <given-names>S.</given-names></name> <name><surname>Larsson</surname> <given-names>S. C.</given-names></name></person-group> (<year>2022</year>). <article-title>Coffee and caffeine consumption and risk of kidney stones: a mendelian randomization study.</article-title> <source><italic>Am. J. Kidney Dis.</italic></source> <volume>79</volume> <fpage>9</fpage>&#x2013;<lpage>14e1</lpage>. <pub-id pub-id-type="doi">10.1053/j.ajkd.2021.04.018</pub-id> <pub-id pub-id-type="pmid">34690004</pub-id></citation></ref>
<ref id="B53"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zeng</surname> <given-names>C.</given-names></name> <name><surname>Zhang</surname> <given-names>C.</given-names></name> <name><surname>He</surname> <given-names>C.</given-names></name> <name><surname>Song</surname> <given-names>H.</given-names></name></person-group> (<year>2023</year>). <article-title>Investigating the causal impact of gut microbiota on glioblastoma: a bidirectional Mendelian randomization study.</article-title> <source><italic>BMC Genomics</italic></source> <volume>24</volume>:<issue>784</issue>. <pub-id pub-id-type="doi">10.1186/s12864-023-09885-2</pub-id> <pub-id pub-id-type="pmid">38110895</pub-id></citation></ref>
<ref id="B54"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zeng</surname> <given-names>Y.</given-names></name> <name><surname>Cao</surname> <given-names>S.</given-names></name> <name><surname>Yang</surname> <given-names>H.</given-names></name></person-group> (<year>2023</year>). <article-title>Roles of gut microbiome in epilepsy risk: a Mendelian randomization study.</article-title> <source><italic>Front. Microbiol.</italic></source> <volume>14</volume>:<issue>1115014</issue>. <pub-id pub-id-type="doi">10.3389/fmicb.2023.1115014</pub-id> <pub-id pub-id-type="pmid">36922970</pub-id></citation></ref>
<ref id="B55"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>B.</given-names></name> <name><surname>Gu</surname> <given-names>J.</given-names></name> <name><surname>Liu</surname> <given-names>J.</given-names></name> <name><surname>Huang</surname> <given-names>B.</given-names></name> <name><surname>Li</surname> <given-names>J.</given-names></name></person-group> (<year>2019</year>). <article-title>Fecal microbiota taxonomic shifts in chinese multiple myeloma patients analyzed by Quantitative Polimerase Chain Reaction (QPCR) and 16S rRNA high-throughput sequencing.</article-title> <source><italic>Med. Sci. Monit.</italic></source> <volume>25</volume> <fpage>8269</fpage>&#x2013;<lpage>8280</lpage>. <pub-id pub-id-type="doi">10.12659/MSM.919988</pub-id> <pub-id pub-id-type="pmid">31678982</pub-id></citation></ref>
<ref id="B56"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>L.</given-names></name> <name><surname>Xiang</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Zhang</surname> <given-names>J.</given-names></name></person-group> (<year>2022</year>). <article-title>Gut microbiome in multiple myeloma: mechanisms of progression and clinical applications.</article-title> <source><italic>Front. Immunol.</italic></source> <volume>13</volume>:<issue>1058272</issue>. <pub-id pub-id-type="doi">10.3389/fimmu.2022.1058272</pub-id> <pub-id pub-id-type="pmid">36569873</pub-id></citation></ref>
<ref id="B57"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname> <given-names>C. B.</given-names></name> <name><surname>Zhou</surname> <given-names>Y. L.</given-names></name> <name><surname>Fang</surname> <given-names>J. Y.</given-names></name></person-group> (<year>2021</year>). <article-title>Gut microbiota in cancer immune response and immunotherapy.</article-title> <source><italic>Trends Cancer</italic></source> <volume>7</volume> <fpage>647</fpage>&#x2013;<lpage>660</lpage>.</citation></ref>
<ref id="B58"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname> <given-names>Y.</given-names></name> <name><surname>Jian</surname> <given-names>X.</given-names></name> <name><surname>Chen</surname> <given-names>S.</given-names></name> <name><surname>An</surname> <given-names>G.</given-names></name> <name><surname>Jiang</surname> <given-names>D.</given-names></name> <name><surname>Yang</surname> <given-names>Q.</given-names></name></person-group> (<year>2024</year>). <article-title>Targeting gut microbial nitrogen recycling and cellular uptake of ammonium to improve bortezomib resistance in multiple myeloma.</article-title> <source><italic>Cell Metab.</italic></source> <volume>36</volume> <fpage>159</fpage>&#x2013;<lpage>175e8</lpage>. <pub-id pub-id-type="doi">10.1016/j.cmet.2023.11.019</pub-id> <pub-id pub-id-type="pmid">38113887</pub-id></citation></ref>
</ref-list>
</back>
</article>
