Studying and modelling dynamic biological processes using time-series gene expression data

GeTeSEPdb: A comprehensive database and online tool for the identification and analysis of gene profiles with temporal-specific expression patterns

Gene expression is dynamic and varies at different stages of processes. The identification of gene profiles with temporal-specific expression patterns can provide valuable insights into ongoing biological processes, such as the cell cycle, cell development, circadian rhythms, or responses to external stimuli such as drug treatments or viral infections. However, currently, no database defines, identifies or archives gene profiles with temporal-specific expression patterns. Here, using a high-throughput regression analysis approach, eight linear and nonlinear parametric models were fitted to gene expression profiles from time-series experiments to identify eight types of gene profiles with temporal-specific expression patterns. We curated 2684 time-series transcriptome datasets and identified 2644,370 gene profiles exhibiting temporal-specific expression patterns. The results were stored in the database GeTeSEPdb (gene profiles with temporal-specific expression patterns database, http://www.inbirg.com/GeTeSEPdb/). Moreover, we implemented an online tool to identify gene profiles with temporal-specific expression patterns from user-submitted data. In summary, GeTeSEPdb is a comprehensive web service that can be used to identify and analyse gene profiles with temporal-specific expression patterns. This approach facilitates the exploration of transcriptional changes and temporal patterns of responses. We firmly believe that GeTeSEPdb will become a valuable resource for biologists and bioinformaticians.

Integrating patients in time series clinical transcriptomics data

MOTIVATION

Analysis of time series transcriptomics data from clinical trials is challenging. Such studies usually profile very few time points from several individuals with varying response patterns and dynamics. Current methods for these datasets are mainly based on linear, global orderings using visit times which do not account for the varying response rates and subgroups within a patient cohort.

RESULTS

We developed a new method that utilizes multi-commodity flow algorithms for trajectory inference in large scale clinical studies. Recovered trajectories satisfy individual-based timing restrictions while integrating data from multiple patients. Testing the method on multiple drug datasets demonstrated an improved performance compared to prior approaches suggested for this task, while identifying novel disease subtypes that correspond to heterogeneous patient response patterns.

AVAILABILITY AND IMPLEMENTATION

The source code and instructions to download the data have been deposited on GitHub at https://github.com/euxhenh/Truffle.

Inferring upstream regulatory genes of FOXP3 in human regulatory T cells from time-series transcriptomic data

The discovery of upstream regulatory genes of a gene of interest still remains challenging. Here we applied a scalable computational method to unbiasedly predict candidate regulatory genes of critical transcription factors by searching the whole genome. We illustrated our approach with a case study on the master regulator FOXP3 of human primary regulatory T cells (Tregs). While target genes of FOXP3 have been identified, its upstream regulatory machinery still remains elusive. Our methodology selected five top-ranked candidates that were tested via proof-of-concept experiments. Following knockdown, three out of five candidates showed significant effects on the mRNA expression of FOXP3 across multiple donors. This provides insights into the regulatory mechanisms modulating FOXP3 transcriptional expression in Tregs. Overall, at the genome level this represents a high level of accuracy in predicting upstream regulatory genes of key genes of interest.

Unveiling gene regulatory networks during cellular state transitions without linkage across time points

Time-stamped cross-sectional data, which lack linkage across time points, are commonly generated in single-cell transcriptional profiling. Many previous methods for inferring gene regulatory networks (GRNs) driving cell-state transitions relied on constructing single-cell temporal ordering. Introducing COSLIR (COvariance restricted Sparse LInear Regression), we presented a direct approach to reconstructing GRNs that govern cell-state transitions, utilizing only the first and second moments of samples between two consecutive time points. Simulations validated COSLIR's perfect accuracy in the oracle case and demonstrated its robust performance in real-world scenarios. When applied to single-cell RT-PCR and RNAseq datasets in developmental biology, COSLIR competed favorably with existing methods. Notably, its running time remained nearly independent of the number of cells. Therefore, COSLIR emerges as a promising addition to GRN reconstruction methods under cell-state transitions, bypassing the single-cell temporal ordering to enhance accuracy and efficiency in single-cell transcriptional profiling.

LogicGep: Boolean networks inference using symbolic regression from time-series transcriptomic profiling data

Reconstructing the topology of gene regulatory network from gene expression data has been extensively studied. With the abundance functional transcriptomic data available, it is now feasible to systematically decipher regulatory interaction dynamics in a logic form such as a Boolean network (BN) framework, which qualitatively indicates how multiple regulators aggregated to affect a common target gene. However, inferring both the network topology and gene interaction dynamics simultaneously is still a challenging problem since gene expression data are typically noisy and data discretization is prone to information loss. We propose a new method for BN inference from time-series transcriptional profiles, called LogicGep. LogicGep formulates the identification of Boolean functions as a symbolic regression problem that learns the Boolean function expression and solve it efficiently through multi-objective optimization using an improved gene expression programming algorithm. To avoid overly emphasizing dynamic characteristics at the expense of topology structure ones, as traditional methods often do, a set of promising Boolean formulas for each target gene is evolved firstly, and a feed-forward neural network trained with continuous expression data is subsequently employed to pick out the final solution. We validated the efficacy of LogicGep using multiple datasets including both synthetic and real-world experimental data. The results elucidate that LogicGep adeptly infers accurate BN models, outperforming other representative BN inference algorithms in both network topology reconstruction and the identification of Boolean functions. Moreover, the execution of LogicGep is hundreds of times faster than other methods, especially in the case of large network inference.

A Causal Regulation Modeling Algorithm for Temporal Events with Application to Escherichia coli's Aerobic to Anaerobic Transition

Time-series experiments are crucial for understanding the transient and dynamic nature of biological phenomena. These experiments, leveraging advanced classification and clustering algorithms, allow for a deep dive into the cellular processes. However, while these approaches effectively identify patterns and trends within data, they often need to improve in elucidating the causal mechanisms behind these changes. Building on this foundation, our study introduces a novel algorithm for temporal causal signaling modeling, integrating established knowledge networks with sequential gene expression data to elucidate signal transduction pathways over time. Focusing on Escherichia coli's (E. coli) aerobic to anaerobic transition (AAT), this research marks a significant leap in understanding the organism's metabolic shifts. By applying our algorithm to a comprehensive E. coli regulatory network and a time-series microarray dataset, we constructed the cross-time point core signaling and regulatory processes of E. coli's AAT. Through gene expression analysis, we validated the primary regulatory interactions governing this process. We identified a novel regulatory scheme wherein environmentally responsive genes, soxR and oxyR, activate fur, modulating the nitrogen metabolism regulators fnr and nac. This regulatory cascade controls the stress regulators ompR and lrhA, ultimately affecting the cell motility gene flhD, unveiling a novel regulatory axis that elucidates the complex regulatory dynamics during the AAT process. Our approach, merging empirical data with prior knowledge, represents a significant advance in modeling cellular signaling processes, offering a deeper understanding of microbial physiology and its applications in biotechnology.

MultiRNAflow: integrated analysis of temporal RNA-seq data with multiple biological conditions

MOTIVATION

The dynamic transcriptional mechanisms that govern eukaryotic cell function can now be analyzed by RNA sequencing. However, the packages currently available for the analysis of raw sequencing data do not provide automatic analysis of complex experimental designs with multiple biological conditions and multiple analysis time-points.

RESULTS

The MultiRNAflow suite combines several packages in a unified framework allowing exploratory and supervised statistical analyses of temporal data for multiple biological conditions.

AVAILABILITY AND IMPLEMENTATION

The R package MultiRNAflow is freely available on Bioconductor (https://bioconductor.org/packages/MultiRNAflow/), and the latest version of the source code is available on a GitHub repository (https://github.com/loubator/MultiRNAflow).

Time-/dose- series transcriptome data analysis and traditional Chinese medicine treatment of pneumoconiosis

Pneumoconiosis' pathogenesis is still unclear and specific drugs for its treatment are lacking. Analysis of series transcriptome data often uses a single comparison method, and there are few reports on using such data to predict the treatment of pneumoconiosis with traditional Chinese medicine (TCM). Here, we proposed a new method for analyzing series transcriptomic data, series difference analysis (SDA), and applied it to pneumoconiosis. By comparison with 5 gene sets including existing pneumoconiosis-related genes and gene set functional enrichment analysis, we demonstrated that the new method was not inferior to two existing traditional analysis methods. Furthermore, based on the TCM-drug target interaction network, we predicted the TCM corresponding to the common pneumoconiosis-related genes obtained by multiple methods, and combined them with the high-frequency TCM for its treatment obtained through literature mining to form a new TCM formula for it. After feeding it to pneumoconiosis modeling mice for two months, compared with the untreated group, the coat color, mental state and tissue sections of the mice in the treated group were markedly improved, indicating that the new TCM formula has a certain efficacy. Our study provides new insights into method development for series transcriptomic data analysis and treatment of pneumoconiosis.

Modeling gene expression cascades during cell state transitions

During cellular processes such as differentiation or response to external stimuli, cells exhibit dynamic changes in their gene expression profiles. Single-cell RNA sequencing (scRNA-seq) can be used to investigate these dynamic changes. To this end, cells are typically ordered along a pseudotemporal trajectory which recapitulates the progression of cells as they transition from one cell state to another. We infer transcriptional dynamics by modeling the gene expression profiles in pseudotemporally ordered cells using a Bayesian inference approach. This enables ordering genes along transcriptional cascades, estimating differences in the timing of gene expression dynamics, and deducing regulatory gene interactions. Here, we apply this approach to scRNA-seq datasets derived from mouse embryonic forebrain and pancreas samples. This analysis demonstrates the utility of the method to derive the ordering of gene dynamics and regulatory relationships critical for proper cellular differentiation and maturation across a variety of developmental contexts.

Powerful and accurate detection of temporal gene expression patterns from multi-sample multi-stage single-cell transcriptomics data with TDEseq

We present a non-parametric statistical method called TDEseq that takes full advantage of smoothing splines basis functions to account for the dependence of multiple time points in scRNA-seq studies, and uses hierarchical structure linear additive mixed models to model the correlated cells within an individual. As a result, TDEseq demonstrates powerful performance in identifying four potential temporal expression patterns within a specific cell type. Extensive simulation studies and the analysis of four published scRNA-seq datasets show that TDEseq can produce well-calibrated p-values and up to 20% power gain over the existing methods for detecting temporal gene expression patterns.

Predicting gene-level sensitivity to JAK-STAT signaling perturbation using a mechanistic-to-machine learning framework

The Janus kinase (JAK)-signal transducer and activator of transcription (STAT) pathway integrates complex cytokine signals via a limited number of molecular components, inspiring numerous efforts to clarify the diversity and specificity of STAT transcription factor function. We developed a computational framework to make global cytokine-induced gene predictions from STAT phosphorylation dynamics, modeling macrophage responses to interleukin (IL)-6 and IL-10, which signal through common STATs, but with distinct temporal dynamics and contrasting functions. Our mechanistic-to-machine learning model identified cytokine-specific genes associated with late pSTAT3 time frames and a preferential pSTAT1 reduction upon JAK2 inhibition. We predicted and validated the impact of JAK2 inhibition on gene expression, identifying genes that were sensitive or insensitive to JAK2 variation. Thus, we successfully linked STAT signaling dynamics to gene expression to support future efforts targeting pathology-associated STAT-driven gene sets. This serves as a first step in developing multi-level prediction models to understand and perturb gene expression outputs from signaling systems. A record of this paper's transparent peer review process is included in the supplemental information.

Inference of Dynamic Growth Regulatory Network in Cancer Using High-Throughput Transcriptomic Data

Growth is regulated by gene expression variation at different developmental stages of biological processes such as cell differentiation, disease progression, or drug response. In cancer, a stage-specific regulatory model constructed to infer the dynamic expression changes in genes contributing to tissue growth or proliferation is referred as a dynamic growth regulatory network (dGRN). Over the past decade, gene expression data has been widely used for reconstructing dGRN by computing correlations between the differentially expressed genes (DEGs). A wide variety of pipelines are available to construct the GRNs using DEGs and the choice of a particular method or tool depends on the nature of the study. In this protocol, we have outlined a step-by-step guide for the analysis of DEGs using RNA-Seq data, beginning from data acquisition, pre-processing, mapping to reference genome, and construction of a correlation-based co-expression network to further downstream analysis. We have also outlined the steps for the inclusion of publicly available interaction/regulation information into the dGRN followed by relevant topological inferences. This tutorial has been designed in a way that early researchers can refer to for an easy and comprehensive glimpse of methodologies used in the inference of dGRN using transcriptomics data.

Genome-Wide Analysis of the MADS-Box Gene Family in Hibiscus syriacus and Their Role in Floral Organ Development

Hibiscus syriacus belongs to the Malvaceae family, and is a plant with medicinal, edible, and greening values. MADS-box transcription factor is a large family of regulatory factors involved in a variety of biological processes in plants. Here, we performed a genome-wide characterization of MADS-box proteins in H. syriacus and investigated gene structure, phylogenetics, cis-acting elements, three-dimensional structure, gene expression, and protein interaction to identify candidate MADS-box genes that mediate petal developmental regulation in H. syriacus. A total of 163 candidate MADS-box genes were found and classified into type I (Mα, Mβ, and Mγ) and type II (MIKC and Mδ). Analysis of cis-acting elements in the promoter region showed that most elements were correlated to plant hormones. The analysis of nine HsMADS expressions of two different H. syriacus cultivars showed that they were differentially expressed between two type flowers. The analysis of protein interaction networks also indicated that MADS proteins played a crucial role in floral organ identification, inflorescence and fruit development, and flowering time. This research is the first to analyze the MADS-box family of H. syriacus and provides an important reference for further study of the biological functions of the MADS-box, especially in flower organ development.

Dynamic changes in transcriptome during orthodontic tooth movement

OBJECTIVES

The objective of this study was to determine global changes in gene expression with next generation sequencing (NGS) in order to assess the biological effects of orthodontic tooth movement (OTM) on alveolar bone in a rat model.

MATERIALS AND METHODS

Thirty-five Wistar rats (age 14 weeks) were used in the study. The OTM was performed using closed coil Nickel-Titanium spring to apply a mesial force on maxillary first molars of 8-10 g. Three hours, 1, 3, 7 and 14 days after the placement of the appliance, rats were killed at each time point respectively. The alveolar bone, around left maxillary first molar, were excised on compression side. The samples were immediately frozen in liquid nitrogen for subsequent RNA extraction. Total RNA samples were prepared for mRNA sequencing using the Illumina kit. RNA-Seq reads were aligned to the rat genomes using the STAR Aligner and bioinformatic analysis was performed.

RESULTS

A total of 18 192 genes were determined. Day 1 has the highest number of differentially expressed genes (DEGs) observed with more upregulated than downregulated genes. A total of 2719 DEGs were identified to use as input for the algorithm. Six distinct clusters of temporal patterns were observed representing proteins that were differentially regulated indicating different expression kinetics. Principal component analysis (PCA) showed distinct clustering by time points and days 3, 7 and 14 share similar gene expression pattern.

CONCLUSIONS

Distinct gene expression pattern was observed at different time points studied. Hypoxia, inflammation and bone remodelling pathways are major mechanisms behind OTM.

Global quantitative understanding of non-equilibrium cell fate decision-making in response to pheromone

Cell-cycle arrest and polarized growth are commonly used to characterize the response of yeast to pheromone. However, the quantitative decision-making processes underlying time-dependent changes in cell fate remain unclear. In this study, we conducted single-cell level experiments to observe multidimensional responses, uncovering diverse fates of yeast cells. Multiple states are revealed, along with the kinetic switching rates and pathways among them, giving rise to a quantitative landscape of mating response. To quantify the experimentally observed cell fates, we developed a theoretical framework based on non-equilibrium landscape and flux theory. Additionally, we performed stochastic simulations of biochemical reactions to elucidate signal transduction and cell growth. Notably, our experimental findings have provided the first global quantitative evidence of the real-time synchronization between intracellular signaling, physiological growth, and morphological functions. These results validate the proposed underlying mechanism governing the emergence of multiple cell fate states. This study introduces an emerging mechanistic approach to understand non-equilibrium cell fate decision-making in response to pheromone.

Principles and challenges of modeling temporal and spatial omics data

Studies with temporal or spatial resolution are crucial to understand the molecular dynamics and spatial dependencies underlying a biological process or system. With advances in high-throughput omic technologies, time- and space-resolved molecular measurements at scale are increasingly accessible, providing new opportunities to study the role of timing or structure in a wide range of biological questions. At the same time, analyses of the data being generated in the context of spatiotemporal studies entail new challenges that need to be considered, including the need to account for temporal and spatial dependencies and compare them across different scales, biological samples or conditions. In this Review, we provide an overview of common principles and challenges in the analysis of temporal and spatial omics data. We discuss statistical concepts to model temporal and spatial dependencies and highlight opportunities for adapting existing analysis methods to data with temporal and spatial dimensions.

Mapping oto-pharyngeal development in a human inner ear organoid model

Inner ear development requires the coordination of cell types from distinct epithelial, mesenchymal and neuronal lineages. Although we have learned much from animal models, many details about human inner ear development remain elusive. We recently developed an in vitro model of human inner ear organogenesis using pluripotent stem cells in a 3D culture, fostering the growth of a sensorineural circuit, including hair cells and neurons. Despite previously characterizing some cell types, many remain undefined. This study aimed to chart the in vitro development timeline of the inner ear organoid to understand the mechanisms at play. Using single-cell RNA sequencing at ten stages during the first 36 days of differentiation, we tracked the evolution from pluripotency to various ear cell types after exposure to specific signaling modulators. Our findings showcase gene expression that influences differentiation, identifying a plethora of ectodermal and mesenchymal cell types. We also discern aspects of the organoid model consistent with in vivo development, while highlighting potential discrepancies. Our study establishes the Inner Ear Organoid Developmental Atlas (IODA), offering deeper insights into human biology and improving inner ear tissue differentiation.

Simultaneous estimation of gene regulatory network structure and RNA kinetics from single cell gene expression

Cells respond to environmental and developmental stimuli by remodeling their transcriptomes through regulation of both mRNA transcription and mRNA decay. A central goal of biology is identifying the global set of regulatory relationships between factors that control mRNA production and degradation and their target transcripts and construct a predictive model of gene expression. Regulatory relationships are typically identified using transcriptome measurements and causal inference algorithms. RNA kinetic parameters are determined experimentally by employing run-on or metabolic labeling (e.g. 4-thiouracil) methods that allow transcription and decay rates to be separately measured. Here, we develop a deep learning model, trained with single-cell RNA-seq data, that both infers causal regulatory relationships and estimates RNA kinetic parameters. The resulting in silico model predicts future gene expression states and can be perturbed to simulate the effect of transcription factor changes. We acquired model training data by sequencing the transcriptomes of 175,000 individual Saccharomyces cerevisiae cells that were subject to an external perturbation and continuously sampled over a one hour period. The rate of change for each transcript was calculated on a per-cell basis to estimate RNA velocity. We then trained a deep learning model with transcriptome and RNA velocity data to calculate time-dependent estimates of mRNA production and decay rates. By separating RNA velocity into transcription and decay rates, we show that rapamycin treatment causes existing ribosomal protein transcripts to be rapidly destabilized, while production of new transcripts gradually slows over the course of an hour. The neural network framework we present is designed to explicitly model causal regulatory relationships between transcription factors and their genes, and shows superior performance to existing models on the basis of recovery of known regulatory relationships. We validated the predictive power of the model by perturbing transcription factors in silico and comparing transcriptome-wide effects with experimental data. Our study represents the first step in constructing a complete, predictive, biophysical model of gene expression regulation.

Transcriptome profiling in rumen, reticulum, omasum, and abomasum tissues during the developmental transition of pre-ruminant to the ruminant in yaks

The development of the four stomachs of yak is closely related to its health and performance, however the underlying molecular mechanisms are largely unknown. Here, we systematically analyzed mRNAs of four stomachs in five growth time points [0 day, 20 days, 60 days, 15 months and 3 years (adult)] of yaks. Overall, the expression patterns of DEmRNAs were unique at 0 d, similar at 20 d and 60 d, and similar at 15 m and adult in four stomachs. The expression pattern in abomasum was markedly different from that in rumen, reticulum and omasum. Short Time-series Expression Miner (STEM) analysis demonstrated that multi-model spectra are drastically enriched over time in four stomachs. All the identified mRNAs in rumen, reticulum, omasum and abomasum were classified into 6, 4, 7, and 5 cluster profiles, respectively. Modules 9, 38, and 41 were the most significant three colored modules. By weighted gene co-expression network analysis (WGCNA), a total of 5,486 genes were categorized into 10 modules. CCKBR, KCNQ1, FER1L6, and A4GNT were the hub genes of the turquoise module, and PAK6, TRIM29, ADGRF4, TGM1, and TMEM79 were the hub genes of the blue module. Furthermore, functional KEGG enrichment analysis suggested that the turquoise module was involved in gastric acid secretion, sphingolipid metabolism, ether lipid metabolism, etc., and the blue module was enriched in pancreatic secretion, pantothenate and CoA biosynthesis, and starch and sucrose metabolism, etc. Our study aims to lay a molecular basis for the study of the physiological functions of rumen, reticulum, omasum and abomasum in yaks. It can further elucidate the important roles of these mRNAs in regulation of growth, development and metabolism in yaks, and to provide a theoretical basis for age-appropriate weaning and supplementary feeding in yaks.

Deep Neural Networks for Predicting Single-Cell Responses and Probability Landscapes

Engineering biology relies on the accurate prediction of cell responses. However, making these predictions is challenging for a variety of reasons, including the stochasticity of biochemical reactions, variability between cells, and incomplete information about underlying biological processes. Machine learning methods, which can model diverse input-output relationships without requiring a priori mechanistic knowledge, are an ideal tool for this task. For example, such approaches can be used to predict gene expression dynamics given time-series data of past expression history. To explore this application, we computationally simulated single-cell responses, incorporating different sources of noise and alternative genetic circuit designs. We showed that deep neural networks trained on these simulated data were able to correctly infer the underlying dynamics of a cell response even in the presence of measurement noise and stochasticity in the biochemical reactions. The training set size and the amount of past data provided as inputs both affected prediction quality, with cascaded genetic circuits that introduce delays requiring more past data. We also tested prediction performance on a bistable auto-activation circuit, finding that our initial method for predicting a single trajectory was fundamentally ill-suited for multimodal dynamics. To address this, we updated the network architecture to predict the entire distribution of future states, showing it could accurately predict bimodal expression distributions. Overall, these methods can be readily applied to the diverse prediction tasks necessary to predict and control a variety of biological circuits, a key aspect of many synthetic biology applications.

From time-series transcriptomics to gene regulatory networks: A review on inference methods

Inference of gene regulatory networks has been an active area of research for around 20 years, leading to the development of sophisticated inference algorithms based on a variety of assumptions and approaches. With the ever increasing demand for more accurate and powerful models, the inference problem remains of broad scientific interest. The abstract representation of biological systems through gene regulatory networks represents a powerful method to study such systems, encoding different amounts and types of information. In this review, we summarize the different types of inference algorithms specifically based on time-series transcriptomics, giving an overview of the main applications of gene regulatory networks in computational biology. This review is intended to give an updated reference of regulatory networks inference tools to biologists and researchers new to the topic and guide them in selecting the appropriate inference method that best fits their questions, aims, and experimental data.

TempShift Reveals the Sequential Development of Human Neocortex and Skewed Developmental Timing of Down Syndrome Brains

Development is a complex process involving precise regulation. Developmental regulation may vary in tissues and individuals, and is often altered in disorders. Currently, the regulation of developmental timing across neocortical areas and developmental changes in Down syndrome (DS) brains remain unclear. The changes in regulation are often accompanied by changes in the gene expression trajectories, which can be divided into two scenarios: (1) changes of gene expression trajectory shape that reflect changes in cell type composition or altered molecular machinery; (2) temporal shift of gene expression trajectories that indicate different regulation of developmental timing. Therefore, we developed an R package TempShift to separates these two scenarios and demonstrated that TempShift can distinguish temporal shift from different shape (DiffShape) of expression trajectories, and can accurately estimate the time difference between multiple trajectories. We applied TempShift to identify sequential gene expression across 11 neocortical areas, which suggested sequential occurrence of synapse formation and axon guidance, as well as reconstructed interneuron migration pathways within neocortex. Comparison between healthy and DS brains revealed increased microglia, shortened neuronal migration process, and delayed synaptogenesis and myelination in DS. These applications also demonstrate the potential of TempShift in understanding gene expression temporal dynamics during different biological processes.

Scalable Gromov-Wasserstein Based Comparison of Biological Time Series

A time series is an extremely abundant data type arising in many areas of scientific research, including the biological sciences. Any method that compares time series data relies on a pairwise distance between trajectories, and the choice of distance measure determines the accuracy and speed of the time series comparison. This paper introduces an optimal transport type distance for comparing time series trajectories that are allowed to lie in spaces of different dimensions and/or with differing numbers of points possibly unequally spaced along each trajectory. The construction is based on a modified Gromov-Wasserstein distance optimization program, reducing the problem to a Wasserstein distance on the real line. The resulting program has a closed-form solution and can be computed quickly due to the scalability of the one-dimensional Wasserstein distance. We discuss theoretical properties of this distance measure, and empirically demonstrate the performance of the proposed distance on several datasets with a range of characteristics commonly found in biologically relevant data. We also use our proposed distance to demonstrate that averaging oscillatory time series trajectories using the recently proposed Fused Gromov-Wasserstein barycenter retains more characteristics in the averaged trajectory when compared to traditional averaging, which demonstrates the applicability of Fused Gromov-Wasserstein barycenters for biological time series. Fast and user friendly software for computing the proposed distance and related applications is provided. The proposed distance allows fast and meaningful comparison of biological time series and can be efficiently used in a wide range of applications.

Systematic review of comparative transcriptomic studies of cellular resistance to genotoxic stress

The development of resistance by tumor cells to various types of therapy is a significant problem that decreases the effectiveness of oncology treatments. For more than two decades, comparative transcriptomic studies of tumor cells with different sensitivities to ionizing radiation and chemotherapeutic agents have been conducted in order to identify the causes and mechanisms underlying this phenomenon. However, the results of such studies have little in common and often contradict each other. We have assumed that a systematic analysis of a large number of such studies will provide new knowledge about the mechanisms of development of therapeutic resistance in tumor cells. Our comparison of 123 differentially expressed gene (DEG) lists published in 98 papers suggests a very low degree of consistency between the study results. Grouping the data by type of genotoxic agent and tumor type did not increase the similarity. The most frequently overexpressed genes were found to be those encoding the transport protein ABCB1 and the antiviral defense protein IFITM1. We put forward a hypothesis that the role played by the overexpression of the latter in the development of resistance may be associated not only with the stimulation of proliferation, but also with the limitation of exosomal communication and, as a result, with a decrease in the bystander effect. Among down regulated DEGs, BNIP3 was observed most frequently. The expression of BNIP3, together with BNIP3L, is often suppressed in cells resistant to non-platinum genotoxic chemotherapeutic agents, whereas it is increased in cells resistant to ionizing radiation. These observations are likely to be mediated by the binary effects of these gene products on survival, and regulation of apoptosis and autophagy. The combined data also show that even such obvious mechanisms as inhibition of apoptosis and increase of proliferation are not universal but show multidirectional changes.

Fibroblasts are not just fibroblasts: clear differences between dermal and pulmonary fibroblasts' response to fibrotic growth factors

Systemic Sclerosis (SSc) hallmark is skin fibrosis, but up to 80% of the patients have fibrotic involvement in the pulmonary system. Antifibrotic drugs which have failed in a general SSc population have now been approved in patients with SSc-associated interstitial lung disease (ILD). This indicates that the fibrotic progression and regulation of fibroblasts likely depend on local factors specific to the tissue type. This study investigated the difference between dermal and pulmonary fibroblasts in a fibrotic setting, mimicking the extracellular matrix. Primary healthy fibroblasts were grown in a crowded environment and stimulated with TGF-β1 and PDGF-AB. The viability, morphology, migration capacity, extracellular matrix formation, and gene expression were assessed: TGF-β1 only increased the viability in the dermal fibroblasts. PDGF-AB increased the migration capacity of dermal fibroblasts while the pulmonary fibroblasts fully migrated. The morphology of the fibroblasts was different without stimulation. TGF-β1 increased the formation of type III collagen in pulmonary fibroblasts, while PDGF-AB increased it in dermal fibroblasts. The gene expression trend of type VI collagen was the opposite after PDGF-AB stimulation. The fibroblasts exhibit different response profiles to TGF-β1 and PDGF-AB; this suggests that drivers of fibrosis are tissue-dependent, which needs to be considered in drug development.

Inferring evolutionary trajectories from cross-sectional transcriptomic data to mirror lung adenocarcinoma progression

Lung adenocarcinoma (LUAD) is a deadly tumor with dynamic evolutionary process. Although much endeavors have been made in identifying the temporal patterns of cancer progression, it remains challenging to infer and interpret the molecular alterations associated with cancer development and progression. To this end, we developed a computational approach to infer the progression trajectory based on cross-sectional transcriptomic data. Analysis of the LUAD data using our approach revealed a linear trajectory with three different branches for malignant progression, and the results showed consistency in three independent cohorts. We used the progression model to elucidate the potential molecular events in LUAD progression. Further analysis showed that overexpression of BUB1B, BUB1 and BUB3 promoted tumor cell proliferation and metastases by disturbing the spindle assembly checkpoint (SAC) in the mitosis. Aberrant mitotic spindle checkpoint signaling appeared to be one of the key factors promoting LUAD progression. We found the inferred cancer trajectory allows to identify LUAD susceptibility genetic variations using genome-wide association analysis. This result shows the opportunity for combining analysis of candidate genetic factors with disease progression. Furthermore, the trajectory showed clear evident mutation accumulation and clonal expansion along with the LUAD progression. Understanding how tumors evolve and identifying mutated genes will help guide cancer management. We investigated the clonal architectures and identified distinct clones and subclones in different LUAD branches. Validation of the model in multiple independent data sets and correlation analysis with clinical results demonstrate that our method is effective and unbiased.

Predicting gene level sensitivity to JAK-STAT signaling perturbation using a mechanistic-to-machine learning framework

The JAK-STAT pathway integrates complex cytokine signals via a limited number of molecular components, inspiring numerous efforts to clarify the diversity and specificity of STAT transcription factor function. We developed a computational workflow to make global cytokine-induced gene predictions from STAT phosphorylation dynamics, modeling macrophage responses to IL-6 and IL-10, which signal through common STATs, but with distinct temporal dynamics and contrasting functions. Our mechanistic-to-machine learning model identified select cytokine-induced gene sets associated with late pSTAT3 timeframes and a preferential pSTAT1 reduction upon JAK2 inhibition. We predicted and validated the impact of JAK2 inhibition on gene expression, identifying dynamically regulated genes that were sensitive or insensitive to JAK2 variation. Thus, we successfully linked STAT signaling dynamics to gene expression to support future efforts targeting pathology-associated STAT-driven gene sets. This serves as a first step in developing multi-level prediction models to understand and perturb gene expression outputs from signaling systems.

Depicting pseudotime-lagged causality across single-cell trajectories for accurate gene-regulatory inference

Identifying the causal interactions in gene-regulatory networks requires an accurate understanding of the time-lagged relationships between transcription factors and their target genes. Here we describe DELAY (short for Depicting Lagged Causality), a convolutional neural network for the inference of gene-regulatory relationships across pseudotime-ordered single-cell trajectories. We show that combining supervised deep learning with joint probability matrices of pseudotime-lagged trajectories allows the network to overcome important limitations of ordinary Granger causality-based methods, for example, the inability to infer cyclic relationships such as feedback loops. Our network outperforms several common methods for inferring gene regulation and, when given partial ground-truth labels, predicts novel regulatory networks from single-cell RNA sequencing (scRNA-seq) and single-cell ATAC sequencing (scATAC-seq) data sets. To validate this approach, we used DELAY to identify important genes and modules in the regulatory network of auditory hair cells, as well as likely DNA-binding partners for two hair cell cofactors (Hist1h1c and Ccnd1) and a novel binding sequence for the hair cell-specific transcription factor Fiz1. We provide an easy-to-use implementation of DELAY under an open-source license at https://github.com/calebclayreagor/DELAY.

AngClust: Angle Feature-Based Clustering for Short Time Series Gene Expression Profiles

When clustering gene expression, it is expected that correlation coefficients of genes in the same clusters are high, and that gene ontology (GO) enrichment analysis of most clusters will be significant. However, existing short-term gene expression clustering algorithms have limitations. To address this problem, we proposed a novel clustering process based on angular features for short-term gene expression. Our method (named AngClust) uses angular features to indicate the change of trend in gene expression levels at two neighboring time points. The changes of angles at multiple time points reflects the change of trend of the overall expression levels. Such changes are used to measure whether the expression trends of different genes are similar. To obtain functionally significant clusters from the clustering results, we evaluated numbers of genes in clusters, average correlation coefficient, fluctuation, and their correlation with GO term enrichment. The efficacy of AngClust outperform two other measures, Euclidean distance (ED) and dynamic time warping of correlation (DTW), on a dataset of yeast gene expression. The ratios of GO and pathway term-enriched of clusters of AngClust is higher than or equal to that of STEM and TMixClust on human, mouse, and yeast time series of gene expression.

Features extracted using tensor decomposition reflect the biological features of the temporal patterns of human blood multimodal metabolome

High-throughput omics technologies have enabled the profiling of entire biological systems. For the biological interpretation of such omics data, two analyses, hypothesis- and data-driven analyses including tensor decomposition, have been used. Both analyses have their own advantages and disadvantages and are mutually complementary; however, a direct comparison of these two analyses for omics data is poorly examined.We applied tensor decomposition (TD) to a dataset representing changes in the concentrations of 562 blood molecules at 14 time points in 20 healthy human subjects after ingestion of 75 g oral glucose. We characterized each molecule by individual dependence (constant or variable) and time dependence (later peak or early peak). Three of the four features extracted by TD were characterized by our previous hypothesis-driven study, indicating that TD can extract some of the same features obtained by hypothesis-driven analysis in a non-biased manner. In contrast to the years taken for our previous hypothesis-driven analysis, the data-driven analysis in this study took days, indicating that TD can extract biological features in a non-biased manner without the time-consuming process of hypothesis generation.

Migraine genetics: Status and road forward

BACKGROUND

Migraine is considered a multifactorial genetic disorder. Different platforms and methods are used to unravel the genetic basis of migraine. Initially, linkage analysis in multigenerational families followed by Sanger sequencing of protein-coding parts (exons) of genes in the genomic region shared by affected family members identified high-effect risk DNA mutations for rare Mendelian forms of migraine, foremost hemiplegic migraine. More recently, genome-wide association studies testing millions of DNA variants in large groups of patients and controls have proven successful in identifying many dozens of low-effect risk DNA variants for the more common forms of migraine with the number of associated DNA variants increasing steadily with larger sample sizes. Currently, next-generation sequencing, utilising whole exome and whole genome sequence data, and other omics data are being used to facilitate their functional interpretation and the discovery of additional risk factors. Various methods and analysis tools, such as genetic correlation and causality analysis, are used to further characterise genetic risk factors.

FINDINGS

We describe recent findings in genome-wide association studies and next-generation sequencing analysis in migraine. We show that the combined results of the two most recent and most powerful migraine genome-wide association studies have identified a total of 178 LD-independent (r² < 0.1) genome-wide significant single nucleotide polymorphisms (SNPs), of which 99 were unique to Hautakangas et al., 11 were unique to Choquet et al., and 68 were identified by both studies. When considering that Choquet et al. also identified three SNPs in a female-specific genome-wide association studies then these two recent studies identified 181 independent SNPs robustly associated with migraine. Cross-trait and causal analyses are beginning to identify and characterise specific biological factors that contribute to migraine risk and its comorbid conditions.

CONCLUSION

This review provides a timely update and overview of recent genetic findings in migraine.

Machine Learning Methods for Cancer Classification Using Gene Expression Data: A Review

Cancer is a term that denotes a group of diseases caused by the abnormal growth of cells that can spread in different parts of the body. According to the World Health Organization (WHO), cancer is the second major cause of death after cardiovascular diseases. Gene expression can play a fundamental role in the early detection of cancer, as it is indicative of the biochemical processes in tissue and cells, as well as the genetic characteristics of an organism. Deoxyribonucleic acid (DNA) microarrays and ribonucleic acid (RNA)-sequencing methods for gene expression data allow quantifying the expression levels of genes and produce valuable data for computational analysis. This study reviews recent progress in gene expression analysis for cancer classification using machine learning methods. Both conventional and deep learning-based approaches are reviewed, with an emphasis on the application of deep learning models due to their comparative advantages for identifying gene patterns that are distinctive for various types of cancers. Relevant works that employ the most commonly used deep neural network architectures are covered, including multi-layer perceptrons, as well as convolutional, recurrent, graph, and transformer networks. This survey also presents an overview of the data collection methods for gene expression analysis and lists important datasets that are commonly used for supervised machine learning for this task. Furthermore, we review pertinent techniques for feature engineering and data preprocessing that are typically used to handle the high dimensionality of gene expression data, caused by a large number of genes present in data samples. The paper concludes with a discussion of future research directions for machine learning-based gene expression analysis for cancer classification.

Time-resolved RNA signatures of CD4+ T cells in Parkinson's disease

Parkinson's disease (PD) emerges as a complex, multifactorial disease. While there is increasing evidence that dysregulated T cells play a central role in PD pathogenesis, elucidation of the pathomechanical changes in related signaling is still in its beginnings. We employed time-resolved RNA expression upon the activation of peripheral CD4+ T cells to track and functionally relate changes on cellular signaling in representative cases of patients at different stages of PD. While only few miRNAs showed time-course related expression changes in PD, we identified groups of genes with significantly altered expression for each different time window. Towards a further understanding of the functional consequences, we highlighted pathways with decreased or increased activity in PD, including the most prominent altered IL-17 pathway. Flow cytometric analyses showed not only an increased prevalence of Th17 cells but also a specific subtype of IL-17 producing γδ-T cells, indicating a previously unknown role in PD pathogenesis.

Inferring gene regulatory networks from single-cell gene expression data via deep multi-view contrastive learning

The inference of gene regulatory networks (GRNs) is of great importance for understanding the complex regulatory mechanisms within cells. The emergence of single-cell RNA-sequencing (scRNA-seq) technologies enables the measure of gene expression levels for individual cells, which promotes the reconstruction of GRNs at single-cell resolution. However, existing network inference methods are mainly designed for data collected from a single data source, which ignores the information provided by multiple related data sources. In this paper, we propose a multi-view contrastive learning (DeepMCL) model to infer GRNs from scRNA-seq data collected from multiple data sources or time points. We first represent each gene pair as a set of histogram images, and then introduce a deep Siamese convolutional neural network with contrastive loss to learn the low-dimensional embedding for each gene pair. Moreover, an attention mechanism is introduced to integrate the embeddings extracted from different data sources and different neighbor gene pairs. Experimental results on synthetic and real-world datasets validate the effectiveness of our contrastive learning and attention mechanisms, demonstrating the effectiveness of our model in integrating multiple data sources for GRN inference.

Comparative network analysis reveals the dynamics of organic acid diversity during fruit ripening in peach (Prunus persica L. Batsch)

BACKGROUND

Organic acids are important components that determine the fruit flavor of peach (Prunus persica L. Batsch). However, the dynamics of organic acid diversity during fruit ripening and the key genes that modulate the organic acids metabolism remain largely unknown in this kind of fruit tree which yield ranks sixth in the world.

RESULTS

In this study, we used 3D transcriptome data containing three dimensions of information, namely time, phenotype and gene expression, from 5 different varieties of peach to construct gene co-expression networks throughout fruit ripening of peach. With the network inferred, the time-ordered network comparative analysis was performed to select high-acid specific gene co-expression network and then clarify the regulatory factors controlling organic acid accumulation. As a result, network modules related to organic acid synthesis and metabolism under high-acid and low-acid comparison conditions were identified for our following research. In addition, we obtained 20 candidate genes as regulatory factors related to organic acid metabolism in peach.

CONCLUSIONS

The study provides new insights into the dynamics of organic acid accumulation during fruit ripening, complements the results of classical co-expression network analysis and establishes a foundation for key genes discovery from time-series multiple species transcriptome data.

FTO/RUNX2 signaling axis promotes cementoblast differentiation under normal and inflammatory condition

N⁶-methyladenosine (m⁶A) is the most prevalent mRNA modification which plays crucial roles in various biological processes, but its role in cementogenesis remains largely unknown. Here, using time-series transcriptomic analysis, we reveal that mRNA m⁶A demethylase Fat mass and obesity-associated protein (FTO) is involved in cementogenesis. Knocking down FTO decreases cementoblast differentiation and mineralization in both OCCM-30 cellular model and murine ectopic bone formation model. Mechanistically, we find that FTO directly binds Runt-related transcription factor 2 (Runx2) mRNA, an important cementogenesis factor, thus protecting it from YTH domain-containing family protein 2 (YTHDF2) mediated degradation, when cementoblasts are differentiating. Knocking down YTHDF2 restores the expression of Runx2 in FTO-knockdown cells. Moreover, under inflammatory conditions, TNF-α inhibits cementoblast differentiation and mineralization partly through FTO/RUNX2 axis. Collectively, our study reveals an important regulatory role of FTO/RUNX2 axis in normal and pathological cementogenesis.

Stakeholder relationships and corporate social goal orientation: Implications for entrepreneurial psychology

As the sensitivity to corporate social responsibility (CSR) continues to grow, the goal of enterprises has expanded beyond the sole pursuit of economic value. Corporate social goal orientation has therefore come to occupy a central position in entrepreneurs' psychology and the transition away from a market-only economy. This study uses secondary data from 4,288 samples of 725 Chinese-listed companies from 2009 to 2020 to explore the driving factors in social goal orientation based on the characteristics of sample companies and their industry groups from the perspective of stakeholder relationships. The results can be summarized as follows: (1) there is an inverted U-shaped relationship between government stakeholder relationships and social goal orientation, and there is a significant positive relationship between financial stakeholder relationships, market stakeholder relationships, and corporate social goal orientation. (2) The correlation between single-dual stakeholder relationships and social goal orientation is not consistent. In light of the nature of the roles of government and the market, the correlation between the government-market dual relationship and corporate social goal orientation is not significant. However, there is a significant correlation between the finance-government dual stakeholder relationship and social goal orientation; that is, the dual stakeholder relationship maintains the existence of non-institutional capital and corporate financial capital. Moreover, there is no significant correlation between the market-finance dual relationship and corporate social goal orientation, and there is substitutability between market and financial stakeholder relationships. With the deepening of our understanding of CSR, the core goal of enterprises is no longer confined to the pursuit of economic value, and their social goal orientation has come to be regarded as a major driving force in sustainable development. This study enriches the research on the relationship between stakeholder relationships and shows that stakeholder relationships also have important significance to both achieving corporate goals and shaping entrepreneurs' psychology.

Efficient reconstruction of directed networks from noisy dynamics using stochastic force inference

We consider coupled network dynamics under uncorrelated noises that fluctuate about the noise-free long-time asymptotic state. Our goal is to reconstruct the directed network only from the time-series data of the dynamics of the nodes. By using the stochastic force inference method with a simple natural choice of linear polynomial basis, we derive a reconstruction scheme of the connection weights and the noise strength of each node. Explicit simulations for directed and undirected random networks with various node dynamics are carried out to demonstrate the good accuracy and high efficiency of the reconstruction scheme. We further consider the case when only a subset of the network and its node dynamics can be observed, and it is demonstrated that the directed weighted connections among the observed nodes can be easily and faithfully reconstructed. In addition, we propose a scheme to infer the number of hidden nodes and their effects on each observed node. The accuracy of these results is illustrated by simulations.

Inferring gene regulation from stochastic transcriptional variation across single cells at steady state

Regulatory relationships between transcription factors (TFs) and their target genes lie at the heart of cellular identity and function; however, uncovering these relationships is often labor-intensive and requires perturbations. Here, we propose a principled framework to systematically infer gene regulation for all TFs simultaneously in cells at steady state by leveraging the intrinsic variation in the transcriptional abundance across single cells. Through modeling and simulations, we characterize how transcriptional bursts of a TF gene are propagated to its target genes, including the expected ranges of time delay and magnitude of maximum covariation. We distinguish these temporal trends from the time-invariant covariation arising from cell states, and we delineate the experimental and technical requirements for leveraging these small but meaningful cofluctuations in the presence of measurement noise. While current technology does not yet allow adequate power for definitively detecting regulatory relationships for all TFs simultaneously in cells at steady state, we investigate a small-scale dataset to inform future experimental design. This study supports the potential value of mapping regulatory connections through stochastic variation, and it motivates further technological development to achieve its full potential.

Constrained Fourier estimation of short-term time-series gene expression data reduces noise and improves clustering and gene regulatory network predictions

BACKGROUND

Biological data suffers from noise that is inherent in the measurements. This is particularly true for time-series gene expression measurements. Nevertheless, in order to to explore cellular dynamics, scientists employ such noisy measurements in predictive and clustering tools. However, noisy data can not only obscure the genes temporal patterns, but applying predictive and clustering tools on noisy data may yield inconsistent, and potentially incorrect, results.

RESULTS

To reduce the noise of short-term (< 48 h) time-series expression data, we relied on the three basic temporal patterns of gene expression: waves, impulses and sustained responses. We constrained the estimation of the true signals to these patterns by estimating the parameters of first and second-order Fourier functions and using the nonlinear least-squares trust-region optimization technique. Our approach lowered the noise in at least 85% of synthetic time-series expression data, significantly more than the spline method ([Formula: see text]). When the data contained a higher signal-to-noise ratio, our method allowed downstream network component analyses to calculate consistent and accurate predictions, particularly when the noise variance was high. Conversely, these tools led to erroneous results from untreated noisy data. Our results suggest that at least 5-7 time points are required to efficiently de-noise logarithmic scaled time-series expression data. Investing in sampling additional time points provides little benefit to clustering and prediction accuracy.

CONCLUSIONS

Our constrained Fourier de-noising method helps to cluster noisy gene expression and interpret dynamic gene networks more accurately. The benefit of noise reduction is large and can constitute the difference between a successful application and a failing one.

Constitutive gene expression differs in three brain regions important for cognition in neophobic and non-neophobic house sparrows (Passer domesticus)

Neophobia (aversion to new objects, food, and environments) is a personality trait that affects the ability of wildlife to adapt to new challenges and opportunities. Despite the ubiquity and importance of this trait, the molecular mechanisms underlying repeatable individual differences in neophobia in wild animals are poorly understood. We evaluated wild-caught house sparrows (Passer domesticus) for neophobia in the lab using novel object tests. We then selected a subset of neophobic and non-neophobic individuals (n = 3 of each, all females) and extracted RNA from four brain regions involved in learning, memory, threat perception, and executive function: striatum, caudal dorsomedial hippocampus, medial ventral arcopallium, and caudolateral nidopallium (NCL). Our analysis of differentially expressed genes (DEGs) used 11,889 gene regions annotated in the house sparrow reference genome for which we had an average of 25.7 million mapped reads/sample. PERMANOVA identified significant effects of brain region, phenotype (neophobic vs. non-neophobic), and a brain region by phenotype interaction. Comparing neophobic and non-neophobic birds revealed constitutive differences in DEGs in three of the four brain regions examined: hippocampus (12% of the transcriptome significantly differentially expressed), striatum (4%) and NCL (3%). DEGs included important known neuroendocrine mediators of learning, memory, executive function, and anxiety behavior, including serotonin receptor 5A, dopamine receptors 1, 2 and 5 (downregulated in neophobic birds), and estrogen receptor beta (upregulated in neophobic birds). These results suggest that some of the behavioral differences between phenotypes may be due to underlying gene expression differences in the brain. The large number of DEGs in neophobic and non-neophobic birds also implies that there are major differences in neural function between the two phenotypes that could affect a wide variety of behavioral traits beyond neophobia.

Fast and accurate inference of gene regulatory networks through robust precision matrix estimation

MOTIVATION

Transcriptional regulation mechanisms allow cells to adapt and respond to external stimuli by altering gene expression. The possible cell transcriptional states are determined by the underlying gene regulatory network (GRN), and reliably inferring such network would be invaluable to understand biological processes and disease progression.

RESULTS

In this article, we present a novel method for the inference of GRNs, called PORTIA, which is based on robust precision matrix estimation, and we show that it positively compares with state-of-the-art methods while being orders of magnitude faster. We extensively validated PORTIA using the DREAM and MERLIN+P datasets as benchmarks. In addition, we propose a novel scoring metric that builds on graph-theoretical concepts.

AVAILABILITY AND IMPLEMENTATION

The code and instructions for data acquisition and full reproduction of our results are available at https://github.com/AntoinePassemiers/PORTIA-Manuscript. PORTIA is available on PyPI as a Python package (portia-grn).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Reconstruction of Gene Regulatory Networks Using Multiple Datasets

MOTIVATION

Laboratory gene regulatory data for a species are sporadic. Despite the abundance of gene regulatory network algorithms that employ single data sets, few algorithms can combine the vast but disperse sources of data and extract the potential information. With a motivation to compensate for this shortage, we developed an algorithm called GENEREF that can accumulate information from multiple types of data sets in an iterative manner, with each iteration boosting the performance of the prediction results.

RESULTS

The algorithm is examined extensively on data extracted from the quintuple DREAM4 networks and DREAM5's Escherichia coli and Saccharomyces cerevisiae networks and sub-networks. Many single-dataset and multi-dataset algorithms were compared to test the performance of the algorithm. Results show that GENEREF surpasses non-ensemble state-of-the-art multi-perturbation algorithms on the selected networks and is competitive to present multiple-dataset algorithms. Specifically, it outperforms dynGENIE3 and is on par with iRafNet. Also, we argued that a scoring method solely based on the AUPR criterion would be more trustworthy than the traditional score.

AVAILABILITY

The Python implementation along with the data sets and results can be downloaded from github.com/msaremi/GENEREF.

Profiling and Functional Analysis of mRNAs during Skeletal Muscle Differentiation in Goats

Skeletal myogenesis is a complicated biological event that involves a succession of tightly controlled gene expressions. In order to identify novel regulators of this process, we performed mRNA-Seq studies of goat skeletal muscle satellite cells (MuSCs) cultured under proliferation (GM) and differentiation (DM1/DM5) conditions. A total of 19,871 goat genes were expressed during these stages, 198 of which represented novel transcripts. Notably, in pairwise comparisons at the different stages, 2551 differentially expressed genes (DEGs) were identified (p < 0.05), including 1560 in GM vs. DM1, 1597 in GM vs. DM5, and 959 in DM1 vs. DM5 DEGs. The time-series expression profile analysis clustered the DEGs into eight gene groups, three of which had significantly upregulated and downregulated patterns (p < 0.05). Functional enrichment analysis showed that DEGs were enriched for essential biological processes such as muscle structure development, muscle contraction, muscle cell development, striated muscle cell differentiation, and myofibril assembly, and were involved in pathways such as the MAPK, Wnt and PPAR signaling pathways. Moreover, the expression of eight DEGs (MYL2, DES, MYOG, FAP, PLK2, ADAM, WWC1, and PRDX1) was validated. These findings offer novel insights into the transcriptional regulation of skeletal myogenesis in goats.

Characterization and identification of cell death dynamics by quantitative phase imaging

SIGNIFICANCE

Investigating cell death dynamics at the single-cell level plays an essential role in biological research. Quantitative phase imaging (QPI), a label-free method without adverse effects of exogenous labels, has been widely used to image many types of cells under various conditions. However, the dynamics of QPI features during cell death have not been thoroughly characterized.

AIM

We aim to develop a label-free technique to quantitatively characterize single-cell dynamics of cellular morphology and intracellular mass distribution of cells undergoing apoptosis and necrosis.

APPROACH

QPI was used to capture time-lapse phase images of apoptotic, necrotic, and normal cells. The dynamics of morphological and QPI features during cell death were fitted by a sigmoid function to quantify both the extent and rate of changes.

RESULTS

The two types of cell death mainly differed from normal cells in the lower phase of the central region and differed from each other in the sharp nuclear boundary shown in apoptotic cells.

CONCLUSIONS

The proposed method characterizes the dynamics of cellular morphology and intracellular mass distributions, which could be applied to studying cells undergoing state transition such as drug response.

Detecting Ultra- and Circadian Activity Rhythms of Dairy Cows in Automatic Milking Systems Using the Degree of Functional Coupling—A Pilot Study

Ultra- and circadian activity rhythms of animals can provide important insights into animal welfare. The consistency of behavioral patterns is characteristic of healthy organisms, while changes in the regularity of behavioral rhythms may indicate health and stress-related challenges. This pilot study aimed to examine whether dairy cows in free-stall barns with an automatic milking system (AMS) and free cow traffic can develop ultra- and circadian activity rhythms. On 4 dairy farms, pedometers recorded the activity of 10 cows each over 28 days. Based on time series calculation, the Degree of Functional Coupling (DFC) was used to determine the cows' activity rhythms. The DFC identified significant rhythmic patterns in sliding 7-day periods and indicated the percentage of activity (0–100%) that was synchronized with the 24-h day-night rhythm. As light is the main factor influencing the sleep-wake cycle of organisms, light intensity was recorded in the AMS, at the feed alley and in the barn of each farm. In addition, feeding and milking management were considered as part of the environmental context. Saliva samples of each cow were taken every 3 h for 1 day to determine the melatonin concentration. The DFC approach was successfully used to detect activity rhythms of dairy cows in commercial housing systems. However, large inter- and intra-individual variations were observed. Due to a high frequency of 0 and 100%, a median split was used to dichotomize into “low” (&lt;72.34%) and “high” (≥72.34%) DFC. Forty percent of the sliding 7-day periods corresponded to a low DFC and 50% to a high DFC. No DFC could be calculated for 10% of the periods, as the cows' activity was not synchronized to 24 h. A generalized linear mixed-effects model revealed that the DFC levels were positively associated with a longer milking interval and a higher amount of daytime activity and negatively associated with higher number of lactations. The DFC is a novel approach to animal behavior monitoring. Due to its automation capability, it represents a promising tool in its further development for the purpose of longitudinal monitoring of animal welfare.

Guidelines for bioinformatics of single-cell sequencing data analysis in Alzheimer's disease: review, recommendation, implementation and application

Alzheimer's disease (AD) is the most common form of dementia, characterized by progressive cognitive impairment and neurodegeneration. Extensive clinical and genomic studies have revealed biomarkers, risk factors, pathways, and targets of AD in the past decade. However, the exact molecular basis of AD development and progression remains elusive. The emerging single-cell sequencing technology can potentially provide cell-level insights into the disease. Here we systematically review the state-of-the-art bioinformatics approaches to analyze single-cell sequencing data and their applications to AD in 14 major directions, including 1) quality control and normalization, 2) dimension reduction and feature extraction, 3) cell clustering analysis, 4) cell type inference and annotation, 5) differential expression, 6) trajectory inference, 7) copy number variation analysis, 8) integration of single-cell multi-omics, 9) epigenomic analysis, 10) gene network inference, 11) prioritization of cell subpopulations, 12) integrative analysis of human and mouse sc-RNA-seq data, 13) spatial transcriptomics, and 14) comparison of single cell AD mouse model studies and single cell human AD studies. We also address challenges in using human postmortem and mouse tissues and outline future developments in single cell sequencing data analysis. Importantly, we have implemented our recommended workflow for each major analytic direction and applied them to a large single nucleus RNA-sequencing (snRNA-seq) dataset in AD. Key analytic results are reported while the scripts and the data are shared with the research community through  GitHub. In summary, this comprehensive review provides insights into various approaches to analyze single cell sequencing data and offers specific guidelines for study design and a variety of analytic directions. The review and the accompanied software tools will serve as a valuable resource for studying cellular and molecular mechanisms of AD, other diseases, or biological systems at the single cell level.