This paper proposes a nonlinear regression model to predict soft tissue deformation after maxillofacial surgery. The feature which served as input in the model is extracted with Finite Element Model (FEM). The output in the model is the facial deformation calculated from the preoperative and postoperative 3D data. After finding the relevance between feature and facial deformation by using the regression model, we establish a general relationship which can be applied to all the patients. As a new patient comes, we predict his/her facial deformation by combining the general relationship and the new patient’s biomechanical properties. Thus, our model is biomechanical relevant and statistical relevant. Validation on eleven patients demonstrates the effectiveness and efficiency of our method.
kernel ridge regression; finite element model; maxillofacial surgery; soft tissue deformation
Recent advances in cone-beam computed tomography (CBCT) have rapidly enabled widepsread applications of dentomaxillofacial imaging and orthodontic practices in the past decades due to its low radiation dose, high spatial resolution, and accessibility. However, low contrast resolution in CBCT image has become its major limitation in building skull models. Intensive hand-segmentation is usually required to reconstruct the skull models. One of the regions affected by this limitation the most is the thin bone images. This paper presents a novel segmentation approach based on wavelet density model (WDM) for a particular interest in the outer surface of anterior wall of maxilla. Nineteen CBCT datasets are used to conduct two experiments. This mode-based segmentation approach is validated and compared with three different segmentation approaches. The results show that the performance of this model-based segmentation approach is better than those of the other approaches. It can achieve 0.25 ± 0.2mm of surface error from ground truth of bone surface.
3D segmentation; Active shape model (ASM); Statistical shape model (SSM); Craniomaxillofacial (CMF) surgeries; Cone-beam computed tomography (CBCT)
The rapid development of next generation sequencing (NGS) technology provides a novel avenue for genomic exploration and research. Single nucleotide variants (SNVs) inferred from next generation sequencing are expected to reveal gene mutations in cancer. However, NGS has lower sequence coverage and poor SNVs detection capability in the regulatory regions of the genome. Post probabilistic based methods are efficient for detection of SNVs in high coverage regions or sequencing data with high depth. However, for data with low sequencing depth, the efficiency of such algorithms remains poor and needs to be improved.
A new tool SNVHMM basing on a discrete hidden Markov model (HMM) was developed to infer the genotype for each position on the genome. We incorporated the mapping quality of each read and the corresponding base quality on the reads into the emission probability of HMM. The context information of the whole observation as well as its confidence were completely utilized to infer the genotype for each position on the genome in study. Therefore, more probability power can be gained over the Bayes based methods, which is very useful for SNVs detection for data with low sequencing depth. Moreover, our model was verified by testing against two sets of lobular breast tumor and Myelodysplastic Syndromes (MDS) data each. Comparing against a recently published SNVs calling algorithm SNVMix2, our model improved the performance of SNVMix2 largely when the sequencing depth is low and also outperformed SNVMix2 when SNVMix2 is well trained by large datasets.
SNVHMM can detect SNVs from NGS cancer data efficiently even if the sequence depth is very low. The training data size can be very small for SNVHMM to work. SNVHMM incorporated the base quality and mapping quality of all observed bases and reads, and also provides the option for users to choose the confidence of the observation for SNVs prediction.
Recently, melanoma has become the most malignant and commonly occurring skin cancer. Melanoma is not only the major source (75%) of deaths related to skin cancer, but also it is hard to be treated by the conventional drugs. Recent research indicated that angiogenesis is an important factor for tumor initiation, expansion, and response to therapy. Thus, we proposed a novel multi-scale agent-based computational model that integrates the angiogenesis into tumor growth to study the response of melanoma cancer under combined drug treatment.
Our multi-scale agent-based model can simulate the melanoma tumor growth with angiogenesis under combined drug treatment. The significant synergistic effects between drug Dox and drug Sunitinib demonstrated the clinical potential to interrupt the communication between melanoma cells and its related vasculatures. Also, the sensitivity analysis of the model revealed that diffusivity related to the micro-vasculatures around tumor tissues closely correlated with the spread, oscillation and destruction of the tumor.
Simulation results showed that the 3D model can represent key features of melanoma growth, angiogenesis, and its related micro-environment. The model can help cancer researchers understand the melanoma developmental mechanism. Drug synergism analysis suggested that interrupting the communications between melanoma cells and the related vasculatures can significantly increase the drug efficacy against tumor cells.
Microenvironment; Drug synergism; Agent-based model; Multi-scale; Melanoma; Anti-angiogenesis
Gene fusions, which result from abnormal chromosome rearrangements, are a pathogenic factor in cancer development. The emerging RNA-Seq technology enables us to detect gene fusions and profile their features.
In this paper, we proposed a novel fusion detection tool, FusionQ, based on paired-end RNA-Seq data. This tool can detect gene fusions, construct the structures of chimerical transcripts, and estimate their abundances. To confirm the read alignment on both sides of a fusion point, we employed a new approach, “residual sequence extension”, which extended the short segments of the reads by aggregating their overlapping reads. We also proposed a list of filters to control the false-positive rate. In addition, we estimated fusion abundance using the Expectation-Maximization algorithm with sparse optimization, and further adopted it to improve the detection accuracy of the fusion transcripts. Simulation was performed by FusionQ and another two stated-of-art fusion detection tools. FusionQ exceeded the other two in both sensitivity and specificity, especially in low coverage fusion detection. Using paired-end RNA-Seq data from breast cancer cell lines, FusionQ detected both the previously reported and new fusions. FusionQ reported the structures of these fusions and provided their expressions. Some highly expressed fusion genes detected by FusionQ are important biomarkers in breast cancer. The performances of FusionQ on cancel line data still showed better specificity and sensitivity in the comparison with another two tools.
FusionQ is a novel tool for fusion detection and quantification based on RNA-Seq data. It has both good specificity and sensitivity performance. FusionQ is free and available at http://www.wakehealth.edu/CTSB/Software/Software.htm.
Fusion detection; chimerical transcripts quantification; EM algorithm
Multiple intergenic single-nucleotide polymorphisms (SNPs) near hedgehog interacting protein (HHIP) on chromosome 4q31 have been strongly associated with pulmonary function levels and moderate-to-severe chronic obstructive pulmonary disease (COPD). However, whether the effects of variants in this region are related to HHIP or another gene has not been proven. We confirmed genetic association of SNPs in the 4q31 COPD genome-wide association study (GWAS) region in a Polish cohort containing severe COPD cases and healthy smoking controls (P = 0.001 to 0.002). We found that HHIP expression at both mRNA and protein levels is reduced in COPD lung tissues. We identified a genomic region located ∼85 kb upstream of HHIP which contains a subset of associated SNPs, interacts with the HHIP promoter through a chromatin loop and functions as an HHIP enhancer. The COPD risk haplotype of two SNPs within this enhancer region (rs6537296A and rs1542725C) was associated with statistically significant reductions in HHIP promoter activity. Moreover, rs1542725 demonstrates differential binding to the transcription factor Sp3; the COPD-associated allele exhibits increased Sp3 binding, which is consistent with Sp3's usual function as a transcriptional repressor. Thus, increased Sp3 binding at a functional SNP within the chromosome 4q31 COPD GWAS locus leads to reduced HHIP expression and increased susceptibility to COPD through distal transcriptional regulation. Together, our findings reveal one mechanism through which SNPs upstream of the HHIP gene modulate the expression of HHIP and functionally implicate reduced HHIP gene expression in the pathogenesis of COPD.
Avian tembusu virus (TMUV), which was first identified in eastern China, is an emerging virus causing serious economic losses in the Chinese poultry industry. Here, we report the complete genome sequence of goose tembusu virus strain JS804, isolated from Jiangnan white geese with severe neurological signs. The genome of JS804 is 10,990 nucleotides (nt) in length and contains a single open reading frame encoding a putative polyprotein of 3,425 amino acids. Research of the whole sequence of tembusu virus will help us to understand further the molecular and evolutionary characteristics and pathogenesis of this virus.
S-glutathionylation, the reversible formation of mixed disulfides between glutathione(GSH) and cysteine residues in proteins, is a specific form of post-translational modification that plays important roles in various biological processes, including signal transduction, redox homeostasis, and metabolism inside cells. Experimentally identifying S-glutathionylation sites is labor-intensive and time consuming, whereas bioinformatics methods provide an alternative way to this problem by predicting S-glutathionylation sites in silico. The bioinformatics approaches give not only candidate sites for further experimental verification but also bio-chemical insights into the mechanism of S-glutathionylation. In this paper, we firstly collect experimentally determined S-glutathionylated proteins and their corresponding modification sites from the literature, and then propose a new method for predicting S-glutathionylation sites by employing machine learning methods based on protein sequence data. Promising results are obtained by our method with an AUC (area under ROC curve) score of 0.879 in 5-fold cross-validation, which demonstrates the predictive power of our proposed method. The datasets used in this work are available at http://csb.shu.edu.cn/SGDB.
Despite successful approaches to preserve organs, tissues, and isolated cells, the maintenance of stem cell viability and function in body fluids during storage for cell distribution and transportation remains unexplored. The aim of this study was to characterize urine-derived stem cells (USCs) after optimal preservation of urine specimens for up to 24 hours. A total of 415 urine specimens were collected from 12 healthy men (age range 20–54 years old). About 6×104 cells shed off from the urinary tract system in 24 hours. At least 100 USC clones were obtained from the stored urine specimens after 24 hours and maintained similar biological features to fresh USCs. The stored USCs had a “rice grain” shape in primary culture, and expressed mesenchymal stem cell surface markers, high telomerase activity, and normal karyotypes. Importantly, the preserved cells retained bipotent differentiation capacity. Differentiated USCs expressed myogenic specific proteins and contractile function when exposed to myogenic differentiation medium, and they expressed urothelial cell-specific markers and barrier function when exposed to urothelial differentiation medium. These data demonstrated that up to 75% of fresh USCs can be safely persevered in urine for 24 hours and that these cells stored in urine retain their original stem cell properties, indicating that preserved USCs could be available for potential use in cell-based therapy or clinical diagnosis.
RNA interference (RNAi) becomes an increasingly important and effective genetic tool to study the function of target genes by suppressing specific genes of interest. This system approach helps identify signaling pathways and cellular phase types by tracking intensity and/or morphological changes of cells. The traditional RNAi screening scheme, in which one siRNA is designed to knockdown one specific mRNA target, needs a large library of siRNAs and turns out to be time-consuming and expensive.
In this paper, we propose a conceptual model, called compressed sensing RNAi (csRNAi), which employs a unique combination of group of small interfering RNAs (siRNAs) to knockdown a much larger size of genes. This strategy is based on the fact that one gene can be partially bound with several small interfering RNAs (siRNAs) and conversely, one siRNA can bind to a few genes with distinct binding affinity. This model constructs a multi-to-multi correspondence between siRNAs and their targets, with siRNAs much fewer than mRNA targets, compared with the conventional scheme. Mathematically this problem involves an underdetermined system of equations (linear or nonlinear), which is ill-posed in general. However, the recently developed compressed sensing (CS) theory can solve this problem. We present a mathematical model to describe the csRNAi system based on both CS theory and biological concerns. To build this model, we first search nucleotide motifs in a target gene set. Then we propose a machine learning based method to find the effective siRNAs with novel features, such as image features and speech features to describe an siRNA sequence. Numerical simulations show that we can reduce the siRNA library to one third of that in the conventional scheme. In addition, the features to describe siRNAs outperform the existing ones substantially.
This csRNAi system is very promising in saving both time and cost for large-scale RNAi screening experiments which may benefit the biological research with respect to cellular processes and pathways.
The epidermal growth factor receptor (EGFR) signaling pathway and angiogenesis in brain cancer act as an engine for tumor initiation, expansion and response to therapy. Since the existing literature does not have any models that investigate the impact of both angiogenesis and molecular signaling pathways on treatment, we propose a novel multi-scale, agent-based computational model that includes both angiogenesis and EGFR modules to study the response of brain cancer under tyrosine kinase inhibitors (TKIs) treatment.
The novel angiogenesis module integrated into the agent-based tumor model is based on a set of reaction–diffusion equations that describe the spatio-temporal evolution of the distributions of micro-environmental factors such as glucose, oxygen, TGFα, VEGF and fibronectin. These molecular species regulate tumor growth during angiogenesis. Each tumor cell is equipped with an EGFR signaling pathway linked to a cell-cycle pathway to determine its phenotype. EGFR TKIs are delivered through the blood vessels of tumor microvasculature and the response to treatment is studied.
Our simulations demonstrated that entire tumor growth profile is a collective behaviour of cells regulated by the EGFR signaling pathway and the cell cycle. We also found that angiogenesis has a dual effect under TKI treatment: on one hand, through neo-vasculature TKIs are delivered to decrease tumor invasion; on the other hand, the neo-vasculature can transport glucose and oxygen to tumor cells to maintain their metabolism, which results in an increase of cell survival rate in the late simulation stages.
Multi-scale; Agent-based modeling; EGFR signaling pathway; Angiogenesis; TKI treatment
In Traditional Chinese Medicine (TCM), the lip diagnosis is an important diagnostic method which has a long history and is applied widely. The lip color of a person is considered as a symptom to reflect the physical conditions of organs in the body. However, the traditional diagnostic approach is mainly based on observation by doctor’s nude eyes, which is non-quantitative and subjective. The non-quantitative approach largely depends on the doctor’s experience and influences accurate the diagnosis and treatment in TCM. Developing new quantification methods to identify the exact syndrome based on the lip diagnosis of TCM becomes urgent and important. In this paper, we design a computer-assisted classification model to provide an automatic and quantitative approach for the diagnosis of TCM based on the lip images.
A computer-assisted classification method is designed and applied for syndrome diagnosis based on the lip images. Our purpose is to classify the lip images into four groups: deep-red, red, purple and pale. The proposed scheme consists of four steps including the lip image preprocessing, image feature extraction, feature selection and classification. The extracted 84 features contain the lip color space component, texture and moment features. Feature subset selection is performed by using SVM-RFE (Support Vector Machine with recursive feature elimination), mRMR (minimum Redundancy Maximum Relevance) and IG (information gain). Classification model is constructed based on the collected lip image features using multi-class SVM and Weighted multi-class SVM (WSVM). In addition, we compare SVM with k-nearest neighbor (kNN) algorithm, Multiple Asymmetric Partial Least Squares Classifier (MAPLSC) and Naïve Bayes for the diagnosis performance comparison. All displayed faces image have obtained consent from the participants.
A total of 257 lip images are collected for the modeling of lip diagnosis in TCM. The feature selection method SVM-RFE selects 9 important features which are composed of 5 color component features, 3 texture features and 1 moment feature. SVM, MAPLSC, Naïve Bayes, kNN showed better classification results based on the 9 selected features than the results obtained from all the 84 features. The total classification accuracy of the five methods is 84%, 81%, 79% and 81%, 77%, respectively. So SVM achieves the best classification accuracy. The classification accuracy of SVM is 81%, 71%, 89% and 86% on Deep-red, Pale Purple, Red and lip image models, respectively. While with the feature selection algorithm mRMR and IG, the total classification accuracy of WSVM achieves the best classification accuracy. Therefore, the results show that the system can achieve best classification accuracy combined with SVM classifiers and SVM-REF feature selection algorithm.
A diagnostic system is proposed, which firstly segments the lip from the original facial image based on the Chan-Vese level set model and Otsu method, then extracts three kinds of features (color space features, Haralick co-occurrence features and Zernike moment features) on the lip image. Meanwhile, SVM-REF is adopted to select the optimal features. Finally, SVM is applied to classify the four classes. Besides, we also compare different feature selection algorithms and classifiers to verify our system. So the developed automatic and quantitative diagnosis system of TCM is effective to distinguish four lip image classes: Deep-red, Purple, Red and Pale. This study puts forward a new method and idea for the quantitative examination on lip diagnosis of TCM, as well as provides a template for objective diagnosis in TCM.
Traditional chinese medicine; Computer-assisted lip diagnosis; Image analysis; Feature selection; Support vector machine
Recent research in cancer biology has suggested the hypothesis that tumors are initiated and driven by a small group of cancer stem cells (CSCs). Furthermore, cancer stem cell niches have been found to be essential in determining fates of CSCs, and several signaling pathways have been proven to play a crucial role in cellular behavior, which could be two important factors in cancer development. To better understand the progression, heterogeneity and treatment response of breast cancer, especially in the context of CSCs, we propose a mathematical model based on the cell compartment method. In this model, three compartments of cellular subpopulations are constructed: CSCs, progenitor cells (PCs), and terminal differentiated cells (TCs). Moreover, 1) the cancer stem cell niche is, considered by modeling its effect on division patterns (symmetric or asymmetric) of CSCs, and 2) the EGFR signaling pathway is integrated by modeling its role in cell proliferation, apoptosis. Our simulation results indicate that 1) a higher probability for symmetric division of CSC may result in a faster expansion of tumor population, and for a larger number of niches, the tumor grows at a slower rate, but the final tumor volume is larger; 2) higher EGFR expression correlates to tumors with larger volumes while a saturation function is observed, and 3) treatments that inhibit tyrosine kinase activity of EGFR may not only repress the tumor volume, but also decrease the CSCs percentages by shifting CSCs from symmetric divisions to asymmetric divisions. These findings suggest that therapies should be designed to effectively control or eliminate the symmetric division of CSCs and to reduce or destroy the CSC niches.
mathematical model; compartment method; signaling pathway; breast cancer; tyrosine kinase inhibitors
Multiscale agent-based modeling (MABM) has been widely used to simulate Glioblastoma Multiforme (GBM) and its progression. At the intracellular level, the MABM approach employs a system of ordinary differential equations to describe quantitatively specific intracellular molecular pathways that determine phenotypic switches among cells (e.g. from migration to proliferation and vice versa). At the intercellular level, MABM describes cell-cell interactions by a discrete module. At the tissue level, partial differential equations are employed to model the diffusion of chemoattractants, which are the input factors of the intracellular molecular pathway. Moreover, multiscale analysis makes it possible to explore the molecules that play important roles in determining the cellular phenotypic switches that in turn drive the whole GBM expansion. However, owing to limited computational resources, MABM is currently a theoretical biological model that uses relatively coarse grids to simulate a few cancer cells in a small slice of brain cancer tissue. In order to improve this theoretical model to simulate and predict actual GBM cancer progression in real time, a graphics processing unit (GPU)-based parallel computing algorithm was developed and combined with the multi-resolution design to speed up the MABM. The simulated results demonstrated that the GPU-based, multi-resolution and multiscale approach can accelerate the previous MABM around 30-fold with relatively fine grids in a large extracellular matrix. Therefore, the new model has great potential for simulating and predicting real-time GBM progression, if real experimental data are incorporated.
Chronic lymphocytic leukemia deletion gene 7 (Clld7) is a candidate tumor suppressor on chromosome 13q14. Clld7 encodes an evolutionarily conserved protein that contains an RCC1 domain plus broad complex, tramtrack, bric-a-brac (BTB) and POZ domains. In this study, we investigated the biological functions of Clld7 protein in inducible osteosarcoma cell lines. Clld7 induction inhibited cell growth, decreased cell viability, and increased gamma-H2AX staining under conditions of caspase inhibition, indicating activation of the DNA damage/repair pathway. Real-time PCR analysis in tumor cells and normal human epithelial cells revealed Clld7 target genes that regulate DNA repair responses. Furthermore, depletion of Clld7 in normal human epithelial cells conferred resistance to apoptosis triggered by DNA damage. Taken together, the biological actions of Clld7 are consistent with those of a tumor suppressor.
tumor suppressor; DNA damage/repair pathway; 13q14; primary human keratinocyte; apoptosis
Motivation: High-throughput protein interaction data, with ever-increasing volume, are becoming the foundation of many biological discoveries, and thus high-quality protein–protein interaction (PPI) maps are critical for a deeper understanding of cellular processes. However, the unreliability and paucity of current available PPI data are key obstacles to the subsequent quantitative studies. It is therefore highly desirable to develop an approach to deal with these issues from the computational perspective. Most previous works for assessing and predicting protein interactions either need supporting evidences from multiple information resources or are severely impacted by the sparseness of PPI networks.
Results: We developed a robust manifold embedding technique for assessing the reliability of interactions and predicting new interactions, which purely utilizes the topological information of PPI networks and can work on a sparse input protein interactome without requiring additional information types. After transforming a given PPI network into a low-dimensional metric space using manifold embedding based on isometric feature mapping (ISOMAP), the problem of assessing and predicting protein interactions is recasted into the form of measuring similarity between points of its metric space. Then a reliability index, a likelihood indicating the interaction of two proteins, is assigned to each protein pair in the PPI networks based on the similarity between the points in the embedded space. Validation of the proposed method is performed with extensive experiments on densely connected and sparse PPI network of yeast, respectively. Results demonstrate that the interactions ranked top by our method have high-functional homogeneity and localization coherence, especially our method is very efficient for large sparse PPI network with which the traditional algorithms fail. Therefore, the proposed algorithm is a much more promising method to detect both false positive and false negative interactions in PPI networks.
Availability: MATLAB code implementing the algorithm is available from the web site http://home.ustc.edu.cn/∼yzh33108/Manifold.htm.
Supplementary information: Supplementary data are available at Bioinformatics online.
High content neuron image processing is considered as an important method for quantitative neurobiological studies. The main goal of analysis in this paper is to provide automatic image processing approaches to process neuron images for studying neuron mechanism in high content screening. In the nuclei channel, all nuclei are segmented and detected by applying the gradient vector field based watershed. Then the neuronal nuclei are selected based on the soma region detected in neurite channel. In neurite images, we propose a novel neurite centerline extraction approach using the improved line-pixel detection technique. The proposed neurite tracing method can detect the curvilinear structure more accurately compared with the current existing methods. An interface called NeuriteIQ based on the proposed algorithms is developed finally for better application in high content screening.
High content screening; Microscopy image; Nuclei segmentation; Neurite outgrowth; Line-pixel detection; Branch area
Cancer cells with active drug-efflux capability are multidrug resistant and pose a significant obstacle for the efficacy of chemotherapy. Moreover, recent evidence suggests that high drug-efflux cancer cells (HDECCs) may be selectively enriched with stem-like cancer cells, which are believed to be the cause for tumor initiation and recurrence. There is a great need for therapeutic reagents that are capable of eliminating HDECCs. We developed an image-based high-content screening (HCS) system to specifically identify and analyze the HDECC population in lung cancer cells. Using the system, we screened 1,280 pharmacologically active compounds which identified twelve potent HDECC inhibitors. It is shown that these inhibitors are able to overcome MDR and sensitize HDECCs to chemotherapeutic drugs, or directly reduce the tumorigenicity of lung cancer cells possibly by affecting stem-like cancer cells. The HCS system we established provides a new approach for identifying therapeutic reagents overcoming MDR. The compounds identified by the screening may potentially be used as potential adjuvant to improve the efficacy of chemotherapeutic drugs.
high drug-efflux cancer cells; multidrug resistance; high content screening; image-based assay
Automated time-lapse microscopy can visualize proliferation of large numbers of individual cells, enabling accurate measurement of the frequency of cell division and the duration of interphase and mitosis. However, extraction of quantitative information by manual inspection of time-lapse movies is too time-consuming to be useful for analysis of large experiments.
Here we present an automated time-series approach that can measure changes in the duration of mitosis and interphase in individual cells expressing fluorescent histone 2B. The approach requires analysis of only 2 features, nuclear area and average intensity. Compared to supervised learning approaches, this method reduces processing time and does not require generation of training data sets. We demonstrate that this method is as sensitive as manual analysis in identifying small changes in interphase or mitotic duration induced by drug or siRNA treatment.
This approach should facilitate automated analysis of high-throughput time-lapse data sets to identify small molecules or gene products that influence timing of cell division.
Recent studies on neuron imaging show that there is a strong relationship between the functional properties of a neuron and its morphology, especially its dendritic spine structures. However, most of the current methods for morphological spine classification only concern features in two-dimensional (2D) space, which consequently decreases the accuracy of dendritic spine analysis. In this paper, we propose a semi-supervised learning (SSL) framework, in which spine phenotypes in three-dimensional (3D) space are considered. With training only on a few pre-classified inputs, the rest of the spines can be identified effectively. We also derived a new scheme using an affinity matrix between features to further improve the accuracy. Our experimental results indicate that a small training dataset is sufficient to classify detected dendritic spines.
dendritic spine; semi-supervised learning; morphological spine classification
In this brief paper we present an overview of the TSC-mTOR pathway and its importance in neurodegenerative disease (ND). We illustrate the influence of ND on dendritic spine morphology. Then we discuss some details of functional gene networks (FGN) and use this information to propose an image driven systems biology approach for the construction of a FGN for ND. We conclude on its importance and the prospective outcome of our study.
In addition to intrinsic genetic alterations, the effects of the extrinsic microenvironment also play a pathological role in cancer development. Altered chemokine/cytokine networks in the tumor microenvironment may contribute to the dysregulation of cellular functions in cancer cells. Anaplastic large cell lymphoma (ALCL) is an aggressive T-cell lymphoma caused by abnormal expression of anaplastic lymphoma kinase due to a chromosomal translocation. Notably, ALCL cells are also characterized by high-level expression of the high-affinity IL-2 receptor subunit CD25 on the cell surface. However, whether the IL-2/IL-2 receptor functions in ALCL cells and how this signaling affects the tumor remain unclear. In this study, we treated cultured ALCL cells with exogenous IL-2 and examined changes in cellular function and signaling pathways. IL-2 stimulated cell growth and augmented activation of the extracellular signal-regulated kinases 1/2 (ERK1/2) pathway. Additionally, IL-2 enhanced lymphoma cell survival by overcoming kinase inhibitor U0126-induced growth arrest and apoptosis. Subsequently, to identify the potential source of IL-2 for lymphoma cells in vivo, we performed gene expression and immunochemical analyses. RT-PCR revealed no IL-2 gene expression in cultured ALCL cells and ruled out the possibility of an IL-2 autocrine loop. Interestingly, immunostaining of lymphoma tumor tissues showed IL-2 protein expression in background cells within tumor tissue, but not in ALCL cells. Our findings demonstrate that IL-2 signaling plays a functional role in ALCL cells, and enhances lymphoma cell survival by increasing activation of the ERK1/2 pathway.
anaplastic large cell lymphoma (ALCL); CD25/IL-2 receptor; IL-2 signaling; extracellular signal-regulated kinases (ERK1/2); tumor microenvironment
Motivation: Prediction of synergistic effects of drug combinations has traditionally been relied on phenotypic response data. However, such methods cannot be used to identify molecular signaling mechanisms of synergistic drug combinations. In this article, we propose an enhanced Petri-Net (EPN) model to recognize the synergistic effects of drug combinations from the molecular response profiles, i.e. drug-treated microarray data.
Methods: We addressed the downstream signaling network of the targets for the two individual drugs used in the pairwise combinations and applied EPN to the identified targeted signaling network. In EPN, drugs and signaling molecules are assigned to different types of places, while drug doses and molecular expressions are denoted by color tokens. The changes of molecular expressions caused by treatments of drugs are simulated by two actions of EPN: firing and blasting. Firing is to transit the drug and molecule tokens from one node or place to another, and blasting is to reduce the number of molecule tokens by drug tokens in a molecule node. The goal of EPN is to mediate the state characterized by control condition without any treatment to that of treatment and to depict the drug effects on molecules by the drug tokens.
Results: We applied EPN to our generated pairwise drug combination microarray data. The synergistic predictions using EPN are consistent with those predicted using phenotypic response data. The molecules responsible for the synergistic effects with their associated feedback loops display the mechanisms of synergism.
Availability: The software implemented in Python 2.7 programming language is available from request.
The development of techniques for sequencing the messenger RNA (RNA-Seq) enables it to study the biological mechanisms such as alternative splicing and gene expression regulation more deeply and accurately. Most existing methods employ RNA-Seq to quantify the expression levels of already annotated isoforms from the reference genome. However, the current reference genome is very incomplete due to the complexity of the transcriptome which hiders the comprehensive investigation of transcriptome using RNA-Seq. Novel study on isoform inference and estimation purely from RNA-Seq without annotation information is desirable.
A Nonnegativity and Sparsity constrained Maximum APosteriori (NSMAP) model has been proposed to estimate the expression levels of isoforms from RNA-Seq data without the annotation information. In contrast to previous methods, NSMAP performs identification of the structures of expressed isoforms and estimation of the expression levels of those expressed isoforms simultaneously, which enables better identification of isoforms. In the simulations parameterized by two real RNA-Seq data sets, more than 77% expressed isoforms are correctly identified and quantified. Then, we apply NSMAP on two RNA-Seq data sets of myelodysplastic syndromes (MDS) samples and one normal sample in order to identify differentially expressed known and novel isoforms in MDS disease.
NSMAP provides a good strategy to identify and quantify novel isoforms without the knowledge of annotated reference genome which can further realize the potential of RNA-Seq technique in transcriptome analysis. NSMAP package is freely available at https://sites.google.com/site/nsmapforrnaseq.
In this paper, we address the issue of peptide ion peak detection for high resolution time-of-flight (TOF) mass spectrometry (MS) data. A novel Bayesian peptide ion peak detection method is proposed for TOF data with resolution of 10 000–15 000 full width at half-maximum (FWHW). MS spectra exhibit distinct characteristics at this resolution, which are captured in a novel parametric model. Based on the proposed parametric model, a Bayesian peak detection algorithm based on Markov chain Monte Carlo (MCMC) sampling is developed. The proposed algorithm is tested on both simulated and real datasets. The results show a significant improvement in detection performance over a commonly employed method. The results also agree with expert’s visual inspection. Moreover, better detection consistency is achieved across MS datasets from patients with identical pathological condition.
Bayesian methods; Markov chain Monte Carlo; mass spectrometry; peptide peak detection; time-of-flight