Gut Microbiota Composition Changes in Constipated Women of Reproductive Age

Introduction

Chronic constipation is a common symptom-based gastrointestinal disease without organic lesions. Its prevalence ranges from 2.6% to 26.9% in the general population, with a median value of 16% in all adults (Mugie and Benninga, 2011; Schmidt, 2014), and disproportionately affecting older adults and women. Chronic constipation often has a striking negative effect on quality of life (Belsey et al., 2010) and has been associated with a reduction in work productivity (Sun et al., 2011). Constipation is not only harmful to the intestinal tract itself but also can affect other diseases, and is known to aggravate cardiovascular and cerebrovascular diseases, leading to mental disorders (Li et al., 2016). Since it is not a life-threatening disease, and it is not a medical emergency, chronic constipation is often overlooked and not reported by patients suffering from it.

Chronic constipation is a multifactorial gastrointestinal disease, and its pathophysiological mechanism is still unclear. Constipation has previously been studied in the field of intestinal function, and only recently have new studies revealed that the gut microbiota of constipated patients differs from that of healthy controls. Some previous studies suggested there is a relationship between constipation and some altered abundance of certain species in the fecal microbiome using culture-based methods (Zoppi et al., 1998; Khalif et al., 2005). New advances in molecular biology methods have replaced traditional culture-based methods and are the current standard approach to analyze gut microbiota. These tools have revealed the tremendous diversity, richness, and functional capacity of human microbiome. Zhu et al. (2014) used 16S rRNA gene sequencing to prove that children with FC (functional constipation, FC) had a significantly decrease of Bacteroidetes, especially Prevotella, and a significant increase of some subgenera of Firmicutes, such as Lactobacillus. – used qPCR to determine that adults with FC had a significant loss in the abundance of Bifidobacterium and Bacteroides. Parthasarathy et al. (2016) used 16S rRNA metagenomics analysis to demonstrate that female patients suffering from chronic constipation or IBS-C (irritable bowel syndrome-constipation, IBS-C) had a significantly higher level of Bacteroidetes in their colonic mucosa.

These inconsistent findings are hard to explain. However, growing evidence indicates that age, sex, and hormone levels can affect the composition of intestinal flora (Bennett et al., 2016). Cross-sectional studies of fecal samples from adult individuals in various age groups suggest that there are age-related changes in gut microbiota composition and diversity. Longitudinal analyses indicate that the intestinal microbial population of healthy individuals is relatively stable for decades (Claesson et al., 2011). However, old age is associated with a more diverse and variable gut microbiome (Faith et al., 2013). With aging, the equilibrium state between intestinal flora and host worsens and gradually reaches a stage of dysbiosis (Kim, 2018; Riaz Rajoka et al., 2018; Xu and Zhu, 2019). In addition, some studies have shown an interaction between the microbiome and the endocrine system. Some bacteria can produce hormones and regulate the host’s hormonal homeostasis by inhibiting gene transcription. Similarly, host hormones may affect bacterial gene expression and bacterial growth and have consequences on host physiology (Rizzetto et al., 2018). The sex of the host also influences the gut microbiome and affects disease susceptibility, and these differences are the result of the actions of sex hormones (Ma, 2019). Importantly, all of these factors known to affect the gut microbiota are also known to influence susceptibility to constipation (Houghton et al., 2016; Lu and Velasco-Benítez, 2017; Shin et al., 2019).

It is conceivable that the large number of gut microbiota residing in the intestinal tract influences intestinal function. Therefore, in this study we examine the overall structure of gut microbiota in constipated women of reproductive age to reduce the effects of age, sex, and hormone levels. We aim to provide a clearer understanding of the role of the microbiome and to facilitate further research on the relationship between chronic constipation and gut microbiota.

Materials and Methods

Human Subjects

All participants were assessed using defecation related scales and a dietary habits questionnaire. Most of them were recruited from among the female staff of reproductive age of Shanxi Bethune Hospital. The constipated group needed to meet the Rome IV diagnosis standards of constipation and the normal group had a healthy lifestyle without any symptoms of constipation. All had a normal body mass index (BMI), indexes ranging from 18 to 25. All participants voluntarily enrolled and signed informed consent forms. The exclusion criteria for the study are summarized as follows: no children, men, pregnant women, or postmenopausal women; no patients with metabolic diseases, neuropsychiatric diseases, cancer, or such conditions as diabetes or Parkinson’s disease; they must have never undergone intestinal surgery; they must not have taken laxatives, antibiotics, probiotics, prebiotics, non-steroidal antiinflammatory drugs (NSAIDs), opioid drugs, traditional Chinese medicines, proton pump inhibitors (PPIs), or histamine receptor antagonists in the month preceding sample collection. All the fecal samples were processed on ice, stored in the Eppendorf tubes, and then transferred to a –80°C freezer within 30 min of collection. The study was approved by the ethics committee of Shanxi Bethune Hospital (Certificate No. XYLL-2019-124).

Detection Methods

Genomic DNA Extraction and PCR Amplification

Microbial DNA of the samples was extracted using SDS (sodium dodecyl sulfate, SDS) method and then purity and concentration of DNA were assessed using agarose gel electrophoresis. An appropriate amount of sample DNA was taken in a centrifuge tube, and the sample was diluted to 1 ng/1g with sterile water. Using the diluted genomic DNA as the template, the V3–V4 regions of the bacterial 16S ribosomal RNA gene were amplified by PCR using the linker primer sequence CCTAYGGGRBGCASCAG, GGACTACNNGGGTATCTAAT, where the barcode is a twelve-base sequence unique to each sample. To ensure the efficiency and accuracy of amplification, Phusion high-fidelity PCR Master Mix with GC Buffer from New England Biolabs was used for PCR.

PCR Amplicon Purification and NovaSeq Sequencing

According to the concentration of duplicate PCR amplicons, the samples were mixed in equal amounts. Barcoded amplicons were recovered using gel recovery kits. TruSeq DNA PCR-free Sample Preparation Kit was used to construct the library. The constructed library was quantified using Qubit and q-PCR. After the library was established, Illumina NovaSeq6000 was used for computer sequencing.

Information Analysis

Sequencing Data Processing

There are a certain proportion of dirty data among the raw data obtained from sequencing. In order to render the results of information analysis more accurate and reliable, the original data were first spliced and filtered to provide clean data. Then chimera filtering was conducted to obtain effective tags for subsequent analysis.

OTU Clustering and Species Annotation

In order to study the species composition of each sample, effective tags from all samples were clustered by OTUs (operational taxonomic units) with 97% identity using Uparse v7.0.1001 (http://www.drive5.com/uparse/).The sequence with the highest frequency of OTUs was selected as the representative sequence of OTUs. According to the OTU clustering results, species annotation was made for the representative sequence of each OTU to obtain the corresponding species information and species-based abundance distribution by using the Mothur method and SSUrRNA(sigmasubunit rRNA) database of SILVA132 (http://www.arb-silva.de/) and by using MUSCLE (version 3.8.31, http://www.drive5.com/muscle/) software for multiple sequence alignment.

Advanced Analysis

The smallest amount of data in the samples was taken as the standard for homogenization of all samples. Based on the homogenized data, the subsequent alpha diversity analysis was performed to find species richness and evenness information. We used Qiime software (version 1.9.1) to calculate Chao1 and Shannon indices; we used R software (version 2.15.3) to draw dilution and Rank abundance curves, all of which were for Alpha diversity index analysis of difference between the groups. Beta diversity analysis was performed to indicate the relative abundance distribution and significant differences of the gut microbiota between the groups. Metastats analysis used R software to perform permutation tests among groups under Phylum Class Order Family Species to determine P values. LEfSe software was used for LEfSe analysis, and the default screening value of LDA score was set to 4, which can be used for comparison of two or more groups. It emphasizes statistical significance and biological correlation, and it can be used to identify statistically different biomarkers between groups. Finally, we used Tax4Fun to predict the functional and metabolic capabilities of microbial communities. Tax4Fun is an open-source R package that predicts the functional or metabolic capabilities of microbial communities based on 16S rRNA datasets. The function of the microbial communities was predicted by constructing a linear relationship between the SILVA (from Latin silva, meaning forest, an improved data processing and web-based tool) classification and the KEGG (Kyoto Encyclopedia of Genes and Genomes) database.

Results

Species Richness and Diversity

The fecal microbiomes of 29 constipated women of reproductive age and 30 healthy controls were analyzed by 16S rRNA gene sequencing, because one patient with constipation failed to collect stool sample. The characteristics of different group subjects are summarized in Table 1. There were no significant differences in age, body mass index, alanine transaminase, aspertate aminotransferase, fasting blood glucose, total bilirubin, albumin, total bile acid, cholesterol, triglyceride, high density lipoprotein, and low density lipoprotein of the two groups. A total of 7,414,293 sequencing reads were obtained from the 59 samples. To assess the alpha diversities of fecal flora in each group, five metrics were calculated: Venn diagrams, rarefaction curve, rank abundance curves, Shannon plots, and Chao1 plots, showing species enrichment and distribution. According to the OTU results and research requirements obtained by clustering, the common and unique OTUs among the groups were analyzed and displayed in a Venn diagram. The two groups had 1,255 common OTUs, and the 196 unique OTUs in the constipated group were lower than 297 in the control group (Figure 1A), indicating that the species diversity in the constipated group was reduced. The rarefaction curve for the healthy group was slightly higher than the constipated group, reflecting the richness in the healthy group. When the curve tended to be flat, it indicates that the sequencing data volume was considered reasonable. The slope of the rank abundance curve reflects the richness and evenness of the species in the sample. The richness of the species is reflected by the width of the curve. In the vertical direction, the smoothness of the curve reflects the uniformity of species in the sample. However, no significant difference was found in either the Shannon or Chao1 index (p > 0.05 for both), suggesting that the community richness and diversity of the two groups were approximately the same.

Table 1 Characteristics of the study groups.

Figure 1 Species richness and diversity of the gut microbiome in two groups. N refers to the normal group; C refers to the constipated group. (A) Species richness and evenness based on the Venn diagram. The species richness in the constipated group is lower than that in the healthy group. (B) The rarefaction curve can directly reflect the rationality of sequencing data and indirectly reflect the richness of species in two groups. (C) The rank abundance curve reflects the richness and evenness of the species in two groups. (D) The Chao1 plot reflects species richness. (E) The Shannon plot reflects species diversity.

Relative Abundance Distributions at Different Levels of the Gut Microbiome in Two Groups

The relative abundance distribution at different levels of the gut microbiome in the two groups summarized in Table 2, which reflects the dominant flora and their relative proportion. The overall microbial composition at the phylum level in each sample is shown in Figure 2A, and the group data are shown in Figure 2B. A total of 25 bacterial phyla were detected in the gut microbiome of the two groups, including four fundamental phyla (Firmicutes, Bacteroidetes, Proteobacteria, and Actinobacteria), and other minor phyla. Firmicutes was the most prominent phylum, taking up 53.49 and 55.76%, respectively, of the gut microbiota in constipated women of reproductive age and healthy control group. The major families were Lachnospiraceae, Ruminococcaceae, Veillonellaceae, Peptostreptococcaceae, Lactobacillaceae, and Streptococcaceae. Bacteroidetes was the second major phyla in both groups, with a proportion of 35.14 and 25.22% in the two groups, including Bacteroidaceae, Prevotellaceae, Rikenellaceae, and Tannerellaceae. Proteobacteria, mostly Enterobacteriaceae, constituted the third most abundant phylum in the healthy control group, with a proportion of 11.49%, and the fourth most abundant phylum in the constipated group, with a proportion of 4.00%. Actinobacteria, mainly containing Bifidobacteriaceae, was the fourth most abundant phylum in the healthy control group and the third most abundant phylum in the constipated group (5.11 vs. 6.41%). According to the above analysis, there was a smaller ratio of Firmicutes/Bacteroidetes (F/B) in the constipated group than that in the healthy control group (1.52 vs. 2.21).

Table 2 The proportion and significant differences in the gut microbiome at different levels in the normal and constipated groups.

Figure 2 Relative abundance distribution at different levels of the gut microbiome in two groups. Here we analyze mainly at the level of Phylum, Family, Genus. (A) The relative abundance histogram of phyla in each sample. (B) The predominant phyla in the two groups. (C) The predominant family in the two groups. (D) The predominant genera in the two groups.

At the family and the genus level, the top 10 most highly abundant bacteria in the healthy control group and constipated women of reproductive age were identifed. As shown in Figures 2C, D, these was mainly included Lachnospiraceae (with a proportion of 24.99 and 18.73%), Ruminococcaceae (21.92 vs. 22.50%), Bacteroidaceae (with a proportion of 16.91 and 23.46%), and Enterobacteriaceae (with a proportion of 10.60 and 3.24%) at the family level in two groups respectively, and mainly consisted of Bacteroidetes (16.91 vs. 23.46%), Feacalibacterium (10.67 vs. 10.97%), Bifidobacterium (5.16 vs. 4.14%), Blautia (6.82 vs. 3.71%), Agathobacter (5.29 vs. 4.68%), and Klebsiella (3.57 vs. 0.77%) at the genus level in two groups, respectively.

Significant Differences in the Microbiota of Constipated Women of Reproductive Age and Healthy Controls

Principal component analysis indicated a remarkable separation in the fecal microbial composition between constipated women of reproductive age and healthy controls, with the first two principal component scores of PC1 (25.35%) and PC2 (12.04%) of explained variance (Figure 3A). The box diagram revealed that there was a significant difference in beta diversity between the two groups as determined by an unweighted t test (Figure 3B, p < 0.01).

Figure 3 Significant differences in fecal microbiota between constipated women of reproductive age and healthy control subjects. (A) Principal component analysis of the two groups. (B) Box plots of beta diversity. (C) Comparison of the top 10 phyla. (D) Comparison of the top 10 families. (E) Comparison of the top 35 genus. (F) Distribution of LDA values displayed as a histogram, with an LDA value of >4 indicating statistical significance between the two groups. (G) Phylogenetic distribution of the microbiome of the two groups.

The p values were obtained by hypothesis testing of species abundance data from different levels using the MetaStat method. A p-value of less than 0.05 or 0.01 was used to select the top 10 species. The data are shown in Table 2 and are displayed as a heatmap with annotations in Figure 3. At the phylum level, Bacteroidetes was higher in the constipated group than that in the healthy control group (p < 0.05); nevertheless, Proteobacteria was lower in the constipated group than that in the healthy control group (p < 0.01). No differences were found in the other top 10 phyla of the two groups (Figure 3C). At the family level, the heatmap showed a significant increase in Bacteroidaceae and Lactobacillaceae (p < 0.05), and a significant decrease in Enterobacteriaceae (p < 0.01). The level of Lachnospiraceae (p = 0.0659) and Prevotellaceae (p = 0.0639) showed marginal values. However, there were no significant differences between the two groups for some common bacteria, such as Ruminococcaceae, Peptostreptococcaceae, Veillonellaceae, and Bifidobacteriaceae (Figure 3D). At the genus level, the top 35 most highly abundant bacteria were selected, as shown in Figure 3E. Bacteroides (p < 0.05) and Parabacteroides (p < 0.01) were more abundant in the constipated group; Fusicatenibacter (p < 0.05), Citrobacter (p < 0.05), and Roseburia (p < 0.01) were less abundant in the constipated group than those in the healthy group.

Using linear discriminant analysis coupled with effect size measurements, we confirmed that the phylum Bacteroidetes, the class Bacteroidia, the order Bacteroidales, the family Bacteroidaceae, and the genus Bacteroides were higher in the gut microbiota in constipated women of reproductive age, while the phylum Proteobacteria, the class Gammaproteobacteria, the order Enterobacteriales, the family Enterobacteriaceae, and the genus unidentified Enterobacteriaceae were higher in the gut microbiota in healthy controls (Figure 3F). The phylogenetic distribution of the microbiome of the two groups show a good correlation from phylum to genus level, suggesting that Bacteroides and Proteobacteria may be the crucial players involved in the pathogenesis of chronic constipation in women of reproductive age (Figure 3G).

Functional Predictions for the Gut Microbiome in Constipated Women of Reproductive Age and Healthy Controls

We used Tax4Fun to predict the functional and metabolic capabilities of microbial communities by constructing a linear relationship between the SILVA classification and the KEGG database. In the KEGG pathway annotations, the relative abundance of microbial genes in both groups was high in carbohydrate metabolism, amino acid metabolism, membrane transport of environmental information processing, translation and replication, and repair of genetic information process (Figure 4A). There were 28 KOs (KEGG Orthologues) be confirmed with significant differences in the fecal microbiomes between the two groups (p < 0.05; Figure 4B). Based on level 3 of KEGG pathway analyses, sphingolipid metabolism, cyanoamino acid metabolism, amino acid metabolism and biotin metabolism, transport, phenylpropanoid biosynthesis, glycosphingolipid biosynthesis, polyketide sugar unit biosynthesis, and various types of N-glycan biosynthesis were more abundance in fecal microbiome of the constipated group. The microbial gene functions involved ribosome biogenesis, glucagon signaling pathways, and naphthalene degradation were more abundant in the gut microbiome of the healthy normal control group (p < 0.05; see Figure 4C).

Figure 4 Functional predictions for the fecal microbiome in constipated women of reproductive age and healthy normal controls. (A) The relative abundance of gene in the KEGG pathway. (B) The KOs in the fecal microbiome of the two groups were shown to have significant differences (p < 0.05), which were identified using Tax4Fun. (C) The significant differences (p < 0.05) in KEGG pathways at level 3 for the fecal microbiome in constipated women of reproductive age and healthy normal individuals.

Discussion

In this cross-sectional study, we compared fecal microbiota from constipated women of reproductive age and healthy controls, based on 16S rRNA sequencing to characterize the overall microbial differences. In contrast to previous studies, which enrolled men and women and a wide range of ages, from 18 to 80 years, we only included women of child-bearing age in our study to reduce the effects of gender, age, and hormone levels. We find that the ecological diversity and richness in the fecal microbiome in constipated women of reproductive age were similar to those in the healthy control individuals, whereas the differences at all taxonomic levels of the fecal microbiome between the two groups were significant, indicating that constipation is associated with an altered microbiome in the gut. As a whole, the fecal microbiome of constipated women of reproductive age exhibited an increased level of Bacteroidetes and decreased level of Proteobacteria, which was mostly explained by the increased abundance in the genus Bacteroides and the reduced numbers of Enterobacteriales. Significantly decreased levels of Proteobacteria and increased levels of Bacteroidetes resulted in a reduced ratio of F/B (Firmicutes to Bacteroidetes). Furthermore, the abundance of some butyrate-producing bacteria was also reduced.

The human colon contains a diverse microbial community, which is inhabited by hundreds of distinct species. Of these, 25% are Bacteroides (Eckburg et al., 2005). Bacteroides are strictly anaerobic, Gram-negative, non-motile, rod-shaped, and non-spore-forming bacteria, and comprises more than 92 species (Smith et al., 2006; Mancabelli et al., 2017). Bacteroides, the predominant genus within the human gut microbiota, usually plays a crucial role in degradation and fermentation of organic matter in the colon and are beneficial symbionts with their hosts (Salyers, 1995). Beneficial symbiosis requires the bacteria to sense changes in the environment so that they can adapt to alterations in their surroundings. The abundance of Bacteroides increased by 40% compared to that of healthy controls in this study, indicating that the increase in Bacteroides numbers may be based on the host’s constipated gut environment. The increased relative abundance of certain taxa in the presence of a gut disorder may not reflect a taxa-specific role in pathogenesis but may be linked to a global alteration of the gut microbiota’s homeostasis (Freitas et al., 2003). However, Bacteroides may simply need to turn on certain genes to change from friendly commensal to dangerous threat (Wexler, 2007). The best support is that Bacteroides also associated with infections, such as colitis and pouchitis (Shepherd et al., 1989; Kmiot et al., 1993), which suggests that Bacteroides may affect mucosal structure. Therefore, the increase of Bacteroidetes may also indicate an increase in pathogenic bacteria. What role Bacteroidetes plays in the constipated women of childbearing age requires further research and exploration.

Interestingly, we found a significant decrease in the abundance of the phylum Proteobacteria, in which Enterobacteriaceae was reduced. This result is similar to a previously published study that found that Enterobacteriaceae increased after constipated mice were given irritant laxatives (Takayama et al., 2019), revealing an association between high levels of the family Enterobacteriaceae and purgative activity. The family Enterobacteriaceae includes diarrheal pathogens such as Shigella and Salmonella. Furthermore, Citrobacter rodentium was found to be reduced in our study, and this bacterium belongs to the family Enterobacteriaceae, a close relative of the human diarrheal pathogen enterohemorrhagic and enteropathogenic Escherichia coli (Tsai et al., 2017; Zhu et al., 2018). Despite the fact that these diarrheal pathogens were not detected in our study, we still observed a reduction in E. coli in the constipated group as determined in the linear discriminant analysis. We assumed that an unclassified genus of the family Enterobacteriaceae might include similar bacteria, the reduction of which is associated with the pathogenesis of constipation. Consistent with the in vitro motility studies using human colon specimen, E. coli strain Nissle 1917 has been reported to promote gastrointestinal motility and muscle cell contractility (Bär et al., 2009).

Additionally, Fusicatenibacter and Roseburia were reduced in constipated women of reproductive age, and these are well-known butyrate-producing bacteria of the Lachnospiraceae family in the Firmicutes phylum. Butyrate producers are an abundant and phylogenetically diverse group of bacteria that are likely to play an important role in maintaining gut health, primarily through the production of butyrate. Studies have shown that the higher the abundance of butyrate-producing bacteria, the faster the colonic transit as a result of the influence of butyrate on gastrointestinal dynamics (Chassard et al., 2012). The butyrate produced by bacteria may accelerate colonic motility by stimulating the release of serotonin or promoting cholinergic pathways (Soret et al., 2010; Reigstad et al., 2015). A recent study shows that an abundance of butyrate-producing gut bacteria relieves constipation symptoms via short-chain fatty acids production and hormone secretion (Zhuang et al., 2019). However, we found no changes in the bacteria of Butyrate-metabolic-capability in the Tax4Fun analysis. This suggests that we need to look at the role of butyric acid-producing bacteria in the future using metabolomics.

The data of functional analysis suggest that the alteration of pathways involved in metabolism, biosynthesis, genetic information process, and environmental information processing is associated with chronic constipation. The 16S rRNA gene sequencing method has its own disadvantages as it could not identify microorganisms at the species level or the strain level or provide direct data on the crucial changes in the functionality of the microbiota (Aßhauer et al., 2015). Future studies should include samples from men and women of various ages, ethnic origins, geographical regions, and dietary differences, and apply metagenomics sequencing analysis techniques to explore the changes in the function of gut microbial genes.

Conclusion

Overall, we find evidence for gut microbiota dysbiosis in constipated women of reproductive age by discussing the fecal microbiota compositional shifts in case groups as compared to healthy controls. The variations may predict the unique interactions between hosts and some certain bacteria, or some specific bacterial metabolic products, which may aid future exploration of the pathogenesis of constipation. Narrowing this study to a small sub-population in order to avoid the influences of age, sex, and hormonal differences had a clear impact on the meaningfulness of the data and made data interpretation more straightforward.

Notably, the taxonomic characters of microorganisms extracted from fecal samples included both the indigenous flora as well as the allochthonous microbiota and thus may not be completely representative of the resident gut microbial population. Nevertheless, analysis of fecal flora is also a promising method for a rapid screen with the aim of identifying biomarkers associated with chronic constipation, since the collection of fecal samples is noninvasive and no special clinical procedures are required. The limitations of the present study should be considered. We did not detect structural changes in colonic mucosal flora, nor did we perform shotgun metagenome analysis to understand intestinal flora function. We did not perform metabolomics analyses to explore the pathogenesis of constipation at a molecular level. Our study preliminarily identified some potential microflora, which still requires further validation based on clinical samples and animal models. Elucidating the differences in the fecal microbiome of female constipated patients of reproductive age may provide the foundation to improve our understanding of the pathogenesis of chronic constipation in specific populations and could support the development of novel therapeutic options aimed at modifying the gut microbiota.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/. SRA accession: PRJNA636012.

Ethics Statement

The studies involving human participants were reviewed and approved by the ethics committee of Shanxi Bethune Hospital. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

XD coordinated the project and conceived of the study. HL and JC recruited patients and conducted the clinical trials. XR and CY participated in its design and coordination. SL and XB collected clinical samples. SS helped refine the design of the clinical trials and helped interpret the data. HL performed the statistical analysis, and prepared and revised the manuscript. All authors had read and approved the final manuscript.

Funding

This research was supported by the Natural Science Foundation of Applied Basic Research of Shanxi Province, China (201801D121199) and the Research Project of Shanxi Provincial Health and Family Planning Commission, China (2018019).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential or actual conflict of interest.

Acknowledgments

We would like to thank Biotree (www.biotree.com.cn) for providing the 16S rRNA gene Sequence and the bioinformatics analysis. We also thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript.

References