RESUMO
Circuit building blocks of gene regulatory networks (GRN) have been identified through the fibration symmetries of the underlying biological graph. Here, we analyse analytically six of these circuits that occur as functional and synchronous building blocks in these networks. Of these, the lock-on, toggle switch, Smolen oscillator, feed-forward fibre and Fibonacci fibre circuits occur in living organisms, notably Escherichia coli; the sixth, the repressilator, is a synthetic GRN. We consider synchronous steady states determined by a fibration symmetry (or balanced colouring) and determine analytic conditions for local bifurcation from such states, which can in principle be either steady-state or Hopf bifurcations. We identify conditions that characterize the first bifurcation, the only one that can be stable near the bifurcation point. We model the state of each gene in terms of two variables: mRNA and protein concentration. We consider all possible 'admissible' models-those compatible with the network structure-and then specialize these general results to simple models based on Hill functions and linear degradation. The results systematically classify using graph symmetries the complexity and dynamics of these circuits, which are relevant to understand the functionality of natural and synthetic cells.
Assuntos
Escherichia coli , Redes Reguladoras de Genes , Modelos Genéticos , Escherichia coli/genética , Escherichia coli/metabolismoRESUMO
Dave Garbers' work significantly contributed to our understanding of sperm's regulated motility, capacitation, and the acrosome reaction. These key sperm functions involve complex multistep signaling pathways engaging numerous finely orchestrated elements. Despite significant progress, many parameters and interactions among these elements remain elusive. Mathematical modeling emerges as a potent tool to study sperm physiology, providing a framework to integrate experimental results and capture functional dynamics considering biochemical, biophysical, and cellular elements. Depending on research objectives, different modeling strategies, broadly categorized into continuous and discrete approaches, reveal valuable insights into cell function. These models allow the exploration of hypotheses regarding molecules, conditions, and pathways, whenever they become challenging to evaluate experimentally. This review presents an overview of current theoretical and experimental efforts to understand sperm motility regulation, capacitation, and the acrosome reaction. We discuss the strengths and weaknesses of different modeling strategies and highlight key findings and unresolved questions. Notable discoveries include the importance of specific ion channels, the role of intracellular molecular heterogeneity in capacitation and the acrosome reaction, and the impact of pH changes on acrosomal exocytosis. Ultimately, this review underscores the crucial importance of mathematical frameworks in advancing our understanding of sperm physiology and guiding future experimental investigations.
Assuntos
Reação Acrossômica , Transdução de Sinais , Capacitação Espermática , Motilidade dos Espermatozoides , Espermatozoides , Masculino , Espermatozoides/metabolismo , Espermatozoides/fisiologia , Humanos , Reação Acrossômica/fisiologia , Capacitação Espermática/fisiologia , Transdução de Sinais/fisiologia , Animais , Motilidade dos Espermatozoides/fisiologia , Modelos Biológicos , Modelos TeóricosRESUMO
Characterization of shared patterns of RNA expression between genes across conditions has led to the discovery of regulatory networks and novel biological functions. However, it is unclear if such coordination extends to translation, a critical step in gene expression. Here, we uniformly analyzed 3,819 ribosome profiling datasets from 117 human and 94 mouse tissues and cell lines. We introduce the concept of Translation Efficiency Covariation (TEC), identifying coordinated translation patterns across cell types. We nominate potential mechanisms driving shared patterns of translation regulation. TEC is conserved across human and mouse cells and helps uncover gene functions. Moreover, our observations indicate that proteins that physically interact are highly enriched for positive covariation at both translational and transcriptional levels. Our findings establish translational covariation as a conserved organizing principle of mammalian transcriptomes.
RESUMO
Single-cell CRISPR screens link genetic perturbations to transcriptional states, but high-throughput methods connecting these induced changes to their regulatory foundations are limited. Here we introduce Multiome Perturb-seq, extending single-cell CRISPR screens to simultaneously measure perturbation-induced changes in gene expression and chromatin accessibility. We apply Multiome Perturb-seq in a CRISPRi screen of 13 chromatin remodelers in human RPE-1 cells, achieving efficient assignment of sgRNA identities to single nuclei via an improved method for capturing barcode transcripts from nuclear RNA. We organize expression and accessibility measurements into coherent programs describing the integrated effects of perturbations on cell state, finding that ARID1A and SUZ12 knockdowns induce programs enriched for developmental features. Pseudotime analysis of perturbations connects accessibility changes to changes in gene expression, highlighting the value of multimodal profiling. Overall, our method provides a scalable and simply implemented system to dissect the regulatory logic underpinning cell state.
RESUMO
BACKGROUND: Lung adenocarcinoma (LUAD) has been observed to have significant sex differences in incidence, prognosis, and response to therapy. However, the molecular mechanisms responsible for these disparities have not been investigated extensively. METHODS: Sample-specific gene regulatory network methods were used to analyze RNA sequencing data from non-cancerous human lung samples from The Genotype Tissue Expression Project (GTEx) and lung adenocarcinoma primary tumor samples from The Cancer Genome Atlas (TCGA); results were validated on independent data. RESULTS: We found that genes associated with key biological pathways including cell proliferation, immune response and drug metabolism are differentially regulated between males and females in both healthy lung tissue and tumor, and that these regulatory differences are further perturbed by tobacco smoking. We also discovered significant sex bias in transcription factor targeting patterns of clinically actionable oncogenes and tumor suppressor genes, including AKT2 and KRAS. Using differentially regulated genes between healthy and tumor samples in conjunction with a drug repurposing tool, we identified several small-molecule drugs that might have sex-biased efficacy as cancer therapeutics and further validated this observation using an independent cell line database. CONCLUSIONS: These findings underscore the importance of including sex as a biological variable and considering gene regulatory processes in developing strategies for disease prevention and management.
Lung adenocarcinoma (LUAD) is a disease that affects males and females differently. Biological sex not only influences chances of developing the disease, but also how the disease progresses and how effective various therapies may be. We analyzed sex-specific gene regulatory networks consisting of transcription factors and the genes they regulate in both healthy lung tissue and in LUAD and identified sex-biased differences. We found that genes associated with cell proliferation, immune response, and drug metabolism are differentially targeted by transcription factors between males and females. We also found that several genes that are drug targets in LUAD, are also regulated differently between males and females. Importantly, these differences are also influenced by an individual's smoking history. Extending our analysis using a drug repurposing tool, we found candidate drugs with evidence that they might work better for one sex or the other. These results demonstrate that considering the differences in gene regulation between males and females will be essential if we are to develop precision medicine strategies for preventing and treating LUAD.
Assuntos
Adenocarcinoma de Pulmão , Redes Reguladoras de Genes , Adenocarcinoma de Pulmão/diagnóstico , Adenocarcinoma de Pulmão/genética , Adenocarcinoma de Pulmão/terapia , Fatores Sexuais , Regulação Neoplásica da Expressão Gênica/genética , Pulmão/metabolismo , Fumar Tabaco/efeitos adversos , Prognóstico , Imunoterapia , Terapia de Alvo Molecular , Linhagem Celular Tumoral , Humanos , Masculino , Feminino , Descoberta de DrogasRESUMO
This paper focuses on inferring a general class of hidden Markov models (HMMs) using data acquired from experts. Expert-acquired data contain decisions/actions made by humans/users for various objectives, such as navigation data reflecting drivers' behavior, cybersecurity data carrying defender decisions, and biological data containing the biologist's actions (e.g., interventions and experiments). Conventional inference methods rely on temporal changes in data without accounting for expert knowledge. This paper incorporates expert knowledge into the inference of HMMs by modeling expert behavior as an imperfect reinforcement learning agent. The proposed method optimally quantifies experts' perceptions about the system model, which, alongside the temporal changes in data, contributes to the inference process. The proposed inference method is derived through a combination of dynamic programming and optimal recursive Bayesian estimation. The applicability of this method is demonstrated to a wide range of inference criteria, such as maximum likelihood and maximum a posteriori. The performance of the proposed method is investigated through a comprehensive numerical experiment using a benchmark problem and biological networks.
RESUMO
Compared to men, women often develop COPD at an earlier age with worse respiratory symptoms despite lower smoking exposure. However, most preventive, and therapeutic strategies ignore biological sex differences in COPD. Our goal was to better understand sex-specific gene regulatory processes in lung tissue and the molecular basis for sex differences in COPD onset and severity. We analyzed lung tissue gene expression and DNA methylation data from 747 individuals in the Lung Tissue Research Consortium (LTRC), and 85 individuals in an independent dataset. We identified sex differences in COPD-associated gene regulation using gene regulatory networks. We used linear regression to test for sex-biased associations of methylation with lung function, emphysema, smoking, and age. Analyzing gene regulatory networks in the control group, we identified that genes involved in the extracellular matrix (ECM) have higher transcriptional factor targeting in females than in males. However, this pattern is reversed in COPD, with males showing stronger regulatory targeting of ECM-related genes than females. Smoking exposure, age, lung function, and emphysema were all associated with sex-specific differential methylation of ECM-related genes. We identified sex-based gene regulatory patterns of ECM-related genes associated with lung function and emphysema. Multiple factors including epigenetics, smoking, aging, and cell heterogeneity influence sex-specific gene regulation in COPD. Our findings underscore the importance of considering sex as a key factor in disease susceptibility and severity.
RESUMO
BACKGROUND: The inflammation phenotypes are often closely related to oxidative stress and autophagy pathway activation, which could be developed as a treatment target. AIMS: The aim of this study was to explore the underlying mechanism of inflammation in chronic obstructive pulmonary disease (COPD). METHODS: The lung tissue single-cell RNA-seq (scRNA-seq) dataset of GSE171541 was downloaded from the Gene Expression Omnibus (GEO) database. The marker genes were obtained from the CellMarker database. "Seurat" and "harmony" R packages were used for single-cell profiling analysis. Then, the "AUCell" R package was employed to calculate the reactive oxygen species (ROS) and autophagy pathway scores. Gene regulation network analysis was performed by applying the "SCENIC" package, followed by conducting correlation analysis with Spearman's rank correlation method. The cigarettes were used to develop a traumatic model in mice, and the expression of relevant genes was measured by qRT-PCR. RESULTS: The scRNA-seq analysis classified 12 cell subgroups in which the contractility of myofibroblasts was closely associated with the progression of COPD. Further analysis showed that ROS and autophagy pathways were significantly activated in myofibroblasts and that the nuclear factor erythroid 2-related factor 2 (NRF2) and its mediated oxidative stress pathway were inhibited in myofibroblasts. In addition, the downregulated NRF2 gene was negatively correlated with the expression of autophagy and ROS activation. In the traumatic mice model, NRF2 was downregulated in COPD mice but further elevated in the COPD+NRF2 mice group. Interestingly, the mRNA levels of Kelchlike ECH-associated protein 1 (Keap1), NADPH oxidase (NOX), and Cathepsin B (CTSB) were upregulated in COPD group in comparison to the control group but they were downregulated by NRF2. These results suggested that low-expressed NFR2 promoted autophagy and ROS pathway activation in myofibroblasts for COPD progression. CONCLUSION: We identified a cell myofibroblast cluster closely associated with COPD progression using the scRNA-seq analysis. The downregulated NFR2, as a key risk factor, mediated myofibroblast death by activating the oxidative stress and autophagy pathway for COPD progression.
RESUMO
Typical 'omic analyses reduce complex biological systems to simple lists of supposedly independent variables, failing to account for changes in the wider transcriptional landscape. In this commentary, we discuss the utility of network approaches for incorporating this wider context into the study of physiological phenomena. We highlight opportunities to build on traditional network tools by utilising cutting-edge techniques to account for higher order interactions (i.e. beyond pairwise associations) within datasets, allowing for more accurate models of complex 'omic systems. Finally, we show examples of previous works utilising network approaches to gain additional insight into their organisms of interest. As 'omics grow in both their popularity and breadth of application, so does the requirement for flexible analytical tools capable of interpreting and synthesising complex datasets.
RESUMO
The genetic architecture of Parkinson's disease (PD) is complex and multiple brain cell subtypes are involved in the neuropathological progression of the disease. Here we aimed to advance our understanding of PD genetic complexity at a cell subtype precision level. Using parallel single-nucleus (sn)RNA-seq and snATAC-seq analyses we simultaneously profiled the transcriptomic and chromatin accessibility landscapes in temporal cortex tissues from 12 PD compared to 12 control subjects at a granular single cell resolution. An integrative bioinformatic pipeline was developed and applied for the analyses of these snMulti-omics datasets. The results identified a subpopulation of cortical glutamatergic excitatory neurons with remarkably altered gene expression in PD, including differentially-expressed genes within PD risk loci identified in genome-wide association studies (GWAS). This was the only neuronal subtype showing significant and robust overexpression of SNCA. Further characterization of this neuronal-subpopulation showed upregulation of specific pathways related to axon guidance, neurite outgrowth and post-synaptic structure, and downregulated pathways involved in presynaptic organization and calcium response. Additionally, we characterized the roles of three molecular mechanisms in governing PD-associated cell subtype-specific dysregulation of gene expression: (1) changes in cis-regulatory element accessibility to transcriptional machinery; (2) changes in the abundance of master transcriptional regulators, including YY1, SP3, and KLF16; (3) candidate regulatory variants in high linkage disequilibrium with PD-GWAS genomic variants impacting transcription factor binding affinities. To our knowledge, this study is the first and the most comprehensive interrogation of the multi-omics landscape of PD at a cell-subtype resolution. Our findings provide new insights into a precise glutamatergic neuronal cell subtype, causal genes, and non-coding regulatory variants underlying the neuropathological progression of PD, paving the way for the development of cell- and gene-targeted therapeutics to halt disease progression as well as genetic biomarkers for early preclinical diagnosis.
Assuntos
Redes Reguladoras de Genes , Neurônios , Doença de Parkinson , Humanos , Doença de Parkinson/genética , Doença de Parkinson/metabolismo , Doença de Parkinson/patologia , Neurônios/metabolismo , Neurônios/patologia , Masculino , Feminino , alfa-Sinucleína/genética , alfa-Sinucleína/metabolismo , Idoso , Fator de Transcrição YY1/genética , Fator de Transcrição YY1/metabolismo , Estudo de Associação Genômica Ampla , Transcriptoma , Análise de Célula Única , Lobo Temporal/metabolismo , Lobo Temporal/patologia , Pessoa de Meia-Idade , Regulação da Expressão Gênica/genética , MultiômicaRESUMO
Inferring gene regulatory networks (GRNs) allows us to obtain a deeper understanding of cellular function and disease pathogenesis. Recent advances in single-cell RNA sequencing (scRNA-seq) technology have improved the accuracy of GRN inference. However, many methods for inferring individual GRNs from scRNA-seq data are limited because they overlook intercellular heterogeneity and similarities between different cell subpopulations, which are often present in the data. Here, we propose a deep learning-based framework, DeepGRNCS, for jointly inferring GRNs across cell subpopulations. We follow the commonly accepted hypothesis that the expression of a target gene can be predicted based on the expression of transcription factors (TFs) due to underlying regulatory relationships. We initially processed scRNA-seq data by discretizing data scattering using the equal-width method. Then, we trained deep learning models to predict target gene expression from TFs. By individually removing each TF from the expression matrix, we used pre-trained deep model predictions to infer regulatory relationships between TFs and genes, thereby constructing the GRN. Our method outperforms existing GRN inference methods for various simulated and real scRNA-seq datasets. Finally, we applied DeepGRNCS to non-small cell lung cancer scRNA-seq data to identify key genes in each cell subpopulation and analyzed their biological relevance. In conclusion, DeepGRNCS effectively predicts cell subpopulation-specific GRNs. The source code is available at https://github.com/Nastume777/DeepGRNCS.
Assuntos
Aprendizado Profundo , Redes Reguladoras de Genes , Análise de Célula Única , Humanos , Análise de Célula Única/métodos , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismo , Biologia Computacional/métodos , Análise de Sequência de RNA/métodos , RNA-Seq/métodosRESUMO
Abscisic acid (ABA) significantly regulates plant growth and development, promoting tuberous root formation in various plants. However, the molecular mechanisms of ABA in the tuberous root development of Pseudostellaria heterophylla are not yet fully understood. This study utilized Illumina sequencing and de novo assembly strategies to obtain a reference transcriptome associated with ABA treatment. Subsequently, integrated transcriptomic and proteomic analyses were used to determine gene expression profiles in P. heterophylla tuberous roots. ABA treatment significantly increases the diameter and shortens the length of tuberous roots. Clustering analysis identified 2,256 differentially expressed genes and 679 differentially abundant proteins regulated by ABA. Gene co-expression and protein interaction networks revealed ABA positively induced 30 vital regulators. Furthermore, we identified and assigned putative functions to transcription factors (PhMYB10, PhbZIP2, PhbZIP, PhSBP) that mediate ABA signaling involved in the regulation of tuberous root development, including those related to cell wall metabolism, cell division, starch synthesis, hormone metabolism. Our findings provide valuable insights into the complex signaling networks of tuberous root development modulated by ABA. It provided potential targets for genetic manipulation to improve the yield and quality of P. heterophylla, which could significantly impact its cultivation and medicinal value.
RESUMO
The Gram-negative pathogen Acinetobacter baumannii is considered an "urgent threat" to human health due to its propensity to become antibiotic resistant. Understanding the distinct regulatory paradigms used by A. baumannii to mitigate cellular stresses may uncover new therapeutic targets. Many γ-proteobacteria use the extracytoplasmic function (ECF) σ factor, RpoE, to invoke envelope homeostasis networks in response to stress. Acinetobacter species contain the poorly characterized ECF "SigAb;" however, it is unclear if SigAb has the same physiological role as RpoE. Here, we show that SigAb is a metal stress-responsive ECF that appears unique to Acinetobacter species and distinct from RpoE. We combine promoter mutagenesis, motif scanning, and ChIP-seq to define the direct SigAb regulon, which consists of sigAb itself, the stringent response mediator, relA, and the uncharacterized small RNA, "sabS." However, RNA-seq of strains overexpressing SigAb revealed a large, indirect regulon containing hundreds of genes. Metal resistance genes are key elements of the indirect regulon, as CRISPRi knockdown of sigAb or sabS resulted in increased copper sensitivity and excess copper induced SigAb-dependent transcription. Further, we found that two uncharacterized genes in the sigAb operon, "aabA" and "aabB", have anti-SigAb activity. Finally, employing a targeted Tn-seq approach that uses CRISPR-associated transposons, we show that sigAb, aabA, and aabB are important for fitness even during optimal growth conditions. Our work reveals new physiological roles for SigAb and SabS, provides a novel approach for assessing gene fitness, and highlights the distinct regulatory architecture of A. baumannii.
RESUMO
Understanding the dynamics of gene regulatory networks (GRNs) across diverse cell types poses a challenge yet holds immense value in unraveling the molecular mechanisms governing cellular processes. Current computational methods, which rely solely on expression changes from bulk RNA-seq and/or scRNA-seq data, often result in high rates of false positives and low precision. Here, we introduce an advanced computational tool, DeepIMAGER, for inferring cell-specific GRNs through deep learning and data integration. DeepIMAGER employs a supervised approach that transforms the co-expression patterns of gene pairs into image-like representations and leverages transcription factor (TF) binding information for model training. It is trained using comprehensive datasets that encompass scRNA-seq profiles and ChIP-seq data, capturing TF-gene pair information across various cell types. Comprehensive validations on six cell lines show DeepIMAGER exhibits superior performance in ten popular GRN inference tools and has remarkable robustness against dropout-zero events. DeepIMAGER was applied to scRNA-seq datasets of multiple myeloma (MM) and detected potential GRNs for TFs of RORC, MITF, and FOXD2 in MM dendritic cells. This technical innovation, combined with its capability to accurately decode GRNs from scRNA-seq, establishes DeepIMAGER as a valuable tool for unraveling complex regulatory networks in various cell types.
Assuntos
Redes Reguladoras de Genes , RNA-Seq , Humanos , Biologia Computacional/métodos , Aprendizado Profundo , Mieloma Múltiplo/genética , RNA-Seq/métodos , Análise da Expressão Gênica de Célula Única , Software , Fatores de Transcrição/metabolismo , Fatores de Transcrição/genéticaRESUMO
Transcription factors (TFs) are crucial in modulating gene expression and sculpting cellular and organismal phenotypes. The identification of TF-target gene interactions is pivotal for comprehending molecular pathways and disease etiologies but has been hindered by the demanding nature of traditional experimental approaches. This paper introduces a novel web application and package utilizing the R program, which predicts TF-target gene relationships and vice versa. Our application integrates the predictive power of various bioinformatic tools, leveraging their combined strengths to provide robust predictions. It merges databases for enhanced precision, incorporates gene expression correlation for accuracy, and employs pan-tissue correlation analysis for context-specific insights. The application also enables the integration of user data with established resources to analyze TF-target gene networks. Despite its current limitation to human data, it provides a platform to explore gene regulatory mechanisms comprehensively. This integrated, systematic approach offers researchers an invaluable tool for dissecting the complexities of gene regulation, with the potential for future expansions to include a broader range of species.
Assuntos
Biologia Computacional , Redes Reguladoras de Genes , Software , Fatores de Transcrição , Humanos , Fatores de Transcrição/metabolismo , Fatores de Transcrição/genética , Biologia Computacional/métodos , Regulação da Expressão Gênica , Bases de Dados GenéticasRESUMO
The development of highly specialized blood cells from hematopoietic stem cells (HSCs) in the bone marrow (BM) is dependent upon a stringently orchestrated network of stage- and lineage-restricted transcription factors (TFs). Thus, the same stem cell can give rise to various types of differentiated blood cells. One of the key regulators of B-lymphocyte development is early B-cell factor 1 (EBF1). This TF belongs to a small, but evolutionary conserved, family of proteins that harbor a Zn-coordinating motif and an IPT/TIG (immunoglobulin-like, plexins, transcription factors/transcription factor immunoglobulin) domain, creating a unique DNA-binding domain (DBD). EBF proteins play critical roles in diverse developmental processes, including body segmentation in the Drosophila melanogaster embryo, and retina formation in mice. While several EBF family members are expressed in neuronal cells, adipocytes, and BM stroma cells, only B-lymphoid cells express EBF1. In the absence of EBF1, hematopoietic progenitor cells (HPCs) fail to activate the B-lineage program. This has been attributed to the ability of EBF1 to act as a pioneering factor with the ability to remodel chromatin, thereby creating a B-lymphoid-specific epigenetic landscape. Conditional inactivation of the Ebf1 gene in B-lineage cells has revealed additional functions of this protein in relation to the control of proliferation and apoptosis. This may explain why EBF1 is frequently targeted by mutations in human leukemia cases. This chapter provides an overview of the biochemical and functional properties of the EBF family proteins, with a focus on the roles of EBF1 in normal and malignant B-lymphocyte development.
Assuntos
Linfócitos B , Linhagem da Célula , Transativadores , Animais , Humanos , Transativadores/genética , Transativadores/metabolismo , Linfócitos B/metabolismo , Linhagem da Célula/genética , Células-Tronco Hematopoéticas/metabolismo , Fatores de Transcrição/metabolismo , Fatores de Transcrição/genéticaRESUMO
Transcriptomic data is a treasure trove in modern molecular biology, as it offers a comprehensive viewpoint into the intricate nuances of gene expression dynamics underlying biological systems. This genetic information must be utilized to infer biomolecular interaction networks that can provide insights into the complex regulatory mechanisms underpinning the dynamic cellular processes. Gene regulatory networks and protein-protein interaction networks are two major classes of such networks. This chapter thoroughly investigates the wide range of methodologies used for distilling insightful revelations from transcriptomic data that include association-based methods (based on correlation among expression vectors), probabilistic models (using Bayesian and Gaussian models), and interologous methods. We reviewed different approaches for evaluating the significance of interactions based on the network topology and biological functions of the interacting molecules and discuss various strategies for the identification of functional modules. The chapter concludes with highlighting network-based techniques of prioritizing key genes, outlining the centrality-based, diffusion- based, and subgraph-based methods. The chapter provides a meticulous framework for investigating transcriptomic data to uncover assembly of complex molecular networks for their adaptable analyses across a broad spectrum of biological domains.
Assuntos
Biologia Computacional , Perfilação da Expressão Gênica , Redes Reguladoras de Genes , Transcriptoma , Humanos , Teorema de Bayes , Biologia Computacional/métodos , Perfilação da Expressão Gênica/métodos , Mapeamento de Interação de Proteínas/métodos , Mapas de Interação de Proteínas/genéticaRESUMO
This article introduces a causal discovery method to learn nonlinear relationships in a directed acyclic graph with correlated Gaussian errors due to confounding. First, we derive model identifiability under the sublinear growth assumption. Then, we propose a novel method, named the Deconfounded Functional Structure Estimation (DeFuSE), consisting of a deconfounding adjustment to remove the confounding effects and a sequential procedure to estimate the causal order of variables. We implement DeFuSE via feedforward neural networks for scalable computation. Moreover, we establish the consistency of DeFuSE under an assumption called the strong causal minimality. In simulations, DeFuSE compares favorably against state-of-the-art competitors that ignore confounding or nonlinearity. Finally, we demonstrate the utility and effectiveness of the proposed approach with an application to gene regulatory network analysis. The Python implementation is available at https://github.com/chunlinli/defuse.
RESUMO
BACKGROUND: Proteins harboring the SPX domain are crucial for the regulation of phosphate (Pi) homeostasis in plants. This study aimed to identify and analyze the entire SPX gene family within the cucumber genome. RESULTS: The cucumber genome encompassed 16 SPX domain-containing genes, which were distributed across six chromosomes and categorized into four distinct subfamilies: SPX, SPX-MFS, SPX-EXS and SPX-RING, based on their structure characteristics. Additionally, gene duplications and synteny analysis were conducted for CsSPXs, revealing that their promoter regions were enriched with a variety of hormone-responsive, biotic/abiotic stress and typical P1BS-related elements. Tissue expression profiling of CsSPX genes revealed that certain members were specifically expressed in particular organs, suggesting essential roles in cucumber growth and development. Under low Pi stress, CsSPX1 and CsSPX2 exhibited a particularly strong response to Pi starvation. It was observed that the cucumber cultivar Xintaimici displayed greater tolerance to low Pi compared to black-spined cucumber under low Pi stress conditions. Protein interaction networks for the 16 CsSPX proteins were predicted, and yeast two-hybrid assay revealed that CsPHR1 interacted with CsSPX2, CsSPX3, CsSPX4 and CsSPX5, implying their involvement in the Pi signaling pathway in conjunction with CsPHR1. CONCLUSION: This research lays the foundation for further exploration of the function of the CsSPX genes in response to low Pi stress and for elucidating the underlying mechanism.
Assuntos
Cucumis sativus , Família Multigênica , Fósforo , Proteínas de Plantas , Cucumis sativus/genética , Proteínas de Plantas/genética , Proteínas de Plantas/metabolismo , Fósforo/metabolismo , Fósforo/deficiência , Genoma de Planta , Genes de Plantas , Regulação da Expressão Gênica de Plantas , FilogeniaRESUMO
Drafting gene regulatory networks (GRNs) requires embryological knowledge pertaining to the cell type families, information on the regulatory genes, causal data from gene knockdown experiments and validations of the identified interactions by cis-regulatory analysis. We use multi-omics involving next-generation sequencing (-seq) to obtain the necessary information for drafting Strongylocentrotus purpuratus posterior gut GRN. Here we present an update to the GRN using i) a single cell RNA-seq derived cell atlas highlighting the 2 day post fertilization (dpf) sea urchin gastrula cell type families, as well as the genes expressed at single cell level, ii) a set of putative cis-regulatory modules and transcription factor (TF) binding sites obtained from chromatin accessibility ATAC-seq data, and iii) interactions directionality obtained from differential bulk RNA-seq following knockdown of the TF Sp-Pdx1, a key regulator of gut patterning in sea urchins. Combining these datasets, we draft the GRN for the hindgut Sp-Pdx1 positive cells in the 2 dpf gastrula embryo. Overall, our data suggests the complex connectivity of the posterior gut GRN and increases the resolution of gene regulatory cascades operating within it.