A Novel Serum Exosomal miRNA Signature in the Early Prediction of Persistent Organ Failure in Patients with Acute Pancreatitis.

OBJECTIVE: Current study aims to investigate whether serum exosomal microRNAs (miRNAs) could be potential biomarkers in predicting APs with POF at early phase. BACKGROUND: Novel biomarkers are sorely needed for early prediction of persistent organ failure (POF) in acute pancreatitis (AP) patients. METHODS: In the discovery stage, exosomal miRNAs were profiled in sera from APs with or without POF (5 vs. 5) using microarrays. POF-associated miRNA signatures then were assessed in training cohort (n=227) and further validated in three independent cohorts (n=516), including one nested case-control cohort. RESULTS: A total of 743 APs were recruited in this large-scale biomarker identification study with a nested case-control study. Data from the discovery cohort demonstrated that 90 exosomal miRNAs were significantly dysregulated in APs with POF compared with controls. One miRNA classifier (Cmi) comprising 3 miRNAs (miR-4265, 1208, 3127-5p) was identified in the training cohort, and was further evaluated in two validation cohorts for their predictive value for POF. AUCs for Cmi ranged from 0.88 to 0.90, which was statistically superior to AUCs of APACHE-II and BISAP, and outperformed BUN and creatinine in POF prediction across all cohorts (P<.05). Higher levels of Cmi indicated increased need for ICU admission, prolonged hospitalization, and elevated mortality rate, thus poor prognosis. In the nested case-control study, Cmi could help identify prediagnostic POF in post-ERCP pancreatitis cases within "golden hours" after ERCP with high efficacy. CONCLUSIONS: Serum exosomal Cmi may be an early predictor for POF in AP, even within "golden hours" after AP onset. TRIAL REGISTRATION: ClinicalTrials.gov (NCT02602808).

Low RCAN1.2 mRNA Expression Is Associated with Poor Prognosis of Patients with Esophageal Squamous Cell Carcinoma.

Background: This study aims to explore the role of RCAN1 in esophageal squamous cell carcinoma (ESCC) cells, determine the mRNA level of three RCAN1 isoforms in ESCC tissue, and evaluate the prognostic value of three RCAN1 isoforms. Methods: Colony-forming assay, Wound-healing assay and Transwell assay were used to evaluate the effect of RCAN1 on cell proliferation, migration and invasion. The mRNA expression of three RCAN1 isoforms was detected in paired tumor and normal tissues from 100 ESCC patients by real-time PCR. Kaplan-Meier survival curves and Cox proportional hazards model were used to evaluate the prognostic value of three RCAN1 isoforms. A nomogram was used to predict the probability of 2-year and 5-year overall survival (OS). Results: In vitro, knockdown of RCAN1 could promote ESCC cell proliferation, migration and invasion abilities. Compared to the paired normal tissues, RCAN1 isoform 1 (RCAN1.1, P=0.0027) and RCAN1 isoform 2 (RCAN1.2, P=0.0006) were significantly decreased in tumor tissues. The low expression of RCAN1.2 mRNA was associated with advanced stage (P=0.0176) and lymph node metastasis (LNM, P=0.0219). ESCC patients with low RCAN1.2 mRNA levels had shorter survival time compared to those with high RCAN1.2 levels (P=0.007). Multivariate COX analysis indicated that RCAN1.2 mRNA level was an independent prognostic indicator of OS of patients with ESCC (hazard ratio=0.5266, P=0.03554). The concordance index of nomogram to predict OS was 0.693 based on LNM, RCAN1.2, tumor stage and patients' age. Conclusion: These findings show that RCAN1 gene play a role in preventing proliferation, migration, and invasive activity of ESCC cells. RCAN1.2 mRNA level is a novel prognostic marker in ESCC, targeting RCAN1.2 may provide a potential therapeutic approach in ESCC.

The Contributions of Trace Elements on Molecular Subtype-Specific Colorectal Cancer.

Purpose: Although growing studies have reported the disturbances of trace elements (TEs) homeostasis was closely associated with the occurrence of colorectal cancer (CRC), the clinical value of TEs in CRC with different molecular subtypes was largely unknown. This study aimed to explore the correlation between KRAS mutations/MSI status and serum TEs levels in patients with CRC. Methods: The serum concentrations of 18 TEs were detected by inductively coupled plasma emission spectrometry (ICP-MS). MSI status (two mononucleotides: BAT25, BAT26, three dinucleotides: D2S123, D5S346, and D17S250), KRAS (G516T, G517A, G518C, G520T, G521A, G522C, and G532A) mutations were detected by the multiplex fluorescent PCR and the real-time fluorescent quantitative PCR, respectively. The correlations among KRAS mutations/MSI status, demographic and clinical characteristics, and TEs were analyzed by Spearman correlation analysis. Results: The propensity score matching (PSM) analysis was adopted to minimize differences between groups. Before PSM, 204 CRC patients were recruited in this study, including 123 KRAS-negative patients and 81 KRAS-positive patients according to the test results of KRAS mutations, and 165 MSS patients and 39 MSI patients based on MSI detection. After PSM, the serum concentration of Mn was significantly lower in CRC patients with KRAS mutations than those without KRAS mutations, and a significant negative correlation was observed between Mn and Pb in the KRAS-positive cases. CRC patients carrying MSI had a significantly lower level of Rb compared to MSS patients. Importantly, Rb was significantly positively correlated with Fe, Mn, Se, and Zn in patients with MSI. Collectively, all our data indicated that the occurrence of different molecular events might be accompanied by different alterations in types and levels of serum TEs. Conclusions: CRC patients with different molecular subtypes presented different alterations in types and levels of serum TEs. Mn was significantly negatively correlated with the KRAS mutations, and Rb was noticeably negatively correlated with the MSI status, indicating certain TEs might contribute to the pathogenesis of molecular subtype-specific colorectal cancer.

Differences of genomic alterations and heavy metals in non-small cell lung cancer with different histological subtypes.

PURPOSE: This study aimed to explore the correlations among heavy metals concentration, histologic subtypes and molecular characteristics in patients with non-small cell lung cancer (NSCLC). METHODS: In this study, an NGS panel of 82 tumor-associated genes was used to identify genomic alternations in 180 newly diagnosed patients with NSCLC. The concentrations of 18 heavy metals in the serum samples were detected by inductively coupled plasma emission spectrometry (ICP-MS). RESULTS: A total of 243 somatic mutations of 25 mutant genes were identified in 115 of 148 patients with LUAD and 45 somatic mutations of 15 mutant genes were found in 24 of 32 patients with LUSC. The genomic alternations, somatic interactions, traditional serum biomarkers, and heavy metals were markedly different between patients with LUAD and LUSC. Moreover, patients with LUSC were significantly positively correlated with Ba, but not LUAD. Lastly, patients with EGFR mutations presented significant negative correlations with Cd and Sr, whereas patients with TP53 mutations showed a significant positive correlation with Pb. CONCLUSION: The genomic alternations, somatic interactions, traditional serum biomarkers, and heavy metals were different between patients with LUAC and LUSC, and heavy metals (e.g., Ba, Pb, and Cd) may contribute to the tumorigenesis of NSCLC with different histological and molecular subtypes.

Effects of N6-Methyladenosine Regulators on LAG3 and Immune Infiltrates in Lung Adenocarcinoma.

Background: Lung adenocarcinoma (LUAD) is the most common histological subtype of lung cancer, which is one of the most commonly diagnosed tumors and the leading causes of death from cancer around the world. Since RNA methylation is a posttranscriptional modification and affects so much biological progress, it is urged to explore the role of N6-methyladenosine (m6A) methylation in LUAD. Methods: We explored the expression of 24 m6A methylation genes, as well as their correlations with LAG3 in 561 LUAD samples from TCGA. Consensus clustering was applied to m6A methylation genes, and two LUAD subgroups were identified. The expression of m6A genes was analyzed by the Wilcoxon test. KEGG and GO enrichment analyses were performed to indicate the pathway affected by differentially expressed genes in the two groups. A prognostic model based on LASSO regression using an eleven-m6A gene signature was constructed according to the expression of these genes. Receiver operating characteristic (ROC) curve was used to confirm the accuracy of the model in the TCGA cohort, as well as in the test cohort from the Gene Expression Omnibus (GEO) database. Results: Compared to cluster 1, cluster 2 showed poorer overall survival (OS) and higher LAG3 expression. In addition, KEGG and GO enrichment analyses indicated that differentially expressed genes are enriched in the immune response. We also observed that the expression of LAG3 is positively correlated with IGF2BP2, CBLL1, and HNRNPA2B1 and negatively correlated with YTHDF2, YTHDF3, and FTO. For patients in the TCGA cohort, the AUC score is 0.7, and the AUC score for the GSE50081 cohort is 0.675. Patients with lower risk scores exhibited better overall survival and lower expression of LAG3 than patients with higher risk scores. Conclusions: In brief, our results indicated the important role of m6 methylation in affecting the tumor immune microenvironment and the survival of patients with LUAD. The m6A methylation gene signatures might serve as promising therapeutic targets and help the immunotherapy of LUAD in the future.

Differences in actionable genomic alterations between brain metastases and non­brain metastases in patients with non­small cell lung cancer.

Brain metastases (BM) have been closely associated with increased morbidity and poor survival outcomes in patients with non­small cell lung cancer (NSCLC). Excluding risk factors in histological subtypes, genomic alterations, including epidermal growth factor receptor mutations and anaplastic lymphoma kinase (ALK) rearrangements have been also regarded as greater risk factors for BM in the aspect of molecular subtypes. In the present study, 69 tumor tissues and 51 peripheral blood samples from patients with NSCLC were analyzed using a hybridization capture­based next­generation sequencing (NGS) panel, including 95 known cancer genes. Among the 90 patients with stage IV NSCLC, 26 cases suffered from BM and 64 cases did not. In total, 174 somatic mutations in 35 mutated genes were identified, and 12 of these genes were concurrently present in the BM group and the non­BM group. Importantly, five mutated genes including ALK, cytidine deaminase (CDA), SMAD family member 4 (SMAD4), superoxide dismutase 2 (SOD2) and Von Hippel­Lindau tumor suppressor (VHL) genes were uniquely detected in the BM group, and they were enriched in the Hippo signaling pathway, pyrimidine metabolism and pantothenate and co­enzyme A (CoA) biosynthesis, as demonstrated using Kyoto Encyclopedia of Genes and Genomes enrichment analysis. RNA polymerase II transcription regulator complex and promyelocytic leukemia nuclear body were the top functional categories according to the Gene Ontology enrichment analysis in the BM group and non­BM group, respectively. Furthermore, 43.33% (13/30) of mutated genes were detected by both tumor tissue deoxyribonucleic acid (DNA) and plasma­derived circulating tumor DNA (ctDNA) in the non­BM group, while this percentage was only limited to 29.41% (5/17) in the BM group. To summarize, significant differences in somatic mutations, somatic interactions, key signaling pathways, functional biological information, and clinical actionability for the therapy of targeted agents were founded between the BM group and the non­BM group, and ctDNA analysis may by applied as a more credible alternative for genomic profiling in patients with advanced NSCLC without BM, due to its higher consistency for genomic profiling between ctDNA analysis and tissue DNA analysis.

Mutational landscape and potential therapeutic targets for sporadic pancreatic neuroendocrine tumors based on target next-generation sequencing.

Pancreatic neuroendocrine tumor (PNET), a heterogenous type of neoplasm with limited treatment options, is relatively rare and to date, the genetic background has remained to be fully elucidated. The present study aimed to determine the mutational landscape of PNET with and without liver metastasis, as well as its clinical application value for treatment. Fresh tumor tissues were collected from 14 patients with PNET following surgery, 4 of whom had developed liver metastasis. Subsequently, targeted next-generation sequencing of 612 cancer-associated genes and comprehensive analysis were performed on the tumor tissues. The results identified 63 somatic mutations in 53 genes in the 14 patients with PNET, amongst which menin 1 was identified as the most recurrently mutated gene. The analysis also identified several novel recurrently mutated genes, including adrenoceptor alpha 2B, ARVCF delta catenin family member, carbamoyl-phosphate synthetase 2, aspartate transcarbamylase, and dihydroorotase and neuregulin 1. Among the 53 mutated genes, 11 were enriched in the PI3K/AKT signaling pathway (adjusted P=7.12x10-5). In addition, 4 patients with PNET with liver metastasis had distinctly different mutational profiles compared with those without liver metastasis; 13 genes were discovered to be exclusively mutated in the liver metastasis group of the patients with PNET, including ATRX chromatin remodeler, thioredoxin reductase 2, anus kinase 3, ARVCF delta catenin family member, integrin subunit alpha V and RAD50 double strand break repair protein. In addition, two potentially actionable alterations in BRCA2 DNA repair-associated (p.Q548Q) and neurofibromin 1 (p.Q1188X) were identified using the OncoKB database. In conclusion, the present study generated a comprehensive mutational profile of 14 patients with PNET and further described the features of patients with liver metastasis, which highlights potential targets for drug development of PNET.

Definition and risk factors of early recurrence based on affecting prognosis of esophageal squamous cell carcinoma patients after radical resection.

Early recurrence after surgery could affect cancerous patients' prognosis, but the definition of early recurrence and its risk factors for esophageal squamous cell carcinoma (ESCC) patients are still unclear. This study analyzed the clinical data of 468 post-surgery recurrent ESCC patients retrospectively. A minimum p-value approach was used to evaluate the optimal cut-off value of recurrence free survival (RFS) to define early recurrence. Risk factors of early recurrence were developed based on a Cox model. The optimal cut-off value of RFS to distinguish early recurrence was 21 months (p <0.001). Independent risk factors for early recurrence included tumor locations (HR=0.562, p <0.001), pathological T stage (HR=1.829, p <0.001), tumor diameter (HR=1.344, p = 0.039), positive lymph nodes (HR=1.361, p <0.001), and total resected lymph nodes (HR=1.271, p = 044). For the late recurrent patients, there was a much more significant survival advantage for recurrence after concurrent chemoradiotherapy than that after sequential chemoradiotherapy and radiotherapy alone (p = 0.0066). In conclusion, this study defined 21 months of RFS as early recurrence and also identified its risk factors. Concurrent chemoradiotherapy was suggested as preferred post-relapse treatment for late recurrent ESCC patients.

MicroRNA-429 acts as a tumor suppressor in colorectal cancer by targeting high mobility group box 3.

Colorectal cancer (CRC) is one of the most common solid tumors worldwide and has an extremely poor prognosis. MicroRNA-429 (miR-429) has been reported to participate in the progression of CRC. However, the pathological mechanisms require further investigation. The aim of the present study was to investigate the association between miR-429 and high mobility group box 3 (HMGB3) in CRC and the associated mechanism. The mRNA expression levels of miR-429 and HMGB3 in 65 paired CRC and adjacent tissues were examined by reverse transcription-quantitative PCR. Furthermore, a dual-luciferase reporter assay was performed to identify the association between miR-429 and HMGB3. Finally, the effects of miR-429 and HMGB3 on the proliferation and apoptosis of CRC cells were detected. As a result, it was identified that miR-429 expression was downregulated and HMGB3 expression was upregulated in CRC tissues compared with in adjacent non-cancer tissues, and the expression levels of miR-429 were negatively associated with those of HMGB3. Notably, HMGB3 was demonstrated to be a direct target of miR-429 by dual-luciferase reporter assay. Furthermore, transfection with a miR-429 mimic significantly inhibited HMGB3 expression and led to decreased proliferation and increased apoptosis of CRC cells. On the other hand, transient overexpression of HMGB3 partially inhibited the antitumor effects of miR-429. To the best of our knowledge, the present study demonstrated for the first time that miR-429 regulated the proliferation and apoptosis of CRC cells via HMGB3, suggesting a specific tumor suppressive function of the miR-429/HMGB3 signaling pathway in CRC.

Clinical Significance of Monitoring Circulating Free DNA and Plasma Heat Shock Protein 90alpha in Patients with Esophageal Squamous Cell Carcinoma.

BACKGROUND: Esophageal squamous cell carcinoma (ESCC) is the predominant histological type of esophageal cancer in China and has an extremely poor prognosis. Circulating free DNA (cfDNA) and plasma heat shock protein 90alpha (Hsp90a) are two novel noninvasive biomarkers for diagnosis and prognostic prediction of several types of cancer. However, to the best of our knowledge, the roles of the two biomarkers in ESCC are still unknown. METHODS: Here, we recruited 93 primary ESCC patients and detected plasma concentrations of the two markers at different time points, including 1-3 days pre-chemotherapy, 1-7 days pre-surgery and 7-14 days post-surgery. Baseline concentrations of the two markers were associated with main characteristics of ESCC patients which were collected at first diagnosis. Correlation between the two markers and traditional serum biomarkers at baseline was also examined. Furthermore, dynamic changes of the cfDNA and Hsp90α concentrations among different time points and the potential clinical significance were assessed. RESULTS: Consequently, there was no significant association between baseline concentrations of the two markers and clinical features. Especially, cfDNA demonstrated stronger correlation with other circulating biomarkers than Hsp90α at baseline level. Importantly, both cfDNA and Hsp90α concentrations were significantly increased after surgery. Kaplan-Meier survival analysis showed that a change in concentration of cfDNA (ΔcfDNA) but not Hsp90α (ΔHSP90ɑ) between pre-surgery and post-surgery had significant effect on the overall survival of surgical patients with ESCC. CONCLUSION: Thus, ΔcfDNA evaluation could be a promising prognostic marker for surgical ESCC patients. Our findings may improve the understanding of the function of cfDNA and Hsp90α in ESCC.

The Landscape of Actionable Genomic Alterations by Next-Generation Sequencing in Tumor Tissue Versus Circulating Tumor DNA in Chinese Patients With Non-Small Cell Lung Cancer.

Background: Circulating tumor DNA (ctDNA) sequence analysis shows great potential in the management of non-small cell lung cancer (NSCLC) and the prediction of drug sensitivity or resistance in many cancers. Here, we drew and compared the somatic mutational profile using ctDNA and tumor tissue sequence analysis in lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC), and assess its potential clinical value. Methods: In this study, 221 tumor tissues and 174 plasma samples from NSCLC patients were analyzed by hybridization capture-based next-generation sequencing (NGS) panel including 95 cancer-associated genes. Tumor response assessments were applied to 137 patients with advanced-stage (III and IV) NSCLC who first received targeted agents. Results: Twenty significantly mutated genes were identified such as TP53, EGFR, RB1, KRAS, PIK3CA, CD3EAP, CTNNB1, ERBB2, APC, BRAF, TERT, FBXW7, and HRAS. Among them, TP53 was the most frequently mutated gene and had a higher mutation probability in male (p = 0.00124) and smoking (p < 0.0001) patients. A total of 48.35% (191/395) of NSCLC patients possessed at least one actionable alteration according to the OncoKB database. Although the sensitivity of genomic profiling from ctDNA was lower than that from tumor tissue DNA, the mutational landscape of target genes from ctDNA is similar to that from tumor tissue DNA, which led to 61.22% (30/49) of mutational concordance in NSCLC. Additionally, the mutational concordance between tissue DNA and ctDNA in LUAD differs from that in LUSC, which is 63.83% versus 46.67%, indicating that NSCLC subtypes influence the specificity of mutation detection in plasma-derived ctDNA. Lastly, patients with EGFR and TP53 co-alterations showed similar responses to Gefitinib and Icotinib, and the co-occurring TP53 mutation was most likely to be a poor prognostic factor for patients receiving Gefitinib, indicating that the distributions and types of TP53 mutations may contribute to the efficacy and prognosis of molecular targeted therapy. Conclusions: As a promising alternative for tumor genomic profiling, ctDNA analysis is more credible in LUAD than in LUSC. Genomic subtyping has strong potential in prognostication and therapeutic decision-making for NSCLC patients, which indicated the necessity for the utility of target NGS in guiding clinical management.

RNA Sequencing Revealed Signals of Evolution From Gallbladder Stone to Gallbladder Carcinoma.

Gallbladder stone is a major risk factor for gallbladder carcinoma (GBC), while there is still a controversy whether period of follow-up since newly diagnoses of asymptomatic gallstones increases the risk of GBC. In this study, 10 GBC patients and 30 patients with gallstones were admitted to our hospital. Patients with gallstones were divided into 3 groups according to the follow-up time, involving 10 patients with follow-up period of 1-3 years (GS3 group), 10 patients with follow-up period of 5-10 years (GS5 group), and 10 patients with follow-up period of more than 10 years (GS10 group). Tumor and para-tumor tissues of GBC patients, and gallbladder tissues of gallstone patients were collected. RNA sequencing was performed on the 50 samples. Besides, 1,704 differentially expressed genes (DEGs) were identified in tumors compared with para-tumor tissues of 10 GBC patients, which were enriched into some well-known cancer-related pathways, such as PI3K-Akt, mitogen-activated protein kinase (MAPK), Ras, and Wnt signaling pathways, and the most significant pathway was neuroactive ligand-receptor interaction. Patients with gallstones with periods of follow-up equal to 1-3 and > 10 years showed to have higher cancer risk than those with 5-10 years. ALPP and GPR87 are potential biomarkers for predicting cancer risk in patients with gallstones. The in vitro results revealed that GPR-87 can promote the proliferation, migration, and invasion of GBC cells. Herein, we explored the relationship between GBC patients and patients with gallstones with different periods of follow-up in transcriptome level.

Exosomal miRNAs are potential diagnostic biomarkers between malignant and benign thyroid nodules based on next-generation sequencing.

An accurate biomarker or method for diagnosis of thyroid nodule with indeterminate fine-needle aspiration result is essential for clinical treatment. Micro RNAs (miRNAs) of exosomes are advantageous in the diagnosis of tumors because they are highly stable, and be protected by a bilayer membrane structure. Exosomes were isolated from 13 papillary thyroid carcinoma (PTC) and 7 nodular goiter (NG) patients' plasma. Small RNA sequencing was performed on exosomes' RNA in next-generation sequencing (NGS) platform. Then, we performed comprehensive analysis on miRNA expression profile in exosome of two groups. One hundred and twenty-nine differentially expressed miRNAs were identified in plasma exosomes between PTC and NG patients. Forty-nine miRNAs were up-regulated, and 80 miRNAs were down-regulated in PTC patients. Receiver operating characteristic (ROC) curves of 129 miRNAs were plotted. Area under curve (AUC) of 129 miRNAs was 0.571-0.951, with distribution peak of 0.82-0.86. AUC of 11 miRNAs was above 0.9, miR-5189-3p had the most optimal performance for diagnosis between PTC and NG, with 0.951 of AUC. Target genes of 129 miRNAs were enriched into 7 cancer-related signaling pathways, including mitogen-activated protein kinase (MAPK), tumor necrosis factor (TNF), NF-kappa B signaling pathway and so on. This study profiled the miRNA signature of exosomes from PTC patients and NG patients. We proposed a group of miRNAs in plasma exosomes as candidate biomarkers for thyroid nodule diagnosis.

Genetic Profiling of Breast Cancer with and Without Preexisting Metabolic Disease.

Breast cancer is the most commonly diagnosed cancer and the leading cause of cancer death among women. Various mechanisms are involved in the initiation and progression of breast cancer. Metabolic dysregulation has been associated with increasing breast cancer incidence and mortality. However, little is known about how metabolic disease regulates the development and progression of breast cancer at the molecular level. Here, using a hybridization capture-based panel including 124 cancer-associated genes, we performed targeted next-generation sequencing of tumor tissues and matched blood samples from 20 postmenopausal patients with primary breast cancer, in which 6 cases suffered from preexisting metabolic disorders including hypertension, type 2 diabetes, and coronary heart disease. We took only the protein-altering variants and identified 170 somatic mutations of 59 genes. Among these, 40 mutated genes were found in the metabolic disease group, and 33 mutated genes were found in the non-metabolic disease group. Importantly, nonsynonymous mutations of 26 genes (MSH3, BRAF, MLH3, MTOR, DDR2, ALK, etc.) were uniquely present in the metabolic disease group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were performed to investigate biological functions and key pathways of somatic mutations. TP53, PIK3CA, and PTEN were the top three commonly mutated genes at a higher frequency compared with the Cancer Genome Atlas (TCGA) data, and several novel but infrequent mutations in other genes were also found. Although further studies are required to validate these variants, our results are the first to suggest a specific molecular profile of breast cancer with preexisting metabolic disease.

Detection of ctDNA in the plasma of patients with papillary thyroid carcinoma.

A total of 15-30% of thyroid nodules that are evaluated by fine-needle aspiration are not clearly determined to be benign or malignant. Gene mutation analysis is recommended for the evaluation of thyroid nodules using clinical guidelines. The detection of circulating tumor DNA (ctDNA) has the potential to aid in the screening, diagnosis and prediction of thyroid cancer prognosis, and can be used when tissues are difficult to obtain. In the present study, whole-exome sequencing (WES) was performed on tumors and matched normal tissues from 10 patients with papillary thyroid carcinoma (PTC) in Quanzhou, China. Hotspot mutations in tumor DNA and cell-free DNA were identified in the validation cohort, which included 59 patients with PTC. BRAF V600E occurred in five samples, and was the most frequent mutation observed. Variation allele frequency (VAF) of BRAF V600E detected by WES was positively correlated with VAF determined using digital PCR (R2=0.9197; P=0.0099). A number of novel mutated genes were identified, including zinc finger protein 717, pleckstrin homology like domain family A member 1, RBMX like 3, lysine methyltransferase 5A and trichohyalin, along with the reported genes BRAF, NRAS and mucin 16, cell surface associated. Somatic mutated genes were enriched in the 'focal adhesion' pathway, as determined by Kyoto Encyclopedia of Genes and Genomes or Gene Ontology analysis. In the validation cohort, 44.07% of tumors were BRAF V600E-positive, and the sensitivity and specificity of BRAF V600E ctDNA were 61.54 and 90.91%, respectively. BRAF V600E was associated with aggressive tumor factors, including lymph node metastasis (P=0.001) and advanced disease stage (P=0.009). The present study investigated the accuracy of ctDNA detection in patients with PTC, and provided evidence that ctDNA can be used as an evaluation of tumor DNA in thyroid nodules.

Whole exome sequencing study of a Chinese concurrent cancer family.

Cancer is one of the leading causes of mortality in China, and poses a threat to public health due to its increasing incidence and mortality rates. Concurrent cancer is defined as one or more organs in the same individual having ≥2 primary malignancies occurring simultaneously or successively; however, concurrent cases are rare and poorly studied. The present study recruited a Chinese family presenting multiple cases of concurrent cancer and performed whole exome sequencing in one unaffected and two affected individuals to identify the causative mutations. DNA was extracted from peripheral blood and tumor tissue samples. Following an exome capture and quality test, the qualified library was sequenced as 100 bp paired-end reads on an Ion Torrent platform. Clean data were obtained by filtering out the low-quality reads. Subsequently, bioinformatics analyses were performed using the clean data. After mapping and annotating in 1000 Genomes Project database, the existing SNP database and the Cancer Gene Census (CGC) database, it was revealed that the NADH:ubiquinone oxidoreductase core subunit S7 gene was a candidate gene with somatic mutations, and a subset of 16 genes were candidate genes with germline mutations. The findings of the present study may improve the understanding of the molecular pathogenesis of concurrent cancer.

Immune-related somatic mutation genes are enriched in PDACs with diabetes.

The bidirectional interaction between pancreatic cancer (PanCa) and diabetes has been confirmed by epidemiological studies, but until now, the underlying molecular mechanisms for this connection is not fully understood yet. Here, we analyzed the clinical and genomic data of 26 pancreatic ductal adenocarcinoma (PDAC) patients without diabetes, and six diabetic PDAC patients, whose tumors underwent targeted next-generation sequencing (551 cancer-related genes included). Ingenuity Pathway Analysis (IPA) was performed to investigate genetic alterations and biological consequences associated with PDACs with or without diabetes. We identified 345 somatic mutations of 153 genes in test cohort and a positive association between diabetes duration and somatic mutation burden. KRAS, TP53, and SMAD4 were the top3 commonly mutated genes at a similar frequency compared to the Cancer Genome Atlas (TCGA) data. Several novel but infrequent mutations in other genes (MDC1, PRB2, and PRB4) were also found. Besides, 13 mutated genes (PIK3CD, SNCAIP, IRF4, HLA-A, NOTCH4, PIM1, ETV6, B2M, CD70, PRDM14, TGFBR1, FLT1, and PARP2) were uniquely found in the diabetic group, mainly involved in immune-related pathways. Further targeted sequencing of these genes in an independent validation cohort (n = 50) revealed significant enrichment in the diabetic group (n = 18, P = 2.6964E-08). Long-standing diabetes (≥3-year duration) may induce increasing somatic mutations with time, facilitating tumor initiation. Gene mutants associated with immune-related pathways could be used to distinguish the diabetic PDAC patients from the non-diabetic cases and allow more selective treatment.

Mutational landscape of gastric cancer and clinical application of genomic profiling based on target next-generation sequencing.

BACKGROUND: Gastric cancer (GC) is a leading cause of cancer deaths, and an increased number of GC patients adopt to next-generation sequencing (NGS) to identify tumor genomic alterations for precision medicine. METHODS: In this study, we established a hybridization capture-based NGS panel including 612 cancer-associated genes, and collected sequencing data of tumors and matched bloods from 153 gastric cancer patients. We performed comprehensive analysis of these sequencing and clinical data. RESULTS: 35 significantly mutated genes were identified such as TP53, AKAP9, DRD2, PTEN, CDH1, LRP2 et al. Among them, 29 genes were novel significantly mutated genes compared with TCGA study. TP53 is the top frequently mutated gene, and tends to mutate in male (p = 0.025) patients and patients whose tumor located in cardia (p = 0.011). High tumor mutation burden (TMB) gathered in TP53 wild-type tumors (p = 0.045). TMB was also significantly associated with DNA damage repair (DDR) genes genotype (p = 0.047), Lauren classification (p = 1.5e-5), differentiation (1.9e-7), and HER2 status (p = 0.023). 38.31% of gastric cancer patients harbored at least one actionable alteration according to OncoKB database. CONCLUSIONS: We drew a comprehensive mutational landscape of 153 gastric tumors and demonstrated utility of target next-generation sequencing to guide clinical management.

Cell adhesion-related gene somatic mutations are enriched in aggressive papillary thyroid microcarcinomas.

BACKGROUND: Approximately half of the documented increases in differentiated thyroid carcinoma is due to identification of papillary thyroid microcarcinomas (PTMCs). Knowing whether PTMC is aggressive is required for proper treatment, but until now, there has been no method for assessing these traits and understanding the underlying mechanisms for aggressiveness. METHODS: We performed whole-exome sequencing of 16 PTMCs and matched normal thyroid tissues and GO/KEGG analysis to study genetic alterations and biological consequences associated with aggressive PTMCs, and then sequenced these genes using a next-generation gene-panel approach in an additional 70 PTMC samples including aggressive (n = 50) and non-aggressive (n = 20) groups. RESULTS: We identified 254 somatic mutations of 234 genes, for which 178 mutations in 168 genes were found in the aggressive group, and 76 mutations in 74 genes were found in the non-aggressive group. Several recurrent mutations in BRAF, VCAN, ALDH1L1, and MUC5B were identified, and many novel but infrequent mutations in other genes were also found. The aggressive cohort had more mutational burdens than the non-aggressive group (P = 0.004). Nonsynonymous mutations of 13 genes (MUC5B, TNN, SSPO, PPFIA1, PCDHGA2, ITGA8, ITGA4, DCHS1, CRNN, ROCK1, RELN, LAMC2, and AEBP1) were involved in cell adhesion, and these were only present in the aggressive group. Targeted sequencing of these genes revealed significant enrichment in the aggressive group (P = 0.000004). CONCLUSION: PTC may have evolved from PTMC due to sharing similar gene mutations, and the accumulation of such mutations promoted the aggressiveness of PTMC. Gene mutants associated with cell adhesion may be used to predict PTMC aggressiveness and allow more selective treatment.

Identification of recurrent and novel mutations by whole­genome sequencing of colorectal tumors from the Han population in Shanghai, eastern China.

Previous studies have identified recurrent oncogenic mutations in colorectal cancer (CRC), but there is limited CRC genomic data from the Chinese Han population. Whole­genome sequencing was performed on 10 primary CRC tumors and matched adjacent normal tissues from patients from the Han population in Shanghai, at an average of 27.8x and 27.9x coverage, respectively. In the 10 tumor samples, 32 significant somatic mutated genes were identified, 13 of which were also reported as CRC mutations in The Cancer Genome Atlas Network. All the mutated genes were enriched in functions associated with channel activity, which has rarely been reported in previous studies investigating CRC. Furthermore, 21 chromosomal rearrangements were detected and 4 rearrangements encoded predicted in­frame fusion proteins, including a fusion of phosphorylase kinase regulatory subunit b and NOTCH2 demonstrated in 2 out of 10 tumors. Chromosome 8 was amplified in 1 tumor and chromosome 20 was amplified in 2 out of 10 CRC patients. The present study produced a genomic mutation profile of CRC, which provides a valuable resource for further insight into the mutations that characterize CRC in patients from the Han population in Shanghai, eastern China.