<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Immunol.</journal-id>
<journal-title>Frontiers in Immunology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Immunol.</abbrev-journal-title>
<issn pub-type="epub">1664-3224</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fimmu.2021.748447</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Immunology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Interaction Between Chronic Endometritis Caused Endometrial Microbiota Disorder and Endometrial Immune Environment Change in Recurrent Implantation Failure</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Peigen</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/977012"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Panyu</given-names>
</name>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Guo</surname>
<given-names>Yingchun</given-names>
</name>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Fang</surname>
<given-names>Cong</given-names>
</name>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Li</surname>
<given-names>Tingting</given-names>
</name>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<institution>Reproductive Medicine Center, The Sixth Affiliated Hospital, Sun Yat-sen University</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Juarez Antonio Sim&#xf5;es Quaresma, Evandro Chagas Institute, Brazil</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Gulrukhsor Haybulloevna Tolibova, D. O. Ott Research Institute of Obstetrics, Gynecology and Reproductology, Russia; Kotaro Kitaya, Reproduction Clinic Osaka, Japan</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Tingting Li, <email xlink:href="mailto:litt33@mail.sysu.edu.cn">litt33@mail.sysu.edu.cn</email>; Cong Fang, <email xlink:href="mailto:fangcong@mail.sysu.edu.cn">fangcong@mail.sysu.edu.cn</email>
</p>
</fn>
<fn fn-type="other" id="fn002">
<p>This article was submitted to Microbial Immunology, a section of the journal Frontiers in Immunology</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>04</day>
<month>10</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>12</volume>
<elocation-id>748447</elocation-id>
<history>
<date date-type="received">
<day>28</day>
<month>07</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>13</day>
<month>09</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2021 Chen, Chen, Guo, Fang and Li</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Chen, Chen, Guo, Fang and Li</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>Objective</title>
<p>To investigate the Interaction between chronic endometritis (CE) caused endometrial microbiota disorder and endometrial immune environment change in recurrent implantation failure (RIF).</p>
</sec>
<sec>
<title>Method</title>
<p>Transcriptome sequencing analysis of the endometrial of 112 patients was preform by using High-Throughput Sequencing. The endometrial microbiota of 43 patients was analyzed by using 16s rRNA sequencing technology.</p>
</sec>
<sec>
<title>Result</title>
<p>In host endometrium, CD4 T cell and macrophage exhibited significant differences abundance between CE and non-CE patients. The enrichment analysis indicated differentially expressed genes mainly enriched in immune-related functional terms. <italic>Phyllobacterium</italic> and <italic>Sphingomonas</italic> were significantly high infiltration in CE patients, and active in pathways related to carbohydrate metabolism and/or fat metabolism. The increased synthesis of lipopolysaccharide, an important immunomodulator, was the result of microbial disorders in the endometrium.</p>
</sec>
<sec>
<title>Conclusion</title>
<p>The composition of endometrial microorganisms in CE and non-CE patients were significantly different. <italic>Phyllobacterium</italic> and <italic>Sphingomonas</italic> mainly regulated immune cells by interfering with the process of carbohydrate metabolism and/or fat metabolism in the endometrium. CE endometrial microorganisms might regulate Th17 response and the ratio of Th1 to Th17 through lipopolysaccharide (LPS).</p>
</sec>
</abstract>
<kwd-group>
<kwd>chronic endometritis</kwd>
<kwd>endometrial microbiota</kwd>
<kwd>immune cells</kwd>
<kwd>recurrent implantation failure</kwd>
<kwd>host-microbiome association</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="44"/>
<page-count count="8"/>
<word-count count="3030"/>
</counts>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<title>Introduction</title>
<p>Recurrent implantation failure (RIF) is usually defined as a woman under 40 who has received at least 4 high-quality embryos in at least three fresh or frozen cycles or has transferred more than 10 high-quality embryos and still cannot obtain a clinical pregnancy (<xref ref-type="bibr" rid="B1">1</xref>&#x2013;<xref ref-type="bibr" rid="B3">3</xref>). Recent studies had confirmed that chronic endometritis (CE) was an important cause of recurrent reproductive failure (included recurrent miscarriage and recurrent implantation failure) (<xref ref-type="bibr" rid="B4">4</xref>). Chronic endometritis (CE) was a persistent inflammation of the endometrium, and CD138 immunohistochemical (IHC) staining of plasma cells was a more accurate and sensitive diagnostic method for the diagnosis of CE (<xref ref-type="bibr" rid="B5">5</xref>&#x2013;<xref ref-type="bibr" rid="B9">9</xref>).</p>
<p>Multiple studies had shown that a good balance of innate and adaptive local and peripheral immune systems was the basis for a successful pregnancy (<xref ref-type="bibr" rid="B10">10</xref>, <xref ref-type="bibr" rid="B11">11</xref>). The study by Yuye Li et&#xa0;al. showed that the abundance of endometrial immune cells in CE patients was increased, and the high abundance of immune cells might be related to decreased endometrial receptivity and repeated pregnancy failures (<xref ref-type="bibr" rid="B4">4</xref>).</p>
<p>It was now clear that microorganisms affect immunity (<xref ref-type="bibr" rid="B12">12</xref>, <xref ref-type="bibr" rid="B13">13</xref>). It also knew that disturbances in the balance between the microbiota and the immune system might lead to inappropriate immune responses or excessive inflammation, or down-regulation of the immune response and pathogenic bacteria&#x2019;s dominance of normal symbionts &#x201c;ecological disorders&#x201d; (<xref ref-type="bibr" rid="B14">14</xref>). However, the mechanism of action of endometrial microbiota in RIF patients with CE remains unclear. The purpose of this study was to explore the regulation of endometrial microbes on the endometrial immune cells of RIF patients with CE.</p>
</sec>
<sec id="s2">
<title>Method</title>
<sec id="s2_1">
<title>Collection of Research Subjects and Ethical Approval</title>
<p>The RIF cohort selected patients who underwent frozen embryo transfer (FET) in Reproductive Medicine Center, The Sixth Affiliated Hospital, Sun Yat-sen University, with the inclusion criteria. (a). At least four good-quality embryos were accepted in at least three fresh or frozen cycles, but still not pregnant. (b). Normal ovarian reserve function (levels of follicle stimulation hormone (FSH) &lt;12 mIU/ml and AMH &gt;1.1&#x2009;ng/ml). (c). Confirmed by the examination (ultrasound or hysteroscopy) confirms that the uterine cavity was in normal shape. (d). Normal chromosome karyotype. (e). No antibiotics were received in the month before collection. (f). All patients were examined for cervical and vaginal secretions one month and one week before the embryo transplant to rule out infectious diseases.</p>
<p>The Ethics Committee approved the study of Sixth Affiliated Hospital of Sun Yat-Sen University. All patients were fully informed and voluntarily participate in the study and sign an informed consent form. At the same time, this study did not cover any research products.</p>
</sec>
<sec id="s2_2">
<title>Sample Collection Procession</title>
<p>Endometrial specimens were collected on LH+7 (natural cycles) or P+5 (artificial hormone cycle) according to the process (described in <xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Data Sheet 1</bold>
</xref>) and stored at -20&#xb0;C.</p>
</sec>
<sec id="s2_3">
<title>Host RNA Extraction and Sequencing</title>
<p>After completing the library construction according to the manufacturer&#x2019;s instructions (as described in <xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Data Sheet 1</bold>
</xref>), pair-end sequencing was performed on the HiSeq 2500 platform (Illumina).</p>
</sec>
<sec id="s2_4">
<title>Host RNA-Seq Analysis</title>
<p>FASTQ files were processed by fastp (<xref ref-type="bibr" rid="B15">15</xref>), Bowite2 (<xref ref-type="bibr" rid="B16">16</xref>), HISAT2 (<xref ref-type="bibr" rid="B17">17</xref>) and StringTie (<xref ref-type="bibr" rid="B18">18</xref>, <xref ref-type="bibr" rid="B19">19</xref>) to get the count matrix (as described in <xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Data Sheet 1</bold>
</xref>). The human genome GRCh38 from Gencode v26 was used as the reference genome. Then the count value of gene expression was normalized to TPM (Transcripts Per Kilobase of exon model per Million mapped reads) value.</p>
</sec>
<sec id="s2_5">
<title>Selection of Differentially Expressed Genes and Functional Enrichment Analysis</title>
<p>The differentially expressed genes between CE and non-CE patients were selected by using DESeq2 (<xref ref-type="bibr" rid="B20">20</xref>) R package with the p value &lt; 0.05 and log fold change &gt;1 as cut-off. Then functional enrichment analysis of DEGs was performed by clusterprofiler R package (<xref ref-type="bibr" rid="B21">21</xref>) and OmicShare tool (<uri xlink:href="http://www.omicshare.com/tools">http://www.omicshare.com/tools</uri>).</p>
</sec>
<sec id="s2_6">
<title>Calculation of Endometrial Immune Cells Infiltration</title>
<p>CIBERSORT is a tool that can estimate the abundances of 22 types of immune cells, including seven T cell types, na&#xef;ve and memory B cells, plasma cells, NK cells, and myeloid subsets (<xref ref-type="bibr" rid="B22">22</xref>). By using CIBERSORT, the infiltration of immune cells in endometrium was calculated. To understand the consistence of T cells, we also used ImmuCellAI web tools which included 18 T-cell subtypes (<xref ref-type="bibr" rid="B23">23</xref>) to calculate the infiltration of subtypes of T-cells.</p>
</sec>
<sec id="s2_7">
<title>16S rRNA Extraction and Sequencing</title>
<p>After completing the library construction according to the manufacturer&#x2019;s instructions with the universal primers 341F (5&#x2032;-CCTACGGGNGGCWGCAG-3&#x2032;)and 806R(5&#x2032;-GGACTACHVGGGTWTCTAAT-3&#x2032;)(as described in <xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Data Sheet 1</bold>
</xref>), NovaSeq 6000 was used for high throughput sequencing.</p>
</sec>
<sec id="s2_8">
<title>Endometrial Microbiota Analysis</title>
<p>Endometrial microbiota analysis was processed following QIIME2 analysis pipeline (as described in <xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Data Sheet 1</bold>
</xref>) (<xref ref-type="bibr" rid="B24">24</xref>).</p>
<p>The differential abundance of endometrial microbiota was selected using LEfSe (Linear discriminant analysis Effect Size) (<xref ref-type="bibr" rid="B25">25</xref>).</p>
<p>PICRUSt2 (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) (<xref ref-type="bibr" rid="B26">26</xref>) was used to predict the functional profiles of the endometrial microbes. Different metacyc (<xref ref-type="bibr" rid="B27">27</xref>) functional profiles were screened by STAMP software (version 2.1.3) (<xref ref-type="bibr" rid="B28">28</xref>).</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<title>Results</title>
<sec id="s3_1">
<title>Endometrial Transcriptome Analysis</title>
<p>A total of 112 subjects (included 32 CE patients and 72 non-CE patients) were recruited in this study. After obtaining the gene expression matrix, we obtained 93 differentially expressed genes through DESeq2 (logFC &#x2265; |1|, <italic>p</italic> &lt; 0,05 [<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>, <xref ref-type="supplementary-material" rid="SF3">
<bold>Supplementary Table 1</bold>
</xref>)]. By functional enrichment analysis, we found that DEGs were mainly enriched in the immune and inflammatory-related GO (Gene Ontology) terms, such as immune response (GO:0006955), immune system process (GO:0002376), and inflammatory response (GO: 0006954) (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1B</bold>
</xref>, <xref ref-type="supplementary-material" rid="SF4">
<bold>Supplementary Table 2</bold>
</xref>). The top 20 enriched pathways were showed in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1D</bold>
</xref> and the interaction network of these pathways was described in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1C</bold>
</xref>. Most of these pathways were also related to inflammation and immune processes in the endometrium.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Functional enrichment analysis of DEGs. <bold>(A)</bold> Volcano map for differential expression analysis; <bold>(B)</bold> Z-score bubble chart of Gene Ontology (GO) enrichment analysis; <bold>(C)</bold> Network diagram of KEGG enrichment analysis; <bold>(D)</bold> Z-score bubble chart of KEGG enrichment analysis.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-12-748447-g001.tif"/>
</fig>
</sec>
<sec id="s3_2">
<title>Calculation of Endometrial Immune Cells Infiltration</title>
<p>To further explore the infiltration of immune cells in the endometrium, we used COBERSORT to calculate the degree of infiltration of 22 immune cells (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>). The infiltration of memory CD4 T cells and Macrophages was a significant difference between the two groups (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>). We also used ImmuCellAI to estimate the infiltration of 18 subtypes of T cells (<xref ref-type="supplementary-material" rid="SF5">
<bold>Supplementary Table 3</bold>
</xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Calculation of endometrial immune cells infiltration. <bold>(A)</bold> The composition of 22 immune cells in the CE group and the non-CE group; <bold>(B)</bold> Comparison of differences in immune cell immersion.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-12-748447-g002.tif"/>
</fig>
</sec>
<sec id="s3_3">
<title>Endometrial Microbiota Analysis</title>
<p>We performed 16s rRNA sequencing on the endometrial microbiota of 43 patients (14 CE patients and 29 non-CE patients). As the rarefaction curve was shown (<xref ref-type="supplementary-material" rid="SF1">
<bold>Supplementary Figure&#xa0;1</bold>
</xref>), the 16s rRNA sequencing in this study had sufficient sequencing depth. The top 15 phylum and top 20 Genus were shown in <xref ref-type="fig" rid="f3">
<bold>Figures&#xa0;3A, B</bold>
</xref>. The &#x3b1; diversity of the endometrium microbiota calculated by Simpson showed that there was no significant difference in species diversity between the two groups of feature level (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>), phylum level (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3D</bold>
</xref>), and genus level (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3E</bold>
</xref>) (T-test).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>The structural characteristics of the endometrial microbiota. Taxonomic classification of the endometrial microbiota of CE (P) and non-CE (N) group at the level of Phylum <bold>(A)</bold> and Genus <bold>(B)</bold>. The &#x3b1; diversity calculated by Simpson index of the endometrium microbiota between CE and non-CE group at the level of feature <bold>(C)</bold>, Phylum <bold>(D)</bold> and Genus <bold>(E)</bold>. <bold>(F)</bold> Linear discriminant analysis of the differential abundance endometrial microbiota between CE and non-CE group.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-12-748447-g003.tif"/>
</fig>
<p>Based on the results of LEfSe (LDA &#x2265; 2), the abundance of <italic>Phyllobacterium</italic> and <italic>Sphingomonas</italic> was significantly increased in CE patients (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3F</bold>
</xref>). Spearman correlation analysis was performed by using the cor. test() function with a two-sided alternative hypothesis to analyze the relationship between the immune cells and genus selected by LEfSe. Compared with other methods (such as Pearson), the normalized count (gene expression) and component data (relative abundance of microbiota) of Spearman&#x2019;s correlation analysis performed better (<xref ref-type="bibr" rid="B29">29</xref>). Cytoscape software (<xref ref-type="bibr" rid="B30">30</xref>) was used to visualize the results of the analysis.</p>
<p>Based on the Metacyc database, PICRUST2 was used to analyze the metabolic function of the endometrial microbiota. Through the STAMP software, we found that the CE group was mainly enriched in the sucrose biosynthesis III pathway (PWY-7347) and sucrose biosynthesis I pathway (SUCSYN-PWY) (Welch&#x2019;s t-test, <italic>p</italic> &lt; 0.05, <xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4A</bold>
</xref> and <xref ref-type="supplementary-material" rid="SF6">
<bold>Supplementary Table 4</bold>
</xref>).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Interaction between endometrial immune cell and microbiota. <bold>(A)</bold> Comparison of differences in metabolic pathway between CE and non-CE group; <bold>(B)</bold> Network visualizing the relationship between genus, immune cell marker genes and immune cells. The red line indicates the positive correlation (red) or negative (blue) relationship between the genera and the gene.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-12-748447-g004.tif"/>
</fig>
</sec>
<sec id="s3_4">
<title>Interaction Between Endometrial Immune Cell and Microbiota</title>
<p>The marker genes of immune cells were obtained from a previous study (<xref ref-type="bibr" rid="B23">23</xref>). The correlation network between the selected Genus and the marker genes of immune cells was constructed by Spearman analysis (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4B</bold>
</xref>). We found that <italic>Phyllobacterium</italic> and <italic>Sphingomonas</italic> had similar relationships with immune cells, and they were mainly positively related to dendritic cells (DC), natural killer (NK) cells, iTreg cells and B cells. These two genera were negatively correlated to macrophage cells.</p>
<p>In the host&#x2019;s pathway analysis, T cell-related signaling pathways show an important position. After calculating the abundance of 18 T cell subtypes using ImmuCellAI, we analyzed the related genus of several subtypes reported in other studies (<xref ref-type="bibr" rid="B14">14</xref>) included Th1, Th17, and T regular (Treg) (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref> and <xref ref-type="supplementary-material" rid="SF7">
<bold>Supplementary Table 5</bold>
</xref>). Endometrial microbiota in patients with high Th1 abundance (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5B</bold>
</xref>) was active in several glycolysis-related pathways, including super pathway of thiamine diphosphate biosynthesis I (THISYN-PWY), reductive TCA cycle I (P23-PWY), super pathway of L-aspartate and L-asparagine biosynthesis (ASPASN-PWY) and purine nucleobases degradation I (P164-PWY) (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5A</bold>
</xref>). The PWY-7332 pathway, which was highly active in the endometrium of patients with high abundance of Th17, was a super pathway for the synthesis of Lipopolysaccharides (LPS), which is an important immunomodulator (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5C, E</bold>
</xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Interaction between T cell subtypes and microbiota. Comparison of differences in metabolic pathway of Th1 <bold>(A)</bold>, Th17 <bold>(E)</bold>, and Treg <bold>(F)</bold>. Different genera related to Th1 <bold>(B)</bold>, Th17 <bold>(C)</bold>, Treg <bold>(D)</bold> abundance.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-12-748447-g005.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>Discussion</title>
<p>Chronic endometritis (CE) is an important factor for recurrent implantation failure (<xref ref-type="bibr" rid="B31">31</xref>). The endometrium of CE patients showed a high abundance of immune cell infiltration, which would cause a decrease in the receptivity of the endometrium, an important feature of recurrent implantation failure (<xref ref-type="bibr" rid="B4">4</xref>). It has been proven that immune adaptation, one of the key mechanisms for the establishment and maintenance of pregnancy, was directly affected by local microorganisms (<xref ref-type="bibr" rid="B32">32</xref>). In this study, we combined the endometrial transcriptome and microbial diversity analysis to illustrate that endometrial microbiota regulated the abundance of endometrial immune cells through metabolic activities.</p>
<p>The functional enrichment analysis of differential genes in the endometrial transcriptome showed that the functions of DEGs were mainly enriched in immune and inflammation-related items. In the network diagram of enriched pathways, cytokine-cytokine receptor interaction and Toll-like receptor signaling pathway were the central pathways in the network. Toll-like receptors (TLRs) are the first line of defense against pathogen invasion and play a key role in inflammation, immune cell regulation, survival, and proliferation (<xref ref-type="bibr" rid="B33">33</xref>). This is the result of the body&#x2019;s defense during endometritis. At the same time, this pathway and cytokine-cytokine receptor interaction played an important role in the regulation of T cell plasticity and differentiation (<xref ref-type="bibr" rid="B33">33</xref>, <xref ref-type="bibr" rid="B34">34</xref>). After calculating the infiltration of endometrial immune cells by CIBERSORT, we found that resting CD4 memory T cells was high infiltration in the non-CE group. This might indicate that in CE, CD4 T cells are significantly activated due to inflammation. This was consistent with the conclusion of another study (<xref ref-type="bibr" rid="B4">4</xref>).</p>
<p>It has been proven that bacteria can affect immunity (<xref ref-type="bibr" rid="B12">12</xref>, <xref ref-type="bibr" rid="B13">13</xref>). We introduced 16s rRNA sequencing technology to analyze the impact of endometrial microbiota on immune cells. The results of LEfSe showed that the abundance of <italic>Phyllobacterium</italic> and <italic>Sphingomonas</italic> in the CE group was significantly higher than that in the non-CE group. In the network diagram of the relationship between different genus and immune cell marker genes, we found that <italic>Phyllobacterium</italic> and <italic>Sphingomonas</italic> had a significant positive correlation with B cells. Since the infiltration of plasma cells in the endometrial stroma was considered a characteristic of endometritis (<xref ref-type="bibr" rid="B35">35</xref>), it confirmed that these two different genera played an important role in endometritis.</p>
<p>At the same time, we noticed a significant positive correlation between <italic>Sphingomonas</italic> and uterus NK (uNK) cells. In addition, in the differential analysis of the transcriptome, the expression of CD16 (Fc&#x3b3;RIII) in the CE group was significantly higher than that in the non-CE group. CD16+ is a characteristic of the uNK cell population that specifically kills infected cells (<xref ref-type="bibr" rid="B36">36</xref>). The abnormal differentiation of uNK cells might lead to embryo implantation failure through the development of endometrial and/or decidual blood vessels and trophoblasts (<xref ref-type="bibr" rid="B37">37</xref>). The study of May-Tal Sauerbrun-Cutler et&#xa0;al. showed that a higher level of CD16 + uNK cells was a risk factor for embryo transfer failure in infertile women before frozen embryo transfer (<xref ref-type="bibr" rid="B38">38</xref>). Macrophage was known as a common anti-epidemic presenting cell (APC). It was the main cytokine producer in the human endometrium (<xref ref-type="bibr" rid="B39">39</xref>) and played an important role in the establishment of endometrial receptivity (<xref ref-type="bibr" rid="B40">40</xref>). The results of this study suggested that <italic>Phyllobacterium</italic> and <italic>Sphingomonas</italic> were significantly negatively correlated with macrophages. This might be one of the important mechanisms of endometrial microbiota affecting endometrial receptivity.</p>
<p>It was known that metabolic activity had an important regulatory effect on immune cells (<xref ref-type="bibr" rid="B41">41</xref>, <xref ref-type="bibr" rid="B42">42</xref>). Endometrial microbiota was metabolically active, and the crosstalk of metabolic activity might be an important way for endometrial microbiota to interact with endometrial immune cells. In the study, two metabolic pathways were highly active in the CE group, PWY-7347, and SUCSYN-PWY. SUCSYN-PWY described the process of generating sucrose by <italic>phosphoglycerol</italic> and <italic>3-phospho-D-glyceric acid</italic>. The intermediate products included <italic>UDP-&#x3b1;-D-glucose</italic>. PWY-7347 described the process of <italic>glucose-6-phosphate</italic> and <italic>UDP-&#x3b1;-D-glucose</italic> to produce sucrose, and phosphate and UDP were produced in the middle. Obviously, these two metabolic pathways were interaction with glucose and lipid metabolism, which had been proved to be closely related to the regulation of T cells and macrophages (<xref ref-type="bibr" rid="B41">41</xref>, <xref ref-type="bibr" rid="B42">42</xref>). These active microbiotas competed for energy materials of the host (including endometrial cells and immune cells) by consuming intermediate products of glycolysis.</p>
<p>Th1 is a subtype of CD4 T cells, and its increased abundance was considered to be unfavorable for embryo implantation (<xref ref-type="bibr" rid="B43">43</xref>). In this study, we found that endometrial microorganisms with high Th1 abundance exhibited significant-high activity in many important metabolic pathways, including L-aspartic acid and L-asparagine biosynthesis, reducing TCA cycle, degradation of purine nucleobases, and biosynthesis of thiamine diphosphate. Obviously, this was related to multiple links in carbohydrate metabolism (including glycolysis). Therefore, the abnormal activities of the microbiota in the endometrium might regulate the conversion of Th1/Th2 by interfering with the process of carbohydrate metabolism and/or fat metabolism. In addition, Th17 and Treg also played an important role in the regulation of immune adaptation during embryo implantation, and the imbalance of Th17/Treg ratio would have an adverse effect on embryo implantation. In the study of intestinal microbes, the researchers confirmed that lipopolysaccharide (LPS) could induce a decrease in the number of Tregs and an increase in the number of Th17 and Th1 (<xref ref-type="bibr" rid="B44">44</xref>). In our study, the endometrial microbiota of patients with high Th17 abundance was significantly active in the relevant pathways of LPS synthesis. It meant that the imbalance in the number and proportion of immune cells in the endometrium may be caused by a similar mechanism.</p>
</sec>
<sec id="s5">
<title>Conclusion</title>
<p>In this study, we combined the transcriptome and 16s rRNA sequencing technology to analyze the interaction between the endometrium microbial disorder caused by chronic endometritis and the immune cells in the endometrium of patients with recurrent implantation failure. We found that the composition of endometrial microorganisms in CE and non-CE patients were significantly different. <italic>Phyllobacterium</italic> and <italic>Sphingomonas</italic> mainly regulated immune cells by interfering with the process of carbohydrate metabolism and/or fat metabolism in the endometrium. The results of this study supplemented part of the mechanism of microbial regulation of the immune environment of the endometrium.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data Availability Statement</title>
<p>The transcriptome sequencing for 112 endometrium samples have been deposited with the National Center for Biotechnology Information (NCBI) under reference number PRJNA747622. The 16s rRNA gene sequencing for 43 endometrial microbiota samples have been deposited with the NCBI under reference number PRJNA732058.</p>
</sec>
<sec id="s7" sec-type="ethics-statement">
<title>Ethics Statement</title>
<p>The studies involving human participants were reviewed and approved by Sixth Affiliated Hospital of Sun Yat-Sen University. The patients/participants provided their written informed consent to participate in this study.</p>
</sec>
<sec id="s8" sec-type="author-contributions">
<title>Author Contributions</title>
<p>CF, TL, and PeC carried out the study. PeC and TL analyzed and interpreted the data and drafted the manuscript. YG and PaC collected and followed up the samples. TL and CF revised the manuscript. CF, TL, and PeC coordinated the study, participated in the design, and reviewed the manuscript. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="s9" sec-type="funding-information">
<title>Funding</title>
<p>This work was supported by the National Natural Science Foundation of China (grant no. 81871214), the National Key R&amp;D Program of China (grant no. 2017YFC1001603), and the National Natural Science Foundation of China (grant no. 81070495).</p>
</sec>
<sec id="s10" 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="s11" 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>
</body>
<back>
<ack>
<title>Acknowledgments</title>
<p>We thank Xukang Medical&#xb7;Yikon Genomics for their assistance in the sequencing process.</p>
</ack>
<sec id="s12" 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/fimmu.2021.748447/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fimmu.2021.748447/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Image_1.tif" id="SF1" mimetype="image/tiff">
<label>Supplementary Figure&#xa0;1</label>
<caption>
<p>The rarefaction curve of 16s rRNA data.</p>
</caption>
</supplementary-material>
  <supplementary-material xlink:href="DataSheet_1.docx" id="SF2" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document">
<label>Supplement Data Sheet 1</label>
<caption>
<p>Detailed description of Methods.</p>
</caption>
</supplementary-material>
  <supplementary-material xlink:href="Table_1.xlsx" id="SF3" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet">
<label>Supplementary Table&#xa0;1</label>
<caption>
<p>Differentially expressed genes of DESeq2.</p>
</caption>
</supplementary-material>
  <supplementary-material xlink:href="Table_2.xlsx" id="SF4" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet">
<label>Supplementary Table&#xa0;2</label>
<caption>
<p>The result table of functional enrichment analysis of DEGs.</p>
</caption>
</supplementary-material>
  <supplementary-material xlink:href="Table_3.xlsx" id="SF5" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet">
<label>Supplementary Table&#xa0;3</label>
<caption>
<p>Infiltration of 18 subtypes of T cells.</p>
</caption>
</supplementary-material>
  <supplementary-material xlink:href="Table_4.xlsx" id="SF6" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet">
<label>Supplementary Table&#xa0;4</label>
<caption>
<p>Comparison Results Table of PICRUST2 Metabolic&#xa0;Pathway.</p>
</caption>
</supplementary-material>
  <supplementary-material xlink:href="Table_5.xlsx" id="SF7" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet">
<label>Supplementary Table&#xa0;5</label>
<caption>
<p>Comparison Results Table of PICRUST2 Metabolic Pathway of Th1, Th17, and Treg.</p>
</caption>
</supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Coughlan</surname> <given-names>C</given-names>
</name>
<name>
<surname>Ledger</surname> <given-names>W</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Q</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>F</given-names>
</name>
<name>
<surname>Demirol</surname> <given-names>A</given-names>
</name>
<name>
<surname>Gurgan</surname> <given-names>T</given-names>
</name>
<etal/>
</person-group>. <article-title>Recurrent Implantation Failure: Definition and Management</article-title>. <source>Reprod BioMed Online</source> (<year>2014</year>) <volume>28</volume>(<issue>1</issue>):<fpage>14</fpage>&#x2013;<lpage>38</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.rbmo.2013.08.011</pub-id>
</citation>
</ref>
<ref id="B2">
<label>2</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Simon</surname> <given-names>A</given-names>
</name>
<name>
<surname>Laufer</surname> <given-names>N</given-names>
</name>
</person-group>. <article-title>Repeated Implantation Failure: Clinical Approach</article-title>. <source>Fertil Steril</source> (<year>2012</year>) <volume>97</volume>(<issue>5</issue>):<page-range>1039&#x2013;43</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.fertnstert.2012.03.010</pub-id>
</citation>
</ref>
<ref id="B3">
<label>3</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thornhill</surname> <given-names>AR</given-names>
</name>
<name>
<surname>deDie-Smulders</surname> <given-names>CE</given-names>
</name>
<name>
<surname>Geraedts</surname> <given-names>JP</given-names>
</name>
<name>
<surname>Harper</surname> <given-names>JC</given-names>
</name>
<name>
<surname>Harton</surname> <given-names>GL</given-names>
</name>
<name>
<surname>Lavery</surname> <given-names>SA</given-names>
</name>
<etal/>
</person-group>. <article-title>ESHRE PGD Consortium 'Best Practice Guidelines for Clinical Preimplantation Genetic Diagnosis (PGD) and Preimplantation Genetic Screening (PGS)'</article-title>. <source>Hum Reprod</source> (<year>2005</year>) <volume>20</volume>(<issue>1</issue>):<fpage>35</fpage>&#x2013;<lpage>48</lpage>. doi: <pub-id pub-id-type="doi">10.1093/humrep/deh579</pub-id>
</citation>
</ref>
<ref id="B4">
<label>4</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>S</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>C</given-names>
</name>
<name>
<surname>Lian</surname> <given-names>R</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>C</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>Evaluation of Peripheral and Uterine Immune Status of Chronic Endometritis in Patients With Recurrent Reproductive Failure</article-title>. <source>Fertil Steril</source> (<year>2020</year>) <volume>113</volume>(<issue>1</issue>):<fpage>187</fpage>&#x2013;<lpage>96.e1</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fertnstert.2019.09.001</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kitaya</surname> <given-names>K</given-names>
</name>
<name>
<surname>Yasuo</surname> <given-names>T</given-names>
</name>
</person-group>. <article-title>Inter-Observer and Intra-Observer Variability in Immunohistochemical Detection of Endometrial Stromal Plasmacytes in Chronic Endometritis</article-title>. <source>Exp Ther Med</source> (<year>2013</year>) <volume>5</volume>(<issue>2</issue>):<page-range>485&#x2013;8</page-range>. doi: <pub-id pub-id-type="doi">10.3892/etm.2012.824</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bayer-Garner</surname> <given-names>IB</given-names>
</name>
<name>
<surname>Korourian</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Plasma Cells in Chronic Endometritis are Easily Identified When Stained With Syndecan-1</article-title>. <source>Mod Pathol</source> (<year>2001</year>) <volume>14</volume>(<issue>9</issue>):<page-range>877&#x2013;9</page-range>. doi: <pub-id pub-id-type="doi">10.1038/modpathol.3880405</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bayer-Garner</surname> <given-names>IB</given-names>
</name>
<name>
<surname>Nickell</surname> <given-names>JA</given-names>
</name>
<name>
<surname>Korourian</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Routine Syndecan-1 Immunohistochemistry Aids in the Diagnosis of Chronic Endometritis</article-title>. <source>Arch Pathol Lab Med</source> (<year>2004</year>) <volume>128</volume>(<issue>9</issue>):<page-range>1000&#x2013;3</page-range>. doi: <pub-id pub-id-type="doi">10.5858/2004-128-1000-RSIAIT</pub-id>
</citation>
</ref>
<ref id="B8">
<label>8</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kitaya</surname> <given-names>K</given-names>
</name>
<name>
<surname>Yasuo</surname> <given-names>T</given-names>
</name>
</person-group>. <article-title>Aberrant Expression of Selectin E, CXCL1, and CXCL13 in Chronic Endometritis</article-title>. <source>Mod Pathol</source> (<year>2010</year>) <volume>23</volume>(<issue>8</issue>):<page-range>1136&#x2013;46</page-range>. doi: <pub-id pub-id-type="doi">10.1038/modpathol.2010.98</pub-id>
</citation>
</ref>
<ref id="B9">
<label>9</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Disep</surname> <given-names>B</given-names>
</name>
<name>
<surname>Innes</surname> <given-names>BA</given-names>
</name>
<name>
<surname>Cochrane</surname> <given-names>HR</given-names>
</name>
<name>
<surname>Tijani</surname> <given-names>S</given-names>
</name>
<name>
<surname>Bulmer</surname> <given-names>JN</given-names>
</name>
</person-group>. <article-title>Immunohistochemical Characterization of Endometrial Leucocytes in Endometritis</article-title>. <source>Histopathology</source> (<year>2004</year>) <volume>45</volume>(<issue>6</issue>):<page-range>625&#x2013;32</page-range>. doi: <pub-id pub-id-type="doi">10.1111/j.1365-2559.2004.02052.x</pub-id>
</citation>
</ref>
<ref id="B10">
<label>10</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aghaeepour</surname> <given-names>N</given-names>
</name>
<name>
<surname>Ganio</surname> <given-names>EA</given-names>
</name>
<name>
<surname>McIlwain</surname> <given-names>D</given-names>
</name>
<name>
<surname>Tsai</surname> <given-names>AS</given-names>
</name>
<name>
<surname>Tingle</surname> <given-names>M</given-names>
</name>
<name>
<surname>Van Gassen</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>An Immune Clock of Human Pregnancy</article-title>. <source>Sci Immunol</source> (<year>2017</year>) <volume>2</volume>(<issue>15</issue>):<fpage>eaan2946</fpage>. doi: <pub-id pub-id-type="doi">10.1126/sciimmunol.aan2946</pub-id>
</citation>
</ref>
<ref id="B11">
<label>11</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>PrabhuDas</surname> <given-names>M</given-names>
</name>
<name>
<surname>Bonney</surname> <given-names>E</given-names>
</name>
<name>
<surname>Caron</surname> <given-names>K</given-names>
</name>
<name>
<surname>Dey</surname> <given-names>S</given-names>
</name>
<name>
<surname>Erlebacher</surname> <given-names>A</given-names>
</name>
<name>
<surname>Fazleabas</surname> <given-names>A</given-names>
</name>
<etal/>
</person-group>. <article-title>Immune Mechanisms at the Maternal-Fetal Interface: Perspectives and Challenges</article-title>. <source>Nat Immunol</source> (<year>2015</year>) <volume>16</volume>(<issue>4</issue>):<page-range>328&#x2013;34</page-range>. doi: <pub-id pub-id-type="doi">10.1038/ni.3131</pub-id>
</citation>
</ref>
<ref id="B12">
<label>12</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hooper</surname> <given-names>LV</given-names>
</name>
<name>
<surname>Littman</surname> <given-names>DR</given-names>
</name>
<name>
<surname>Macpherson</surname> <given-names>AJ</given-names>
</name>
</person-group>. <article-title>Interactions Between the Microbiota and the Immune System</article-title>. <source>Science</source> (<year>2012</year>) <volume>336</volume>(<issue>6086</issue>):<page-range>1268&#x2013;73</page-range>. doi: <pub-id pub-id-type="doi">10.1126/science.1223490</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Littman</surname> <given-names>DR</given-names>
</name>
<name>
<surname>Pamer</surname> <given-names>EG</given-names>
</name>
</person-group>. <article-title>Role of the Commensal Microbiota in Normal and Pathogenic Host Immune Responses</article-title>. <source>Cell Host Microbe</source> (<year>2011</year>) <volume>10</volume>(<issue>4</issue>):<page-range>311&#x2013;23</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.chom.2011.10.004</pub-id>
</citation>
</ref>
<ref id="B14">
<label>14</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Al-Nasiry</surname> <given-names>S</given-names>
</name>
<name>
<surname>Ambrosino</surname> <given-names>E</given-names>
</name>
<name>
<surname>Schlaepfer</surname> <given-names>M</given-names>
</name>
<name>
<surname>Morre</surname> <given-names>SA</given-names>
</name>
<name>
<surname>Wieten</surname> <given-names>L</given-names>
</name>
<name>
<surname>Voncken</surname> <given-names>JW</given-names>
</name>
<etal/>
</person-group>. <article-title>The Interplay Between Reproductive Tract Microbiota and Immunological System in Human Reproduction</article-title>. <source>Front Immunol</source> (<year>2020</year>) <volume>11</volume>:<elocation-id>378</elocation-id>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2020.00378</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>S</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Gu</surname> <given-names>J</given-names>
</name>
</person-group>. <article-title>Fastp: An Ultra-Fast All-in-One FASTQ Preprocessor</article-title>. <source>Bioinformatics</source> (<year>2018</year>) <volume>34</volume>(<issue>17</issue>):<page-range>i884&#x2013;i90</page-range>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/bty560</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Langmead</surname> <given-names>B</given-names>
</name>
<name>
<surname>Salzberg</surname> <given-names>SL</given-names>
</name>
</person-group>. <article-title>Fast Gapped-Read Alignment With Bowtie 2</article-title>. <source>Nat Methods</source> (<year>2012</year>) <volume>9</volume>(<issue>4</issue>):<page-range>357&#x2013;9</page-range>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.1923</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname> <given-names>D</given-names>
</name>
<name>
<surname>Langmead</surname> <given-names>B</given-names>
</name>
<name>
<surname>Salzberg</surname> <given-names>SL</given-names>
</name>
</person-group>. <article-title>HISAT: A Fast Spliced Aligner With Low Memory Requirements</article-title>. <source>Nat Methods</source> (<year>2015</year>) <volume>12</volume>(<issue>4</issue>):<page-range>357&#x2013;60</page-range>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.3317</pub-id>
</citation>
</ref>
<ref id="B18">
<label>18</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pertea</surname> <given-names>M</given-names>
</name>
<name>
<surname>Pertea</surname> <given-names>GM</given-names>
</name>
<name>
<surname>Antonescu</surname> <given-names>CM</given-names>
</name>
<name>
<surname>Chang</surname> <given-names>TC</given-names>
</name>
<name>
<surname>Mendell</surname> <given-names>JT</given-names>
</name>
<name>
<surname>Salzberg</surname> <given-names>SL</given-names>
</name>
</person-group>. <article-title>StringTie Enables Improved Reconstruction of a Transcriptome From RNA-Seq Reads</article-title>. <source>Nat Biotechnol</source> (<year>2015</year>) <volume>33</volume>(<issue>3</issue>):<page-range>290&#x2013;5</page-range>. doi: <pub-id pub-id-type="doi">10.1038/nbt.3122</pub-id>
</citation>
</ref>
<ref id="B19">
<label>19</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pertea</surname> <given-names>M</given-names>
</name>
<name>
<surname>Kim</surname> <given-names>D</given-names>
</name>
<name>
<surname>Pertea</surname> <given-names>GM</given-names>
</name>
<name>
<surname>Leek</surname> <given-names>JT</given-names>
</name>
<name>
<surname>Salzberg</surname> <given-names>SL</given-names>
</name>
</person-group>. <article-title>Transcript-Level Expression Analysis of RNA-Seq Experiments With HISAT, StringTie and Ballgown</article-title>. <source>Nat Protoc</source> (<year>2016</year>) <volume>11</volume>(<issue>9</issue>):<page-range>1650&#x2013;67</page-range>. doi: <pub-id pub-id-type="doi">10.1038/nprot.2016.095</pub-id>
</citation>
</ref>
<ref id="B20">
<label>20</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Love</surname> <given-names>MI</given-names>
</name>
<name>
<surname>Huber</surname> <given-names>W</given-names>
</name>
<name>
<surname>Anders</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2</article-title>. <source>Genome Biol</source> (<year>2014</year>) <volume>15</volume>(<issue>12</issue>):<fpage>550</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s13059-014-0550-8</pub-id>
</citation>
</ref>
<ref id="B21">
<label>21</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yu</surname> <given-names>G</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>LG</given-names>
</name>
<name>
<surname>Han</surname> <given-names>Y</given-names>
</name>
<name>
<surname>He</surname> <given-names>QY</given-names>
</name>
</person-group>. <article-title>Clusterprofiler: An R Package for Comparing Biological Themes Among Gene Clusters</article-title>. <source>OMICS</source> (<year>2012</year>) <volume>16</volume>(<issue>5</issue>):<page-range>284&#x2013;7</page-range>. doi: <pub-id pub-id-type="doi">10.1089/omi.2011.0118</pub-id>
</citation>
</ref>
<ref id="B22">
<label>22</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Newman</surname> <given-names>AM</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>CL</given-names>
</name>
<name>
<surname>Green</surname> <given-names>MR</given-names>
</name>
<name>
<surname>Gentles</surname> <given-names>AJ</given-names>
</name>
<name>
<surname>Feng</surname> <given-names>W</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>Y</given-names>
</name>
<etal/>
</person-group>. <article-title>Robust Enumeration of Cell Subsets From Tissue Expression Profiles</article-title>. <source>Nat Methods</source> (<year>2015</year>) <volume>12</volume>(<issue>5</issue>):<page-range>453&#x2013;7</page-range>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.3337</pub-id>
</citation>
</ref>
<ref id="B23">
<label>23</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Miao</surname> <given-names>YR</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Q</given-names>
</name>
<name>
<surname>Lei</surname> <given-names>Q</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>M</given-names>
</name>
<name>
<surname>Xie</surname> <given-names>GY</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H</given-names>
</name>
<etal/>
</person-group>. <article-title>ImmuCellAI: A Unique Method for Comprehensive T-Cell Subsets Abundance Prediction and its Application in Cancer Immunotherapy</article-title>. <source>Adv Sci (Weinh)</source> (<year>2020</year>) <volume>7</volume>(<issue>7</issue>):<fpage>1902880</fpage>. doi: <pub-id pub-id-type="doi">10.1002/advs.201902880</pub-id>
</citation>
</ref>
<ref id="B24">
<label>24</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bolyen</surname> <given-names>E</given-names>
</name>
<name>
<surname>Rideout</surname> <given-names>JR</given-names>
</name>
<name>
<surname>Dillon</surname> <given-names>MR</given-names>
</name>
<name>
<surname>Bokulich</surname> <given-names>NA</given-names>
</name>
<name>
<surname>Abnet</surname> <given-names>CC</given-names>
</name>
<name>
<surname>Al-Ghalith</surname> <given-names>GA</given-names>
</name>
<etal/>
</person-group>. <article-title>Reproducible, Interactive, Scalable and Extensible Microbiome Data Science Using QIIME 2</article-title>. <source>Nat Biotechnol</source> (<year>2019</year>) <volume>37</volume>(<issue>8</issue>):<page-range>852&#x2013;7</page-range>. doi: <pub-id pub-id-type="doi">10.1038/s41587-019-0209-9</pub-id>
</citation>
</ref>
<ref id="B25">
<label>25</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Segata</surname> <given-names>N</given-names>
</name>
<name>
<surname>Izard</surname> <given-names>J</given-names>
</name>
<name>
<surname>Waldron</surname> <given-names>L</given-names>
</name>
<name>
<surname>Gevers</surname> <given-names>D</given-names>
</name>
<name>
<surname>Miropolsky</surname> <given-names>L</given-names>
</name>
<name>
<surname>Garrett</surname> <given-names>WS</given-names>
</name>
<etal/>
</person-group>. <article-title>Metagenomic Biomarker Discovery and Explanation</article-title>. <source>Genome Biol</source> (<year>2011</year>) <volume>12</volume>(<issue>6</issue>):<fpage>R60</fpage>. doi: <pub-id pub-id-type="doi">10.1186/gb-2011-12-6-r60</pub-id>
</citation>
</ref>
<ref id="B26">
<label>26</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Douglas</surname> <given-names>GM</given-names>
</name>
<name>
<surname>Maffei</surname> <given-names>VJ</given-names>
</name>
<name>
<surname>Zaneveld</surname> <given-names>JR</given-names>
</name>
<name>
<surname>Yurgel</surname> <given-names>SN</given-names>
</name>
<name>
<surname>Brown</surname> <given-names>JR</given-names>
</name>
<name>
<surname>Taylor</surname> <given-names>CM</given-names>
</name>
<etal/>
</person-group>. <article-title>PICRUSt2 for Prediction of Metagenome Functions</article-title>. <source>Nat Biotechnol</source> (<year>2020</year>) <volume>38</volume>(<issue>6</issue>):<page-range>685&#x2013;8</page-range>. doi: <pub-id pub-id-type="doi">10.1038/s41587-020-0548-6</pub-id>
</citation>
</ref>
<ref id="B27">
<label>27</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Caspi</surname> <given-names>R</given-names>
</name>
<name>
<surname>Billington</surname> <given-names>R</given-names>
</name>
<name>
<surname>Keseler</surname> <given-names>IM</given-names>
</name>
<name>
<surname>Kothari</surname> <given-names>A</given-names>
</name>
<name>
<surname>Krummenacker</surname> <given-names>M</given-names>
</name>
<name>
<surname>Midford</surname> <given-names>PE</given-names>
</name>
<etal/>
</person-group>. <article-title>The MetaCyc Database of Metabolic Pathways and Enzymes - a 2019 Update</article-title>. <source>Nucleic Acids Res</source> (<year>2020</year>) <volume>48</volume>(<issue>D1</issue>):<page-range>D445&#x2013;D53</page-range>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkz862</pub-id>
</citation>
</ref>
<ref id="B28">
<label>28</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Parks</surname> <given-names>DH</given-names>
</name>
<name>
<surname>Tyson</surname> <given-names>GW</given-names>
</name>
<name>
<surname>Hugenholtz</surname> <given-names>P</given-names>
</name>
<name>
<surname>Beiko</surname> <given-names>RG</given-names>
</name>
</person-group>. <article-title>STAMP: Statistical Analysis of Taxonomic and Functional Profiles</article-title>. <source>Bioinformatics</source> (<year>2014</year>) <volume>30</volume>(<issue>21</issue>):<page-range>3123&#x2013;4</page-range>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btu494</pub-id>
</citation>
</ref>
<ref id="B29">
<label>29</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Weiss</surname> <given-names>S</given-names>
</name>
<name>
<surname>Van Treuren</surname> <given-names>W</given-names>
</name>
<name>
<surname>Lozupone</surname> <given-names>C</given-names>
</name>
<name>
<surname>Faust</surname> <given-names>K</given-names>
</name>
<name>
<surname>Friedman</surname> <given-names>J</given-names>
</name>
<name>
<surname>Deng</surname> <given-names>Y</given-names>
</name>
<etal/>
</person-group>. <article-title>Correlation Detection Strategies in Microbial Data Sets Vary Widely in Sensitivity and Precision</article-title>. <source>ISME J</source> (<year>2016</year>) <volume>10</volume>(<issue>7</issue>):<page-range>1669&#x2013;81</page-range>. doi: <pub-id pub-id-type="doi">10.1038/ismej.2015.235</pub-id>
</citation>
</ref>
<ref id="B30">
<label>30</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shannon</surname> <given-names>P</given-names>
</name>
<name>
<surname>Markiel</surname> <given-names>A</given-names>
</name>
<name>
<surname>Ozier</surname> <given-names>O</given-names>
</name>
<name>
<surname>Baliga</surname> <given-names>NS</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>JT</given-names>
</name>
<name>
<surname>Ramage</surname> <given-names>D</given-names>
</name>
<etal/>
</person-group>. <article-title>Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks</article-title>. <source>Genome Res</source> (<year>2003</year>) <volume>13</volume>(<issue>11</issue>):<page-range>2498&#x2013;504</page-range>. doi: <pub-id pub-id-type="doi">10.1101/gr.1239303</pub-id>
</citation>
</ref>
<ref id="B31">
<label>31</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tan</surname> <given-names>BK</given-names>
</name>
<name>
<surname>Vandekerckhove</surname> <given-names>P</given-names>
</name>
<name>
<surname>Kennedy</surname> <given-names>R</given-names>
</name>
<name>
<surname>Keay</surname> <given-names>SD</given-names>
</name>
</person-group>. <article-title>Investigation and Current Management of Recurrent IVF Treatment Failure in the UK</article-title>. <source>BJOG</source> (<year>2005</year>) <volume>112</volume>(<issue>6</issue>):<page-range>773&#x2013;80</page-range>. doi: <pub-id pub-id-type="doi">10.1111/j.1471-0528.2005.00523.x</pub-id>
</citation>
</ref>
<ref id="B32">
<label>32</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Franasiak</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Scott</surname> <given-names>RT</given-names>
</name>
</person-group>. <article-title>Endometrial Microbiome</article-title>. <source>Curr Opin Obstet Gynecol</source> (<year>2017</year>) <volume>29</volume>(<issue>3</issue>):<page-range>146&#x2013;52</page-range>. doi: <pub-id pub-id-type="doi">10.1097/GCO.0000000000000357</pub-id>
</citation>
</ref>
<ref id="B33">
<label>33</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kawai</surname> <given-names>T</given-names>
</name>
<name>
<surname>Akira</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>The Role of Pattern-Recognition Receptors in Innate Immunity: Update on Toll-Like Receptors</article-title>. <source>Nat Immunol</source> (<year>2010</year>) <volume>11</volume>(<issue>5</issue>):<page-range>373&#x2013;84</page-range>. doi: <pub-id pub-id-type="doi">10.1038/ni.1863</pub-id>
</citation>
</ref>
<ref id="B34">
<label>34</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dong</surname> <given-names>C</given-names>
</name>
</person-group>. <article-title>Cytokine Regulation and Function in T Cells</article-title>. <source>Annu Rev Immunol</source> (<year>2021</year>) <volume>39</volume>:<fpage>51</fpage>&#x2013;<lpage>76</lpage>. doi: <pub-id pub-id-type="doi">10.1146/annurev-immunol-061020-053702</pub-id>
</citation>
</ref>
<ref id="B35">
<label>35</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Greenwood</surname> <given-names>SM</given-names>
</name>
<name>
<surname>Moran</surname> <given-names>JJ</given-names>
</name>
</person-group>. <article-title>Chronic Endometritis: Morphologic and Clinical Observations</article-title>. <source>Obstet Gynecol</source> (<year>1981</year>) <volume>58</volume>(<issue>2</issue>):<page-range>176&#x2013;84</page-range>. doi: <pub-id pub-id-type="doi">10.1016/0378-5122(81)90041-4</pub-id>
</citation>
</ref>
<ref id="B36">
<label>36</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hanna</surname> <given-names>J</given-names>
</name>
<name>
<surname>Mandelboim</surname> <given-names>O</given-names>
</name>
</person-group>. <article-title>When Killers Become Helpers</article-title>. <source>Trends Immunol</source> (<year>2007</year>) <volume>28</volume>(<issue>5</issue>):<page-range>201&#x2013;6</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.it.2007.03.005</pub-id>
</citation>
</ref>
<ref id="B37">
<label>37</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>El-Azzamy</surname> <given-names>H</given-names>
</name>
<name>
<surname>Dambaeva</surname> <given-names>SV</given-names>
</name>
<name>
<surname>Katukurundage</surname> <given-names>D</given-names>
</name>
<name>
<surname>Salazar Garcia</surname> <given-names>MD</given-names>
</name>
<name>
<surname>Skariah</surname> <given-names>A</given-names>
</name>
<name>
<surname>Hussein</surname> <given-names>Y</given-names>
</name>
<etal/>
</person-group>. <article-title>Dysregulated Uterine Natural Killer Cells and Vascular Remodeling in Women With Recurrent Pregnancy Losses</article-title>. <source>Am J Reprod Immunol</source> (<year>2018</year>) <volume>80</volume>(<issue>4</issue>):<fpage>e13024</fpage>. doi: <pub-id pub-id-type="doi">10.1111/aji.13024</pub-id>
</citation>
</ref>
<ref id="B38">
<label>38</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sauerbrun-Cutler</surname> <given-names>MT</given-names>
</name>
<name>
<surname>Huber</surname> <given-names>WJ</given-names>
</name>
<name>
<surname>Krueger</surname> <given-names>PM</given-names>
</name>
<name>
<surname>Sung</surname> <given-names>CJ</given-names>
</name>
<name>
<surname>Has</surname> <given-names>P</given-names>
</name>
<name>
<surname>Sharma</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Do Endometrial Natural Killer and Regulatory T Cells Differ in Infertile and Clinical Pregnancy Patients? An Analysis in Patients Undergoing Frozen Embryo Transfer Cycles</article-title>. <source>Am J Reprod Immunol</source> (<year>2021</year>) <volume>85</volume>(<issue>6</issue>):<fpage>e13393</fpage>. doi: <pub-id pub-id-type="doi">10.1111/aji.13393</pub-id>
</citation>
</ref>
<ref id="B39">
<label>39</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kammerer</surname> <given-names>U</given-names>
</name>
<name>
<surname>von Wolff</surname> <given-names>M</given-names>
</name>
<name>
<surname>Markert</surname> <given-names>UR</given-names>
</name>
</person-group>. <article-title>Immunology of Human Endometrium</article-title>. <source>Immunobiology</source> (<year>2004</year>) <volume>209</volume>(<issue>7</issue>):<page-range>569&#x2013;74</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.imbio.2004.04.009</pub-id>
</citation>
</ref>
<ref id="B40">
<label>40</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robertson</surname> <given-names>SA</given-names>
</name>
<name>
<surname>Moldenhauer</surname> <given-names>LM</given-names>
</name>
</person-group>. <article-title>Immunological Determinants of Implantation Success</article-title>. <source>Int J Dev Biol</source> (<year>2014</year>) <volume>58</volume>(<issue>2-4</issue>):<page-range>205&#x2013;17</page-range>. doi: <pub-id pub-id-type="doi">10.1387/ijdb.140096sr</pub-id>
</citation>
</ref>
<ref id="B41">
<label>41</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yan</surname> <given-names>J</given-names>
</name>
<name>
<surname>Horng</surname> <given-names>T</given-names>
</name>
</person-group>. <article-title>Lipid Metabolism in Regulation of Macrophage Functions</article-title>. <source>Trends Cell Biol</source> (<year>2020</year>) <volume>30</volume>(<issue>12</issue>):<page-range>979&#x2013;89</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.tcb.2020.09.006</pub-id>
</citation>
</ref>
<ref id="B42">
<label>42</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kempkes</surname> <given-names>RWM</given-names>
</name>
<name>
<surname>Joosten</surname> <given-names>I</given-names>
</name>
<name>
<surname>Koenen</surname> <given-names>H</given-names>
</name>
<name>
<surname>He</surname> <given-names>X</given-names>
</name>
</person-group>. <article-title>Metabolic Pathways Involved in Regulatory T Cell Functionality</article-title>. <source>Front Immunol</source> (<year>2019</year>) <volume>10</volume>:<elocation-id>2839</elocation-id>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2019.02839</pub-id>
</citation>
</ref>
<ref id="B43">
<label>43</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chaouat</surname> <given-names>G</given-names>
</name>
</person-group>. <article-title>The Th1/Th2 Paradigm: Still Important in Pregnancy</article-title>? <source>Semin Immunopathol</source> (<year>2007</year>) <volume>29</volume>(<issue>2</issue>):<fpage>95</fpage>&#x2013;<lpage>113</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00281-007-0069-0</pub-id>
</citation>
</ref>
<ref id="B44">
<label>44</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Candelli</surname> <given-names>M</given-names>
</name>
<name>
<surname>Franza</surname> <given-names>L</given-names>
</name>
<name>
<surname>Pignataro</surname> <given-names>G</given-names>
</name>
<name>
<surname>Ojetti</surname> <given-names>V</given-names>
</name>
<name>
<surname>Covino</surname> <given-names>M</given-names>
</name>
<name>
<surname>Piccioni</surname> <given-names>A</given-names>
</name>
<etal/>
</person-group>. <article-title>Interaction Between Lipopolysaccharide and Gut Microbiota in Inflammatory Bowel Diseases</article-title>. <source>Int J Mol Sci</source> (<year>2021</year>) <volume>22</volume>(<issue>12</issue>). doi: <pub-id pub-id-type="doi">10.3390/ijms22126242</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>