Search

Abstract ARID1A, encoding a subunit of SWI/SNF chromatin remodeling complex, is widely recognized as a tumor suppressor gene in multiple tumor types including liver cancer. Previous studies have demonstrated that ARID1A deficiency can cause liver cancer metastasis, possibly due to the altered chromatin organization, however the underlying mechanisms remain poorly understood. To address the effect of Arid1a deficiency on chromatin organization, we generated chromatin interaction matrices, and exploited the conformation changes upon Arid1a depletion in hepatocytes. Our results demonstrated that Arid1a deficiency induced A/B compartment switching, topologically associated domain (TAD) remodeling, and decrease of chromatin loops. Further mechanism studies revealed that ATPase BRG1 of SWI/SNF complex could physically interact with RAD21, a structural subunit of chromatin architectural element cohesin; whereas ARID1A deficiency significantly diminished the coupled BRG1-RAD21. Interestingly, the tumor-associated genes within the switched compartments were differentially expressed depending upon Arid1a depletion or not. As a consequence of ARID1A deficiency-induced conformational alteration, the dysregulation of some genes such as PMP22 and GSC, promoted the invasion capacity of liver cancer cells. This study provides an insight into liver cancer tumorigenesis and progression related to ARID1A mutations.

Abstract Genetic heterogeneity of tumor is closely related to its clonal evolution, phenotypic diversity and treatment resistance, and such heterogeneity has only been characterized at single-cell sub-chromosomal scale in liver cancer. Here we reconstructed the single-variant resolution clonal evolution in human liver cancer based on single-cell mutational profiles. The results indicated that key genetic events occurred early during tumorigenesis, and an early metastasis followed by independent evolution was observed in primary liver tumor and intrahepatic metastatic portal vein tumor thrombus. By parallel single-cell RNA-Seq, the transcriptomic phenotype of HCC was found to be related with genetic heterogeneity. For the first time we reconstructed the single-cell and single-variant clonal evolution in human liver cancer, and dissection of both genetic and phenotypic heterogeneity will facilitate better understanding of their relationship.

Abstract Background High-throughput sequencing technology has yielded reliable and ultra-fast sequencing for DNA and RNA. For tumor cells of cancer patients, when combining the results of DNA and RNA sequencing, one can identify potential neoantigens that stimulate the immune response of the T cell. However, when the somatic mutations are abundant, it is computationally challenging to efficiently prioritize the identified neoantigen candidates according to their ability of activating the T cell immuno-response. Methods Numerous prioritization or prediction approaches have been proposed to address this issue but none of them considers the original DNA loci of the neoantigens from the perspective of 3D genome. Based on our previous discoveries, we propose to investigate the distribution of neoantigens with different immunogenicity abilities in 3D genome and propose to adopt this important information into neoantigen prediction. Results We retrospect the DNA origins of the immuno-positive and immuno-negative neoantigens in the context of 3D genome and discovered that DNA loci of the immuno-positive neoantigens and immuno-negative neoantigens have very different distribution pattern. Specifically, comparing to the background 3D genome, DNA loci of the immuno-positive neoantigens tend to locate at specific regions in the 3D genome. We thus used this information into neoantigen prediction and demonstrated the effectiveness of this approach. Conclusion We believe that the 3D genome information will help to increase the precision of neoantigen prioritization and discovery and eventually benefit precision and personalized medicine in cancer immunotherapy.

Glioblastoma multiforme (GBM) is the most common and malignant brain tumor, and almost half of the patients carrying EGFR-driven tumor with PTEN deficiency are resistant to EGFR-targeted therapy. EGFR amplification and/or mutation is reported in various epithelial tumors. This series of studies aimed to identify a potent compound against EGFR-driven tumor. We screened a chemical library containing over 600 individual compounds purified from Traditional Chinese Medicine against GBM cells with EGFR amplification and found that cinobufagin, the major active ingredient of Chansu, inhibited the proliferation of EGFR amplified GBM cells and PTEN deficiency enhanced its anti-proliferation effects. Cinobufagin also strongly inhibited the proliferation of carcinoma cell lines with wild-type or mutant EGFR expression. In contrast, the compound only weakly inhibited the proliferation of cancer cells with low or without EGFR expression. Cinobufagin blocked EGFR phosphorylation and its downstream signaling, which additionally induced apoptosis and cytotoxicity in EGFR amplified cancer cells. In vivo, cinobufagin blocked EGFR signaling, inhibited cell proliferation, and elicited apoptosis, thereby suppressing tumor growth in both subcutaneous and intracranial U87MG-EGFR xenograft mouse models and increasing the median survival of nude mice bearing intracranial U87MG-EGFR tumors. Cinobufagin is a potential therapeutic agent for treating malignant glioma and other human cancers expressing EGFR.

Background: Although significant progress has been made in understanding the mechanisms of steatosis and insulin resistance, the physiological functions of the epigenetic regulators in these processes remain largely elusive. Methods: Hepatocyte-specific Arid1a knockout mice were administrated with high-fat diet (HFD) for 12 weeks, then insulin sensitivity was assessed by glucose tolerance test (GTT) and insulin tolerance test (ITT). The metabolism-related indicators were determined by employing a variety of biological methods, including histology, real-time PCR, enzyme-linked immunosorbent assay (ELISA), Western blotting assay, Chromatin immunoprecipitation (ChIP), RNA-seq and assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq). Findings: Hepatocyte-specific Arid1a deletion significantly increases susceptibility to develop hepatic steatosis, insulin resistance and inflammation in mice fed a HFD. In vitro, Arid1a deletion in isolated hepatocytes directly leads to free fatty acid-induced lipid accumulation and insulin resistance. Mechanically, Arid1a deficiency impairs fatty acid oxidation by downregulating PPARα and altering the epigenetic landscape of some metabolism genes. Interpretation: These findings reveal that targeting Arid1a might be a promising therapeutic strategy for liver steatosis and insulin resistance. Fund: This work was supported by National Natural Science Foundation of China (81672772 and 81472621), China National Science and Technology Major Project for Prevention and Treatment of Infectious Diseases (No.2017ZX 10203207) and National Program on Key Research Project of China (grant no. 2016YFC0902701). Keywords: Arid1a, Steatosis, Insulin resistance, Fatty acid oxidation, PPARα

Abstract Background A major facilitator superfamily transporter Dehp2 was recently shown to be playing an important role in transport and biodegradation of haloacids in Paraburkholderia caribensis MBA4, and Dehp2 is phylogenetically conserved in Burkholderia sensu lato. Results We designed both Burkholderia sensu stricto-specific and Paraburkholderia-specific qPCR assays based on dehp2 and 16S rRNA, and validated the qPCR assays in 12 bacterial strains. The qPCR assays could detect single species of Burkholderia sensu stricto or Paraburkholderia with high sensitivity and discriminate them in mixtures with high specificity over a wide dynamic range of relative concentrations. At relatively lower cost compared with sequencing-based approach, the qPCR assays will facilitate discrimination of Burkholderia sensu stricto and Paraburkholderia in a large number of samples. Conclusions For the first time, we report the utilization of a haloacids transporter gene for discriminative purpose in Burkholderia sensu lato. This enables not only quick decision on proper handling of putative pathogenic samples in Burkholderia sensu stricto group but also future exploitation of relevant species in Paraburkholderia group for haloacids biodegradation purposes.

The conserved zona pellucida (ZP) domain is found in hundreds of extracellular proteins that are expressed in various organs and play a variety of roles as structural components, receptors and tumor suppressors. A liver-specific zona pellucida domain-containing protein (LZP), also named OIT3, has been shown to be mainly expressed in human and mouse hepatocytes; however, the physiological function of LZP in the liver remains unclear. Here, we show that Lzp deletion inhibited very low-density lipoprotein (VLDL) secretion, leading to hepatic TG accumulation and lower serum TG levels in mice. The apolipoprotein B (apoB) levels were significantly decreased in the liver, serum, and VLDL particles of LZP-deficient mice. In the presence of LZP, which is localized to the endoplasmic reticulum (ER) and Golgi apparatus, the ER-associated degradation (ERAD) of apoB was attenuated; in contrast, in the absence of LZP, apoB was ubiquitinated by AMFR, a known E3 ubiquitin ligase specific for apoB, and was subsequently degraded, leading to lower hepatic apoB levels and inhibited VLDL secretion. Interestingly, hepatic LZP levels were elevated in mice challenged with a high-fat diet and humans with simple hepatic steatosis, suggesting that LZP contributes to the physiological regulation of hepatic TG homeostasis. In general, our data establish an essential role for LZP in hepatic TG transportation and VLDL secretion by preventing the AMFR-mediated ubiquitination and degradation of apoB and therefore provide insight into the molecular function of LZP in hepatic lipid metabolism.

Upstream pathways regulating Brg1 stability and their role in carcinogenesis are unknown. Here they show Brg1 to be phosphorylated by CK1δ to promote its ubiquitination by SCFFBW7 (FBW7), Brg1 stabilization to promote gastric cancer metastasis, and suggest targeting Brg1 in FBW7 compromised gastric cancer.

Abstract Background With the developments of DNA sequencing technology, large amounts of sequencing data have been produced that provides unprecedented opportunities for advanced association studies between somatic mutations and cancer types/subtypes which further contributes to more accurate somatic mutation based cancer typing (SMCT). In existing SMCT methods however, the absence of high-level feature extraction is a major obstacle in improving the classification performance. Results We propose DeepCNA, an advanced convolutional neural network (CNN) based classifier, which utilizes copy number aberrations (CNAs) and HiC data, to address this issue. DeepCNA first pre-process the CNA data by clipping, zero padding and reshaping. Then, the processed data is fed into a CNN classifier, which extracts high-level features for accurate classification. Experimental results on the COSMIC CNA dataset indicate that 2D CNN with both cell lines of HiC data lead to the best performance. We further compare DeepCNA with three widely adopted classifiers, and demonstrate that DeepCNA has at least 78% improvement of performance. Conclusions This paper demonstrates the advantages and potential of the proposed DeepCNA model for processing of somatic point mutation based gene data, and proposes that its usage may be extended to other complex genotype-phenotype association studies.

Abstract Background The differentiation and maturation trajectories of fetal liver stem/progenitor cells (LSPCs) are not fully understood at single-cell resolution, and a priori knowledge of limited biomarkers could restrict trajectory tracking. Results We employed marker-free single-cell RNA-Seq to characterize comprehensive transcriptional profiles of 507 cells randomly selected from seven stages between embryonic day 11.5 and postnatal day 2.5 during mouse liver development, and also 52 Epcam-positive cholangiocytes from postnatal day 3.25 mouse livers. LSPCs in developing mouse livers were identified via marker-free transcriptomic profiling. Single-cell resolution dynamic developmental trajectories of LSPCs exhibited contiguous but discrete genetic control through transcription factors and signaling pathways. The gene expression profiles of cholangiocytes were more close to that of embryonic day 11.5 rather than other later staged LSPCs, cuing the fate decision stage of LSPCs. Our marker-free approach also allows systematic assessment and prediction of isolation biomarkers for LSPCs. Conclusions Our data provide not only a valuable resource but also novel insights into the fate decision and transcriptional control of self-renewal, differentiation and maturation of LSPCs.

Innate lymphoid cells (ILC) are thought to direct immune responses, but little is known about the development and maintenance of ILC subsets. Here the authors show that WASH maintains the pool of NKp46+ ILC3s by recruiting Arid1a to the aryl hydrocarbon receptor promoter and inducing its expression.

ARID1A, encoding a subunit of chromatin remodeling SWI/SNF complexes, has recently been considered as a new type of tumor suppressor gene for its somatic mutations frequently found in various human tumors, including hepatocellular carcinoma (HCC). However, the role and mechanism of inactivated ARID1A mutations in tumorigenesis remain unclear. To investigate the role of ARID1A inactivation in HCC pathogenesis, we generated hepatocyte-specific Arid1a knockout (Arid1aLKO) mice by crossing mice carrying loxP-flanked Arid1a exon 8 alleles (Arid1af/f) with albumin promoter-Cre transgenic mice. Significantly, the hepatocyte-specific Arid1a deficiency results in mouse steatohepatitis and HCC development. In Arid1aLKO mice, we found that innate immune cells, including F4/80+ macrophages and CD11c+ neutrophil cells, infiltrate into the liver parenchyma, accompanied by the increased tumor necrosis factor (TNF)-α and interleukin (IL)-6, and activation of STAT3 and NF-κB pathways. In conclusion, hepatocyte-specific Arid1a deficiency could lead to mouse steatohepatitis and HCC development. This study provides an alternative mechanism by which Arid1a deficiency contributes to HCC tumorigenesis.

PURPOSE: Hepatocellular carcinoma (HCC) is the sixth most common solid tumor worldwide and the third leading cause of cancer-related death. HCC is a particularly serious threat to the Chinese population. Although many molecular alterations are known to be involved in the tumorigenesis of hepatocytes, no systemic survey has examined the somatic mutations in HCC samples from Chinese patients. Our goal was to elucidate somatic mutations in Chinese HCC patients and investigate the possible molecular mechanisms involved in tumorigenesis. EXPERIMENTAL DESIGN: A total of 110 hepatitis B virus (HBV)-positive HCC samples and 46 HBV-negative HCC samples were genotyped for hot-spot mutations in the CSF1R, CTNNB1, KRAS, BRAF, NRAS, ERBB2, MET, PIK3CA, JAK1, and SMO genes. The transcriptomes of the CTNNB1 mutation-positive HCC samples from the HBV-positive patients (CB+ HCC) were compared to adjacent non-cancerous livers, and significantly altered genes were functionally validated in vitro. RESULTS: CTNNB1 mutations accounted for the majority of the mutations detected in our study. A slightly higher mutation rate was found in the HBV-positive patients than in their negative counterparts. A distinct pattern of CTNNB1 mutation was detected in these two populations, and drastic changes at the transcriptomic level were detected in the CB+ tumors compared to adjacent non-cancerous livers. Potential tumor suppressors (FoxA3 and Onecut1) and oncogenes (MAFG and SSX1) were functionally validated. CONCLUSIONS: Our work is the first systemic characterization of oncogenic mutations in HCC samples from Chinese patients. Targeting the Wnt-β-catenin pathway may represent a valid treatment option for Chinese HCC patients. Our work also suggests that targeting ONECUT1, FOXA3, SSX1, and MAFG may be a valid treatment option for CTNNB1 mutation positive HCC patients.

Poor prognosis of hepatocellular carcinoma (HCC) associated with late diagnosis necessitates the development of early diagnostic biomarkers. We have previously delineated the landscape of DNA methylation in HCC patients unraveling the importance of promoter hypomethylation in activation of cancer- and metastasis-driving genes. The purpose of the present study was to test the feasibility that genes that are hypomethylated in HCC could serve as candidate diagnostic markers. We use high resolution melting analysis (HRM) as a simple translatable PCR-based method to define methylation states in clinical samples. We tested seven regions selected from the shortlist of genes hypomethylated in HCC and showed that HRM analysis of several of them distinguishes methylation states in liver cancer specimens from normal adjacent liver and chronic hepatitis in the Shanghai area. Such regions were identified within promoters of neuronal membrane glycoprotein M6-B (GPM6B) and melanoma antigen family A12 (MAGEA12) genes. Differences in HRM in the immunoglobulin superfamily Fc receptor (FCRL1) separated invasive tumors from less invasive HCC. The identified biomarkers differentiated HCC from chronic hepatitis in another set of samples from Dhaka. Although the main thrust in DNA methylation diagnostics in cancer is on hypermethylated genes, our study for the first time illustrates the potential use of hypomethylated genes as markers for solid tumors. After further validation in a larger cohort, the identified DNA hypomethylated regions can become important candidate biomarkers for liver cancer diagnosis and prognosis, especially in populations with high risk for HCC development.

Both hepatitis B virus (HBV) and aflatoxin B1 (AFB1) exposure can cause liver damage as well as increase the probability of hepatocellular carcinoma (HCC). To investigate the underlying genetic changes that may influence development of HCC associated with HBV infection and AFB1 exposure, HCC patients were subdivided into 4 groups depending upon HBV and AFB1 exposure status: (HBV(+)/AFB1(+), HBV(+)/AFB1(-), HBV(-)/AFB1(+), HBV(-)/AFB1(-)). Genetic abnormalities and protein expression profiles were analyzed by array-based comparative genomic hybridization and isobaric tagging for quantitation. A total of 573 chromosomal aberrations (CNAs) including 184 increased and 389 decreased were detected in our study population. Twenty-five recurrently altered regions (RARs; chromosomal alterations observed in ≥10 patients) in chromosomes were identified. Loss of 4q13.3-q35.2, 13q12.1-q21.2 and gain of 7q11.2-q35 were observed with a higher frequency in the HBV(+)/AFB1(+), HBV(+)/AFB1(-) and HBV(-)/AFB1(+) groups compared to the HBV(-)/AFB(-) group. Loss of 8p12-p23.2 was associated with high TNM stage tumors (P = 0.038) and was an unfavorable prognostic factor for tumor-free survival (P =0.045). A total of 133 differentially expressed proteins were identified in iTRAQ proteomics analysis, 69 (51.8%) of which mapped within identified RARs. The most common biological processes affected by HBV and AFB1 status in HCC tumorigenesis were detoxification and drug metabolism pathways, antigen processing and anti-apoptosis pathways. Expression of AKR1B10 was increased significantly in the HBV(+)/AFB1(+) and HBV(-)/AFB1(+) groups. A significant correlation between the expression of AKR1B10 mRNA and protein levels as well as AKR1B10 copy number was observered, which suggest that AKR1B10 may play a role in AFB1-related hepatocarcinogenesis. In summary, a number of genetic and gene expression alterations were found to be associated with HBV and AFB1- related HCC. The possible synergistic effects of HBV and AFB1 in hepatocarcinogenesis warrant further investigations.

The tumor suppressor p53 controls multiple cellular functions including DNA repair, cell cycle arrest and apoptosis. MDM2-mediated p53 ubiquitination affects both degradation and cytoplasmic localization of p53. Several cofactors are known to modulate MDM2-mediated p53 ubiquitination and proteasomal degradation. Here we show that IRTKS, a novel IRSp53-like protein inhibited p53-induced apoptosis and depressed its transcription activity. IRTKS bound directly to p53 and increased p53 ubiquitination and cytoplasmic localization. Further studies revealed that IRTKS interacted with MDM2 and promoted low levels of MDM2-mediated p53 ubiquitination in vitro and in vivo. In unstressed cells with low levels of MDM2, IRTKS was found to stabilize the interaction of p53 and MDM2. In stressed cells, IRTKS dissociated from p53, and high levels of MDM2 induced by p53 activation mediate IRTKS poly-ubiquitination and subsequent proteasomal degradation. These data suggest that IRTKS is a novel regulator of p53, modulating low level of MDM2-mediated p53 ubiquitination in unstressed cells.

BackgroundHuman schistosomiasis is one of the most prevalent and serious parasitic diseases worldwide. Schistosoma japonicum is one of important pathogens of this disease. MicroRNAs (miRNAs) are a large group of non-coding RNAs that play important roles in regulating gene expression and protein translation in animals. Genome-wide identification of miRNAs in a given organism is a critical step to facilitating our understanding of genome organization, genome biology, evolution, and posttranscriptional regulation.Methodology/principal findingsWe sequenced two small RNA libraries prepared from different stages of the life cycle of S. japonicum, immature schistosomula and mature pairing adults, through a deep DNA sequencing approach, which yielded approximately 12 million high-quality short sequence reads containing a total of approximately 2 million non-redundant tags. Based on a bioinformatics pipeline, we identified 176 new S. japonicum miRNAs, of which some exhibited a differential pattern of expression between the two stages. Although 21 S. japonicum miRNAs are orthologs of known miRNAs within the metazoans, some nucleotides at many positions of Schistosoma miRNAs, such as miR-8, let-7, miR-10, miR-31, miR-92, miR-124, and miR-125, are indeed significantly distinct from other bilaterian orthologs. In addition, both miR-71 and some miR-2 family members in tandem are found to be clustered in a reversal direction model on two genomic loci, and two pairs of novel S. japonicum miRNAs were derived from sense and antisense DNA strands at the same genomic loci.Conclusions/significanceThe collection of S. japonicum miRNAs could be used as a new platform to study the genomic structure, gene regulation and networks, evolutionary processes, development, and host-parasite interactions. Some S. japonicum miRNAs and their clusters could represent the ancestral forms of the conserved orthologues and a model for the genesis of novel miRNAs.

Schistosomiasis remains a serious public health problem with an estimated 200 million people infected in 76 countries. Here we isolated ~ 8,400 potential protein-encoding cDNA contigs from Schistosoma japonicum after sequencing circa 84,000 expressed sequence tags. In tandem, we undertook a high-throughput proteomics approach to characterize the protein expression profiles of a number of developmental stages (cercariae, hepatic schistosomula, female and male adults, eggs, and miracidia) and tissues at the host-parasite interface (eggshell and tegument) by interrogating the protein database deduced from the contigs. Comparative analysis of these transcriptomic and proteomic data, the latter including 3,260 proteins with putative identities, revealed differential expression of genes among the various developmental stages and sexes of S. japonicum and localization of putative secretory and membrane antigens, enzymes, and other gene products on the adult tegument and eggshell, many of which displayed genetic polymorphisms. Numerous S. japonicum genes exhibited high levels of identity with those of their mammalian hosts, whereas many others appeared to be conserved only across the genus Schistosoma or Phylum Platyhelminthes. These findings are expected to provide new insights into the pathophysiology of schistosomiasis and for the development of improved interventions for disease control and will facilitate a more fundamental understanding of schistosome biology, evolution, and the host-parasite interplay.

The human genome sequence defines our inherent biological potential; the realization of the biology encoded therein requires knowledge of the function of each gene. Currently, our knowledge in this area is still limited. Several lines of investigation have been used to elucidate the structure and function of the genes in the human genome. Even so, gene prediction remains a difficult task, as the varieties of transcripts of a gene may vary to a great extent. We thus performed an exhaustive integrative characterization of 41,118 full-length cDNAs that capture the gene transcripts as complete functional cassettes, providing an unequivocal report of structural and functional diversity at the gene level. Our international collaboration has validated 21,037 human gene candidates by analysis of high-quality full-length cDNA clones through curation using unified criteria. This led to the identification of 5,155 new gene candidates. It also manifested the most reliable way to control the quality of the cDNA clones. We have developed a human gene database, called the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/). It provides the following: integrative annotation of human genes, description of gene structures, details of novel alternative splicing isoforms, non-protein-coding RNAs, functional domains, subcellular localizations, metabolic pathways, predictions of protein three-dimensional structure, mapping of known single nucleotide polymorphisms (SNPs), identification of polymorphic microsatellite repeats within human genes, and comparative results with mouse full-length cDNAs. The H-InvDB analysis has shown that up to 4% of the human genome sequence (National Center for Biotechnology Information build 34 assembly) may contain misassembled or missing regions. We found that 6.5% of the human gene candidates (1,377 loci) did not have a good protein-coding open reading frame, of which 296 loci are strong candidates for non-protein-coding RNA genes. In addition, among 72,027 uniquely mapped SNPs and insertions/deletions localized within human genes, 13,215 nonsynonymous SNPs, 315 nonsense SNPs, and 452 indels occurred in coding regions. Together with 25 polymorphic microsatellite repeats present in coding regions, they may alter protein structure, causing phenotypic effects or resulting in disease. The H-InvDB platform represents a substantial contribution to resources needed for the exploration of human biology and pathology.