Postazacitidine clone size predicts long-term outcome of patients with myelodysplastic syndromes and related myeloid neoplasms.

Azacitidine is a mainstay of therapy for myelodysplastic syndrome (MDS)-related diseases. The purpose of our study is to elucidate the effect of gene mutations on hematological response and overall survival (OS), particularly focusing on their posttreatment clone size. We enrolled a total of 449 patients with MDS or related myeloid neoplasms. They were analyzed for gene mutations in pretreatment (n = 449) and posttreatment (n = 289) bone marrow samples using targeted-capture sequencing to assess the impact of gene mutations and their posttreatment clone size on treatment outcomes. In Cox proportional hazard modeling, multihit TP53 mutation (hazard ratio [HR], 2.03; 95% confidence interval [CI], 1.42-2.91; P < .001), EZH2 mutation (HR, 1.71; 95% CI, 1.14-2.54; P = .009), and DDX41 mutation (HR, 0.33; 95% CI, 0.17-0.62; P < .001), together with age, high-risk karyotypes, low platelets, and high blast counts, independently predicted OS. Posttreatment clone size accounting for all drivers significantly correlated with International Working Group (IWG) response (P < .001, using trend test), except for that of DDX41-mutated clones, which did not predict IWG response. Combined, IWG response and posttreatment clone size further improved the prediction of the original model and even that of a recently proposed molecular prediction model, the molecular International Prognostic Scoring System (IPSS-M; c-index, 0.653 vs 0.688; P < .001, using likelihood ratio test). In conclusion, evaluation of posttreatment clone size, together with the pretreatment mutational profile as well as the IWG response play a role in better prognostication of azacitidine-treated patients with myelodysplasia.

Amplified EPOR/JAK2 Genes Define a Unique Subtype of Acute Erythroid Leukemia.

Acute erythroid leukemia (AEL) is a unique subtype of acute myeloid leukemia characterized by prominent erythroid proliferation whose molecular basis is poorly understood. To elucidate the underlying mechanism of erythroid proliferation, we analyzed 121 AEL using whole-genome, whole-exome, and/or targeted-capture sequencing, together with transcriptome analysis of 21 AEL samples. Combining publicly available sequencing data, we found a high frequency of gains and amplifications involving EPOR/JAK2 in TP53-mutated cases, particularly those having >80% erythroblasts designated as pure erythroid leukemia (10/13). These cases were frequently accompanied by gains and amplifications of ERG/ETS2 and associated with a very poor prognosis, even compared with other TP53-mutated AEL. In addition to activation of the STAT5 pathway, a common feature across all AEL cases, these AEL cases exhibited enhanced cell proliferation and heme metabolism and often showed high sensitivity to ruxolitinib in vitro and in xenograft models, highlighting a potential role of JAK2 inhibition in therapeutics of AEL. SIGNIFICANCE: This study reveals the major role of gains, amplifications, and mutations of EPOR and JAK2 in the pathogenesis of pure erythroleukemia. Their frequent response to ruxolitinib in patient-derived xenograft and cell culture models highlights a possible therapeutic role of JAK2 inhibition for erythroleukemia with EPOR/JAK2-involving lesions. This article is highlighted in the In This Issue feature, p. 369.

Clonal evolution and clinical implications of genetic abnormalities in blastic transformation of chronic myeloid leukaemia.

Blast crisis (BC) predicts dismal outcomes in patients with chronic myeloid leukaemia (CML). Although additional genetic alterations play a central role in BC, the landscape and prognostic impact of these alterations remain elusive. Here, we comprehensively investigate genetic abnormalities in 136 BC and 148 chronic phase (CP) samples obtained from 216 CML patients using exome and targeted sequencing. One or more genetic abnormalities are found in 126 (92.6%) out of the 136 BC patients, including the RUNX1-ETS2 fusion and NBEAL2 mutations. The number of genetic alterations increase during the transition from CP to BC, which is markedly suppressed by tyrosine kinase inhibitors (TKIs). The lineage of the BC and prior use of TKIs correlate with distinct molecular profiles. Notably, genetic alterations, rather than clinical variables, contribute to a better prediction of BC prognosis. In conclusion, genetic abnormalities can help predict clinical outcomes and can guide clinical decisions in CML.

Discovery of thieno[2,3-d]pyrimidin-4(3H)-one derivatives as a new class of ROCK inhibitors.

Herein, we report the discovery of a series of thieno[2,3-d]pyrimidin-4(3H)-one derivatives as a new class of ROCK inhibitors. Structure-activity relationship studies of these compounds led to the identification of the most potent compound, 3-(3-methoxybenzyl)-6-(1H-pyrrolo[2,3-b]pyridin-4-yl)thieno[2,3-d]pyrimidin-4(3H)-one (8k), which showed IC50 values of 0.004 µM and 0.001 µM against ROCK Ⅰ and ROCK Ⅱ, respectively. In vitro, 8k significantly reduced the phosphorylation level of ROCK downstream signaling protein and induce changes in cell morphology and migration. Overall, this study provides a promising lead compound for drug discovery targeting ROCKs.

Discovery of 4H-Chromen-4-one Derivatives as a New Class of Selective Rho Kinase (ROCK) Inhibitors, which Showed Potent Activity in ex Vivo Diabetic Retinopathy Models.

Diabetic retinopathy (DR) is a major cause of blindness, and there is a lack of effective treatment at present. Rho-associated coiled-coil containing serine/threonine protein kinases (ROCKs) have recently been suggested as potential targets for the DR treatment. We herein report the discovery of 4H-chromen-4-one derivatives as a new class of ROCK inhibitors. Structure-activity relationship analyses led to the identification of the most active compound, 4-(dimethylamino)-N-(3-{2-[(4-oxo-4H-chromen-7-yl)oxy]acetamido}phenyl) (12j). This compound showed excellent kinase selectivity for ROCK I and ROCK II against 387 other kinases. In retinal explants, compound 12j protected retinal neurons from high glucose-induced oxidative stress and apoptosis-mediated cell death. Furthermore, 12j administration suppressed the improper proliferation of Müller cells and promoted the regression of vascular vessels in retinal explants cultured in a high glucose microenvironment. Collectively, our data suggest that 12j could be a potential lead compound for the treatment of DR, hence deserving further in-depth studies.

Molecular heterogeneity in peripheral T-cell lymphoma, not otherwise specified revealed by comprehensive genetic profiling.

Peripheral T-cell lymphoma, not otherwise specified (PTCL, NOS) is a diagnosis of exclusion, being the most common entity in mature T-cell neoplasms, and its molecular pathogenesis remains significantly understudied. Here, combining whole-exome and targeted-capture sequencing, gene-expression profiling, and immunohistochemical analysis of tumor samples from 133 cases, we have delineated the entire landscape of somatic alterations, and discovered frequently affected driver pathways in PTCL, NOS, with and without a T-follicular helper (TFH) cell phenotype. In addition to previously reported mutational targets, we identified a number of novel recurrently altered genes, such as KMT2C, SETD1B, YTHDF2, and PDCD1. We integrated these genetic drivers using hierarchical clustering and identified a previously undescribed molecular subtype characterized by TP53 and/or CDKN2A mutations and deletions in non-TFH PTCL, NOS. This subtype exhibited different prognosis and unique genetic features associated with extensive chromosomal instability, which preferentially affected molecules involved in immune escape and transcriptional regulation, such as HLA-A/B and IKZF2. Taken together, our findings provide novel insights into the molecular pathogenesis of PTCL, NOS by highlighting their genetic heterogeneity. These results should help to devise a novel molecular classification of PTCLs and to exploit a new therapeutic strategy for this group of aggressive malignancies.

Description and initial characterization of metatranscriptomic nidovirus-like genomes from the proposed new family Abyssoviridae, and from a sister group to the Coronavirinae, the proposed genus Alphaletovirus.

Transcriptomics has the potential to discover new RNA virus genomes by sequencing total intracellular RNA pools. In this study, we have searched publicly available transcriptomes for sequences similar to viruses of the Nidovirales order. We report two potential nidovirus genomes, a highly divergent 35.9 kb likely complete genome from the California sea hare Aplysia californica, which we assign to a nidovirus named Aplysia abyssovirus 1 (AAbV), and a coronavirus-like 22.3 kb partial genome from the ornamented pygmy frog Microhyla fissipes, which we assign to a nidovirus named Microhyla alphaletovirus 1 (MLeV). AAbV was shown to encode a functional main proteinase, and a translational readthrough signal. Phylogenetic analysis suggested that AAbV represents a new family, proposed here as Abyssoviridae. MLeV represents a sister group to the other known coronaviruses. The importance of MLeV and AAbV for understanding nidovirus evolution, and the origin of terrestrial nidoviruses are discussed.

De novo transcriptome assembly for the lung of the ornamented pygmy frog (Microhyla fissipes).

Microhyla fissipes, belonging to Neobatrachia, is a new model organism to study developmental biology, adaptive mechanisms from aquatic to terrestrial life, environmental toxicology, and human disease. M. fissipes use of lungs soon after hatching makes it extremely valuable for the study of lung function and development mechanisms. However, our knowledge of genes and pathways associated with lung development in M. fissipes is very limited. In this study, we conducted de novo transcriptome assembly for the lung of M. fissipes using the Illumina HiSeq4000 platform. We obtained approximately 9.0 GB clean data from the lung of the stage 28 tadpole with lung inflation (NCBI accession numbers: SRP107055). De novo transcriptome assembly identified 209,358 transcripts and 93,813 unigenes. In addition, BLASTX against NR, NT, KO, SwissProt, PFAM, GO and KOG databases were used to annotate all the 93,813 unigenes. This study provides the transcriptome and functional annotation of genes in M. fissipes lung development, which will be useful for comparative transcriptome analyses and promote research into mechanism of lung development in anuran.

Significant association of RNF213 p.R4810K, a moyamoya susceptibility variant, with coronary artery disease.

BACKGROUND: The genetic architecture of coronary artery disease has not been fully elucidated, especially in Asian countries. Moyamoya disease is a progressive cerebrovascular disease that is reported to be complicated by coronary artery disease. Because most Japanese patients with moyamoya disease carry the p.R4810K variant of the ring finger 213 gene (RNF213), this may also be a risk factor for coronary artery disease; however, this possibility has never been tested. METHODS AND RESULTS: We genotyped the RNF213 p.R4810K variant in 956 coronary artery disease patients and 716 controls and tested the association between p.R4810K and coronary artery disease. We also validated the association in an independent population of 311 coronary artery disease patients and 494 controls. In the replication study, the p.R4810K genotypes were imputed from genome-wide genotyping data based on the 1000 Genomes Project. We used multivariate logistic regression analyses to adjust for well-known risk factors such as dyslipidemia and smoking habits. In the primary study population, the frequency of the minor variant allele was significantly higher in patients with coronary artery disease than in controls (2.04% vs. 0.98%), with an odds ratio of 2.11 (p = 0.017). Under a dominant model, after adjustment for risk factors, the association remained significant, with an odds ratio of 2.90 (95% confidence interval: 1.37-6.61; p = 0.005). In the replication study, the association was significant after adjustment for age and sex (odds ratio = 4.99; 95% confidence interval: 1.16-21.53; p = 0.031), although it did not reach statistical significance when further adjusted for risk factors (odds ratio = 3.82; 95% confidence interval: 0.87-16.77; p = 0.076). CONCLUSIONS: The RNF213 p.R4810K variant appears to be significantly associated with coronary artery disease in the Japanese population.

Research proceedings on amphibian model organisms.

Model organisms have long been important in biology and medicine due to their specific characteristics. Amphibians, especially Xenopus, play key roles in answering fundamental questions on developmental biology, regeneration, genetics, and toxicology due to their large and abundant eggs, as well as their versatile embryos, which can be readily manipulated and developed in vivo. Furthermore, amphibians have also proven to be of considerable benefit in human disease research due to their conserved cellular developmental and genomic organization. This review gives a brief introduction on the progress and limitations of these animal models in biology and human disease research, and discusses the potential and challenge of Microhyla fissipes as a new model organism.

Transcriptome profiles of metamorphosis in the ornamented pygmy frog Microhyla fissipes clarify the functions of thyroid hormone receptors in metamorphosis.

Anuran metamorphosis is an excellent system in which to study postembryonic development. Studies on Xenopus (Mesobatrachia) show that thyroid hormone receptors (TRs) regulate metamorphosis in a ligand-dependent manner by coordinating the action of hundreds of genes. However, whether this mechanism is conserved among amphibians is still unknown. To understand the molecular mechanism of this universal phenomenon, we report the transcriptional profiles of the three key developmental stages in Microhyla fissipes (Neobatrachia): premetamorphosis (PM), metamorphic climax (MC) and completion of metamorphosis (CM). In total, 2,293 differentially expressed genes were identified from comparisons of transcriptomes, and these genes showed stage-specific expression patterns. Unexpectedly, we found that TRα was highly expressed in Xenopus laevis and Bufo gargarizans at premetamorphosis but showed low expression in M. fissipes. In contrast, TRß was highly expressed during metamorphosis in M. fissipes and X. laevis. This result may imply that TRß is essential for initiating metamorphosis, at least in M. fissipes. Thus, our work not only identifies genes that are likely to be involved in Neobatrachia metamorphosis but also clarifies the roles of unliganded TRα in regulating tadpole growth and timing of metamorphosis, which may be conserved in anurans, and the role of liganded TRß in launching metamorphosis.

Complete mitochondrial genomes of two relatively closed species from Gymnocypris (Cypriniformes: Cyprinidae): genome characterization and phylogenetic considerations.

The complete mitochondrial genome for Gymnocypris eckloni and Gymnocypris przewalskii ganzihonensis is 16,686 and 16,682 bp in size, respectively. Both of the mitochondrial genomes have identical genomic organization and structure including 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes and 1 control region as found in other vertebrates, with the exception of 86 bp non-coding nucleotides between tRNA(Thr) and tRNA(Pro). The phylogenetic analyses based on the concatenated nucleotide sequence of 12 PCGs on the heavy strand supported a relatively closed relationship among G. eckloni, G. p. przewalskii and G. p. ganzihonensis. The sequence divergence between G. eckloni and G. przewalskii, as well as within G. przewalskii, is relatively low compared to the average sequence divergence among other cyprinid fishes, suggesting that they are the most closely related species resulted from a relatively recent speciation event. The estimated molecular divergence time is coincided with a relatively recent speciation event occurring during the separation of the Yellow River and Qinghai Lake.

[cDNA cloning and expression analysis of MSTN gene from Schizopygopsis pylzovi].

Myostatin (MSTN) is a member of the TGF-ß superfamily that acts as a negative regulator of skeletal muscle growth. A full-length, 2 180 bp, cDNA sequence of the myostatin gene from Schizopygopisis pylzovi was cloned with RT-PCR,5'-RACE and 3'-RACE and the cDNA clone included a 1 128 bp ORF, encoding a 375 amino acid peptide. Using PCR, we obtained the sequences of two introns of the MSTN gene and found that its structure in Schizopygopsis pylzovi was similar to that of other vertebrates, including three exons and two introns. Likewise, the putative MSTN peptide of Schizopygopsis pylzovi contains a conserved RXXR proteolytic cleavage domain, and 8 conserved cysteine residues in the C terminal of the protein, similar to other vertebrates. Phylogenetic analysis showed that the MSTN of Schizopygopsis pylzovi has high homology with other cyprinid fishes, but a low homology with mammals and birds. In the 9 examined tissues, the MSTN gene was highly expressed in heart, kidney, intestine and spermary, while weakly expressed in muscle, brain, fat, gill and hepatopancreas. Quantitative real-time PCR analysis showed that the expression of MSTN gene was different during embryo development, suggesting that the fish MSTN may not only play roles in muscle development but also contribute to other biological functions.

Convergent, parallel and correlated evolution of trophic morphologies in the subfamily schizothoracinae from the Qinghai-Tibetan plateau.

Schizothoracine fishes distributed in the water system of the Qinghai-Tibetan plateau (QTP) and adjacent areas are characterized by being highly adaptive to the cold and hypoxic environment of the plateau, as well as by a high degree of diversity in trophic morphology due to resource polymorphisms. Although convergent and parallel evolution are prevalent in the organisms of the QTP, it remains unknown whether similar evolutionary patterns have occurred in the schizothoracine fishes. Here, we constructed for the first time a tentative molecular phylogeny of the schizothoracine fishes based on the complete sequences of the cytochrome b gene. We employed this molecular phylogenetic framework to examine the evolution of trophic morphologies. We used Pagel's maximum likelihood method to estimate the evolutionary associations of trophic morphologies and food resource use. Our results showed that the molecular and published morphological phylogenies of Schizothoracinae are partially incongruent with respect to some intergeneric relationships. The phylogenetic results revealed that four character states of five trophic morphologies and of food resource use evolved at least twice during the diversification of the subfamily. State transitions are the result of evolutionary patterns including either convergence or parallelism or both. Furthermore, our analyses indicate that some characters of trophic morphologies in the Schizothoracinae have undergone correlated evolution, which are somewhat correlated with different food resource uses. Collectively, our results reveal new examples of convergent and parallel evolution in the organisms of the QTP. The adaptation to different trophic niches through the modification of trophic morphologies and feeding behaviour as found in the schizothoracine fishes may account for the formation and maintenance of the high degree of diversity and radiations in fish communities endemic to QTP.

[Construction and expression of recombinant human serum albumin-EPO fusion protein].

UNLABELLED: OBJECTIVE To construct the recombinant plasmid pCI-HLE encoding human serum album-EPO (HSA-EPO) fusion protein and to express it in CHO cell. METHODS: The cDNA encoding human serum album and EPO were amplified by PCR, and then spliced with the synsitic DNA fragment encoding GS (GGGGS), by overlap PCR extension to form LEPO. After BamH I digestion, the HSA and LEPO was ligated to generate the fusion HSA-EPO gene and was then cloned into the expression vector pCI-neo to generate the recombinant plasmid pCI-HLE. The plasmid pCI-HLE was transfected into CHO cell by liposome protocol. Then, the recombinant cells were screened by G418 and identified by PCR and Western blot. Expression of fusion protein was evaluated by Enzyme Linked Immunosorbent Assay (ELISA). RESULTS: Restrictive enzymes digestion and DNA sequencing revealed that HSA-EPO fusion gene was cloned into expression vector pCI-neo successfully. PCR and Western blot analysis confirmed that the fusion gene was integrated in the genome of CHO cells and expressed successfully. The HSA-EPO production varied from 86 Iu/(mL x 10(6) x 72 h) to 637 IU/(mLx 10(6) x 72 h). CONCLUSION: The results confirmed that HSA-EPO fusion gene can be expressed in the CHO cells, with EPO immunogenicity, which could serve as foundation for the development of long-lasting recombinant HSA-EPO protein.

Sonic hedgehog is involved in formation of the ventral optic cup by limiting Bmp4 expression to the dorsal domain.

Accumulating evidence suggests that Sonic hedgehog (Shh) signaling plays a crucial role in eye vesicle patterning in vertebrates. Shh promotes expression of Pax2 in the optic stalk and represses expression of Pax6 in the optic cup. Shh signaling contributes to establishment of both proximal-distal and dorsal-ventral axes by activating Vax1, Vax2, and Pax2. In the dorsal part of the developing retina, Bmp4 is expressed and antagonizes the ventralizing effects of Shh signaling through the activation of Tbx5 expression in chick and Xenopus. To examine the roles of Shh signaling in optic cup formation and optic stalk development, we utilized the Smoothened (Smo) conditional knockout (CKO) mouse line. Smo is a membrane protein which mediates Shh signaling into inside of cells. Cre expression was driven by Fgf15 enhancer. The ventral evagination of the optic cup deteriorated from E10 in the Smo-CKO, whereas the dorsal optic cup and optic stalk develop normally until E11. We analyzed expression of various genes such as Pax family (Pax2/Pax6), Vax family (Vax1/Vax2) and Bmp4. Bmp4 expression was greatly upregulated in the optic vesicle by the 21-somite stage. Then Vax1/2 expression was decreased at the 20- to 24-somite stages. Pax2/6 expression was affected at the 27- to 32-somite stages. Our data suggest that the effects of the absence of Shh signaling on Vax1/Vax2 are mediated through increased Bmp4 expression throughout the optic cup. Also unchanged patterns of Raldh2 and Raldh3 suggest that retinoic acid is not the downstream to Shh signaling to control the ventral optic cup morphology.

[Construction of recombinant adenovirus vector carrying human TIMP-1 cDNA and its expression in vitro].

The full-length cDNA of hTIMP-1 was obtained from a surgical patient with HCC by the method of RT-PCR. Then it was cloned into the adenoviral shuttle plasmid pAdTrack-CMV, and subsequently cotransformed into competent BJ5183 cells with the adenoviral backbone plasmid pAdEasy-1. Thereupon, a recombinant adenoviral plasmid containing full-length cDNA of hTIMP-1 was generated by homologous recombination in E. coli. The adenoviruses (AdhTIMP-1) were packaged and amplified in adenoviral packaging cells HEK 293. Then the viral titer was checked by green fluorescent protein (GFP), and the expression of hTIMP-1 was detected by the techniques of Western blot and RT-PCR. The recombinant adenovirus vector carrying human TIMP-1 was successfully constructed and expressed in vitro and may pave the way for further application in liver gene therapy.

Effects of multidrug resistance, antisense RNA on the chemosensitivity of hepatocellular carcinoma cells.

BACKGROUND: Multidrug resistance is a major obstacle in cancer chemotherapy. We examined whether the antisense RNA of multidrug resistance gene 1 (mdr1) could reverse multidrug resistance in the human hepatocellular carcinoma (HCC) cell line SMMC7721/ADM. METHODS: The recombinant adenoviruses pAdEasy-GFP-ASmdr1 product was produced by the adenoviral vector AdEasy system, which can express antisense RNA against the mdr1 gene. Following that, the recombinant adenovirus was transfected into the P-glycoprotein-producing multidrug resistance cell line, SMMC7721/ADM human HCC cells resistant to adriamycin (ADM) and daunorubicin (DNR). In order to investigate the reversal of multidrug resistance phenotype, we measured the expression of mdr1 mRNA by RT-PCR and the production of P-glycoprotein by flow cytometry. The sensitivities for ADM and DNR SMMC7721/ADM cells were examined by [3-(4, 5-dimethylthi-azol-2-yl)-2,5 diphenyl-terazolium bromide] (MTT) analysis. RESULTS: The low-level expression of mdr1 mRNA and P-glycoprotein production were observed in parental sensitive cells SMMC/7721 in addition to the overexpression of mdr1 mRNA and P-glycoprotein in SMMC7721/ADM cells. The transfection of antisense-RNA into SMMC7721/ADM cells resulted in decreases of mdr1 mRNA and P-glycoprotein, but increase of drug sensitivities. The sensitivities of transfected SMMC7721/ADM cells to ADM and DNR in IC50 reduced by 31.25% and 62.96% respectively. CONCLUSIONS: Mdr1 antisense RNA can increase the sensitivities of SMMC7721/ADM cells to anticancer drug by decreasing the expression of the mdr1 gene and inhibiting P-glycoprotein expression. This strategy may be applicable to cancer patients with P-glycoprotein mediated multidrug resistance.

Overexpression of TIMP-1 mediated by recombinant adenovirus in hepatocellular carcinoma cells inhibits proliferation and invasion in vitro.

BACKGROUND: Matrix metalloproteinases (MMPs) and its natural tissue inhibitors of metalloproteinases (TIMPs) are involved in cancer progression. This study was undertaken to determine the effects of overexpression of TIMP-1 on human hepatocellular carcinoma (HCC) cell growth, proliferation, and invasion. METHODS: Employing the efficient AdEasy(TM) system, recombinant adenovirus AdTIMP-1 containing full-length cDNA of TIMP-1 was generated by homologous recombination and amplified in 293 cells. Then, human HCC cell line (HepG2) underwent gene transfection to overexpress TIMP-1 (so-called HepG-T cells). The mRNA and protein expressions of TIMP-1 were detected with RT-PCR and Western blotting, respectively. The ultrastructure was observed with a transmission electron microscope and the proliferation of HepG-T cells was determined by MTT assay and growth curve. The potential of in vitro invasion was measured with Millicell Chamber. RESULTS: The resulting AdTIMP-1 and HepG-T cells were generated and the expression of TIMP-1 was detected in vitro. The cell proliferation curves and MTT assay showed HepG-T cells' growth, and proliferation were obviously inhibited. The invasion across Matrigel-coated filters was significantly decreased compared with controls. The suppression rate of HepG-2 cells with AdhTIMP-1 transfection was 50%, and AdhTIMP-1 transfection inhibited by more than 91.6% of the invasion into the Matrigel-coated filter (P<0.01). CONCLUSIONS: TIMP-1 overexpression results in the suppression of proliferative and invasive potential of HepG2 cells in vitro. This study demonstrates the potential role of TIMP-1 as a target for liver cancer gene therapy and has laid a foundation for further study on its anticancer function.