Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gene regulatory network”

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 55 records · Page 3

Engineering activatable promoters for scalable and multi-input CRISPRa/i circuits

Dynamic, multi-input gene regulatory networks (GRNs) are ubiquitous in nature. Multilayer CRISPR-based genetic circuits hold great promise for building GRNs akin to those found in naturally occurring biological systems. We develop an approach for creating high-performing activatable promoters that can be assembled into deep, wide, and multi-input CRISPR-activation and -interference (CRISPRa/i) GRNs. By integrating sequence-based design and in vivo screening, we engineer activatable promoters that achieve up to 1,000-fold dynamic range in an Escherichia coli-based cell-free system. These components enable CRISPRa GRNs that are six layers deep and four branches wide. We show the generalizability of the promoter engineering workflow by improving the dynamic range of the light-dependent EL222 optogenetic system from 6-fold to 34-fold. Additionally, high dynamic range promoters enable CRISPRa systems mediated by small molecules and protein–protein interactions. We apply these tools to build input-responsive CRISPRa/i GRNs, including feedback loops, logic gates, multilayer cascades, and dynamic pulse modulators. Our work provides a generalizable approach for the design of high dynamic range activatable promoters and enables classes of gene regulatory functions in cell-free systems.

42 ENGINEERING↗

Differential Etv2 threshold requirement for endothelial and erythropoietic development

Endothelial and erythropoietic lineages arise from a common developmental progenitor. Etv2 is a master transcriptional regulator required for the development of both lineages. However, the mechanisms through which Etv2 initiates the gene-regulatory networks (GRNs) for endothelial and erythropoietic specification and how the two GRNs diverge downstream of Etv2 remain incompletely understood. Here, by analyzing a hypomorphic Etv2 mutant, we demonstrate different threshold requirements for initiation of the downstream GRNs for endothelial and erythropoietic development. We show that Etv2 functions directly in a coherent feedforward transcriptional network for vascular endothelial development, and a low level of Etv2 expression is sufficient to induce and sustain the endothelial GRN. In contrast, Etv2 induces the erythropoietic GRN indirectly via activation of Tal1, which requires a significantly higher threshold of Etv2 to initiate and sustain erythropoietic development. These results provide important mechanistic insight into the divergence of the endothelial and erythropoietic lineages.

59 BASIC BIOLOGICAL SCIENCES↗

Single-cell and spatial omics in plants: from cellular atlases to regulatory mechanisms

Single-cell RNA sequencing (scRNA-seq) has transformed transcriptomic studies by enabling gene expression profiling at the resolution of individual cells within and across a broad range of tissue types, revealing cellular heterogeneity that is obscured in bulk tissue transcriptomes. Over the past decade, improvements in microfluidics and library preparation have drastically increased throughput, allowing tens of thousands of cells to be assayed in a single experiment. Although initially developed in animal systems, scRNA-seq has rapidly emerged as a powerful and widely adopted approach in plant biology. Beyond transcriptomics, the integration of single-cell data with chromatin accessibility, proteomics, metabolomics, and spatial omics is enabling a system-level understanding of plant gene regulation and cellular organization. Network-based analytical frameworks further support the reconstruction of gene regulatory networks and the interpretation of complex single-cell data. In this review, we summarize the current technological landscape of plant single-cell studies, discuss key experimental and analytical challenges, and review emerging strategies for validating single-cell discoveries. We also discuss future directions in applying single-cell technologies to woody perennials plants and bioenergy-relevant crops, emphasizing their potential to accelerate the discovery of cell type-specific regulatory mechanisms underlying growth, stress resilience, and biomass production.

Li, Miaomiao [ORNL] (ORCID:0000000321326168)↗

Global transcriptome profile of the developmental principles of in vitro iPSC-to-motor neuron differentiation

Abstract Background Human induced pluripotent stem cells (iPSC) have opened new avenues for regenerative medicine. Consequently, iPSC-derived motor neurons have emerged as potentially viable therapies for spinal cord injuries and neurodegenerative disorders including Amyotrophic Lateral Sclerosis. However, direct clinical application of iPSC bears in itself the risk of tumorigenesis and other unforeseeable genetic or epigenetic abnormalities. Results Employing RNA-seq technology, we identified and characterized gene regulatory networks triggered by in vitro chemical reprogramming of iPSC into cells with the molecular features of motor neurons (MNs) whose function in vivo is to innervate effector organs. We present meta-transcriptome signatures of 5 cell types: iPSCs, neural stem cells, motor neuron progenitors, early motor neurons, and mature motor neurons. In strict response to the chemical stimuli, along the MN differentiation axis we observed temporal downregulation of tumor growth factor-β signaling pathway and consistent activation of sonic hedgehog, Wnt/β-catenin, and Notch signaling. Together with gene networks defining neuronal differentiation (neurogenin 2, microtubule-associated protein 2, Pax6, and neuropilin-1), we observed steady accumulation of motor neuron-specific regulatory genes, including Islet-1 and homeobox protein HB9. Interestingly, transcriptome profiling of the differentiation process showed that Ca 2+ signaling through cAMP and LPC was downregulated during the conversion of the iPSC to neural stem cells and key regulatory gene activity of the pathway remained inhibited until later stages of motor neuron formation. Pathways shaping the neuronal development and function were well-represented in the early motor neuron cells including, neuroactive ligand-receptor interactions, axon guidance, and the cholinergic synapse formation. A notable hallmark of our in vitro motor neuron maturation in monoculture was the activation of genes encoding G-coupled muscarinic acetylcholine receptors and downregulation of the ionotropic nicotinic acetylcholine receptors expression. We observed the formation of functional neuronal networks as spontaneous oscillations in the extracellular action potentials recorded on multi-electrode array chip after 20 days of differentiation. Conclusions Detailed transcriptome profile of each developmental step from iPSC to motor neuron driven by chemical induction provides the guidelines to novel therapeutic approaches in the re-construction efforts of muscle innervation.

59 BASIC BIOLOGICAL SCIENCES↗

Widespread misregulation of inter-species hybrid transcriptomes due to sex-specific and sex-chromosome regulatory evolution

When gene regulatory networks diverge between species, their dysfunctional expression in inter-species hybrid individuals can create genetic incompatibilities that generate the developmental defects responsible for intrinsic post-zygotic reproductive isolation. Both cis- and trans-acting regulatory divergence can be hastened by directional selection through adaptation, sexual selection, and inter-sexual conflict, in addition to cryptic evolution under stabilizing selection. Dysfunctional sex-biased gene expression, in particular, may provide an important source of sexually-dimorphic genetic incompatibilities. Here, we characterize and compare male and female/hermaphrodite transcriptome profiles for sibling nematode species Caenorhabditis briggsae and C. nigoni, along with allele-specific expression in their F 1 hybrids, to deconvolve features of expression divergence and regulatory dysfunction. Despite evidence of widespread stabilizing selection on gene expression, misexpression of sex-biased genes pervades F 1 hybrids of both sexes. This finding implicates greater fragility of male genetic networks to produce dysfunctional organismal phenotypes. Spermatogenesis genes are especially prone to high divergence in both expression and coding sequences, consistent with a “faster male” model for Haldane’s rule and elevated sterility of hybrid males. Moreover, underdominant expression pervades male-biased genes compared to female-biased and sex-neutral genes and an excess of cis-trans compensatory regulatory divergence for X-linked genes underscores a “large-X effect” for hybrid male expression dysfunction. Extensive regulatory divergence in sex determination pathway genes likely contributes to demasculinization of XX hybrids. The evolution of genetic incompatibilities due to regulatory versus coding sequence divergence, however, are expected to arise in an uncorrelated fashion. This study identifies important differences between the sexes in how regulatory networks diverge to contribute to sex-biases in how genetic incompatibilities manifest during the speciation process.

59 BASIC BIOLOGICAL SCIENCES↗

Towards the Reconstruction of Integrated Genome-Scale Models of Metabolism and Gene Expression

The reconstruction of integrated genome-scale models of metabolism and gene expression has been a challenge for a while now. In fact, various methods that allow integrating reconstructions of Transcriptional Regulatory Networks, gene expression data or both into Genome-Scale Metabolic Models have been proposed. Several of these methods are surveyed in this article, which allowed identifying their strengths and weaknesses concerning the reconstruction of integrated models for multiple prokaryotic organisms. Additionally, the main resources of regulatory information were also surveyed, as the existence of novel sources of regulatory information and gene expression data may contribute for the improvement of methodologies referred herein.

Databases of regulatory information↗

Disrupting autorepression circuitry generates “open-loop lethality” to yield escape-resistant antiviral agents

Across biological scales, gene-regulatory networks employ autorepression (negative feedback) to maintain homeostasis and minimize failure from aberrant expression. In this study, we present a proof of concept that disrupting transcriptional negative feedback dysregulates viral gene expression to therapeutically inhibit replication and confers a high evolutionary barrier to resistance. We find that nucleic-acid decoys mimicking cis-regulatory sites act as “feedback disruptors,” break homeostasis, and increase viral transcription factors to cytotoxic levels (termed “open-loop lethality”). Feedback disruptors against herpesviruses reduced viral replication >2-logs without activating innate immunity, showed sub-nM IC 50 , synergized with standard-of-care antivirals, and inhibited virus replication in mice. In contrast to approved antivirals where resistance rapidly emerged, no feedback-disruptor escape mutants evolved in long-term cultures. For SARS-CoV-2, disruption of a putative feedback circuit also generated open-loop lethality, reducing viral titers by >1-log. These results demonstrate that generating open-loop lethality, via negative-feedback disruption, may yield a class of antimicrobials with a high genetic barrier to resistance.

59 BASIC BIOLOGICAL SCIENCES↗

Exploring Camelina sativa lipid metabolism regulation by combining gene co-expression and DNA affinity purification analyses

Camelina (Camelina sativa) is an annual oilseed plant that is gaining momentum as a biofuel cover crop. Understanding gene regulatory networks is essential to deciphering plant metabolic pathways, including lipid metabolism. Furthermore, we take advantage of a growing collection of gene expression datasets to predict transcription factors (TFs) associated with the control of Camelina lipid metabolism. We identified approximately 350 TFs highly co-expressed with lipid-related genes (LRGs). These TFs are highly represented in the MYB, AP2/ERF, bZIP, and bHLH families, including a significant number of homologs of well-known Arabidopsis lipid and seed developmental regulators. After prioritizing the top 22 TFs for further validation, we identified DNA-binding sites and predicted target genes for 16 out of the 22 TFs tested using DNA affinity purification followed by sequencing (DAP-seq). Enrichment analyses of targets supported the co-expression prediction for most TF candidates, and the comparison to Arabidopsis revealed some common themes, but also aspects unique to Camelina. Within the top potential lipid regulators, we identified CsaMYB1, CsaABI3AVP1-2, CsaHB1, CsaNAC2, CsaMYB3, and CsaNAC1 as likely involved in the control of seed fatty acid elongation and CsaABI3AVP1-2 and CsabZIP1 as potential regulators of the synthesis and degradation of triacylglycerols (TAGs), respectively. Altogether, the integration of co-expression data and DNA-binding assays permitted us to generate a high-confidence and short list of Camelina TFs involved in the control of lipid metabolism during seed development.

59 BASIC BIOLOGICAL SCIENCES↗

Bacterium-enabled transient gene activation by artificial transcription factors for resolving gene regulation in maize

Abstract Understanding gene regulatory networks is essential to elucidate developmental processes and environmental responses. Here, we studied regulation of a maize (Zea mays) transcription factor gene using designer transcription activator-like effectors (dTALes), which are synthetic Type III TALes of the bacterial genus Xanthomonas and serve as inducers of disease susceptibility gene transcription in host cells. The maize pathogen Xanthomonas vasicola pv. vasculorum was used to introduce 2 independent dTALes into maize cells to induced expression of the gene glossy3 (gl3), which encodes a MYB transcription factor involved in biosynthesis of cuticular wax. RNA-seq analysis of leaf samples identified, in addition to gl3, 146 genes altered in expression by the 2 dTALes. Nine of the 10 genes known to be involved in cuticular wax biosynthesis were upregulated by at least 1 of the 2 dTALes. A gene previously unknown to be associated with gl3, Zm00001d017418, which encodes aldehyde dehydrogenase, was also expressed in a dTALe-dependent manner. A chemically induced mutant and a CRISPR-Cas9 mutant of Zm00001d017418 both exhibited glossy leaf phenotypes, indicating that Zm00001d017418 is involved in biosynthesis of cuticular waxes. Bacterial protein delivery of dTALes proved to be a straightforward and practical approach for the analysis and discovery of pathway-specific genes in maize.

59 BASIC BIOLOGICAL SCIENCES↗

Sorghum bicolor BTx623 Nitrogen Grown Conditions Set2 Gene Expression Profiling

Dhurrin, a cyanogenic glucoside, plays an important role in Sorghum bicolor physiology and defense. The concentration of dhurrin in sorghum is influenced by both nitrogen status and stage of plant organ development. While nitrogen resupply activates the expression of genes for dhurrin biosynthesis, the molecular mechanisms underlying this regulation remain unclear. In this study, we investigated the transcriptional response of sorghum to nitrogen resupply following growth under nitrogen-limiting conditions. Using a time-course design, we measured hydrogen cyanide potential (HCNp), growth, and nitrate content at 0-, 2-, 6-, 12-, 24-, 36-, 48-, and 60-h after resupply and collected tissue for RNAseq analysis in parallel for analysis of gene expression and construction of gene regulatory networks (GRNs). HCNp (mg g−1 DW) increased significantly in leaf and stem tissues following nitrogen resupply, with increases in the leaf partially driven by continued declines in controls under ongoing nitrogen stress. Expression of the dhurrin pathway genes was upregulated in leaves from 24 h after nitrogen resupply, with diel expression patterns observable over the remaining time points. No upregulation was observed in roots or stems, suggesting that developmental context overrides environmental cues. GRN analysis identified candidate transcription factors regulating dhurrin biosynthesis genes, including members of the MYB, bZIP, and GARP-type transcription factor families. Some of these candidate transcription factors may be involved in relieving senescence-associated suppression of dhurrin biosynthesis and link nitrogen signaling to pathway activation. These findings provide new insight into the nitrogen-responsive regulation of dhurrin in sorghum, highlighting candidate regulators for future functional characterization.

cyanogenic glucoside↗

Sorghum bicolor BTx623 Nitrogen Grown Conditions Gene Expression Profiling

Dhurrin, a cyanogenic glucoside, plays an important role in Sorghum bicolor physiology and defense. The concentration of dhurrin in sorghum is influenced by both nitrogen status and stage of plant organ development. While nitrogen resupply activates the expression of genes for dhurrin biosynthesis, the molecular mechanisms underlying this regulation remain unclear. In this study, we investigated the transcriptional response of sorghum to nitrogen resupply following growth under nitrogen-limiting conditions. Using a time-course design, we measured hydrogen cyanide potential (HCNp), growth, and nitrate content at 0-, 2-, 6-, 12-, 24-, 36-, 48-, and 60-h after resupply and collected tissue for RNAseq analysis in parallel for analysis of gene expression and construction of gene regulatory networks (GRNs). HCNp (mg g−1 DW) increased significantly in leaf and stem tissues following nitrogen resupply, with increases in the leaf partially driven by continued declines in controls under ongoing nitrogen stress. Expression of the dhurrin pathway genes was upregulated in leaves from 24 h after nitrogen resupply, with diel expression patterns observable over the remaining time points. No upregulation was observed in roots or stems, suggesting that developmental context overrides environmental cues. GRN analysis identified candidate transcription factors regulating dhurrin biosynthesis genes, including members of the MYB, bZIP, and GARP-type transcription factor families. Some of these candidate transcription factors may be involved in relieving senescence-associated suppression of dhurrin biosynthesis and link nitrogen signaling to pathway activation. These findings provide new insight into the nitrogen-responsive regulation of dhurrin in sorghum, highlighting candidate regulators for future functional characterization.

cyanogenic glucoside↗

Nitrogen Status Rewires Transcriptional Regulation of Dhurrin, a Dual‐Purpose Defense Metabolite in Sorghum bicolor

Dhurrin, a cyanogenic glucoside, plays an important role in Sorghum bicolor physiology and defense. The concentration of dhurrin in sorghum is influenced by both nitrogen status and stage of plant organ development. While nitrogen resupply activates the expression of genes for dhurrin biosynthesis, the molecular mechanisms underlying this regulation remain unclear. In this study, we investigated the transcriptional response of sorghum to nitrogen resupply following growth under nitrogen-limiting conditions. Using a time-course design, we measured hydrogen cyanide potential (HCNp), growth, and nitrate content at 0-, 2-, 6-, 12-, 24-, 36-, 48-, and 60-h after resupply and collected tissue for RNAseq analysis in parallel for analysis of gene expression and construction of gene regulatory networks (GRNs). HCNp (mg g −1 DW) increased significantly in leaf and stem tissues following nitrogen resupply, with increases in the leaf partially driven by continued declines in controls under ongoing nitrogen stress. Expression of the dhurrin pathway genes was upregulated in leaves from 24 h after nitrogen resupply, with diel expression patterns observable over the remaining time points. No upregulation was observed in roots or stems, suggesting that developmental context overrides environmental cues. GRN analysis identified candidate transcription factors regulating dhurrin biosynthesis genes, including members of the MYB, bZIP, and GARP-type transcription factor families. Some of these candidate transcription factors may be involved in relieving senescence-associated suppression of dhurrin biosynthesis and link nitrogen signaling to pathway activation. These findings provide new insight into the nitrogen-responsive regulation of dhurrin in sorghum, highlighting candidate regulators for future functional characterization.

S. bicolor↗

Disentangling direct from indirect relationships in association networks

Networks are vital tools for understanding and modeling interactions in complex systems in science and engineering, and direct and indirect interactions are pervasive in all types of networks. However, quantitatively disentangling direct and indirect relationships in networks remains a formidable task. Here, we present a framework, called iDIRECT (Inference of Direct and Indirect Relationships with Effective Copula-based Transitivity), for quantitatively inferring direct dependencies in association networks. Using copula-based transitivity, iDIRECT eliminates/ameliorates several challenging mathematical problems, including ill-conditioning, self-looping, and interaction strength overflow. With simulation data as benchmark examples, iDIRECT showed high prediction accuracies. Application of iDIRECT to reconstruct gene regulatory networks in Escherichia coli also revealed considerably higher prediction power than the best-performing approaches in the DREAM5 (Dialogue on Reverse Engineering Assessment and Methods project, #5) Network Inference Challenge. In addition, applying iDIRECT to highly diverse grassland soil microbial communities in response to climate warming showed that the iDIRECT-processed networks were significantly different from the original networks, with considerably fewer nodes, links, and connectivity, but higher relative modularity. Further analysis revealed that the iDIRECT-processed network was more complex under warming than the control and more robust to both random and target species removal ( P < 0.001). As a general approach, iDIRECT has great advantages for network inference, and it should be widely applicable to infer direct relationships in association networks across diverse disciplines in science and engineering.

59 BASIC BIOLOGICAL SCIENCES↗

Parallel measurement of transcriptomes and proteomes from same single cells using nanodroplet splitting

Single-cell multiomics provides comprehensive insights into gene regulatory networks, cellular diversity, and temporal dynamics. While tools for co-profiling single-cell genomes, transcriptomes, and epigenomes are available, accessing proteomes in parallel is more challenging. We developed nanoSPLITS (nanodroplet SPlitting for Linked-multimodal Investigations of Trace Samples), an integrated platform that enables global profiling of the transcriptome and proteome from same single cells using RNA sequencing and mass spectrometry-based proteomics, respectively. nanoSPLITS can precisely quantify over 5000 genes, 2000 proteins, and 140 phosphopeptides per single cell and identify candidate cell markers from these modalities. By exploring Cdk1-mediated cell cycle arrest, we demonstrate how nanoSPLITS single-cell multiomics can provide comprehensive cellular characterization with insights into covarying protein/gene clusters, unique phosphorylation events, and mitotic pathways.

59 BASIC BIOLOGICAL SCIENCES↗

Learning genetic perturbation effects with variational causal inference

Advances in sequencing technologies have enhanced the understanding of gene regulation in cells. In particular, Perturb-seq has enabled high-resolution profiling of the transcriptomic response to genetic perturbations at the single-cell level. This understanding has implications in functional genomics and potentially for identifying therapeutic targets. Various computational models have been developed to predict perturbational effects. While deep learning models excel at interpolating observed perturbational data, they tend to overfit in the lack of enough data and may not generalize well to unseen perturbations. In contrast, mechanistic models, such as linear causal models based on gene regulatory networks, hold greater potential for extrapolation, as they encapsulate regulatory information that can predict responses to unseen perturbations. However, their application has been limited to small studies due to overly simplistic assumptions, making them less effective in handling noisy, large-scale single-cell data. We propose a hybrid approach that combines a mechanistic causal model with variational deep learning, termed Single Cell Causal Variational Autoencoder (SCCVAE). The mechanistic model employs a learned regulatory network to represent perturbational changes as shift interventions that propagate through the learned network. SCCVAE integrates this mechanistic causal model into a variational autoencoder, generating rich, comprehensive transcriptomic responses. Our results indicate that SCCVAE exhibits superior performance over current state-of-the-art baselines for extrapolating to predict unseen perturbational responses. Additionally, for the observed perturbations, the latent space learned by SCCVAE allows for the identification of functional perturbation modules and simulation of single-gene knockdown experiments of varying penetrance, presenting a robust tool for interpreting and interpolating perturbational responses at the single-cell level.

59 BASIC BIOLOGICAL SCIENCES↗

SOX2-driven enhancer landscape defines the transcriptional architecture of retinogenesis

Retinal neurogenesis is mediated by the coordinated activities of a complex gene regulatory network (GRN) of transcription factors (TFs) in multipotent retinal progenitor cells (RPCs). How this GRN mechanistically guides neural competence remains poorly understood. In this study, we present integrated transcriptional, genetic and genomic analyses to uncover the regulatory mechanisms of SOX2, a key factor in establishing neural identity in RPCs. We show that SOX2 is preferentially enriched in the RPC-specific enhancer landscape associated with essential regulators of retinogenesis. Disruption of SOX2 expression impairs retinogenesis, marked by a selective loss of enhancer activity near genes essential for RPC proliferation and lineage specification. We identified the RPC transcription factor VSX2 as a binding partner for SOX2 and, together, SOX2 and VSX2 co-target a core, retina-specific chromatin repertoire characterized by enhanced TF binding and robust chromatin accessibility. This cooperative binding establishes a shared SOX2-VSX2 transcriptional code that promotes the expression of crucial regulators of neurogenesis while repressing the acquisition of alternative lineage cell fate. Our data illuminate fundamental biological insights on how transcription factors act in concert to drive chromatin-based genetic programs underlying retinal neural identity.

Chromatin↗

Canalization of Phenotypes—When the Transcriptome is Constantly but Weakly Perturbed

Abstract Recent studies have increasingly pointed to microRNAs (miRNAs) as the agent of gene regulatory network (GRN) stabilization as well as developmental canalization against constant but small environmental perturbations. To analyze mild perturbations, we construct a Dicer-1 knockdown line (dcr-1 KD) in Drosophila that modestly reduces all miRNAs by, on average, ∼20%. The defining characteristic of stabilizers is that, when their capacity is compromised, GRNs do not change their short-term behaviors. Indeed, even with such broad reductions across all miRNAs, the changes in the transcriptome are very modest during development in stable environment. By comparison, broad knockdowns of other regulatory genes (esp. transcription factors) by the same method should lead to drastic changes in the GRNs. The consequence of destabilization may thus be in long-term development as postulated by the theory of canalization. Flies with modest miRNA reductions may gradually deviate from the developmental norm, resulting in late-stage failures such as shortened longevity. In the optimal culture condition, the survival to adulthood is indeed normal in the dcr-1 KD line but, importantly, adult longevity is reduced by ∼90%. When flies are stressed by high temperature, dcr-1 KD induces lethality earlier in late pupation and, as the perturbations are shifted earlier, the affected stages are shifted correspondingly. Hence, in late stages of development with deviations piling up, GRN would be increasingly in need of stabilization. In conclusion, miRNAs appear to be a solution to weak but constant environmental perturbations.

Lu, Guang-An↗

Atlas of Transcription Factor Binding Sites from ENCODE DNase Hypersensitivity Data across 27 Tissue Types

Characterizing the tissue-specific binding sites of transcription factors (TFs) is essential to reconstruct gene regulatory networks and predict functions for non-coding genetic variation. DNase-seq footprinting enables the prediction of genome-wide binding sites for hundreds of TFs simultaneously. Despite the public availability of high-quality DNase-seq data from hundreds of samples, a comprehensive, up-to-date resource for the locations of genomic footprints is lacking. Here, we develop a scalable footprinting workflow using two state-of-the-art algorithms: Wellington and HINT. We apply our workflow to detect footprints in 192 ENCODE DNase-seq experiments and predict the genomic occupancy of 1,515 human TFs in 27 human tissues. We validate that these footprints overlap true-positive TF binding sites from ChIP-seq. We demonstrate that the locations, depth, and tissue specificity of footprints predict effects of genetic variants on gene expression and capture a substantial proportion of genetic risk for complex traits.

59 BASIC BIOLOGICAL SCIENCES↗