A commentary has been posted on this article:
Commentary: Genomic Analysis Reveals Heterogeneity Between Lesions in Synchronous Primary Right-Sided and Left-Sided Colon Cancer
Hanqing Hu1†
Qian Zhang2†Rui Huang1†Zhifeng Gao3
Ziming Yuan1Qingchao Tang1Feng Gao1Meng Wang2Weiyuan Zhang1Tianyi Ma1Tianyu Qiao1Yinghu Jin1
Guiyu Wang2*- 1Colorectal Cancer Surgery Department, The Second Affiliated Hospital of Harbin Medical University, Harbin, China
- 2Colorectal Cancer Surgery Department, Cancer Hospital of the University of Chinese Academy of Sciences (Zhejiang Cancer Hospital), Hangzhou, China
- 3Department of Surgery, The Second Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
Background: The synchronous primary right-sided and left-sided colon cancer (sRL-CC) is a peculiar subtype of colorectal cancer. However, the genomic landscape of sRL-CC remains elusive.
Methods: Twenty-eight paired tumor samples and their corresponding normal mucosa samples from 14 patients were collected from the Second Affiliated Hospital of Harbin Medical University from 2011 to 2018. The clinical–pathological data were obtained, and whole-exome sequencing was performed based on formalin-fixed and paraffin-embedded samples of these patients, and then, comprehensive bioinformatic analyses were conducted.
Results: Both the lesions of sRL-CC presented dissimilar histological grade and differentiation. Based on sequencing data, few overlapping SNV signatures, onco-driver gene mutations, and SMGs were identified. Moreover, the paired lesions harbored a different distribution of copy number variants (CNVs) and loss of heterozygosity. The clonal architecture analysis demonstrated the polyclonal origin of sRL-CC and inter-cancerous heterogeneity between two lesions.
Conclusion: Our work provides evidence that lesions of sRL-CC share few overlapping mutational signatures and CNVs, and may originate from different clones.
Introduction
Synchronous multiple primary colorectal cancer (SM-CRC) refers to at least two primary lesions identified in a single patient. Synchronous primary right-sided and left-sided colon cancer (sRL-CC) is a peculiar subtype of SM-CRC. Previous studies have shown that genomic molecular aberration is a contributor to synchronous neoplasia (Ogino et al., 2006). However, the genomic landscape of sRL-CC remains to be elucidated.
Recently, many works have shown that solitary left-sided colon cancer and right-sided colon cancer have different biological behaviors. Clinical studies revealed that patients with right colon cancer (RCC) have a worse overall survival (Ishihara et al., 2018). Patients with metastatic left colon cancer (LCC) respond more effectively to cetuximab-based target therapy compared to those with metastatic RCC (Elez et al., 2015; Arnold et al., 2017; Tejpar et al., 2017). In addition to clinical trials, molecular studies have shown different genomic phenotypes in LCC and RCC. TP53 and APC are mutated more frequently in LCC, whereas PIK3CA, CTNNB1, ATM, PTEN, and BRCA1 are mutated more frequently in RCC (Yamauchi et al., 2012; Takahashi et al., 2016; Salem et al., 2017). The consensus molecular subtype (CMS) according to the transcriptome has divided colorectal cancer into four subtypes: CMS1 (immune activation and JAK-STAT activation), CMS2 (WNT activation, MYC activation, EGFR or SRC activation, and VEGF or VEGFR activation), CMS3 (DNA damage repair, glutaminolysis, and lipidogenesis), and CMS4 (mesenchymal transition and complement activation immunosuppression) (Guinney et al., 2015). CMS1 and CMS3 are frequent in RCC, whereas CMS2 and CMS4 are frequent in LCC (Dienstmann et al., 2017). Moreover, LCC and RCC also present different immune landscapes. For RCC, decreased infiltration of CD8+ T cells and Th1 cells were identified. For LCC, infiltration of CD56high natural killer cells and activation of IFN-α signaling were identified (Zhang et al., 2018). Despite continuous anatomy, left-sided and right-sided colon cancer could be two sides of a coin in solitary colorectal cancer, which provokes us to explore the molecular phenotype of sRL-CC.
Previous studies had reported the heterogeneity and independent genetic origin of synchronous colorectal cancer (Cereda et al., 2016; Wang et al., 2018), but these studies have not compared the genetic phenotype according to the tumor location. In our present work, we analyzed 14 cases of sRL-CC to investigate single nucleotide variation, somatic mutation, and copy number alteration in sRL-CC patients based on whole-exome sequencing data.
Materials and Methods
Sample Collection
Nineteen sRL-CC patients were enrolled from March 2011 to October 2018 and five patients were excluded due to the small sample volume of the tumor tissue. The diagnosis was confirmed by two experts in the Department of Pathology from the Second Affiliated Hospital of Harbin Medical University. Fourteen patients denied hereditary history of CRC and were diagnosed with sporadic CRC. All the samples were formalin-fixed and paraffin-embedded. The splenic flexure, descending colon, and sigmoid colon were classified as the left-sided colon, while the caecum, hepatic flexure, and ascending colon were classified as the right-sided colon (Missiaglia et al., 2014). The clinicopathological data of the 14 sRL-CC were all available. All the samples were acquired with the approval of the ethics committee from the Second Affiliated Hospital of Harbin Medical University and written informed consent was obtained from all the participants.
DNA Extraction and Whole-Exome Sequencing
The DNA was extracted from FFPE samples using a QIAamp DNA FFPE Tissue Kit (Qiagen, China), according to the manufacturer’s protocol, after each slide was reviewed by two pathological experts to ensure tumor purity was greater than 50%. The quality of DNA and contamination was evaluated on 1% agarose gels and the concentration of DNA was calculated by a Qubit DNA Assay Kit in a Qubit 2.0 Fluorometer (Invitrogen, China). Two micrograms of genomic DNA was used to prepare the captured libraries by an Agilent SureSelect Human All Exon V5 kit (Agilent Technologies, China), following the manufacturer’s recommendations. The reads library was sequenced on the Illumina Novaseq 6000 platform.
Raw Data Processing for Calling Single Nucleotide Variants and Indels
The clean data were obtained after discarding the adapter and paired reads of the raw data from the Novaseq 6000 platform. Next, quality control was performed including reads number, error rate, and percentage of reads with average quality (>Q20 and >Q30). Burrows–Wheeler Aligner (BWA) software was used to map the paired-end clean reads to the reference genome (UCSC Human Genome Reference hg38) (Li and Durbin, 2009). Aligned reads were processed in terms of marked duplicates, realignment of indels, and base recalibration by the Genome Analysis Toolkit (GATK) (McKenna et al., 2010). Variants were identified in accordance to dbSNP (Sherry et al., 2001) and the 1000 Genomes database (Genomes Project et al., 2012), which was annotated by ANNOVAR (Wang et al., 2010). Next, SNVs and indels were identified by MuTech (Cibulskis et al., 2013) and Strelka (Saunders et al., 2012).
Based on these data, cluster analysis was conducted on 96 somatic mutational nucleotide types through nonnegative matrix factorization (Nik-Zainal et al., 2012), and three signatures were identified (Alexandrov et al., 2013a). These signatures were clustered based on 30 known signatures according to the COSMIC database to reveal the biological process of the signatures (Alexandrov et al., 2013b). The degree of similarity was assessed by the cosine similarity coefficient index.
Mutation Signature Analysis
To investigate the predisposing genes, we analyzed germline mutation in the normal tissue. Compared with the Cancer Gene Census (CGC) database (Gerlinger et al., 2014) using an in-house algorithm, the predisposing genes were identified. To determine driver mutation genes in carcinogenesis, the mutation spectrum was aligned with published driver mutations via OncodriveCLUSTL (Arnedo-Pac et al., 2019) and OncodriveFM (Gonzalez-Perez and Lopez-Bigas, 2012) software referring to four databases, the CGC513 database (Futreal et al., 2004), 125 mutation genes reported by Bert Vogelstein (Vogelstein et al., 2013), SMG127 (Kandoth et al., 2013), and Comprehensive435 (Tamborero et al., 2013). A significantly mutated genes test was applied to define the SMGs in the tumor, and a mutation relation rest was applied to explore the relationship between SMGs (Dees et al., 2012).
Phylogenetic Tress and Clonal Architecture Analysis
Phylogenetic trees analysis was performed based on the WES data. The branch and trunk lengths reflected the number of nonsynonymous mutations, which were also marked beside the branch and trunk (Gerlinger et al., 2014). In accordance with the Pyclone algorithm, clonal architecture analysis was performed (Roth et al., 2014).
Copy Number Analysis
The copy number variants (CNVs) were evaluated by control-FREEC software based on the WES data (Boeva et al., 2012). The profiles of CNVs were calculated with alignment to the BAM data. Next, the CNVs were normalized to obtain the number of CNVs in different regions on the chromosome.
Statistical Analyses
Statistical analyses were conducted by SPSS software (version 25.0). The comparison between the left-sided and right-sided lesions with SNP signature and TMB was performed by a two-tailed paired t test. The differences were considered to be significant at p < 0.05. The p values are described in the corresponding figure legend.
Results
Clinical Characteristics of Synchronous Primary Right-Sided and Left-Sided Colon Cancer
We collected 9,876 colorectal cancer patients, who were confirmed by pathological examination, in the Second Affiliated Hospital of Harbin Medical University from 2011 to 2018. A total of 7,369 patients received radical colectomy, and the samples were preserved by the formalin-fixed paraffin-embedded method. Among these patients, 67 were diagnosed with SM-CRC and 19 SM-CRC patients were identified as sRL-CC. The tumor location was evaluated by colonoscopy or computed tomography. During the surgery, the tumor location was confirmed by two experienced surgeons and recorded in the pathological reports. At last, 14 patients were enrolled for further analysis and 5 patients were excluded due to the small size of the tumor. The demographic characteristics of the 14 patients are presented in Table 1. As shown in Table 1, the median age of these patients was 69 years. Of the 14 patients, 4 were women and 10 were men. The CEA and CA-199 level was elevated (cutoff value, CEA: 5 ng/ml, CA-199: 30 U/ml) in most cases, suggesting a heavy tumor burden. As for the two lesions in most patients, the left lesions were located at the sigmoid colon and the right lesions were located at the right colon. However, in two patients, the left lesions were located at the left colon and right lesions were located at the right colon. Mucinous adenocarcinoma was identified in right-sided lesions but not in left-sided lesions. The right-sided lesions showed an advanced disease compared to the left-sided lesions according to the pathological stages and indicators, which suggested that left-sided and right-sided lesions may be at different stages of CRC development.
TABLE 1
TABLE 1. Demographic characteristic of the patients in our study.
Single Nucleotide Variants Analysis of Synchronous Primary Right-Sided and Left-Sided Colon Cancer
We analyzed the single nucleotide variants (SNVs) based on the WES data of 14 paired tumor specimen and their corresponding normal mucosa. Six types of SNVs were identified in sRL-CC patients, and C>T/G>A transition was preponderant nucleotide substitution in all the samples (Supplementary Figure S1A). Then, we evaluated the percentage of each type, which were similar in left and right lesions (Supplemental Figure S1B). To explore the signatures of the SNV spectrum, we clustered SNVs to obtain three signatures (Figure 1A). In signature A, C>T/G>A and T>C/A>G transition was abundant, signature B featured C>A/G>T and C>T/G>A transition, and signature C was defined as a high level of C>T/G>A transition (Figure 1A). To investigate the function of the three signatures, we performed cosine similarity analysis according to the COSMIC database (Blokzijl et al., 2018). It revealed that signature A was mostly related to signature 6 (correlation coefficient: 0.837), which was associated with defective DNA mismatch repair. Signature B was mostly related to signature 1 (correlation coefficient: 0.725), which was associated with the biological process initiated by spontaneous deamination of 5-methylcytosine. Signature C was also mostly related to signature 6 (correlation coefficient: 0.920) (Figure 1B). To evaluate the composition of signatures in paired samples, we analyzed the left-sided and right-sided lesions as a whole. In spite of no statistical significance resulting from the small number of patients, the average percentage of signature A and C contribution was higher in right lesions, but the average percentage of signature B contribution was higher in left lesions (Figure 1C). Then, we evaluated the signatures in a single patient, respectively. For instance, signature B was dominated in L2, but signatures A and C were dominated in R2. In patient 9, signature A was dominated in the left lesion and signature C was dominated in the right lesions. There were still some special cases, such as patients 5–8 and 13, in whom the dominated signature was similar (Figure 1D). This finding suggests that distinctive mutational processes happen in the right or left lesion during carcinogenesis.
FIGURE 1
FIGURE 1. Single nucleotide variants analysis of sRL-CC. (A) Three SNV signatures are obtained from WES data in 14 sRL-CC patients. (B) Correlation analysis of three SNV signatures with known 30 signatures. The bar represents the coefficient index. (C) The box plot shows the percentage of the three signatures in left-sided and right-sided lesions when analyzed as a whole. Signature A: p = 0.722, two-tailed paired t test; signature B: p = 0.539, two-tailed paired t test; and signature C: p = 0.589, two-tailed paired t test. (D) The distribution of the three signatures in two lesions of each patient is presented by the heatmap.
Mutation Signature Between Lesions of Synchronous Primary Right-Sided and Left-Sided Colon Cancer
We analyzed the germline mutation in the normal tissues of these patients to define predisposing genes using the in-house method and several frequent predisposing genes were identified. The nonsense mutation of PDE4DIP was found in all of 14 patients. The in-frame insertion of BPTF was found in 13 patients except patient 6. The in-frame deletion of MAP3K1/4 and ZNF384 was found in the majority of the patients. The missense mutation of MUC20 was found in 11 out of 14 patients (Figure 2A). Cancer driver gene mutation plays a significant role in carcinogenesis, and we explored the oncodriver mutation in sRL-CC. The completely matched lesions failed to be identified in all 14 patients. Although similar in a small part, the features of cancer driver gene mutation between lesions were diversified in most patients. For example, frameshift insertion of APC and frameshift deletion of TP53 was detected in L2 but not in R2, whereas missense mutation of SOX9, frameshift deletion of ZFHX3, and missense mutation of SYNE1 were detected in R2 but not in L2 (Figure 2B). The significant mutation genes (SMGs) reflect the mutational phenotype of diseases; and then, we wanted to investigate the SMGs in all the lesions. When treated as a whole, the tumor mutation burden (TMB) was alike between left- and right-sided lesions (Figure 2C). When analyzed separately, TMB was different in patients 2, 3, 4, 8, 9, 10, and 14 and was similar in patients 5, 6, 7, 11, 12, and 13. Particularly, the TMB was almost the same in patient 1 (L: 4.1537674 vs. R: 4.1569451) (Figure 2D). Among these significant mutation genes, the APC mutation tended to happen in the left lesions, while BAX, KRAS, and SOX9 mutations tended to happen in right-sided lesions. The PI3KCA mutation was also distributed unevenly, which was found in L2, L6, L9, R1, and R7 but not in the corresponding lesions on the other side (Figure 2E). The frequent mutation genes identified in sRL-CCs were similar to those identified in solitary colon cancer. Two studies identified eight frequent mutation genes (APC, TP53, SMAD4, PIK3CA, KRAS, ARID1A, SOX9, and FAM123B). APC, TP53, KRAS, and SOX9 were the top four mutation genes in sRL-CCs (Cancer Genome Atlas, 2012; Vasaikar et al., 2019). The ACVR2A mutation was found in hypermutation colon cancers (Vasaikar et al., 2019), which was also identified in sRL-CCs. Mutation genes in a sample may have a synergistic effect or a mutually exclusive effect. Through mutation relation rest analysis (MRT), a synergistic and mutually exclusive relationship was identified. We found that ACVR2A had a synergistic effect with PLEC, SLC4A11, and CEL (Supplementary Figure S2A). Besides, we found that the APC mutation was mutually exclusive with EPS8L1, PPP1R12C, SETD1B, RNF43, and SLC4A11. The KRAS mutation was mutually exclusive with TP53, SLC1A7, and PIK3CA (Supplementary Figure S2B). These works showed the different mutational landscape in two lesions of a single sRL-CC patient.
FIGURE 2
FIGURE 2. The mutational analysis of sRL-CC. (A) The frequent predisposing genes are presented by the heatmap. The tumor mutation burden is marked on the top of the heatmap, the name of predisposing genes is marked on the left, and the number of mutated samples is marked on the right. (B) The oncodriver gene mutations are presented by the heatmap. The tumor mutation burden is marked on the top of the heatmap, the name of predisposing genes is marked on the left, and the number of mutated samples is marked on the right. (C) The tumor mutation burden is calculated in the left-sided and right-sided lesions. (D) The tumor mutation burden is calculated in both lesions of each patient. (E) The significant mutated genes are presented by the heatmap. The tumor mutation burden is marked on the top of the heatmap.
Copy Number Variant Analysis of Synchronous Primary Right-Sided and Left-Sided Colon Cancer
Subsequently, we assessed the CNVs in both lesions of sRL-CCs. The GISTIC curve showed that amplification of the chromosome was predominant and more deletion events were identified in left-sided lesions (Figure 3A). The distribution of somatic CNVs in left-sided lesions was also different from those in right-sided lesions (Figure 3B). Loss of heterozygosity (LOH) was evaluated by β-allelic frequency (BAF). The right-sided and left-sided lesions shared different BAF (Figure 3B). When analyzed separately, the heterogeneity stood out more clearly. As shown in the Circos plot (Figure 3C), the copy number was almost normal in L2 and R4 lesions, whereas the copy number was amplified in R2 and L4 lesions. Similar phenomena were also found in other patients of our cohort (Supplementary Figure S3). These data imply that the lesions share different patterns with respect to copy number and different stages during carcinogenesis.
FIGURE 3
FIGURE 3. The copy number variants analysis of sRL-CC. (A) GISTIC analysis of left-sided and right-sided lesions. The G-score represents the degree of CNV amplification (upper panel) and deletion (lower panel). L-group means left-sided lesions and R-group means right-sided lesions. (B) The upper chart illustrates the distribution of copy number variants. Red presents the CNV gain, green presents normal CNV, and blue presents CNV loss. The lower chart illustrates the distribution of β allelic frequency (BAF). Orange means the normal distribution of allele and blue means loss of heterozygosity. (C) Circos plots illustrate the molecular landscape of genome in patients 2 and 4 in terms of copy number variants in each patient. The first circle means the depth of sequencing, the second circle means the density of SNP insert and deletion, and the third circle means the distribution of CNVs. In the third circle, red means the CNV gains, green means normal copy number, and blue means CNV loss.
Clonal Architecture and Evolution Analysis of Synchronous Primary Right-Sided and Left-Sided Colon Cancer
To acquire insights into the origin of both lesions shaping sRL-CC tumorigenesis, phylogenetic trees analyses were constructed to evaluate the ancestral relationship of individual lesions, significantly overlapping variant sets failed to be identified in most cases in our cohort. None of the overlapping variant sets were identified in patients 3, 4, 5, 7, 8, 9, 10, 11, 13, and 14, and few overlapping variant sets were identified in patients 1, 2, and 6. But as an exception, both the lesions in patient 12 shared 53 overlapping variant sets (Figure 4A and Supplementary Figure S4A). Clonal architecture analysis revealed different clusters between two lesions. Most of the patients shared few overlapping clusters. In patients 3 and 12, relatively more sharing clusters were identified with some independent clusters. However, in patient 9, a single shared cluster was identified (Figure 4B and Supplementary Figure S4B). These data may suggest different origins and clonal architectures in the lesion of sRL-CCs.
FIGURE 4
FIGURE 4. The clonal architecture and evolution analysis of sRL-CC. (A) The phylogenetic tress analysis of patients 1–6. The trunk represents the common mutations and the branch represents unique mutations in each lesion. The number of mutations is marked. (B) The scatter plot demonstrates cancer cell fraction (CCF) of the mutations in both lesions of patients 1–6. Different colors represent different clusters. The clusters located in the middle area of the plot mean shared clusters by both lesions, and those located near the X-axis and Y-axis mean unique clusters in each lesion.
Discussion
In this study, we performed an unprecedented molecular characterization of both lesions and corresponding normal adjacent tissue in sRL-CCs with comprehensive data from WES. Our work revealed that heterogeneity occurred in both lesions in the same patient and each lesion was in a different stage during carcinogenesis. Our study is in line with the multi-omics analysis in single primary colorectal cancer. Vasaikar et al. reported some significant mutations such as APC, TP53, KRAS, SOX9, and ACVR2A, which were also found in our cohort (Vasaikar et al., 2019). The mutation of mismatch repair genes is a key molecular event in colorectal cancer (Sahin et al., 2019; Schrock et al., 2019). Our work has also identified that SNV signature is associated with defective DNA mismatch repair. All this evidence confirms that our data are accurate enough to reflect the biological process in sRL-CCs.
Recently, more and more attention has been paid to the multiple primary tumors. Grolleman et al. has reported the mutation of the base excision repair gene NTHL1 that can trigger the development of malignancy in many organs (Grolleman et al., 2019). Our work also showed the mutational signature involved in the dysfunction of DNA mismatch repair in sRL-CCs. Our study has also highlighted predisposing genes in sRL-CC. The nonsense mutation of PDE4DIP was identified in all the cases in our cohort. The PDE4DIP gene encodes protein myomegalin and exerts its function as an anchor to sequester components to Golgi and/or centrosomes (Wang et al., 2014; Wu et al., 2016; Bouguenina et al., 2017; Yang et al., 2017). However, the role of PDE4DIP in cancers is poorly understood. Thus, prospective studies are needed to explore the role of nonsense mutation in PDE4DIP of sRL-CCs.
Above all, our study reveals heterogeneity of both lesions in sRL-CCs. We have identified different copy numbers in paired right-sided and left-sided lesions, which indicates that the lesions are in different stages of carcinogenesis. LOH is an important mechanism for the disability of tumor suppressor genes. Through the control-FREEC method (Boeva et al., 2012), the regions of LOH have been identified. The paired lesions share unmatched distribution of LOH on the chromosomes. This evidence suggests that inconsistent degrees of genome instability happen in the both lesions. Combined with the phylogenetic trees and clonal architecture analysis, we have proved the multi-origins of the paired lesions and different clonal fractions in the paired lesions. Studies based on multiple primary cancer in other organs have drawn a similar conclusion. Ma et al. reported that multicentric lesions harbor distinct oncogenic alterations and genomic heterogeneity (Ma et al., 2017). In synchronous bilateral renal cancer, the lesions in the kidneys also originate from separated clones (Linehan, 2009). In line with these previous studies, our work has shown the polyclonal origin of sRL-CC.
Limitations do exist in our study. The sample size is small and only 14 patients are included. Besides, our study is a single center, retrospective study and all the samples have to be preserved by FFPE due to the long time span. Thus, a multi-center, prospective study could collect a larger number of sRL-CC patients and fresh tumor tissue in a shorter time. Moreover, genomic analysis has trouble in reflecting the whole scenarios of sRL-CC. Systematic studies including the transcriptome, epi-transcriptome, proteome, and metabolome can provide a wider horizon for sRL-CC.
In conclusion, we performed whole-exome sequencing analyses, which are suggestive of heterogeneity between lesions and the polyclonal origin of sRL-CC. Moreover, we illustrated the genomic landscape of sRL-CC and provided an insight into the molecular pattern of sRL-CC, which could make treatment more precise and effective.
Data Availability Statement
The datasets presented in this article are not readily available because the Chinese laws forbid the publication the genomic data based on the Chinese population. Requests to access the datasets should be directed to the corresponding author GW, Z3VpeXdhbmdAMTYzLmNvbQ==.
Ethics Statement
The studies involving human participants were reviewed and approved by the ethics committee from the Second Affiliated Hospital of Harbin Medical University. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
HH, RH, and GW: study design, data interpretation, preparation, and submission of the manuscript. ZY, QZ, QT, and FG: clinical and pathological data collection, DNA extraction, and DNA quality control. MW, WZ, TM, TQ, and YJ: whole-exsome sequencing and raw data collection.
Funding
This work was supported by the China Postdoctoral Science Foundation (Grant number: 2018M641856), Heilongjiang Postdoctoral Fund (Grant number: LBH-Z18200), Heilongjiang Applied and Development Program (Grant number: GA19C003), Natural Science Foundation of Heilongjiang Province (Grant number: H2017016), Postgraduate Research & Practice Innovation Program of Harbin Medical University (Grant number: YJSKYCX 2019-51HYD), and National Natural Science Foundation of China (Grant number: 81872034).
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 conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.689466/full#supplementary-material
References
Alexandrov, L. B., Nik-Zainal, S., Wedge, D. C., Aparicio, S. A., Behjati, S., Biankin, A. V., et al. (2013a). Signatures of Mutational Processes in Human Cancer. Nature 500, 415–421. doi:10.1038/nature12477
Alexandrov, L. B., Nik-Zainal, S., Wedge, D. C., Campbell, P. J., and Stratton, M. R. (2013b). Deciphering Signatures of Mutational Processes Operative in Human Cancer. Cell Rep. 3, 246–259. doi:10.1016/j.celrep.2012.12.008
Arnedo-Pac, C., Mularoni, L., Muiños, F., Gonzalez-Perez, A., and Lopez-Bigas, N. (2019). OncodriveCLUSTL: A Sequence-Based Clustering Method to Identify Cancer Drivers. Bioinformatics 35, 4788–4790. doi:10.1093/bioinformatics/btz501
Arnold, D., Lueza, B., Douillard, J.-Y., Peeters, M., Lenz, H.-J., Venook, A., et al. (2017). Prognostic and Predictive Value of Primary Tumour Side in Patients With RAS Wild-Type Metastatic Colorectal Cancer Treated with Chemotherapy and EGFR Directed Antibodies in Six Randomized Trials. Ann. Oncol. 28, 1713–1729. doi:10.1093/annonc/mdx175
Blokzijl, F., Janssen, R., van Boxtel, R., and Cuppen, E. (2018). Mutational Patterns: Comprehensive Genome-Wide Analysis of Mutational Processes. Genome Med. 10, 33. doi:10.1186/s13073-018-0539-0
Boeva, V., Popova, T., Bleakley, K., Chiche, P., Cappo, J., Schleiermacher, G., et al. (2012). Control-FREEC: A Tool for Assessing Copy Number and Allelic Content Using Next-Generation Sequencing Data. Bioinformatics 28, 423–425. doi:10.1093/bioinformatics/btr670
Bouguenina, H., Salaun, D., Mangon, A., Muller, L., Baudelet, E., Camoin, L., et al. (2017). EB1-Binding-Myomegalin Protein Complex Promotes Centrosomal Microtubules Functions. Proc. Natl. Acad. Sci. U.S.A. 114, E10687–E10696. doi:10.1073/pnas.1705682114
Cancer Genome Atlas (2012). Comprehensive Molecular Characterization of Human colon and Rectal Cancer. Nature 487, 330–337. doi:10.1038/nature11252
Cereda, M., Gambardella, G., Benedetti, L., Iannelli, F., Patel, D., Basso, G., et al. (2016). Patients with Genetically Heterogeneous Synchronous Colorectal Cancer Carry Rare Damaging Germline Mutations in Immune-Related Genes. Nat. Commun. 7, 12072. doi:10.1038/ncomms12072
Cibulskis, K., Lawrence, M. S., Carter, S. L., Sivachenko, A., Jaffe, D., Sougnez, C., et al. (2013). Sensitive Detection of Somatic Point Mutations in Impure and Heterogeneous Cancer Samples. Nat. Biotechnol. 31, 213–219. doi:10.1038/nbt.2514
Dees, N. D., Zhang, Q., Kandoth, C., Wendl, M. C., Schierding, W., Koboldt, D. C., et al. (2012). MuSiC: Identifying Mutational Significance in Cancer Genomes. Genome Res. 22, 1589–1598. doi:10.1101/gr.134635.111
Dienstmann, R., Vermeulen, L., Guinney, J., Kopetz, S., Tejpar, S., and Tabernero, J. (2017). Consensus Molecular Subtypes and the Evolution of Precision Medicine in Colorectal Cancer. Nat. Rev. Cancer 17, 79–92. doi:10.1038/nrc.2016.126
Elez, E., Argiles, G., and Tabernero, J. (2015). First-Line Treatment of Metastatic Colorectal Cancer: Interpreting FIRE-3, PEAK, and CALGB/SWOG 80405. Curr. Treat Options Oncol. 16, 52. doi:10.1007/s11864-015-0369-x
Futreal, P. A., Coin, L., Marshall, M., Down, T., Hubbard, T., Wooster, R., et al. (2004). A Census of Human Cancer Genes. Nat. Rev. Cancer. 4, 177–183. doi:10.1038/nrc1299
Genomes Project, C., Abecasis, G. R., Auton, A., Brooks, L. D., DePristo, M. A., Durbin, R. M., et al. (2012). An Integrated Map of Genetic Variation From 1,092 Human Genomes. Nature 491, 56–65. doi:10.1038/nature11632
Gerlinger, M., Horswell, S., Larkin, J., Rowan, A. J., Salm, M. P., Varela, I., et al. (2014). Genomic Architecture and Evolution of Clear Cell Renal Cell Carcinomas Defined by Multiregion Sequencing. Nat. Genet. 46, 225–233. doi:10.1038/ng.2891
Gonzalez-Perez, A., and Lopez-Bigas, N. (2012). Functional Impact Bias Reveals Cancer Drivers. Nucleic Acids Res. 40, e169. doi:10.1093/nar/gks743
Grolleman, J. E., de Voer, R. M., Elsayed, F. A., Nielsen, M., Weren, R. D. A., Palles, C., et al. (2019). Mutational Signature Analysis Reveals NTHL1 Deficiency to Cause a Multi-Tumor Phenotype. Cancer Cell 35, 256–266. doi:10.1016/j.ccell.2018.12.011
Guinney, J., Dienstmann, R., Wang, X., de Reyniès, A., Schlicker, A., Soneson, C., et al. (2015). The Consensus Molecular Subtypes of Colorectal Cancer. Nat. Med. 21, 1350–1356. doi:10.1038/nm.3967
Ishihara, S., Murono, K., Sasaki, K., Yasuda, K., Otani, K., Nishikawa, T., et al. (2018). Impact of Primary Tumor Location on Postoperative Recurrence and Subsequent Prognosis in Nonmetastatic Colon Cancers. Ann. Surg. 267, 917–921. doi:10.1097/sla.0000000000002206
Kandoth, C., McLellan, M. D., Vandin, F., Ye, K., Niu, B., Lu, C., et al. (2013). Mutational Landscape and Significance Across 12 Major Cancer Types. Nature 502, 333–339. doi:10.1038/nature12634
Li, H., and Durbin, R. (2009). Fast and Accurate Short Read Alignment With Burrows-Wheeler Transform. Bioinformatics 25, 1754–1760. doi:10.1093/bioinformatics/btp324
Linehan, W. M. (2009). Genetic Basis of Bilateral Renal Cancer: Implications for Evaluation and Management. J. Clin. Oncol. 27, 3731–3733. doi:10.1200/jco.2009.23.0045
Ma, P., Fu, Y., Cai, M. C., Yan, Y., Jing, Y., Zhang, S., et al. (2017). Simultaneous Evolutionary Expansion and Constraint of Genomic Heterogeneity in Multifocal Lung Cancer. Nat. Commun. 8, 823. doi:10.1038/s41467-017-00963-0
McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The Genome Analysis Toolkit: A MapReduce Framework for Analyzing Next-Generation DNA Sequencing Data. Genome Res. 20, 1297–1303. doi:10.1101/gr.107524.110
Missiaglia, E., Jacobs, B., D'Ario, G., Di Narzo, A. F., Soneson, C., Budinska, E., et al. (2014). Distal and Proximal colon Cancers Differ in Terms of Molecular, Pathological, and Clinical Features. Ann. Oncol. 25, 1995–2001. doi:10.1093/annonc/mdu275
Nik-Zainal, S., Alexandrov, L. B., Wedge, D. C., Van Loo, P., Greenman, C. D., Raine, K., et al. (2012). Mutational Processes Molding the Genomes of 21 Breast Cancers. Cell 149, 979–993. doi:10.1016/j.cell.2012.04.024
Ogino, S., Brahmandam, M., Kawasaki, T., Kirkner, G. J., Loda, M., and Fuchs, C. S. (2006). Epigenetic Profiling of Synchronous Colorectal Neoplasias by Quantitative DNA Methylation Analysis. Mod. Pathol. 19, 1083–1090. doi:10.1038/modpathol.3800618
Roth, A., Khattra, J., Yap, D., Wan, A., Laks, E., Biele, J., et al. (2014). PyClone: Statistical Inference of Clonal Population Structure in Cancer. Nat. Methods 11, 396–398. doi:10.1038/nmeth.2883
Sahin, I. H., Akce, M., Alese, O., Shaib, W., Lesinski, G. B., El-Rayes, B., et al. (2019). Immune Checkpoint Inhibitors for the Treatment of MSI-H/MMR-D Colorectal Cancer and a Perspective on Resistance Mechanisms. Br. J. Cancer 121, 809–818. doi:10.1038/s41416-019-0599-y
Salem, M. E., Weinberg, B. A., Xiu, J., El-Deiry, W. S., Hwang, J. J., Gatalica, Z., et al. (2017). Comparative Molecular Analyses of Left-Sided Colon, Right-Sided colon, and Rectal Cancers. Oncotarget 8, 86356–86368. doi:10.18632/oncotarget.21169
Saunders, C. T., Wong, W. S. W., Swamy, S., Becq, J., Murray, L. J., and Cheetham, R. K. (2012). Strelka: Accurate Somatic Small-Variant Calling From Sequenced Tumor-normal Sample Pairs. Bioinformatics 28, 1811–1817. doi:10.1093/bioinformatics/bts271
Schrock, A. B., Ouyang, C., Sandhu, J., Sokol, E., Jin, D., Ross, J. S., et al. (2019). Tumor Mutational Burden is Predictive of Response to Immune Checkpoint Inhibitors in MSI-High Metastatic Colorectal Cancer. Ann. Oncol. 30, 1096–1103. doi:10.1093/annonc/mdz134
Sherry, S. T., Ward, M. H., Kholodov, M., Baker, J., Phan, L., Smigielski, E. M., et al. (2001). dbSNP: The NCBI Database of Genetic Variation. Nucleic. Acids Res. 29, 308–311. doi:10.1093/nar/29.1.308
Takahashi, Y., Sugai, T., Habano, W., Ishida, K., Eizuka, M., Otsuka, K., et al. (2016). Molecular Differences in the Microsatellite Stable Phenotype Between Left-Sided and Right-Sided Colorectal Cancer. Int. J. Cancer. 139, 2493–2501. doi:10.1002/ijc.30377
Tamborero, D., Gonzalez-Perez, A., Perez-Llamas, C., Deu-Pons, J., Kandoth, C., Reimand, J., et al. (2013). Comprehensive Identification of Mutational Cancer Driver Genes Across 12 Tumor Types. Sci. Rep. 3, 2650. doi:10.1038/srep02952
Tejpar, S., Stintzing, S., Ciardiello, F., Tabernero, J., Van Cutsem, E., Beier, F., et al. (2017). Prognostic and Predictive Relevance of Primary Tumor Location in Patients With RAS Wild-Type Metastatic Colorectal Cancer. JAMA Oncol. 3, 194–201. doi:10.1001/jamaoncol.2016.3797
Vasaikar, S., Huang, C., Wang, X., Petyuk, V. A., Savage, S. R., Wen, B., et al. (2019). Proteogenomic Analysis of Human Colon Cancer Reveals New Therapeutic Opportunities. Cell 177, 1035–1049. doi:10.1016/j.cell.2019.03.030
Vogelstein, B., Papadopoulos, N., Velculescu, V. E., Zhou, S., Diaz, L. A., and Kinzler, K. W. (2013). Cancer Genome Landscapes. Science 339, 1546–1558. doi:10.1126/science.1235122
Wang, K., Li, M., and Hakonarson, H. (2010). ANNOVAR: Functional Annotation of Genetic Variants From High-Throughput Sequencing Data. Nucleic. Acids Res. 38, e164. doi:10.1093/nar/gkq603
Wang, X., Fang, H., Cheng, Y., Li, L., Sun, X., Fu, T., et al. (2018). The Molecular Landscape of Synchronous Colorectal Cancer Reveals Genetic Heterogeneity. Carcinogenesis 39, 708–718. doi:10.1093/carcin/bgy040
Wang, Z., Zhang, C., and Qi, R. Z. (2014). A Newly Identified Myomegalin Isoform Functions in Golgi Microtubule Organization and ER-Golgi Transport. J. Cell Sci. 127, 4904–4917. doi:10.1242/jcs.155408
Wu, J., de Heus, C., Liu, Q., Bouchet, B. P., Noordstra, I., Jiang, K., et al. (2016). Molecular Pathway of Microtubule Organization at the Golgi Apparatus. Dev. Cell. 39, 44–60. doi:10.1016/j.devcel.2016.08.009
Yamauchi, M., Morikawa, T., Kuchiba, A., Imamura, Y., Qian, Z. R., Nishihara, R., et al. (2012). Assessment of Colorectal Cancer Molecular Features Along Bowel Subsites Challenges the Conception of Distinct Dichotomy of Proximal Versus Distal Colorectum. Gut 61, 847–854. doi:10.1136/gutjnl-2011-300865
Yang, C., Wu, J., de Heus, C., Grigoriev, I., Liv, N., Yao, Y., et al. (2017). EB1 and EB3 Regulate Microtubule Minus End Organization and Golgi Morphology. J. Cell Biol. 216, 3179–3198. doi:10.1083/jcb.201701024
Keywords: synchronous multiple primary cancer, genome, right-sided colon cancer, left-sided colon cancer, heterogeneity
Citation: Hu H, Zhang Q, Huang R, Gao Z, Yuan Z, Tang Q, Gao F, Wang M, Zhang W, Ma T, Qiao T, Jin Y and Wang G (2021) Genomic Analysis Reveals Heterogeneity Between Lesions in Synchronous Primary Right-Sided and Left-Sided Colon Cancer. Front. Mol. Biosci. 8:689466. doi: 10.3389/fmolb.2021.689466
Received: 31 March 2021; Accepted: 07 June 2021;
Published: 04 August 2021.
Edited by:
Xin Wang, City University of Hong Kong, SAR ChinaReviewed by:
Feng Gao, The Sixth Affiliated Hospital of Sun Yat-sen University, ChinaWei Wang, Huzhou Maternity and Child Health Care Hospital, China
Copyright © 2021 Hu, Zhang, Huang, Gao, Yuan, Tang, Gao, Wang, Zhang, Ma, Qiao, Jin and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guiyu Wang, Z3VpeXdhbmdAMTYzLmNvbQ==
†These authors have contributed equally to this work and shared first authorship