Receptor tyrosine kinase fusions act as a significant alternative driver of the serrated pathway in colorectal cancer development.

Serrated polyps are a clinically and molecularly heterogeneous group of lesions that can contribute to the development of colorectal cancers (CRCs). However, the molecular mechanism underlying the development of serrated lesions is still not well understood. Here, we combined multiple approaches to analyze the genetic alterations in 86 colorectal adenomas (including 35 sessile serrated lesions, 15 traditional adenomas, and 36 conventional adenomatous polyps). We also investigated the in vitro and in vivo oncogenic properties of a novel variant of the NCOA4-RET fusion gene. Molecular profiling revealed that sessile serrated lesions and traditional serrated adenomas have distinct clinicopathological and molecular features. Moreover, we identified receptor tyrosine kinase translocations exclusively in sessile serrated lesions (17%), and the observation was validated in a separate cohort of 34 sessile serrated lesions (15%). The kinase fusions as well as the BRAF and KRAS mutations were mutually exclusive to each other. Ectopic expression of a novel variant of the NCOA4-RET fusion gene promoted cell proliferation in vitro and in vivo, and the proliferation was significantly suppressed by RET kinase inhibitors. All of these underscored the importance of mitogen-activated protein kinase (MAPK) pathway activation in the serrated pathway of colorectal tumorigenesis. In addition, we demonstrated that the kinase fusion may occur early in the precursor lesion and subsequent loss of TP53 may drives the transformation to carcinoma during serrated tumorigenesis. In conclusion, we identified kinase fusions as a significant alternative driver of the serrated pathway in colorectal cancer development, and detecting their presence may serve as a biomarker for the diagnosis of sessile serrated lesions. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

WITER: a powerful method for estimation of cancer-driver genes using a weighted iterative regression modelling background mutation counts.

Genomic identification of driver mutations and genes in cancer cells are critical for precision medicine. Due to difficulty in modelling distribution of background mutation counts, existing statistical methods are often underpowered to discriminate cancer-driver genes from passenger genes. Here we propose a novel statistical approach, weighted iterative zero-truncated negative-binomial regression (WITER, http://grass.cgs.hku.hk/limx/witer or KGGSeq,http://grass.cgs.hku.hk/limx/kggseq/), to detect cancer-driver genes showing an excess of somatic mutations. By fitting the distribution of background mutation counts properly, this approach works well even in small or moderate samples. Compared to alternative methods, it detected more significant and cancer-consensus genes in most tested cancers. Applying this approach, we estimated 229 driver genes in 26 different types of cancers. In silico validation confirmed 78% of predicted genes as likely known drivers and many other genes as very likely new drivers for corresponding cancers. The technical advances of WITER enable the detection of driver genes in TCGA datasets as small as 30 subjects and rescue of more genes missed by alternative tools in moderate or small samples.

EBV-encoded miRNAs target ATM-mediated response in nasopharyngeal carcinoma.

Nasopharyngeal carcinoma (NPC) is a highly invasive epithelial malignancy that is prevalent in southern China and Southeast Asia. It is consistently associated with latent Epstein-Barr virus (EBV) infection. In NPC, miR-BARTs, the EBV-encoded miRNAs derived from BamH1-A rightward transcripts, are abundantly expressed and contribute to cancer development by targeting various cellular and viral genes. In this study, we establish a comprehensive transcriptional profile of EBV-encoded miRNAs in a panel of NPC patient-derived xenografts and an EBV-positive NPC cell line by small RNA sequencing. Among the 40 miR-BARTs, predominant expression of 22 miRNAs was consistently detected in these tumors. Among the abundantly expressed EBV-miRNAs, BART5-5p, BART7-3p, BART9-3p, and BART14-3p could negatively regulate the expression of a key DNA double-strand break (DSB) repair gene, ataxia telangiectasia mutated (ATM), by binding to multiple sites on its 3'-UTR. Notably, the expression of these four miR-BARTs represented more than 10% of all EBV-encoded miRNAs in tumor cells, while downregulation of ATM expression was commonly detected in all of our tested sequenced samples. In addition, downregulation of ATM was also observed in primary NPC tissues in both qRT-PCR (16 NP and 45 NPC cases) and immunohistochemical staining (35 NP and 46 NPC cases) analysis. Modulation of ATM expression by BART5-5p, BART7-3p, BART9-3p, and BART14-3p was demonstrated in the transient transfection assays. These findings suggest that EBV uses miRNA machinery as a key mechanism to control the ATM signaling pathway in NPC cells. By suppressing these endogenous miR-BARTs in EBV-positive NPC cells, we further demonstrated the novel function of miR-BARTs in inhibiting Zta-induced lytic reactivation. These findings imply that the four viral miRNAs work co-operatively to modulate ATM activity in response to DNA damage and to maintain viral latency, contributing to the tumorigenesis of NPC. © 2017 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

PIEZO1 functions as a potential oncogene by promoting cell proliferation and migration in gastric carcinogenesis.

Gastric cancer (GC) is one of the most common malignancies in Asian areas. PIEZO1 has been implied to regulate epithelial homeostasis to play an important role in physiological processes. It is also related to tumor initiation and progression. However, the role of PIEZO1 has not yet been explored in GC. The expression of PIEZO1 in GC cell lines and primary samples was determined by qRT-PCR and Western blot. The clinical correlation of PIEZO1 in GC was evaluated by immunohistochemistry on tissue microarray. The oncogenic role of PIEZO1 was demonstrated in gastric tumorigenesis through a series of functional assays, including cell proliferation, cell invasion, and flow cytometry analysis. Drug sensitivity was also assessed by PIEZO1 knockdown experiment. PIEZO1 exhibited an upregulation in most of the GC cell lines and primary samples compared with non-tumorous gastric epithelial tissues. Increase of PIEZO1 was associated with poor disease specific survival. PIEZO1 knockdown led to inhibitory effect by suppressing cell proliferation and invasion and inhibiting xenograft formation. Decreased PIEZO1 enhanced the sensitivity of Cisplatin or 5-FU treatment. Morphological alteration was also observed in siPIEZO1 treated cells. GTP-Rac1 showed accumulated activated form, while total-RhoA was decreased accompanied with PIEZO1 knockdown. In the present study, PIEZO1 is required for cell proliferation, migration and invasion to promote GC progression. PIEZO1 also maintains cellular morphology related to GC cell motility by regulating the activity of Rho GTPase family members. Our data not only suggested a novel prognostic marker, but also provided a useful clinical therapeutic target for GC.

Assessment of programmed cell death ligand-1 expression by 4 diagnostic assays and its clinicopathological correlation in a large cohort of surgical resected non-small cell lung carcinoma.

Immune checkpoint blockade targeting the PD-1/PD-L1 axis has recently demonstrated efficacy and promise in cancer treatment. Appropriate biomarker selection is therefore essential for improving treatment efficacy. However, the establishment of PD-L1 assay in pathology laboratories is complicated by the presence of multiple testing platforms using different scoring systems. Here we assessed the PD-L1 expression in 713 consecutive non-small cell lung carcinomas by four commercially available PD-L1 immunohistochemical assays, namely, 22C3, 28-8, SP142 and SP263. The analytical performances of the four assays and diagnostic performances across clinically relevant cutoffs were evaluated. The prevalence of PD-L1 (22C3) expression was 21% with a ≥50% cutoff and 56% with a ≥1% cutoff. High PD-L1 expression (using a ≥50% cutoff) was significantly associated with male sex (P = 0.001), ever smoking history (P < 0.001), squamous cell carcinoma (P = 0.001), large cell carcinoma (P < 0.001), lymphoepithelioma-like carcinoma (P = 0.006), sarcomatoid carcinoma (P < 0.001), mutant KRAS (P = 0.005) and wild-type EGFR (P = 0.003). Elevated PD-L1 expression was also significantly associated with shorter survival in patients with adenocarcinoma (log-rank P = 0.026) and remained an independent prognostic factor by multivariable analysis. Among the four assays, 22C3, 28-8 and SP263 were highly concordant for tumor cell scoring. With a cutoff of ≥50% (i.e., the threshold for first-line patient selection), inter-rater agreement was high among the three assays with percentage agreement >97%. In conclusion, three PD-L1 assays showed good analytical performance and a high agreement with each other, but not all cases were correctly classified using the same clinical cutoff. Further studies comparing the predictive value of these assays are required to address the interchangeability of these assays for clinical use.

Inheritance-mode specific pathogenicity prioritization (ISPP) for human protein coding genes.

MOTIVATION: Exome sequencing studies have facilitated the detection of causal genetic variants in yet-unsolved Mendelian diseases. However, the identification of disease causal genes among a list of candidates in an exome sequencing study is still not fully settled, and it is often difficult to prioritize candidate genes for follow-up studies. The inheritance mode provides crucial information for understanding Mendelian diseases, but none of the existing gene prioritization tools fully utilize this information. RESULTS: We examined the characteristics of Mendelian disease genes under different inheritance modes. The results suggest that Mendelian disease genes with autosomal dominant (AD) inheritance mode are more haploinsufficiency and de novo mutation sensitive, whereas those autosomal recessive (AR) genes have significantly more non-synonymous variants and regulatory transcript isoforms. In addition, the X-linked (XL) Mendelian disease genes have fewer non-synonymous and synonymous variants. As a result, we derived a new scoring system for prioritizing candidate genes for Mendelian diseases according to the inheritance mode. Our scoring system assigned to each annotated protein-coding gene (N = 18 859) three pathogenic scores according to the inheritance mode (AD, AR and XL). This inheritance mode-specific framework achieved higher accuracy (area under curve = 0.84) in XL mode. CONCLUSION: The inheritance-mode specific pathogenicity prioritization (ISPP) outperformed other well-known methods including Haploinsufficiency, Recessive, Network centrality, Genic Intolerance, Gene Damage Index and Gene Constraint scores. This systematic study suggests that genes manifesting disease inheritance modes tend to have unique characteristics. AVAILABILITY AND IMPLEMENTATION: ISPP is included in KGGSeq v1.0 (http://grass.cgs.hku.hk/limx/kggseq/), and source code is available from (https://github.com/jacobhsu35/ISPP.git). CONTACT: mxli@hku.hkSupplementary information: Supplementary data are available at Bioinformatics online.

Meta-analysis of genome-wide association studies identifies two loci associated with circulating osteoprotegerin levels.

Osteoprotegerin (OPG) is involved in bone homeostasis and tumor cell survival. Circulating OPG levels are also important biomarkers of various clinical traits, such as cancers and atherosclerosis. OPG levels were measured in serum or in plasma. In a meta-analysis of genome-wide association studies in up to 10 336 individuals from European and Asian origin, we discovered that variants >100 kb upstream of the TNFRSF11B gene encoding OPG and another new locus on chromosome 17q11.2 were significantly associated with OPG variation. We also identified a suggestive locus on chromosome 14q21.2 associated with the trait. Moreover, we estimated that over half of the heritability of OPG levels could be explained by all variants examined in our study. Our findings provide further insight into the genetic regulation of circulating OPG levels.

Predicting mendelian disease-causing non-synonymous single nucleotide variants in exome sequencing studies.

Exome sequencing is becoming a standard tool for mapping Mendelian disease-causing (or pathogenic) non-synonymous single nucleotide variants (nsSNVs). Minor allele frequency (MAF) filtering approach and functional prediction methods are commonly used to identify candidate pathogenic mutations in these studies. Combining multiple functional prediction methods may increase accuracy in prediction. Here, we propose to use a logit model to combine multiple prediction methods and compute an unbiased probability of a rare variant being pathogenic. Also, for the first time we assess the predictive power of seven prediction methods (including SIFT, PolyPhen2, CONDEL, and logit) in predicting pathogenic nsSNVs from other rare variants, which reflects the situation after MAF filtering is done in exome-sequencing studies. We found that a logit model combining all or some original prediction methods outperforms other methods examined, but is unable to discriminate between autosomal dominant and autosomal recessive disease mutations. Finally, based on the predictions of the logit model, we estimate that an individual has around 5% of rare nsSNVs that are pathogenic and carries ~22 pathogenic derived alleles at least, which if made homozygous by consanguineous marriages may lead to recessive diseases.

HYST: a hybrid set-based test for genome-wide association studies, with application to protein-protein interaction-based association analysis.

The extended Simes' test (known as GATES) and scaled chi-square test were proposed to combine a set of dependent genome-wide association signals at multiple single-nucleotide polymorphisms (SNPs) for assessing the overall significance of association at the gene or pathway levels. The two tests use different strategies to combine association p values and can outperform each other when the number of and linkage disequilibrium between SNPs vary. In this paper, we introduce a hybrid set-based test (HYST) combining the two tests for genome-wide association studies (GWASs). We describe how HYST can be used to evaluate statistical significance for association at the protein-protein interaction (PPI) level in order to increase power for detecting disease-susceptibility genes of moderate effect size. Computer simulations demonstrated that HYST had a reasonable type 1 error rate and was generally more powerful than its parents and other alternative tests to detect a PPI pair where both genes are associated with the disease of interest. We applied the method to three complex disease GWAS data sets in the public domain; the method detected a number of highly connected significant PPI pairs involving multiple confirmed disease-susceptibility genes not found in the SNP- and gene-based association analyses. These results indicate that HYST can be effectively used to examine a collection of predefined SNP sets based on prior biological knowledge for revealing additional disease-predisposing genes of modest effects in GWASs.

Two-stage genome-wide association study identifies variants in CAMSAP1L1 as susceptibility loci for epilepsy in Chinese.

In the majority of patients, epilepsy is a complex disorder with multiple susceptibility genes interacting with environmental factors. However, we understand little about its genetic risks. Here, we report the first genome-wide association study (GWAS) to identify common susceptibility variants of epilepsy in Chinese. This two-stage GWAS included a total of 1087 patients and 3444 matched controls. In the combined analysis of the two stages, the strongest signals were observed with two highly correlated variants, rs2292096 [G] [P= 1.0 × 10(-8), odds ratio (OR) = 0.63] and rs6660197 [T] (P= 9.9 × 10(-7), OR = 0.69), with the former reaching genome-wide significance, on 1q32.1 in the CAMSAP1L1 gene, which encodes a cytoskeletal protein. We also refined a previously reported association with rs9390754 (P= 1.7 × 10(-5)) on 6q21 in the GRIK2 gene, which encodes a glutamate receptor, and identified several other loci in genes involved in neurotransmission or neuronal networking that warrant further investigation. Our results suggest that common genetic variants may increase the susceptibility to epilepsy in Chinese.

Post-genome wide association studies and functional analyses identify association of MPP7 gene variants with site-specific bone mineral density.

Our previous genome-wide association study (GWAS) in a Hong Kong Southern Chinese population with extreme bone mineral density (BMD) scores revealed suggestive association with MPP7, which ranked second after JAG1 as a candidate gene for BMD. To follow-up this suggestive signal, we replicated the top single-nucleotide polymorphism rs4317882 of MPP7 in three additional independent Asian-descent samples (n= 2684). The association of rs4317882 reached the genome-wide significance in the meta-analysis of all available subjects (P(meta)= 4.58 × 10(-8), n= 4204). Site heterogeneity was observed, with a larger effect on spine than hip BMD. Further functional studies in a zebrafish model revealed that vertebral bone mass was lower in an mpp7 knock-down model compared with the wide-type (P= 9.64 × 10(-4), n= 21). In addition, MPP7 was found to have constitutive expression in human bone-derived cells during osteogenesis. Immunostaining of murine MC3T3-E1 cells revealed that the Mpp7 protein is localized in the plasma membrane and intracytoplasmic compartment of osteoblasts. In an assessment of the function of identified variants, an electrophoretic mobility shift assay demonstrated the binding of transcriptional factor GATA2 to the risk allele 'A' but not the 'G' allele of rs4317882. An mRNA expression study in human peripheral blood mononuclear cells confirmed that the low BMD-related allele 'A' of rs4317882 was associated with lower MPP7 expression (P= 9.07 × 10(-3), n= 135). Our data suggest a genetic and functional association of MPP7 with BMD variation.

GATES: a rapid and powerful gene-based association test using extended Simes procedure.

The gene has been proposed as an attractive unit of analysis for association studies, but a simple yet valid, powerful, and sufficiently fast method of evaluating the statistical significance of all genes in large, genome-wide datasets has been lacking. Here we propose the use of an extended Simes test that integrates functional information and association evidence to combine the p values of the single nucleotide polymorphisms within a gene to obtain an overall p value for the association of the entire gene. Our computer simulations demonstrate that this test is more powerful than the SNP-based test, offers effective control of the type 1 error rate regardless of gene size and linkage-disequilibrium pattern among markers, and does not need permutation or simulation to evaluate empirical significance. Its statistical power in simulated data is at least comparable, and often superior, to that of several alternative gene-based tests. When applied to real genome-wide association study (GWAS) datasets on Crohn disease, the test detected more significant genes than SNP-based tests and alternative gene-based tests. The proposed test, implemented in an open-source package, has the potential to identify additional novel disease-susceptibility genes for complex diseases from large GWAS datasets.

Risk prediction of complex diseases from family history and known susceptibility loci, with applications for cancer screening.

Risk prediction based on genomic profiles has raised a lot of attention recently. However, family history is usually ignored in genetic risk prediction. In this study we proposed a statistical framework for risk prediction given an individual's genotype profile and family history. Genotype information about the relatives can also be incorporated. We allow risk prediction given the current age and follow-up period and consider competing risks of mortality. The framework allows easy extension to any family size and structure. In addition, the predicted risk at any percentile and the risk distribution graphs can be computed analytically. We applied the method to risk prediction for breast and prostate cancers by using known susceptibility loci from genome-wide association studies. For breast cancer, in the population the 10-year risk at age 50 ranged from 1.1% at the 5th percentile to 4.7% at the 95th percentile. If we consider the average 10-year risk at age 50 (2.39%) as the threshold for screening, the screening age ranged from 62 at the 20th percentile to 38 at the 95th percentile (and some never reach the threshold). For women with one affected first-degree relative, the 10-year risks ranged from 2.6% (at the 5th percentile) to 8.1% (at the 95th percentile). For prostate cancer, the corresponding 10-year risks at age 60 varied from 1.8% to 14.9% in the population and from 4.2% to 23.2% in those with an affected first-degree relative. We suggest that for some diseases genetic testing that incorporates family history can stratify people into diverse risk categories and might be useful in targeted prevention and screening.

A comprehensive framework for prioritizing variants in exome sequencing studies of Mendelian diseases.

Exome sequencing strategy is promising for finding novel mutations of human monogenic disorders. However, pinpointing the casual mutation in a small number of samples is still a big challenge. Here, we propose a three-level filtration and prioritization framework to identify the casual mutation(s) in exome sequencing studies. This efficient and comprehensive framework successfully narrowed down whole exome variants to very small numbers of candidate variants in the proof-of-concept examples. The proposed framework, implemented in a user-friendly software package, named KGGSeq (http://statgenpro.psychiatry.hku.hk/kggseq), will play a very useful role in exome sequencing-based discovery of human Mendelian disease genes.

H. pylori-induced NF-κB-PIEZO1-YAP1-CTGF axis drives gastric cancer progression and cancer-associated fibroblast-mediated tumour microenvironment remodelling.

BACKGROUND: Gastric cancer (GC) is one of the most common tumours in East Asia countries and is associated with Helicobacter pylori infection. H. pylori utilizes virulence factors, CagA and VacA, to up-regulate pro-inflammatory cytokines and activate NF-κB signaling. Meanwhile, the PIEZO1 upregulation and cancer-associated fibroblast (CAF) enrichment were found in GC progression. However, the mechanisms of PIEZO1 upregulation and its involvement in GC progression have not been fully elucidated. METHODS: The CAF enrichment and clinical significance were investigated in animal models and primary samples. The expression of NF-κB and PIEZO1 in GC was confirmed by immunohistochemistry staining, and expression correlation was analysed in multiple GC datasets. GSEA and Western blot analysis revealed the YAP1-CTGF axis regulation by PIEZO1. The stimulatory effects of CTGF on CAFs were validated by the co-culture system and animal studies. Patient-derived organoid and peritoneal dissemination models were employed to confirm the role of the PIEZO1-YAP1-CTGF cascade in GC. RESULTS: Both CAF signature and PIEZO1 were positively correlated with H. pylori infection. PIEZO1, a mechanosensor, was confirmed as a direct downstream of NF-κB to promote the transformation from intestinal metaplasia to GC. Mechanistic studies revealed that PIEZO1 transduced the oncogenic signal from NF-κB into YAP1 signaling, a well-documented oncogenic pathway in GC progression. PIEZO1 expression was positively correlated with the YAP1 signature (CTGF, CYR61, and c-Myc, etc.) in primary samples. The secreted CTGF by cancer cells stimulated the CAF infiltration to form a stiffened collagen-enrichment microenvironment, thus activating PIEZO1 to form a positive feedback loop. Both PIEZO1 depletion by shRNA and CTGF inhibition by Procyanidin C1 enhanced the efficacy of 5-FU in suppressing the GC cell peritoneal metastasis. CONCLUSION: This study elucidates a novel driving PIEZO1-YAP1-CTGF force, which opens a novel therapeutic avenue to block the transformation from precancerous lesions to GC. H. pylori-NF-κB activates the PIEZO1-YAP1-CTGF axis to remodel the GC microenvironment by promoting CAF infiltration. Targeting PIEZO1-YAP1-CTGF plus chemotherapy might serve as a potential therapeutic option to block GC progression and peritoneal metastasis.

Heritability of serum osteoprotegerin.

Serum osteoprotegerin (OPG) level is a key biomarker for numerous traits of clinical importance like diabetes, coronary artery disease, blood pressure, lipid profile, and cancers, but its genetic basis remains poorly understood. We estimated the heritability (h(2)) of serum OPG level in 1442 southern Chinese subjects from 306 families. The h(2) for unadjusted OPG was 0.62 for females and 0.17 for males; and for age-adjusted OPG, 0.75 for females and 0.37 for males. Adjustment for lifestyle factors including calcium and phytoestrogen intake, exercise, smoking, and alcohol consumption exerted only a modest effect on the h(2). In conclusion, we confirmed that circulating OPG is a heritable trait and there is a significant difference in heritability between sexes.

A simple bias correction in linear regression for quantitative trait association under two-tail extreme selection.

Selective genotyping can increase power in quantitative trait association. One example of selective genotyping is two-tail extreme selection, but simple linear regression analysis gives a biased genetic effect estimate. Here, we present a simple correction for the bias.

Whole-genome profiling of nasopharyngeal carcinoma reveals viral-host co-operation in inflammatory NF-κB activation and immune escape.

Interplay between EBV infection and acquired genetic alterations during nasopharyngeal carcinoma (NPC) development remains vague. Here we report a comprehensive genomic analysis of 70 NPCs, combining whole-genome sequencing (WGS) of microdissected tumor cells with EBV oncogene expression to reveal multiple aspects of cellular-viral co-operation in tumorigenesis. Genomic aberrations along with EBV-encoded LMP1 expression underpin constitutive NF-κB activation in 90% of NPCs. A similar spectrum of somatic aberrations and viral gene expression undermine innate immunity in 79% of cases and adaptive immunity in 47% of cases; mechanisms by which NPC may evade immune surveillance despite its pro-inflammatory phenotype. Additionally, genomic changes impairing TGFBR2 promote oncogenesis and stabilize EBV infection in tumor cells. Fine-mapping of CDKN2A/CDKN2B deletion breakpoints reveals homozygous MTAP deletions in 32-34% of NPCs that confer marked sensitivity to MAT2A inhibition. Our work concludes that NPC is a homogeneously NF-κB-driven and immune-protected, yet potentially druggable, cancer.

Novel sib pair selection strategy increases power in quantitative association analysis.

Quantitative-trait association studies have been widely used in search for genetic loci for complex traits in recent years. Yet, fiscal constraints still prohibit many on-going research projects from recruiting a large number of individuals for genotyping to reach a desired level of statistical power. Accordingly, in this article, we describe a novel sib pair sampling strategy for genotyping in QTL association studies. With the use of phenotypic scores (and IBD allele-sharing probabilities if available), the genetic effect of a biallelic additive trait locus can be properly modelled within the maximum-likelihood variance components framework proposed by Fulker et al. (Am J Hum Genet 64(1):259-267, 1999) and sib pairs can be rank-ordered by use of informativeness indices. The performance of our method was investigated using simulation. The power of our approach was shown to be higher when compared with other phenotypic selection schemes. An R-script implementing all the selection approaches (including the traditional phenotype-based ones) used in the simulation is available at http://statgen.hku.hk/jshkwan .

The Landscape of Actionable Molecular Alterations in Immunomarker-Defined Large-Cell Carcinoma of the Lung.

INTRODUCTION: Patients with pulmonary large-cell carcinoma (LCC) have poor prognosis and limited treatment options. The identification of clinically actionable molecular alterations helps to guide personalized cancer treatment decisions. PATIENTS AND METHODS: A consecutive cohort of 789 resected NSCLC cases were reviewed. Fifty-nine NSCLC cases lacking morphologic differentiation, accounting for 7.5% of all resected NSCLCs, were identified and further characterized by immunohistochemistry according to the 2015 WHO lung tumor classification. Molecular alterations were investigated by multiple technologies including target capture sequencing, immunohistochemistry, and fluorescence in situ hybridizations. RESULTS: Of 59 NSCLC cases lacking morphologic differentiation, 20 (33.9%) were reclassified as adenocarcinoma (LCC-AD), 14 (23.7%) as squamous cell carcinoma (LCC-SqCC), and 25 (42.4%) as LCC-Null. Approximately 92% of LCC-Null, 95% of LCC-AD, and 86% of LCC-SqCC harbored clinically relevant alterations. Alterations characteristic of adenocarcinoma (EGFR, KRAS, ALK receptor tyrosine kinase [ALK], ROS1, and serine/threonine kinase 11 [STK11]) were detected in the LCC-AD subgroup but not in LCC-SqCC, whereas squamous-lineage alterations (phosphatidylinositol-4,5-biphosphate 3-kinase catalytic subunit alpha [PIK3CA], SRY-box 2 [SOX2], fibroblast growth factor receptor 1 [FGFR1], and AKT1) were detected in the LCC-SqCC subgroup but not in the LCC-AD group. Although some LCC-Null tumors displayed a genetic profile similar to either adenocarcinoma or squamous-cell carcinoma, more than half of the LCC-Null group were completely devoid of recognizable lineage-specific genetic profiles. High programmed death ligand 1 expression and high frequency of cell cycle regulatory gene alterations were found in the LCC-Null group offering alternative options of targeted therapy. CONCLUSIONS: This comprehensive molecular study provided further insight into the genetic architecture of LCC. The presence of clinically actionable alterations in a majority of the tumors allowed personalized treatment to emerge.