Multiple primary lung cancers (MPLCs) harbour various genetic profiles among the tumours, even from individuals with same non-intrinsic risk factors. Paired mutational analyses were performed to obtain a census of mutational events in MPLC and assess their relationship with non-intrinsic risk factors. Thirty-eight surgical specimens from 17 patients diagnosed as MPLC were used. Extracted DNAs were sequenced for somatic mutations in 409 cancer-associated genes from a comprehensive cancer panel. We statistically analysed the correlation between each driver mutation frequency and non-intrinsic risk factors using Fisher's exact test, and whether genetic mutations occurred concomitantly or randomly in MPLC using an exact test. Comprehensive genetic analyses suggested different mutation profiles in tumours within the same individuals, with some exceptions. EGFR, KRAS, TP53, or PARP1 mutations were concomitantly detected in some MPLC cases. EGFR mutations were significantly more frequent in never or light smokers and females. Concomitant EGFR or KRAS mutations in MPLCs were significantly more frequent than expected by chance (P = .0023 and .0049, respectively) suggesting a more prominent role of non-intrinsic risk factors in EGFR and KRAS mutations than other mutations, which occurred more randomly. Concomitant EGFR or KRAS mutations were particularly prominent in never or light smokers and males.

Multiple myeloma is a plasma cell neoplasm with an extremely variable clinical course. Animal models are needed to better understand its pathophysiology and for preclinical testing of potential therapeutic agents. Hematopoietic cells expressing the hypermorphic Rad50(s) allele show hematopoietic failure, which can be mitigated by the lack of a transcription factor, Mef/Elf4. However, we find that 70% of Mef(-/-)Rad50(s/s) mice die from multiple myeloma or other plasma cell neoplasms. These mice initially show an abnormal plasma cell proliferation and monoclonal protein production, and then develop anemia and a decreased bone mineral density. Tumor cells can be serially transplanted and according to array CGH and whole exome sequencing, the pathogenesis of plasma cell neoplasms in these mice is not linked to activation of a specific oncogene, or inactivation of a specific tumor suppressor. This model recapitulates the systemic manifestations of human plasma cell neoplasms, and implicates cooperativity between the Rad50(s) and Mef/Elf4 pathways in initiating myelomagenic mutations that promote plasma cell transformation.

Alcohol consumption and smoking pose a significant risk for esophageal squamous cell neoplasia (ESCN) development in males; however, ESCN is often diagnosed in non-drinking and non-smoking females. The mechanisms underlying these differences remain elusive, and understanding them can potentially identify novel pathways involved in ESCN development. We performed short-read sequencing to identify somatic variants on a cancer panel targeting 409 genes using DNA extracted from the superficial squamous cell carcinoma (ESCC) tissues and adjacent non-neoplastic epithelium (NE), and immunohistochemical staining of the protein encoded by the target gene. All male patients (n = 117) were drinkers or smokers, whereas 45% of the female patients (n = 33) were not. Somatic variants were compared among three age-matched groups: 13 female ESCC patients with smoking and drinking habits (known-risk group, F-KR), 13 female ESCC patients without these habits (unknown-risk group, F-UR), and 27 males with ESCC and smoking and drinking habits (M-KR). In the NE, the frequencies of CDKN2A variants were significantly higher in F-UR than in F-KR and M-KR. In both ESCC and NE, p14ARF was significantly overexpressed in F-UR than in the other groups. In conclusion, CDKN2A might be important in ESCC development, independent of known risk factors.

Protein phosphorylation is an important type of post-translational modification that is involved in a variety of biological activities. Most phosphorylation events occur on serine, threonine and tyrosine residues in eukaryotes. In recent years, many phosphorylation sites have been identified as a result of advances in mass-spectrometric techniques. However, a large percentage of phosphorylation sites may be non-functional. Systematically prioritizing functional sites from a large number of phosphorylation sites will be increasingly important for the study of their biological roles. This study focused on exploring the intrinsic features of functional phosphorylation sites to predict whether a phosphosite is likely to be functional. We found significant differences in the distribution of evolutionary conservation, kinase association, disorder score, and secondary structure between known functional and background phosphorylation datasets. We built four different types of classifiers based on the most representative features and found that their performances were similar. We also prioritized 213,837 human phosphorylation sites from a variety of phosphorylation databases, which will be helpful for subsequent functional studies. All predicted results are available for query and download on our website (Predict Functional Phosphosites, PFP, http://pfp.biosino.org/).

Multiple myeloma (MM) was characterized by frequent mutations in KRAS/NRAS/BRAF within the EGFR pathway that could induce resistance to EGFR inhibitors. We here report that EGFR inhibition solely exhibited moderate inhibition in KRAS/NRAS/BRAF wildtype (triple-WT) MM cells, whilst had no effect in myeloma cells with any of the mutated genes. The moderate inhibitory effect was conferred by induction of pentose phosphate pathway (PPP) when cells were treated with Gefitinib, the EGFR inhibitor. Combination of Gefitinib with PPP inhibitor 6AN effected synergistically in triple-WT cells. The inhibition could be restored by addition of NADPH. Dual EGFR/ERBB2 inhibitor Afatinib also exhibited similar effects. Further genetic silencing of EGFR, ERBB2 and mTOR indicated that major effect conferred by ERBB2 was via convergence to EGFR pathway in MM. Our results contributed to the individualized targeted therapy with EGFR inhibitors in MM.

Liver cancer is associated with genetic mutations caused by environmental exposures, including occupational exposure to alpha radiation emitted by plutonium. We used whole exome sequencing (WES) to characterize somatic mutations in 3 histologically distinct primary liver tumors (angiosarcoma of the liver (ASL), cholangiocarcinoma (CCA) and hepatocellular carcinoma (HCC)) from Mayak worker subjects occupationally exposed to ionizing radiation (IR) to investigate the contribution of IR to the mutational landscape of liver cancer. DNA sequence analysis revealed these tumors harbor an excess of deletions, with a deletions:substitutions ratio similar to that previously reported in radiation-associated tumors. These tumors were also enriched for clustered mutations, a signature of radiation exposure. Multiple tumors displayed similarities in abrogated gene pathways including actin cytoskeletal signaling and DNA double-strand break (DSB) repair. WES identified novel candidate driver genes in ASL involved in angiogenesis and PIK3CA/AKT/mTOR signaling. We confirmed known driver genes of CCA, and identified candidate driver genes involved in chromatin remodeling. In HCC tumors we validated known driver genes, and identified novel putative driver genes involved in Wnt/β-catenin signaling, chromatin remodeling, PIK3CA/AKT/mTOR signaling, and angiogenesis. This pilot study identifies several novel candidate driver mutations that are likely to be caused by IR exposure, and provides the first data on the mutational landscape of liver cancer after IR exposure.

Canonical ultraviolet (UV) mutation type and spectra are traditionally defined by direct sequencing-based approaches to map mutations in a limited number of representative DNA elements. To obtain an unbiased view of genome wide UV mutation features, we performed whole exome-sequencing (WES) to profile single nucleotide substitutions in UVB-irradiated primary human keratinocytes. Cross comparison of UV mutation profiles under different UVB radiation conditions revealed that T > C transition was highly prevalent in addition to C > T transition. We also identified 5'-ACG-3' as a common sequence motif of C > T transition. Furthermore, our analyses uncovered several recurring UV mutations following acute UVB radiation affecting multiple genes including HRNR, TRIOBP, KCNJ12, and KMT2C, which are frequently mutated in skin cancers, indicating their potential role as founding mutations in UV-induced skin tumorigenesis. Pretreatment with trichostatin A, a pan-histone deacetylase inhibitor that renders chromatin decondensation, significantly decreased the number of mutations in UVB-irradiated keratinocytes. Unexpectedly, we found trichostatin A to be a mutagen that caused DNA damage and mutagenesis at least partly through increased reactive oxidation. In summary, our study reveals new UV mutation features following acute UVB radiation and identifies novel UV mutation hotspots that may potentially represent founding driver mutations in skin cancer development.

Multiple primary malignant neoplasms are rare entities in the clinical setting, but represent an important issue in the clinical management of patients since they could be expression of a genetic predisposition to malignancy. A high resolution genome wide array CGH led us to identify the first case of a de novo constitutional deletion confined to the FBXW7 gene, a well known tumor suppressor, in a patient with a syndromic phenotype characterized by focal segmental glomerulosclerosis and multiple primary early/atypical onset tumors, including Hodgkin's lymphoma, Wilms tumor and breast cancer. Other genetic defects may be associated with patient's phenotype. In this light, constitutional mutations at BRCA1, BRCA2, TP53, PALB2 and WT1 genes were excluded by performing sequencing and MLPA analysis; similarly, we ruled out constitutional abnormalities at the imprinted 11p15 region by methylation specific -MLPA assay. Our observations sustain the role of FBXW7 as cancer predisposition gene and expand the spectrum of its possible associated diseases.

The most common histological subtypes of cutaneous melanoma include superficial spreading and nodular melanoma. However, the spectrum of somatic mutations developed in those lesions and all potential druggable targets have not yet been fully elucidated. We present the results of a sequence capture NGS analysis of 114 primary nodular and superficial spreading melanomas identifying driver mutations using biostatistical, immunohistochemical and/or functional approach. The spectrum and frequency of pathogenic or likely pathogenic variants were identified across 54 evaluated genes, including 59 novel mutations, and the newly identified TP53 loss-of-function mutations p.(L194P) and p.(R280K). Frequently mutated genes most commonly affected the MAPK pathway, followed by chromatin remodeling, and cell cycle regulation. Frequent aberrations were also detected in the genes coding for proteins involved in DNA repair and the regulation and modification of cellular tight junctions. Furthermore, relatively frequent mutations were described in KDR and MET, which represent potential clinically important targets. Those results suggest that with the development of new therapeutic possibilities, not only BRAF testing, but complex molecular testing of cutaneous melanoma may become an integral part of the decision process concerning the treatment of patients with melanoma.

Multiple myeloma (MM) is a B-cell neoplasm that is characterized by the accumulation of malignant plasma cells in the bone marrow. The transcription factor PRDM1 is a master regulator of plasma cell development and is considered to be an oncosuppressor in several lymphoid neoplasms. The PRDM1β isoform is an alternative promoter of the PRDM1 gene that may interfere with the normal role of the PRDM1α isoform. To explain the induction of the PRDM1β isoform in MM and to offer potential therapeutic strategies to modulate its expression, we characterized the cis regulatory elements and epigenetic status of its promoter. We observed unexpected patterns of hypermethylation and hypomethylation at the PRDM1α and PRDM1β promoters, respectively, and prominent H3K4me1 and H3K9me2 enrichment at the PRDM1β promoter in non-expressing cell lines compared to PRDM1β-expressing cell lines. After treatment with drugs that inhibit DNA methylation, we were able to modify the activity of the PRDM1β promoter but not that of the PRDM1α promoter. Epigenetic drugs may offer the ability to control the expression of the PRDM1α/PRDM1β promoters as components of novel therapeutic approaches.

Cancer is a disease of complex genetic alterations, and comprehensive genetic diagnosis is beneficial to match each patient to appropriate therapy. However, acquisition of representative tumor samples is invasive and sometimes impossible. Circulating tumor DNA (ctDNA) is a promising tool to use as a non-invasive biomarker for cancer mutation profiling. Here we implemented targeted next generation sequencing (NGS) with a customized gene panel of 382 cancer-relevant genes on 605 ctDNA samples in multiple cancer types. Overall, tumor-specific mutations were identified in 87% of ctDNA samples, with mutation spectra highly concordant with their matched tumor tissues. 71% of patients had at least one clinically-actionable mutation, 76% of which have suggested drugs approved or in clinical trials. In particular, our study reveals a unique mutation spectrum in Chinese lung cancer patients which could be used to guide treatment decisions and monitor drug-resistant mutations. Taken together, our study demonstrated the feasibility of clinically-useful targeted NGS-based ctDNA mutation profiling to guide treatment decisions in cancer.

The gastrointestinal tract's most commonly occurring primary mesenchymal tumor is the gastrointestinal stromal tumor (GIST). However, few cases worldwide were reported associated with renal cell carcinoma (RCC). Therefore, we aimed to identify the association of genitourinary tumors in patients with GIST in our tertiary care hospital in Saudi Arabia and compare it to the literature. We identified all patients in the pathology department database with the diagnosis of GIST. We excluded duplicate and recurrent cases. We examined patients' files for the presence of RCC, adrenal tumors, or other genitourinary cancer. A systematic review of the association was conducted. From 2003 to 2020, 170 patients had a histopathologic diagnosis of primary GIST, 100 men and 70 women, median age of 57 (range 9-91) years at the time of diagnosis. The site of primary GIST was gastric 103, small bowel 43, mesenteric 5, omentum/peritoneum 7, abdomen 4, isolated adrenal 1, and other 7. Six patients had associated primary genitourinary cancer. Three patients had RCC (two clear cell RCC and one radiologic diagnosis only), and three had adrenal tumors (one adrenal carcinoma, one an isolated adrenal GIST, and one pheochromocytoma). In addition, two patients had a tumor invading the urinary bladder. Although the cohort included 63 men aged 60 or above (median 71 ± 8.7 years, range 60-94), none demonstrated clinical prostatic carcinoma. Data was compared to 69 systematic review articles. We report the rare association between GIST tumors and primary genitourinary cancer, mainly RCC and adrenal tumors. Also, we identified a secondary invasion of the urinary bladder. Unlike the reported series, none of the older male patients had clinical prostate cancer.

Thyroid hormones (THs) are critical regulators of cellular processes, while changes in their levels impact all the hallmarks of cancer. Disturbed expression of type 3 deiodinase (DIO3), the main TH-inactivating enzyme, occurs in several human neoplasms and has been associated with adverse outcomes. Here, we investigated the patterns of DIO3 expression and its prognostic significance in breast cancer. DIO3 expression was evaluated by immunohistochemistry in a primary cohort of patients with breast cancer and validated in a second cohort using RNA sequencing data from the TCGA database. DNA methylation data were obtained from the same database. DIO3 expression was present in normal and tumoral breast tissue. Low levels of DIO3 expression were associated with increased mortality in the primary cohort. Accordingly, low DIO3 mRNA levels were associated with an increased risk of death in a multivariate model in the validation cohort. DNA methylation analysis revealed that the DIO3 gene promoter is hypermethylated in tumors when compared to normal tissue. In conclusion, DIO3 is expressed in normal and tumoral breast tissue, while decreased expression relates to poor overall survival in breast cancer patients. Finally, loss of DIO3 expression is associated with hypermethylation of the gene promoter and might have therapeutic implications.

Rectal cancer (RC) appears to behave differently compared with colon cancer. We aimed to analyze existence of different subtypes of RC depending on distinct features (age of onset and the presence of synchronous primary malignant neoplasms). We compared the clinicopathological, familial and molecular features of three different populations diagnosed with RC (early-onset RC [EORC], late-onset RC, and synchronous RC [SRC]). Eighty-five RCs were identified and were evaluated according to their microsatellite instability, CpG Island Methylator Phenotype (CIMP) and chromosomal instability, as assessed by Next Generation Sequencing and microarray-based comparative genomic hybridization approaches. The results were subjected to cluster analysis. SRCs displayed the most specific characteristics including a trend for the development of multiple malignant neoplasms, a greater proportion of CIMP-High tumors (75%) and more frequent genomic alterations. These findings were confirmed by a clustering analysis that stratified RCs according to their genomic alterations. We also found that EORCs exhibited their own features including an important familial cancer component and a remarkable rate of mutations in TP53 (53%). Together, heterogeneity in RC characteristics by age of disease-onset and SRC warrants further study to optimize tailored prevention, detection and intervention strategies-particularly among young adults.

Circulating cell-free DNA (cfDNA) has the potential to be a specific biomarker for the therapeutic management of lung cancer patients. Here, a new sequencing error-reduction method based on molecular amplification pools (MAPs) was utilized to analyze cfDNA in lung cancer patients. We determined the accuracy of MAPs plasma sequencing with respect to droplet digital polymerase chain reaction assays (ddPCR), and tested whether actionable mutation discovery is improved by next-generation sequencing (NGS) in a clinical setting. This study reports data from 356 lung cancer patients receiving plasma testing as part of routine clinical management. Sequencing of cfDNA via MAPs had a sensitivity of 98.5% and specificity 98.9%. The ddPCR assay was used as the reference, since it is an established, accurate assay that can be performed contemporaneously on the same plasma sample. MAPs sequencing detected somatic variants in 261 of 356 samples (73%). Non-actionable clonal hematopoiesis-associated variants were identified via sequencing in 21% of samples. The accuracy of this cfDNA sequencing approach was similar to that of ddPCR assays in a clinical setting, down to an allele frequency of 0.1%. Due to broader coverage and high sensitivity for insertions and deletions, sequencing via MAPs afforded important detection of additional actionable mutations.

Mutational signatures' association with replication timing (RT) has been studied in cancer samples, but the RT distribution of somatic mutations in non-cancerous cells was only minimally explored. Here, we performed comprehensive analyses of mutational signatures in 2.9 million somatic mutations across multiple non-cancerous tissues, stratified by early and late RT regions. We found that many mutational processes are active mainly or solely in early RT, such as SBS16 in hepatocytes and SBS88 in the colon, or in late RT, such as SBS4 in lung and hepatocytes, and SBS18 across many tissues. The two ubiquitous signatures, SBS1 and SBS5, showed late and early bias, respectively, across multiple tissues and in mutations representing germ cells. We also performed a direct comparison with cancer samples in 4 matched tissue-cancer types. Unexpectedly, while for most signatures the RT bias was consistent in normal tissue and in cancer, we found that SBS1's late RT bias is lost in cancer.

Identification of the associations between microRNA molecules and human diseases from large-scale heterogeneous biological data is an important step for understanding the pathogenesis of diseases in microRNA level. However, experimental verification of microRNA-disease associations is expensive and time-consuming. To overcome the drawbacks of conventional experimental methods, we presented a combinatorial prioritization algorithm to predict the microRNA-disease associations. Importantly, our method can be used to predict microRNAs (diseases) associated with the diseases (microRNAs) without the known associated microRNAs (diseases). The predictive performance of our proposed approach was evaluated and verified by the internal cross-validations and external independent validations based on standard association datasets. The results demonstrate that our proposed method achieves the impressive performance for predicting the microRNA-disease association with the Area Under receiver operation characteristic Curve (AUC), 86.93%, which is indeed outperform the previous prediction methods. Particularly, we observed that the ensemble-based method by integrating the predictions of multiple algorithms can give more reliable and robust prediction than the single algorithm, with the AUC score improved to 92.26%. We applied our combinatorial prioritization algorithm to lung neoplasms and breast neoplasms, and revealed their top 30 microRNA candidates, which are in consistent with the published literatures and databases.

The FGF receptor signaling pathway is recurrently involved in the leukemogenic processes. Oncogenic fusions of FGFR1 with various fusion partners were described in myeloid proliferative neoplasms, and overexpression and mutations of FGFR3 are common in multiple myeloma. In addition, fibroblast growth factors are abundant in the bone marrow, and they were shown to enhance the survival of acute myeloid leukemia cells. Here we investigate the effect of FGFR stimulation on pediatric BCP-ALL cells in vitro, and search for mutations with deep targeted next-generation sequencing of mutational hotspots in FGFR1, FGFR2, and FGFR3. In 481 primary BCP-ALL cases, 28 samples from 19 unique relapsed BCP-ALL cases, and twelve BCP-ALL cell lines we found that mutations are rare (4/481 = 0.8%, 0/28 and 0/12) and do not affect codons which are frequently mutated in other malignancies. However, recombinant ligand FGF2 reduced the response to prednisolone in several BCP-ALL cell lines in vitro. We therefore conclude that FGFR signaling can contribute to prednisolone resistance in BCP-ALL cells, but that activating mutations in this receptor tyrosine kinase family are very rare.

An emergent area of cancer genomics is the identification of driver genes. Driver genes confer a selective growth advantage to the cell. While several driver genes have been discovered, many remain undiscovered, especially those mutated at a low frequency across samples. This study defines new features and builds a pan-cancer model, cTaG, to identify new driver genes. The features capture the functional impact of the mutations as well as their recurrence across samples, which helps build a model unbiased to genes with low frequency. The model classifies genes into the functional categories of driver genes, tumour suppressor genes (TSGs) and oncogenes (OGs), having distinct mutation type profiles. We overcome overfitting and show that certain mutation types, such as nonsense mutations, are more important for classification. Further, cTaG was employed to identify tissue-specific driver genes. Some known cancer driver genes predicted by cTaG as TSGs with high probability are ARID1A, TP53, and RB1. In addition to these known genes, potential driver genes predicted are CD36, ZNF750 and ARHGAP35 as TSGs and TAB3 as an oncogene. Overall, our approach surmounts the issue of low recall and bias towards genes with high mutation rates and predicts potential new driver genes for further experimental screening. cTaG is available at https://github.com/RamanLab/cTaG .

We describe a novel ERBB1/EGFR somatic mutation (p. C329R; c.985 T > C) identified in a patient with JAK2V617F Polycythaemia Vera (PV). This substitution affects a conserved cysteine residue in EGFR domain 2 and leads to the formation of a ligand-independent covalent receptor dimer, associated with increased transforming potential. Aberrant signalling from the EGFRC329R receptor is cell type-dependent and in the TF1.8 erythroid cell line expression of this mutant suppresses EPO-induced differentiation. Clonal analysis shows that the dominant JAK2V617F-positive clone in this PV patient harbors EGFRC329R, thus this mutation may contribute to clonal expansion. Somatic mutations affecting other ERBB and related receptor tyrosine kinases are observed in myeloproliferative neoplasms (MPN), and we show elevated EGFR levels in MPN samples, consistent with previous reports. Thus activation of this group of receptors, via multiple mechanisms, may contribute to clonal growth and survival of the JAK2V617F disease clone in MPN.