Fecal gene detection based on next generation sequencing for colorectal cancer diagnosis

Si-Yu He, Ying-Chun Li, Yong Wang, Hai-Lin Peng, Cheng-Lin Zhou, Chuan-Meng Zhang, Sheng-Lan Chen,Jian-Feng Yin,Mei Lin

Abstract

Key Words: Colorectal cancer; Feces; Next generation sequencing; Diagnosis; Gene

INTRODUCTION

Colorectal cancer (CRC), as a health issue of great concern worldwide, poses a serious threat to human health. According to the latest data released by the International Agency for Research on Cancer, there were approximately 19.3 million new cancer cases and nearly 10 million cancer deaths globally in 2020.Among them, new CRC cases amounted to about 1.93 million (10%), ranking third, while CRC deaths amounted to about 0.94 million (9.4%), ranking second. For female patients, there were approximately 860000 new cases of CRC (9.4%), second only to breast cancer, and nearly 870000 CRC deaths (9.5%),ranking third. As for male population, the numbers of new cases and deaths both ranked third for CRC,which were 930000 (10.6%) and 50000 (9.3%), respectively[1].

In recent years, the morbidity and mortality of CRC have been on the rise in China, and the onset age of patients has gradually become younger. Early detection of CRC is difficult due to the insidious onset,long progressive course, and unobvious symptoms and signs of patients at the onset stage. As the tumor grows and gradually compresses the intestinal cavity, such symptoms as altered bowel habit and abnormal stool color or traits may occur. At this time, the condition often already progresses to advanced stage. Studies have shown that the 5-year survival of patients with early CRC can reach 90%,while is only about 12.5% for advanced stage patients[2]. Thus, early diagnosis is conducive to early intervention and treatment of CRC, which is of great significance for improving the efficacy and prognosis[3,4]. In response, the American Cancer Society recommended the CRC screening since the age of 45 for adults in 2018, which is 5 years earlier than the previously recommended age[5].

The occurrence and progression are a multi-factor, multi-stage complex process[6-8]. It is generally believed that CRC is caused by the mutations and malignant proliferation of colorectal mucosal epithelial cells. Other studies have shown that CRC can develop from inflammatory bowel disease(IBD). IBD is mainly characterized by intestinal mucosal barrier dysfunction, because the changes of cytokines aggravate tissue damage and eventually lead to tumorigenesis. It is reported that patients with ulcerative colitis have a higher risk of developing CRC than the general population[9]. In addition,microRNA and tumor microenvironment also have potential pathogenic effects on CRC[10,11].Currently, the clinically recognized gold standard for CRC detection remains endoscopy combined with histopathological analysis. Given the heterogeneity of tumor tissues and the invasiveness of endoscopy,the same patient cannot be sampled multiple times, so that the genetic testing results cannot fully display the mutation spectrum of oncogenic genes[12]. According to the ChineseExpert Consensus on Early Diagnosis and Screening Strategies for Colorectal Tumorspublished in October 2018, fecal detection is a promising method for early diagnosis[13]. Its theoretical basis is that the tumor cells in the intestinal cavity of CRC patients adhere less to the base and can be continuously exfoliated from the colorectal mucosa. Studies have shown that tumor-associated DNA can stably exist in the stools, which has high specificity for the detection of CRC and is the most suitable biomarker in fecal detection[14-16]. By extracting the DNA of shed tumor cells in stools and analyzing their mutations, the traits of intestinal tumor cells can be reflected. Fecal DNA detection, as an emerging, noninvasive technique for screening CRC, has attracted great attention in recent years, which has the advantages of noninvasiveness,convenient access to materials, good patient compliance, and continuous dynamic monitoring. It is in line with the concept and conforms to the trend of "precision medicine", which is an exceptionally promising detection method.

With the continuous development of molecular biotechnology, next generation sequencing (NGS) has emerged to gradually become a hot research topic worldwide. It is also known as high-throughput sequencing, which operates on the following principles: Four different deoxynucleotide triphosphates(dNTPs) are labeled with different fluorescent colors. Conforming to the complementary base pairing rule, the bound base releases the corresponding fluorescence by the action of DNA polymerase every time a dNTP is added. Finally, the base is identified according to the fluorescent signal. Compared to the first-generation Sanger sequencing, NGS is a set of technologies that completely change the standard concept of nucleic acid sequencing. Capable of detecting millions of DNA molecules simultaneously, it achieves concurrent synthesis and sequencing, which can greatly shorten the detection time, and has the advantages of high throughput, high efficiency, and rich information[17,18]. NGS can more comprehensively display the full spectrum of tumor-associated genes, which can reflect the changes in tumorspecific genes, expand the understanding of tumor-specific genes, and improve the predictability of targeted therapy, suggesting its significance in achieving personalized tumor treatment and developing new therapeutic strategies[19-21]. Currently, NGS has been widely applied in the research of various solid tumors, including CRC. Kraus and his colleagues performed NGS analysis on the tissue paraffin sections and blood samples of 152 CRC patients based on 18 CRC-related genes. Their results showed that NGS had a high detection rate and could effectively identify the high-risk patients with hereditary CRC[22]. In a study by Liet al[23], the NGS (Ion Torrent PGM platform) technology was applied to analyze 22 tumor-associated mutation hotspots in 526 CRC patients, finding thatTP53andKRASmutations were the most common in tumor tissues, and confirming that NGS was a highly sensitive detection technique.

This study performed NGS-based genetic testing of 50 tumor-associated mutation hotspots in the stools of CRC patients, and compared them with the genetic testing results of corresponding tumor tissues, with an aim to find the fecal genes or gene combinations with high specificity and sensitivity,and to establish a technological platform for CRC screening, diagnosis, and efficacy monitoring through fecal DNA detection. The findings herein will provide support and basis for personalized diagnosis and treatment of CRC.

MATERIALS AND METHODS

Subjects

Fifty tumor tissues, 27 preoperative stools, and 19 postoperative stools (stools formed over 7 d postoperatively) collected from patients with CRC were included in an experimental group. Twenty stool samples collected from healthy subjects were included in a normal control group. Seventeen stool samples of patients with benign intestinal polyps confirmed by digestive endoscopy were collected, as well as three stool samples from patients with ulcerative colitis, totaling 20 samples, were included in a benign control group. Postoperative pathological reports of polyp patients revealed inflammatory or adenomatous polyps, none of which had progressed to the CRC stage.

This study was approved by the Ethics Committee of Taizhou People's Hospital in Jiangsu, and the patients provided signed informed consent for sample collection.

Main reagents

Fecal DNA sample preservation tubes, tissue nucleic acid preservation tubes, fecal DNA extraction kit(immunomagnetic bead method), universal columnar genomic nucleic acid extraction kit, fast DNA library preparation kit were provided by Cowin Bio., Jiangsu, while high-throughput kit (MGISEQ-2000RS) was provided by MGI, Shenzhen.

Samples sources

The subjects of this project are patients who were diagnosed with CRC for the first time and underwent surgery at Taizhou People's Hospital from January 2019 to January 2021. A total of 54 tumor tissue samples, 45 preoperative stool samples, and 41 postoperative stool samples were collected. The normal control group consisted of healthy physical examination takers without malignancies (including CRC),intestinal polyps, or IBD, from whom 20 stool samples were collected. The benign control group comprised 20 patients who were diagnosed with intestinal polyps or IBD by digestive endoscopy, and admitted to the Department of Gastroenterology, Taizhou People's Hospital for polypectomy or antiinflammatory treatment. Twenty formed stool samples were collected from patients who had not taken laxatives preoperatively or from patients before receiving anti-inflammatory therapies. The above samples were tested by NGS, unreliable results were filtered during analysis of genetic testing results based on biological information, and poor quality samples were eliminated (heterogeneous tumor tissues, atypical fecal specimens, and little or degraded DNA in stools). Finally, effectively detected samples included 50 tumor tissues, 27 preoperative stools, 19 postoperative stools, 20 stools in the normal control group, and 20 stools in the benign control group.

Sample collection and processing

Stool samples: Patients were instructed to use the fecal DNA preservation tube correctly (tubes contained protective solution to prevent the degradation of nucleic acid, while steel beads played a mixing role). They were asked to take samples immediately after natural defecation. Using sampling spoon located on the tube cap, typical formed stool about the size of soybean was excavated while avoiding spillover of protective solution and steel beads in the tube. After tightening the tube cap, the patients were asked to shake the tubes for approximately 30 s for uniform mixing. The samples were stored in a -80 °C refrigerator for subsequent testing.

Tissue samples: The freshly isolated tumor tissues were quickly cut into pieces (< 0.5 cm in thickness)using sterile tissue scissors, and then the tissue pieces were soaked in the protective solution of preservation tubes (if the volume was excessively large, the protective solution would not easily penetrate the tissues, which could not protect the nucleic acid from degradation). After tightening the tube caps, the samples were stored in a 4–8 °C refrigerator for subsequent testing.

DNA extraction from stool samples

The DNA in stool samples should be extracted strictly in accordance with the kit instructions (immunomagnetic bead method). Initially, the samples were homogenized, and then the precipitates (fecal particles) were removed, and the supernatants were retained for subsequent use. Different reagents were added to the corresponding wells of deep 96-well plates asperthe instruction procedures, and then the deep well plates and magnetic sleeves containing the reagents were placed at the corresponding positions of nucleic acid extractor, followed by running of the stool extraction program. Under high salinity condition, DNA can bind to the surfaces of silica-coated Magbeads. After rinsing, highpurity DNA was eluted and stored in Buffer GE. The DNA-containing Buffer GE was transferred to an EP tube and stored under refrigeration for subsequent testing.

Table 1 details the reaction system. Relevant reagents were added to the corresponding wells of deep 96-well plates (wells 1–6 corresponded to one sample, wells 7–12 corresponded to another sample, and so forth).

The deep 96-well plates and corresponding magnetic sleeves were placed into the corresponding positions of an automatic nucleic acid extractor, and then the fecal DNA extraction program was initiated. After running for about 20 min, the program was paused, and each 220 μL of pre-prepared mixture (Magbeads PN 20 μL, isopropanol 200 μL) was added to the first well of each sample, followed by continuation of the program. The product in the last well of each sample was transferred to the centrifuge tube and stored at low temperature for subsequent testing.

DNA extraction from tissue samples

DNA was extracted from tissue samples strictly following the instructions of universal columnar genomic nucleic acid extraction kit. Initially, fresh tumor tissues were ground to small pieces and added with lysis buffer and proteinase K for complete digestion and lysis. Then, RNase was added to remove RNA. Since DNA was insoluble in absolute ethanol, the samples added with absolute ethanol were centrifuged to adsorb the DNA on the column. The residual impurities dissolved in absolute ethanol were discarded along with the liquid waste in the collection tube. Next, the collection tube was replaced with a new one, and the centrifugation was repeated several times. The DNA-containing adsorption column was placed at room temperature, and let stand to dry the absolute ethanol. The elution buffer was suspended, added to the adsorption column with a new collection tube, and centrifuged, and then the DNA-containing eluate in the collection tube was collected and refrigerated for subsequent testing.

The nucleic acid extraction was carried out in Cowin Bio., Jiangsu.

Library construction and NGS detection

In accordance with the instructions of fast DNA library preparation kit, the hotspot mutation panels of 50 tumor-associated genes designed for specific genomic regions or targets were added initially based on the multiplex polymerase chain reaction (PCR) library construction technology, and then multiplex PCR enrichment was performed on specific gene region fragments. Next, the obtained target gene fragments were subjected to end repair, barcode adapter ligation, and PCR amplification (Eastwin Scientific Equipment, Suzhou) to prepare a library. Afterwards, the library was quantified with Qubit 3.0 fluorescence quantifier (Thermo Fisher Scientific, Shanghai), and the fragment size of the library was analyzed with the Bioptic Qsep 100 automatic nucleic acid analyzer (Bioptic Inc., Jiangsu) to evaluate the library quality. Finally, a high-throughput sequencer (MGI, Shenzhen) was used for sequencing, and the data splitting, denoising, and comparison were accomplished in the Linux system to obtain the gene mutation results of various samples.

The library construction and NGS were carried out in Cowin Bio., Jiangsu.

Statistical analysis

Data were analyzed using SPSS 26.0 software. The gene mutation frequencies are expressed as the numbers of cases or percentages. Inter-group differences were examined byχ2and Fisher exact tests.P<0.05 was considered statistically significant.

RESULTS

Clinical data analysis of CRC patients

Among the CRC patients corresponding to 50 tumor tissue samples, 24 were female (48.00%) and 26 were male (52.00%). The age span was large (35–84 years), and the median age was 59.5 years.According to the histopathological results of tumor tissues, there were 29 cases of rectal cancer (58.00%),11 cases of left colon cancer + sigmoid colon cancer (22.00%), and 10 cases of right colon cancer (20.00%).Twelve cases (24.00%) were well-differentiated + well-to-moderately differentiated, 25 (50.00%) were moderately differentiated, and 12 (24.00%) were moderately-to-poorly differentiated + poorly differentiated (the degree of differentiation was not determined in 1 of 50 colorectal cancer tumor tissues).Among 19 cases (38.00%) of lymph node metastasis, only 1 had distant metastasis (metastasis to the liver). Regarding pathological staging (according to eighth edition of AJCC TNM staging system), 33 cases (66.00%) were at stages I + II, and 17 (34.00%) were at stages III + IV (Table 2).

Gene mutations in various samples

As is clear from the NGS testing results of CRC specimens in Table 3, genes likeTP53,APC,KRAS,PTEN,MET,HRAS,andPDGFRAhad high incidences of mutations in tumor tissues and preoperative stools, while in postoperative stools, the mutation frequencies ofTP53,HRAS,andPDGFRAwere high.Despite high mutation frequencies ofHRASandPDGFRAin the three groups of samples, their mutations were mostly same sense mutations, meaning that there were only base mutations, without changes in amino acids, which had no clinical significance.

Table 1 The reaction system

Table 2 Clinical data analysis of colorectal cancer patients

As shown in Table 4, thePDGFRA,HRAS, andKITmutations occurred in only 3 of 20 stool samples in the normal control group. In the benign control group, only 4 stool samples hadPDGFRA,HRAS,KIT,orSTK11mutations.

Since every gene has numerous mutation sites, different site mutations have different meanings. We classified and summarized the results of gene sequencing, and queried them online (https://www.ncbi.nlm.nih.gov/clinvar/). After combining the start and stop positions of each gene mutation site, various gene mutation sites were classified into four types depending on their pathogenicity:"Pathogenic", "benign", "unknown clinical significance", and "undetected in the system" (hereinafter referred to as "undetected"). "Benign" gene mutation sites imply a large number of mutations that can occur multiple times in the normal population, albeit less than 1% incidences. Mutation sites of"unknown clinical significance" indicate that such mutations are neither pathogenic nor benign, on whom the experimental conclusions in the literature reports are inconsistent. There is insufficient evidence to clearly classify them, and further exploration is required. "Undetected" mutation sites are identified by comprehensively checking multiple items such as the exons, the types of amino acid mutations, and the start and end positions of mutations. In the case of non-conformity of any item, the corresponding mutation site is regarded as "undetected". In Table 3, the mutations of thePTENgene in tumor tissues were quite frequent, all of which were mutation sites of "unknown clinical significance".Among the 27 cases of preoperative stools, 1 had pathogenic mutation. Given the small case number, the possible cause of tumor heterogeneity was not ruled out. TheMETgene exhibited high mutation frequencies in both tumor tissues and preoperative stools, all of which were "benign" mutations upon query. The fourMETmutations in the postoperative stools were all "undetected" mutations. In Table 4,the mutation sites ofKITandSTK11in the benign control group were benign mutation sites, while the rest were all same sense mutations and had no clinical significance.

Table 3 Summary of gene mutations in colorectal cancer patients

This study focused on analyzing the pathogenic mutation sites of various genes.

Comparison between preoperative stool and tumor tissue results in CRC patients

As shown in Table 5, theTP53,APC, andKRASgenes had high incidences of pathogenic mutations in tumor tissues and preoperative stools among CRC patients. Despite high mutation frequency ofPIK3CAin tumor tissues, its pathogenic mutation frequency in preoperative stools was only 3.70% (1/27). The possible causes of tumor heterogeneity or individual disparity could not be ruled out. Although other genes had pathogenic mutations, their mutation frequencies were rather low.

Comparison between preoperative/postoperative stool results of CRC patients and control stool results

For CRC patients, the pathogenic mutation incidences ofAPC,TP53,andKRASin the preoperative stools were 11.11% (3/27), 37.04% (10/27), and 25.93% (7/27), respectively. There were mutations ofPIK3CA,BRAF,andMPLgenes in the postoperative stools. Although they were pathogenic, only onecase was found for each gene. The possible cause of tumor heterogeneity could not be ruled out. No pathogenic gene mutation sites were detected in stool samples from the normal and benign control groups (Tables 6–8).

Table 4 Fecal gene mutation results in the control groups

Table 5 Comparison of pathogenic mutation sites in preoperative stools vs tumor tissues

According to a combination of Tables 5–8, theTP53,APC,andKRASgenes all had high mutation frequencies in the preoperative stools and tumor tissues of CRC patients, while exhibiting no pathogenic mutations in the postoperative stools, or in the stool samples of the normal or benign control group.This indicates that tumor-specific DNA can be detected in the preoperative stools of CRC patients, and that the preoperative fecal expression of tumor-associated genes can reflect the gene mutations in tumor tissue to some extent. Compared to the postoperative stools and the stool samples of the two control groups, the pathogenic mutation frequencies ofTP53andKRASwere significantly higher in the preoperative stools (P< 0.05), suggesting that fecalTP53andKRASgenes can be used for the screening,diagnosis, and prognostic prediction of CRC. Contrastively, the pathogenic mutation frequency of theAPCgene in the preoperative stools differed insignificantly from that in the postoperative stool or the two control groups (P> 0.05), which was probably associated with the excessively small sample size.

Table 6 Comparison of pathogenic mutation sites in preoperative stools vs normal control stools

Mutation site analysis for TP53, KRAS, and APC genes

As shown in Table 9,TP53mutations occurred in 31 of 50 tumor tissues, and the total mutation frequency was 40 times, of which pathogenic mutation sites accounted for 60.00% (24/40), with exon 5 R175H, exon 7 R248Q/W, and exon 8 R273H/C being the most common. The proportion of pathogenic mutation sites in 29 tumor tissues withKRASmutations was 90.00% (27/30), most of which were located in exon 2, with G12D/V/S being predominant, followed by G13D/C. Additionally, there was 1 case of exon 3 Q61H mutation. As for theAPCgene, all its pathogenic mutation sites were located in exon 17.

The 12 of 27TP53mutations in preoperative stools were pathogenic mutations, which were mostly exon 7 R248Q and exon 8 R273H/C. The pathogenic mutation sites of theKRASgene were primarily exon 2 G12D/V. As for theAPCgene, all its pathogenic mutation sites were located in exon 17(Table 10).

As are clear from Tables 9 and 10, the aforementioned pathogenic mutation sites occurred in both preoperative stools and tumor tissues, suggesting the possible correlation between preoperative stool and tumor tissue in CRC patients.

Combined TP53–KRAS detection

Among the 27 preoperative stools of CRC patients, 16 had pathogenic mutations of theTP53gene, with a true positive rate of 59.26%, and 10 had pathogenic mutations of theKRASgene, with a true positive rate of 37.04%. Pathogenic mutations ofTP53orKRASoccurred in 18 preoperative stools, showing a true positive rate of 66.67% (18/27). According to Table 11, the sensitivity and negative predictive value of the combined TP53–KRAS detection were 66.67% and 68.97%, respectively, which were higher than that those ofTP53orKRASmutation detection alone. This suggests that the combined testing ofTP53andKRASgenes can improve the detection rate.

Distribution of mutation sites of "unknown clinical significance"

Mutation sites of "unknown clinical significance" occurred for such genes as TP53,PTEN,KRAS,BRAF,andAKT1in the tumor tissues of CRC patients. The five mutation sites of "unknown clinical significance" forTP53were respectively exon 5 A159V/V172G/G154D/P153S and exon 7 S241T. In the case ofPTEN, 9 of 11 mutations were of "unknown clinical significance", all of which were exon 6 V133I.For theKRASgene, the mutation sites of "unknown clinical significance" accounted for 10.00% (3/30),all of which were exon 4 A146T. In preoperative stools, the mutation sites of "unknown clinical significance" were scattered inTP53,KRAS,AKT1,andSTK11. In the case ofTP53, such mutation sites were exon 5 V172G and exon 7 I251T/S241T, which were identical to the mutation locations in tumor tissues. For theKRASgene, the mutation site of "unknown clinical significance" was exon 4 A146T. As forAKT1andSTK11, such mutation sites were located separately in exon 3 E17K and exon 4 E199K.

Distribution of "undetected" mutation sites

As shown in Table 12, genes likeTP53,APC,EGFR,FBXW7,andNOTCH1in tumor tissues all had"undetected" mutation sites, among which exon 4 A84G ofTP53occurred more frequently (4/40). For theAPCgene, "undetected" mutation sites were all found in exon 17. The pathogenic mutation sites ofAPCin Tables 9 and 10 were also located in exon 17, suggesting that the "undetected" mutation sites ofAPCin Table 12 may be unreported new sites. They are probably associated with the carcinogenesis and progression of CRC, which require further exploration. Table 13 shows that 6 of 27TP53mutation sites in the preoperative stools were "undetected" sites. They were distributed in exon 4, exon 5, exon 7, and exon 8, among which there were three exon 7 N247D (mutation start and stop positions were both 7577542). Additionally, genes likeCDKN2A,HRAS,EGFR,IDH1,andPDGFRAalso had "undetected"mutations sites.

Table 7 Comparison of pathogenic mutation sites in preoperative stools vs benign control group

Table 8 Comparison of pathogenic mutation sites in preoperative stools vs postoperative stools

DISCUSSION

Studies have shown that the genetic factors coupled with prolonged unfavorable external factors induce the disturbance of homeostasis, which leads to a series of changes including angiogenesis, cell hyperproliferation, apoptosis evasion, and enhanced invasion capacity, ultimately resulting in the occurrence of malignancies. It is generally believed that CRC often develops slowly from adenomas,during which multiple steps are involved, including the activation of proto-oncogenes and theinactivation of tumor suppressor genes. The mutation statuses of genes likeTP53,APC,andKRAShave been reported to be closely associated with the carcinogenesis and progression of CRC[24,25]. With the introduction of the concept of "precision medicine", molecular targeted therapy has become a hot research topic in recent years. This NGS-based study investigated the value of fecal genetic testing in the CRC diagnosis.

Table 9 Analysis of pathogenic mutation sites in tumor tissues

As a tumor suppressor gene,TP53regulates downstream genes when cells are under stress, which induces apoptosis and cell cycle arrestviaa series of processes, thereby preventing the growth of tumor cells[26,27]. Its mutation is a key step in facilitating the adenoma–adenocarcinoma transition. In malignancies, the p53 protein encoded byTP53is mutant. The mutant p53 protein can lose its tumor suppressor effect and acquire new oncogenic functions to promote the invasion and metastasis of tumor cells[28]. In this study,TP53had 60.00% (24/40) and 59.26% (16/27) of mutations in the tumor tissues and preoperative stools of CRC patients, respectively, among which pathogenic mutation sites accounted for 60.00% (24/40) and 40.74% (11/27) of total mutation sites. Such pathogenic mutation sites as exon 5 C176Y, exon 6 R196*, exon 7 R248Q, exon 7 G245C/S, and exon 8 R273H/C existed in both types of samples, showing consistency. Such mutation hotspots asR273H,R248Q,andR282Whave been reported as typical "hotspots" in the DNA binding domain. Over 95% of missenseTP53mutations occur in the DNA binding domain, which can affect the ability of proteins to bind to their target DNA sequences[29-31]. Loet al[32] screened a patient with CRC lymphatic metastasis for tumor infiltrating lymphocytes. Since some tumor patients had TP53 R175H mutations during genetic testing, and HLA-A*0201 could limit the recognition of mutant TP53 R175H, they identified the T cell receptors capable of recognizing the TP53 R175H mutation hotspots, thereby finding patients with different tumor types who shared the same mutation hotspots. In their opinion, the mutantTP53could serve as an ideal common antigen target, and these specific hotspot mutations had potential as valuable therapeutic targets for a variety of tumor cells. A study by Olszewskiet al[33] suggested that changes in theoncogenic activity of TP53 led to acquisition of new functions for the mutated tumor suppressor p53 protein, and that the overexpression of such mutation hotspots as codons 175, 248, or 273 might be associated with the chemotherapeutic resistance or invasive changes. Their analysis ofTP53mutation database in human somatic cells revealed that patients carrying TP53 R248Q allele had earlier tumor onsets, but lower frequencies of metastasis.

Table 10 Analysis of pathogenic mutation sites in preoperative feces

Table 11 Colorectal cancer diagnosis results by TP53 and KRAS mutations in preoperative stools, n (%)

The 21 kDa protein encoded byKRAS, a member of the ras proto-oncogene family, can participate in the G protein signaling pathway and regulate the proliferation and differentiation of cells.KRASmutations can continuously activate the RAS-RAF-MAPK pathway, causing uncontrolled cellular proliferation and differentiation, thereby inducing CRC[14]. Research has demonstrated that the patients'resistance to EGFR-targeted drugs was closely associated with the activation of theKRASgene[34]. In the present study, the incidences ofKRASmutations in tumor tissue and preoperative stools of CRC patients were 58.00% (29/50) and 37.04% (10/27), respectively. Among them, pathogenic mutation sites existed in 27 tissues and 7 preoperative stools, with the most common ones being G12D/V and G13D located in exon 2. Armengolet al[35] explored the gene mutations in stool samples of 52 CRC patients.They found 12 cases ofKRASgene mutations, of which 58.3% (7/12) were located in codons 12 and 13 of exon 2. Knight argued that when there was a KRAS G12D carcinogenic mutation in the body, the formation or growth of intestinal adenomas was not inhibited by rapamycin[36]. Compared to KRAS G12D, patients with KRAS A146T mutation could have stronger resistance to EGFR tolerance, and better survival rates[37]. Mullaet al[38] analyzed the RAS mutations in 51 CRC patients, finding that 43% of them had RAS mutations, 91% of which wereKRASmutations. The majority ofKRASmutations were located in G12D, followed by G13D. They also found that the incidence ofKRASmutations was slightly higher among individuals under 50 years of age than that of the elderly (> 50 years). Ottaianoet al[39]reported that mutantKRASwas a negative influencing factor of prognosis in CRC patients. Among these patients, those with G12C/S mutation showed the shortest survival and the worst prognosis.Varshaviet al[40] investigated howKRASmutations in different codons affected the metabolic pathways. According to their findings, metabolic differences existed between variousKRASmutation sites, which might be one of the reasons for the varying efficacies of anti-cancer therapies across patients. Hence, more effective, individualized molecular targeted drugs can be searched based on such metabolic differences.

Table 12 "Undetected" gene mutation sites in tumor tissues

Table 13 "Undetected" gene mutation sites in preoperative feces

Mutations ofAPCwere linked to familial adenomatous polyposis, and its inactivation and the activation of Wnt signaling pathway also played pivotal roles in the occurrence of CRC, which could affect the adhesion between tumor cells to accelerate the cellular migration[41,42]. It has been reported that in sporadic CRC, most of theAPCmutations occurred around the mutation cluster region,especially between codons 1286 and 1513[43]. In this study, the incidences ofAPCmutations in tumor tissues and preoperative stools were 24.00% (12/50) and 11.11% (3/27), respectively. As shown in Tables 5 and 6, the pathogenic mutation frequency of theAPCgene in preoperative stools differed insignificantly from that in control stools (P= 0.370), which thus cannot serve as a specific gene for CRC screening and diagnosis. This may be attributed to the small sample size of this study, which remains further exploration. We found that all the pathogenic mutation sites ofAPCwere located in exon 17,most of which used the mutations of certain amino acids as stop codons, thus allowing premature termination of mRNA translation (Tables 8 and 9). Luo's team performed NGS-based gene sequencing on the blood and tissue samples of 22 CRC patients and 21 patients with intestinal polyps, in order to understand gene mutations. Through screening, they found that theAPCmutations were most common in both types of samples, which were mostly pathogenic mutation sites. These pathogenic mutation sites of theAPCgene accounted for 77.3% (17/22) of the total pathogenic mutation sites in all samples[44].

In addition to the pathogenic mutation sites of various samples, this study also found multiple mutation sites of "unknown clinical significance", which were scattered in such genes asTP53,PTEN,KRAS,BRAF,AKT1, andPIK3CA. There were also some gene mutation sites that had not been found in the Clinvar database.TP53exon4 A84G andEGFRexon20 I821T were present in both the preoperative stools and tumor tissues of CRC patients. The mutation start and stop positions were both 7579436 for the former, while 55249164 for the latter.

CONCLUSION

Conclusively, NGS-based fecal DNA detection can be used as a complementary technique for the CRC diagnosis. FecalTP53andKRAScan serve as a specific gene combination for the screening, diagnosis,prognostic prediction, and recurrence monitoring of CRC. Moreover, the combined TP53–KRAS testing can improve the CRC detection rate, which provides a molecular reference for developing personalized treatment programs. The mutation sites of "unknown clinical significance" and the "undetected"mutation sites may be unreported new mutation sites associated with CRC. Whether these genes are linked to the CRC diagnosis and prognosis remains to be further explored and researched.

ARTICLE HIGHLIGHTS

Research conclusions

NGS-based fecal genetic testing can be used as a complementary technique for the CRC diagnosis. FecalTP53andKRAScan be used as specific genes for the screening, diagnosis, prognostic prediction, and recurrence monitoring of CRC. Moreover, the combined testing ofTP53andKRASgenes can improve the CRC detection rate.

Research perspectives

Fecal genetic detection is a new method for CRC diagnosis, which has the advantages of noninvasiveness, convenient sampling, and dynamic monitoring. Although the sensitivity of fecal genetic test in CRC screening is low, it is certain that it has great potential and broad prospects in the diagnosis and prognosis assessment of CRC. In addition, the "undetected" mutation sites in preoperative stools of CRC patients and the "unknown clinical significance" mutation sites are related to the occurrence and development of CRC, which requires further research and exploration.

FOOTNOTES

Author contributions:He SY performed most of the experiments, analyzed some data, and drafted the manuscript; Li YC and Wang Y provided specimens for the study; Peng HL and Zhou CL gave constructive guidance on the study;Zhang CM contributed to the statistics; Chen SL and Yin JF performed some of the experiments; Lin M designed and supervised the study and edited the manuscript; all authors approved the final version of the article.

Supported byTaizhou Social Development Plan, No. TS202004; Natural Science Foundation of Nanjing University of Chinese Medicine China, No. XZR2020093; and Taizhou People's Hospital Medical Innovation Team Foundation, No.CXTDA201901.

Institutional review board statement:This study was reviewed and approved by the Ethics Committee of Taizhou people's Hospital in Jiangsu Province (No. KY201912501).

Conflict-of-interest statement:All the authors report no relevant conflicts of interest for this article.

Data sharing statement:No additional data are available.

Open-Access:This article is an open-access article that was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution NonCommercial (CC BYNC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is noncommercial. See: https://creativecommons.org/Licenses/by-nc/4.0/

Country/Territory of origin:China

ORCID number:Si-Yu He 0000-0003-3411-8616; Ying-Chun Li 0000-0003-0650-531X; Yong Wang 0000-0002-6348-4562;Hai-Lin Peng 0000-0001-8332-1704; Cheng-Lin Zhou 0000-0003-2661-2067; Chuan-Meng Zhang 0000-0002-5642-1512;Sheng-Lan Chen 0000-0002-8366-451X; Jian-Feng Yin 0000-0002-5638-263X; Mei Lin 0000-0001-9815-2966.

S-Editor:Fan JR

L-Editor:Wang TQ

P-Editor:Fan JR