Genomic progression for local invasion of cutaneous squamous cell carcinoma from the superficial to the deep portion.

Cutaneous squamous cell carcinoma (cSCC) is the second most common skin cancer. While many treatments exist, our understanding of its genomic progression, especially from the epidermis to the deep dermis, remains limited. This study aims to identify genetic mutations associated with the progression of cSCC into the deep dermis, providing insights into its aggressive behavior and high-risk features. We performed high-depth whole-exome sequencing on 12 cSCC tissues, along with paired normal tissues from six patients, using microdissection techniques. The mutational analysis focused on identifying alterations enriched during cSCC progression. Gene Ontology enrichment analysis, immunohistochemical assays, and external single-cell RNA data were utilized for validation. A total of 8863 non-synonymous somatic mutations were identified in 4092 genes across the superficial and deep portions of cSCCs. Analysis of deep portion mutations revealed a significant correlation with gene ontology biological processes, particularly cell junction organization, and cell-cell adhesion. Clonal mutations in these processes were more prevalent in the deep portions, indicating their impact on the cSCC mutation landscape. Genetic evolution analysis identified 29 causal genes associated with dermal invasion in cSCC. We highlight somatic mutations in cSCC, revealing heterogeneity between superficial and deep regions. Altered genes in cell junction organization and cell-cell adhesion emerged as pivotal in dermal invasion. We identified 29 causal genes primarily in deep tumor regions. Our findings emphasize analyzing multiple tumor regions to capture varied mutational landscapes. These insights advance our understanding of cSCC progression, emphasizing genetic and cellular changes during tumor evolution.

Risk Allele Frequency Analysis and Risk Prediction of Single-Nucleotide Polymorphisms for Prostate Cancer.

The incidence of prostate cancer (PCa) varies by ethnicity. This study aimed to provide insights into the genetic cause of PCa, which can result in differences in incidence among individuals of diverse ancestry. We collected data on PCa-associated single-nucleotide polymorphisms (SNPs) from a genome-wide association study catalog. Fisher's exact tests were used to analyze the significance of enrichment or depletion of the effect on the allele at a given SNP. A network analysis was performed based on PCa-related SNPs that showed significant differences among ethnicities. The SNP-based polygenic risk score (PRS) was calculated, and its correlation with PCa incidence was evaluated. European, African, and East Asian populations had different heatmap patterns. Calculated PRS from the allele frequencies of PCa was the highest among Africans, followed by Europeans, and was the lowest among East Asians. PRS was positively correlated with the incidence and mortality of PCa. Network analysis revealed that AR, CDKN1B, and MAD1L1 are genes related to ethnic differences in PCa. The incidence and mortality of PCa showed a strong correlation with PRS according to ethnicity, which may suggest the effect of genetic factors, such as the AR gene, on PCa pathogenesis.

Monitoring Circulating Tumor DNA in Untreated Non-Small-Cell Lung Cancer Patients.

Circulating tumor DNA (ctDNA) has been utilized to monitor the clinical course of patients of non-small-cell lung cancer (NSCLC) who receive therapies targeting druggable mutations. However, despite providing valuable information on how NSCLC would naturally progress, the clinical utility of ctDNA for clinical-course monitoring and prediction of treatment-naïve NSCLC patients without druggable mutations remain unknown. We longitudinally followed a total of 12 treatment-naïve NSCLC patients, who did not harbor EGFR and ALK mutations, by collecting clinical information, radiological data, and plasma samples. Changes in ctDNA levels and tumor burden (TB) were compared with each other. New metastasis development, volume doubling time (VDT), and overall survival (OS) were analyzed regarding ctDNA detection at diagnosis. ctDNA was detected in the plasma of seven (58.3%) patients. Changes in ctDNA levels correlated with those in TB in a substantial fraction (57.1%) of patients and was also associated with brain metastasis, tumor necrosis, or pneumonia in other patients. All patients with ctDNA detection developed new metastasis during follow-ups in the organs that had been devoid of metastasis at diagnosis. The patients without ctDNA detection did not develop new metastasis (median duration of follow-ups: 9.8 months). In addition, patients with ctDNA detection had shorter VDT (p = 0.039) and worse OS (p = 0.019) than those without ctDNA detection. The natural course of NSCLC progression can be monitored by measuring ctDNA levels. Detection of ctDNA at diagnosis can predict development of new metastasis, rapid tumor growth and poor survival of NSCLC patients.

Long-term efficacy and safety of intravenous injection of clonal mesenchymal stem cells derived from bone marrow in five adults with moderate to severe atopic dermatitis.

Atopic dermatitis is a chronic and relapsing inflammatory skin disease that is treated with immunosuppressants. However, long-term use of immunosuppressants may cause toxicity and severe side-effects. To confirm the long-term efficacy and safety of clonal mesenchymal stem cell therapy, we performed investigator-initiated clinical trials and long-term observation in five adult patients with moderate to severe atopic dermatitis that was refractory to conventional treatments. The clinical response assessment values such as Eczema Area and Severity Index (EASI) improved significantly at 16 weeks, and 80% (4/5) of the patients achieved EASI-50 after one or two treatment cycles. Patients were observed for long-term efficacy and safety for an average of 38 weeks (range, 16-86) and showed no serious side-effects. Among the cytokines tested, CCL-17, interleukin (IL)-13, and IL-22 significantly decreased at the end-point of the five participants, two patients who maintained good clinical response over 84 weeks showed increased IL-17 cytokine levels in the blood.

Dysregulation of cancer genes by recurrent intergenic fusions.

BACKGROUND: Gene fusions have been studied extensively, as frequent drivers of tumorigenesis as well as potential therapeutic targets. In many well-known cases, breakpoints occur at two intragenic positions, leading to in-frame gene-gene fusions that generate chimeric mRNAs. However, fusions often occur with intergenic breakpoints, and the role of such fusions has not been carefully examined. RESULTS: We analyze whole-genome sequencing data from 268 patients to catalog gene-intergenic and intergenic-intergenic fusions and characterize their impact. First, we discover that, in contrast to the common assumption, chimeric oncogenic transcripts-such as those involving ETV4, ERG, RSPO3, and PIK3CA-can be generated by gene-intergenic fusions through splicing of the intervening region. Second, we find that over-expression of an upstream or downstream gene by a fusion-mediated repositioning of a regulatory sequence is much more common than previously suspected, with enhancers sometimes located megabases away. We detect a number of recurrent fusions, such as those involving ANO3, RGS9, FUT5, CHI3L1, OR1D4, and LIPG in breast; IGF2 in colon; ETV1 in prostate; and IGF2BP3 and SIX2 in thyroid cancers. CONCLUSION: Our findings elucidate the potential oncogenic function of intergenic fusions and highlight the wide-ranging consequences of structural rearrangements in cancer genomes.

Lineage-dependent gene expression programs influence the immune landscape of colorectal cancer.

Immunotherapy for metastatic colorectal cancer is effective only for mismatch repair-deficient tumors with high microsatellite instability that demonstrate immune infiltration, suggesting that tumor cells can determine their immune microenvironment. To understand this cross-talk, we analyzed the transcriptome of 91,103 unsorted single cells from 23 Korean and 6 Belgian patients. Cancer cells displayed transcriptional features reminiscent of normal differentiation programs, and genetic alterations that apparently fostered immunosuppressive microenvironments directed by regulatory T cells, myofibroblasts and myeloid cells. Intercellular network reconstruction supported the association between cancer cell signatures and specific stromal or immune cell populations. Our collective view of the cellular landscape and intercellular interactions in colorectal cancer provide mechanistic information for the design of efficient immuno-oncology treatment strategies.

Junction Location Identifier (JuLI): Accurate Detection of DNA Fusions in Clinical Sequencing for Precision Oncology.

Accurate detection of genomic fusions by high-throughput sequencing in clinical samples with inadequate tumor purity and formalin-fixed, paraffin-embedded tissue is an essential task in precise oncology. We developed the fusion detection algorithm Junction Location Identifier (JuLI) for optimization of high-depth clinical sequencing. Novel filtering steps were implemented to minimize false positives in the clinical setting. The algorithm was comprehensively validated using high-depth sequencing data from cancer cell lines and clinical samples and genome sequencing data from NA12878. JuLI showed improved performance mainly in positive predictive value over state-of-the-art fusion callers in cases with high-depth clinical sequencing and rescued a driver fusion from false negative in plasma cell-free DNA using joint calling.

Evolutionary processes of melanomas from giant congenital melanocytic nevi.

Melanoma can develop in a congenital melanocytic nevus (CMN). In fact, a large CMN is associated with a high risk of developing melanoma. Although melanomas arising from CMNs are thought to have a pathogenesis distinct from conventional melanomas, no studies have been conducted on the evolution or tumor heterogeneity of CMN melanomas. We applied multi-region whole-exome sequencing to investigate the clonal nature of driver events and evolutionary processes in CMNs and melanomas arising from CMNs. In two patients, we observed an independent subclonal evolution in cancerized fields of CMNs and chromosome 8q amplification in both melanomas arising from CMNs. The amplification of MYC, located in chromosome 8q, was correlated with the percentage of tumor cells expressing high levels of MYC protein detected in melanoma cells by immunohistochemistry. Our analysis suggests that each CMN cell may evolve sporadically and that amplification of MYC might be a key event for melanoma development in CMNs.

Assessment of intratumoral heterogeneity with mutations and gene expression profiles.

Intratumoral heterogeneity (ITH) refers to the presence of distinct tumor cell populations. It provides vital information for the clinical prognosis, drug responsiveness, and personalized treatment of cancer patients. As genomic ITH in various cancers affects the expression patterns of genes, the expression profile could be utilized for determining ITH level. Herein, we present a novel approach to directly detect high ITH defined as a larger number of subclones from the gene expression pattern through machine learning approaches. We examined associations between gene expression profile and ITH of 12 cancer types from The Cancer Genome Atlas (TCGA) database. Using stomach adenocarcinoma (STAD) showing high association, we evaluated the performance of our method in predicting ITH by employing three machine learning algorithms using gene expression profile data. We classified tumors into high and low heterogeneity groups using the learning model through the selection of LASSO feature. The result showed that support vector machines (SVMs) outperformed other algorithms (AUC = 0.84 in SVMs and 0.82 in Naïve Bayes) and we were able to improve predictive power by using both combined data from mutation and expression. Furthermore, we evaluated the prediction ability of each model using simulation data generated by mixing cell lines of the Cancer Cell Line Encyclopedia (CCLE), and obtained consistent results with using real dataset. Our approach could be utilized for discriminating tumors with heterogeneous cell populations to characterize ITH.

The concept of nail matrix onychodermis (onychomatricodermis) in the nail unit: Histology and elastin immunohistochemistry.

BACKGROUNDS: We previously demonstrated the presence of onychodermis below nail matrix and nail bed. Because nail matrix is a producer of nail plate, we hypothesized that onychodermis below nail matrix could be the nail counterpart of follicular dermal papilla. In this study, we sought to further characterize histologic, histochemical, and immunohistochemical features of nail matrix onychodermis. METHODS AND RESULTS: Hematoxylin and eosin slides of 10 polydactyly nail units and 10 nail matrix biopsies from children and adults were reviewed. In polydactyly nail units, the onychodermis beneath nail matrix was characterized by onychofibroblasts showing abundant cytoplasm, and this area was slightly separated from the undersurface of the nail matrix. Nail matrix biopsy specimens also showed similar histology in the nail matrix onychodermis. Alcian blue stain demonstrated mucin deposition in onychofibroblasts within the nail matrix onychodermis. Immunohistochemically, elastin was rarely expressed in the nail matrix onychodermis while it was strongly expressed in the dermis of other areas of polydactyly nail units. Elastin was not expressed in follicular dermal papilla of terminal hair follicles of the scalp. CONCLUSION: Our findings demonstrate the presence and localization of nail matrix onychodermis (onychomatricodermis). Our study also demonstrates similar elastin expression patterns in the onychomatricodermis and follicular dermal papilla.

Intratumor heterogeneity inferred from targeted deep sequencing as a prognostic indicator.

Tumor genetic heterogeneity may underlie poor clinical outcomes because diverse subclones could be comprised of metastatic and drug resistant cells. Targeted deep sequencing has been used widely as a diagnostic tool to identify actionable mutations in cancer patients. In this study, we evaluated the clinical utility of estimating tumor heterogeneity using targeted panel sequencing data. We investigated the prognostic impact of a tumor heterogeneity (TH) index on clinical outcomes, using mutational profiles from targeted deep sequencing data acquired from 1,352 patients across 8 cancer types. The TH index tended to be increased in high pathological stage disease in several cancer types, indicating clonal expansion of cancer cells as tumor progression proceeds. In colorectal cancer patients, TH index values also correlated significantly with clinical prognosis. Integration of the TH index with genomic and clinical features could improve the power of risk prediction for clinical outcomes. In conclusion, deep sequencing to determine the TH index could serve as a promising prognostic indicator in cancer patients.

Clinical Relevance of Genomic Changes in Recurrent Pediatric Solid Tumors.

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 (<10%). Patients were classified into three categories based on the mutation patterns after cancer treatment. A significant association between the major mutation patterns and clinical outcome was observed. Patients whose relapsed tumor had fewer mutations than the diagnostic sample tended to be older, had longer progression-free survival, and achieved complete remission after relapse. Contrastingly, patients whose genetic profile only had concordant mutations without any change had the worst outcome. CONCLUSIONS: We characterized genomic changes in recurrent pediatric solid tumors. These findings could help to understand the biology of relapsed childhood cancer and to develop personalized treatment based on their genetic profile.

Integrin ß3 Inhibition Enhances the Antitumor Activity of ALK Inhibitor in ALK-Rearranged NSCLC.

Purpose: Anaplastic lymphoma kinase (ALK)-positive cancers are sensitive to small-molecule ALK kinase inhibitors, but most cases experience failure following treatment. Hence, additional drug targets and combination therapeutic treatments are needed. We investigated gene expression that is regulated by the expression of ALK and explored its roles in cancer progression and therapeutic implication.Experimental Design: We screened ALK-rearranged non-small cell lung cancer (NSCLC) cases using immunohistochemistry and fluorescence in situ hybridization and then conducted multiplex gene expression analysis. We also performed a clinicopathologic analysis to validate the findings. Additional cellular experiments, including inhibition and migration assays, and in vivo lung cancer model studies were performed.Results: Among patients with ALK-rearranged NSCLC, integrin ß3 (ITGB3) was one of the overexpressed genes in comparison with that in ALK-negative NSCLC (P = 0.0003). ALK and integrin ß3 expression were positively correlated, and we discovered that high integrin ß3 mRNA expression was associated with metastasis and more advanced tumor stages (P < 0.005; P < 0.05). Furthermore, we found that inhibition of both ALK and integrin ß3 led to increased drug sensitivity in vitro and in vivo (both P < 0.05).Conclusions: We discovered a positive correlation between ALK and integrin ß3 expression levels in ALK-rearranged NSCLC. Our findings suggest that high integrin ß3 expression in ALK-rearranged NSCLC is associated with tumor progression and a worse prognosis. This finding demonstrates the prognostic value of integrin ß3 and provides a rationale for combination treatment with ALK and integrin ß3 inhibitors in patients with ALK-rearranged NSCLC. Clin Cancer Res; 24(17); 4162-74. ©2018 AACR.

Genomic alterations of ground-glass nodular lung adenocarcinoma.

In-depth molecular pathogenesis of ground-glass nodular lung adenocarcinoma has not been well understood. The objectives of this study were to identify genomic alterations in ground-glass nodular lung adenocarcinomas and to investigate whether viral transcripts were detected in these tumors. Nine patients with pure (n = 4) and part-solid (n = 5) ground-glass nodular adenocarcinomas were included. Six were females with a median age of 58 years. We performed targeted exon sequencing and RNA sequencing. EGFR (n = 10), IDH2 (n = 2), TP53 (n = 1), PTEN (n = 1), EPHB4 (n = 1), and BRAF (n = 1) were identified as driver mutations by targeted exon sequencing. Vasculogenesis-associated genes including NOTCH4 and TGFBR3 expression were significantly downregulated in adenocarcinoma tissue versus normal tissue (adjusted P values < 0.001 for both NOTCH4 and TGFBR3). In addition, five novel fusion gene loci were identified in four lung adenocarcinomas. However, no significant virus-associated transcripts were detected in tumors. In conclusions, EGFR, IDH2, TP53, PTEN, EPHB4, and BRAF were identified as putative driver mutations of ground-glass nodular adenocarcinomas. Five novel fusion genes were also identified in four tumors. Viruses do not appear to be involved in the tumorigenesis of ground-glass nodular lung adenocarcinoma.

Prevalence and detection of low-allele-fraction variants in clinical cancer samples.

Accurate detection of genomic alterations using high-throughput sequencing is an essential component of precision cancer medicine. We characterize the variant allele fractions (VAFs) of somatic single nucleotide variants and indels across 5095 clinical samples profiled using a custom panel, CancerSCAN. Our results demonstrate that a significant fraction of clinically actionable variants have low VAFs, often due to low tumor purity and treatment-induced mutations. The percentages of mutations under 5% VAF across hotspots in EGFR, KRAS, PIK3CA, and BRAF are 16%, 11%, 12%, and 10%, respectively, with 24% for EGFR T790M and 17% for PIK3CA E545. For clinical relevance, we describe two patients for whom targeted therapy achieved remission despite low VAF mutations. We also characterize the read depths necessary to achieve sensitivity and specificity comparable to current laboratory assays. These results show that capturing low VAF mutations at hotspots by sufficient sequencing coverage and carefully tuned algorithms is imperative for a clinical assay.

Characterization of background noise in capture-based targeted sequencing data.

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.

Molecular breakdown: a comprehensive view of anaplastic lymphoma kinase (ALK)-rearranged non-small cell lung cancer.

Most anaplastic lymphoma kinase (ALK)-rearranged non-small cell lung cancers (NSCLCs) show good clinical response to ALK inhibitors. However, some ALK-rearranged NSCLC patients show various primary responses with unknown reasons. Previous studies focused on the clinical aspects of ALK fusions in small cohorts, or were conducted in vitro and/or in vivo to investigate the function of ALK. One of the suggested theories describes how echinoderm microtubule-associated protein-like 4 (EML4)-ALK variants play a role towards different sensitivities in ALK inhibitors. Until now, there has been no integrated comprehensive study that dissects ALK at the molecular level in a large scale. Here, we report the largest extensive molecular analysis of 158 ALK-rearranged NSCLCs and have investigated these findings in a cell line construct experiment. We discovered that NSCLCs with EML4-ALK short forms (variant 3/others) had more advanced stage and frequent metastases than cases with the long forms (variant 1/others) (p = 0.057, p < 0.05). In vitro experiments revealed that EML4-ALK short forms show lower sensitivity to ALK inhibitors than do long forms. Clinical analysis also showed a trend for the short forms showing worse PFS. Interestingly, we found that breakpoints of ALK are evenly distributed mainly in intron 19 and almost all of them undergo a non-homologous end-joining repair to generate ALK fusions. We also discovered four novel somatic ALK mutations in NSCLC (T1151R, R1192P, A1280V, and L1535Q) that confer primary resistance; all of them showed strong resistance to ALK inhibitors, as G1202R does. Through targeted deep sequencing, we discovered three novel ALK fusion partners (GCC2, LMO7, and PHACTR1), and different ALK fusion partners showed different intracellular localization. With our findings that the EML4-ALK variants, new ALK somatic mutations, and novel ALK-fusion partners may affect sensitivity to ALK inhibitors, we stress the importance of targeted therapy to take the ALK molecular profiling into consideration. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Mutational profiling of acral melanomas in Korean populations.

The proportion of acral melanoma (AM) is much higher in Asian populations than in Caucasian populations. Although mutational profiles associated with AM have been discovered in Caucasian populations, knowledge of its genetic alterations in Asian populations is limited. To describe the molecular nature of AM in Korean patients, we performed mutational profiling of AM and matched normal tissues in patients. Fifty-one formalin-fixed paraffin-embedded AM samples and 32 matched pairs from patients' saliva DNA were analysed by next-generation sequencing. Only mutations confirmed via digital droplet PCR or in BRAF, KIT and NRAS, the most frequently altered cancer genes in cutaneous melanoma, were considered as positive. The relationship between mutational status and clinicopathological features were examined. Of the 47 AM patients screened, alteration of BRAF, NRAS and KIT genes was observed in 6.4%, 4.3% and 8.5%, respectively. We also tested matched normal tissues of patients to identify tumor-specific mutations. Examination of the mutational profile in a cohort of 28 primary melanomas and matched normal controls found BRAF mutations in two cases (7.1%), KIT mutations in three cases (10.7%) and CTNNB1 mutations in one case (3.6%). The BRAF, NRAS and KIT mutation status did not correlate with clinicopathological characteristics. Our results show that KIT, NRAS and BRAF hotspot mutations occur at a low frequency in Korean populations. We also observed a case with the CTNNB1 mutation, which raises the possibility that other pathways are associated with AM development.

Recurrent mutations of MAPK pathway genes in multiple myeloma but not in amyloid light-chain amyloidosis.

Clinically applicable platforms revealing actionable genomic alterations may improve the treatment efficacy of myeloma patients. In this pilot study, we used a high depth targeted sequencing panel containing 83 anti-cancer drug target genes to sequence genomic DNAs extracted from bone marrow aspirates of 23 patients with myeloma and 12 patients with amyloid light-chain (AL) amyloidosis. Mutation analysis revealed NRAS as the most commonly mutated gene (30%, 7/23) in myeloma patients followed by KRAS (26%, 6/23) and BRAF (22%, 5/23). However, no significant mutations were found in the 12 patients with AL amyloidosis. Notably, 6 of the 23 myeloma patients showed multi-site and/or multi-gene mutations in NRAS, KRAS, or BRAF, indicating compound aberrations in the Mitogen activated protein kinase (MAPK) pathway. Gene panel sequencing also revealed cytogenetic abnormalities associated with prognosis in myeloma patients. In conclusion, our pilot study suggests that targeted gene sequencing may have an important prognostic value for myeloma patients for the identification of actionable genomic alterations and cytogenetic aberrations.