ITDetect: a method to detect internal tandem duplication of FMS-like tyrosine kinase (FLT3) from next-generation sequencing data with high sensitivity and clinical application.

Internal tandem duplication (ITD) of the FMS-like tyrosine kinase (FLT3) gene is associated with poor clinical outcomes in patients with acute myeloid leukemia. Although recent methods for detecting FLT3-ITD from next-generation sequencing (NGS) data have replaced traditional ITD detection approaches such as conventional PCR or fragment analysis, their use in the clinical field is still limited and requires further information. Here, we introduce ITDetect, an efficient FLT3-ITD detection approach that uses NGS data. Our proposed method allows for more precise detection and provides more detailed information than existing in silico methods. Further, it enables FLT3-ITD detection from exome sequencing or targeted panel sequencing data, thereby improving its clinical application. We validated the performance of ITDetect using NGS-based and experimental ITD detection methods and successfully demonstrated that ITDetect provides the highest concordance with the experimental methods. The program and data underlying this study are available in a public repository.

The molecular pathogenesis of Trichilemmal carcinoma.

BACKGROUND: Trichilemmal carcinoma (TC) is an extremely rare hair follicle tumor. We aimed to explore the genetic abnormalities involved in TC to gain insight into its molecular pathogenesis. METHODS: Data from patients diagnosed with TC within a 12-year period were retrospectively reviewed. Genomic DNA isolated from a formalin-fixed paraffin-embedded (FFPE) tumor tissue block was sequenced and explored for a panel of cancer genes. RESULTS: DNA was extracted from the FFPE tissue of four patients (50% female; mean age, 51.5 years) diagnosed with TC for analysis. The tumor was located in the head and neck of three patients and in the shoulder of one patient. TP53 mutations (p.Arg213*, p.Arg249Trp, and p.Arg248Gln) were found in three patients. Fusions previously identified in melanoma were detected in two patients (TACC3-FGFR3 and ROS1-GOPC fusions). Other mutations found included NF1-truncating mutation (Arg1362*), NRAS mutation (p.Gln61Lys), TOP1 amplification, and PTEN deletion. Overall, genetic changes found in TC resemble that of other skin cancers, suggesting similar pathogenesis. All patients with TP53 mutations had aggressive clinical course, two who died (OS 93 and 36 months), and one who experienced recurrent relapse. CONCLUSIONS: We reported the genomic variations found in TC, which may give insight into the molecular pathogenesis. Overall, genetic changes found in TC resembled that of other skin cancers, suggesting similar pathogenesis. TP53 mutations was were identified in patients who had an aggressive clinical course. Genetic alterations identified may further suggest the potential treatment options of TC.

NPM1 as a potential therapeutic target for atypical teratoid/rhabdoid tumors.

BACKGROUND: Atypical teratoid/rhabdoid tumors (AT/RTs) are highly malignant brain tumors with inactivation of the SMARCB1 gene, which play a critical role in genomic transcriptional control. In this study, we analyzed the genomic and transcriptomic profiles of human AT/RTs to discover new druggable targets. METHODS: Multiplanar sequencing analyses, including whole exome sequencing (WES), single nucleotide polymorphism (SNP) arrays, array comparative genomic hybridization (aCGH), and whole transcriptome sequencing (RNA-Seq), were performed on 4 AT/RT tissues. Validation of a druggable target was conducted using AT/RT cell lines. RESULTS: WES revealed that the AT/RT genome is extremely stable except for the inactivation of SMARCB1. However, we identified 897 significantly upregulated genes and 523 significantly downregulated genes identified using RNA-Seq, indicating that the transcriptional profiles of the AT/RT tissues changed substantially. Gene set enrichment assays revealed genes related to the canonical pathways of cancers, and nucleophosmin (NPM1) was the most significantly upregulated gene in the AT/RT samples. An NPM1 inhibitor (NSC348884) effectively suppressed the viability of 7 AT/RT cell lines. Network analyses showed that genes associated with NPM1 are mainly involved in cell cycle regulation. Upon treatment with an NPM1 inhibitor, cell cycle arrest at G1 phase was observed in AT/RT cells. CONCLUSIONS: We propose that NPM1 is a novel therapeutic target for AT/RTs.

Genetic Characteristics of Polycythemia Vera and Essential Thrombocythemia in Korean Patients.

BACKGROUND: Despite recent advances in the investigation of myeloproliferative neoplasms (MPN), the impact of genetic heterogeneity on its molecular pathogenesis has not been fully elucidated. Thus, in this study, we aim to characterize the genetic complexity in Korean patients with polycythemia vera (PV) and essential thrombocythemia (ET). METHODS: We conducted association studies using 84 single-nucleotide polymorphisms (SNPs) in 229 patients (96 with PV and 133 with ET) and 170 controls. Further, whole-genome sequencing was performed in six patients (two with JAK2 V617F and four with wild-type JAK2), and putative somatic mutations were validated in a further 69 ET patients. Clinical and laboratory characteristics were also analyzed. RESULTS: Several germline SNPs and the 46 haplotype were significantly associated with PV and ET. Three somatic mutations in MPDZ, IQCH, and CALR genes were selected and validated. The frequency of the CALR mutation was 58.0% (40/69) in ET patients, who did not carry JAK2/MPL mutations. Moreover, compared with JAK2 V617F-positive patients, those with CALR mutations showed lower hemoglobin and hematocrit levels (P = 0.004 and P = 0.002, respectively), higher platelet counts (P =0.008), and a lower frequency of cytoreductive therapy (P = 0.014). CONCLUSION: This study was the first comprehensive investigation of the genetic characteristics of Korean patients with PV and ET. We found that somatic mutations and the 46 haplotype contribute to PV and ET pathogenesis in Korean patients.

Combination of acalabrutinib with lenalidomide and rituximab in relapsed/refractory aggressive B-cell non-Hodgkin lymphoma: a single-arm phase II trial.

Potential synergism between Bruton's tyrosine kinase (BTK) inhibitor and lenalidomide in treating aggressive B-cell lymphoma has been suggested. Here, the authors report a single-arm phase II clinical trial of combination of acalabrutinib, lenalidomide and rituximab (R2A) in patients with aggressive relapsed/refractory aggressive (R/R) B-cell non-Hodgkin lymphoma (NHL). The primary endpoint of this study is objective response rate (ORR), and the secondary endpoints are complete remission (CR) rate, duration of response (DoR), progression-free survival (PFS) and overall survival (OS). A total of 66 patients are enrolled mostly with diffuse large B-cell lymphoma. The ORR is 54.5% and CR rate is 31.8% meeting the primary end point. The median DoR is 12.9 months, and 1-year PFS and OS rate is 33.1% and 67.5% respectively. Adverse events (AE) are manageable with the most frequent AE being neutropenia (31.8%). Patients with MYD88 mutations, subtypes known for NF-κB activation, and high BTK expression by immunohistochemistry respond well. Overall, these results show a significant efficacy of the R2A regimen in patients with aggressive R/R B-cell NHL, with exploratory biomarkers suggesting potential associations with response. (ClinicalTrials.gov 51 identifier: NCT04094142).

Human nuclear clusterin mediates apoptosis by interacting with Bcl-XL through C-terminal coiled coil domain.

Clusterin (CLU), a glycoprotein, is involved in apoptosis, producing two alternatively spliced isoforms in various cell types. The pro-apoptotic CLU appears to be a nuclear isoform (nuclear clusterin; nCLU), and the secretory CLU (sCLU) is thought to be anti-apoptotic. The detailed molecular mechanism of nCLU as a pro-apoptotic molecule has not yet been clear. In the current study, overexpressed nCLU induced apoptosis in human kidney cells. Biochemical studies revealed that nCLU sequestered Bcl-XL via a putative BH3 motif in the C-terminal coiled coil (CC2) domain, releasing Bax, and promoted apoptosis accompanied by activation of caspase-3 and cytochrome c release. These results suggest a novel mechanism of apoptosis mediated by nCLU as a pro-apoptotic molecule.

DynaMod: dynamic functional modularity analysis.

A comprehensive analysis of enriched functional categories in differentially expressed genes is important to extract the underlying biological processes of genome-wide expression profiles. Moreover, identification of the network of significant functional modules in these dynamic processes is an interesting challenge. This study introduces DynaMod, a web-based application that identifies significant functional modules reflecting the change of modularity and differential expressions that are correlated with gene expression profiles under different conditions. DynaMod allows the inspection of a wide variety of functional modules such as the biological pathways, transcriptional factor-target gene groups, microRNA-target gene groups, protein complexes and hub networks involved in protein interactome. The statistical significance of dynamic functional modularity is scored based on Z-statistics from the average of mutual information (MI) changes of involved gene pairs under different conditions. Significantly correlated gene pairs among the functional modules are used to generate a correlated network of functional categories. In addition to these main goals, this scoring strategy supports better performance to detect significant genes in microarray analyses, as the scores of correlated genes show the superior characteristics of the significance analysis compared with those of individual genes. DynaMod also offers cross-comparison between different analysis outputs. DynaMod is freely accessible at http://piech.kaist.ac.kr/dynamod.

COFECO: composite function annotation enriched by protein complex data.

COFECO is a web-based tool for a composite annotation of protein complexes, KEGG pathways and Gene Ontology (GO) terms within a class of genes and their orthologs under study. Widely used functional enrichment tools using GO and KEGG pathways create large list of annotations that make it difficult to derive consolidated information and often include over-generalized terms. The interrelationship of annotation terms can be more clearly delineated by integrating the information of physically interacting proteins with biological pathways and GO terms. COFECO has the following advanced characteristics: (i) The composite annotation sets of correlated functions and cellular processes for a given gene set can be identified in a more comprehensive and specified way by the employment of protein complex data together with GO and KEGG pathways as annotation resources. (ii) Orthology based integrative annotations among different species complement the defective annotations in an individual genome and provide the information of evolutionary conserved correlations. (iii) A term filtering feature enables users to collect the specified annotations enriched with selected function terms. (iv) A cross-comparison of annotation results between two different datasets is possible. In addition, COFECO provides a web-based GO hierarchical viewer and KEGG pathway viewer where the enrichment results can be summarized and further explored. COFECO is freely accessible at http://piech.kaist.ac.kr/cofeco.

Clonal hematopoiesis is associated with risk of severe Covid-19.

Acquired somatic mutations in hematopoietic stem and progenitor cells (clonal hematopoiesis or CH) are associated with advanced age, increased risk of cardiovascular and malignant diseases, and decreased overall survival. These adverse sequelae may be mediated by altered inflammatory profiles observed in patients with CH. A pro-inflammatory immunologic profile is also associated with worse outcomes of certain infections, including SARS-CoV-2 and its associated disease Covid-19. Whether CH predisposes to severe Covid-19 or other infections is unknown. Among 525 individuals with Covid-19 from Memorial Sloan Kettering (MSK) and the Korean Clonal Hematopoiesis (KoCH) consortia, we show that CH is associated with severe Covid-19 outcomes (OR = 1.85, 95%=1.15-2.99, p = 0.01), in particular CH characterized by non-cancer driver mutations (OR = 2.01, 95% CI = 1.15-3.50, p = 0.01). We further explore the relationship between CH and risk of other infections in 14,211 solid tumor patients at MSK. CH is significantly associated with risk of Clostridium Difficile (HR = 2.01, 95% CI: 1.22-3.30, p = 6×10-3) and Streptococcus/Enterococcus infections (HR = 1.56, 95% CI = 1.15-2.13, p = 5×10-3). These findings suggest a relationship between CH and risk of severe infections that warrants further investigation.

FiGS: a filter-based gene selection workbench for microarray data.

BACKGROUND: The selection of genes that discriminate disease classes from microarray data is widely used for the identification of diagnostic biomarkers. Although various gene selection methods are currently available and some of them have shown excellent performance, no single method can retain the best performance for all types of microarray datasets. It is desirable to use a comparative approach to find the best gene selection result after rigorous test of different methodological strategies for a given microarray dataset. RESULTS: FiGS is a web-based workbench that automatically compares various gene selection procedures and provides the optimal gene selection result for an input microarray dataset. FiGS builds up diverse gene selection procedures by aligning different feature selection techniques and classifiers. In addition to the highly reputed techniques, FiGS diversifies the gene selection procedures by incorporating gene clustering options in the feature selection step and different data pre-processing options in classifier training step. All candidate gene selection procedures are evaluated by the .632+ bootstrap errors and listed with their classification accuracies and selected gene sets. FiGS runs on parallelized computing nodes that capacitate heavy computations. FiGS is freely accessible at http://gexp.kaist.ac.kr/figs. CONCLUSION: FiGS is an web-based application that automates an extensive search for the optimized gene selection analysis for a microarray dataset in a parallel computing environment. FiGS will provide both an efficient and comprehensive means of acquiring optimal gene sets that discriminate disease states from microarray datasets.

Enhancement of TREK1 channel surface expression by protein-protein interaction with beta-COP.

TREK1 belongs to a family of two-pore-domain K(+) (K(2P)) channels and produce background currents that regulate cell excitability. In the present study, we identified a vesicle transport protein, beta-COP, as an interacting partner by yeast two-hybrid screening of a human brain cDNA library with N-terminal region of TREK1 (TREK1-N) as bait. Several in vitro and in vivo binding assays confirmed the protein-protein interaction between TREK1 and beta-COP. We also found that beta-COP was associated with TREK1 in native condition at the PC3 cells. When RFP-beta-COP was co-transfected with GFP-TREK1 into COS-7 cells, both proteins were found localized to the plasma membrane. In addition, the channel activity and surface expression of GFP-TREK1 increased dramatically by co-transfection with RFP-beta-COP. Surface expression of the TREK1 channel was also clearly reduced with the addition of beta-COP-specific shRNA. Collectively, these data suggest that beta-COP plays a critical role in the forward transport of TREK1 channel to the plasma membrane.

MED-TMA: A clinical decision support tool for differential diagnosis of TMA with enhanced accuracy using an ensemble method.

Considering difficulties in on-site ADAMTS13 testing and the performance instability of PLASMIC score according to ethnicity, we developed a prediction tool, MED-TMA (machine learning (ML) method for differential diagnosis (DDx) of thrombotic microangiopathy (TMA)) to support clinical decision. Data from 319 patients visiting 31 hospitals in Korea clinically diagnosed with primary TMA was randomly separated by 2:1 into two groups - the development dataset (D-set, n = 212), the validation dataset (V-set, n = 107). Feature elimination was conducted to select optimal clinical predictors. We developed the model with the selected features using ML methods, verifying using V-set. After the feature elimination using 19 clinical variables, five variables were selected with high importance value. Among nine ML methods, four ML methods were chosen considering the Area Under the Curves (AUC) and the correlation between the methods using D-set. We developed MED-TMA based on an optimized ensemble model with the selected four ML methods resulting in AUC values of 0.945 and 0.924 in D-set and V-set, respectively. In addition to the binary outcome, MED-TMA was capable of providing a probability for DDx of TMA. The ensemble approach driven MED-TMA showed comparable accurate and intuitive decision support for DDx of TMA to that of the existing models based on a single ML method. We provide a web-based nomogram for the appropriate use of effective but costly therapeutics to treat TMA patients (http://hematology.snu.ac.kr/medtma/).

Clonal hematopoiesis is associated with risk of severe Covid-19.

Acquired somatic mutations in hematopoietic stem and progenitor cells (clonal hematopoiesis or CH) are associated with advanced age, increased risk of cardiovascular and malignant diseases, and decreased overall survival. 1-4 These adverse sequelae may be mediated by altered inflammatory profiles observed in patients with CH. 2,5,6 A pro-inflammatory immunologic profile is also associated with worse outcomes of certain infections, including SARS-CoV-2 and its associated disease Covid-19. 7,8 Whether CH predisposes to severe Covid-19 or other infections is unknown. Among 515 individuals with Covid-19 from Memorial Sloan Kettering (MSK) and the Korean Clonal Hematopoiesis (KoCH) consortia, we found that CH was associated with severe Covid-19 outcomes (OR=1.9, 95%=1.2-2.9, p=0.01). We further explored the relationship between CH and risk of other infections in 14,211 solid tumor patients at MSK. CH was significantly associated with risk of Clostridium Difficile (HR=2.0, 95% CI: 1.2-3.3, p=6×10 -3 ) and Streptococcus/Enterococcus infections (HR=1.5, 95% CI=1.1-2.1, p=5×10 -3 ). These findings suggest a relationship between CH and risk of severe infections that warrants further investigation.

Gpx3-dependent responses against oxidative stress in Saccharomyces cerevisiae.

The yeast Saccharomyces cerevisiae has defense mechanisms identical to higher eukaryotes. It offers the potential for genome-wide experimental approaches owing to its smaller genome size and the availability of the complete sequence. It therefore represents an ideal eukaryotic model for studying cellular redox control and oxidative stress responses. S. cerevisiae Yap1 is a well-known transcription factor that is required for H2O2-dependent stress responses. Yap1 is involved in various signaling pathways in an oxidative stress response. The Gpx3 (Orp1/PHGpx3) protein is one of the factors related to these signaling pathways. It plays the role of a transducer that transfers the hydroperoxide signal to Yap1. In this study, using extensive proteomic and bioinformatics analyses, the function of the Gpx3 protein in an adaptive response against oxidative stress was investigated in wild-type, gpx3-deletion mutant, and gpx3-deletion mutant overexpressing Gpx3 protein strains. We identified 30 proteins that are related to the Gpx3- dependent oxidative stress responses and 17 proteins that are changed in a Gpx3-dependent manner regardless of oxidative stress. As expected, H2O2-responsive Gpx3-dependent proteins include a number of antioxidants related with cell rescue and defense. In addition, they contain a variety of proteins related to energy and carbohydrate metabolism, transcription, and protein fate. Based upon the experimental results, it is suggested that Gpx3-dependent stress adaptive response includes the regulation of genes related to the capacity to detoxify oxidants and repair oxidative stress-induced damages affected by Yap1 as well as metabolism and protein fate independent from Yap1.

A scientific treatment approach for acute mast cell leukemia: using a strategy based on next-generation sequencing data.

BACKGROUND: Mast cell leukemia (MCL) is the most aggressive form of systemic mastocytosis disorders. Owing to its rarity, neither pathogenesis nor standard treatment is established for this orphan disease. Hence, we tried to treat a patient with MCL based on the exome and transcriptome sequencing results of the patient's own DNA and RNA. METHODS: First, tumor DNA and RNA were extracted from bone marrow at the time of diagnosis. Germline DNA was extracted from the patient's saliva 45 days after induction chemotherapy and used as a control. Then, we performed whole-exome sequencing (WES) using the DNA and whole transcriptome sequencing (WTS) using the RNA. Single nucleotide variants (SNVs) were called using MuTect and GATK. Samtools, FusionMap, and Gene Set Enrichment Analysis were utilized to analyze WTS results. RESULTS: WES and WTS results revealed mutation in KIT S476I. Fusion analysis was performed using WTS data, which suggested a possible RARα-B2M fusion. When RNA expression analysis was performed using WTS data, upregulation of PIK3/AKT pathway, downstream of KIT and mTOR, was observed. Based on our WES and WTS results, we first administered all-trans retinoic acid, then dasatinib, and finally, an mTOR inhibitor. CONCLUSION: We present a case of orphan disease where we used a targeted approach using WES and WTS data of the patient. Even though our treatment was not successful, use of our approach warrants further validation.

Integrated genomic analyses identify frequent gene fusion events and VHL inactivation in gastrointestinal stromal tumors.

Gastrointestinal stromal tumors (GISTs) are the most common mesenchymal tumors of the gastrointestinal tract. We sequenced nine exomes and transcriptomes, and two genomes of GISTs for integrated analyses. We detected 306 somatic variants in nine GISTs and recurrent protein-altering mutations in 29 genes. Transcriptome sequencing revealed 328 gene fusions, and the most frequently involved fusion events were associated with IGF2 fused to several partner genes including CCND1, FUS, and LASP1. We additionally identified three recurrent read-through fusion transcripts: POLA2-CDC42EP2, C8orf42-FBXO25, and STX16-NPEPL1. Notably, we found intragenic deletions in one of three exons of the VHL gene and increased mRNAs of VEGF, PDGF-ß, and IGF-1/2 in 56% of GISTs, suggesting a mechanistic link between VHL inactivation and overexpression of hypoxia-inducible factor target genes in the absence of hypoxia. We also identified copy number gain and increased mRNA expression of AMACR, CRIM1, SKP2, and CACNA1E. Mapping of copy number and gene expression results to the KEGG pathways revealed activation of the JAK-STAT pathway in small intestinal GISTs and the MAPK pathway in wild-type GISTs. These observations will allow us to determine the genetic basis of GISTs and will facilitate further investigation to develop new therapeutic options.

Genomic dynamics associated with malignant transformation in IDH1 mutated gliomas.

The genomic mechanism responsible for malignant transformation remains an open question for glioma researchers, where differing conclusions have been drawn based on diverse study conditions. Therefore, it is essential to secure direct evidence using longitudinal samples from the same patient. Moreover, malignant transformation of IDH1-mutated gliomas is of potential interest, as its genomic mechanism under influence of oncometabolite remains unclear, and even higher rate of malignant transformation was reported in IDH1-mutated low grade gliomas than in wild-type IDH1 tumors. We have analyzed genomic data using next-generation sequencing technology for longitudinal samples from 3 patients with IDH1-mutated gliomas whose disease had progressed from a low grade to a high grade phenotype. Comprehensive analysis included chromosomal aberrations as well as whole exome and transcriptome sequencing, and the candidate driver genes for malignant transformation were validated with public database. Integrated analysis of genomic dynamics in clonal evolution during the malignant transformation revealed alterations in the machinery regulating gene expression, including the spliceosome complex (U2AF2), transcription factors (TCF12), and chromatin remodelers (ARID1A). Moreover, consequential expression changes implied the activation of genes associated with the restoration of the stemness of cancer cells. The alterations in genetic regulatory mechanisms may be the key factor for the major phenotypic changes in IDH1 mutated gliomas. Despite being limited to a small number of cases, this analysis provides a direct example of the genomic changes responsible for malignant transformation in gliomas.

Identification and integrative analysis of 28 novel genes specifically expressed and developmentally regulated in murine spermatogenic cells.

Mammalian spermatogenesis is a highly ordered process that occurs in mitotic, meiotic, and postmeiotic phases. The unique mechanisms responsible for this tightly regulated developmental process suggest the presence of an intrinsic genetic program composed of spermatogenic cell-specific genes. In this study, we analyzed the mouse round spermatid UniGene library currently containing 2124 gene-oriented transcript clusters, predicting that 467 of them are testis-specific genes, and systematically identified 28 novel genes with evident testis-specific expression by in silico and in vitro approaches. We analyzed these genes by Northern blot hybridization and cDNA cloning, demonstrating the presence of additional transcript sequences in five genes and multiple transcript isoforms in six genes. Genomic analysis revealed lack of human orthologues for 10 genes, implying a relationship between these genes and male reproduction unique to mouse. We found that all of the novel genes are expressed in developmentally regulated and stage-specific patterns, suggesting that they are primary regulators of male germ cell development. Using computational bioinformatics tools, we found that 20 gene products are potentially involved in various processes during spermatogenesis or fertilization. Taken together, we predict that over 20% of the genes from the round spermatid library are testis-specific, have discovered the 28 authentic, novel genes with probable spermatogenic cell-specific expression by the integrative approach, and provide new and thorough information about the novel genes by various in vitro and in silico analyses. Thus, the study establishes on a comprehensive scale a new basis for studies to uncover molecular mechanisms underlying the reproductive process.