Genomic and transcriptomic profiles associated with response to eribulin and nivolumab combination in HER-2-negative metastatic breast cancer.

BACKGROUND: Biomarkers for predicting response to the immunotherapy and chemotherapy combination in breast cancer patients are not established. In this study, we report exploratory genomic and transcriptomic analyses of pretreatment tumor tissues from patients enrolled in phase II clinical trial of a combination of eribulin and nivolumab for HER-2-negative metastatic breast cancer (MBC) (KORNELIA trial, NCT04061863). METHODS: We analyzed associations between tumor molecular profiles based on genomic (n = 76) and transcriptomic data (n = 58) and therapeutic efficacy. Patients who achieved progression-free survival (PFS) ≥ 6 months were defined as PFS6-responders and PFS6-nonresponders otherwise. FINDINGS: Analyses on tumor mutation burden (TMB) showed a tendency toward a favorable effect on efficacy, while several analyses related to homologous recombination deficiency (HRD) indicated a potentially negative impact on efficacy. Patients harboring TP53 mutations showed significantly poor PFS6 rate and PFS, which correlated with the enrichment of cell cycle-related signatures in PFS6-nonresponders. High antigen presentation gene set enrichment scores (≥ median) were significantly associated with longer PFS. Naïve B-cell and plasma cell proportions were considerably higher in long responders (≥ 18 months). INTERPRETATION: Genomic features including TMB, HRD, and TP53 mutations and transcriptomic features related to immune cell profiles and cell cycle may distinguish responders. Our findings provide insights for further exploring the combination regimen and its biomarkers in these tumors.

Comparative RNA-Seq Analysis Revealed Tissue-Specific Splicing Variations during the Generation of the PDX Model.

Tissue-specific gene expression generates fundamental differences in the function of each tissue and affects the characteristics of the tumors that are created as a result. However, it is unclear how much the tissue specificity is conserved during grafting of the primary tumor into an immune-compromised mouse model. Here, we performed a comparative RNA-seq analysis of four different primary-patient derived xenograft (PDX) tumors. The analysis revealed a conserved RNA biotype distribution of primary-PDX pairs, as revealed by previous works. Interestingly, we detected significant changes in the splicing pattern of PDX, which was mainly comprised of skipped exons. This was confirmed by splicing variant-specific RT-PCR analysis. On the other hand, the correlation analysis for the tissue-specific genes indicated overall strong positive correlations between the primary and PDX tumor pairs, with the exception of gastric cancer cases, which showed an inverse correlation. These data propose a tissue-specific change in splicing events during PDX formation as a variable factor that affects primary-PDX integrity.

Genomic hypomethylation in cell-free DNA predicts responses to checkpoint blockade in lung and breast cancer.

Genomic hypomethylation has recently been identified as a determinant of therapeutic responses to immune checkpoint blockade (ICB). However, it remains unclear whether this approach can be applied to cell-free DNA (cfDNA) and whether it can address the issue of low tumor purity encountered in tissue-based methylation profiling. In this study, we developed an assay named iMethyl, designed to estimate the genomic hypomethylation status from cfDNA. This was achieved through deep targeted sequencing of young LINE-1 elements with > 400,000 reads per sample. iMethyl was applied to a total of 653 ICB samples encompassing lung cancer (cfDNA n = 167; tissue n = 137; cfDNA early during treatment n = 40), breast cancer (cfDNA n = 91; tissue n = 50; PBMC n = 50; cfDNA at progression n = 44), and ovarian cancer (tissue n = 74). iMethyl-liquid predicted ICB responses accurately regardless of the tumor purity of tissue samples. iMethyl-liquid was also able to monitor therapeutic responses early during treatment (3 or 6 weeks after initiation of ICB) and detect progressive hypomethylation accompanying tumor progression. iMethyl-tissue had better predictive power than tumor mutation burden and PD-L1 expression. In conclusion, our iMethyl-liquid method allows for reliable noninvasive prediction, early evaluation, and monitoring of clinical responses to ICB therapy.

DeepNeo: a webserver for predicting immunogenic neoantigens.

Non-self epitopes, whether originated from foreign substances or somatic mutations, trigger immune responses when presented by major histocompatibility complex (MHC) molecules and recognized by T cells. Identification of immunogenically active neoepitopes has significant implications in cancer and virus medicine. However, current methods are mostly limited to predicting physical binding of mutant peptides and MHCs. We previously developed a deep-learning based model, DeepNeo, to identify immunogenic neoepitopes by capturing the structural properties of peptide-MHC pairs with T cell reactivity. Here, we upgraded our DeepNeo model with up-to-date training data. The upgraded model (DeepNeo-v2) was improved in evaluation metrics and showed prediction score distribution that better fits known neoantigen behavior. The immunogenic neoantigen prediction can be conducted at https://deepneo.net.

MHC II immunogenicity shapes the neoepitope landscape in human tumors.

Despite advances in predicting physical peptide-major histocompatibility complex I (pMHC I) binding, it remains challenging to identify functionally immunogenic neoepitopes, especially for MHC II. By using the results of >36,000 immunogenicity assay, we developed a method to identify pMHC whose structural alignment facilitates T cell reaction. Our method predicted neoepitopes for MHC II and MHC I that were responsive to checkpoint blockade when applied to >1,200 samples of various tumor types. To investigate selection by spontaneous immunity at the single epitope level, we analyzed the frequency spectrum of >25 million mutations in >9,000 treatment-naive tumors with >100 immune phenotypes. MHC II immunogenicity specifically lowered variant frequencies in tumors under high immune pressure, particularly with high TCR clonality and MHC II expression. A similar trend was shown for MHC I neoepitopes, but only in particular tissue types. In summary, we report immune selection imposed by MHC II-restricted natural or therapeutic T cell reactivity.

Efficacy and Tolerability of 14-Day Tegoprazan- versus Rabeprazole-Based Triple Therapy for Eradication of Helicobacter pylori: A Real-World Evidence Study.

Background/Aims: Tegoprazan, a new, fast, and strong potassium-competitive acid blocker, has been approved for the treatment of gastric acid-related diseases in Korea. However, real-world clinical data regarding this drug are scarce. We aimed to compare the Helicobacter pylori eradication rates of tegoprazan- and rabeprazole-based triple therapy. Methods: We retrospectively reviewed data from patients who received first-line treatment for H. pylori infection using tegoprazan- or rabeprazole-based triple therapy for 2 weeks (50 mg tegoprazan or 20 mg rabeprazole+1,000 mg amoxicillin+500 mg clarithromycin twice daily). The primary endpoint was the eradication rate as determined by intention-to-treat analysis. Results: Of the 677 patients included in our study, 344 and 333 received tegoprazan-based and rabeprazole-based triple therapy, respectively. The eradication rate from intention-to-treat analysis was 76.7% (95% confidence interval [CI], 72.1% to 81.0%) for tegoprazan-based triple therapy and 75.4% (95% CI, 70.5% to 79.8%) for rabeprazole-based triple therapy. There was no significant difference in the eradication rates between the two groups (p>0.999). Per-protocol analysis also revealed no significant difference between the eradication rates of the two groups (tegoprazan 83.4% [95% CI, 79.0% to 87.2%] vs rabeprazole 83.5% [79.0% to 87.4%], p>0.999). Furthermore, there was no significant difference in adverse event rates between the two groups (tegoprazan, 27.6%; rabeprazole, 25.8%; p=0.604). Conclusions: The eradication rate of tegoprazan-based triple therapy was similar to that of rabeprazole-based triple therapy. Further studies on the dose-escalation effect of tegoprazan for H. pylori eradication and the efficacy of tegoprazan in regimens other than conventional triple therapy are needed.

7-day versus 14-day tegoprazan-based triple therapy to treat Helicobacter pylori infection: Real-world evidence.

BACKGROUND AND AIM: Potassium-competitive acid blockers (P-CABs) can be used to eradicate Helicobacter pylori infection. We aimed to evaluate the impact of treatment duration (7 vs 14 days) on successful H. pylori eradication with P-CAB-based triple therapy in Korea, where clarithromycin resistance rate is high. METHODS: We retrospectively reviewed the data of patients who received first-line treatment for H. pylori infection with tegoprazan-based triple therapy (50 mg tegoprazan + 1000 mg amoxicillin + 500 mg clarithromycin twice daily for 1 or 2 weeks). The primary endpoint was the eradication rate in intention-to-treat (ITT) analysis. RESULTS: Of the 948 patients included in the study, 435 and 513 received 7-day and 14-day tegoprazan-based triple therapy, respectively. The eradication rate was higher in the 14-day therapy group than in the 7-day therapy group (ITT, 63.9%; 95% confidence interval [CI], 59.3-68.3%] vs 78.6% [95% CI, 74.9-81.9%], respectively, P < 0.001; per-protocol, 70.5% [95% CI, 65.8-74.8%] vs 85.1% [81.7-88.1%], respectively, P < 0.001). Overall adverse event rates did not differ between the two groups. Although six patients in the 14-day treatment group discontinued the prescribed medications due to adverse events, four of them (67%) discontinued the medication within 4 days. CONCLUSIONS: The 14-day tegoprazan-based triple therapy showed a superior eradication rate and acceptable adverse events compared with the 7-day tegoprazan-based triple therapy. A 14-day treatment regimen may be required when H. pylori infection is treated with tegoprazan-based triple therapy in regions with high clarithromycin resistance.

Use of dual-electron probes reveals the role of ferritin as an iron depot in ex vivo erythropoiesis.

In the finely regulated process of mammalian erythropoiesis, the path of the labile iron pool into mitochondria for heme production is not well understood. Existing models for erythropoiesis do not include a central role for the ubiquitous iron storage protein ferritin; one model proposes that incoming endosomal Fe3+ bound to transferrin enters the cytoplasm through an ion transporter after reduction to Fe2+ and is taken up into mitochondria through mitoferrin-1 transporter. Here, we apply a dual three-dimensional imaging and spectroscopic technique, based on scanned electron probes, to measure Fe3+ in ex vivo human hematopoietic stem cells. After seven days in culture, we observe cells displaying a highly specialized architecture with anchored clustering of mitochondria and massive accumulation of nanoparticles containing high iron concentrations localized to lysosomal storage depots, identified as ferritin. We hypothesize that lysosomal ferritin iron depots enable continued heme production after expulsion of most of the cellular machinery.

Frequency of hereditary neuropathy with liability to pressure palsies (HNPP) due to 17p11.2 deletion in a Korean newborn population.

Hereditary neuropathy with liability to pressure palsies (HNPP) is an autosomal dominant disorder mainly due to a deletion of chromosome 17p11.2 including PMP22 (PMP22 Del HNPP). The prevalence of HNPP is estimated to be 0.84 to 16 per 100,000, but could be underestimated because of the mild symptoms of HNPP. In this study, we estimated the prevalence of PMP22 Del HNPP in a Korean newborn population who underwent next-generation sequencing (NGS)-based copy number variation (CNV) analysis. Of the 11,885 newborns tested by NGS-based CNV analysis, 17p11.2 deletions were found in seven samples. The prevalence of PMP22 Del HNPP was estimated to be 58.9 per 100,000 (95% confidence interval (CI), 25.8-116.5) or 1 in 1698 (95% CI, 1/909-1/5000). Our data suggest that PMP22 Del HNPP might not be uncommon at least in the Korean population.

Transcriptome and metabolome analysis in shoot and root of Valeriana fauriei.

BACKGROUND: Valeriana fauriei is commonly used in the treatment of cardiovascular diseases in many countries. Several constituents with various pharmacological properties are present in the roots of Valeriana species. Although many researches on V. fauriei have been done since a long time, further studies in the discipline make a limit due to inadequate genomic information. Hence, Illumina HiSeq 2500 system was conducted to obtain the transcriptome data from shoot and root of V. fauriei. RESULTS: A total of 97,595 unigenes were noticed from 346,771,454 raw reads after preprocessing and assembly. Of these, 47,760 unigens were annotated with Uniprot BLAST hits and mapped to COG, GO and KEGG pathway. Also, 70,013 and 88,827 transcripts were expressed in root and shoot of V. fauriei, respectively. Among the secondary metabolite biosynthesis, terpenoid backbone and phenylpropanoid biosynthesis were large groups, where transcripts was involved. To characterize the molecular basis of terpenoid, carotenoid, and phenylpropanoid biosynthesis, the levels of transcription were determined by qRT-PCR. Also, secondary metabolites content were measured using GC/MS and HPLC analysis for that gene expression correlated with its accumulation respectively between shoot and root of V. fauriei. CONCLUSIONS: We have identified the transcriptome using Illumina HiSeq system in shoot and root of V. fauriei. Also, we have demonstrated gene expressions associated with secondary metabolism such as terpenoid, carotenoid, and phenylpropanoid.

Complete mitochondrial genome of Pacific abalone (Haliotis discus hannai) from Korea.

We have sequenced and characterized the complete mitochondrial DNA of an economically and ecologically important Pacific abalone, Haliotis discus hannai (Haliotidae, Gastropoda). The mitogenome of the Pacific abalone is 16,886 nt total length with a 39.6% G+C composition. Thirty-seven genes were identified including 13 protein-coding, 2 rRNA and 22 tRNA genes. We compared the mitogenome of the Pacific abalone to a putative relative species, H. rubra.

Small non-coding RNA profiling and the role of piRNA pathway genes in the protection of chicken primordial germ cells.

BACKGROUND: Genes, RNAs, and proteins play important roles during germline development. However, the functions of non-coding RNAs (ncRNAs) on germline development remain unclear in avian species. Recent high-throughput techniques have identified several classes of ncRNAs, including micro RNAs (miRNAs), small-interfering RNAs (siRNAs), and PIWI-interacting RNAs (piRNAs). These ncRNAs are functionally important in the genome, however, the identification and annotation of ncRNAs in a genome is challenging. The aim of this study was to identify different types of small ncRNAs particularly piRNAs, and the role of piRNA pathway genes in the protection of chicken primordial germ cells (PGCs). RESULTS: At first, we performed next-generation sequencing to identify ncRNAs in chicken PGCs, and we performed ab initio predictive analysis to identify putative piRNAs in PGCs. Then, we examined the expression of three repetitive sequence-linked piRNAs and 14 genic-transcript-linked piRNAs along with their linked genes using real-time PCR. All piRNAs and their linked genes were highly expressed in PGCs. Subsequently, we knocked down two known piRNA pathway genes of chicken, PIWI-like protein 1 (CIWI) and 2 (CILI), in PGCs using siRNAs. After knockdown of CIWI and CILI, we examined their effects on the expression of six putative piRNA-linked genes and DNA double-strand breakage in PGCs. The knockdown of CIWI and CILI upregulated chicken repetitive 1 (CR1) element and RAP2B, a member of RAS oncogene family, and increased DNA double-strand breakage in PGCs. CONCLUSIONS: Our results increase the understanding of PGC-expressed piRNAs and the role of piRNA pathway genes in the protection of germ cells.

Estimation of the genome sizes of the chigger mites Leptotrombidium pallidum and Leptotrombidium scutellare based on quantitative PCR and k-mer analysis.

BACKGROUND: Leptotrombidium pallidum and Leptotrombidium scutellare are the major vector mites for Orientia tsutsugamushi, the causative agent of scrub typhus. Before these organisms can be subjected to whole-genome sequencing, it is necessary to estimate their genome sizes to obtain basic information for establishing the strategies that should be used for genome sequencing and assembly. METHOD: The genome sizes of L. pallidum and L. scutellare were estimated by a method based on quantitative real-time PCR. In addition, a k-mer analysis of the whole-genome sequences obtained through Illumina sequencing was conducted to verify the mutual compatibility and reliability of the results. RESULTS: The genome sizes estimated using qPCR were 191 ± 7 Mb for L. pallidum and 262 ± 13 Mb for L. scutellare. The k-mer analysis-based genome lengths were estimated to be 175 Mb for L. pallidum and 286 Mb for L. scutellare. The estimates from these two independent methods were mutually complementary and within a similar range to those of other Acariform mites. CONCLUSIONS: The estimation method based on qPCR appears to be a useful alternative when the standard methods, such as flow cytometry, are impractical. The relatively small estimated genome sizes should facilitate whole-genome analysis, which could contribute to our understanding of Arachnida genome evolution and provide key information for scrub typhus prevention and mite vector competence.

Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species.

Hot pepper (Capsicum annuum), one of the oldest domesticated crops in the Americas, is the most widely grown spice crop in the world. We report whole-genome sequencing and assembly of the hot pepper (Mexican landrace of Capsicum annuum cv. CM334) at 186.6× coverage. We also report resequencing of two cultivated peppers and de novo sequencing of the wild species Capsicum chinense. The genome size of the hot pepper was approximately fourfold larger than that of its close relative tomato, and the genome showed an accumulation of Gypsy and Caulimoviridae family elements. Integrative genomic and transcriptomic analyses suggested that change in gene expression and neofunctionalization of capsaicin synthase have shaped capsaicinoid biosynthesis. We found differential molecular patterns of ripening regulators and ethylene synthesis in hot pepper and tomato. The reference genome will serve as a platform for improving the nutritional and medicinal values of Capsicum species.

A synthetic model of human beta-thalassemia erythropoiesis using CD34+ cells from healthy adult donors.

Based upon the lack of clinical samples available for research in many laboratories worldwide, a significant gap exists between basic and clinical studies of beta-thalassemia major. To bridge this gap, we developed an artificially engineered model for human beta thalassemia by knocking down beta-globin gene and protein expression in cultured CD34+ cells obtained from healthy adults. Lentiviral-mediated transduction of beta-globin shRNA (beta-KD) caused imbalanced globin chain production. Beta-globin mRNA was reduced by 90% compared to controls, while alpha-globin mRNA levels were maintained. HPLC analyses revealed a 96% reduction in HbA with only a minor increase in HbF. During the terminal phases of differentiation (culture days 14-21), beta-KD cells demonstrated increased levels of insoluble alpha-globin, as well as activated caspase-3. The majority of the beta-KD cells underwent apoptosis around the polychromatophilic stage of maturation. GDF15, a marker of ineffective erythropoiesis in humans with thalassemia, was significantly increased in the culture supernatants from the beta-KD cells. Knockdown of beta-globin expression in cultured primary human erythroblasts provides a robust ex vivo model for beta-thalassemia.

LIN28B-mediated expression of fetal hemoglobin and production of fetal-like erythrocytes from adult human erythroblasts ex vivo.

Reactivation of fetal hemoglobin (HbF) holds therapeutic potential for sickle cell disease and ß-thalassemias. In human erythroid cells and hematopoietic organs, LIN28B and its targeted let-7 microRNA family, demonstrate regulated expression during the fetal-to-adult developmental transition. To explore the effects of LIN28B in human erythroid cell development, lentiviral transduction was used to knockdown LIN28B expression in erythroblasts cultured from human umbilical cord CD34+ cells. The subsequent reduction in LIN28B expression caused increased expression of let-7 and significantly reduced HbF expression. Conversely, LIN28B overexpression in cultured adult erythroblasts reduced the expression of let-7 and significantly increased HbF expression. Cellular maturation was maintained including enucleation. LIN28B expression in adult erythroblasts increased the expression of γ-globin, and the HbF content of the cells rose to levels >30% of their hemoglobin. Expression of carbonic anhydrase I, glucosaminyl (N-acetyl) transferase 2, and miR-96 (three additional genes marking the transition from fetal-to-adult erythropoiesis) were reduced by LIN28B expression. The transcription factor BCL11A, a well-characterized repressor of γ-globin expression, was significantly down-regulated. Independent of LIN28B, experimental suppression of let-7 also reduced BCL11A expression and significantly increased HbF expression. LIN28B expression regulates HbF levels and causes adult human erythroblasts to differentiate with a more fetal-like phenotype.

Let-7 microRNAs are developmentally regulated in circulating human erythroid cells.

BACKGROUND: MicroRNAs are approximately 22nt-long small non-coding RNAs that negatively regulate protein expression through mRNA degradation or translational repression in eukaryotic cells. Based upon their importance in regulating development and terminal differentiation in model systems, erythrocyte microRNA profiles were examined at birth and in adults to determine if changes in their abundance coincide with the developmental phenomenon of hemoglobin switching. METHODS: Expression profiling of microRNA was performed using total RNA from four adult peripheral blood samples compared to four cord blood samples after depletion of plasma, platelets, and nucleated cells. Labeled RNAs were hybridized to custom spotted arrays containing 474 human microRNA species (miRBase release 9.1). Total RNA from Epstein-Barr virus (EBV)-transformed lymphoblastoid cell lines provided a hybridization reference for all samples to generate microRNA abundance profile for each sample. RESULTS: Among 206 detected miRNAs, 79% of the microRNAs were present at equivalent levels in both cord and adult cells. By comparison, 37 microRNAs were up-regulated and 4 microRNAs were down-regulated in adult erythroid cells (fold change > 2; p < 0.01). Among the up-regulated subset, the let-7 miRNA family consistently demonstrated increased abundance in the adult samples by array-based analyses that were confirmed by quantitative PCR (4.5 to 18.4 fold increases in 6 of 8 let-7 miRNA). Profiling studies of messenger RNA (mRNA) in these cells additionally demonstrated down-regulation of ten let-7 target genes in the adult cells. CONCLUSION: These data suggest that a consistent pattern of up-regulation among let-7 miRNA in circulating erythroid cells occurs in association with hemoglobin switching during the fetal-to-adult developmental transition in humans.

Cytokine-mediated increases in fetal hemoglobin are associated with globin gene histone modification and transcription factor reprogramming.

Therapeutic regulation of globin genes is a primary goal of translational research aimed toward hemoglobinopathies. Signal transduction was used to identify chromatin modifications and transcription factor expression patterns that are associated with globin gene regulation. Histone modification and transcriptome profiling were performed using adult primary CD34(+) cells cultured with cytokine combinations that produced low versus high levels of gamma-globin mRNA and fetal hemoglobin (HbF). Embryonic, fetal, and adult globin transcript and protein expression patterns were determined for comparison. Chromatin immunoprecipitation assays revealed RNA polymerase II occupancy and histone tail modifications consistent with transcriptional activation only in the high-HbF culture condition. Transcriptome profiling studies demonstrated reproducible changes in expression of nuclear transcription factors associated with high HbF. Among the 13 genes that demonstrated differential transcript levels, 8 demonstrated nuclear protein expression levels that were significantly changed by cytokine signal transduction. Five of the 8 genes are recognized regulators of erythropoiesis or globin genes (MAFF, ID2, HHEX, SOX6, and EGR1). Thus, cytokine-mediated signal transduction in adult erythroid cells causes significant changes in the pattern of globin gene and protein expression that are associated with distinct histone modifications as well as nuclear reprogramming of erythroid transcription factors.

Identification of TWSG1 as a second novel erythroid regulator of hepcidin expression in murine and human cells.

In thalassemia and other iron loading anemias, ineffective erythropoiesis and erythroid signaling molecules are thought to cause inappropriate suppression of a small peptide produced by hepatocytes named hepcidin. Previously, it was reported that the erythrokine GDF15 is expressed at very high levels in thalassemia and suppresses hepcidin expression. In this study, erythroblast expression of a second molecule named twisted gastrulation (TWSG1) was explored as a potential erythroid regulator of hepcidin. Transcriptome analyses suggest TWSG1 is produced during the earlier stages of erythropoiesis. Hepcidin suppression assays demonstrated inhibition by TWSG1 as measured by quantitative polymerase chain reaction (PCR) in dosed assays (1-1000 ng/mL TWSG1). In human cells, TWSG1 suppressed hepcidin indirectly by inhibiting the signaling effects and associated hepcidin up-regulation by bone morphogenic proteins 2 and 4 (BMP2/BMP4). In murine hepatocytes, hepcidin expression was inhibited by murine Twsg1 in the absence of additional BMP. In vivo studies of Twsg1 expression were performed in healthy and thalassemic mice. Twsg1 expression was significantly increased in the spleen, bone marrow, and liver of the thalassemic animals. These data demonstrate that twisted gastrulation protein interferes with BMP-mediated hepcidin expression and may act with GDF15 to dysregulate iron homeostasis in thalassemia syndromes.

A combined approach for the classification of G protein-coupled receptors and its application to detect GPCR splice variants.

G protein-coupled receptors (GPCRs) constitute the largest family of cell surface receptors and play a central role in cellular signaling pathways. The importance of GPCRs has led to their becoming the targets of more than 50% of prescription drugs. However, drug compounds that do not differentiate between receptor subtypes can have considerable side effects and efficacy problems. An accurate classification of GPCRs can solve the side effect problems and raise the efficacy of drugs. Here, we introduce an approach that combines a fingerprint method, statistical profiles and physicochemical properties of transmembrane (TM) domains for a highly accurate classification of the receptors. The approach allows both the recognition and classification for GPCRs at the subfamily and subtype level, and allows the identification of splice variants. We found that the approach demonstrates an overall accuracy of 97.88% for subfamily classification, and 94.57% for subtype classification.