Mutational evolution after chemotherapy-progression in metastatic colorectal cancer revealed by circulating tumor DNA analysis.

Emerging new mutations after treatment can provide clues to acquired resistant mechanisms. Circulating tumor DNA (ctDNA) sequencing has enabled noninvasive repeated tumor mutational profiling. We aimed to investigate newly emerging mutations in ctDNA after disease progression in metastatic colorectal cancer (mCRC). Blood samples were prospectively collected from mCRC patients receiving palliative chemotherapy before treatment and at radiological evaluations. ctDNA from pretreatment and progressive disease (PD) samples were sequenced with a next-generation sequencing panel targeting 106 genes. A total of 712 samples from 326 patients were analyzed, and 381 pretreatment and PD pairs (163 first-line, 85 second-line and 133 later-line [≥third-line]) were compared. New mutations in PD samples (mean 2.75 mutations/sample) were observed in 49.6% (189/381) of treatments. ctDNA samples from later-line had more baseline mutations (P = .002) and were more likely to have new PD mutations (adjusted odds ratio [OR] 2.27, 95% confidence interval [CI]: 1.40-3.69) compared to first-line. RAS/BRAF wild-type tumors were more likely to develop PD mutations (adjusted OR 1.87, 95% CI: 1.22-2.87), independent of cetuximab treatment. The majority of new PD mutations (68.5%) were minor clones, suggesting an increasing clonal heterogeneity after treatment. Pathways involved by PD mutations differed by the treatment received: MAPK cascade (Gene Ontology [GO]: 0000165) in cetuximab and regulation of kinase activity (GO: 0043549) in regorafenib. The number of mutations revealed by ctDNA sequencing increased during disease progression in mCRC. Clonal heterogeneity increased after chemotherapy progression, and pathways involved were affected by chemotherapy regimens.

Prevalence and founder effect of DRC1 exon 1-4 deletion in Korean patients with primary ciliary dyskinesia.

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous disorder affecting ciliary structure and function. PCD exhibiting dynein regulatory complex subunit 1 (DRC1) exon 1-4 deletion has been reported in several Japanese PCD patients; however, no large scale studies have been performed. Here, we aimed to determine the prevalence and founder effect of this variant in the Korean population. Using an in-house copy number variation tool, we screened for DRC1 exon 1-4 deletion in 20 patients with PCD and exome data of 1435 patients in the Seoul National University Hospital repository. In cases of suspected DRC1 deletion, confirmatory gap-PCR was performed. In a PCD cohort, three of 20 (15%) patients were positive for DRC1 exon 1-4 deletion (NM_145038.5(DRC1): c.1-3952_540 + 1331del27748-bp) while pathogenic variants were found in CCDC39 (N = 1), DNAAF6 (N = 1), DNAH9 (N = 1). In the 1,435-sample exome data, seven patients (0.49%) were confirmed to have DRC1 exon 1-4 deletion. A chimeric sequence including the junction was searched from the 1000 Genomes Project data repository. One Japanese patient (0.96%) was found to have the same DRC1 exon 1-4 deletion, which was absent in other populations. This study demonstrated that the DRC1 exon 1-4 deletion is a founder mutation based on haplotype analysis. In summary, the prevalence of PCD based on DRC1 exon 1-4 deletion is particularly high in Korean and Japanese populations, which is attributed to the founder effect. Genetic testing for DRC1 exon 1-4 deletion should be considered as an initial screening tool for Korean and Japanese patients with PCD.

Dynamic changes in longitudinal circulating tumour DNA profile during metastatic colorectal cancer treatment.

BACKGROUND: Circulating tumour DNA (ctDNA) has been spotlighted as an attractive biomarker because of its easy accessibility and real-time representation of tumour genetic profile. However, the clinical utility of longitudinal ctDNA monitoring has not been clearly defined. METHODS: Serial blood samples were obtained from metastatic colorectal cancer patients undergoing first-line chemotherapy. ctDNA was sequenced using a targeted next-generation sequencing platform which included 106 genes. Changes in ctDNA profile and treatment outcome were comprehensively analysed. RESULTS: A total of 272 samples from 62 patients were analysed. In all, 90.3% of patients had detectable ctDNA mutation before treatment. ctDNA clearance after chemotherapy was associated with longer progression-free survival which was independent of radiological response (adjusted hazard ratio 0.22, 95% confidence interval 0.10-0.46). Longitudinal monitoring was able to detect ctDNA progression which preceded radiological progressive disease (PD) in 58.1% (median 3.3 months). Diverse resistant mutations (34.9%) and gene amplification (7.0%) at the time of PD were discovered. For 16.3% of the PD patients, the newly identified mutations could be potential candidates of targeted therapy or clinical trial. CONCLUSION: ctDNA profile provided a more accurate landscape of tumour and dynamic changes compared to radiological evaluation. Longitudinal ctDNA monitoring may improve personalised treatment decision-making.

A highly annotated whole-genome sequence of a Korean individual.

Recent advances in sequencing technologies have initiated an era of personal genome sequences. To date, human genome sequences have been reported for individuals with ancestry in three distinct geographical regions: a Yoruba African, two individuals of northwest European origin, and a person from China. Here we provide a highly annotated, whole-genome sequence for a Korean individual, known as AK1. The genome of AK1 was determined by an exacting, combined approach that included whole-genome shotgun sequencing (27.8x coverage), targeted bacterial artificial chromosome sequencing, and high-resolution comparative genomic hybridization using custom microarrays featuring more than 24 million probes. Alignment to the NCBI reference, a composite of several ethnic clades, disclosed nearly 3.45 million single nucleotide polymorphisms (SNPs), including 10,162 non-synonymous SNPs, and 170,202 deletion or insertion polymorphisms (indels). SNP and indel densities were strongly correlated genome-wide. Applying very conservative criteria yielded highly reliable copy number variants for clinical considerations. Potential medical phenotypes were annotated for non-synonymous SNPs, coding domain indels, and structural variants. The integration of several human whole-genome sequences derived from several ethnic groups will assist in understanding genetic ancestry, migration patterns and population bottlenecks.

Pursuing dynamics of minimal residual leukemic subclones in relapsed and refractory acute myeloid leukemia during conventional therapy.

BACKGROUND: Acute myeloid leukemia (AML) is characterized by clonal heterogeneity, leading to frequent relapses and drug resistance despite intensive clinical therapy. Although AML's clonal architecture has been addressed in many studies, practical monitoring of dynamic changes in those subclones during relapse and treatment is still understudied. METHOD: Fifteen longitudinal bone marrow (BM) samples were collected from three relapsed and refractory (R/R) AML patients. Using droplet digital polymerase chain reaction (ddPCR), the frequencies of patient's leukemic variants were assessed in seven cell populations that were isolated from each BM sample based on cellular phenotypes. By quantifying mutant clones at the diagnosis, remission, and relapse stages, the distribution of AML subclones was sequentially monitored. RESULTS: Minimal residual (MR) leukemic subclones exhibit heterogeneous distribution among BM cell populations, including mature leukocyte populations. During AML progression, these subclones undergo active phenotypic transitions and repopulate into distinct cell population regardless of normal hematopoiesis hierarchic order. Of these, MR subclones in progenitor populations of patient BM predominantly carry MR leukemic properties, leading to more robust expansion and stubborn persistence than those in mature populations. Moreover, a minor subset of MR leukemic subclones could be sustained at an extremely low frequency without clonal expansion during relapse. CONCLUSIONS: In this study, we observed treatment persistent MR leukemic subclones and their phenotypic changes during the treatment process of R/R AML patients. This underscores the importance of preemptive inhibition of clonal promiscuity in R/R AML, proposing a practical method for monitoring AML MR subclones.

Diagnostic uplift through the implementation of short tandem repeat analysis using exome sequencing.

To date, approximately 50 short tandem repeat (STR) disorders have been identified; yet, clinical laboratories rarely conduct STR analysis on exomes. To assess its diagnostic value, we analyzed STRs in 6099 exomes from 2510 families with mostly suspected neurogenetic disorders. We employed ExpansionHunter and REViewer to detect pathogenic repeat expansions, confirming them using orthogonal methods. Genotype-phenotype correlations led to the diagnosis of thirteen individuals in seven previously undiagnosed families, identifying three autosomal dominant disorders: dentatorubral-pallidoluysian atrophy (n = 3), spinocerebellar ataxia type 7 (n = 2), and myotonic dystrophy type 1 (n = 2), resulting in a diagnostic gain of 0.28% (7/2510). Additionally, we found expanded ATXN1 alleles (≥39 repeats) with varying patterns of CAT interruptions in twelve individuals, accounting for approximately 0.19% in the Korean population. Our study underscores the importance of integrating STR analysis into exome sequencing pipeline, broadening the application of exome sequencing for STR assessments.

Broadening the scope of multigene panel analysis for adult epilepsy patients.

OBJECTIVE: Epilepsy is a suitable target for gene panel sequencing because a considerable portion of epilepsy is now explained by genetic components, especially in syndromic cases. However, previous gene panel studies on epilepsy have mostly focused on pediatric patients. METHODS: We enrolled adult epilepsy patients meeting any of the following criteria: family history of epilepsy, seizure onset age ≤ 19 years, neuronal migration disorder, and seizure freedom not achieved by dual anti-seizure medications. We sequenced the exonic regions of 211 epilepsy genes in these patients. To confirm the pathogenicity of a novel MTOR truncating variant, we electroporated vectors with different MTOR variants into developing mouse brains. RESULTS: A total of 92 probands and 4 affected relatives were tested, and the proportion of intellectual disability (ID) and/or developmental disability (DD) was 21.7%. As a result, twelve probands (13.0%) had pathogenic or likely pathogenic variants in the following genes or regions: DEPDC5, 15q12-q13 duplication (n = 2), SLC6A1, SYNGAP1, EEF1A2, LGI1, MTOR, KCNQ2, MEF2C, and TSC1 (n = 1). We confirmed the functional impact of a novel truncating mutation in the MTOR gene (c.7570C > T, p.Gln2524Ter) that disrupted neuronal migration in a mouse model. The diagnostic yield was higher in patients with ID/DD or childhood-onset seizures. We also identified additional candidate variants in 20 patients that could be reassessed by further studies. SIGNIFICANCE: Our findings underscore the clinical utility of gene panel sequencing in adult epilepsy patients suspected of having genetic etiology, especially those with ID/DD or early-onset seizures. Gene panel sequencing could not only lead to genetic diagnosis in a substantial portion of adult epilepsy patients but also inform more precise therapeutic decisions based on their genetic background. PLAIN LANGUAGE SUMMARY: This study demonstrated the effectiveness of gene panel sequencing in adults with epilepsy, revealing pathogenic or likely pathogenic variants in 13.0% of patients. Higher diagnostic yields were observed in those with neurodevelopmental disorders or childhood-onset seizures. Additionally, we have shown that expanding genetic studies into adult patients would uncover new types of pathogenic variants for epilepsy, contributing to the advancement of precision medicine for individuals with epilepsy. In conclusion, our results highlight the practical value of employing gene panel sequencing in adult epilepsy patients, particularly when genetic etiology is clinically suspected.

Prevalence and Characterization of NOTCH2NLC GGC Repeat Expansions in Koreans: From a Hospital Cohort Analysis to a Population-Wide Study.

Background and Objectives: GGC repeat expansions in the NOTCH2NLC gene are associated with a broad spectrum of progressive neurologic disorders, notably, neuronal intranuclear inclusion disease (NIID). We aimed to investigate the population-wide prevalence and clinical manifestations of NOTCH2NLC-related disorders in Koreans. Methods: We conducted a study using 2 different cohorts from the Korean population. Patients with available brain MRI scans from Seoul National University Hospital (SNUH) were thoroughly reviewed, and NIID-suspected patients presenting the zigzag edging signs underwent genetic evaluation for NOTCH2NLC repeats by Cas9-mediated nanopore sequencing. In addition, we analyzed whole-genome sequencing data from 3,887 individuals in the Korea Biobank cohort to estimate the distribution of the repeat counts in Koreans and to identify putative patients with expanded alleles and neurologic phenotypes. Results: In the SNUH cohort, among 90 adult-onset leukoencephalopathy patients with unknown etiologies, we found 20 patients with zigzag edging signs. Except for 2 diagnosed with fragile X-associated tremor/ataxia syndrome and 2 with unavailable samples, all 16 patients (17.8%) were diagnosed with NIID (repeat range: 87-217). By analyzing the Korea Biobank cohort, we estimated the distribution of repeat counts and threshold (>64) for Koreans, identifying 6 potential patients with NIID. Furthermore, long-read sequencing enabled the elucidation of transmission and epigenetic patterns of NOTCH2NLC repeats within a family affected by pediatric-onset NIID. Discussion: This study presents the population-wide distribution of NOTCH2NLC repeats and the estimated prevalence of NIID in Koreans, providing valuable insights into the association between repeat counts and disease manifestations in diverse neurologic disorders.

TIARA: a database for accurate analysis of multiple personal genomes based on cross-technology.

High-throughput genomic technologies have been used to explore personal human genomes for the past few years. Although the integration of technologies is important for high-accuracy detection of personal genomic variations, no databases have been prepared to systematically archive genomes and to facilitate the comparison of personal genomic data sets prepared using a variety of experimental platforms. We describe here the Total Integrated Archive of Short-Read and Array (TIARA; http://tiara.gmi.ac.kr) database, which contains personal genomic information obtained from next generation sequencing (NGS) techniques and ultra-high-resolution comparative genomic hybridization (CGH) arrays. This database improves the accuracy of detecting personal genomic variations, such as SNPs, short indels and structural variants (SVs). At present, 36 individual genomes have been archived and may be displayed in the database. TIARA supports a user-friendly genome browser, which retrieves read-depths (RDs) and log2 ratios from NGS and CGH arrays, respectively. In addition, this database provides information on all genomic variants and the raw data, including short reads and feature-level CGH data, through anonymous file transfer protocol. More personal genomes will be archived as more individuals are analyzed by NGS or CGH array. TIARA provides a new approach to the accurate interpretation of personal genomes for genome research.

Utilizing Plasma Circulating Tumor DNA Sequencing for Precision Medicine in the Management of Solid Cancers.

Plasma circulating tumor DNA (ctDNA) sequencing has demonstrated clinical utility for tumor molecular profiling at initial diagnosis or tumor progression in advanced solid cancers and is being rapidly incorporated into the clinical practice guidelines, including non-small cell lung and breast cancer. Despite relatively low sensitivity, plasma ctDNA sequencing has several advantages over tissue-based assays, including ease of sampling, rapid turnaround time, repeatability, and the ability to overcome spatial heterogeneity, which makes it ideal for investigating acquired resistance and monitoring tumor evolution and dynamics. With technological advancement and declining costs, the clinical application of plasma ctDNA is expanding, and numerous ongoing clinical trials are examining its potential to guide the management of advanced, localized, and even preclinical cancers of various tumor types. The ability of plasma ctDNA analysis to detect minimal residual disease following curative treatment in the absence of clinical disease is among its most promising attributes. Plasma ctDNA sequencing can also facilitate the conduct of clinical trials and drug development, particularly in immunotherapy. In order to incorporate plasma ctDNA sequencing for clinical decision-making, it is important to understand the preanalytical and analytical factors that may affect its sensitivity and reliability.

Efficacy of Olaparib in Treatment-Refractory, Metastatic Breast Cancer with Uncommon Somatic BRCA Mutations Detected in Circulating Tumor DNA.

Poly(ADP-ribose) polymerase inhibitors have been shown dramatic responses in patients with BRCAness. However, clinical studies have been limited to breast cancer patients with germline mutations. Here, we describe a patient with metastatic breast cancer who had a rare BRCA1 somatic mutation (BRCA1 c.4336G>T (p.E1446*)) detected by cell-free DNA analysis after failing standard therapies. This tier III variant of unknown significance was predicted to be a pathogenic variant in our assessment, leading us to consider off-label treatment with olaparib. The patient responded well to olaparib for several months, with a decrease in allele frequency of this BRCA1 somatic mutation in cell-free DNA. Olaparib resistance subsequently developed with an increase in the allele frequency and new BRCA1 reversion mutations. To our knowledge, this is the first report confirming BRCA1 c.4336G>T (p.E1446*) as a mutation sensitive to olaparib in breast cancer and describing the dynamic changes in the associated mutations using liquid biopsy.

Expansion of clinico-genetic spectrum of PRDX3 disease: a literature review with two additional cases.

Graphical Abstract.

Clinical validity of oncogenic driver genes detected from circulating tumor DNA in the blood of lung cancer patients.

Background: This study aimed to evaluate the concordance of oncogenic driver mutations between tumor tissues and circulating tumor DNA (ctDNA) in patients with lung cancer. In addition, this study attempted to reveal the clinical utility of ctDNA in lung cancer treatment. Methods: Recurrent or metastatic non-small cell lung cancer (NSCLC) patients were prospectively enrolled in this study. Tumor tissue and serial blood samples were obtained from newly diagnosed patients (Cohort A) or patients treated with targeted therapy (Cohort B) and targeted gene panel sequencing was conducted to identify tumor mutational profiles. Results: At the time of diagnosis, patients in Cohort A with a high cell-free DNA (cfDNA) concentration had poorer overall survival than those with a low cfDNA concentration. The sensitivity and precision of ctDNA analysis in pre-treatment patients compared with those of tissue sequencing were 58.4% and 61.5%, respectively. Known lung cancer-associated variants of oncogenic driver genes, including EGFR and KRAS, and tumor suppressor genes, including TP53 and APC, were frequently detected in the ctDNA of the patients (76.9%). An association between smoking and TP53 mutation status was observed in both tissues and ctDNA (P=0.005 and 0.037, respectively). In addition, the EGFR T790M resistance mutation was detected solely from the ctDNA of two patients after treatment with an EGFR tyrosine kinase inhibitor. Conclusions: ctDNA may be a reliable prognostic biomarker with an additional role in treating patients with lung cancer. Further analyses are necessary to understand the properties of ctDNA and widen its clinical use.

Clinical Practice Recommendations for the Use of Next-Generation Sequencing in Patients with Solid Cancer: A Joint Report from KSMO and KSP.

In recent years, next-generation sequencing (NGS)-based genetic testing has become crucial in cancer care. While its primary objective is to identify actionable genetic alterations to guide treatment decisions, its scope has broadened to encompass aiding in pathological diagnosis and exploring resistance mechanisms. With the ongoing expansion in NGS application and reliance, a compelling necessity arises for expert consensus on its application in solid cancers. To address this demand, the forthcoming recommendations not only provide pragmatic guidance for the clinical use of NGS but also systematically classify actionable genes based on specific cancer types. Additionally, these recommendations will incorporate expert perspectives on crucial biomarkers, ensuring informed decisions regarding circulating tumor DNA panel testing.

Reference-unbiased copy number variant analysis using CGH microarrays.

Comparative genomic hybridization (CGH) microarrays have been used to determine copy number variations (CNVs) and their effects on complex diseases. Detection of absolute CNVs independent of genomic variants of an arbitrary reference sample has been a critical issue in CGH array experiments. Whole genome analysis using massively parallel sequencing with multiple ultra-high resolution CGH arrays provides an opportunity to catalog highly accurate genomic variants of the reference DNA (NA10851). Using information on variants, we developed a new method, the CGH array reference-free algorithm (CARA), which can determine reference-unbiased absolute CNVs from any CGH array platform. The algorithm enables the removal and rescue of false positive and false negative CNVs, respectively, which appear due to the effects of genomic variants of the reference sample in raw CGH array experiments. We found that the CARA remarkably enhanced the accuracy of CGH array in determining absolute CNVs. Our method thus provides a new approach to interpret CGH array data for personalized medicine.

Patient-derived organoids as a preclinical platform for precision medicine in colorectal cancer.

Patient-derived organoids are being considered as models that can help guide personalized therapy through in vitro anticancer drug response evaluation. However, attempts to quantify in vitro drug responses in organoids and compare them with responses in matched patients remain inadequate. In this study, we investigated whether drug responses of organoids correlate with clinical responses of matched patients and disease progression of patients. Organoids were established from 54 patients with colorectal cancer who (except for one patient) did not receive any form of therapy before, and tumor organoids were assessed through whole-exome sequencing. For comparisons of in vitro drug responses in matched patients, we developed an 'organoid score' based on the variable anticancer treatment responses observed in organoids. Very interestingly, a higher organoid score was significantly correlated with a lower tumor regression rate after the standard-of-care treatment in matched patients. Additionally, we confirmed that patients with a higher organoid score (≥ 2.5) had poorer progression-free survival compared with those with a lower organoid score (< 2.5). Furthermore, to assess potential drug repurposing using an FDA-approved drug library, ten tumor organoids derived from patients with disease progression were applied to a simulation platform. Taken together, organoids and organoid scores can facilitate the prediction of anticancer therapy efficacy, and they can be used as a simulation model to determine the next therapeutic options through drug screening. Organoids will be an attractive platform to enable the implementation of personalized therapy for colorectal cancer patients.

Circulating tumor DNA sequencing in colorectal cancer patients treated with first-line chemotherapy with anti-EGFR.

Circulating tumor DNA (ctDNA) may reveal dynamic tumor status during therapy. We conducted serial ctDNA analysis to investigate potential association with clinical outcome in metastatic colorectal cancer (mCRC) patients receiving chemotherapy. Tissue KRAS/NRAS wild-type mCRC patients were enrolled and treated with first-line cetuximab-containing chemotherapy. ctDNA isolated from plasma were analyzed by next generation sequencing (NGS) with 16 targeted gene panel. Among 93 patients, 84 (90.3%) had at least 1 somatic mutation in baseline ctDNA samples (average 2.74). Five patients with KRAS or NRAS hotspot mutation in the ctDNA showed significantly worse progression-free survival (PFS) (p = 0.029). Changes in average variant allele frequency (VAF) in ctDNA showed significant correlation with tumor size change at the time of first response evaluation (p = 0.020) and progressive disease (PD) (p = 0.042). Patients whose average VAF decreased below cutoff (< 1%) at the first evaluation showed significantly better PFS (p < 0.001), and the average VAF change further discriminated the PFS in the patients in partial response (p = 0.018). At the time of PD, 54 new mutations including KRAS and MAP2K1 emerged in ctDNA. ctDNA sequencing can provide mutation profile that could better reflect tumor mutation status and predict treatment outcome.

Characterization and Prognosis of Secondary Acute Myeloid Leukemia in an Asian Population: AML With Antecedent Hematological Disease Confers Worst Outcomes, Irrespective of Cytogenetic Risk.

BACKGROUND/AIM: Acute myeloid leukemia (AML) is a biologically heterogeneous disease that can be classified into de novo AML and secondary AML. Secondary AML can be further divided into therapy-related AML (t-AML) or AML evolving from antecedent hematological disorder (AHD-AML). This study evaluated the characteristics and prognosis of secondary AML in a homogeneous East Asian population who are often under-represented. PATIENTS AND METHODS: This was a retrospective, longitudinal cohort study of Korean AML patients over 18 years old treated between January 2000 and December 2013. A total of 437 de novo AML (80.3%), 41 t-AML (7.5%), and 66 AHD-AML (12.1%) were evaluated. RESULTS: First, we found that secondary AML constituted about 19.7% of all AML cases, and t-AML was more prevalent than AHD-AML. Second, we determined AHD-AML as a prognostic factor for inferior survival, independent of other risk factors (HR=2.137, 95%CI=1.534-2.977, p<0.001). The induction response rates correlated well with the overall survival. Furthermore, AHD-AML was associated with worst treatment outcomes and prognosis regardless of cytogenetic risk or age. Interestingly, t-AML was generally associated with better outcomes compared to AHD-AML despite the similarities in treatment schema. CONCLUSION: Secondary AML represents a broad spectrum of diseases and t-AML should be addressed separately from AHD-AML.

Extensive genomic and transcriptional diversity identified through massively parallel DNA and RNA sequencing of eighteen Korean individuals.

Massively parallel sequencing technologies have identified a broad spectrum of human genome diversity. Here we deep sequenced and correlated 18 genomes and 17 transcriptomes of unrelated Korean individuals. This has allowed us to construct a genome-wide map of common and rare variants and also identify variants formed during DNA-RNA transcription. We identified 9.56 million genomic variants, 23.2% of which appear to be previously unidentified. From transcriptome sequencing, we discovered 4,414 transcripts not previously annotated. Finally, we revealed 1,809 sites of transcriptional base modification, where the transcriptional landscape is different from the corresponding genomic sequences, and 580 sites of allele-specific expression. Our findings suggest that a considerable number of unexplored genomic variants still remain to be identified in the human genome, and that the integrated analysis of genome and transcriptome sequencing is powerful for understanding the diversity and functional aspects of human genomic variants.