Predicting DNA Methylation State of CpG Dinucleotide Using Genome Topological Features and Deep Networks

The hypo- or hyper-methylation of the human genome is one of the epigenetic features of leukemia. However, experimental approaches have only determined the methylation state of a small portion of the human genome. We developed deep learning based (stacked denoising autoencoders, or SdAs) software named DeepMethyl to predict the methylation state of DNA CpG dinucleotides using features inferred from three-dimensional genome topology (based on Hi-C) and DNA sequence patterns. We used the experimental data from immortalised myelogenous leukemia (K562) and healthy lymphoblastoid (GM12878) cell lines to train the learning models and assess prediction performance. We have tested various SdA architectures with different configurations of hidden layer(s) and amount of pre-training data and compared the performance of deep networks relative to support vector machines (SVMs). Using the methylation states of sequentially neighboring regions as one of the learning features, an SdA achieved a blind test accuracy of 89.7% for GM12878 and 88.6% for K562. When the methylation states of sequentially neighboring regions are unknown, the accuracies are 84.82% for GM12878 and 72.01% for K562. We also analyzed the contribution of genome topological features inferred from Hi-C. DeepMethyl can be accessed at http://dna.cs.usm.edu/deepmethyl/

TADKB:Family Classification and a Knowledge Base of Topologically Associating Domains

Background: Topologically associating domains (TADs) are considered the structural and functional units of the genome. However, there is a lack of an integrated resource for TADs in the literature where researchers can obtain family classifications and detailed information about TADs. Results: We built an online knowledge base TADKB integrating knowledge for TADs in eleven cell types of human and mouse. For each TAD, TADKB provides the predicted three-dimensional (3D) structures of chromosomes and TADs, and detailed annotations about the protein-coding genes and long non-coding RNAs (lncRNAs) existent in each TAD. Besides the 3D chromosomal structures inferred by population Hi-C, the single-cell haplotype-resolved chromosomal 3D structures of 17 GM12878 cells are also integrated in TADKB. A user can submit query gene/lncRNA ID/sequence to search for the TAD(s) that contain(s) the query gene or lncRNA. We also classified TADs into families. To achieve that, we used the TM-scores between reconstructed 3D structures of TADs as structural similarities and the Pearson’s correlation coefficients between the fold enrichment of chromatin states as functional similarities. All of the TADs in one cell type were clustered based on structural and functional similarities respectively using the spectral clustering algorithm with various predefined numbers of clusters. We have compared the overlapping TADs from structural and functional clusters and found that most of the TADs in the functional clusters with depleted chromatin states are clustered into one or two structural clusters. This novel finding indicates a connection between the 3D structures of TADs and their DNA functions in terms of chromatin states. Conclusion: TADKB is available at href= http://dna.cs.miami.edu/TADKB/ target= _blank \u3ehttp://dna.cs.miami.edu/TADKB

Comprehensive Network Analysis Reveals Alternative Splicing-Related lncRNAs in Hepatocellular Carcinoma

© Copyright © 2020 Wang, Wang, Bhat, Chen, Xu, Mo, Yi and Zhou. It is increasingly appreciated that long non-coding RNAs (lncRNAs) associated with alternative splicing (AS) could be involved in aggressive hepatocellular carcinoma. Although many recent studies show the alteration of RNA alternative splicing by deregulated lncRNAs in cancer, the extent to which and how lncRNAs impact alternative splicing at the genome scale remains largely elusive. We analyzed RNA-seq data obtained from 369 hepatocellular carcinomas (HCCs) and 160 normal liver tissues, quantified 198,619 isoform transcripts, and identified a total of 1,375 significant AS events in liver cancer. In order to predict novel AS-associated lncRNAs, we performed an integration of co-expression, protein-protein interaction (PPI) and epigenetic interaction networks that links lncRNA modulators (such as splicing factors, transcript factors, and miRNAs) along with their targeted AS genes in HCC. We developed a random walk-based multi-graphic (RWMG) model algorithm that prioritizes functional lncRNAs with their associated AS targets to computationally model the heterogeneous networks in HCC. RWMG shows a good performance evaluated by the ROC curve based on cross-validation and bootstrapping strategies. As a conclusion, our robust network-based framework has derived 31 AS-related lncRNAs that not only validates known cancer-associated cases MALAT1 and HOXA11-AS, but also reveals new players such as DNM1P35 and DLX6-AS1with potential functional implications. Survival analysis further provides insights into the clinical significance of identified lncRNAs

Rapid exchange of mammalian topoisomerase IIα at kinetochores and chromosome arms in mitosis

Astable cell line (GT2-LPk) derived from LLC-Pk was created in which endogenous DNA topoisomerase IIα (topoIIα) protein was downregulated and replaced by the expression of topoIIα fused with enhanced green fluorescent protein (EGFP–topoIIα). The EGFP–topoIIα faithfully mimicked the distribution of the endogenous protein in both interphase and mitosis. In early stages of mitosis, EGFP–topoIIα accumulated at kinetochores and in axial lines extending along the chromosome arms. During anaphase, EGFP–topoIIα diminished at kinetochores and increased in the cytoplasm with a portion accumulating into large circular foci that were mobile and appeared to fuse with the reforming nuclei. These cytoplasmic foci appearing at anaphase were coincident with precursor organelles of the reforming nucleolus called nucleolus-derived foci (NDF). Photobleaching of EGFP–topoIIα associated with kinetochores and chromosome arms showed that the majority of the protein rapidly exchanges (t1/2 of 16 s). Catalytic activity of topoIIα was essential for rapid dynamics, as ICRF-187, an inhibitor of topoIIα, blocked recovery after photobleaching. Although some topoIIα may be stably associated with chromosomes, these studies indicate that the majority undergoes rapid dynamic exchange. Rapid mobility of topoIIα in chromosomes may be essential to resolve strain imparted during chromosome condensation and segregation

Long Non-Coding RNAs As Prognostic Markers In Human Breast Cancer

Long non-coding RNAs (lncRNAs) have been recently shown to play an important role in gene regulation and normal cellular functions, and disease processes. However, despite the overwhelming number of lncRNAs identified to date, little is known about their role in cancer for vast majority of them. The present study aims to determine whether lncRNAs can serve as prognostic markers in human breast cancer. We interrogated the breast invasive carcinoma dataset of the Cancer Genome Atlas (TCGA) at the cBioPortal consisting of ~ 1,000 cases. Among 2,730 lncRNAs analyzed, 577 lncRNAs had alterations ranging from 1% to 32% frequency, which include mutations, alterations of copy number and RNA expression. We found that deregulation of 11 lncRNAs, primarily due to copy number alteration, is associated with poor overall survival. At RNA expression level, upregulation of 4 lncRNAs (LINC00657, LINC00346, LINC00654 and HCG11) was associated with poor overall survival. A third signature consists of 9 lncRNAs (LINC00705, LINC00310, LINC00704, LINC00574, FAM74A3, UMODL1-AS1, ARRDC1-AS1, HAR1A, and LINC00323) and their upregulation can predict recurrence. Finally, we selected LINC00657 to determine their role in breast cancer, and found that LINC00657 knockout significantly suppresses tumor cell growth and proliferation, suggesting that it plays an oncogenic role. Together, these results highlight the clinical significance of lncRNAs, and thus, these lncRNAs may serve as prognostic markers for breast cancer

Establishment of an Efficient and Flexible Genetic Manipulation Platform Based on a Fosmid Library for Rapid Generation of Recombinant Pseudorabies Virus

Conventional genetic engineering of pseudorabies virus (PRV) is essentially based on homologous recombination or bacterial artificial chromosome. However, these techniques require multiple plaque purification, which is labor-intensive and time-consuming. The aim of the present study was to develop an efficient, direct, and flexible genetic manipulation platform for PRV. To this end, the PRV genomic DNA was extracted from purified PRV virions and sheared into approximately 30–45-kb DNA fragments. After end-blunting and phosphorylation, the DNA fragments were separated by pulsed-field gel electrophoresis, the recovered DNA fragments were inserted into the cloning-ready fosmids. The fosmids were then transformed into Escherichia coli and selected clones were end-sequenced for full-length genome assembly. Overlapping fosmid combinations that cover the complete genome of PRV were directly transfected into Vero cells and PRV was rescued. The morphology and one-step growth curve of the rescued virus were indistinguishable from those of the parent virus. Based on this system, a recombinant PRV expressing enhanced green fluorescent protein fused with the VP26 gene was generated within 2 weeks, and this recombinant virus can be used to observe the capsid transport in axons. The new genetic manipulation platform developed in the present study is an efficient, flexible, and stable method for the study of the PRV life cycle and development of novel vaccines

Generation of shRNAs from randomized oligonucleotides

Suppression of gene expression by small interfering RNA (siRNA) has proved to be a gene-specific and cost effective alternative to other gene suppression technologies. Short hairpin RNAs (shRNAs) generated from the vector-based expression are believed to be processed into functional siRNAs in vivo, leading to gene silencing. Since an shRNA library carries a large pool of potential siRNAs, such a library makes it possible to knock down gene expression at the genome wide scale. Although much of research has been focused on generating shRNA libraries from either individually made gene specific sequences or cDNA libraries, there is no report on constructing randomized shRNA libraries, which could provide a good alternative to these existing libraries. We have developed a method of constructing shRNAs from randomized oligonucleotides. Through this method, one can generate a partially or fully randomized shRNA library for various functional analyses. We validated this procedure by constructing a p53-specific shRNA. Western blot revealed that the p53-shRNA successfully suppressed expression of the endogenous p53 in MCF-7 cells. We then made a partially randomized shRNA library. Sequencing of 15 randomly picked cloned confirmed the randomness of the library. Therefore, the library can be used for various functional assays, such as target validation when a suitable screening or selection method is available

Adjustment of Prescription by ORAC assay to Enhance the Effect of Medicinal Liqueur

Objective: The ingredient herbs of JING liqueur were studied, and prescription was adjusted by oxygen radical absorbance capacity (ORAC) assay to enhance the effect of this medicinal liqueur. Method: Values of very herb's ORAC were measured, the proportions of high ORAC value herbs were increased, while the proportions of low ORAC value herbs were decreased. Adjusted JING liqueur was compared with primary product by the test of anti-fatigue and enhance immunity experiment. Result: Compared with primary product, the anti-fatigue effect of adjusted JING liqueur was significantly increased, and the enhance immunity effect of adjusted JING liqueur was obviously improved. Conclusion: The effect of JING liqueur could enhanced by adjust proportions of every ingredients across their ORAC, increased the proportions of high ORAC value herb could obviously improved the anti-fatigue and enhance immunity efficacy

Genome-Wide and Differential Proteomic Analysis of Hepatitis B Virus and Aflatoxin B1 Related Hepatocellular Carcinoma in Guangxi, China

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 observed, 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

N-myc downstream regulated gene 1 modulates Wnt-β-catenin signalling and pleiotropically suppresses metastasis

Wnt signalling has pivotal roles in tumour progression and metastasis; however, the exact molecular mechanism of Wnt signalling in the metastatic process is as yet poorly defined. Here we demonstrate that the tumour metastasis suppressor gene, NDRG1, interacts with the Wnt receptor, LRP6, followed by blocking of the Wnt signalling, and therefore, orchestrates a cellular network that impairs the metastatic progression of tumour cells. Importantly, restoring NDRG1 expression by a small molecule compound significantly suppressed the capability of otherwise highly metastatic tumour cells to thrive in circulation and distant organs in animal models. In addition, our analysis of clinical cohorts data indicate that Wnt+/NDRG−/LRP+ signature has a strong predictable value for recurrence-free survival of cancer patients. Collectively, we have identified NDRG1 as a novel negative master regulator of Wnt signalling during the metastatic progression, which opens an opportunity to define a potential therapeutic target for metastatic disease