Altered Fecal Microbiota Signatures in Patients With Anxiety and Depression in the Gastrointestinal Cancer Screening: A Case-Control Study

Background: Increasing attention has been devoted to cancer screening and microbiota in recent decades, but currently there is less focus on microbiota characterization among screeners and its relationship to anxiety and depression. Methods: We characterized the microbial communities of fecal samples collected through the FOBT card from anxiety and depression screeners and paired controls in Henan, China (1:2, N = 69). DNA was extracted using the MOBIO PowerSoil kit. The V4 region of the 16S rRNA gene was sequenced using MiniSeq and processed using QIIME1. LEfSe was used to identify differentially abundant microbes, the Wilcoxon rank-sum test was used to test alpha diversity differences, and permutational multivariate analysis of variance was used to test for differences in beta diversity. Results: Similar fecal microbiota signatures in composition were found among screeners. The intestinal microbial environments by phylum were all composed primarily of Firmicutes, Bacteroidetes, and Proteobacteria, and the corresponding top genera were Faecalibacterium, Roseburia, and Prevotella. Compared with controls, the ranking of the top five genera in the anxiety and depression group changed, and the dominant genus was Prevotella in the anxiety and depression group and Faecalibacterium in the control group. There was a lower relative abundance of Gemmiger (1.4 vs. 2.3%, P = 0.025), Ruminococcus (0.6 vs. 0.8%, P = 0.037), and Veillonella (0.6 vs. 1.3%, P = 0.020). This may be linked to the lower alpha diversity in participants with anxiety and depression (Observed OTUs: 122.35 vs. 143.24; Chao1: 127.35 vs. 149.98), although no significant differences were observed. Distinct clustering in microbial composition between the two groups was detected for the Jaccard distance (P = 0.011). Conclusions: Our study showed differing microbial characterization among participants with anxiety and depression in the endoscopic screening of upper gastrointestinal cancer. Gemmiger, Ruminococcus, and Veillonella were informative and have potential clinical implications, which need to be confirmed by large-scale, prospective cohort studies and biological mechanism research.


INTRODUCTION
In excess of 100 trillion microorganisms colonize the human intestinal plot, which assumes a vital part in human wellbeing and illness conditions (1). The normal intestinal microbiota act as significant functions in host metabolism, xenobiotics, integrity maintenance of the intestinal mucosal barrier, immunomodulation, and assurance against microorganisms (2,3). The microbiota-gut-brain axis, a research hotspot, refers to the bidirectional communication between the microorganisms residing in the gut and our brain function, behavior, and emotion (4,5). Adequate evidence highlighted the multifaceted role of the intestinal microbiota in carcinogenesis (e.g., gastrointestinal cancer) and psychological distress (e.g., anxiety and depression disorders) (6)(7)(8)(9). Previous research showed that anxiety and depression patients were characterized by a higher abundance of proinflammatory species (Enterobacteriaceae and Desulfovibrio), lower microbiota diversity, and a lower abundance of short-chain fatty acid-producing species (Faecalibacterium) (10).
Current national cancer screening recommendations referenced the potential harm of mental health owing to cancer screening (11). As people increase their emphasis on health problems, studies on cancer screening and psychology represent a growing field. Invasive endoscopic screening for gastrointestinal cancer is often accompanied by negative psychosocial consequences to participants (anxiety and depression symptoms) (12,13). Considering the psychological distress in cancer screening and the microbiota-gut-brain axis, the relationship of microbiota and anxiety and depression among screeners becomes interesting. Discovering the microbial characteristics of screeners and microbiota diversity and characteristic genera affecting anxiety and depression would be valuable to optimize the strategy of cancer screening and reduce the negative psychological effects (anxiety and depression) caused by cancer screening. However, so far reliable evidence on microbiota characterization among screeners is limited and insufficient. No known research has investigated the relationship between intestinal microbiota and anxiety and depression among screeners. Therefore, the study aimed to explore the microbial characterization of participants in endoscopic screening and to identify psychological distress-associated gut microbiota.

Study Participants and Sample Collection
Based on the endoscopic screening of the National Cohort of Esophageal Cancer (NCEC) project in China, we retrospectively recruited permanent residents aged 40-69 years in August 2019. Fecal samples were collected before endoscopic screening of upper gastrointestinal cancer in Linzhou Cancer Hospital, Henan Province. Fecal sample collection process: (1) The fecal collection kit, including a fecal collection box and fecal occult blood test card (FOBT) for smearing feces, was prepared in advance. (2) The kit was distributed to the participants, and the sampling box was directly placed in the squatting stool by themselves. The fresh fecal collection was completed before endoscopy. (3) After defecation, the FOBT card was opened, the stool collection stick was used to pick up a small number of feces and smear them on the two panes of the FOBT card, and then the card was closed and the FOBT card was placed in the sealed bags. (4) The sealed bags were stored in the −80 • C refrigerator in the biobank of Linzhou Cancer Hospital and shipped to the laboratory with dry ice.
Only people with both anxiety and depression symptoms were regarded as anxiety and depression screeners. Those who had neither anxiety nor depression symptoms were regarded as the control group (screeners without anxiety and depression). The control group was matched by age and sex (1:2). A total of 69 participants were included, with 23 anxiety and depression screeners and 46 paired screeners in the control group.
All participants signed written informed consent. This study was approved by the Institutional Review Board of the Cancer Hospital of the Chinese Academy of Medical Sciences (No. 16-171/1250). Participants' sociodemographic information was gathered by trained staff via a uniform questionnaire.
Polymerase chain reaction (PCR) mixtures contained 1 µL of forward and reverse primer, 1 µL of template DNA, 4 µL of deoxyribonucleoside-triphosphates (dNTPs), 5 µL of 10× EasyPfu Buffer, 1 µL of EasyPfu DNA Polymerase, and 1 µL of double distilled water into a 50 µL total reaction volume. The PCR amplicons were quantified using the Qubit dsDNA HS Assay Kit (Thermo Fisher/Invitrogen Cat. no. Q32854, Waltham, USA) following the manufacturer's instructions. All sequencing was acted in a solitary MiniSeq run and exported in the FASTQ format. Illumina MiniSeq Reporter was carried out to remove adapter and primer sequences. All specimens collected were successfully amplified and sequenced. Sequencing data were performed with the Quantitative Insights into Microbial Ecology (QIIME2) platform (16). The raw sequences were processed to remove low-quality reads, under strict quality control and feature table construction using the Divisive Amplicon Denoising Algorithm 2 (DADA2) algorithm (17). A similarity threshold of 97% was matched. The taxonomic assignment of the sequence variants was assigned using the Greengenes 13_8 (18). The Shannon index rarefaction curve was represented in Supplementary Figure 1. A total of 23 positive anxiety and depression screeners with a mean of 51,272 reads and 46 non-anxiety and depression screeners with a mean of 60,697 reads were included in the analysis. Then we generated alpha diversity metrics and beta diversity metrics using QIIME.

Measurement of Anxiety and Depression
The anxiety symptoms were evaluated by the seven-item Generalized Anxiety Disorder (GAD-7), a widely used and acknowledged measurement tool worldwide. Good psychometrics of GAD-7 has been proved in primary medical care (19). The reliability of internal consistency of GAD-7 in the study was strong (Cronbach's alpha = 0.888). GAD-7 was used to identify anxiety symptoms of individuals in the past 2 weeks, with seven items and four responses (0 = never; 1 = sometimes; 2 ≥ half of day; 3 = almost every day). The anxiety score was calculated by adding the answers for each of the items and ranges from 0 to 21. The higher the score, the worse anxious symptoms. A result of five was regarded as the threshold for positive anxiety symptoms (20). The nine-item Patient Health Questionnaire (PHQ-9), one of the most well-known self-reported tools for assessing depression symptoms (21), has shown good performance for evaluating depressive disorder (22,23). Cronbach's alpha coefficients of PHA-9 in our study were 0.896. PHQ-9 was used to identify depressive symptoms of individuals in the past 2 weeks, with nine items and four responses (similar to GAD-7). The PHQ-9 score was the sum of each item. The higher the score, the worse the depression symptoms. People with a PHQ-9 score of 5 or higher were considered positive for depression symptoms (23).

Statistical Analysis
Chi-square tests and T-tests were used to compare basic characteristics between anxiety and depression screeners and the controls. Evenness index, Shannon index, Observed OTUs, and Chao1 index were used to reflect the alpha diversity. Differences in alpha diversity were analyzed between the anxiety and depression group and paired control group by Wilcoxon rank-sum test. Permutational Multivariate Analysis of Variance (PERMANOVA, R-vegan function adonis) was used to explore whether the flora composition differed by anxiety and depression status (beta diversity). Principal coordinate analysis (PCoA) was used to visualize clustering and find discrepancy among the independent β diversity matrices, based on Bary-Curtis dissimilarity, Jaccard distance, and weighted and unweighted unifrac distances. High relative abundance (≥0.01) genera were compared between the two groups by the Wilcoxon ranksum test.
Linear discriminant analysis effect size (LEfSe) was used to identify microbes associated with anxiety and depression symptoms (24). Using Wilcoxon rank-sum test, LEfSe detects microbiota with significant differences between the two groups. Microbiota, with linear discriminant analysis scores (LDA) ≥ 2.032, were identified as potential characteristic flora associated with anxiety and depression symptoms. P < 0.05 were considered statistically significant. All analyses were conducted using the software program R Studio (Version 1.1.456).

Baseline Characteristics of Participants
The average age of screeners in the anxiety and depression group and control group was 55.30 (SD = 8.20) and 55.63 (SD = 7.78), respectively. No statistically significant differences were observed for the baseline characteristics, including BMI, marital status, highest education level, household income, smoking status, alcohol drinking, hot food, and life satisfaction between the anxiety and depression group and the control group (see Table 1).

Alpha Diversity
As shown in Figure 1, the alpha diversity in the anxiety and depression group was decreased compared with the control group, although no significant differences were observed (Observed OTUs: 122. 35

Beta Diversity
When considering microbial community composition (i.e., beta diversity), significant clustering was found for the Jaccard distance (P = 0.011) between the anxiety and depression group and the control group but not for the Bray-Curtis dissimilarity, Weighted and Unweighted Unifrac distance (Figure 2).

Microbial Composition
Microbial relative abundances at the phylum, family, genus, and species levels for anxiety and depression group and control group were shown in Figure 3. Similar fecal microbiota signatures in composition were found between the two groups. At the phylum level, the intestinal microbial environments of the two groups were all comprised primarily of Firmicutes (66.6 vs. 73.4%, P = 0.240), Bacteroidetes (15.5 vs. 14.1%, P = 0.620), Proteobacteria (7.8 vs. 4.7%, P = 0.100), and Actinobacteria (2.2 vs.3.1%, P = 0.541). The Bacteroidetes and Proteobacteria in feces of patients with anxiety and depression increased, while the Firmicutes and Actinobacteria decreased, although there was no significant difference.
At the genus level, compared with controls, the prevotella and Blautia in feces of anxious and depressed patients showed an increasing trend, while the Faecalibacterium and Roseburia showed a downward trend, although the difference is not statistically significant. The relative abundance of prevotella (10.9 vs. 10.0%, P = 0.600), Faecalibacterium (10.5 vs. 11.8%, P = 0.470), Blautia (10.2 vs. 9.8%, P = 0.720), Roseburia (9.2 vs. 10.9%, P = 0.600), and Escherichia (3.9 vs. 2.6, P = 0.679) were the top five in both groups, but the ranking of the top five genera changed for subjects with anxiety and depression symptoms. The top one genera were Prevotella in the anxiety and depression group and Faecalibacterium in the control group (Figure 3 and Supplementary Table 1). Microbial relative abundances of each sample for the anxiety and depression group and control group were shown in Supplementary Figure 2.

Microbial Diversity and Characteristic Genus
The microbiota associated with anxiety and depression symptoms from LEfSe is shown in

DISCUSSION
In this study, we investigated the microbial characterization of participants in gastrointestinal cancer screening and identified psychological distress-associated gut microbiota. The microbial environments of screeners all comprised primarily of Firmicutes, Bacteroidetes, and Proteobacteria at the phylum level and Faecalibacterium, Roseburia, and Prevotella at the genus level. Compared with the controls, the microbial characterization of screeners was distinct among participants with anxiety and depression, and the ranking of the top five genera in the anxiety and depression group changed. There was a lower relative abundance of Gemmiger, Ruminococcus, and Veillonella in participants with anxiety and depression, which was also reflected by the decreased alpha diversity in screeners who suffered psychological distress, although the difference was not significant. The findings filled the gap in the field of screening and contribute to a better understanding of endoscopic screening for gastrointestinal cancer and psychological distress, which would provide innovative strategies for relieving anxiety and depression and the optimization and implementation of endoscopic screening programs for upper gastrointestinal cancer in China. It was evidently clear from the findings that participation in endoscopic screening may increase screeners' anxiety and depression symptoms. Current national cancer screening recommendations referenced the potential mental harm owing to cancer screening (11). The role of psychological status on screening has been seriously underestimated (12,13). The psychological impact of the screening procedure itself is a common problem in all cancer screenings, but the psychological problem caused by endoscopic screening is more prominent due to its invasive nature, which presents a challenge to screeners' psychological states and emotions, increasing anxiety and depression levels. Taking esophageal cancer as an example, on the one hand, waiting for an invasive endoscopic examination may trigger or increase anxiety and depression levels. On the other hand, screeners are worried about screening results. Even low-risk grade health states (e.g., mild dysplasia and moderate dysplasia) are screened and diagnosed, and the risk of esophageal cancer is nearly 3-10 times higher than that of normal people. In this situation, patients may be scared of malignant deterioration and metastasis. Low-risk grade health states, such as moderate dysplasia, had ∼28 times higher esophageal cancer incidence than normal individuals (25). In this case, it is difficult for patients to accept the fact in a short time, which may act as a serious stressor and stimulation of life-stress events, especially for patients who have been screened for EC and precancerous lesions. Considering cancer progression, recurrence, and prognosis, they are prone to distress (26).
Although debate is ongoing, the microbiota-gut-brain axis is becoming as significant as the microbiota for monitoring bidirectional gut-brain communication pathways (27). Reliable evidence has demonstrated the relationship among brain cognitive function, mood, and intestinal flora (4,5). Reviews have shown that germ-free animals and animals with pathogenic bacterial infections played a vital role in the intestinal microbiota in the modulation of mood and cognition (27). Growing evidence indicates that the gastrointestinal microbiota is connected with anxiety and depression disorders. A wide range of studies consistently proposed that anxiety and depression impaired microbial characterization and diversity (8,9). The study found that subjects with anxiety and depression had lower alpha diversity, although no significant differences were observed. The results were consistent with findings from previous research that the diversity and abundance of intestinal flora in patients with anxiety and depression decreased overall (9,28). An important systematic review showed that α and β diversity were inconsistent. It indicated that the difference of bacterial taxa related to disorders may be manifested in a higher abundance of prion-flammatory species (e.g., Enterobacteriaceae) and lower bacteria producing short-chain fatty acid (e.g., Faecalibacterium) (10). Strong and consistent evidence has shown that Firmicutes and Bacteroidetes are dominant in human intestinal microbial flora (29). The change was inconspicuous in our study due to the relatively small sample size.
Several studies found consistent taxonomic differences among participants with generalizing anxiety disorders or depression relative to healthy controls, including higher Bacteroidetes and Proteobacteria and lower Firmicutes at the phylum level (10,28,30) and higher Prevotella and lower Faecalibacterium at the genus level (31), which is consistent with our results. Although these studies have found that the fecal flora of depressed patients is different from that of healthy individuals, the specific difference may vary, which may be related to the diagnostic criteria of research, inclusion criteria, and fecal detection methods (32). Animal experiments also found that the flora composition of depressed animals was similar to that of depressed patients, such as increased Bacteroidetes and decreased Firmicutes (33). However, several studies found the Lactobacillus and Bifidobacterium were decreased in depressed patients or animal models, which was not found in our study. It may be affected by the selection of subjects, inclusion criteria, and sample size.
In addition, we found that the specific genera Gemmiger, Ruminococcus, and Veillonella decreased in participants with anxiety and depression. Song et al. (34) proposed that a higher abundance of Bacteroides was linked with a higher fear of cancer recurrence. The relative abundance of Gemmiger in the anxiety and depression group was similar to that in another study related to diarrhea-predominant irritable bowel syndrome (2.4%) (35). In addition, Aranaz's study showed that the abundance of Gemmiger decreased in subjects with a higher inflammatory index (36). Similar results were found in a systematic review in which a reduced abundance of Ruminococcus was observed in depressed people (37). Evidence has shown that the high abundance of Parabacteroide, Oscillibacter, Paraprevotella, Veillonella, Klebsiella, and Desulfovibrio in patients with depression may demonstrate the role of flora in the emergence of depression (38). These results indicated that Gemmiger, Ruminococcus, and Veillonella may be the characteristic and specific genus of high-risk groups and vulnerable participants of anxiety and depression.
A growing body of evidence intriguingly suggests that the microbiota composition of individuals may affect their susceptibility to anxiety and depression (27). A key study presented the role of mouse microbiota transplantation in detecting the microbiota-gut-brain axis (39). A landmark study showed that sterile mice changed the function of the hypothalamic-pituitary-adrenal axis (HPA), which can be reversed by inhabiting specific bacterial strains early in life (40). The mechanism of how microbiota influences gut-brain signaling may be associated with the pathophysiology of anxiety and depression by delivering peripheral inflammation to the central nerve (28). These mechanisms may include modulating microbial composition, activating immunity, transducing vagal signals, alternating tryptophan metabolism, and producing specific microbial neuroactive metabolites (30).
In fact, due to the lack of longitudinal investigation in this study, we do not know how long psychological distressassociated gut microbiota would persist, and longitudinal investigations are sparse. First, we measured the symptoms of anxiety and depression with the GAD-7 and PHQ-9 in the past 2 weeks. Second, the intestinal flora was greatly influenced by diet, lifestyle, geography, and age, and the composition was dynamic and fluctuating (41,42). Largescale population studies found that antibiotics used in antiinfective treatment significantly increased the risk of individual psychological diseases such as anxiety and depression (43,44). Once the microbiota becomes unbalanced, alteration may occur to the microenvironment and then lead to gastrointestinal diseases and even cancer (6,7). Experimental evidence has shown that the human gut flora affects the occurrence and progression of gastrointestinal tumors by activating carcinogenic signaling pathways, producing tumor-promoting metabolites and inhibiting antitumor immune responses (45). In addition, evidence from saliva and tissues showed that oral flora may act as potential risk factors for oral and gastrointestinal cancer (7,46), but there were differences in flora changes among various studies. Fusobacterium nucleatum mainly inhabits the oral cavity and causes periodontal disease, which may promote the aggressive behaviors of tumors by activating chemokines (e.g., CCL20) in esophageal cancer tissues (47). Therefore, given the important role of the gut microbiome in maintaining homeostasis, a better understanding of the microbiome in cancer screeners is increasingly important. We can develop innovative cancer prevention and therapeutic strategies by targeting the gut microbiota. The human intestinal microbiome plays an important role in cancer screening, especially for gastrointestinal cancer (48). Meta-analysis indicated that fecal bacteria and oral flora act as promising biomarkers for the noninvasive diagnosis of gastrointestinal cancer (49).
Since the gastrointestinal microbiota is altered through the rational use of prebiotics, probiotics, and antibiotics (50), the relationship between mental disorders, the microbiota, and tumors may be of clinical significance and implications. The novel concept of the microbiota-gut-brain axis indicated that regulation of the intestinal flora may be a feasible strategy to develop innovative therapeutics for psychological distress.
These mechanisms may clear the way for microbial-based psychotherapies. Mind-altering microorganisms refer to the gut microbiota that could alter the brain and behavior. An important study showed that the potential probiotic could regulate behaviors related to anxiety and depressive disorders and alter central levels of γ-aminobutyric acid receptors (51). Compared with the placebo control group, probiotic Bifidobacterium longum and Lactobacillus helveticus reduced depression scores and altered the brain activity of patients with anxiety and depression (52). An important study consistently confirmed that probiotics help prevent and relieve depression disorders (31). In addition, clinical studies found that taking prebiotics daily for 3 weeks reduced the activation of the cortex caused by negative information, thereby reducing anxiety-like and depressive-like behaviors (53).
The pattern of intestinal microbiota changed significantly in patients with anxiety and depression. Further evidence is needed to translate microbiome findings into innovative clinical treatments to improve therapeutic effects in patients with mental disorders. Anxiety and depression have the characteristics of a low treatment rate, poor compliance, and recurrence, and this study provides a new promising research direction to improve psychological distress, namely mindaltering microorganisms. This suggests that we can attempt to explore probiotics, prebiotics, and other microecological agents to regulate the balance and homeostasis of intestinal flora in cancer screening progress, which contributes to optimizing screening and maximizing the net benefits of cancer screening. However, little clinical and large-scale research on probiotics and prebiotics has been used in treatment strategies. To date, the selection of probiotics is relatively random and difficult to predict, and it is hard to show stable efficacy in different groups. There are some shortcomings in the existing studies, such as a small sample size and poor contrast in the use of probiotics. Therefore, more randomized controlled trials are needed to further verify the efficacy of probiotics and prebiotics.
Adequate evidence has shown the bidirectional relationship between cancer and intestinal microbiota (54), such as colorectal cancer and malignant gastrointestinal diseases. Evidence highlighted the multifaceted role of the intestinal microbiota in cancer. The occurrence of cancer is usually accompanied by inflammation and leads to microbial alteration and disorder of the intestinal microbial environment, such as increased abundance of Escherichia coli and Fusobacterium nucleatum in colorectal cancer (54,55). Conversely, the unbalanced microbiota could also cause a proinflammatory microenvironment and DNA damage, increasing the risk of cancer and deteriorating the prognosis, such as Helicobacter pylori, and invasive Escherichia coli. Further studies are needed to decipher whether there is a synergistic effect of microbial and psychological distress in tumor promotion, which may be exploited therapeutically in the future.
To our knowledge, first, this is the first study to explore the microbial characterization of participants in gastrointestinal cancer screening, filling the gap in the field of screening. This perspective is innovative and provides new ideas and optimization strategies for the comprehensive evaluation of cancer screening. Second, the included subjects were from the same regions, with similar dietary patterns, and the controls were matched by age and sex, which makes the findings more objective and credible to some extent. Some limitations of our work should be acknowledged. First, the stratified analysis of anxiety and depression separately was confined by the relatively small sample size. This may lead to no significant difference on some microbial characterizations and diversity, but the trend of our results provided scientific references in this field. Our exploratory study largely provides clues from a novel perspective. This study is an attempt to explore the microbial characterizations of screeners. Larger studies will be needed to reduce the uncertainty and confirm our associations. Second, we did not have specimens from nonscreeners as bank controls. Third, selection bias may exist, and the sample may not be representative of all screeners. Further multicenter, large-scale, and prospective cohort studies, randomized controlled trials, and clinical trials are needed to validate the results. Fourth, anxiety and depression symptoms were evaluated in the study, not a clinical diagnosis of anxiety and depression disorders. In the future, psychiatrists could be considered in the screening process. Finally, considering the goal of this study and other factors (e.g., economy, efficiency, and data processing), we chose 16S rRNA gene sequencing. In the future, full metagenomics could be applied for further exploration of mechanisms, pathway in-depth, and functional prediction analysis.

CONCLUSIONS
Differing microbial characterization among participants with anxiety and depression was found in the endoscopic screening of upper gastrointestinal cancer in China. Gemmiger, Ruminococcus, and Veillonella were informative for psychological distress in cancer screening and have potential clinical implications for mental disorders, which provides references for optimizing cancer screening and minimizing psychological harm. The results should be explained cautiously, and more large-scale, prospective cohort studies are needed in the future to validate the results and further explore biological mechanisms and the relationship among gut microbiota, psychological distress, and cancer risk in cancer screening.

DATA AVAILABILITY STATEMENT
The data presented in the study are deposited in the the Genome Sequence Archive (GSA) repository, accession number CRA005126.