Parkinson's disease is associated with clonal hematopoiesis with TET2 mutation.

Clonal hematopoiesis of indeterminate potential (CHIP), a premalignant expansion of mutated hematopoietic stem cells, is linked to immune alterations. Given the role of neuroinflammation and immune dysfunction in Parkinson's disease (PD), we hypothesized a connection between CHIP and PD. We analyzed peripheral blood DNA from 341 PD, 92 isolated REM sleep behavior disorder (iRBD) patients, and 5003 controls using targeted sequencing of 24 genes associated with hematologic neoplasms. PD cases were classified by clinical progression mode: fast, slow, and typical. Using multivariable logistic regression models, CHIP prevalence was assessed against controls with a 1.0% variant allele fraction threshold. CHIP with TET2 mutations was more prevalent in PD than controls (aOR 1.75, 95% CI 1.11-2.77, p = 0.017), particularly in the fast motor progression subgroup (aOR 3.19, p = 0.004). No distinct associations were observed with iRBD. PD is linked to increased odds of CHIP with TET2 mutations, suggesting immune dysregulation in PD pathophysiology.

Clonal hematopoiesis with DNMT3A mutations is associated with multiple system atrophy.

BACKGROUND: Clonal hematopoiesis of indeterminate potential (CHIP) is associated with cardiovascular diseases and other disorders, possibly via inflammation. Recent research suggests a connection of CHIP with neurodegenerative disorders. OBJECTIVE: We aimed to investigate the association between multiple system atrophy (MSA) and CHIP. METHODS: We included 100 patients with MSA and 4457 controls. Targeted sequencing of peripheral blood DNA samples was performed, focusing on a panel of 25 genes commonly. LINKED TO CHIP: The prevalence of CHIP in patients with MSA was assessed against controls at variant allele frequency (VAF) thresholds of 1.5 % and 2.0 %. RESULTS: DNMT3A mutation rates were significantly higher in patients with MSA, with a VAF of 1.5 %, which remained significant after adjusting for age and sex (adjusted odds ratio, 1.848; 95 % CI, 1.024-3.335; p = 0.0416). CONCLUSION: Our results suggest an association between DNMT3A mutations and MSA.

Clonal haematopoiesis of indeterminate potential and atrial fibrillation: an east Asian cohort study.

BACKGROUND AND AIMS: Both clonal haematopoiesis of indeterminate potential (CHIP) and atrial fibrillation (AF) are age-related conditions. This study investigated the potential role of CHIP in the development and progression of AF. METHODS: Deep-targeted sequencing of 24 CHIP mutations (a mean depth of coverage = 1000×) was performed in 1004 patients with AF and 3341 non-AF healthy subjects. Variant allele fraction ≥ 2.0% indicated the presence of CHIP mutations. The association between CHIP and AF was evaluated by the comparison of (i) the prevalence of CHIP mutations between AF and non-AF subjects and (ii) clinical characteristics discriminated by CHIP mutations within AF patients. Furthermore, the risk of clinical outcomes-the composite of heart failure, ischaemic stroke, or death-according to the presence of CHIP mutations in AF was investigated from the UK Biobank cohort. RESULTS: The mean age was 67.6 ± 6.9 vs. 58.5 ± 6.5 years in AF (paroxysmal, 39.0%; persistent, 61.0%) and non-AF cohorts, respectively. CHIP mutations with a variant allele fraction of ≥2.0% were found in 237 (23.6%) AF patients (DNMT3A, 13.5%; TET2, 6.6%; and ASXL1, 1.5%) and were more prevalent than non-AF subjects [356 (10.7%); P < .001] across the age. After multivariable adjustment (age, sex, smoking, body mass index, diabetes, and hypertension), CHIP mutations were 1.4-fold higher in AF [adjusted odds ratio (OR) 1.38; 95% confidence interval 1.10-1.74, P < .01]. The ORs of CHIP mutations were the highest in the long-standing persistent AF (adjusted OR 1.50; 95% confidence interval 1.14-1.99, P = .004) followed by persistent (adjusted OR 1.44) and paroxysmal (adjusted OR 1.33) AF. In gene-specific analyses, TET2 somatic mutation presented the highest association with AF (adjusted OR 1.65; 95% confidence interval 1.05-2.60, P = .030). AF patients with CHIP mutations were older and had a higher prevalence of diabetes, a longer AF duration, a higher E/E', and a more severely enlarged left atrium than those without CHIP mutations (all P < .05). In UK Biobank analysis of 21 286 AF subjects (1297 with CHIP and 19 989 without CHIP), the CHIP mutation in AF is associated with a 1.32-fold higher risk of a composite clinical event (heart failure, ischaemic stroke, or death). CONCLUSIONS: CHIP mutations, primarily DNMT3A or TET2, are more prevalent in patients with AF than non-AF subjects whilst their presence is associated with a more progressive nature of AF and unfavourable clinical outcomes.

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.

Consistent count region-copy number variation (CCR-CNV): an expandable and robust tool for clinical diagnosis of copy number variation at the exon level using next-generation sequencing data.

PURPOSE: Despite the importance of exonic copy number variations (CNVs) in human genetic diseases, reliable next-generation sequencing-based methods for detecting them are unavailable. We developed an expandable and robust exonic CNV detection tool called consistent count region (CCR)-CNV. METHODS: In total, about 1000 samples of the truth set were used for validating CCR-CNV. We compared CCR-CNV performance with 2 well-known CNV tools. Finally, to overcome the limitations of CCR-CNV, we devised a combined approach. RESULTS: The mean sensitivity and specificity of CCR-CNV alone were above 95%, which was superior to that of other CNV tools, such as DECoN and Atlas-CNV. However, low covered region and positive predictive value and high false discovery rate act as obstacles to its use in clinical settings. The combined approach showed much improved performance than CCR-CNV alone. CONCLUSION: In this study, we present a novel diagnostic tool that allows the identification of exonic CNVs with high confidence using various reagents and clinical next-generation sequencing platforms. We validated this method using the largest multiple ligation-dependent probe amplification-confirmed data set, including sufficient copy normal control data. The approach, combined with existing CNV tools, allows the implementation of CCR-CNV in clinical settings.

Noninvasive prenatal test of single-gene disorders by linked-read direct haplotyping: application in various diseases.

Direct haplotyping enables noninvasive prenatal testing (NIPT) without analyzing proband, which is a promising strategy for pregnancies at risk of an inherited single-gene disorder. Here, we aimed to expand the scope of single-gene disorders that NIPT using linked-read direct haplotyping would be applicable to. Three families at risk of myotonic dystrophy type 1, lipoid congenital adrenal hyperplasia, and Fukuyama congenital muscular dystrophy were recruited. All cases exhibited distinct characteristics that are often encountered as hurdles (i.e., repeat expansion, identical variants in both parents, and novel variants with retrotransposon insertion) in the universal clinical application of NIPT. Direct haplotyping of parental genomes was performed by linked-read sequencing, combined with allele-specific PCR, if necessary. Target DMPK, STAR, and FKTN genes in the maternal plasma DNA were sequenced. Posterior risk calculations and an Anderson-Darling test were performed to deduce the maternal and paternal inheritance, respectively. In all cases, we could predict the inheritance of maternal mutant allele with > 99.9% confidence, while paternal mutant alleles were not predicted to be inherited. Our study indicates that direct haplotyping and posterior risk calculation can be applied with subtle modifications to NIPT for the detection of an expanded range of diseases.

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/).

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.

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.

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.

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.

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.

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.