Search

Understanding the mechanisms that lead to lung adenocarcinoma metastasis is important for identifying new therapeutics. Here, the authors document the changes in the transcriptome of human lung adenocarcinoma using single-cell sequencing and link cancer cell signatures to immune cell dynamics.

Abstract Background Tumor cell-intrinsic mechanisms and complex interactions with the tumor microenvironment contribute to therapeutic failure via tumor evolution. It may be possible to overcome treatment resistance by developing a personalized approach against relapsing cancers based on a comprehensive analysis of cell type-specific transcriptomic changes over the clinical course of the disease using single-cell RNA sequencing (scRNA-seq). Methods Here, we used scRNA-seq to depict the tumor landscape of a single case of chemo-resistant metastatic, muscle-invasive urothelial bladder cancer (MIUBC) addicted to an activating Harvey rat sarcoma viral oncogene homolog (HRAS) mutation. In order to analyze tumor evolution and microenvironmental changes upon treatment, we also applied scRNA-seq to the corresponding patient-derived xenograft (PDX) before and after treatment with tipifarnib, a HRAS-targeting agent under clinical evaluation. Results In the parallel analysis of the human MIUBC and the PDX, diverse stromal and immune cell populations recapitulated the cellular composition in the human and mouse tumor microenvironment. Treatment with tipifarnib showed dramatic anticancer effects but was unable to achieve a complete response. Importantly, the comparative scRNA-seq analysis between pre- and post-tipifarnib-treated PDX revealed the nature of tipifarnib-refractory tumor cells and the tumor-supporting microenvironment. Based on the upregulation of programmed death-ligand 1 (PD-L1) in surviving tumor cells, and the accumulation of multiple immune-suppressive subsets from post-tipifarnib-treated PDX, a PD-L1 inhibitor, atezolizumab, was clinically applied; this resulted in a favorable response from the patient with acquired resistance to tipifarnib. Conclusion We presented a single case report demonstrating the power of scRNA-seq for visualizing the tumor microenvironment and identifying molecular and cellular therapeutic targets in a treatment-refractory cancer patient.

Adult T-cell Leukemia/Lymphoma (ATLL) is a rare aggressive T-cell malignancy caused by human T-cell leukemia virus type 1 (HTLV-1) infection. However, little is known about the underlying activated molecular pathways at the single cell level. Moreover, the intercellular communications between the tumor microenvironment (TME) and tumor cells in this malignancy are currently unknown. Difficulties in harvesting fresh tissue in a clinical setting have hampered our deeper understanding of this malignancy. Herein, we examined ATLL using archived fresh frozen tissue after biopsy using single-cell RNA sequencing (scRNA-seq) with T-cell receptor (TCR) clonal analysis. Highly clonal tumor cells showed multiple activating pathways, suggesting dynamic evolution of the malignancy. By dissecting diverse cell types comprising the TME, we identified a novel subset of cancer-associated fibroblast, which showed enriched epidermal growth factor receptor (EGFR)-related transcripts including early growth response 1 and 2 (EGR1 and EGR2). Cancer associated fibroblasts (CAFs) of ATLL play an important role for CD4 T-cell proliferation via FGF7-FGF1 and PDGFA-PDGFRA/B signaling, and CAFs, particularly EGR-enriched, are also associated with CD8 and NKT expansion by EGFR. These findings suggest a potential targeted therapeutic pathway to better treat this neoplasm.

Summary: Despite the worldwide effect of the coronavirus disease 2019 (COVID-19) pandemic, the underlying mechanisms of fatal viral pneumonia remain elusive. Here, we show that critical COVID-19 is associated with enhanced eosinophil-mediated inflammation when compared to non-critical cases. In addition, we confirm increased T helper (Th)2-biased adaptive immune responses, accompanying overt complement activation, in the critical group. Moreover, enhanced antibody responses and complement activation are associated with disease pathogenesis as evidenced by formation of immune complexes and membrane attack complexes in airways and vasculature of lung biopsies from six fatal cases, as well as by enhanced hallmark gene set signatures of Fcγ receptor (FcγR) signaling and complement activation in myeloid cells of respiratory specimens from critical COVID-19 patients. These results suggest that SARS-CoV-2 infection may drive specific innate immune responses, including eosinophil-mediated inflammation, and subsequent pulmonary pathogenesis via enhanced Th2-biased immune responses, which might be crucial drivers of critical disease in COVID-19 patients.

Abstract Background Gynecologic malignancy is one of the leading causes of mortality in female adults worldwide. Comprehensive genomic analysis has revealed a list of molecular aberrations that are essential to tumorigenesis, progression, and metastasis of gynecologic tumors. However, targeting such alterations has frequently led to treatment failures due to underlying genomic complexity and simultaneous activation of various tumor cell survival pathway molecules. A compilation of molecular characterization of tumors with pharmacological drug response is the next step toward clinical application of patient-tailored treatment regimens. Results Toward this goal, we establish a library of 139 gynecologic tumors including epithelial ovarian cancers (EOCs), cervical, endometrial tumors, and uterine sarcomas that are genomically and/or pharmacologically annotated and explore dynamic pharmacogenomic associations against 37 molecularly targeted drugs. We discover lineage-specific drug sensitivities based on subcategorization of gynecologic tumors and identify TP53 mutation as a molecular determinant that elicits therapeutic response to poly (ADP-Ribose) polymerase (PARP) inhibitor. We further identify transcriptome expression of inhibitor of DNA biding 2 (ID2) as a potential predictive biomarker for treatment response to olaparib. Conclusions Together, our results demonstrate the potential utility of rapid drug screening combined with genomic profiling for precision treatment of gynecologic cancers.

Abstract Background The immunogenomic changes triggered by concurrent chemoradiation therapy (CCRT), a standard neoadjuvant treatment for locally advanced esophageal squamous cell carcinoma (ESCC), are unknown. We aimed to analyze the early immunogenomic changes in ESCC induced by CCRT and to correlate them with clinical outcomes. Methods We collected biopsy samples from 40 patients with ESCC and the surgical candidates were treated with 5-fluorouracil (5-FU)/Cisplatin and concurrent radiation therapy. Endoscopic biopsy was performed before and after one treatment cycle of 5-FU/Cisplatin and 5 to 18 fractions of radiation. We analyzed immunogenomic changes using paired whole-exome sequencing (n = 29) and paired whole-transcriptome sequencing (WTS, n = 23). Multiplex immunohistochemistry (IHC) was conducted in four representative pair samples. Results Fourteen out of 23 WTS samples (60.8%) showed increased immune scores after CCRT, as calculated by ESTIMATE. The rate of progression-free survival was higher in patients with increased immune scores compared with the remaining patients (83.1% vs. 57.1%, p = 0.25). Tumor mutation burden and neoantigen load were significantly reduced after CCRT (p

Objective The aim of this study was to investigate the efficacy of globus pallidus interna deep brain stimulation (GPi-DBS) for treating dystonia due to the GNAL mutation. Methods We provide the first report of a dystonia patient with a genetically confirmed GNAL mutation in the Korean population and reviewed the literature on patients with the GNAL mutation who underwent GPi-DBS. We compared the effectiveness of DBS in patients with the GNAL mutation compared to that in patients with DYT1 and DYT6 in a previous study. Results Patients with the GNAL mutation and those with DYT1 had higher early responder rates (GNAL, 5/5, 100%; DYT1, 7/7, 100%) than did patients with DYT6 (p = 0.047). The responder rates at late follow-up did not differ statistically among the three groups (p = 0.278). The decrease in the dystonia motor scale score in the GNAL group was 46.9% at early follow-up and 63.4% at late follow-up. Conclusion GPi-DBS would be an effective treatment option for dystonia patients with the GNAL mutation who are resistant to medication or botulinum toxin treatment.

Abstract Background Target enrichment is a critical component of targeted deep next-generation sequencing for the cost-effective and sensitive detection of mutations, which is predominantly performed by either hybrid selection or PCR. Despite the advantages of efficient enrichment, PCR-based methods preclude the identification of PCR duplicates and their subsequent removal. Recently, this limitation was overcome by assigning a unique molecular identifier(UMI) to each template molecule. Currently, several commercial library construction kits based on PCR enrichment are available for UMIs, but there have been no systematic studies to compare their performances. In this study, we evaluated and compared the performances of five commercial library kits from four vendors: the Archer® Reveal ctDNA™ 28 Kit, NEBNext Direct® Cancer HotSpot Panel, Nugen Ovation® Custom Target Enrichment System, Qiagen Human Comprehensive Cancer Panel(HCCP), and Qiagen Human Actionable Solid Tumor Panel(HASTP). Results We evaluated and compared the performances of the five kits using 50 ng of genomic DNA for the library construction in terms of the library complexity, coverage uniformity, and errors in the UMIs. While the duplicate rates for all kits were dramatically decreased by identifying unique molecules with UMIs, the Qiagen HASTP achieved the highest library complexity based on the depth of unique coverage indicating superb library construction efficiency. Regarding the coverage uniformity, the kits from Nugen and NEB performed the best followed by the kits from Qiagen. We also analyzed the UMIs, including errors, which allowed us to adjust the depth of unique coverage and the length required for sufficient complexity. Based on these comparisons, we selected the Qiagen HASTP for further performance evaluations. The targeted deep sequencing method based on PCR target enrichment combined with UMI tagging sensitively detected mutations present at a frequency as low as 1% using 6.25 ng of human genomic DNA as the starting material. Conclusion This study is the first systematic evaluation of commercial library construction kits for PCR-based targeted deep sequencing utilizing UMIs. Because the kits displayed significant variability in different quality metrics, our study offers a practical guideline for researchers to choose appropriate options for PCR-based targeted sequencing and useful benchmark data for evaluating new kits.

Primary refractory/relapsed diffuse large B-cell lymphoma (rrDLBCL) is an unresolved issue for DLBCL treatment and new treatments to overcome resistance is required. To explore the genetic mechanisms underlying treatment resistance in rrDLBCL and to identify candidate genes, we performed targeted deep sequencing of 430 lymphoma-related genes from 58 patients diagnosed with rrDLBCL. Genetic alterations found between the initial biopsy and biopsy at recurrence or refractory disease were investigated. The genes most frequently altered (> 20%) were (in decreasing order of frequency) CDKN2A, PIM1, CD79B, TP53, MYD88, MYC, BTG2, BTG1, CDKN2B, DTX1, CD58, ETV6, and IRF4. Genes mutation of which in pretreatment sample were associated with poor overall survival included NOTCH1, FGFR2, BCL7A, BCL10, SPEN and TP53 (P < 0.05). FGFR2, BCL2, BCL6, BCL10, and TP53 were associated with poor progression-free survival (P < 0.05). Most mutations were truncal and were maintained in both the initial biopsy and post-treatment biopsy with high dynamics of subclones. Immune-evasion genes showed increased overall mutation frequency (CD58, B2M) and variant allele fraction (CD58), and decreased copy number (B2M, CD70) at the post-treatment biopsy. Using the established mutational profiles and integrative analysis of mutational evolution, we identified information about candidate genes that may be useful for the development of future treatment strategies.

Microcephaly is a prevalent phenotype in patients with neurodevelopmental problems, often with genetic causes. We comprehensively investigated the clinical phenotypes and genetic background of microcephaly in 40 Korean patients. We analyzed their clinical phenotypes and radiologic images and conducted whole exome sequencing (WES) and analysis of copy number variation (CNV). Infantile hypotonia and developmental delay were present in all patients. Thirty-four patients (85%) showed primary microcephaly. The diagnostic yield from the WES and CNV analyses was 47.5%. With WES, we detected pathogenic or likely pathogenic variants that were previously associated with microcephaly in 12 patients (30%); nine of these were de novo variants with autosomal dominant inheritance. Two unrelated patients had mutations in the KMT2A gene. In 10 other patients, we found mutations in the GNB1, GNAO1, TCF4, ASXL1, SMC1A, VPS13B, ACTG1, EP300, and KMT2D genes. Seven patients (17.5%) were diagnosed with pathogenic CNVs. Korean patients with microcephaly show a genetic spectrum that is different from that of patients with microcephaly of other ethnicities. WES along with CNV analysis represents an effective approach for diagnosis of the underlying causes of microcephaly.

Background Tumor mutational burden (TMB) measurement is limited by low tumor purity of samples, which can influence prediction of the immunotherapy response, particularly when using whole-exome sequencing-based TMB (wTMB). This issue could be overcome by targeted panel sequencing-based TMB (pTMB) with higher depth of coverage, which remains unexplored.Methods We comprehensively reanalyzed four public datasets of immune checkpoint inhibitor (ICI)-treated cohorts (adopting pTMB or wTMB) to test each biomarker’s predictive ability for low purity samples (cut-off: 30%). For validation, paired genomic profiling with the same tumor specimens was performed to directly compare wTMB and pTMB in patients with breast cancer (paired-BRCA, n=165) and ICI-treated patients with advanced non-small-cell lung cancer (paired-NSCLC, n=156).Results Low tumor purity was common (range 30%–45%) in real-world samples from ICI-treated patients. In the survival analyzes of public cohorts, wTMB could not predict the clinical benefit of immunotherapy when tumor purity was low (log-rank p=0.874), whereas pTMB could effectively stratify the survival outcome (log-rank p=0.020). In the paired-BRCA and paired-NSCLC cohorts, pTMB was less affected by tumor purity, with significantly more somatic variants identified at low allele frequency (p

Summary: Exome and transcriptome analyses of clinically homogeneous early-stage never-smoker female patients with lung adenocarcinoma were performed to understand tumor-T cell interactions and immune escape points. Using our novel gene panels of eight functional categories in the cancer-immunity cycle, three distinct subgroups were identified in this immune checkpoint blockade-refractory cohort by defective gene expression in two major domains, i.e., type I interferon production/signaling pathway and antigen-presenting machinery. Our approach could play a critical role in understanding immune evasion mechanisms, developing a method for effective selection of rare immune checkpoint blockade responders, and finding new treatment strategies.

BACKGROUND: Prostate cancers exhibit intratumor heterogeneity (ITH), like other cancer types. The ITH may affect diverse phenotypes such as treatment response, drug resistance, and clinical outcomes. It is crucial to consider ITH to understand tumorigenesis. METHODS: Genomic and transcriptomic profiles of prostate cancer patients were investigated to determine which markers are correlated with the degree of tumor heterogeneity. In addition, the correlation between the immune activity and clonality of tumors was examined. RESULTS: Tumor heterogeneity across all prostate cancer samples was variable. However, ITH events were dependent on genomic and clinical features. Interestingly, prostate-specific antigen score increased in tumors with multiple subclones, indicating high-grade tumor heterogeneity. On the other hand, CD8-positive T-cell activation decreased in highly heterogeneous tumors. Intriguingly, PTEN deletion was prominently enriched in high heterogeneity groups, with a strong association with heterozygous loss. Expression of major genes including PTEN, CDC42EP5, RNLS, GP2, NETO2, and AMPD3 was closely related to tumor heterogeneity in association with PTEN deletion. CONCLUSIONS: In prostate cancer, ITH, a potential factor affecting tumor progression, is associated with PTEN deletion and cytotoxic T cell inactivation.

PURPOSE: Relapsed/refractory pediatric cancers show poor prognosis; however, their genomic patterns remain unknown. To investigate the genetic mechanisms of tumor relapse and therapy resistance, we characterized genomic alterations in diagnostic and relapsed lesions in patients with relapsed/refractory pediatric solid tumors using targeted deep sequencing. PATIENTS AND METHODS: A targeted sequencing panel covering the exons of 381 cancer genes was used to characterize 19 paired diagnostic and relapsed samples from patients with relapsed/refractory pediatric solid tumors. RESULTS: The mean coverage for all samples was 930.6× (SD = 213.8). Among the 381 genes, 173 single nucleotide variations (SNVs)/insertion-deletions (InDels), 100 copy number alterations, and 1 structural variation were detected. A total of 72.6% of SNVs in primary tumors were also found in recurrent lesions, and 27.2% of SNVs in recurrent tumors had newly occurred. Among SNVs/InDels detected only in recurrent lesions, 71% had a low variant allele fraction (

Abstract Background While auditory neuropathy spectrum disorder (ANSD) is a heterogeneous disorder and its management quite varies depending upon the etiology, even including self-resolution, OTOF is an important molecular etiology of prelingual ANSD and has emerged as an attractive target for implementation of precision medicine in terms of timing and prognosis prediction of auditory rehabilitation. However, to date, the literature is lacking in the genotype–phenotype relationship of this gene as well as efficient molecular testing strategy in the clinic in many populations and to make things more complicated in Koreans, the most prevalent variant p.Arg1939Gln among Korean ANSD children frequently evaded detection by next generation sequencing (NGS), resulting in delayed genetic diagnosis and late cochlear implantation (CI). The aims of this study are to document the mutational and phenotypic spectrum of OTOF-related ANSD (DFNB9) in the Korean population, further establishing genotype–phenotype correlation and proposing a set of the most commonly found OTOF variants to be screened first. Methods Genetic diagnosis through the NGS-based sequencing was made on patients with ANSD in two tertiary hospitals. Genotype and phenotypes of eleven DFNB9 patients were reviewed. For data analysis, Mann–Whitney test and Fisher’s exact test were applied. Results This study disclosed four prevalent variants in Koreans: p.Arg1939Gln with an allele frequency of 40.9%, p.Glu841Lys (13.6%), p.Leu1011Pro and p.Arg1856Trp (9.1%). Three novel variants (c.4227 + 5G > C, p.Gly1845Glu, and p.Pro1931Thr) were identified. Interestingly, a significant association of p.Arg1939Gln with worse ASSR thresholds was observed despite consistently no ABR response. Ten of 11 DFNB9 patients received CI for auditory rehabilitation, showing favorable outcomes with more rapid improvement on early-CI group (age at CI ≤ 18 mo.) than late-CI group. Conclusions This study included the largest Korean DFNB9 cohort to date and proposed a set of the most frequent four OTOF variants, allowing the potential prioritization of exons during Sanger sequencing. Further, a significant association of p.Arg1939Gln homozygotes with poor residual hearing was observed. We may have to suspect p.Arg1939Gln homozygosity in cases of poor auditory thresholds in ANSD children with putative negative OTOF variants solely screened by NGS. Reciprocal feedback between bench and clinics regarding DFNB9 would complement each other.

Abstract Spatial mapping of genomic data to tissue context in a high-throughput and high-resolution manner has been challenging due to technical limitations. Here, we describe PHLI-seq, a novel approach that enables high-throughput isolation and genome-wide sequence analysis of single cells or small numbers of cells to construct genomic maps within cancer tissue in relation to the images or phenotypes of the cells. By applying PHLI-seq, we reveal the heterogeneity of breast cancer tissues at a high resolution and map the genomic landscape of the cells to their corresponding spatial locations and phenotypes in the 3D tumor mass.

BACKGROUND: Signet-ring cell carcinoma (SRCC) is a very rare subtype of colorectal adenocarcinoma (COAD) with a poor clinical prognosis. Although understanding key mechanisms of tumor progression in SRCCs is critical for precise treatment, a comprehensive view of genomic alterations is lacking. MATERIALS AND METHODS: We performed whole-exome sequencing of tumors and matched normal blood as well as RNA sequencing of tumors and matched normal colonic tissues from five patients with SRCC. RESULTS: We identified major somatic alterations and characterized transcriptional changes at the gene and pathway level. Based on high-throughput sequencing, the pattern of mutations and copy number variations was overall similar to that of COAD. Transcriptome analysis revealed that major transcription factors, such as SRF, HNF4A, ZEB1, and RUNX1, with potential regulatory roles in key pathways, including focal adhesion, the PI3K-Akt signaling pathway, and the MAPK signaling pathway, may play a role in the tumorigenesis of SRCC. Furthermore, significantly upregulated genes in SRCCs were enriched for epithelial-mesenchymal transition genes, and accumulation of mucin in intracytoplasm was associated with the overexpression of MUC2. CONCLUSION: The results indicate that the molecular basis of colorectal SRCC exhibits key differences from that of consensus COAD. Our findings clarify important genetic features of particular abnormalities in SRCCs.

While breast cancer incidence in the Asia Pacific region is rising, the molecular basis remains poorly characterized. Here the authors perform genomic screening of 187 Korean breast cancer patients and find differences in molecular subtype distribution, mutation pattern and prevalence, and gene expression signature when compared to TCGA.

Abstract Background MYO15A variants, except those in the N-terminal domain, have been shown to be associated with congenital or pre-lingual severe-to-profound hearing loss (DFNB3), which ultimately requires cochlear implantation in early childhood. Recently, such variants have also been shown to possibly cause moderate-to-severe hearing loss. Herein, we also demonstrate that some MYO15A mutant alleles can cause postlingual onset of progressive partial deafness. Methods Two multiplex Korean families (SB246 and SB224), manifesting postlingual, progressive, partial deafness in an autosomal recessive fashion, were recruited. Molecular genetics testing was performed in two different pipelines, in a parallel fashion, for the SB246 family: targeted exome sequencing (TES) of 129 known deafness genes from the proband and whole exome sequencing (WES) of all affected subjects. Only the former pipeline was performed for the SB224 family. Rigorous bioinformatics analyses encompassing structural variations were executed to investigate any causative variants. Results In the SB246 family, two different molecular diagnostic pipelines provided exactly the same candidate variants: c.5504G > A (p.R1835H) in the motor domain and c.10245_10247delCTC (p.S3417del) in the FERM domain of MYO15A. In the SB224 family, c.9790C > T (p.Q3264X) and c.10263C > G (p.I3421M) in the FERM domain were detected as candidate variants. Conclusions Some recessive MYO15A variants can cause postlingual onset of progressive partial deafness. The phenotypic spectrum of DFNB3 should be extended to include such partial deafness. The mechanism for a milder phenotype could be due to the milder pathogenic potential from hypomorphic alleles of MYO15A or the presence of modifier genes. This merits further investigation.

PURPOSE: Phyllodes tumors (PTs) of the breast are rare, accounting for less than 1% of all breast tumors. Among PTs, malignant PTs (MPTs) have malignant characteristics and distant metastases occur in about 20% to 30% of MPTs. However, there is no effective treatment for MPTs with distant metastasis, resulting in an abject prognosis. We performed targeted deep sequencing on PTs to identify the associations between genetic alterations and clinical prognosis. METHODS: We performed targeted deep sequencing to evaluate the genetic characteristics of PTs and analyzed the relationships between clinical and genetic characteristics. RESULTS: A total of 17 PTs were collected between 2001 and 2012. Histologic review was performed by pathologists. The samples included three benign PTs, one borderline PT, and 13 MPTs. The most frequently detected genetic alteration occurred in the TERT promoter region (70.6%), followed by MED12 (64.7%). EGFR amplification and TP53 alteration were detected in four MPTs without genetic alterations in MED12 and TERT promoter regions. Genetic alterations of RARA and ZNF703 were repeatedly found in PTs with local recurrence, and genetic alterations of SETD2, BRCA2, and TSC1 were detected in PTs with distant metastasis. Especially, MPT harboring PTEN and RB1 copy number deletion showed rapid disease progression. CONCLUSIONS: In this study, we provide genetic characterization and potential therapeutic target for this rare, potentially lethal disease. Further large-scale comprehensive genetic study and functional validation are warranted.

High-throughput sequencing is used to identify somatic variants in cancer patients. Here, the authors perform panel-based profiling of 5095 clinical samples and demonstrate that many clinically-actionable variants have low variant allele fractions, requiring assays with high detection sensitivity.

Purpose Neurotropic tropomyosin receptor kinase (NTRK) fusions have been identified in a variety of cancers, and tyrosine kinase inhibitors targeting the tropomyosin receptor kinase (TRK) receptor are currently in clinical trials. However, no reports are available on the effects of NTRK gene amplification. Methods Samples from patients enrolled in the sequencing program were analyzed using a next-generation sequencing (NGS) cancer panel. For cases in which NTRK amplification (defined as ≥ 4.0 copies) was identified, panTRK immunohistochemical (IHC) staining of tissue microarrays was performed. Results A total of 1,250 tumor specimens collected between February 2014 and January 2016 were analyzed using the NGS cancer panel. NTRK amplification was detected in 28 cases of various types of cancer. Among 27 cases, only four were positive for pan-TRK IHC. These four cases were melanoma, sarcoma, lung cancer, and gastric cancer. We found that 2.2% of cancer patients showed NTRK amplification using NGS cancer panel and NTRK amplification resulted in protein overexpression in 14.8% of these patients. Conclusion Patients with NTRK amplification and increased TRK protein expression may be considered for inclusion in clinical trials for NTRK inhibitors.

Abstract Background Patient-derived xenograft (PDX) models are important tools in precision medicine and for the development of targeted therapies to treat cancer patients. This study aimed to evaluate our precision medicine strategy that integrates genomic profiling and preclinical drug-screening platforms, in order to personalize cancer treatments using PDX models. Methods We performed array-comparative genomic hybridization, microarray, and targeted next-generation sequencing analyses, in order to determine the oncogenic driver mutations. PDX cells were obtained from PDXs and subsequently screened in vitro with 17 targeted agents. Results PDX tumors recapitulated the histopathologic and genetic features of the patient tumors. Among the samples from lung cancer patients that were molecularly-profiled, copy number analysis identified unique focal MET amplification in one sample, 033 T, without RTK/RAS/RAF oncogene mutations. Although HER2 amplification in 033 T was not detected in the cancer panel, the selection of HER2-amplified clones was found in PDXs and PDX cells. Additionally, MET and HER2 overexpression were found in patient tumors, PDXs, and PDX cells. Crizotinib or EGFR tyrosine kinase inhibitor treatments significantly inhibited cell growth and impaired tumor sphere formation in 033 T PDX cells. Conclusions We established PDX cell models using surgical samples from lung cancer patients, and investigated their preclinical and clinical implications for personalized targeted therapy. Additionally, we suggest that MET and EGFR inhibitor-based therapy can be used to treat MET and HER2-overexpressing lung cancers, without receptor tyrosine kinase /RAS/RAF pathway alterations.

Abstract Background Targeted deep sequencing is increasingly used to detect low-allelic fraction variants; it is therefore essential that errors that constitute baseline noise and impose a practical limit on detection are characterized. In the present study, we systematically evaluate the extent to which errors are incurred during specific steps of the capture-based targeted sequencing process. Results We removed most sequencing artifacts by filtering out low-quality bases and then analyze the remaining background noise. By recognizing that plasma DNA is naturally fragmented to be of a size comparable to that of mono-nucleosomal DNA, we were able to identify and characterize errors that are specifically associated with acoustic shearing. Two-thirds of C:G > A:T errors and one quarter of C:G > G:C errors were attributed to the oxidation of guanine during acoustic shearing, and this was further validated by comparative experiments conducted under different shearing conditions. The acoustic shearing step also causes A > G and A > T substitutions localized to the end bases of sheared DNA fragments, indicating a probable association of these errors with DNA breakage. Finally, the hybrid selection step contributes to one-third of the remaining C:G > A:T and one-fifth of the C > T errors. Conclusions The results of this study provide a comprehensive summary of various errors incurred during targeted deep sequencing, and their underlying causes. This information will be invaluable to drive technical improvements in this sequencing method, and may increase the future usage of targeted deep sequencing methods for low-allelic fraction variant detection.

Some cancers with DNA mismatch repair deficiency display microsatellite instability. Here the authors analyse twenty three cancer types at the exome and whole-genome level, and identify loci with recurrent microsatellite instability that could be used to identify patients who would benefit from immunotherapy.

PURPOSE: Although targeting angiogenesis with tyrosine kinase inhibitors (TKIs) has become standard of care in the treatment of clear cell renal cell carcinoma (RCC), resistance mechanism are not fully understood, and there is a need to develop new therapeutic options overcoming them. METHODS AND MATERIALS: To develop a preclinical model that predicts clinical activity of novel agents in 19 RCC patients, we established patient-derived cell (PDC) and xenograft (PDX) models derived from malignant effusions or surgical specimen. RESULTS: Successful PDCs, defined as cells that maintained growth following two passages, were established in 5 of 15 malignant effusions and 1 of 4 surgical specimens. One PDC, clinically refractory to TKIs, was implanted and engrafted in mice, resulting in a comparable histology to the primary tumor. The PDC-PDX model also showed similar genomic features when tested using targeted sequencing of cancer-related genes. When we examined the drug effects of the PDX model, the tumor cells showed resistance to TKIs and everolimus in vitro. CONCLUSION: The results suggest that the PDC-PDX preclinical model we developed using malignant effusions can be a useful preclinical model to interrogate sensitivity to targeted agents based on genomic alterations.

Abstract Despite efforts to interrogate human genome variation through large-scale databases, systematic preference toward populations of Caucasian descendants has resulted in unintended reduction of power in studying non-Caucasians. Here we report a compilation of coding variants from 1,055 healthy Korean individuals (KOVA; Korean Variant Archive). The samples were sequenced to a mean depth of 75x, yielding 101 singleton variants per individual. Population genetics analysis demonstrates that the Korean population is a distinct ethnic group comparable to other discrete ethnic groups in Africa and Europe, providing a rationale for such independent genomic datasets. Indeed, KOVA conferred 22.8% increased variant filtering power in addition to Exome Aggregation Consortium (ExAC) when used on Korean exomes. Functional assessment of nonsynonymous variant supported the presence of purifying selection in Koreans. Analysis of copy number variants detected 5.2 deletions and 10.3 amplifications per individual with an increased fraction of novel variants among smaller and rarer copy number variable segments. We also report a list of germline variants that are associated with increased tumor susceptibility. This catalog can function as a critical addition to the pre-existing variant databases in pursuing genetic studies of Korean individuals.

Abstract Somatic mutations in cancer genomes often show allelic imbalance (AI) of mutation abundance between the genome and transcriptome, but there is not yet a systematic understanding of AI. In this study, we performed large-scale DNA and RNA AI analyses of >100,000 somatic mutations in >2,000 cancer specimens across five tumor types using the exome and transcriptome sequencing data of the Cancer Genome Atlas consortium. First, AI analysis of nonsense mutations and frameshift indels revealed that nonsense-mediated decay is typical in cancer genomes, and we identified the relationship between the extent of AI and the location of mutations in addition to the well-recognized 50-nt rules. Second, the AI with splice site mutations may reflect the extent of intron retention and is frequently observed in known tumor suppressor genes. For missense mutations, we observed that mutations frequently subject to AI are enriched to genes related to cancer, especially those of apoptosis and the extracellular matrix, and C:G > A:T transversions. Our results suggest that mutations in known cancer-related genes and their transcripts are subjected to different levels of transcriptional or posttranscriptional regulation compared to wildtype alleles and may add an additional regulatory layer to the functions of cancer-relevant genes.

Genetic heterogeneity in breast cancer has been demonstrated at a single-cell resolution with high levels of genome coverage. Here, the authors perform transcriptome analysis of 515 single cells from 11 patients and define core gene expression signatures for subtype-specific single breast cancer cells and tumour-infiltrating immune cells.

Next-generation sequencing (NGS) has recently emerged as an essential component of personalized cancer medicine due to its high throughput and low per-base cost. However, no sufficient guidelines for implementing NGS as a clinical molecular pathology test are established in Korea. To ensure clinical grade quality without inhibiting adoption of NGS, a taskforce team assembled by the Korean Society of Pathologists developed laboratory guidelines for NGS cancer panel testing procedures and requirements for clinical implementation of NGS. This consensus standard proposal consists of two parts: laboratory guidelines and requirements for clinical NGS laboratories. The laboratory guidelines part addressed several important issues across multistep NGS cancer panel tests including choice of gene panel and platform, sample handling, nucleic acid management, sample identity tracking, library preparation, sequencing, analysis and reporting. Requirements for clinical NGS tests were summarized in terms of documentation, validation, quality management, and other required written policies. Together with appropriate pathologist training and international laboratory standards, these laboratory standards would help molecular pathology laboratories to successfully implement NGS cancer panel tests in clinic. In this way, the oncology community would be able to help patients to benefit more from personalized cancer medicine.

Glioblastoma multiforme (GBM) is clinically highly aggressive as a result of evolutionary dynamics induced by cross-talk between cancer cells and a heterogeneous group of immune cells in tumor microenvironment. The brain harbors limited numbers of immune cells with few lymphocytes and macrophages; thus, innate-like lymphocytes, such as γδ T cells, have important roles in antitumor immunity. Here, we characterized GBM-infiltrating γδ T cells, which may have roles in regulating the GBM tumor microenvironment and cancer cell gene expression. V(D)J repertoires of tumor-infiltrating and blood-circulating γδ T cells from four patients were analyzed by next-generation sequencing-based T-cell receptor (TCR) sequencing in addition to mutation and immune profiles in four GBM cases. In all tumor tissues, abundant innate and effector/memory lymphocytes were detected, accompanied by large numbers of tumor-associated macrophages and closely located tumor-infiltrating γδ T cells, which appear to have anti-tumor activity. The immune-related gene expression analysis using the TCGA database showed that the signature gene expression extent of γδ T cells were more associated with those of cytotoxic T and Th1 cells and M1 macrophages than those of Th2 cells and M2 macrophages. Although the most abundant γδ T cells were Vγ9Vδ2 T cells in both tumor tissues and blood, the repertoire of intratumoral Vγ9Vδ2 T cells was distinct from that of peripheral blood Vγ9Vδ2 T cells and was dominated by Vγ9Jγ2 sequences, not by canonical Vγ9JγP sequences that are mostly commonly found in blood γδ T cells. Collectively, unique GBM-specific TCR clonotypes were identified by comparing TCR repertoires of peripheral blood and intra-tumoral γδ T cells. These findings will be helpful for the elucidation of tumor-specific antigens and development of anticancer immunotherapies using tumor-infiltrating γδ T cells.

Background: Epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs) are considered the first line treatment for a subset of EGFR-mutated non-small cell lung cancer (NSCLC) patients. Although transformation to small cell lung cancer (SCLC) is one of the known mechanisms of resistance to EGFR TKIs, it is not certain whether transformation to SCLC is exclusively found as a mechanism of TKI resistance in EGFR-mutant tumors. Methods: We identified six patients with primary lung adenocarcinoma that showed transformation to SCLC on second biopsy (n = 401) during a 6-year period. Clinicopathologic information was analyzed and EGFR mutation results were compared between initial and second biopsy samples. Results: Six patients showed transformation from adenocarcinoma to SCLC, of which four were pure SCLCs and two were combined adenocarcinoma and SCLCs. Clinically, four cases were EGFR-mutant tumors from non-smoking females who underwent TKI treatment, and the EGFR mutation was retained in the transformed SCLC tumors. The remaining two adenocarcinomas were EGFR wild-type, and one of these patients received EGFR TKI treatment. Conclusions: NSCLC can acquire a neuroendocrine phenotype with or without EGFR TKI treatment.

Background: Biliary tract cancers (BTCs) are rare and heterogeneous group of tumors classified anatomically into intrahepatic and extrahepatic bile ducts and gallbladder adenocarcinomas. Patient-derived tumor cell (PDC) models with genome analysis can be a valuable platform to develop a method to overcome the clinical barrier on BTCs. Material and Methods: Between January 2012 and June 2015, 40 BTC patients’ samples were collected. PDCs were isolated and cultured from surgical specimens, biopsy tissues, or malignant effusions including ascites and pleural fluid. Genome analysis using targeted panel sequencing as well as digital multiplexed gene analysis was applied to PDCs as well as primary tumors. Results: Extrahepatic cholangiocarcinoma (N = 15, 37.5%), intrahepatic cholangiocarcinoma (N = 10, 25.0%), gallbladder cancer (N = 14, 35.0%), and ampulla of Vater cancer (N = 1, 2.5%) were included. We identified 15 mutations with diverse genetic alterations in 19 cases of BTC from primary tumor specimens. The most common molecular alterations were in TP53 (8/19, 42.1%), including missense mutations such as C242Y, E285K, G112S, P19T, R148T, R248Q, and R273L. We also detected two NRAS mutations (G12C and Q61L), two KRAS mutations (G12A and G12S), two ERBB2 mutations (V777L and pM774delinsMA) and amplification, and three PIK3CA mutations (N345K, E545K, and E521K). PDC models were successfully established in 27 of 40 samples (67.5%), including 22/24 from body fluids (91.7%) and 5/16 from tissue specimens (31.3%). Conclusions: PDC models are promising tools for uncovering driver mutations and identifying rational therapeutic strategies in BTC. Application of this model is expected to inform clinical trials of drugs for molecular-based targeted therapy.

BACKGROUND: Although pazopanib treatment has become the standard chemotherapy in salvage setting for metastatic sarcoma patients, most patients progress after pazopanib treatment in 4 to 6 months. After failure to pazopanib, patients have limited options for treatment. Therefore, subsequent therapy in patients who failed to pazopanib is urgently needed and the use of patient derived cells or patient derived tumors for accompanying testing with various pharmacological inhibitors could offer additional treatment options for these patients. METHODS: Patient derived tumor cells were collected from ascites at the time of progression to pazopanib and a 13-drug panel was tested for drug sensitivity. We confirmed the results using in vitro cell viability assay and immunoblot assay. We also performed the genomic profiling of PDX model. RESULTS: The growth of patient derived tumor cells was significantly reduced by exposure to 1.0 μM AZD2014 compared with control (control versus AZD2014, mean growth = 100.0% vs 16.04%, difference = 83.96%, 95% CI = 70.01% to 97.92%, P = .0435). Similarly, 1.0 μM BEZ235 profoundly inhibited tumor cell growth in vitro when compared to control (control versus BEZ235, mean growth = 100.0% vs 7.308%, difference = 92.69%, 95% CI = 78.87% to 106.5%, P < .0001). Despite the presence of CDK4 amplification in the patient-derived tumor cells, LEE011 did not considerably inhibit cell proliferation when compared with control (control vs LEE011, mean growth = 100.0% vs 80.23%, difference = 19.77%, 95% CI = 1.828% to 37.72%, P = .0377). The immunoblot analysis showed that BEZ235 treatment decreased pAKT, pmTOR and pERK whereas AZD2014 decreased only pmTOR. CONCLUSION: Taken together, upregulation of mTOR/AKT pathway in sarcoma patient derived cells was considerably inhibited by the treatment of AZD2014 and BEZ235 with downregulation of AKT pathway (greater extent for BEZ235). These molecules may be considered as treatment option in STS patient who have failed to pazopanib in the context of clinical trials.