BAP1 is a novel regulator of HIF-1α.

BAP1 is a powerful tumor suppressor gene characterized by haplo insufficiency. Individuals carrying germline BAP1 mutations often develop mesothelioma, an aggressive malignancy of the serosal layers covering the lungs, pericardium, and abdominal cavity. Intriguingly, mesotheliomas developing in carriers of germline BAP1 mutations are less aggressive, and these patients have significantly improved survival. We investigated the apparent paradox of a tumor suppressor gene that, when mutated, causes less aggressive mesotheliomas. We discovered that mesothelioma biopsies with biallelic BAP1 mutations showed loss of nuclear HIF-1α staining. We demonstrated that during hypoxia, BAP1 binds, deubiquitylates, and stabilizes HIF-1α, the master regulator of the hypoxia response and tumor cell invasion. Moreover, primary cells from individuals carrying germline BAP1 mutations and primary cells in which BAP1 was silenced using siRNA had reduced HIF-1α protein levels in hypoxia. Computational modeling and co-immunoprecipitation experiments revealed that mutations of BAP1 residues I675, F678, I679, and L691 -encompassing the C-terminal domain-nuclear localization signal- to A, abolished the interaction with HIF-1α. We found that BAP1 binds to the N-terminal region of HIF-1α, where HIF-1α binds DNA and dimerizes with HIF-1ß forming the heterodimeric transactivating complex HIF. Our data identify BAP1 as a key positive regulator of HIF-1α in hypoxia. We propose that the significant reduction of HIF-1α activity in mesothelioma cells carrying biallelic BAP1 mutations, accompanied by the significant reduction of HIF-1α activity in hypoxic tissues containing germline BAP1 mutations, contributes to the reduced aggressiveness and improved survival of mesotheliomas developing in carriers of germline BAP1 mutations.

Circulating CD8 lymphocytes predict response to atezolizumab-bevacizumab in hepatocellular carcinoma.

Due to the lack of biomarkers predictive of response to atezolizumab-bevacizumab, the standard of care for advanced HCC, we analyzed baseline and early on-treatment variation of peripheral lymphocyte populations of 37 prospective patients treated by atezolizumab-bevacizumab and in 15 prospective patients treated by sorafenib or lenvatinib (TKIs). RNAseq analysis followed by RT-PCR validation on patients-derived PBMC was also performed. At first imaging, re-evaluation 13 patients receiving atezolizumab-bevacizumab, showed an objective response, 17 stable disease, while 7 were nonresponders. Baseline CD8+ and CD8+PD-L1+ peripheral lymphocytes were lower in responders versus nonresponders (T-test, p = 0.012 and 0.004, respectively). At 3 weeks, 28 of 30 responders displayed a rise of CD8+PD1+ lymphocytes with a positive mean fold change of 4.35 (±5.6 SD), whereas 6 of 7 nonresponders displayed a negative fold change of 0.89 (±0.84 SD). These changes were not observed in patients treated by TKIs. TRIM56, TRIM16, TRIM64, and Ki67 mRNAs were validated as upregulated in responders versus nonresponders after 3 weeks after treatment start, providing possible evidence of immune activation. Baseline CD8+ and CD8+PD-L1+ peripheral lymphocytes and early changes in CD8+PD1+ lymphocytes predict response to atezolizumab-bevacizumab providing noninvasive markers to complement clinical practice in the very early phases of treatment of HCC patients.

Near-Complete Response to Osimertinib for Advanced Non-Small-Cell Lung Cancer in a Pretreated Patient Bearing Rare Compound Exon 20 Mutation (S768I + V774M): A Case Report.

Third-generation tyrosine kinase inhibitors are the first-line gold standard in treating advanced non-small-cell lung cancer bearing common EGFR mutations, but data documenting clinical efficacy in uncommon mutations are currently limited. In this paper, we describe the case of a patient bearing uncommon compound EGFR mutations in exon 20, who experienced a near-complete response to third-line Osimertinib, with metabolic complete response of pulmonary, nodal and ostheolytic lesions. This radiological assessment corresponded to an ECOG PS improvement (from three to one) and a substantial clinical benefit for the patients. Out of two mutations, S768I was associated with poor response to third-generation TKI and V774M had unknown clinical significance, highlighting the complexity of the correct management of these kinds of mutations. We reviewed the literature to document the up-to-date preclinical and clinical data concerning third-generation tyrosine kinase inhibitors for the treatment of patients bearing uncommon EGFR mutations.

Detection of disease-causing mutations in prostate cancer by NGS sequencing.

Gene mutations may affect the fate of many tumors including prostate cancer (PCa); therefore, the research of specific mutations associated with tumor outcomes might help the urologist to identify the best therapy for PCa patients such as surgical resection, adjuvant therapy or active surveillance. Genomic DNA (gDNA) was extracted from 48 paraffin-embedded PCa samples and normal paired tissues. Next, gDNA was amplified and analyzed by next-generation sequencing (NGS) using a specific gene panel for PCa. Raw data were refined to exclude false-positive mutations; thus, variants with coverage and frequency lower than 100× and 5%, respectively were removed. Mutation significance was processed by Genomic Evolutionary Rate Profiling, ClinVar, and Varsome tools. Most of 3000 mutations (80%) were single nucleotide variants and the remaining 20% indels. After raw data elaboration, 312 variants were selected. Most mutated genes were KMT2D (26.45%), FOXA1 (16.13%), ATM (15.81%), ZFHX3 (9.35%), TP53 (8.06%), and APC (5.48%). Hot spot mutations in FOXA1, ATM, ZFHX3, SPOP, and MED12 were also found. Truncating mutations of ATM, lesions lying in hot spot regions of SPOP and FOXA1 as well as mutations of TP53 correlated with poor prognosis. Importantly, we have also found some germline mutations associated with hereditary cancer-predisposing syndrome. gDNA sequencing of 48 cancer tissues by NGS allowed to detect new tumor variants as well as confirmed lesions in genes linked to prostate cancer. Overall, somatic and germline mutations linked to good/poor prognosis could represent new prognostic tools to improve the management of PCa patients.

Human-rat integrated microRNAs profiling identified a new neonatal cerebral hypoxic-ischemic pathway melatonin-sensitive.

Neonatal encephalopathy (NE) is a pathological condition affecting long-term neurodevelopmental outcomes. Hypothermia is the only therapeutic option, but does not always improve outcomes; hence, researchers continue to hunt for pharmaceutical compounds. Melatonin treatment has benefitted neonates with hypoxic-ischemic (HI) brain injury. However, unlike animal models that enable the study of the brain and the pathophysiologic cascade, only blood is available from human subjects. Therefore, due to the unavailability of neonatal brain tissue, assumptions about the pathophysiology in pathways and cascades are made in human subjects with NE. We analyzed animal and human specimens to improve our understanding of the pathophysiology in human neonates. A neonate with NE who underwent hypothermia and enrolled in a melatonin pharmacokinetic study was compared to HI rats treated/untreated with melatonin. MicroRNA (miRNA) analyses provided profiles of the neonate's plasma, rat plasma, and rat brain cortexes. We compared these profiles through a bioinformatics tool, identifying Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways common to HI brain injury and melatonin treatment. After evaluating the resulting pathways and the literature, to validate the method, the key proteins expressed in HI brain injury were investigated using cerebral cortexes. The upregulated miRNAs in human neonate and rat plasma helped identify two KEGG pathways, glioma and long-term potentiation, common to HI injury and melatonin treatment. A unified neonatal cerebral melatonin-sensitive HI pathway was designed and validated by assessing the expression of protein kinase Cα (PKCα), phospho (p)-Akt, and p-ERK proteins in rat brain cortexes. PKCα increased in HI-injured rats and further increased with melatonin. p-Akt and p-ERK returned phosphorylated to their basal level with melatonin treatment after HI injury. The bioinformatics analyses validated by key protein expression identified pathways common to HI brain injury and melatonin treatment. This approach helped complete pathways in neonates with NE by integrating information from animal models of HI brain injury.

Gene Expression Landscape of Chronic Myeloid Leukemia K562 Cells Overexpressing the Tumor Suppressor Gene PTPRG.

This study concerns the analysis of the modulation of Chronic Myeloid Leukemia (CML) cell model K562 transcriptome following transfection with the tumor suppressor gene encoding for Protein Tyrosine Phosphatase Receptor Type G (PTPRG) and treatment with the tyrosine kinase inhibitor (TKI) Imatinib. Specifically, we aimed at identifying genes whose level of expression is altered by PTPRG modulation and Imatinib concentration. Statistical tests as differential expression analysis (DEA) supported by gene set enrichment analysis (GSEA) and modern methods of ontological term analysis are presented along with some results of current interest for forthcoming experimental research in the field of the transcriptomic landscape of CML. In particular, we present two methods that differ in the order of the analysis steps. After a gene selection based on fold-change value thresholding, we applied statistical tests to select differentially expressed genes. Therefore, we applied two different methods on the set of differentially expressed genes. With the first method (Method 1), we implemented GSEA, followed by the identification of transcription factors. With the second method (Method 2), we first selected the transcription factors from the set of differentially expressed genes and implemented GSEA on this set. Method 1 is a standard method commonly used in this type of analysis, while Method 2 is unconventional and is motivated by the intention to identify transcription factors more specifically involved in biological processes relevant to the CML condition. Both methods have been equipped in ontological knowledge mining and word cloud analysis, as elements of novelty in our analytical procedure. Data analysis identified RARG and CD36 as a potential PTPRG up-regulated genes, suggesting a possible induction of cell differentiation toward an erithromyeloid phenotype. The prediction was confirmed at the mRNA and protein level, further validating the approach and identifying a new molecular mechanism of tumor suppression governed by PTPRG in a CML context.

Monitoring Somatic Genetic Alterations in Circulating Cell-Free DNA/RNA of Patients with "Oncogene-Addicted" Advanced Lung Adenocarcinoma: A Real-World Clinical Study.

Liquid biopsy has advantages over tissue biopsy, but also some technical limitations that hinder its wide use in clinical applications. In this study, we aimed to evaluate the usefulness of liquid biopsy for the clinical management of patients with advanced-stage oncogene-addicted non-small-cell lung adenocarcinomas. The investigation was conducted on a series of cases-641 plasma samples from 57 patients-collected in a prospective consecutive manner, which allowed us to assess the benefits and limitations of the approach in a real-world clinical context. Thirteen samples were collected at diagnosis, and the additional samples during the periodic follow-up visits. At diagnosis, we detected mutations in ctDNA in 10 of the 13 cases (77%). During follow-up, 36 patients progressed. In this subset of patients, molecular analyses of plasma DNA/RNA at progression revealed the appearance of mutations in 29 patients (80.6%). Mutations in ctDNA/RNA were typically detected an average of 80 days earlier than disease progression assessed by RECIST or clinical evaluations. Among the cases positive for mutations, we observed 13 de novo mutations, responsible for the development of resistance to therapy. This study allowed us to highlight the advantages and disadvantages of liquid biopsy, which led to suggesting algorithms for the use of liquid biopsy analyses at diagnosis and during monitoring of therapy response.

In chronic lymphocytic leukaemia, SLAMF1 deregulation is associated with genomic complexity and independently predicts a worse outcome.

In a series of 349 patients with chronic lymphocytic leukaemia (CLL), we found lower levels of signalling lymphocytic activation molecule family member 1 (SLAMF1) expression in cases with highly complex karyotypes, as defined by the presence of five or more chromosomal abnormalities (CK5; P < 0·001) and with major chromosomal structural abnormalities (P < 0·001). SLAMF1 downregulation was significantly associated with advanced Binet Stage (P = 0·001), CD38 positivity (P < 0·001), high ß2 -microglobulin levels (P < 0·001), immunoglobulin heavy chain variable region gene (IGHV) unmutated status (P < 0·001), 11q deletion (P < 0·001), tumour protein p53 (TP53) disruption (P = 0·011) and higher risk CLL International Prognostic Index categories (P < 0·001). Multivariate analysis showed that downregulated SLAMF1 levels had independent negative prognostic impact on time-to-first treatment (P < 0·001) and overall survival (P < 0·001).

Molecular biomarkers predicting early development of endometrial carcinoma: A pilot study.

OBJECTIVE: Endometrial carcinoma represents the most common gynaecological cancer and the sixth most frequent cancer among women worldwide. The 5-year survival of patients with stage I endometrial carcinoma is 75%-88% versus 50% for stage III or 15% for stage IV disease. Therefore, early detection could improve survival rates. Specifically, in the most prevalent, type 1 endometrial cancer develops from hyperplastic endometrium. The aim of the study was to evaluate the utility of cancer gene mutations from endometrial biopsies towards predicting synchronous or metachronous development of malignant lesions. The aim of the study was to evaluate whether endometrial biopsies could already carry mutations in cancer genes useful for predicting or anticipating subsequent cancer development. METHODS: Patients with a previous endometrial biopsy negative for cancer, followed by a subsequent biopsy positive for cancer, were included in the study. A fifty cancer genes targeted next-generation sequencing panel were used to investigate mutations in matched non-cancerous and malignant samples. RESULTS: All biopsies from cancer tissues harboured mutations in one or more of the following genes: APC, CTNNB1, FBXW7, HNF1A, KRAS, MTOR, NRAS, PIK3CA, PTEN, RB1 and TP53. Additionally, 50% of the biopsies from matched non-cancerous tissues exhibited mutations in PTEN, KRAS or PIK3CA genes. CONCLUSIONS: These results suggest that detecting pathogenic mutations in oncogenes or tumour suppressor genes in an otherwise benign condition is associated with a risk of developing a malignant disease. Given the identification of mutations several months or years before the appearance of a malignancy, our finding suggests that a closer monitoring of patients who present such molecular alterations in non-cancerous uterine mass is warranted.

In chronic lymphocytic leukaemia with complex karyotype, major structural abnormalities identify a subset of patients with inferior outcome and distinct biological characteristics.

Complex karyotype (CK) is a negative prognostic factor in chronic lymphocytic leukaemia (CLL). However, CK is a heterogeneous cytogenetic category. Unbalanced rearrangements were present in 73·3% of 90 CLL patients with CK (i.e. ≥3 chromosome aberrations in the same clone), and were associated with a shorter overall survival (P = 0·025) and a shorter time to first treatment (P = 0·043) by multivariate analysis. Patients with unbalanced rearrangements presented a distinct mRNA expression profile. In conclusion, CLL patients with unbalanced rearrangements might represent a subset of very high-risk CLL patients with distinct clinical and biological characteristics.

Differential expression of hsa-miR-221, hsa-miR-21, hsa-miR-135b, and hsa-miR-29c suggests a field effect in oral cancer.

BACKGROUND: The theory of field effect suggests that the tumor-adjacent area, besides histopathologically normal, undergoes genetic and epigenetic changes that can eventually affect epithelial homeostasis, predisposing the patient to cancer development. One of the many molecular changes described in cancer are microRNAs (miRNAs), which regulates the expression of important genes during carcinogenesis. Thus, the aim of this study was to investigate the field effect in oral cancer. METHODS: We investigated the differential expression profile of four miRNAs (hsa-miR-221, hsa-miR-21, hsa-miR-135b, and hsa-miR-29c) in cancerous oral tissue, in tumor-adjacent tissue and and in non-cancerous tissue samples from healthy volunteers. RESULTS: Our results showed significant overexpression profiles of all four studied miRNAs in cancerous oral tissue compared to non-cancerous samples, as well as in tumor-adjacent tissue compared to cancer-free tissue. No significant difference was found when comparing the expression profile of cancerous and tissue-adjacent tissue groups. We found a negative correlation between the expression of hsa-miR-21 expression and STAT3 in oral squamous cell carcinoma. CONCLUSION: These results suggest that the tissue adjacent to cancer cannot be considered a normal tissue because its molecular aspects are significantly altered. Our data corroborates the hypothesis of field cancerization.

Wnt signalling modulates transcribed-ultraconserved regions in hepatobiliary cancers.

OBJECTIVE: Transcribed-ultraconserved regions (T-UCR) are long non-coding RNAs which are conserved across species and are involved in carcinogenesis. We studied T-UCRs downstream of the Wnt/ß-catenin pathway in liver cancer. DESIGN: Hypomorphic Apc mice (Apcfl/fl) and thiocetamide (TAA)-treated rats developed Wnt/ß-catenin dependent hepatocarcinoma (HCC) and cholangiocarcinoma (CCA), respectively. T-UCR expression was assessed by microarray, real-time PCR and in situ hybridisation. RESULTS: Overexpression of the T-UCR uc.158- could differentiate Wnt/ß-catenin dependent HCC from normal liver and from ß-catenin negative diethylnitrosamine (DEN)-induced HCC. uc.158- was overexpressed in human HepG2 versus Huh7 cells in line with activation of the Wnt pathway. In vitro modulation of ß-catenin altered uc.158- expression in human malignant hepatocytes. uc.158- expression was increased in CTNNB1-mutated human HCCs compared with non-mutated human HCCs, and in human HCC with nuclear localisation of ß-catenin. uc.158- was increased in TAA rat CCA and reduced after treatment with Wnt/ß-catenin inhibitors. uc.158- expression was negative in human normal liver and biliary epithelia, while it was increased in human CCA in two different cohorts. Locked nucleic acid-mediated inhibition of uc.158- reduced anchorage cell growth, 3D-spheroid formation and spheroid-based cell migration, and increased apoptosis in HepG2 and SW1 cells. miR-193b was predicted to have binding sites within the uc.158- sequence. Modulation of uc.158- changed miR-193b expression in human malignant hepatocytes. Co-transfection of uc.158- inhibitor and anti-miR-193b rescued the effect of uc.158- inhibition on cell viability. CONCLUSIONS: We showed that uc.158- is activated by the Wnt pathway in liver cancers and drives their growth. Thus, it may represent a promising target for the development of novel therapeutics.

Characterisation of peripheral blood mononuclear cell microRNA in early onset psoriatic arthritis.

OBJECTIVES: To evaluate the micro-RNA (miRNA) expression profile in patients with early psoriatic arthritis (PsA) in order to assess the role of miRNAs as potential PsA biomarkers. METHODS: The expression of 723 mature miRNAs in peripheral blood mononuclear cells of early PsA patients in comparison with early-rheumatoid arthritis (ERA) patients and healthy controls (HC) was evaluated using a miRNA microarray. All patients had active disease and were naïve from treatment. The results were validated for a specific miRNA (miR-21-5p) in the entire series of patients plus an additional group of early PsA, ERA and HC using droplet digital PCR. RESULTS: In PsA, microarray analysis revealed a distinct pattern of 19 (vs. HC) and 48 (vs. ERA) deregulated miRNAs (p<0.05). The significant up-regulation of miR-21-5p both in early PsA and in ERA in comparison with HC was validated and confirmed. In PsA, miR-21-5p was found significantly down regulated after 12 weeks of therapy, which significantly correlated with the reduction of DAPSA score. CONCLUSIONS: In early PsA, a 19- (vs. HC) and 48- (vs. ERA) miRNA signature was identified. A differential expression of miRNAs in patients with active disease makes them attractive biomarkers of psoriatic disease. MiR-21-5p was found up-regulated both in early PsA and ERA, a finding which highlights its role in the inflammatory process in general and its potential role as a therapeutic target in different inflammatory disorders. A potential role of miR-21-5p as a response to treatment biomarker in early PsA has been identified.

CCAT2, a novel noncoding RNA mapping to 8q24, underlies metastatic progression and chromosomal instability in colon cancer.

The functional roles of SNPs within the 8q24 gene desert in the cancer phenotype are not yet well understood. Here, we report that CCAT2, a novel long noncoding RNA transcript (lncRNA) encompassing the rs6983267 SNP, is highly overexpressed in microsatellite-stable colorectal cancer and promotes tumor growth, metastasis, and chromosomal instability. We demonstrate that MYC, miR-17-5p, and miR-20a are up-regulated by CCAT2 through TCF7L2-mediated transcriptional regulation. We further identify the physical interaction between CCAT2 and TCF7L2 resulting in an enhancement of WNT signaling activity. We show that CCAT2 is itself a WNT downstream target, which suggests the existence of a feedback loop. Finally, we demonstrate that the SNP status affects CCAT2 expression and the risk allele G produces more CCAT2 transcript. Our results support a new mechanism of MYC and WNT regulation by the novel lncRNA CCAT2 in colorectal cancer pathogenesis, and provide an alternative explanation of the SNP-conferred cancer risk.

Mesenchymal progenitors aging highlights a miR-196 switch targeting HOXB7 as master regulator of proliferation and osteogenesis.

Human aging is associated with a decrease in tissue functions combined with a decline in stem cells frequency and activity followed by a loss of regenerative capacity. The molecular mechanisms behind this senescence remain largely obscure, precluding targeted approaches to counteract aging. Focusing on mesenchymal stromal/stem cells (MSC) as known adult progenitors, we identified a specific switch in miRNA expression during aging, revealing a miR-196a upregulation which was inversely correlated with MSC proliferation through HOXB7 targeting. A forced HOXB7 expression was associated with an improved cell growth, a reduction of senescence, and an improved osteogenesis linked to a dramatic increase of autocrine basic fibroblast growth factor secretion. These findings, along with the progressive decrease of HOXB7 levels observed during skeletal aging in mice, indicate HOXB7 as a master factor driving progenitors behavior lifetime, providing a better understanding of bone senescence and leading to an optimization of MSC performance.

Circulating Non-coding RNA as Biomarkers in Colorectal Cancer.

Recent studies suggested that colorectal cancer influences the types and quantity of nucleic acids - especially microRNAs - detected in the bloodstream. Concentration of circulating (cell-free) microRNAs, and possibly of other non-coding RNAs, could therefore serve as valuable colorectal cancer biomarker and could deliver insight into the disease process. This chapter addresses the recent discoveries on circulating microRNA and long non-coding RNA as diagnostic or prognostic biomarkers in colorectal cancer.

Chromosome aberrations detected by conventional karyotyping using novel mitogens in chronic lymphocytic leukemia: Clinical and biologic correlations.

To clarify whether karyotype aberrations (KA) involving regions not covered by the standard fluorescence in situ hybridization (FISH) panel have independent prognostic relevance, we evaluated KA by conventional cytogenetics in a learning cohort (LC; n = 166) and a validation cohort (VC; n = 250) of untreated chronic lymphocytic leukemia (CLL) patients. In the VC, novel mitogens were used to improve metaphase generation and TP53, NOTCH1, and SF3B1 mutations were assessed. KA undetected by FISH were found in 35 and 35% of the cases in the LC and VC, respectively. In addition to FISH, KA allowed reclassification of 23 and 26% of cases in the LC and VC, respectively, into a higher cytogenetic risk group. By multivariate analysis, both in the LC and VC, KA other than isolated 13q deletion correlated with a shorter time to first treatment (TFT; P < 0.001 and 0.003, respectively), while a complex karyotype predicted a worse overall survival (OS, P = 0.015 and 0.010, respectively). In the VC, where a comprehensive biologic assessment was performed, a shorter TFT was also predicted by stage (P < 0.001), IGHV mutational status (P = 0.05), and del(17p)/TP53 mutations (P = 0.033) while stage (P = 0.023) and del(17p)/TP53 mutations (P = 0.024) independently predicted a shorter OS. FISH results did not independently impact on TFT and OS, in the LC and VC cohorts; this was also the case for NOTCH1 and SF3B1 mutations in the VC. We suggest that in CLL, conventional karyotyping with novel mitogens could be more effective than FISH for the detection of KA allowing for a more precise refinement of prognosis.

Gene expression changes in progression of cervical neoplasia revealed by microarray analysis of cervical neoplastic keratinocytes.

To evaluate the gene expression changes involved in neoplastic progression of cervical intraepithelial neoplasia. Using microarray analysis, large-scale gene expression profile was carried out on HPV16-CIN2, HPV16-CIN3, and normal cervical keratinocytes derived from two HPV16-CIN2, two HPV-CIN3 lesions, and two corresponding normal cervical tissues, respectively. Differentially expressed genes were analyzed in normal cervical keratinocytes compared with HPV16-CIN2 keratinocytes and in HPV16-CIN2 keratinocytes compared with HPV16-CIN3 keratinocytes; 37 candidate genes with continuously increasing or decreasing expression during CIN progression were identified. One of these genes, phosphoglycerate dehydrogenase, was chosen for further characterization. Quantitative reverse transcription-polymerase chain reaction and immunohistochemical analysis confirmed that expression of phosphoglycerate dehydrogenase consistently increases during progression of CIN toward cancer. Gene expression changes occurring during CIN progression were investigated using microarray analysis, for the first time, in CIN2 and CIN3 keratinocytes naturally infected with HPV16. Phosphoglycerate dehydrogenase is likely to be associated with tumorigenesis and may be a potential prognostic marker for CIN progression.

Prediction of response to anti-EGFR antibody-based therapies by multigene sequencing in colorectal cancer patients.

BACKGROUND: The anti-epidermal growth factor receptor (EGFR) monoclonal antibodies (moAbs) cetuximab or panitumumab are administered to colorectal cancer (CRC) patients who harbor wild-type RAS proto-oncogenes. However, a percentage of patients do not respond to this treatment. In addition to mutations in the RAS genes, mutations in other genes, such as BRAF, PI3KCA, or PTEN, could be involved in the resistance to anti-EGFR moAb therapy. METHODS: In order to develop a comprehensive approach for the detection of mutations and to eventually identify other genes responsible for resistance to anti-EGFR moAbs, we investigated a panel of 21 genes by parallel sequencing on the Ion Torrent Personal Genome Machine platform. We sequenced 65 CRCs that were treated with cetuximab or panitumumab. Among these, 37 samples were responsive and 28 were resistant. RESULTS: We confirmed that mutations in EGFR-pathway genes (KRAS, NRAS, BRAF, PI3KCA) were relevant for conferring resistance to therapy and could predict response (p = 0.001). After exclusion of KRAS, NRAS, BRAF and PI3KCA combined mutations could still significantly associate to resistant phenotype (p = 0.045, by Fisher exact test). In addition, mutations in FBXW7 and SMAD4 were prevalent in cases that were non-responsive to anti-EGFR moAb. After we combined the mutations of all genes (excluding KRAS), the ability to predict response to therapy improved significantly (p = 0.002, by Fisher exact test). CONCLUSIONS: The combination of mutations at KRAS and at the five gene panel demonstrates the usefulness and feasibility of multigene sequencing to assess response to anti-EGFR moAbs. The application of parallel sequencing technology in clinical practice, in addition to its innate ability to simultaneously examine the genetic status of several cancer genes, proved to be more accurate and sensitive than the presently in use traditional approaches.