Time-course Gene Expression Data Research Articles

Network Activity Evaluation Reveals Significant Gene Regulatory Architectures during SARS-CoV-2 Viral Infection from Dynamic scRNA-seq Data.

The key to understand COVID-19 caused by SARS-CoV-2, which has caused massive deaths worldwide, is to reveal the gene activities at molecular level. Single-cell RNA-sequencing (scRNA-seq) technology allows us to capture gene expression at high resolution, thereby delineating cell-specific gene regulatory network (GRN). Network activity refers to the degree of consistency between GRN architectures and gene expression profiles in a specific condition or cellular microenvironment. Currently, numerous experimentally determined molecular interactions, including regulatory relationships closely related to SARS-CoV-2 infection, are documented in knowledge-bases. However, GRN activity is closely related to the cell dynamic environment and the heterogeneity of cell clusters. Therefore, to evaluate the consistency of GRN with gene expression profiles, we propose a single-cell Network Activity Evaluation framework, called scNAE. First, scNAE performs ODE modeling of time-course gene expression data. Then, the loss function with regularization penalty terms is constructed for formulating GRN inference rules from transcriptomic data. Furthermore, we have devised a rapid-convergence alternating direction method of multipliers to solve the regularized and constrained programs. Finally, an empirical P-value is derived based on a permutation statistical testing procedure to quantify the likelihood significance of the network matching with the data. The efficiency and advantage of scNAE have also been demonstrated by extensive numerical experiments, which can clearly depict the dynamic responses underlying GRN architectures triggered by the infection of SARS-CoV-2 in cells. The code and data of scNAE are available at https://github.com/zpliulab/scNAE.

Chromosome-level genome assembly of Quercus variabilis provides insights into the molecular mechanism of cork thickness

Quercus variabilis is a deciduous woody species with high ecological and economic value, and is a major source of cork in East Asia. Cork from thick softwood sheets have higher commercial value than those from thin sheets. It is extremely difficult to genetically improve Q. variabilis to produce high quality softwood due to the lack of genomic information. Here, we present a high-quality chromosomal genome assembly for Q. variabilis with length of 791,89 Mb and 54,606 predicted genes. Comparative analysis of protein sequences of Q. variabilis with 11 other species revealed that specific and expanded gene families were significantly enriched in the "fatty acid biosynthesis" pathway in Q. variabilis, which may contribute to the formation of its unique cork. Based on weighted correlation network analysis of time-course (i.e., five important developmental ages) gene expression data in thick-cork versus thin-cork genotypes of Q. variabilis, we identified one co-expression gene module associated with the thick-cork trait. Within this co-expression gene module, 10 hub genes were associated with suberin biosynthesis. Furthermore, we identified a total of 198 suberin biosynthesis-related new candidate genes that were up-regulated in trees with a thick cork layer relative to those with a thin cork layer. Also, we found that some genes related to cell expansion and cell division were highly expressed in trees with a thick cork layer. Collectively, our results revealed that two metabolic pathways (i.e., suberin biosynthesis, fatty acid biosynthesis), along with other genes involved in cell expansion, cell division, and transcriptional regulation, were associated with the thick-cork trait in Q. variabilis, providing insights into the molecular basis of cork development and knowledge for informing genetic improvement of cork thickness in Q. variabilis and closely related species.

Estimation of Ordinary Differential Equation Models for Gene Regulatory Networks Through Data Cloning.

Ordinary differential equations (ODEs) are widely used for elucidating dynamic processes in various fields. One of the applications of ODEs is to describe dynamics of gene regulatory networks (GRNs), which is a critical step in understanding disease mechanisms. However, estimation of ODE models for GRNs is challenging because of inflexibility of the model and noisy data with complex error structures such as heteroscedasticity, correlations between genes, and time dependency. In addition, either a likelihood or Bayesian approach is commonly used for estimation of ODE models, but both approaches have benefits and drawbacks in their own right. Data cloning is a maximum likelihood (ML) estimation method through the Bayesian framework. Since it works in the Bayesian framework, it is free from local optimum problems that are common drawbacks of ML methods. Also, its inference is invariant for the selection of prior distributions, which is a major issue in Bayesian methods. This study proposes an estimation method of ODE models for GRNs through data cloning. The proposed method is demonstrated through simulation and it is applied to real gene expression time-course data.

DynDeepDRIM: a dynamic deep learning model to infer direct regulatory interactions using time-course single-cell gene expression data.

Time-course single-cell RNA sequencing (scRNA-seq) data have been widely used to explore dynamic changes in gene expression of transcription factors (TFs) and their target genes. This information is useful to reconstruct cell-type-specific gene regulatory networks (GRNs). However, the existing tools are commonly designed to analyze either time-course bulk gene expression data or static scRNA-seq data via pseudo-time cell ordering. A few methods successfully utilize the information from multiple time points while also considering the characteristics of scRNA-seq data. We proposed dynDeepDRIM, a novel deep learning model to reconstruct GRNs using time-course scRNA-seq data. It represents the joint expression of a gene pair as an image and utilizes the image of the target TF-gene pair and the ones of the potential neighbors to reconstruct GRNs from time-course scRNA-seq data. dynDeepDRIM can effectively remove the transitive TF-gene interactions by considering neighborhood context and model the gene expression dynamics using high-dimensional tensors. We compared dynDeepDRIM with six GRN reconstruction methods on both simulation and four real time-course scRNA-seq data. dynDeepDRIM achieved substantially better performance than the other methods in inferring TF-gene interactions and eliminated the false positives effectively. We also applied dynDeepDRIM to annotate gene functions and found it achieved evidently better performance than the other tools due to considering the neighbor genes.

APRIL Drives a Coordinated but Diverse Response as a Foundation for Plasma Cell Longevity.

Ab-secreting cells survive in niche microenvironments, but cellular responses driven by particular niche signals are incompletely defined. The TNF superfamily member a proliferation-inducing ligand (APRIL) can support the maturation of transitory plasmablasts into long-lived plasma cells. In this study, we explore the biological programs established by APRIL in human plasmablasts. Under conditions allowing the maturation of ex vivo- or in vitro-generated plasmablasts, we find that APRIL drives activation of ERK, p38, and JNK, accompanied by a classical NF-κB response and activation of the AKT/FOXO1 pathway. Time-course gene expression data resolve coordinated transcriptional responses propagated via immediate early genes and NF-κB targets and converging onto modules of genes enriched for MYC targets and metabolism/cell growth-related pathways. This response is shared between APRIL and an alternate TNF superfamily member CD40L but is not a feature of alternative niche signals delivered by IFN-α or SDF1. However, APRIL and CD40L responses also diverge. CD40L drives expression of genes related to the activated B cell state whereas APRIL does not. Thus, APRIL establishes a broad foundation for plasma cell longevity with features of cellular refueling while being uncoupled from support of the B cell state.

Generating dynamic gene expression patterns without the need for regulatory circuits.

Synthetic biology has successfully advanced our ability to design and implement complex, time-varying genetic circuits to control the expression of recombinant proteins. However, these circuits typically require the production of regulatory genes whose only purpose is to coordinate expression of other genes. When designing very small genetic constructs, such as viral genomes, we may want to avoid introducing such auxiliary gene products while nevertheless encoding complex expression dynamics. To this end, here we demonstrate that varying only the placement and strengths of promoters, terminators, and RNase cleavage sites in a computational model of a bacteriophage genome is sufficient to achieve solutions to a variety of basic gene expression patterns. We discover these genetic solutions by computationally evolving genomes to reproduce desired gene expression time-course data. Our approach shows that non-trivial patterns can be evolved, including patterns where the relative ordering of genes by abundance changes over time. We find that some patterns are easier to evolve than others, and comparable expression patterns can be achieved via different genetic architectures. Our work opens up a novel avenue to genome engineering via fine-tuning the balance of gene expression and gene degradation rates.

Biclustering analysis of functionals via penalized fusion

In biomedical data analysis, clustering is commonly conducted. Biclustering analysis conducts clustering in both the sample and covariate dimensions and can more comprehensively describe data heterogeneity. In most of the existing biclustering analyses, scalar measurements are considered. In this study, motivated by time-course gene expression data and other examples, we take the “natural next step” and consider the biclustering analysis of functionals under which, for each covariate of each sample, a function (to be exact, its values at discrete measurement points) is present. We develop a doubly penalized fusion approach, which includes a smoothness penalty for estimating functionals and, more importantly, a fusion penalty for clustering. Statistical properties are rigorously established, providing the proposed approach a strong ground. We also develop an effective ADMM algorithm and accompanying R code. Numerical analysis, including simulations, comparisons, and the analysis of two time-course gene expression data, demonstrates the practical effectiveness of the proposed approach.

A novel approach for the analysis of time-course gene expression data based on computing with words

In this paper, a novel approach is proposed for the analysis of time-course gene expression data based on the path-breaking work of Zadeh, Computing with Words. This method can automatically discover the patterns of temporal gene expression profile in terms of two distinguishing descriptions: linguistic description that is understandable and interpretable for human inference; and type-2 fuzzy description that is suitable for robust machine inference in the presence of uncertainty. In contrast to conventional static data mining methods which focus on the steady-state gene expression levels, the proposed scheme is a new time-series pattern mining technique for dynamical modeling of gene expression. To evaluate the performance of this paradigm, it is applied to a case study dataset from Gene Expression Omnibus (GEO) which includes the temporal transcriptional profile of human colon cancer cells. The goal is to investigate the pharmacodynamics of two anticancer drugs. The transient and steady-state analysis of the transcriptional response clearly demonstrates the ability of the proposed approach to reveal the pharmacodynamical effects of drug type and dosage on the expression of genes.

Integrative Modeling of Gene Expression and Metabolic Networks of Arabidopsis Embryos for Identification of Seed Oil Causal Genes.

Fatty acids in crop seeds are a major source for both vegetable oils and industrial applications. Genetic improvement of fatty acid composition and oil content is critical to meet the current and future demands of plant-based renewable seed oils. Addressing this challenge can be approached by network modeling to capture key contributors of seed metabolism and to identify underpinning genetic targets for engineering the traits associated with seed oil composition and content. Here, we present a dynamic model, using an Ordinary Differential Equations model and integrated time-course gene expression data, to describe metabolic networks during Arabidopsis thaliana seed development. Through in silico perturbation of genes, targets were predicted in seed oil traits. Validation and supporting evidence were obtained for several of these predictions using published reports in the scientific literature. Furthermore, we investigated two predicted targets using omics datasets for both gene expression and metabolites from the seed embryo, and demonstrated the applicability of this network-based model. This work highlights that integration of dynamic gene expression atlases generates informative models which can be explored to dissect metabolic pathways and lead to the identification of causal genes associated with seed oil traits.

Inferring latent temporal progression and regulatory networks from cross-sectional transcriptomic data of cancer samples.

Unraveling molecular regulatory networks underlying disease progression is critically important for understanding disease mechanisms and identifying drug targets. The existing methods for inferring gene regulatory networks (GRNs) rely mainly on time-course gene expression data. However, most available omics data from cross-sectional studies of cancer patients often lack sufficient temporal information, leading to a key challenge for GRN inference. Through quantifying the latent progression using random walks-based manifold distance, we propose a latent-temporal progression-based Bayesian method, PROB, for inferring GRNs from the cross-sectional transcriptomic data of tumor samples. The robustness of PROB to the measurement variabilities in the data is mathematically proved and numerically verified. Performance evaluation on real data indicates that PROB outperforms other methods in both pseudotime inference and GRN inference. Applications to bladder cancer and breast cancer demonstrate that our method is effective to identify key regulators of cancer progression or drug targets. The identified ACSS1 is experimentally validated to promote epithelial-to-mesenchymal transition of bladder cancer cells, and the predicted FOXM1-targets interactions are verified and are predictive of relapse in breast cancer. Our study suggests new effective ways to clinical transcriptomic data modeling for characterizing cancer progression and facilitates the translation of regulatory network-based approaches into precision medicine.

Bayesian Functional Mixed-effects Models with Grouped Smoothness for Analyzing Time-course Gene Expression Data

Objective: As a result of the development of microarray technologies, gene expression levels of thousands of genes involved in a given biological process can be measured simultaneously, and it is important to study their temporal behavior to understand their mechanisms. Since the dependence between gene expression levels over time for a given gene is often too complicated to model parametrically, sparse functional data analysis has received an increasing amount of attention for analyzing such data. Methods: We propose a new functional mixed-effects model for analyzing time-course gene expression data. Specifically, the model groups individual functions with heterogeneous smoothness. The proposed method utilizes the mixed-effects model representation of penalized splines for both the mean function and the individual functions. Given noninformative or weakly informative priors, Bayesian inference on the proposed models was developed, and Bayesian computation was implemented by using Markov chain Monte Carlo methods. Results: The performance of our new model was studied by two simulation studies and illustrated using a yeast cell cycle gene expression dataset. Simulation results suggest that our proposed methods can outperform the previously used methods in terms of the mean integrated squared error. The yeast gene expression data application suggests that the proposed model with two latent groups should be used on this dataset. Conclusion: The new Bayesian functional mixed-effects model that assumes multiple groups of functions with different smoothing parameters provides an enhanced approach to analyzing timecourse gene expression data.

A pipeline to analyse time-course gene expression data

The phenotypic diversity of cells is governed by a complex equilibrium between their genetic identity and their environmental interactions: Understanding the dynamics of gene expression is a fundamental question of biology. However, analysing time-course transcriptomic data raises unique challenging statistical and computational questions, requiring the development of novel methods and software. This workflow provides a step-by-step tutorial of the methodology used to analyse time-course data: (1) quality control and normalization of the dataset; (2) differential expression analysis using functional data analysis; (3) clustering of time-course data; (4) interpreting clusters with GO term and KEGG pathway enrichment analysis. As a case study, we apply this workflow to time-course transcriptomic data from mice exposed to four strains of influenza to showcase every step of the pipeline.

An informatics research platform to make public gene expression time-course datasets reusable for more scientific discoveries.

The exponential growth of genomic/genetic data in the era of Big Data demands new solutions for making these data findable, accessible, interoperable and reusable. In this article, we present a web-based platform named Gene Expression Time-Course Research (GETc) Platform that enables the discovery and visualization of time-course gene expression data and analytical results from the NIH/NCBI-sponsored Gene Expression Omnibus (GEO). The analytical results are produced from an analytic pipeline based on the ordinary differential equation model. Furthermore, in order to extract scientific insights from these results and disseminate the scientific findings, close and efficient collaborations between domain-specific experts from biomedical and scientific fields and data scientists is required. Therefore, GETc provides several recommendation functions and tools to facilitate effective collaborations. GETc platform is a very useful tool for researchers from the biomedical genomics community to present and communicate large numbers of analysis results from GEO. It is generalizable and broadly applicable across different biomedical research areas. GETc is a user-friendly and efficient web-based platform freely accessible at http://genestudy.org/

A recursive framework for predicting the time-course of drug sensitivity

The biological processes involved in a drug’s mechanisms of action are oftentimes dynamic, complex and difficult to discern. Time-course gene expression data is a rich source of information that can be used to unravel these complex processes, identify biomarkers of drug sensitivity and predict the response to a drug. However, the majority of previous work has not fully utilized this temporal dimension. In these studies, the gene expression data is either considered at one time-point (before the administration of the drug) or two time-points (before and after the administration of the drug). This is clearly inadequate in modeling dynamic gene–drug interactions, especially for applications such as long-term drug therapy. In this work, we present a novel REcursive Prediction (REP) framework for drug response prediction by taking advantage of time-course gene expression data. Our goal is to predict drug response values at every stage of a long-term treatment, given the expression levels of genes collected in the previous time-points. To this end, REP employs a built-in recursive structure that exploits the intrinsic time-course nature of the data and integrates past values of drug responses for subsequent predictions. It also incorporates tensor completion that can not only alleviate the impact of noise and missing data, but also predict unseen gene expression levels (GEXs). These advantages enable REP to estimate drug response at any stage of a given treatment from some GEXs measured in the beginning of the treatment. Extensive experiments on two datasets corresponding to multiple sclerosis patients treated with interferon are included to showcase the effectiveness of REP.

Can ODE gene regulatory models neglect time lag or measurement scaling?

Many ordinary differential equation (ODE) models have been introduced to replace linear regression models for inferring gene regulatory relationships from time-course gene expression data. But, since the observed data are usually not direct measurements of the gene products or there is an unknown time lag in gene regulation, it is problematic to directly apply traditional ODE models or linear regression models. We introduce a lagged ODE model to infer lagged gene regulatory relationships from time-course measurements, which are modeled as linear transformation of the gene products. A time-course microarray dataset from a yeast cell-cycle study is used for simulation assessment of the methods and real data analysis. The results show that our method, by considering both time lag and measurement scaling, performs much better than other linear and ODE models. It indicates the necessity of explicitly modeling the time lag and measurement scaling in ODE gene regulatory models. R code is available at https://www.sta.cuhk.edu.hk/xfan/share/lagODE.zip.

Core transcriptional regulatory circuits in prion diseases

Complex diseases involve dynamic perturbations of pathophysiological processes during disease progression. Transcriptional programs underlying such perturbations are unknown in many diseases. Here, we present core transcriptional regulatory circuits underlying early and late perturbations in prion disease. We first identified cellular processes perturbed early and late using time-course gene expression data from three prion-infected mouse strains. We then built a transcriptional regulatory network (TRN) describing regulation of early and late processes. We found over-represented feed-forward loops (FFLs) comprising transcription factor (TF) pairs and target genes in the TRN. Using gene expression data of brain cell types, we further selected active FFLs where TF pairs and target genes were expressed in the same cell type and showed correlated temporal expression changes in the brain. We finally determined core transcriptional regulatory circuits by combining these active FFLs. These circuits provide insights into transcriptional programs for early and late pathophysiological processes in prion disease.

TGCnA: temporal gene coexpression network analysis using a low-rank plus sparse framework

ABSTRACTVarious gene network models with distinct physical nature have been widely used in biological studies. For temporal transcriptomic studies, the current dynamic models either ignore the temporal variation in the network structure or fail to scale up to a large number of genes due to severe computational bottlenecks and sample size limitation. Although the correlation-based gene networks are computationally affordable, they have limitations after being applied to gene expression time-course data. We proposed Temporal Gene Coexpression Network Analysis (TGCnA) framework for the transcriptomic time-course data. The mathematical nature of TGCnA is the joint modeling of multiple covariance matrices across time points using a ‘low-rank plus sparse’ framework, in which the network similarity across time points is explicitly modeled in the low-rank component. We demonstrated the advantage of TGCnA in covariance matrix estimation and gene module discovery using both simulation data and real transcriptomic data. The code is available at https://github.com/QiZhangStat/TGCnA.

Comparative Genomics Analysis Provides New Insight Into Molecular Basis of Stomatal Movement in Kalanchoë fedtschenkoi.

CO2 uptake and water loss in plants are regulated by microscopic pores on the surface of leaves, called stomata. This enablement of gas exchange by the opening and closing of stomata is one of the most essential processes in plant photosynthesis and transpiration, affecting water-use efficiency (WUE) and thus drought susceptibility. In plant species with crassulacean acid metabolism (CAM) photosynthesis, diel stomatal movement pattern is inverted relative to C3 and C4 photosynthesis species, resulting in much higher WUE and drought tolerance. However, little is known about the molecular basis of stomatal movement in CAM species. The goal of this study is to identify candidate genes that could play a role in stomatal movement in an obligate CAM species, Kalanchoë fedtschenkoi. By way of a text-mining approach, proteins were identified in various plant species, spanning C3, C4, and CAM photosynthetic types, which are orthologous to proteins known to be involved in stomatal movement. A comparative analysis of diel time-course gene expression data was performed between K. fedtschenkoi and two C3 species (i.e., Arabidopsis thaliana and Solanum lycopersicum) to identify differential gene expression between the dusk and dawn phases of the 24-h cycle. A rescheduled catalase gene known to be involved in stomatal movement was identified, suggesting a role for H2O2 in CAM-like stomatal movement. Overall, these results provide new insights into the molecular regulation of stomatal movement in CAM plants, facilitating genetic improvement of drought resistance in agricultural crops through manipulation of stomata-related genes.

Identifying Windows of Susceptibility by Temporal Gene Analysis

Increased understanding of developmental disorders of the brain has shown that genetic mutations, environmental toxins and biological insults typically act during developmental windows of susceptibility. Identifying these vulnerable periods is a necessary and vital step for safeguarding women and their fetuses against disease causing agents during pregnancy and for developing timely interventions and treatments for neurodevelopmental disorders. We analyzed developmental time-course gene expression data derived from human pluripotent stem cells, with disease association, pathway, and protein interaction databases to identify windows of disease susceptibility during development and the time periods for productive interventions. The results are displayed as interactive Susceptibility Windows Ontological Transcriptome (SWOT) Clocks illustrating disease susceptibility over developmental time. Using this method, we determine the likely windows of susceptibility for multiple neurological disorders using known disease associated genes and genes derived from RNA-sequencing studies including autism spectrum disorder, schizophrenia, and Zika virus induced microcephaly. SWOT clocks provide a valuable tool for integrating data from multiple databases in a developmental context with data generated from next-generation sequencing to help identify windows of susceptibility.

Conservation and Diversification of Circadian Rhythmicity Between a Model Crassulacean Acid Metabolism Plant Kalanchoë fedtschenkoi and a Model C3 Photosynthesis Plant Arabidopsis thaliana.

Crassulacean acid metabolism (CAM) improves photosynthetic efficiency under limited water availability relative to C3 photosynthesis. It is widely accepted that CAM plants have evolved from C3 plants and it is hypothesized that CAM is under the control of the internal circadian clock. However, the role that the circadian clock plays in the evolution of CAM is not well understood. To identify the molecular basis of circadian control over CAM evolution, rhythmic gene sets were identified in a CAM model plant species (Kalanchoë fedtschenkoi) and a C3 model plant species (Arabidopsis thaliana) through analysis of diel time-course gene expression data using multiple periodicity detection algorithms. Based on protein sequences, ortholog groups were constructed containing genes from each of these two species. The ortholog groups were categorized into five gene sets based on conservation and diversification of rhythmic gene expression. Interestingly, minimal functional overlap was observed when comparing the rhythmic gene sets of each species. Specifcally, metabolic processes were enriched in the gene set under circadian control in K. fedtschenkoi and numerous genes were found to have retained or gained rhythmic expression in K. fedtsechenkoi. Additonally, several rhythmic orthologs, including CAM-related orthologs, displayed phase shifts between species. Results of this analysis point to several mechanisms by which the circadian clock plays a role in the evolution of CAM. These genes provide a set of testable hypotheses for future experiments.