Reconstruction Of Gene Regulatory Networks Research Articles

Scalable identification of lineage-specific gene regulatory networks from metacells with NetID

The identification of gene regulatory networks (GRNs) is crucial for understanding cellular differentiation. Single-cell RNA sequencing data encode gene-level covariations at high resolution, yet data sparsity and high dimensionality hamper accurate and scalable GRN reconstruction. To overcome these challenges, we introduce NetID leveraging homogenous metacells while avoiding spurious gene–gene correlations. Benchmarking demonstrates superior performance of NetID compared to imputation-based methods. By incorporating cell fate probability information, NetID facilitates the prediction of lineage-specific GRNs and recovers known network motifs governing bone marrow hematopoiesis, making it a powerful toolkit for deciphering gene regulatory control of cellular differentiation from large-scale single-cell transcriptome data.

GRNMOPT: Inference of gene regulatory networks based on a multi-objective optimization approach

Background and objectiveThe reconstruction of gene regulatory networks (GRNs) stands as a vital approach in deciphering complex biological processes. The application of nonlinear ordinary differential equations (ODEs) models has demonstrated considerable efficacy in predicting GRNs. Notably, the decay rate and time delay are pivotal in authentic gene regulation, yet their systematic determination in ODEs models remains underexplored. The development of a comprehensive optimization framework for the effective estimation of these key parameters is essential for accurate GRN inference. MethodThis study introduces GRNMOPT, an innovative methodology for inferring GRNs from time-series and steady-state data. GRNMOPT employs a combined use of decay rate and time delay in constructing ODEs models to authentically represent gene regulatory processes. It incorporates a multi-objective optimization approach, optimizing decay rate and time delay concurrently to derive Pareto optimal sets for these factors, thereby maximizing accuracy metrics such as AUROC (Area Under the Receiver Operating Characteristic curve) and AUPR (Area Under the Precision-Recall curve). Additionally, the use of XGBoost for calculating feature importance aids in identifying potential regulatory gene links. ResultsComprehensive experimental evaluations on two simulated datasets from DREAM4 and three real gene expression datasets (Yeast, In vivo Reverse-engineering and Modeling Assessment [IRMA], and Escherichia coli [E. coli]) reveal that GRNMOPT performs commendably across varying network scales. Furthermore, cross-validation experiments substantiate the robustness of GRNMOPT. ConclusionWe propose a novel approach called GRNMOPT to infer GRNs based on a multi-objective optimization framework, which effectively improves inference accuracy and provides a powerful tool for GRNs inference.

A multi-layered integrative analysis reveals a cholesterol metabolic program in outer radial glia with implications for human brain evolution.

The definition of molecular and cellular mechanisms contributing to brain ontogenetic trajectories is essential to investigate the evolution of our species. Yet their functional dissection at an appropriate level of granularity remains challenging. Capitalizing on recent efforts that have extensively profiled neural stem cells from the developing human cortex, we develop an integrative computational framework to perform trajectory inference and gene regulatory network reconstruction, (pseudo)time-informed non-negative matrix factorization for learning the dynamics of gene expression programs, and paleogenomic analysis for a higher-resolution mapping of derived regulatory variants in our species in comparison with our closest relatives. We provide evidence for cell type-specific regulation of gene expression programs during indirect neurogenesis. In particular, our analysis uncovers a key role for a cholesterol program in outer radial glia, regulated by zinc-finger transcription factor KLF6. A cartography of the regulatory landscape impacted by Homo sapiens-derived variants reveals signals of selection clustering around regulatory regions associated with GLI3, a well-known regulator of radial glial cell cycle, and impacting KLF6 regulation. Our study contributes to the evidence of significant changes in metabolic pathways in recent human brain evolution.

GRAMP: A gene ranking and model prioritisation framework for building consensus genetic networks

Despite significant recent advancements in computational and statistical methods, different methods have specific strengths and weaknesses in the accurate reconstruction of gene regulatory networks (GRNs), making it difficult to determine the best method for each specific problem. To overcome these challenges, ensemble approaches, which combine the strengths of individual inference methods, are valuable. However, existing ensemble methods for GRN inference lack a sophisticated network aggregation method and generally rely solely on ranking approaches. These ensemble methods have no reliable mechanisms to identify highly performing inference methods specific to a given problem. They therefore tend to aggregate weak methods, diminishing the overall accuracy of the approach. Thus, developing a reliable mechanism to identify the most effective methods for specific problems and prioritize them in consensus network building is important. This paper presents a novel ensemble approach for reconstructing GRNs by integrating previously developed diverse GRN inference approaches. A novel network aggregation method called GRAMP, Gene Ranking And Model Prioritisation framework was developed, taking into consideration both local and global gene ranking and the performance of different inference approaches on a specific network. The proposed ensemble approach demonstrated performance superior to those of other state-of-the-art methods, as evidenced by results from simulated datasets and a real-world gene expression dataset.

Reconstruction of gene regulatory networks for Caenorhabditis elegans using tree-shaped gene expression data.

Constructing gene regulatory networks is a widely adopted approach for investigating gene regulation, offering diverse applications in biology and medicine. A great deal of research focuses on using time series data or single-cell RNA-sequencing data to infer gene regulatory networks. However, such gene expression data lack either cellular or temporal information. Fortunately, the advent of time-lapse confocal laser microscopy enables biologists to obtain tree-shaped gene expression data of Caenorhabditis elegans, achieving both cellular and temporal resolution. Although such tree-shaped data provide abundant knowledge, they pose challenges like non-pairwise time series, laying the inaccuracy of downstream analysis. To address this issue, a comprehensive framework for data integration and a novel Bayesian approach based on Boolean network with time delay are proposed. The pre-screening process and Markov Chain Monte Carlo algorithm are applied to obtain the parameter estimates. Simulation studies show that our method outperforms existing Boolean network inference algorithms. Leveraging the proposed approach, gene regulatory networks for five subtrees are reconstructed based on the real tree-shaped datatsets of Caenorhabditis elegans, where some gene regulatory relationships confirmed in previous genetic studies are recovered. Also, heterogeneity of regulatory relationships in different cell lineage subtrees is detected. Furthermore, the exploration of potential gene regulatory relationships that bear importance in human diseases is undertaken. All source code is available at the GitHub repository https://github.com/edawu11/BBTD.git.

Graph attention network with convolutional layer for predicting gene regulations from single-cell ribonucleic acid sequence data

Reconstructing gene regulatory networks (GRNs) is an important task to reveal the regulatory relationship between genes and understand the mechanism of intracellular gene expression regulation. With the development of single-cell ribonucleic acid sequencing (scRNA-seq) technology, researchers begin to attempt to infer GRN within cells. In this paper, we propose graph attention network with convolutional layer (GATCL) to infer the latent interactions between transcription factors (TFs) and target genes in GRN. Firstly, GATCL uses graph attention network (GAT), which can effectively extract information about genes and TFs. Secondly, we combine multi-head attention layer with one-head attention layer and propose a new method of using convolution instead of weight matrix. Thirdly, we use the exponential linear unit (ELU) activation function to replace the leaky rectified linear unit (LReLU) commonly used in GAT, which further improves the accuracy of GATCL. The AUROC of our method on seven scRNA-seq datasets with four types of ground-truth networks reached an average of 0.827, which is higher than other state-of-the-art models. GATCL applies a supervised deep learning algorithm to solve the problems existing in the inference of GRN from scRNA-seq in the engineering field, and the effectiveness of this method is verified by a large number of experiments.

CRISPR-GEM: A Novel Machine Learning Model for CRISPR Genetic Target Discovery and Evaluation.

CRISPR gene editing strategies are shaping cell therapies through precise and tunable control over gene expression. However, achieving reliable therapeutic effects with improved safety and efficacy requires informed target gene selection. This depends on a thorough understanding of the involvement of target genes in gene regulatory networks (GRNs) that regulate cell phenotype and function. Machine learning models have been previously used for GRN reconstruction using RNA-seq data, but current techniques are limited to single cell types and focus mainly on transcription factors. This restriction overlooks many potential CRISPR target genes, such as those encoding extracellular matrix components, growth factors, and signaling molecules, thus limiting the applicability of these models for CRISPR strategies. To address these limitations, we have developed CRISPR-GEM, a multi-layer perceptron (MLP)-based synthetic GRN constructed to accurately predict the downstream effects of CRISPR gene editing. First, input and output nodes are identified as differentially expressed genes between defined experimental and target cell/tissue types respectively. Then, MLP training learns regulatory relationships in a black-box approach allowing accurate prediction of output gene expression using only input gene expression. Finally, CRISPR-mimetic perturbations are made to each input gene individually and the resulting model predictions are compared to those for the target group to score and assess each input gene as a CRISPR candidate. The top scoring genes provided by CRISPR-GEM therefore best modulate experimental group GRNs to motivate transcriptomic shifts towards a target group phenotype. This machine learning model is the first of its kind for predicting optimal CRISPR target genes and serves as a powerful tool for enhanced CRISPR strategies across a range of cell therapies.

SFINN: inferring gene regulatory network from single-cell and spatial transcriptomic data with shared factor neighborhood and integrated neural network.

The rise of single-cell RNA sequencing (scRNA-seq) technology presents new opportunities for constructing detailed cell type-specific gene regulatory networks (GRNs) to study cell heterogeneity. However, challenges caused by noises, technical errors, and dropout phenomena in scRNA-seq data pose significant obstacles to GRN inference, making the design of accurate GRN inference algorithms still essential. The recent growth of both single-cell and spatial transcriptomic sequencing data enables the development of supervised deep learning methods to infer GRNs on these diverse single-cell datasets. In this study, we introduce a novel deep learning framework based on shared factor neighborhood and integrated neural network (SFINN) for inferring potential interactions and causalities between transcription factors and target genes from single-cell and spatial transcriptomic data. SFINN utilizes shared factor neighborhood to construct cellular neighborhood network based on gene expression data and additionally integrates cellular network generated from spatial location information. Subsequently, the cell adjacency matrix and gene pair expression are fed into an integrated neural network framework consisting of a graph convolutional neural network and a fully-connected neural network to determine whether the genes interact. Performance evaluation in the tasks of gene interaction and causality prediction against the existing GRN reconstruction algorithms demonstrates the usability and competitiveness of SFINN across different kinds of data. SFINN can be applied to infer GRNs from conventional single-cell sequencing data and spatial transcriptomic data. SFINN can be accessed at GitHub: https://github.com/JGuan-lab/SFINN.

Reconstruction of gene regulatory networks using graph neural networks

Gene regulatory network (GRN) inference, a longstanding challenge in computational biology, aims to construct GRNs from genomic data. Graph Neural Networks (GNNs) are well-suited for this task due to their ability to leverage both node features and topological relationships. This research systematically evaluated various GNN variants, gradually narrowing the focus through a filtering process. The study considered multiple design aspects, including layers, epochs, decoders, activation functions, graph structures, aggregation methods, skip connections, dropout, and hidden dimensions. Ultimately, two promising models emerged, one based on the Chebyshev spectral graph convolutional operator and the other on the Hypergraph convolutional operator, demonstrating state-of-the-art performance. Notably, hypergraphs demonstrated superior performance on real datasets with higher-order dependencies, while the Chebyshev model showed greater generalization across both simulated and real datasets. The code for this research is available online at https://github.com/EmmaDPaul/GRN-inference-using-GNN.

Network-based prioritization and validation of regulators of vascular smooth muscle cell proliferation in disease

Aberrant vascular smooth muscle cell (VSMC) homeostasis and proliferation characterize vascular diseases causing heart attack and stroke. Here we elucidate molecular determinants governing VSMC proliferation by reconstructing gene regulatory networks from single-cell transcriptomics and epigenetic profiling. We detect widespread activation of enhancers at disease-relevant loci in proliferation-predisposed VSMCs. We compared gene regulatory network rewiring between injury-responsive and nonresponsive VSMCs, which suggested shared transcription factors but differing target loci between VSMC states. Through in silico perturbation analysis, we identified and prioritized previously unrecognized regulators of proliferation, including RUNX1 and TIMP1. Moreover, we showed that the pioneer transcription factor RUNX1 increased VSMC responsiveness and that TIMP1 feeds back to promote VSMC proliferation through CD74-mediated STAT3 signaling. Both RUNX1 and the TIMP1–CD74 axis were expressed in human VSMCs, showing low levels in normal arteries and increased expression in disease, suggesting clinical relevance and potential as vascular disease targets.

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.

Distributed Heterogeneous Transfer Learning

Transfer learning has proved to be effective for building predictive models even in complex conditions with a low amount of available labeled data, by constructing a predictive model for a target domain also using the knowledge coming from a separate domain, called source domain. However, several existing transfer learning methods assume identical feature spaces between the source and the target domains. This assumption limits the possible real-world applications of such methods, since two separate, although related, domains could be described by totally different feature spaces. Heterogeneous transfer learning methods aim to overcome this limitation, but they usually i) make other assumptions on the features, such as requiring the same number of features, ii) are not generally able to distribute the workload over multiple computational nodes, iii) cannot work in the Positive-Unlabeled (PU) learning setting, which we also considered in this study, or iv) their applicability is limited to specific application domains, i.e., they are not general-purpose methods.In this manuscript, we present a novel distributed heterogeneous transfer learning method, implemented in Apache Spark, that overcomes all the above-mentioned limitations. Specifically, it is able to work also in the PU learning setting by resorting to a clustering-based approach, and can align totally heterogeneous feature spaces, without exploiting peculiarities of specific application domains. Moreover, our distributed approach allows us to process large source and target datasets.Our experimental evaluation was performed in three different application domains that can benefit from transfer learning approaches, namely the reconstruction of the human gene regulatory network, the prediction of cerebral stroke in hospital patients, and the prediction of customer energy consumption in power grids. The results show that the proposed approach is able to outperform 4 state-of-the-art heterogeneous transfer learning approaches and 3 baselines, and exhibits ideal performances in terms of scalability.

Abstract PO2-01-08: Differential HOXB13 gene expression and promoter methylation analysis in breast cancer

Abstract Background: The Breast Cancer Index (BCI) is a gene expression-based signature comprising two functional biomarker panels: the Molecular Grade Index (MGI) and BCI (H/I), which is a ratio of the HOXB13 and IL17BR gene expression. Several studies have shown that BCI (H/I) is a predictive biomarker for extended adjuvant endocrine therapy benefit in hormone receptor-positive (HR+) early-stage breast cancer. Molecular mechanisms underlying endocrine responsiveness remain to be fully explored. The objective of this analysis is to evaluate the correlation between gene expression and methylation in breast cancer patients displaying different HOXB13 mRNA expression levels from The Cancer Genome Atlas (TCGA) project. Methods: Data from methylation microarrays and mRNA sequencing were examined in combination with clinical metadata from 1,095 TCGA breast cancer patients. After removing outliers, HR+ samples were ordered based on their HOXB13 expression levels and the top 10% of samples (HOXB13-high; n=95) and bottom 10% (HOXB13-low; n=95) were selected for further analysis. Differentially expressed genes (DEGs) were determined by comparing all protein-coding genes between the HOXB13-high and HOXB13-low sets, based on an absolute log-fold change above 2 and a Benjamini-Hochberg adjusted p-value threshold of 0.01. Among the subset of samples with matching methylation data (n=60 for the HOXB13-high group and n=48 for the HOXB13-low group), hypo-methylated and hyper-methylated differentially methylated probes (DMPs) were identified based on at least 10% change of β (i.e., average promoter methylation) and a Benjamini-Hochberg adjusted p-value threshold of 0.01. Motif and transcription factor (TF) enrichment analyses were identified at a minimum incidence of 10 and a lower boundary of 1.5 for the 95% confidence interval of the odds ratio. Results: Our analysis identified a total of 613 DEGs between the HOXB13-high and HOXB13-low groups. Gene ontology analysis of DEGs revealed a statistically significant enrichment of genes associated with PTK6 regulated cell cycle, oncogene-induced cellular senescence, and regulation of T cell activation. Analysis of protein-protein interactions of differentially expressed genes revealed differences in key biological processes including changes in cell cycle regulation, p53 pathway, and IL-17 signaling. A total of 5,295 hyper-methylated and 580 hypo-methylated DMPs were discovered when comparing HOXB13-high and HOXB13-low samples. Changes in global methylation patterns enabled classification of breast cancer samples with high or low HOXB13 expression with an accuracy of 0.90. Reconstruction of gene regulatory networks from DNA methylation and transcriptome profiles elucidated significant pairs of DMPs and differentially expressed genes. Specifically, a total of 197 differentially expressed genes were attributed to hyper-methylated DMPs, including genes involved in response to estrogen, activation of HOX genes during differentiation, and drug-mediated inhibition of CDK4/6 activity. Additionally, 8 differentially expressed genes were attributed to hypo-methylated DMPs and were associated with proteasomal protein catabolic processes. Lastly, motif analysis of transcriptional factors for the hyper-methylated DMPs revealed enrichment of binding motifs attributed to FOXA1, ESR1, XBR1, and FOXP1. Conclusion: A comprehensive and comparative analysis of mRNA expression levels and methylation patterns between TCGA breast cancer samples with high and low expression of HOXB13 revealed key signaling pathways and biological processes that may provide insights on the molecular mechanism of HOXB13 in the regulation of response to endocrine therapy in HR+ breast cancer. Citation Format: Natalia Siuliukina, Yi Zhang, Kai Treuner, Mark Pegram. Differential HOXB13 gene expression and promoter methylation analysis in breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO2-01-08.

S100a9 might act as a modulator of the Toll-like receptor 4 transduction pathway in chronic rhinosinusitis with nasal polyps

Chronic rhinosinusitis with nasal polyp (CRSwNP) is a highly prevalent disorder characterized by persistent nasal and sinus mucosa inflammation. Despite significant morbidity and decreased quality of life, there are limited effective treatment options for such a disease. Therefore, identifying causal genes and dysregulated pathways paves the way for novel therapeutic interventions. In the current study, a three-way interaction approach was used to detect dynamic co-expression interactions involved in CRSwNP. In this approach, the internal evolution of the co-expression relation between a pair of genes (X, Y) was captured under a change in the expression profile of a third gene (Z), named the switch gene. Subsequently, the biological relevancy of the statistically significant triplets was confirmed using both gene set enrichment analysis and gene regulatory network reconstruction. Finally, the importance of identified switch genes was confirmed using a random forest model. The results suggested four dysregulated pathways in CRSwNP, including “positive regulation of intracellular signal transduction”, “arachidonic acid metabolic process”, “spermatogenesis” and “negative regulation of cellular protein metabolic process”. Additionally, the S100a9 as a switch gene together with the gene pair {Cd14, Tpd52l1} form a biologically relevant triplet. More specifically, we suggested that S100a9 might act as a potential upstream modulator in toll-like receptor 4 transduction pathway in the major CRSwNP pathologies.

Organ-delimited gene regulatory networks provide high accuracy in candidate transcription factor selection across diverse processes

Organ-specific gene expression datasets that include hundreds to thousands of experiments allow the reconstruction of organ-level gene regulatory networks (GRNs). However, creating such datasets is greatly hampered by the requirements of extensive and tedious manual curation. Here, we trained a supervised classification model that can accurately classify the organ-of-origin for a plant transcriptome. This K-Nearest Neighbor-based multiclass classifier was used to create organ-specific gene expression datasets for the leaf, root, shoot, flower, and seed in Arabidopsis thaliana. A GRN inference approach was used to determine the: i. influential transcription factors (TFs) in each organ and, ii. most influential TFs for specific biological processes in that organ. These genome-wide, organ-delimited GRNs (OD-GRNs), recalled many known regulators of organ development and processes operating in those organs. Importantly, many previously unknown TF regulators were uncovered as potential regulators of these processes. As a proof-of-concept, we focused on experimentally validating the predicted TF regulators of lipid biosynthesis in seeds, an important food and biofuel trait. Of the top 20 predicted TFs, eight are known regulators of seed oil content, e.g., WRI1, LEC1, FUS3. Importantly, we validated our prediction of MybS2, TGA4, SPL12, AGL18, and DiV2 as regulators of seed lipid biosynthesis. We elucidated the molecular mechanism of MybS2 and show that it induces purple acid phosphatase family genes and lipid synthesis genes to enhance seed lipid content. This general approach has the potential to be extended to any species with sufficiently large gene expression datasets to find unique regulators of any trait-of-interest.

BHLHE40/41 regulate microglia and peripheral macrophage responses associated with Alzheimer’s disease and other disorders of lipid-rich tissues

Genetic and experimental evidence suggests that Alzheimer’s disease (AD) risk alleles and genes may influence disease susceptibility by altering the transcriptional and cellular responses of macrophages, including microglia, to damage of lipid-rich tissues like the brain. Recently, sc/nRNA sequencing studies identified similar transcriptional activation states in subpopulations of macrophages in aging and degenerating brains and in other diseased lipid-rich tissues. We collectively refer to these subpopulations of microglia and peripheral macrophages as DLAMs. Using macrophage sc/nRNA-seq data from healthy and diseased human and mouse lipid-rich tissues, we reconstructed gene regulatory networks and identified 11 strong candidate transcriptional regulators of the DLAM response across species. Loss or reduction of two of these transcription factors, BHLHE40/41, in iPSC-derived microglia and human THP-1 macrophages as well as loss of Bhlhe40/41 in mouse microglia, resulted in increased expression of DLAM genes involved in cholesterol clearance and lysosomal processing, increased cholesterol efflux and storage, and increased lysosomal mass and degradative capacity. These findings provide targets for therapeutic modulation of macrophage/microglial function in AD and other disorders affecting lipid-rich tissues.

Population-level comparisons of gene regulatory networks modeled on high-throughput single-cell transcriptomics data.

Single-cell technologies enable high-resolution studies of phenotype-defining molecular mechanisms. However, data sparsity and cellular heterogeneity make modeling biological variability across single-cell samples difficult. Here we present SCORPION, a tool that uses a message-passing algorithm to reconstruct comparable gene regulatory networks from single-cell/nuclei RNA-sequencing data that are suitable for population-level comparisons by leveraging the same baseline priors. Using synthetic data, we found that SCORPION outperformed 12 existing gene regulatory network reconstruction techniques. Using supervised experiments, we show that SCORPION can accurately identify differences in regulatory networks between wild-type and transcription factor-perturbed cells. We demonstrate SCORPION's scalability to population-level analyses using a single-cell RNA-sequencing atlas containing 200,436 cells from colorectal cancer and adjacent healthy tissues. The differences between tumor regions detected by SCORPION are consistent across multiple cohorts as well as with our understanding of disease progression, and elucidate phenotypic regulators that may impact patient survival.

SPREd: a simulation-supervised neural network tool for gene regulatory network reconstruction.

Reconstruction of gene regulatory networks (GRNs) from expression data is a significant open problem. Common approaches train a machine learning (ML) model to predict a gene's expression using transcription factors' (TFs') expression as features and designate important features/TFs as regulators of the gene. Here, we present an entirely different paradigm, where GRN edges are directly predicted by the ML model. The new approach, named "SPREd," is a simulation-supervised neural network for GRN inference. Its inputs comprise expression relationships (e.g. correlation, mutual information) between the target gene and each TF and between pairs of TFs. The output includes binary labels indicating whether each TF regulates the target gene. We train the neural network model using synthetic expression data generated by a biophysics-inspired simulation model that incorporates linear as well as non-linear TF-gene relationships and diverse GRN configurations. We show SPREd to outperform state-of-the-art GRN reconstruction tools GENIE3, ENNET, PORTIA, and TIGRESS on synthetic datasets with high co-expression among TFs, similar to that seen in real data. A key advantage of the new approach is its robustness to relatively small numbers of conditions (columns) in the expression matrix, which is a common problem faced by existing methods. Finally, we evaluate SPREd on real data sets in yeast that represent gold-standard benchmarks of GRN reconstruction and show it to perform significantly better than or comparably to existing methods. In addition to its high accuracy and speed, SPREd marks a first step toward incorporating biophysics principles of gene regulation into ML-based approaches to GRN reconstruction. Data and code are available from https://github.com/iiiime/SPREd.

Discovery of novel plant biomass conversion associated fungal transcription factors using a network-based approach

Fungal plant biomass conversion (FPBC) is an important component of the global carbon cycle and has been widely applied for the production of biofuels, enzymes and biochemicals. Identification of transcription factors (TFs) governing FPBC is crucial for genetic engineering of industrial fungi towards sustainable production of high-value bioproducts from renewable lignocellulose. Here, we developed a bioinformatics framework for the identification of FPBC related TFs based on reconstructed gene regulatory networks and enrichment analysis of manually curatedFPBC gene sets. Applying this approach to model fungi Aspergillus niger and Neurospora crassa, we successfully identified both known TFs and promising candidates. The function of one identified TF, HapX, has been experimentally validated, and several candidates were supported by literature, transcriptome data and initial growth analysis. Our new approach will accelerate the identification of novel TFs involved in FPBC, and facilitate the further improvement of fungal cell factories.

A review on gene regulatory network reconstruction algorithms based on single cell RNA sequencing.

Understanding gene regulatory networks (GRNs) is essential for unraveling the molecular mechanisms governing cellular behavior. With the advent of high-throughput transcriptome measurement technology, researchers have aimed to reverse engineer the biological systems, extracting gene regulatory rules from their outputs, which represented by gene expression data. Bulk RNA sequencing, a widely used method for measuring gene expression, has been employed for GRN reconstruction. However, it falls short in capturing dynamic changes in gene expression at the level of individual cells since it averages gene expression across mixed cell populations. In this review, we provide an overview of 15 GRN reconstruction tools and discuss their respective strengths and limitations, particularly in the context of single cell RNA sequencing (scRNA-seq). Recent advancements in scRNA-seq break new ground of GRN reconstruction. They offer snapshots of the individual cell transcriptomes and capturing dynamic changes. We emphasize how these technological breakthroughs have enhanced GRN reconstruction. GRN reconstructors can be classified based on their requirement for cellular trajectory, which represents a dynamical cellular process including differentiation, aging, or disease progression. Benchmarking studies support the superiority of GRN reconstructors that do not require trajectory analysis in identifying regulator-target relationships. However, methods equipped with trajectory analysis demonstrate better performance in identifying key regulatory factors. In conclusion, researchers should select a suitable GRN reconstructor based on their specific research objectives.