Functional clustering of time series gene expression data by Granger causality

Time-Series Clustering of Single-Cell Trajectories in Collective Cell Migration

Simple Summary

In this study, we normalized trajectories containing both mesenchymal and epithelial cells to remove the effect of cell location on clustering, and performed a dimensionality reduction on the time series data before clustering. When the clustering results were superimposed on the trajectories prior to normalization, the results still showed similarities in location, indicating that this method can find cells with similar migration patterns. These data highlight the reliability of this method in identifying consistent migration patterns in collective cell migration.

Abstract

Collective invasion drives multicellular cancer cells to spread to surrounding normal tissues. To fully comprehend metastasis, the methodology of analysis of individual cell migration in tissue should be well developed. Extracting and classifying cells with similar migratory characteristics in a colony would facilitate an understanding of complex cell migration patterns. Here, we used electrospun fibers as the extracellular matrix for the in vitro modeling of collective cell migration, clustering of mesenchymal and epithelial cells based on trajectories, and analysis of collective migration patterns based on trajectory similarity. We normalized the trajectories to eliminate the effect of cell location on clustering and used uniform manifold approximation and projection to perform dimensionality reduction on the time-series data before clustering. When the clustering results were superimposed on the trajectories before normalization, the results still exhibited positional similarity, thereby demonstrating that this method can identify cells with similar migration patterns. The same cluster contained both mesenchymal and epithelial cells, and this result was related to cell location and cell division. These data highlight the reliability of this method in identifying consistent migration patterns during collective cell migration. This provides new insights into the epithelial–mesenchymal interactions that affect migration patterns.

Dense time-course gene expression profiling of the Drosophila melanogaster innate immune response

BACKGROUND

Immune responses need to be initiated rapidly, and maintained as needed, to prevent establishment and growth of infections. At the same time, resources need to be balanced with other physiological processes. On the level of transcription, studies have shown that this balancing act is reflected in tight control of the initiation kinetics and shutdown dynamics of specific immune genes.

RESULTS

To investigate genome-wide expression dynamics and trade-offs after infection at a high temporal resolution, we performed an RNA-seq time course on D. melanogaster with 20 time points post Imd stimulation. A combination of methods, including spline fitting, cluster analysis, and Granger causality inference, allowed detailed dissection of expression profiles, lead-lag interactions, and functional annotation of genes through guilt-by-association. We identified Imd-responsive genes and co-expressed, less well characterized genes, with an immediate-early response and sustained up-regulation up to 5 days after stimulation. In contrast, stress response and Toll-responsive genes, among which were Bomanins, demonstrated early and transient responses. We further observed a strong trade-off with metabolic genes, which strikingly recovered to pre-infection levels before the immune response was fully resolved.

CONCLUSIONS

This high-dimensional dataset enabled the comprehensive study of immune response dynamics through the parallel application of multiple temporal data analysis methods. The well annotated data set should also serve as a useful resource for further investigation of the D. melanogaster innate immune response, and for the development of methods for analysis of a post-stress transcriptional response time-series at whole-genome scale.

Discovering key transcriptomic regulators in pancreatic ductal adenocarcinoma using Dirichlet process Gaussian mixture model

Pancreatic Ductal Adenocarcinoma (PDAC) is the most lethal type of pancreatic cancer, late detection leading to its therapeutic failure. This study aims to determine the key regulatory genes and their impacts on the disease’s progression, helping the disease’s etiology, which is still mostly unknown. We leverage the landmark advantages of time-series gene expression data of this disease and thereby identified the key regulators that capture the characteristics of gene activity patterns in the cancer progression. We have identified the key gene modules and predicted the functions of top genes from a reconstructed gene association network (GAN). A variation of the partial correlation method is utilized to analyze the GAN, followed by a gene function prediction task. Moreover, we have identified regulators for each target gene by gene regulatory network inference using the dynamical GENIE3 (dynGENIE3) algorithm. The Dirichlet process Gaussian process mixture model and cubic spline regression model (splineTimeR) are employed to identify the key gene modules and differentially expressed genes, respectively. Our analysis demonstrates a panel of key regulators and gene modules that are crucial for PDAC disease progression.

Clustering of Bacterial Growth Dynamics in Response to Growth Media by Dynamic Time Warping

Bacterial growth curves, representing population dynamics, are still poorly understood. The growth curves are commonly analyzed by model-based theoretical fitting, which is limited to typical S-shape fittings and does not elucidate the dynamics in their entirety. Thus, whether a certain growth condition results in any particular pattern of growth curve remains unclear. To address this question, up-to-date data mining techniques were applied to bacterial growth analysis for the first time. Dynamic time warping (DTW) and derivative DTW (DDTW) were used to compare the similarity among 1015 growth curves of 28 Escherichia coli strains growing in three different media. In the similarity evaluation, agglomerative hierarchical clustering, assessed with four statistic benchmarks, successfully categorized the growth curves into three clusters, roughly corresponding to the three media. Furthermore, a simple benchmark was newly proposed, providing a highly improved accuracy (~99%) in clustering the growth curves corresponding to the growth media. The biologically reasonable categorization of growth curves suggested that DTW and DDTW are applicable for bacterial growth analysis. The bottom-up clustering results indicate that the growth media determine some specific patterns of population dynamics, regardless of genomic variation, and thus have a higher priority of shaping the growth curves than the genomes do.

A Composite Mode Differential Gene Regulatory Architecture based on Temporal Expression Profiles

Exploring the complex interactive mechanism in a Gene Regulatory Network (GRN) developed using transcriptome data obtained from standard microarray and/or RNA-seq experiments helps us to understand the triggering factors in cancer research. The Transcription Factor (TF) genes generate protein complexes which affect the transcription of various target genes. However, considering the mode of regulation in a time frame such transcriptional activities are dependent on some specific activation time points only. It is also crucial to check whether the regulating capabilities are uniform across varied stages, especially when periodicity is a big issue. In this context, we propose an algorithm called RIFT which helps to monitor the temporal differential regulatory pattern of a Differentially Expressed (DE) target gene either by a TF gene or a group of TF genes from a large time series (TS) data. We have tested our algorithm on HeLa cell cycle data and compared the result with its most advanced state of the art counterpart proposed so far. As our algorithm yields up stringent mode and target specific significant valid TF genes for a DE gene, we can expect to have new forms of genetic interactions.

Neonatal brain resting-state functional connectivity imaging modalities

Infancy is the most critical period in human brain development. Studies demonstrate that subtle brain abnormalities during this state of life may greatly affect the developmental processes of the newborn infants. One of the rapidly developing methods for early characterization of abnormal brain development is functional connectivity of the brain at rest. While the majority of resting-state studies have been conducted using magnetic resonance imaging (MRI), there is clear evidence that resting-state functional connectivity (rs-FC) can also be evaluated using other imaging modalities. The aim of this review is to compare the advantages and limitations of different modalities used for the mapping of infants' brain functional connectivity at rest. In addition, we introduce photoacoustic tomography, a novel functional neuroimaging modality, as a complementary modality for functional mapping of infants' brain.

DNA Damage Is a Potential Marker for TP53 Mutation in Colorectal Carcinogenesis

PURPOSE

The ability to measure oxidative DNA damage in a tissue allows establishment of the relationship between DNA damage and mutations in normal and neoplastic cells. It is well known that TP53 is a key inhibitor of tumor development and preserves the genome integrity in each cell. The aim of the present study was to investigate the relationship between DNA damage and TP53 mutation in colorectal adenoma and adenocarcinoma, and the value of DNA damage as potential marker of TP53 mutation in non-tumor tissues adjacent to colon malignant lesions.

METHODS

Tissue samples were obtained by colonoscopy from patients with adenoma and/or adenocarcinoma and from healthy volunteers. Diagnosis was defined by histopathology. Immunohistochemistry with computer-assisted image analysis was performed to quantify TP53 mutation. Oxidative DNA damage was determined by comet assay. Statistical analyses were performed with 5 % of significance level.

RESULTS

The TP53 level was higher in non-tumor tissues from tumor patients than in normal tissues from healthy volunteers (p = 0.01). Likewise, higher TP53 levels were observed in tumor tissues compared with the non-tumor tissues (p = 0.00). Oxidative DNA damage levels were higher in tumor tissues than in non-tumor tissues (p = 0.00). The amount of TP53 (p = 0.00) and oxidative DNA damage (p = 0.00) in normal and tumor tissue was related. The relationship between oxidative DNA damage and TP53 mutation was demonstrated in all samples (p = 0.00).

CONCLUSION

Oxidative DNA damage is an intervening variable for TP53 mutation in colorectal adenoma-carcinoma. Our data suggests that oxidative DNA damage is a potential marker of TP53 mutation in colorectal carcinogenesis.

A tutorial to identify nonlinear associations in gene expression time series data

The study of gene regulatory networks is the basis to understand the biological complexity of several diseases and/or cell states. It has become the core of research in the field of systems biology. Several mathematical methods have been developed in the last decade, especially in the analysis of time series gene expression data derived from microarrays and sequencing-based methods. Most of the models available in the literature assumes linear associations among genes and do not infer directionality in these connections or uses a priori biological knowledge to set the directionality. However, in several cases, a priori biological information is not available. In this context, we describe a statistical method, namely nonlinear vector autoregressive model to estimate nonlinear relationships and also to infer directionality at the edges of the network by using the temporal information of the time series gene expression data without a priori biological information.

OLYMPUS: an automated hybrid clustering method in time series gene expression. Case study: host response after Influenza A (H1N1) infection

The increasing flow of short time series microarray experiments for the study of dynamic cellular processes poses the need for efficient clustering tools. These tools must deal with three primary issues: first, to consider the multi-functionality of genes; second, to evaluate the similarity of the relative change of amplitude in the time domain rather than the absolute values; third, to cope with the constraints of conventional clustering algorithms such as the assignment of the appropriate cluster number. To address these, we propose OLYMPUS, a novel unsupervised clustering algorithm that integrates Differential Evolution (DE) method into Fuzzy Short Time Series (FSTS) algorithm with the scope to utilize efficiently the information of population of the first and enhance the performance of the latter. Our hybrid approach provides sets of genes that enable the deciphering of distinct phases in dynamic cellular processes. We proved the efficiency of OLYMPUS on synthetic as well as on experimental data. The discriminative power of OLYMPUS provided clusters, which refined the so far perspective of the dynamics of host response mechanisms to Influenza A (H1N1). Our kinetic model sets a timeline for several pathways and cell populations, implicated to participate in host response; yet no timeline was assigned to them (e.g. cell cycle, homeostasis). Regarding the activity of B cells, our approach revealed that some antibody-related mechanisms remain activated until day 60 post infection. The Matlab codes for implementing OLYMPUS, as well as example datasets, are freely accessible via the Web (http://biosignal.med.upatras.gr/wordpress/biosignal/).