Application of Multigroup Technology in Non-Small-Cell Lung Cancer with Qi Stagnation and Blood Stasis Syndrome

Objective To explore the basic characteristics of intestinal flora, metabolomics, and proteomics of non-small cell lung cancer (NSCLC) in patients with Qi stagnation and blood stasis syndrome. Methods Twelve NSCLC patients with Qi stagnation and blood stasis syndrome were selected for the QZXY group and 15 healthy volunteers were selected for the control group. Fecal samples from the two groups were collected to evaluate intestinal microecology using the 16s rDNA technique. Serum samples were collected to compare the differences in metabolomics and proteomics between the two groups using liquid chromatography-mass spectrometry (LC-MS). Another 34 NSCLC patients with other syndromes were selected for the nQZXY group and their serum samples were collected. Metabolomics differences between the QZXY and nQZXY groups were compared using LC-MS, and four metabolites with the most obvious differences were selected for receiver operation characteristic curve representation. Finally, multigroup results were analyzed using the WGCNA software. Results There were two significantly different types of bacteria (Aerococcaceae and Abiotrophia), 11 different proteins (six upregulated and five downregulated), and 38 different metabolites (nine upregulated, 29 downregulated) between the QZXY and control groups. There was a correlation between differential bacteria, proteins, and metabolites. The conjoint analysis found that the different substances were related to MAPK, PI3K/Akt, Ras signaling pathway, cancer pathways, and cytokine-cytokine receptor interaction. There were four significant differences in metabolites (Pseudouridine, phenlacetyl-C0A, L-glutamic, and phospho-anandamide) between the QZXY and nQZXY groups. Conclusions NSCLC with Qi stagnation and blood stasis syndrome had specific intestinal flora and protein and metabolites, which were closely related to the occurrence and development of tumors.


Introduction
Lung cancer ranks first in the incidence and mortality rates of malignant tumors in China. Many studies [1][2][3] have shown that the application of traditional Chinese medicine (TCM) can control tumor progression, improve quality of life, and prolong survival time during treatment of nonsmall cell lung cancer (NSCLC). TCM may be a strong therapeutic candidate for tumors. Syndrome differentiation is one of the important principles of TCM targeted at recognizing and dealing with diseases. It is based on a comprehensive analysis of clinical information collected using the main TCM diagnosis procedures: observation, listening, questioning, and pulse analyses. In recent years, research on the standardization, objectification, and quantification of TCM syndromes has become popular.
Systems biology is an emerging discipline based on the development of modern and continuous biological experiment technology. Its research methods include transcriptomics, proteomics, metabolomics, and intestinal microecology. e concept of the system is similar to the holistic concept of TCM. TCM syndromes describe the functional state of the body at a certain stage as part of TCM theory. It is a high level generalization of the physiological and pathological state at a certain stage during the occurrence and evolution of the disease. e syndromes can be classified as instantaneous, dynamic, and integral.
Using systems biology methods to objectively study TCM syndromes is conducive to revealing their essence. e TCM pathological properties of lung cancer are as follows: deficiency in origin and excess in excess, the healthy Qi hurts first, and the evil Qi takes the opportunity to enter the lung. Qi stagnation and blood stasis are important pathogenesis processes of lung cancer [4] and [5], and the prescription of regulating Qi and activating blood has a good curative effect [6]. In this study, 16S rDNA and liquid chromatographymass spectrometry (LC-MS) technology were used to study the proteomics, metabolomics, and intestinal microecology of NSCLC with Qi stagnation and blood stasis syndrome, aiming to provide a biological basis for the diagnosis of TCM syndromes in NSCLC, and an experimental basis for the future of anticancer treatment using Chinese medicine.

General Information.
e subjects of this study were 46 patients from the thoracic tumor surgery outpatient department of Zhejiang Cancer Hospital between March 1, 2019, and September 30, 2019, who planned to undergo radical resection of lung cancer without any treatment before enrollment. Twelve patients were identified to have the QZXY syndrome by two vice directors or TCM physicians at the same time. e remaining 34 patients were identified as participants with other syndromes (nine patients had spleen deficiency and dampness stagnation (PXSZ) syndrome, 11 patients had a deficiency of Qi and Yin (QYLX) syndrome, 14 patients had Yin deficiency with internal heat (YXNR) syndrome). In addition, 15 healthy volunteers were included.
is clinical study was approved by the Ethics Committee of Zhejiang Cancer Hospital (Approval no. IRB-2018-219).

Diagnostic Criteria of the Patients
e patients were confirmed to have NSCLC by histopathology and/or cytopathology [7]. e patients met the following diagnostic criteria of the QZXY syndrome [8,9], which included a cough that was not smooth, shortness of breath, chest pain as cone or thorns, dry mouth, constipation, hemoptysis or phlegm that is bloody and dark red, lip dark purple tongue, tongue ecchymosis, thin yellow moss, and pulse string or astringent.

NSCLC Patients
Patients meet the diagnostic criteria Patients had no history of other malignancies Patient ages ranged from 18 to 75 years Patients signed the informed consent form in person

Volunteers
Volunteers had no history of malignancies or major diseases, such as cardiovascular, cerebrovascular, liver, kidney, and endocrine diseases Volunteer ages ranged from 18 to 75 years Volunteers signed the informed consent form in person

Exclusion Criteria of Patients
Patients were mentally ill and incapacitated Patients or volunteers had been treated with antibiotics for nearly 3 weeks

Intestinal Microecological Analysis.
Fecal samples (200 mg each) were mixed with 30 mL of buffer solution. After centrifugation at 1,000 rpm for 5 min, the precipitation was collected and resuspended in 10 mL of buffer solution. Fecal microbial DNA was extracted using the QIAAMP Fast DNA Stool Mini Kit. After extraction, the total DNA was inspected using agarose gel electrophoresis and ermo Nanodrop 2000 UV microspectrophotometer. e samples were amplified in 16S V3-V4 region. e primer information was as follows: forward primer: CCTACGGGNGGCWGCAG; reverse primer: GACTACHVGGGTATCTAATCC. e diluted genomic DNA was used as a template for PCR using Phanta Max Master Mix (Vazyme) high fidelity enzyme. e reaction conditions of PCR were as follows: predenaturation at 95°C for 3°min; 95°C denaturation for 30 s, 55°C annealing for 30s, 72 extension for 45s, and repeat 25 cycles. e final incubation was carried out at 72°C for 5 min. e PCR reaction system included the following: 25 μL of 2 × Phanta Max Master Mix (Vazyme), 2 μL of each primer (10 μM), and DDH2O added to make up a total of 50 μL. PCR product (1.5 μL) was analyzed using 120 V continuous electrophoresis for 20 min in 2% agarose gel. Subsequent UV imaging was performed using the gel imaging system.
After the quality inspection of the library, Qubit was used to quantify the concentration of the mixed library pool, Mixing was carried out according to the volume requirements of each sample. Finally, the Miseq sequencing program was run. Clean reads with the same sequence were sorted according to their abundance, while singletons were filtered out. USearch was used for clustering at a similarity degree of 0.97. After chimeric filtering of the clustered e concentrate was digested in a wet chamber at 37°C for 12 h with trypsin (enzymatic to protein ratio: 1 : 100, 50 M MABC) and then centrifuged. e digestive solution was collected and rinsed with 50 μL of 0.5 M sodium chloride. e resulting solution was combined with 10% trifluoroacetic acid, acidified, and then desalted using C18Tips (Pierce, ermo Science Titi C). Activation, equalization, peptide elution, and elution were performed according to the manufacturer's instructions. TIPS were activated with two rinses in 50% acetonitrile (ACN) and equilibrated with two rinses in 0.1% trifluoroacetic acid (TFA) in water. e peptide-rich solution was slowly aspirated and distributed using TIP for 10 cycles. Two washing desalination processes were performed with 0.1% TFA. e desalinated peptide was eluted with 50 μL of 50% can. e eluted product was dried in a vacuum and then suspended in 30 μL of 0.1% formic acid. e solution containing the desalinized peptide (3 μL) was separated using an EASY-Nano-LC system ( ermo Fisher Scientific) and online nanomobile liquid chromatography. e LC was connected to a 2 cm precolumn with an inner diameter of 100 μm and filled with 5 μm of C18 resin ( ermo Fisher Scientific). After the precolumn reaction, the 75 μm × 15 cm capillary column was filled with 3 μm of C18 resin ( ermo Fisher Scientific). e polypeptides were eluted from the column front and analyzed on the column using step-by-step gradient elution, followed by a search of all the original files against the UniProt human protein sequence database in MaxQuant (version 1.6). Perseus was used to calculate the folding changes of each subgroup. Proteomes with more than twofold changes in the data frame (>2 or < 0.5) were labeled as significantly altered proteins.

Combined Analysis of Intestinal Microecology and
Metabolomics. Spearman's correlation analysis was performed on the secondary differential metabolites screened using metabolomics and significant differential level bacteria obtained via 16S sequencing analysis in order to identify the relationship between microflora and metabolites. Based on the calculation results, appropriate screening conditions were selected to obtain the final correlation and network diagram.

Combined Proteomics-Metabolomics Analysis.
By combining metabolomics and proteomics data using the KEGG pathway analysis, proteins and metabolites that are involved in significant changes in the same biological process (KEGG pathway) can be identified and used to quickly target key genes. is process was divided into the following steps: (1) check whether proteome and metabolome are related through KEGG metabolic pathway; (2) data screening for significant difference data and regulatory relationship; (3) GO and KEGG enrichment analysis; (4) clustering of expression patterns for linked data; and (5) creation of linked data network diagram.

Statistical Method.
e intestinal floras with significant differences among different groups were identified using a rank-sum test. e threshold value for significance screening was a P value < 0.05. All p values were verified using FDR. In the analysis of differential proteins, the differentially expressed protein needed to meet the following conditions: the coincidence difference ratio was 1.2 times (up-down) and P value < 0.05. e t-test method was limited to the comparison between two groups of samples with biological duplicates. When conducting differential metabolite analysis between two groups of samples, univariate analysis can intuitively show the significance of metabolite changes between two samples, thus helping to screen for potential marker metabolites (usually FC > 1.5 and P value < 0.05 as the screening criteria). Spearman's correlation coefficient was used for joint analysis. It is a linear correlation coefficient, which is a statistical quantity used to reflect the degree of linear correlation between two variables. e correlation coefficient is expressed using r, where n is the sample size. It is the observed value and mean value of the two variables, respectively. R describes the degree of linear correlation between two variables. e greater the absolute value of R, the stronger the correlation. e correlation analysis of Evidence-Based Complementary and Alternative Medicine significantly different metabolites between the QZXY and nQZXY groups was carried out using the binary logistic regression stepwise method. A receiver operating characteristic curve was used to evaluate the efficacy of differential metabolites in differentiating the two.

Baseline Subject Data.
Twelve NSCLC patients with QZXY syndrome were included in the QZXY group. A total of 34 NSCLC patients with other syndromes (PXSZ, QYLX, and YXNR) were included in the nQZXY group. Fifteen healthy volunteers formed the control group. e above information is presented in Table 1.

Differences between Microecological Groups
e 16 S sequencing analysis showed that a total of 249 OTUs were generated from 27 samples ( Figure 1).

Differential Species Analysis.
Linear discriminant analysis (LDA) effect size (LEFSE) analysis used linear discriminant analysis (LDA) to estimate the influence of the abundance of each component (species) on the difference effect and to determine the communities or species that had a significantly different influence on sample division. e LDA scores obtained using linear regression analysis for the significant microbial groups are shown in Figure 2.

Analysis of Spearman's Correlation Coefficient in
Different Species. LEFSE was used at each level or the ranksum test was used at the genus level (or at a specific level) to select the species with different abundance TOP30. Spearman's correlation heat map for species was created using the Corrplot package of R software. e specific results are shown in Figure 3.

Statistical Analysis of Upregulation and Downregulation of Differential Proteins.
e up-and downregulation frequency statistics for differential proteins were used to determine the number of differential proteins under different experimental conditions. e horizontal axis represents the comparison group information, and the vertical axis represents the number of proteins. Red identifies the upregulated proteins, blue shows the downregulated proteins, and numbers represent the up-and downregulated proteins ( Figure 4).

Functional Analysis of Differential Protein Genes.
e term is the basic unit of GO analysis, and each term corresponds to an attribute. e GO significance enrichment analysis first mapped all of the differentially expressed genes to each term in the GO database, calculated the number of genes in each term, and then applied a hypergeometric test to identify the GO terms significantly enriched in the differentially expressed genes compared to the entire genome background. e results are shown in Figure 5.

KEGG Pathway Analysis of Differential Protein.
e KEGG pathway was used as a unit of pathway enrichment analysis. e hypergeometric test was used to identify the pathways that were significantly enriched in the differentially expressed genes compared to the overall genome background. e results are shown in Figure 6. Statistical analysis indicated that the baseline data were comparable. e abundance level in the QZXY group was significantly lower than that in the control group. 4 Evidence-Based Complementary and Alternative Medicine

Differences between Metabolomics
Groups. In this study, univariate statistical analysis and multivariate statistical analysis were used to screen for significant difference metabolites. MetDNA and OSI/SMMS software were used for the structure identification of metabolites with significant differences. e specific identification results are shown in Table 2 and Figure 7. e obtained differential metabolites were submitted to the Metaboanalyst 4.0 website for metabolic pathway analysis. e results are shown in Table 3.      Figures 8 and 9, where red identifies a positive correlation and blue signifies a negative correlation. Figures 8 and 9 show that Peptoniphilus was positively correlated with phenylacetyl-L-glutamine and indoxyl sulfate. Prevotella was negatively correlated with phenylacetyl-L-glutamine and 16-hydroxypalmitate. Bifidobacteria and 18-hydroxyoleate were also negatively correlated.

Differential Metabolites and Microflora Network
Regulation Analysis. Different nodes in the figure represent different flora (round) or metabolites (triangular). e line between the flora and the metabolites represents the correlation between the two, where a solid line identifies a positive correlation and a dashed line represents a negative correlation. e results are shown in Figures 10 and 11.

Combined Proteomics-Metabolomics Analysis.
e results of the combined proteomics-metabolite analysis are shown in Figure 12.

Metabolite Analysis of Difference between the QZXY
Group and the nQZXY Group. ROC curves were used to analyze the different metabolites between the two groups. e results are shown in Figure 13.

Discussion
At present, the NSCLC treatment mainly involves different comprehensive treatment strategies according to pathological tumor types and clinical stages, which can somewhat improve the clinical symptoms and prolong survival time. However, individual differences lead to great variations in prognosis [10]. Increasing evidence [11,12] has indicated that TCM with syndrome differentiation can improve the current situation. However, its validity has always been questioned because of the undefined molecular mechanism of syndromes. e present report attempted to provide a scientific explanation for the syndrome by observing the biological basis for the syndrome's type of Qi stagnation and blood stasis in NSCLC patients. e present study found that peptone bacteria were positively correlated with glutamine, indolophenol, and other metabolites through network regulation analysis of differential metabolites and flora. e upregulated metabolite glutamine was found to inhibit the migration of human lung adenocarcinoma A549 cells, possibly by downregulating the expression of TNF-α and NF-ΚB p65 protein   Figure 6: Pathways most associated with outcomes in the QZXY group included MAPK signaling pathway, Phagosome, Ras signaling pathway, and cytokine-cytokine receptor interaction. Tumor-related pathways were also analyzed.     Evidence-Based Complementary and Alternative Medicine [13]. Indophenol can be produced by isopentenylation of tryptophan derivatives. Tryptophan has been found to be metabolized rapidly in many tumors. Peptoneophilus is positively correlated with metabolites, such as glutamine and indophenol. It has been found to be enriched in endometrial carcinoma tissues [14]. No correlation has been found between peptones and lung cancer. is suggests that peptones may affect the formation and development of lung cancer by regulating a variety of metabolic pathways. Prevotella and Bifidobacteria were negatively correlated with lipid metabolites (palmitic acid, oleate ester). Palmitic acid can promote the epithelial-mesenchymal transformation of lung tumor cells and thus enhance their invasion. Microflora has been shown to be involved in host regulation of metabolism and in the production of hormones and bile acids [15]. e present study found that many bacterial genera were negatively correlated with lipid metabolism. is suggested that there may be a certain mechanism between bacterial flora and lipid metabolism to regulate the occurrence and development of lung cancer. NSCLC patients with QZXY * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 5a-Pregnane-3,20-dione syndrome often exhibit blood stasis signs, such as stasis spots on the tip of the tongue. Studies [16] have shown that application of Danshensu, an extract of Salvia miltiorrhiza with the effect of promoting blood circulation and removing blood stasis, can effectively regulate dyslipidemia in hyperlipidemia rats, providing a theoretical basis for the correlation between abnormal lipid metabolism and QZXY syndrome. e present study explored the changes in signaling pathways in NSCLC patients with QZXY at the multiomics level. Five signaling pathways, including MAPK, PI3K-AKT, Ras, cancer pathways, and cytokine-cytokine receptor interaction, were screened using the linked data network diagram. As a traditional Chinese medicine for promoting blood circulation and removing blood stasis, the active ingredients in Panax notoginseng can act on EGFR, MAPK14, and other targets, regulating cancer pathways, cytokine-cytokine receptor interactions, and other signaling pathways and exert anti-inflammatory effects [17]. It was also found that corydus can exert the effect of promoting blood circulation by regulating the cytokine-cytokine receptor interaction and other signaling pathways [18]. e combination of Astragalus and Zedoary (toning Qi and activating blood circulation) can downregulate the expression of MAPK and inhibit the growth and metastasis of transplanted Lewis lung cancer tumors in tumor-bearing mice [19]. To summarize, we speculated that the above five pathways were closely related to lung cancer. Further studies may help to explore how to use TCM theories combined with biomarkers to deeply explore the mechanism of TCM treatment for lung cancer, which is particularly important and can provide new targets for clinical treatment.
We also attempted to analyze the metabolomic differences between the QZXY and nQZXY groups. Four metabolites with ROC > 0.9 and high sensitivity and specificity were found. Pseudouridine is a component of tRNA synthesized by the human body, which can significantly reduce  Figure 10: Results of differential metabolite and flora network regulation analysis in anion mode show that Abiotrophia and Peptoniphilus are closely related to many metabolites. is suggests that they may be the key bacteria in the development of lung cancer. the ability of RNA to stimulate the immune system. It may help lung cancer cells to avoid an immune response. L-Glutamine has also been analyzed. e study results indicated that compared to NSCLC patients with other syndromes, there were significant differences in metabolites in patients with QZXY syndrome.
is meant that QZXY syndrome could be diagnosed and differentiated from other syndromes, suggesting that there was microscopic evidence of its presence. However, the sample size in the present study was small, potentially creating a source of bias. erefore, these results only provided a limited reference value for future investigations. One prior study [20] has found that the syndromic types had different effects on tumor progression by establishing the "Zheng-first" and "Tumor-first" models of subcutaneous pancreatic cancer in mice with different syndromic types. ese results may provide a direction for follow-up studies. e study of TCM syndromes using systems biology methods is helpful to provide an objective basis for syndrome differentiation and treatment and to suggest highly effective and low-toxicity individualized treatment for tumor patients. In this study, combined multiomics analysis was conducted in NSCLC patients with QZXY.
is investigation not only promoted the study of lung cancer mechanisms and determine its pathogenic targets, but also provided theoretical support and technical foundation for the prevention and treatment of lung cancer using traditional Chinese medicine. Evidence-Based Complementary and Alternative Medicine 13

Conclusions
NSCLC with Qi stagnation and blood stasis syndrome had specific intestinal flora and protein and metabolites, which were closely related to the occurrence and development of tumors.

Data Availability
e data used to support the findings of this study are available from the author upon request.

Conflicts of Interest
e authors declare no competing interests.