Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Expression systems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

MEPHESTO: Modeling Energy-Performance in Heterogeneous SoCs and Their Trade-Offs

Integrated shared memory heterogeneous architectures are pervasive because they satisfy the diverse needs of mobile, autonomous, and edge computing platforms. Although specialized processing units (PUs) that share a unified system memory improve performance and energy efficiency by reducing data movement, they also increase contention for this memory since the PUs interact with each other. Prior work has investigated performance degradation due to memory contention, but few have studied the relationship of power and energy to memory contention. Moreover, a comprehensive solution that models memory contention for kernel placement on contemporary heterogeneous systems on chip (SoCs) in response to energy and performance has been largely unaddressed.This paper presents MEPHESTO, a novel and holistic approach for managing this balance. The authors characterize applications and PUs in terms of two memory contention factors - time factors and power factors - to achieve the desired trade-off between energy and performance for collocated kernel execution on heterogeneous systems. The authors believe that this investigation is the first to combine all of these factors and present a simple knob-based approach that expresses the target trade-off. The approach is evaluated on a diverse integrated shared memory heterogeneous system with a CPU, GPU, and programmable vision accelerator. By using an empirical model for memory contention that provides up to 92% accuracy, the kernel collocation approach can provide a near-optimal ordering and placement based on the user-defined, energy-performance trade-off parameter. Moreover, the dynamic programming-based heuristics provide up to 30% better energy or 20% performance benefits when compared with the greedy approaches commonly employed by previous studies.

Alaul haque monil, Mohammad↗

Synergy and antagonism in a genome-scale model of metabolic hijacking by bacteriophages

Bacteriophage auxiliary metabolic genes (AMGs) alter host metabolism by hijacking reactions, but previous studies mostly inferred their roles from annotations, ignoring system-wide impacts and phage production. Here we integrate AMGs and phage assembly into a genome-scale metabolic model of Prochloroccocus marinus MED4 infected by P-HM2. We show that 17 directly hijacked reactions substantially affect more than 30% of the reactions in MED4 metabolism, including carbon fixation, photosynthesis, and nucleotide synthesis, distinguishing these AMGs as either phage aligned—shifting feasible reaction velocities in accordance with maximal phage production—or phage antialigned. Pareto optimization reveals that phage-aligned reactions alter phage-host growth trade-offs, while phage-antialigned reactions do not. We experimentally validate our predictions of system-level AMG impacts by measuring the N-dependent effect of P-HM2 cp12 expression on growth in a model relative of the genetically intractable MED4, Synechococcus elongatus. We also show that AMGs’ indirect impacts are synergistically and antagonistically coupled, providing systems-level insight into AMG perturbations and highlighting how nontrivial cascading effects shape host metabolism.

Rozum, Jordan C. [Pacific Northwest National Labor↗

Metabolic multistability and hysteresis in a model aerobe-anaerobe microbiome community

Major changes in the microbiome are associated with health and disease. Some microbiome states persist despite seemingly unfavorable conditions, such as the proliferation of aerobe-anaerobe communities in oxygen-exposed environments in wound infections or small intestinal bacterial overgrowth. Mechanisms underlying transitions into and persistence of these states remain unclear. Using two microbial taxa relevant to the human microbiome, we combine genome-scale mathematical modeling, bioreactor experiments, transcriptomics, and dynamical systems theory to show that multistability and hysteresis (MSH) is a mechanism describing the shift from an aerobe-dominated state to a resilient, paradoxically persistent aerobe-anaerobe state. We examine the impact of changing oxygen and nutrient regimes and identify changes in metabolism and gene expression that lead to MSH and associated multi-stable states. In such systems, conceptual causation-correlation connections break and MSH must be used for analysis. Using MSH to analyze microbiome dynamics will improve our conceptual understanding of stability of microbiome states and transitions between states.

59 BASIC BIOLOGICAL SCIENCES↗

The role of CenKR in the coordination of Rhodobacter sphaeroides cell elongation and division

ABSTRACT Cell elongation and division are essential aspects of the bacterial life cycle that must be coordinated for viability and replication. The impact of misregulation of these processes is not well understood as these systems are often not amenable to traditional genetic manipulation. Recently, we reported on the CenKR two-component system (TCS) in the Gram-negative bacterium Rhodobacter sphaeroides that is genetically tractable, widely conserved in α-proteobacteria, and directly regulates the expression of components crucial for cell elongation and division, including genes encoding subunit of the Tol-Pal complex. In this work, we show that overexpression of cenK results in cell filamentation and chaining. Using cryo-electron microscopy (cryo-EM) and cryo-electron tomography (cryo-ET), we generated high-resolution two-dimensional (2D) images and three-dimensional (3D) volumes of the cell envelope and division septum of wild-type cells and a cenK overexpression strain finding that these morphological changes stem from defects in outer membrane (OM) and peptidoglycan (PG) constriction. By monitoring the localization of Pal, PG biosynthesis, and the bacterial cytoskeletal proteins MreB and FtsZ, we developed a model for how increased CenKR activity leads to changes in cell elongation and division. This model predicts that increased CenKR activity decreases the mobility of Pal, delaying OM constriction, and ultimately disrupting the midcell positioning of MreB and FtsZ and interfering with the spatial regulation of PG synthesis and remodeling. IMPORTANCE By coordinating cell elongation and division, bacteria maintain their shape, support critical envelope functions, and orchestrate division. Regulatory and assembly systems have been implicated in these processes in some well-studied Gram-negative bacteria. However, we lack information on these processes and their conservation across the bacterial phylogeny. In R. sphaeroides and other α-proteobacteria, CenKR is an essential two-component system (TCS) that regulates the expression of genes known or predicted to function in cell envelope biosynthesis, elongation, and/or division. Here, we leverage unique features of CenKR to understand how increasing its activity impacts cell elongation/division and use antibiotics to identify how modulating the activity of this TCS leads to changes in cell morphology. Our results provide new insight into how CenKR activity controls the structure and function of the bacterial envelope, the localization of cell elongation and division machinery, and cellular processes in organisms with importance in health, host-microbe interactions, and biotechnology.

59 BASIC BIOLOGICAL SCIENCES↗

Language-Theoretic Data Analysis to Support ICS Protocol Baselining

Critical infrastructure stakeholders need to baseline their systems to understand expected protocol communications. Baseline behaviors may vary based on operational context. Expected operations during a maintenance window, for example, may be different from normal operations. Furthermore, constructing system baselines for Industrial Control Systems (ICS) is difficult and time-consuming. ICS processes generate artifacts expressed across heterogeneous data sources such as network traffic and device logs. This paper explores the hypothesis that such ICS artifacts form a language in the language-theoretic sense. From a theoretical perspective, the variety of implementations of ICS protocols and constrained environment of OT networks provide a rich application domain for language-theoretic approaches. We present several use cases related to the practical construction of system baselines: grammars for data fusion, language dialects for device fingerprinting, and security automata for system baselining

24 POWER TRANSMISSION AND DISTRIBUTION↗

γ-Aminobutyric acid modulates terpene biosynthesis through the ATG8a-mediated pathway

Terpenes play crucial roles in plant growth, development, and stress responses. The biosynthesis of terpenes is influenced by abiotic stress factors, such as drought, temperature, or light. Here, however, the molecular network underlying how terpenes are regulated in response to environmental stimuli remains largely unknown. Here, we identified the autophagy protein SgATG8a as a key mediator of GABA-regulated terpene production and drought tolerance in Sindora glabra. SgATG8a, evolutionarily related to the animal GABA receptor-associated protein (GABARAP) subfamily, localizes in both the nucleus and cytoplasm. Exogenous GABA treatment not only increased the expression level of terpene synthase genes (SgTPSs) but also led to enhanced accumulation of six main terpene components in Sindora glabra. In addition, GABA alleviated the photosynthesis damage and enhanced leaf biomass under drought conditions. Consistently, overexpression of SgATG8a in Arabidopsis increased terpene synthase gene (SgTPS) expression, leading to the enhanced production of four major terpenes and improved the tolerance of transgenic plants to drought stress by regulating reactive oxygen species (ROS) scavenging systems. Moreover, the transcription factors SgWRKY13 and SgERF4 were identified as interacting partners of SgATG8a, activating SgTPS3 expression. Lectin receptor-like kinase (LecRK1) is involved in the GABA-mediated pathway by interacting with the SgWRKY13/SgERF4-SgATG8a proteins, and the LecRK1-SgWRKY13/SgERF4 phosphorylation module fine-tunes the transcription of the downstream SgTPS3 gene. Taken together, these findings reveal a novel role for GABA in regulating terpene biosynthesis and drought tolerance, providing insights into the molecular mechanism underlying GABA-mediated terpene production.

59 BASIC BIOLOGICAL SCIENCES↗

Too dim, too bright, and just right: Systems analysis of the Chlamydomonas diurnal program under limiting and excess light

Photosynthetic organisms coordinate their metabolism and growth with diurnal light, which can range in intensity from limiting to excessive. Little is known about how light intensity impacts the diurnal program in Chlamydomonas reinhardtii, or how diurnal rhythms in gene expression and metabolism shape photoprotective responses at different times of day. To address these questions, we performed a systems analysis of synchronized Chlamydomonas populations acclimated to low, moderate, and high diurnal light. Transcriptomic and proteomic data revealed that the Chlamydomonas rhythmic gene expression program is resilient to limiting and excess light: genome-wide, waves of transcripts, and proteins peak at the same times in populations acclimated to stressful light intensities as in populations acclimated to moderate light. Yet, diurnal photoacclimation gives rise to hundreds of gene expression changes, even at night. Time course measurements of photosynthetic efficiency and pigments responsive to excess light showed that high light-acclimated cells partially overcome photodamage in the latter half of the day prior to cell division. Although gene expression and photodamage are dynamic over the diurnal cycle, Chlamydomonas populations acclimated to low and high diurnal light maintain altered photosystem abundance, thylakoid architecture, and non-photochemical quenching capacity through the night phase. This suggests that cells remember or anticipate the light intensities that they have typically encountered during the day. The integrated data constitute an excellent resource for understanding photoacclimation in eukaryotes under environmentally relevant conditions.

Dupuis, Sunnyjoy↗

Engineering quorum-sensing circuits in Synechococcus elongatus PCC 7942 towards self-inducible systems

Despite significant potential for cyanobacteria as sustainable bioproduction chases, there are limited examples of scaled cyanobacterial bioproduction. In part, this is because most cyanobacterial species are poorly adapted to bioreactor cultivation conditions and lack features that facilitate biomass growth and harvesting at scale. We explored quorum sensing (QS) pathways derived from heterotrophic microbes as a method for autoinduction of gene expression circuits coordinated to population density in cyanobacteria. Here, we integrated genetic modules designed to produce and detect the diffusible QS signal, acyl-homoserine lactones (AHLs), in the cyanobacterial model, Synechococcus elongatus PCC 7942 (S. elongatus). We demonstrate that S. elongatus heterologously produces sufficient AHL signals to activate gene expression in a dose-dependent and population density-responsive manner. A hybrid combination of AHL synthesis enzyme from Vibrio fischeri (Lux system) with the transcription factor receiver from Pseudomonas aeruginosa (Las system) provides an ideal activation ratio and mitigates toxicity observed with some AHL systems. As a proof of concept, we coupled the QS pathway to the expression of a cell division inhibitory gene, cdv3, facilitating late-phase cell elongation, cell sedimentation, and improved biomass recovery. Our findings provide a foundation for the development of auto-induction systems leverageable to improve cyanobacterial biotechnology applications.

Acyl homoserine lactones (AHL)↗

Gapless Criterion for Crystals from Effective Axion Field

Gapless criteria that can efficiently determine whether a crystal is gapless or not are particularly useful for identifying topological semimetals. In this work, we propose a sufficient gapless criterion for three-dimensional noninteracting crystals, based on the simplified expressions for the bulk average value of the static axion field. The brief logic is that two different simplified expressions give the same value in an insulator, and thus the gapless phase can be detected by the mismatch of them. We demonstrate the effectiveness of the gapless criterion in the magnetic systems with space groups 26 and 13, where mirror, glide, and inversion symmetries provide the simplified expressions. In particular, the gapless criterion can identify gapless phases that are missed by the symmetry-representation approach, as illustrated by space group 26. Our proposal serves as a guiding principle for future discovery of topological semimetals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Improving Arsenic Tolerance of Pyrococcus furiosus by Heterologous Expression of a Respiratory Arsenate Reductase

Arsenate is a notorious toxicant that is known to disrupt multiple biochemical pathways. Many microorganisms have developed mechanisms to detoxify arsenate using the ArsC-type arsenate reductase, and some even use arsenate as a terminal electron acceptor for respiration involving arsenate respiratory reductase (Arr). ArsC-type reductases have been studied extensively, but the phylogenetically unrelated Arr system is less investigated and has not been characterized from Archaea. Here, we heterologously expressed the genes encoding Arr from the crenarchaeon Pyrobaculum aerophilum in the euryarchaeon Pyrococcus furiosus, both of which grow optimally near 100°C. Recombinant P. furiosus was grown on molybdenum (Mo)- or tungsten (W)-containing medium, and two types of recombinant Arr enzymes were purified, one containing Mo (Arr-Mo) and one containing W (Arr-W). Purified Arr-Mo had a 140-fold higher specific activity in arsenate [As(V)] reduction than Arr-W, and Arr-Mo also reduced arsenite [As(III)]. The P. furiosus strain expressing Arr-Mo (the Arr strain) was able to use arsenate as a terminal electron acceptor during growth on peptides. In addition, the Arr strain had increased tolerance compared to that of the parent strain to arsenate and also, surprisingly, to arsenite. Compared to the parent, the Arr strain accumulated intracellularly almost an order of magnitude more arsenic when cells were grown in the presence of arsenite. Furthermore, x-ray absorption spectroscopy (XAS) results suggest that the Arr strain of P. furiosus improves its tolerance to arsenite by increasing production of less-toxic arsenate and nontoxic methylated arsenicals compared to that by the parent.

59 BASIC BIOLOGICAL SCIENCES↗

The changing mouse embryo transcriptome at whole tissue and single-cell resolution

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system 1 . Here we systematically quantified mouse polyA-RNA from day 10.5 of embryonic development to birth, sampling 17 tissues and organs. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets that were further characterized by the transcription factor motif codes of their promoters. We decomposed the tissue-level transcriptome using single-cell RNA-seq (sequencing of RNA reverse transcribed into cDNA) and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for one-third of differential gene expression and more than 40% of identified cell types. By integrating promoter sequence motifs with companion ENCODE epigenomic profiles, we identified a prominent promoter de-repression mechanism in neuronal expression clusters that was attributable to known and novel repressors. Focusing on the developing limb, single-cell RNA data identified 25 candidate cell types that included progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell-type transcription factor networks and complementary sets of candidate enhancer elements by using single-cell RNA-seq to decompose integrative cis -element (IDEAS) models that were derived from whole-tissue epigenome chromatin data. These ENCODE reference data, computed network components and IDEAS chromatin segmentations are companion resources to the matching epigenomic developmental matrix, and are available for researchers to further mine and integrate.

59 BASIC BIOLOGICAL SCIENCES↗

A Mathematical Model of the Dynamics of Cytokine Expression and Human Immune Cell Activation in Response to the Pathogen Staphylococcus aureus

Cell-based mathematical models have previously been developed to simulate the immune system in response to pathogens. Mathematical modeling papers which study the human immune response to pathogens have predicted concentrations of a variety of cells, including activated and resting macrophages, plasma cells, and antibodies. This study aims to create a comprehensive mathematical model that can predict cytokine levels in response to a gram-positive bacterium, S. aureus by coupling previous models. To accomplish this, the cytokines Tumor Necrosis Factor Alpha (TNF-α), Interleukin 6 (IL-6), Interleukin 8 (IL-8), and Interleukin 10 (IL-10) are included to quantify the relationship between cytokine release from macrophages and the concentration of the pathogen, S. aureus, ex vivo. Partial differential equations (PDEs) are used to model cellular response and ordinary differential equations (ODEs) are used to model cytokine response, and interactions between both components produce a more robust and more complete systems-level understanding of immune activation. In the coupled cellular and cytokine model outlined in this paper, a low concentration of S. aureus is used to stimulate the measured cellular response and cytokine expression. Results show that our cellular activation and cytokine expression model characterizing septic conditions can predict ex vivo mechanisms in response to gram-negative and gram-positive bacteria. Our simulations provide new insights into how the human immune system responds to infections from different pathogens. Novel applications of these insights help in the development of more powerful tools and protocols in infection biology.

60 APPLIED LIFE SCIENCES↗

A hybrid cyt c maturation system enhances the bioelectrical performance of engineered Escherichia coli by improving the rate-limiting step

Bioelectronic devices can use electron flux to enable communication between biotic components and abiotic electrodes. We have modified Escherichia coli to electrically interact with electrodes by expressing the cytochrome c from Shewanella oneidensis MR-1. However, we observe inefficient electrical performance, which we hypothesize is due to the limited compatibility of the E. coli cytochrome c maturation (Ccm) systems with MR-1 cytochrome c. Here we test whether the bioelectronic performance of E. coli can be improved by constructing hybrid Ccm systems containing protein domains from both E. coli and S. oneidensis MR-1. The hybrid CcmH increased cytochrome c expression by increasing the abundance of CymA 60%, while only slightly changing the abundance of the other cytochromes c. Electrochemical measurements showed that the overall current from the hybrid ccm strain increased 121% relative to the wildtype ccm strain, with an electron flux per cell of 12.3 ± 0.3 fA·cell -1 . Additionally, the hybrid ccm strain doubled its electrical response with the addition of exogenous flavin, and quantitative analysis of this demonstrates CymA is the rate-limiting step in this electron conduit. These results demonstrate that this hybrid Ccm system can enhance the bioelectrical performance of the cyt c expressing E. coli, allowing the construction of more efficient bioelectronic devices.

59 BASIC BIOLOGICAL SCIENCES↗

SolarAPP+ Performance Review (2024 Data)

The Solar Automated Permit Processing Plus (SolarAPP+) platform is an online portal to facilitate and expedite rooftop solar photovoltaic (PV) and battery storage permitting processes. SolarAPP+ allows PV contractors to upload system specifications, have that information automatically reviewed for code compliance, and receive instant approval for code-compliant systems, reducing authority having jurisdiction (AHJ) staff time needed for review. SolarAPP+ also provides inspection checklists to verify installation practices and adherence to approved designs. This report is part of an ongoing series of reviews of SolarAPP+ performance. Consistent with previous performance reviews, we summarize SolarAPP+ adoption trends to date and compare various metrics for PV systems permitted through SolarAPP+ versus systems permitted through traditional AHJ permitting processes. As of the end of 2024, 799 AHJs had expressed interest in the platform, with 264 fully adopting (215) or piloting (49) the platform. In 2024, 861 installers submitted 37,393 permits through the SolarAPP+ platform, including 27,375 permits for PV+storage systems. SolarAPP+ permits accounted for around 43% of all permits issued in all participating AHJs, and more than 60% of all permits in several participating AHJs. We compare permitting timelines through SolarAPP+ to traditional AHJ permitting processes to assess the platform's performance. Consistent with previous SolarAPP+ performance reviews, we find that permitting timelines are significantly shorter for SolarAPP+ projects. Based on median timelines, a typical SolarAPP+ project is permitted and inspected 12 business days sooner than traditional projects. We estimate that automatic SolarAPP+ permitting saved around 18,400 hours of AHJ staff time in 2024. Finally, we estimate that SolarAPP+ eliminated over 100,000 business days in permitting-related delays in 2024.

14 SOLAR ENERGY↗

A KBase Case Study on Genome-wide Transcriptomics and Plant Primary Metabolism in Sorghum

A better understanding of the genetic and metabolic mechanisms that confer stress resistance and tolerance in plants is key to engineering new crops through advanced breeding technologies. This requires a systems biology approach that builds on a genome-wide understanding of the regulation of gene expression, plant metabolism, physiology and growth. In this study, we examine the response to drought stress in Sorghum, as we leverage the tools for transcriptomics and plant metabolic modeling we have implemented at the U.S. Department of Energy Systems Biology Knowledgebase (KBase). KBase enables researchers worldwide to collaborate and advance research by allowing them to upload private or public data into the KBase Narrative Interface, empowering them to analyze this data using a rich, extensible array of computational and data-analytics tools, and allows them to securely share scientific workflows and conclusions. We demonstrate how to use the current RNA-seq tools in KBase, applicable to both plants and microbes, to assemble and quantify long transcripts and identify differentially expressed genes effectively. More specifically, we demonstrate the utility of the platform by identifying key genes that are differentially expressed during drought-stress in Sorghum bicolor, which is an important sustainable production crop plant. We then show how to use KBase tools to predict the membership of genes in metabolic pathways and examine expression data in the context of metabolic subsystems. We demonstrate the power of the platform by making the data analysis and interpretation available to the biologists in the reproducible, re-usable, point-and-click format of a KBase Narrative thus promoting FAIR (Findable, Accessible, Interoperable and Reusable) guiding principles for scientific data management and stewardship.

59 BASIC BIOLOGICAL SCIENCES↗

Local SU(2) × U(1) quark flavor symmetry in the RS bulk

We propose a model of quark flavor based on an additional SU(2) × U(1) local symmetry in a warped extra dimensional bulk. In contrast to other works, we break the additional gauge symmetry in the bulk via two complex scalars which acquire bulk vevs, rather than relying on brane-localized symmetry breaking. A gauge-covariant Kaluza-Klein decomposition of a theory with a bulk spontaneously broken gauge symmetry is performed, and exact expressions for the bulk profiles of all physical particles in such systems are given. The SM quark masses and mixings are then recreated using gauge-covariant bulk quark mass terms and Yukawa-like couplings to the new bulk scalars. A numerical sampling of points in the model parameter space that recreate the quark masses and mixings is performed at a KK scale of M KK = 5 TeV. We then compute the ΔF = 2 4-quark operators arising from our new flavor gauge bosons and scalars, and those arising from Kaluza-Klein modes of SM gauge bosons. By decoupling one of our bulk scalar fields to all quark fields except the right-handed up-like sector, we find that it is possible to greatly suppress tree-level contributions to the highly constrained Kaon mixing parameters. Instead, the dominant constraints on the model emerge from neutral Bd and D meson mixing. These constraints are explored with our numerical sampling of the model parameter space, and the specific contribution of the new flavor gauge bosons and scalars is discussed. We find that for a significant range of realistic flavor gauge couplings, the new gauge bosons compete with the normally dominant gluon flavor-changing currents, but flavor-changing operators emerging from the bulk scalar fields are highly suppressed. Finally, we briefly comment on flavor constraints that are independent of the flavor gauge sector arising from the Z$\overline{b}$ L b L coupling and rare top decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The effective number of parameters in kernel density estimation

We devise a new formula for measuring the effective degrees of freedom (EDoF) in kernel density estimation (KDE). Starting from the orthogonal polynomial sequence (OPS) expansion for the ratio of the empirical to the oracle density, we show how convolution with the kernel leads to a new OPS with respect to which one may express the resulting KDE. The expansion coefficients of the two OPS systems can then be related via a kernel sensitivity matrix, which leads to a natural oracle definition of EDoF through the trace operator. Asymptotic properties of the (empirical) plug-in EDoF are worked out through influence functions, and connections with other empirical EDoFs are established. Minimization of Kullback-Leibler divergence is investigated as an alternative to integrated squared error based bandwidth selection rules, yielding a new normal scale rule. The methodology, which arises from a proper oracle formulation and is not restricted to convolution kernels, suggests the possibility of a new bandwidth selection rule based on an information criterion such as AIC.

bandwidth selection↗

Identifying strengths and weaknesses of methods for computational network inference from single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) offers unparalleled insight into the transcriptional programs of different cellular states by measuring the transcriptome of thousands of individual cells. An emerging problem in the analysis of scRNA-seq is the inference of transcriptional gene regulatory networks and a number of methods with different learning frameworks have been developed to address this problem. Here, we present an expanded benchmarking study of eleven recent network inference methods on seven published scRNA-seq datasets in human, mouse, and yeast considering different types of gold standard networks and evaluation metrics. We evaluate methods based on their computing requirements as well as on their ability to recover the network structure. We find that, while most methods have a modest recovery of experimentally derived interactions based on global metrics such as Area Under the Precision Recall curve, methods are able to capture targets of regulators that are relevant to the system under study. Among the top performing methods that use only expression were SCENIC, PIDC, MERLIN or Correlation. Addition of prior biological knowledge and the estimation of transcription factor activities resulted in the best overall performance with the Inferelator and MERLIN methods that use prior knowledge outperforming methods that use expression alone. We found that imputation for network inference did not improve network inference accuracy and could be detrimental. Comparisons of inferred networks for comparable bulk conditions showed that the networks inferred from scRNA-seq datasets are often better or at par with the networks inferred from bulk datasets. Our analysis should be beneficial in selecting methods for network inference. At the same time, this highlights the need for improved methods and better gold standards for regulatory network inference from scRNAseq datasets.

59 BASIC BIOLOGICAL SCIENCES↗