Computational Inferences of the Functions of Alternative/Noncanonical Splice Isoforms Specific to HER2+/ER-/PR- Breast Cancers, a Chromosome 17 C-HPP Study

Identification and validation of hub genes and molecular classifications associated with chronic myeloid leukemia

Background

Chronic myeloid leukemia (CML) is a kind of malignant blood tumor, which is prone to drug resistance and relapse. This study aimed to identify novel diagnostic and therapeutic targets for CML.

Methods

Differentially expressed genes (DEGs) were obtained by differential analysis of the CML cohort in the GEO database. Weighted gene co-expression network analysis (WGCNA) was used to identify CML-related co-expressed genes. Least absolute shrinkage and selection operator (LASSO) regression analysis was used to screen hub genes and construct a risk score model based on hub genes. Consensus clustering algorithm was used for the identification of molecular subtypes. Clinical samples and in vitro experiments were used to verify the expression and biological function of hub genes.

Results

A total of 378 DEGs were identified by differential analysis. 369 CML-related genes were identified by WGCNA analysis, which were mainly enriched in metabolism-related signaling pathways. In addition, CML-related genes are mainly involved in immune regulation and anti-tumor immunity, suggesting that CML has some immunodeficiency. Immune infiltration analysis confirmed the reduced infiltration of immune killer cells such as CD8+ T cells in CML samples. 6 hub genes (LINC01268, NME8, DMXL2, CXXC5, SCD and FBN1) were identified by LASSO regression analysis. The receiver operating characteristic (ROC) curve confirmed the high diagnostic value of the hub genes in the analysis and validation cohorts, and the risk score model further improved the diagnostic accuracy. hub genes were also associated with cell proliferation, cycle, and metabolic pathway activity. Two molecular subtypes, Cluster A and Cluster B, were identified based on hub gene expression. Cluster B has a lower risk score, higher levels of CD8+ T cell and activated dendritic cell infiltration, and immune checkpoint expression, and is more sensitive to commonly used tyrosine kinase inhibitors. Finally, our clinical samples validated the expression and diagnostic efficacy of hub genes, and the knockdown of LINC01268 inhibited the proliferation of CML cells, and promoted apoptosis.

Conclusion

Through WGCNA analysis and LASSO regression analysis, our study provides a new target for CML diagnosis and treatment, and provides a basis for further CML research.

Short- and long-term effects of radiation exposure at low dose and low dose rate in normal human VH10 fibroblasts

Introduction

Experimental studies complement epidemiological data on the biological effects of low doses and dose rates of ionizing radiation and help in determining the dose and dose rate effectiveness factor.

Methods

Human VH10 skin fibroblasts exposed to 25, 50, and 100 mGy of 137Cs gamma radiation at 1.6, 8, 12 mGy/h, and at a high dose rate of 23.4 Gy/h, were analyzed for radiation-induced short- and long-term effects. Two sample cohorts, i.e., discovery (n = 30) and validation (n = 12), were subjected to RNA sequencing. The pool of the results from those six experiments with shared conditions (1.6 mGy/h; 24 h), together with an earlier time point (0 h), constituted a third cohort (n = 12).

Results

The 100 mGy-exposed cells at all abovementioned dose rates, harvested at 0/24 h and 21 days after exposure, showed no strong gene expression changes. DMXL2, involved in the regulation of the NOTCH signaling pathway, presented a consistent upregulation among both the discovery and validation cohorts, and was validated by qPCR. Gene set enrichment analysis revealed that the NOTCH pathway was upregulated in the pooled cohort (p = 0.76, normalized enrichment score (NES) = 0.86). Apart from upregulated apical junction and downregulated DNA repair, few pathways were consistently changed across exposed cohorts. Concurringly, cell viability assays, performed 1, 3, and 6 days post irradiation, and colony forming assay, seeded just after exposure, did not reveal any statistically significant early effects on cell growth or survival patterns. Tendencies of increased viability (day 6) and reduced colony size (day 21) were observed at 12 mGy/h and 23.4 Gy/min. Furthermore, no long-term changes were observed in cell growth curves generated up to 70 days after exposure.

Discussion

In conclusion, low doses of gamma radiation given at low dose rates had no strong cytotoxic effects on radioresistant VH10 cells.

Cellular, transcriptomic and isoform heterogeneity of breast cancer cell line revealed by full-length single-cell RNA sequencing

Tumor heterogeneity is generated through a combination of genetic and epigenetic mechanisms, the latter of which plays an important role in the generation of stem like cells responsible for tumor formation and metastasis. Although the development of single cell transcriptomic technologies holds promise to deconvolute this complexity, a number of these techniques have limitations including drop-out and uneven coverage, which challenge the further delineation of tumor heterogeneity. We adopted deep and full-length single-cell RNA sequencing on Fluidigm's Polaris platform to reveal the cellular, transcriptomic, and isoform heterogeneity of SUM149, a triple negative breast cancer (TNBC) cell line. We first validate the quality of the TNBC sequencing data with the sequencing data from erythroleukemia K562 cell line as control. We next scrutinized well-defined marker genes for cancer stem-like cell to identify different cell populations. We then profile the isoform expression data to investigate the heterogeneity of alternative splicing patterns. Though classified as triple-negative breast cancer, the SUM149 stem cells show heterogeneous expression of marker receptors (ER, PR, and HER2) across the cells. We identified three cell populations that express patterns of stemness: epithelial-mesenchymal transition (EMT) cancer stem cells (CSCs), mesenchymal-epithelial transition (MET) CSCs and Dual-EMT-MET CSCs. These cells also manifested a high level of heterogeneity in alternative splicing patterns. For example, CSCs have shown different expression patterns of the CD44v6 exon, as well as different levels of truncated EGFR transcripts, which may suggest different potentials for proliferation and invasion among cancer stem cells. Our study identified features of the landscape of previously underestimated cellular, transcriptomic, and isoform heterogeneity of cancer stem cells in triple-negative breast cancers.

ThreaDomEx: a unified platform for predicting continuous and discontinuous protein domains by multiple-threading and segment assembly

We develop a hierarchical pipeline, ThreaDomEx, for both continuous domain (CD) and discontinuous domain (DCD) structure predictions. Starting from a query sequence, ThreaDomEx first threads it through the PDB to identify multiple structure templates, where a profile of domain conservation score (DC-score) is derived for domain-segment assignment. To further detect DCDs that consist of separated segments along the sequence, a boundary-clustering algorithm is used to refine the DCD-linker locations. In case that the templates do not contain DCDs, a domain-segment assembly process, guided by symmetry comparison, is applied for further DCD detections. ThreaDomEx was tested a set of 1111 proteins and achieved a normalized domain overlap score of 89.3% compared to experimental data, which is significantly higher than other state-of-the-art methods. It also recalls 26.7% of DCDs with 72.7% precision on the proteins for which threading failed to detect any DCDs. The server provides facilities for users to interactively refine the domain models by adjusting DC-score threshold, deleting and adding domain linkers, and assembling domain segments, which are particularly helpful for the hard targets for which current methods have a low accuracy while human-expert knowledge and experimental insights can be used for refining models. ThreaDomEX server is available at http://zhanglab.ccmb.med.umich.edu/ThreaDomEx.

Structure and Protein Interaction-Based Gene Ontology Annotations Reveal Likely Functions of Uncharacterized Proteins on Human Chromosome 17

Understanding the function of human proteins is essential to decipher the molecular mechanisms of human diseases and phenotypes. Of the 17 470 human protein coding genes in the neXtProt 2018-01-17 database with unequivocal protein existence evidence (PE1), 1260 proteins do not have characterized functions. To reveal the function of poorly annotated human proteins, we developed a hybrid pipeline that creates protein structure prediction using I-TASSER and infers functional insights for the target protein from the functional templates recognized by COFACTOR. As a case study, the pipeline was applied to all 66 PE1 proteins with unknown or insufficiently specific function (uPE1) on human chromosome 17 as of neXtProt 2017-07-01. Benchmark testing on a control set of 100 well-characterized proteins randomly selected from the same chromosome shows high Gene Ontology (GO) term prediction accuracies of 0.69, 0.57, and 0.67 for molecular function (MF), biological process (BP), and cellular component (CC), respectively. Three pipelines of function annotations (homology detection, protein-protein interaction network inference, and structure template identification) have been exploited by COFACTOR. Detailed analyses show that structure template detection based on low-resolution protein structure prediction made the major contribution to the enhancement of the sensitivity and precision of the annotation predictions, especially for cases that do not have sequence-level homologous templates. For the chromosome 17 uPE1 proteins, the I-TASSER/COFACTOR pipeline confidently assigned MF, BP, and CC for 13, 33, and 49 proteins, respectively, with predicted functions ranging from sphingosine N-acyltransferase activity and sugar transmembrane transporter to cytoskeleton constitution. We highlight the 13 proteins with confident MF predictions; 11 of these are among the 33 proteins with confident BP predictions and 12 are among the 49 proteins with confident CC. This study demonstrates a novel computational approach to systematically annotate protein function in the human proteome and provides useful insights to guide experimental design and follow-up validation studies of these uncharacterized proteins.

Alternative Splicing as a Target for Cancer Treatment

Alternative splicing is a key mechanism determinant for gene expression in metazoan. During alternative splicing, non-coding sequences are removed to generate different mature messenger RNAs due to a combination of sequence elements and cellular factors that contribute to splicing regulation. A different combination of splicing sites, exonic or intronic sequences, mutually exclusive exons or retained introns could be selected during alternative splicing to generate different mature mRNAs that could in turn produce distinct protein products. Alternative splicing is the main source of protein diversity responsible for 90% of human gene expression, and it has recently become a hallmark for cancer with a full potential as a prognostic and therapeutic tool. Currently, more than 15,000 alternative splicing events have been associated to different aspects of cancer biology, including cell proliferation and invasion, apoptosis resistance and susceptibility to different chemotherapeutic drugs. Here, we present well established and newly discovered splicing events that occur in different cancer-related genes, their modification by several approaches and the current status of key tools developed to target alternative splicing with diagnostic and therapeutic purposes.

Advances in the Chromosome-Centric Human Proteome Project: looking to the future

INTRODUCTION

The mission of the Chromosome-Centric Human Proteome Project (C-HPP), is to map and annotate the entire predicted human protein set (~20,000 proteins) encoded by each chromosome. The initial steps of the project are focused on 'missing proteins (MPs)', which lacked documented evidence for existence at protein level. In addition to remaining 2,579 MPs, we also target those annotated proteins having unknown functions, uPE1 proteins, alternative splice isoforms and post-translational modifications. We also consider how to investigate various protein functions involved in cis-regulatory phenomena, amplicons lncRNAs and smORFs. Areas covered: We will cover the scope, historic background, progress, challenges and future prospects of C-HPP. This review also addresses the question of how we can best improve the methodological approaches, select the optimal biological samples, and recommend stringent protocols for the identification and characterization of MPs. A new strategy for functional analysis of some of those annotated proteins having unknown function will also be discussed. Expert commentary: If the project moves well by reshaping the original goals, the current working modules and team work in the proposed extended planning period, it is anticipated that a progressively more detailed draft of an accurate chromosome-based proteome map will become available with functional information.

Pathology-Driven Comprehensive Proteomic Profiling of the Prostate Cancer Tumor Microenvironment

Prostate cancer is the second most common cancer in men worldwide. Gleason grading is an important predictor of prostate cancer outcomes and is influential in determining patient treatment options. Clinical decisions based on a Gleason score of 7 are difficult as the prognosis for individuals diagnosed with Gleason 4+3 cancer is much worse than for those diagnosed with Gleason 3+4 cancer. Laser capture microdissection (LCM) is a highly precise method to isolate specific cell populations or discrete microregions from tissues. This report undertook a detailed molecular characterization of the tumor microenvironment in prostate cancer to define the proteome in the epithelial and stromal regions from tumor foci of Gleason grades 3 and 4. Tissue regions of interest were isolated from several Gleason 3+3 and Gleason 4+4 tumors using telepathology to leverage specialized pathology expertise to support LCM. Over 2,000 proteins were identified following liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis of all regions of interest. Statistical analysis revealed significant differences in protein expression (>100 proteins) between Gleason 3 and Gleason 4 regions-in both stromal and epithelial compartments. A subset of these proteins has had prior strong association with prostate cancer, thereby providing evidence for the authenticity of the approach. Finally, validation of these proteins by immunohistochemistry has been obtained using an independent cohort of prostate cancer tumor specimens.Implications: This unbiased strategy provides a strong foundation for the development of biomarker protein panels with significant diagnostic and prognostic potential. Mol Cancer Res; 15(3); 281-93. ©2017 AACR.

Annotation of Alternatively Spliced Proteins and Transcripts with Protein-Folding Algorithms and Isoform-Level Functional Networks

Tens of thousands of splice isoforms of proteins have been catalogued as predicted sequences from transcripts in humans and other species. Relatively few have been characterized biochemically or structurally. With the extensive development of protein bioinformatics, the characterization and modeling of isoform features, isoform functions, and isoform-level networks have advanced notably. Here we present applications of the I-TASSER family of algorithms for folding and functional predictions and the IsoFunc, MIsoMine, and Hisonet data resources for isoform-level analyses of network and pathway-based functional predictions and protein-protein interactions. Hopefully, predictions and insights from protein bioinformatics will stimulate many experimental validation studies.

The intricacy of biomarker complexity-the identification of a genuine proteomic biomarker is more complicated than believed

Several reasons have been put forward to explain the irreproducibility of proteomic biomarker search. However, these reasons pertain to almost every part of biomarker search across the entire analytical workflow but are entirely experimental or methodological. However, in this article we point out that there is a further cause of such irreproducibility. This is not an additional methodological or experimental cause but arises directly from the biology of protein expression. It arises from the fact that disease changes the diversity within protein families. This cause of irreproducibility has been very little studied in relation to proteomic biomarker search. Gene expression is highly variable even in healthy people. Therefore, multiple proteoforms are also to be expected when gene expression is disrupted by disease, proteoforms that may be differently altered by pathology. In consequence, it is illogical to expect that the whole protein family produces a reliably usable biomarker. It is more reasonable to expect that a specific proteoform fulfills this role. Appropriate sample pre-fractionation methods and data analyses could help to identify this version, carrying the modification or the epitope required.

Genome-Wide Functional Annotation of Human Protein-Coding Splice Variants Using Multiple Instance Learning

The vast majority of human multiexon genes undergo alternative splicing and produce a variety of splice variant transcripts and proteins, which can perform different functions. These protein-coding splice variants (PCSVs) greatly increase the functional diversity of proteins. Most functional annotation algorithms have been developed at the gene level; the lack of isoform-level gold standards is an important intellectual limitation for currently available machine learning algorithms. The accumulation of a large amount of RNA-seq data in the public domain greatly increases our ability to examine the functional annotation of genes at isoform level. In the present study, we used a multiple instance learning (MIL)-based approach for predicting the function of PCSVs. We used transcript-level expression values and gene-level functional associations from the Gene Ontology database. A support vector machine (SVM)-based 5-fold cross-validation technique was applied. Comparatively, genes with multiple PCSVs performed better than single PCSV genes, and performance also improved when more examples were available to train the models. We demonstrated our predictions using literature evidence of ADAM15, LMNA/C, and DMXL2 genes. All predictions have been implemented in a web resource called "IsoFunc", which is freely available for the global scientific community through http://guanlab.ccmb.med.umich.edu/isofunc .