Whole-exome sequencing to identify novel somatic mutations in squamous cell lung cancers

Uncovering novel mutational signatures by de novo extraction with SigProfilerExtractor

Mutational signature analysis is commonly performed in cancer genomic studies. Here, we present SigProfilerExtractor, an automated tool for de novo extraction of mutational signatures, and benchmark it against another 13 bioinformatics tools by using 34 scenarios encompassing 2,500 simulated signatures found in 60,000 synthetic genomes and 20,000 synthetic exomes. For simulations with 5% noise, reflecting high-quality datasets, SigProfilerExtractor outperforms other approaches by elucidating between 20% and 50% more true-positive signatures while yielding 5-fold less false-positive signatures. Applying SigProfilerExtractor to 4,643 whole-genome- and 19,184 whole-exome-sequenced cancers reveals four novel signatures. Two of the signatures are confirmed in independent cohorts, and one of these signatures is associated with tobacco smoking. In summary, this report provides a reference tool for analysis of mutational signatures, a comprehensive benchmarking of bioinformatics tools for extracting signatures, and several novel mutational signatures, including one putatively attributed to direct tobacco smoking mutagenesis in bladder tissues.

Whole Exome Sequencing of Biliary Tubulopapillary Neoplasms Reveals Common Mutations in Chromatin Remodeling Genes

Simple Summary

Intraductal tubulopapillary neoplasms (ITPN) have recently been described as rare precursor lesions of pancreatic ductal adenocarcinoma and cholangiocarcinoma. Despite a high number of associated invasive adenocarcinomas at the time of diagnosis, patients with ITPN tend to have a much better clinical outcome than those with classical pancreato-biliary adenocarcinoma. Furthermore, rare molecular studies of ITPN show an unexpected lack of hotspot mutations in common driver genes of pancreato-biliary adenocarcinoma, including KRAS. This article reports the first large comprehensive and comparative molecular study of pancreato-biliary ITPN. In the absence of KRAS mutations, we found a high genetic heterogeneity with enrichment in core signaling pathways, including putative actionable genomic targets in one-third of the cases. Whereas, pancreatic ITPN demonstrates a highly distinct genetic profile, differing from classical pancreatic carcinogenesis, biliary ITPN and classical cholangiocarcinoma share common alterations in key genes of the chromatin remodeling pathway, and therefore, appear more closely related than pancreatic ITPN and classical pancreatic ductal adenocarcinoma PDAC.

Abstract

The molecular carcinogenesis of intraductal tubulopapillary neoplasms (ITPN), recently described as rare neoplasms in the pancreato-biliary tract with a favorable prognosis despite a high incidence of associated pancreato-biliary adenocarcinoma, is still poorly understood. To identify driver genes, chromosomal gains and losses, mutational signatures, key signaling pathways, and potential therapeutic targets, the molecular profile of 11 biliary and 6 pancreatic ITPNs, associated with invasive adenocarcinoma in 14/17 cases, are studied by whole exome sequencing (WES). The WES of 17 ITPNs reveals common copy number variants (CNVs) broadly distributed across the genome, with recurrent chromosomal deletions primarily in 1p36 and 9p21 affecting the tumor suppressors CHD5 and CDKN2A, respectively, and gains in 1q affecting the prominent oncogene AKT3. The identified somatic nucleotide variants (SNVs) involve few core signaling pathways despite high genetic heterogeneity with diverse mutational spectra: Chromatin remodeling, the cell cycle, and DNA damage/repair. An OncoKB search identifies putative actionable genomic targets in 35% of the cases (6/17), including recurrent missense mutations of the FGFR2 gene in biliary ITPNs (2/11, 18%). Our results show that somatic SNV in classical cancer genes, typically associated with pancreato-biliary carcinogenesis, were absent (KRAS, IDH1/2, GNAS, and others) to rare (TP53 and SMAD4, 6%, respectively) in ITPNs. Mutational signature pattern analysis reveals a predominance of an age-related pattern. Our findings highlight that biliary ITPN and classical cholangiocarcinoma display commonalities, in particular mutations in genes of the chromatin remodeling pathway, and appear, therefore, more closely related than pancreatic ITPN and classical pancreatic ductal adenocarcinoma.

A Density Peak-Based Method to Detect Copy Number Variations From Next-Generation Sequencing Data

Copy number variation (CNV) is a common type of structural variations in human genome and confers biological meanings to human complex diseases. Detection of CNVs is an important step for a systematic analysis of CNVs in medical research of complex diseases. The recent development of next-generation sequencing (NGS) platforms provides unprecedented opportunities for the detection of CNVs at a base-level resolution. However, due to the intrinsic characteristics behind NGS data, accurate detection of CNVs is still a challenging task. In this article, we propose a new density peak-based method, called dpCNV, for the detection of CNVs from NGS data. The algorithm of dpCNV is designed based on density peak clustering algorithm. It extracts two features, i.e., local density and minimum distance, from sequencing read depth (RD) profile and generates a two-dimensional data. Based on the generated data, a two-dimensional null distribution is constructed to test the significance of each genome bin and then the significant genome bins are declared as CNVs. We test the performance of the dpCNV method on a number of simulated datasets and make comparison with several existing methods. The experimental results demonstrate that our proposed method outperforms others in terms of sensitivity and F1-score. We further apply it to a set of real sequencing samples and the results demonstrate the validity of dpCNV. Therefore, we expect that dpCNV can be used as a supplementary to existing methods and may become a routine tool in the field of genome mutation analysis.

The Emerging Role of Adhesion GPCRs in Cancer

Aberrant expression, function, and mutation of G protein-coupled receptors (GPCRs) and their signaling partners, G proteins, have been well documented in many forms of cancer. These cell surface receptors and their endogenous ligands are implicated in all aspects of cancer including proliferation, angiogenesis, invasion, and metastasis. Adhesion GPCRs (aGPCRs) form the second largest family of GPCRs, most of which are orphan receptors with unknown physiological functions. This is mainly due to our limited insight into their structure, natural ligands, signaling pathways, and tissue expression profiles. Nevertheless, recent studies show that aGPCRs play important roles in cell adhesion to the extracellular matrix and cell-cell communication, processes that are dysregulated in cancer. Emerging evidence suggests that aGPCRs are implicated in migration, proliferation, and survival of tumor cells. We here review the role of aGPCRs in the five most common types of cancer (lung, breast, colorectal, prostate, and gastric) and emphasize the importance of further translational studies in this field.

A genome-wide association study of cocaine use disorder accounting for phenotypic heterogeneity and gene–environment interaction

BACKGROUND

Phenotypic heterogeneity and complicated gene–environment interplay in etiology are among the primary factors that hinder the identification of genetic variants associated with cocaine use disorder.

METHODS

To detect novel genetic variants associated with cocaine use disorder, we derived disease traits with reduced phenotypic heterogeneity using cluster analysis of a study sample (n = 9965). We then used these traits in genome-wide association tests, performed separately for 2070 African Americans and 1570 European Americans, using a new mixed model that accounted for the moderating effects of 5 childhood environmental factors. We used an independent sample (918 African Americans, 1382 European Americans) for replication.

RESULTS

The cluster analysis yielded 5 cocaine use disorder subtypes, of which subtypes 4 (n = 3258) and 5 (n = 1916) comprised heavy cocaine users, had high heritability estimates (h2 = 0.66 and 0.64, respectively) and were used in association tests. Seven of the 13 identified genetic loci in the discovery phase were available in the replication sample. In African Americans, rs114492924 (discovery p = 1.23 × E−8), a single nucleotide polymorphism in LINC01411, was replicated in the replication sample (p = 3.63 × E−3). In a meta-analysis that combined the discovery and replication results, 3 loci in African Americans were significant genome-wide: rs10188036 in TRAK2 (p = 2.95 × E−8), del-1:15511771 in TMEM51 (p = 9.11 × E−10) and rs149843442 near LPHN2 (p = 3.50 × E−8).

LIMITATIONS

Lack of data prevented us from replicating 6 of the 13 identified loci.

CONCLUSION

Our results demonstrate the importance of considering phenotypic heterogeneity and gene–environment interplay in detecting genetic variations that contribute to cocaine use disorder, because new genetic loci have been identified using our novel analytic method.

A curated target gene pool assisting disease prediction and patient-specific biomarker selection for lung squamous cell carcinoma

There have been hundreds of genes demonstrated to be associated with lung squamous cell carcinoma (LSCC), presenting various degrees of association with this disease. In the present study, gene vectors were investigated as genetic biomarkers for the diagnosis and personalized treatment of LSCC. A LSCC genetic database (LSCC_GD) was developed through literature-associated data analysis, where 260 LSCC target genes were curated. Subsequently, numerous associations between these genes and LSCC were studied. Following this, a sparse representation-based variable selection (SRVS) was employed for gene selection from two LSCC gene expression datasets, followed by a case/control classification. Results were compared using analysis of variance (ANOVA)-based gene selection approaches. Using SRVS, a gene vector was selected from each dataset, resulting in significantly higher classification accuracy (CR), compared with randomly selected genes (For datasets GSE18842 and GSE1987, CR=100 and 100% and permutation P=5.0×10⁻⁴ and 1.8×10⁻³, respectively). The SRVS method outperformed ANOVA in terms of the classification ratio. The results indicated that, for a given dataset, there may be a gene vector from the 260 curated LSCC genes that possesses significant prediction power. SRVS is effective in identifying the optimum gene subset target for personalized treatment.

Familial Esophageal Squamous Cell Carcinoma with damaging rare/germline mutations in KCNJ12/KCNJ18 and GPRIN2 genes

In Iran, esophageal cancer is the fourth common cancers in women and sixth common cancers in men. Here we evaluated the importance of familial risk factors and the role of genetic predisposition in Esophageal Squamous Cell Carcinoma (ESCC) using Whole-Exome Sequencing (WES). Germline damaging mutations were identified in WES data from 9 probands of 9 unrelated ESCC pedigrees. Mutations were confirmed with Sanger sequencing and evaluated amplification-refractory mutation system-Polymerase Chain Reaction (ARMS-PCR) in 50 non-related ethnically matched samples and in complete genomics database. Sixteen candidate variants were detected in ESCC 9 probands. Four of these 16 variants were rare damaging mutations including novel mutations in KCNJ12/KCNJ18, and GPRIN2 genes. This WES study in Iranian patients with ESCC, provides insight into the identification of novel germline mutations in familial ESCC. Our data suggest an association between specific mutations and increased risk of ESCC.

Functions of bromodomain-containing proteins and their roles in homeostasis and cancer

Bromodomains (BRDs) are evolutionarily conserved protein-protein interaction modules that are found in a wide range of proteins with diverse catalytic and scaffolding functions and are present in most tissues. BRDs selectively recognize and bind to acetylated Lys residues - particularly in histones - and thereby have important roles in the regulation of gene expression. BRD-containing proteins are frequently dysregulated in cancer, they participate in gene fusions that generate diverse, frequently oncogenic proteins, and many cancer-causing mutations have been mapped to the BRDs themselves. Importantly, BRDs can be targeted by small-molecule inhibitors, which has stimulated many translational research projects that seek to attenuate the aberrant functions of BRD-containing proteins in disease.

Identification of informative genes and pathways using an improved penalized support vector machine with a weighting scheme

Incorporation of pathway knowledge into microarray analysis has brought better biological interpretation of the analysis outcome. However, most pathway data are manually curated without specific biological context. Non-informative genes could be included when the pathway data is used for analysis of context specific data like cancer microarray data. Therefore, efficient identification of informative genes is inevitable. Embedded methods like penalized classifiers have been used for microarray analysis due to their embedded gene selection. This paper proposes an improved penalized support vector machine with absolute t-test weighting scheme to identify informative genes and pathways. Experiments are done on four microarray data sets. The results are compared with previous methods using 10-fold cross validation in terms of accuracy, sensitivity, specificity and F-score. Our method shows consistent improvement over the previous methods and biological validation has been done to elucidate the relation of the selected genes and pathway with the phenotype under study.

G protein-coupled receptors as promising cancer targets

G protein-coupled receptors (GPCRs) regulate an array of fundamental biological processes, such as growth, metabolism and homeostasis. Specifically, GPCRs are involved in cancer initiation and progression. However, compared with the involvement of the epidermal growth factor receptor in cancer, that of GPCRs have been largely ignored. Recent findings have implicated many GPCRs in tumorigenesis, tumor progression, invasion and metastasis. Moreover, GPCRs contribute to the establishment and maintenance of a microenvironment which is permissive for tumor formation and growth, including effects upon surrounding blood vessels, signaling molecules and the extracellular matrix. Thus, GPCRs are considered to be among the most useful drug targets against many solid cancers. Development of selective ligands targeting GPCRs may provide novel and effective treatment strategies against cancer and some anticancer compounds are now in clinical trials. Here, we focus on tumor related GPCRs, such as G protein-coupled receptor 30, the lysophosphatidic acid receptor, angiotensin receptors 1 and 2, the sphingosine 1-phosphate receptors and gastrin releasing peptide receptor. We also summarize their tissue distributions, activation and roles in tumorigenesis and discuss the potential use of GPCR agonists and antagonists in cancer therapy.

Lung Cancer Genomics in the Era of Accelerated Targeted Drug Development

Lung cancer is the leading cause of cancer-related deaths in the United States and the 5-year overall survival outlook for a patient has not improved in several decades. Recently, however, molecular and genomic profiling of the lung tumors has revealed recurring somatic mutations. As a result the therapeutic landscape of lung cancer is undergoing a paradigm shift from a purely histology-based understanding of the disease to subtype distinctions based on tumor genetics, which has launched cancer-specific, mechanism-based targeted therapies with clear benefit to patients. While targeted therapy advancements are being made at an ever increasing rate, a new challenge in the form of drug resistance has also emerged. This review summarizes the current literature for these issues.

Aberrant Epigenetic Modifications of LPHN2 Function as a Potential Cisplatin-Specific Biomarker for Human Gastrointestinal Cancer

Purpose

Epigenetic alterations of specific genes have recently been identified as diagnostic biomarkers for human cancers. However, there are currently no standardized epigenetic biomarkers for drug sensitivity in human gastrointestinal cancer. Therefore, the aim of this study is to identify a novel epigenetic biomarker in gastrointestinal cancer.

Materials and Methods

Using bisulfite sequencing and pyrosequencing analysis, DNA methylation patterns of gastric, colon primary tissues and their cancer cells were analyzed, and histone modifications were analyzed using chromatin immunoprecipitation assay. In addition, cancer cells were exposed to cisplatin and treated with a DNA methyltransferase inhibitor.

Results

We report that in human gastric and colon cancers, latrophilin 2 (LPHN2) is silenced by epigenetic modifications, including CpG island methylation and aberrant histone modifications. We also confirmed that LPHN2 was silenced by DNA hypermethylation in primary gastric and colon tumor tissues compared to their normal counterparts. Interestingly, we found that cancer cells with methylated LPHN2 showed higher sensitivity to cisplatin. Also, 5-aza- 2′-deoxycytidine combined with cisplatin decreased the cytotoxicity of cisplatin in cancer cells with methylated LPHN2. In addition, LPHN2 knockdown in cancer cells with high LPHN2 expression sensitized these cells to the anti-proliferative effects of cisplatin.

Conclusion

In human gastrointestinal cancer, we found that LPHN2 is regulated by epigenetic modifications, and that cancer cells with lower LPHN2 expression show higher sensitivity to cisplatin. Therefore, the methylation status of LPHN2 is a potential novel epigenetic biomarker for cisplatin treatment in human gastric and colon cancers.

Deletion of olfactomedin 2 induces changes in the AMPA receptor complex and impairs visual, olfactory, and motor functions in mice

Olfactomedin 2 (Olfm2) is a secretory glycoprotein belonging to the family of olfactomedin domain-containing proteins. A previous study has shown that a mutation in OLFM2 is associated with primary open angle glaucoma in Japanese patients. In the present study, we generated Olfm2 deficient mice by replacing the Olfm2 gene with the LacZ gene. The loss of Olfm2 resulted in no gross abnormalities. However, Olfm2 null mice showed reduced exploration, locomotion, olfactory sensitivity, abnormal motor coordination, and anxiety related behavior. The pattern of the Olfm2 gene expression was studied in the brain and eye using β-galactosidase staining. In the brain, Olfm2 was mainly expressed in the olfactory bulb, cortex, piriform cortex, olfactory trabeculae, and inferior and superior colliculus. In the eye expression was detected mainly in retinal ganglion cells. In Olfm2 null mice, the amplitude of the first negative wave in the visual evoked potential test was significantly reduced as compared with wild-type littermates. Olfm2, similar to Olfm1, interacted with the GluR2 subunit of the AMPAR complexes and Olfm2 co-segregated with the AMPA receptor subunit GluR2 and other synaptic proteins in the synaptosomal membrane fraction upon biochemical fractionation of the adult mice cortex and retina. Immunoprecipitation from the synaptosomal membrane fraction of the Olfm2 null mouse brain cortex using the GluR2 antibody showed reduced levels of several components of the AMPAR complex in the immunoprecipitates including Olfm1, PSD95 and CNIH2. These results suggest that heterodimers of Olfm1 and Olfm2 interact with AMPAR more efficiently than Olfm2 homodimers and that Olfm2 plays a role in the organization of the AMPA receptor complexes.

Targeting the phosphatidylinositol 3-kinase (PI3K)/AKT/mammalian target of rapamycin (mTOR) pathway: an emerging treatment strategy for squamous cell lung carcinoma

Squamous cell lung carcinoma accounts for approximately 30% of all non-small cell lung cancers (NSCLCs). Despite progress in the understanding of the biology of cancer, cytotoxic chemotherapy remains the standard of care for patients with squamous cell lung carcinoma, but the prognosis is generally poor. The phosphatidylinositol 3-kinase (PI3K)/AKT/mammalian target of rapamycin (mTOR) pathway is one of the most commonly activated signaling pathways in cancer, leading to cell proliferation, survival, and differentiation. It has therefore become a major focus of clinical research. Various alterations in the PI3K/AKT/mTOR pathway have been identified in squamous cell lung carcinoma and a number of agents targeting these alterations are in clinical development for use as single agents and in combination with other targeted and conventional treatments. These include pan-PI3K inhibitors, isoform-specific PI3K inhibitors, AKT inhibitors, mTOR inhibitors, and dual PI3K/mTOR inhibitors. These agents have demonstrated antitumor activity in preclinical models of NSCLC and preliminary clinical evidence is also available for some agents. This review will discuss the role of the PI3K/AKT/mTOR pathway in cancer and how the discovery of genetic alterations in this pathway in patients with squamous cell lung carcinoma can inform the development of targeted therapies for this disease. An overview of ongoing clinical trials investigating PI3K/AKT/mTOR pathway inhibitors in squamous cell lung carcinoma will also be included.

[Research progress of thyroid transcription factor-1 as molecular marker in lung carcinoma]

肺癌是全球癌症相关死亡的主要病因之一。多数的肺癌患者确诊时已处于晚期，中位生存期为1年左右，5年生存率不足16%，预后极差。近年来的研究热点主要集中于涉及致癌过程中新的生物因子的潜在作用。最近研究表明，甲状腺转录因子（thyroid transcription factor-1, TTF-1）是肺癌的一种特异性的谱系生存癌基因。在肿瘤进展过程中，TTF-1的生物学活性和临床功能表现出相反的作用。在此，本文总结了TTF-1在肺癌的发生、发展、诊断及预后方面的作用，以期为TTF-1作为肺癌新型生物标记物的作用提供见解。

Oligogenic germline mutations identified in early non-smokers lung adenocarcinoma patients

OBJECTIVES

A polygenic model is commonly assumed for the predisposition to common cancers. With respect to lung cancer, Genome Wide Association Studies (GWAS) have identified three loci at 15q25, 5p15.33, and 6p21. However, the relative risks associated with alleles at these loci are low; in addition, the data are limited to smokers, and have not been quite reproducible.

MATERIALS AND METHODS

In order to investigate genetic susceptibility we have adopted an entirely novel patient selection strategy. First, we have selected for adenocarcinoma (ADCA) histology only; second, we have selected non-smokers; third we have selected patients who developed ADCA of lung before the age of 60 and who had an older unaffected sib: we have identified 31 such sib-pairs. Among them, we selected two patients with very early age at disease onset (37- and 49-years old), and having a healthy sibling available for genome comparison older than at least 7 years.

RESULTS

On germline DNA samples of four subjects of two such pairs we have carried out whole exome sequencing. Truncating mutations were detected in 8 'cancer genes' in one affected, and in 5 cancer genes in the other affected subject: but none in the two healthy sibs (p=0.0026). Some of these mutant genes (such as BAG6, SPEN and WISP3) are recognized as major cancer players in lung tumors; others have been previously identified in other human cancers (JAK2, TCEB3C, NELFE, TAF1B, EBLN2), in mouse models (GON4L, NOP58, and RBMX) or in genome-wide association studies (KIAA2018, ZNF311).

CONCLUSIONS

This study identifies for the first time in non-smokers with lung adenocarcinoma specific sets of germline mutations that, together, may predispose to this tumor.