Genetic Control of Expression and Splicing in Developing Human Brain Informs Disease Mechanisms.

Tissue-specific regulatory regions harbor substantial genetic risk for disease. Because brain development is a critical epoch for neuropsychiatric disease susceptibility, we characterized the genetic control of the transcriptome in 201 mid-gestational human brains, identifying 7,962 expression quantitative trait loci (eQTL) and 4,635 spliceQTL (sQTL), including several thousand prenatal-specific regulatory regions. We show that significant genetic liability for neuropsychiatric disease lies within prenatal eQTL and sQTL. Integration of eQTL and sQTL with genome-wide association studies (GWAS) via transcriptome-wide association identified dozens of novel candidate risk genes, highlighting shared and stage-specific mechanisms in schizophrenia (SCZ). Gene network analysis revealed that SCZ and autism spectrum disorder (ASD) affect distinct developmental gene co-expression modules. Yet, in each disorder, common and rare genetic variation converges within modules, which in ASD implicates superficial cortical neurons. More broadly, these data, available as a web browser and our analyses, demonstrate the genetic mechanisms by which developmental events have a widespread influence on adult anatomical and behavioral phenotypes.

Gene Essentiality Profiling Reveals Gene Networks and Synthetic Lethal Interactions with Oncogenic Ras.

The genetic dependencies of human cancers widely vary. Here, we catalog this heterogeneity and use it to identify functional gene interactions and genotype-dependent liabilities in cancer. By using genome-wide CRISPR-based screens, we generate a gene essentiality dataset across 14 human acute myeloid leukemia (AML) cell lines. Sets of genes with correlated patterns of essentiality across the lines reveal new gene relationships, the essential substrates of enzymes, and the molecular functions of uncharacterized proteins. Comparisons of differentially essential genes between Ras-dependent and -independent lines uncover synthetic lethal partners of oncogenic Ras. Screens in both human AML and engineered mouse pro-B cells converge on a surprisingly small number of genes in the Ras processing and MAPK pathways and pinpoint PREX1 as an AML-specific activator of MAPK signaling. Our findings suggest general strategies for defining mammalian gene networks and synthetic lethal interactions by exploiting the natural genetic and epigenetic diversity of human cancer cells.

Dissecting Immune Circuits by Linking CRISPR-Pooled Screens with Single-Cell RNA-Seq.

In multicellular organisms, dedicated regulatory circuits control cell type diversity and responses. The crosstalk and redundancies within these circuits and substantial cellular heterogeneity pose a major research challenge. Here, we present CRISP-seq, an integrated method for massively parallel single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-pooled screens. We show that profiling the genomic perturbation and transcriptome in the same cell enables us to simultaneously elucidate the function of multiple factors and their interactions. We applied CRISP-seq to probe regulatory circuits of innate immunity. By sampling tens of thousands of perturbed cells in vitro and in mice, we identified interactions and redundancies between developmental and signaling-dependent factors. These include opposing effects of Cebpb and Irf8 in regulating the monocyte/macrophage versus dendritic cell lineages and differential functions for Rela and Stat1/2 in monocyte versus dendritic cell responses to pathogens. This study establishes CRISP-seq as a broadly applicable, comprehensive, and unbiased approach for elucidating mammalian regulatory circuits.

Pioneer Transcription Factors Initiating Gene Network Changes.

Pioneer transcription factors have the intrinsic biochemical ability to scan partial DNA sequence motifs that are exposed on the surface of a nucleosome and thus access silent genes that are inaccessible to other transcription factors. Pioneer factors subsequently enable other transcription factors, nucleosome remodeling complexes, and histone modifiers to engage chromatin, thereby initiating the formation of an activating or repressive regulatory sequence. Thus, pioneer factors endow the competence for fate changes in embryonic development, are essential for cellular reprogramming, and rewire gene networks in cancer cells. Recent studies with reconstituted nucleosomes in vitro and chromatin binding in vivo reveal that pioneer factors can directly perturb nucleosome structure and chromatin accessibility in different ways. This review focuses on our current understanding of the mechanisms by which pioneer factors initiate gene network changes and will ultimately contribute to our ability to control cell fates at will.

Small RNAs, Large Networks: Posttranscriptional Regulons in Gram-Negative Bacteria.

Small regulatory RNA (sRNAs) are key mediators of posttranscriptional gene control in bacteria. Assisted by RNA-binding proteins, a single sRNA often modulates the expression of dozens of genes, and thus sRNAs frequently adopt central roles in regulatory networks. Posttranscriptional regulation by sRNAs comes with several unique features that cannot be achieved by transcriptional regulators. However, for optimal network performance, transcriptional and posttranscriptional control mechanisms typically go hand-in-hand. This view is reflected by the ever-growing class of mixed network motifs involving sRNAs and transcription factors, which are ubiquitous in biology and whose regulatory properties we are beginning to understand. In addition, sRNA activity can be antagonized by base-pairing with sponge RNAs, adding yet another layer of complexity to these networks. In this article, we summarize the regulatory concepts underlying sRNA-mediated gene control in bacteria and discuss how sRNAs shape the output of a network, focusing on several key examples.

A multistep computational approach reveals a neuro-mesenchymal cell population in the embryonic hematopoietic stem cell niche.

The first hematopoietic stem and progenitor cells (HSPCs) emerge in the Aorta-Gonad-Mesonephros (AGM) region of the mid-gestation mouse embryo. However, the precise nature of their supportive mesenchymal microenvironment remains largely unexplored. Here, we profiled transcriptomes of laser micro-dissected aortic tissues at three developmental stages and individual AGM cells. Computational analyses allowed the identification of several cell subpopulations within the E11.5 AGM mesenchyme, with the presence of a yet unidentified subpopulation characterized by the dual expression of genes implicated in adhesive or neuronal functions. We confirmed the identity of this cell subset as a neuro-mesenchymal population, through morphological and lineage tracing assays. Loss of function in the zebrafish confirmed that Decorin, a characteristic extracellular matrix component of the neuro-mesenchyme, is essential for HSPC development. We further demonstrated that this cell population is not merely derived from the neural crest, and hence, is a bona fide novel subpopulation of the AGM mesenchyme.

A Journey Through Genetics to Biology.

Although my engagement with human genetics emerged gradually, and sometimes serendipitously, it has held me spellbound for decades. Without my teachers, students, postdocs, colleagues, and collaborators, I would not be writing this review of my scientific adventures. Early gene and disease mapping was a satisfying puzzle-solving exercise, but building biological insight was my main goal. The project trajectory was hugely influenced by the evolutionarily conserved nature of the implicated genes and by the pace of progress in genetic technologies. The rich detail of clinical observations, particularly in eye disease, makes humans an excellent model, especially when complemented by the use of multiple other animal species for experimental validation. The contributions of collaborators and rivals also influenced our approach. We are very fortunate to work in this era of unprecedented progress in genetics and genomics.

Blood transcriptome profiling reveals distinct gene networks induced by mRNA vaccination against COVID-19.

Messenger RNA (mRNA) vaccines represent a new class of vaccines that has been shown to be highly effective during the COVID-19 pandemic and that holds great potential for other preventative and therapeutic applications. While it is known that the transcriptional activity of various genes is altered following mRNA vaccination, identifying and studying gene networks could reveal important scientific insights that might inform future vaccine designs. In this study, we conducted an in-depth weighted gene correlation network analysis of the blood transcriptome before and 24 h after the second and third vaccination with licensed mRNA vaccines against COVID-19 in humans, following a prime vaccination with either mRNA or ChAdOx1 vaccines. Utilizing this unsupervised gene network analysis approach, we identified distinct modular networks of co-varying genes characterized by either an expressional up- or downregulation in response to vaccination. Downregulated networks were associated with cell metabolic processes and regulation of transcription factors, while upregulated networks were associated with myeloid differentiation, antigen presentation, and antiviral, interferon-driven pathways. Within this interferon-associated network, we identified highly connected hub genes such as STAT2 and RIGI and associated upstream transcription factors, potentially playing important regulatory roles in the vaccine-induced immune response. The expression profile of this network significantly correlated with S1-specific IgG levels at the follow-up visit in vaccinated individuals. Those findings could be corroborated in a second, independent cohort of mRNA vaccine recipients. Collectively, results from this modular gene network analysis enhance the understanding of mRNA vaccines from a systems immunology perspective. Influencing specific gene networks could lead to optimized vaccines that elicit augmented vaccine responses.

Regeneration among animals: An evolutionary hypothesis related to aquatic versus terrestrial environment.

The present hypothesis tries to explain animal regeneration in relation to the life cycles and environment of different animals. Regeneration is a basic phenomenon present since the origin of life in the sea, as testimonial in lower or more complex extant marine animals. Aquatic animals that evolved an indirect development, forming larvae and transiting into the adult stage through metamorphosis, use gene networks present in their genome for these transformations. In case of injury or organ loss as adults, they can re-utilize most or part of the gene networks previously activated during larval growth and metamorphosis. In contrast, terrestrial animals that evolved life cycles with the elimination of larvae and metamorphosis for the adaptation to land conditions lost some of the genes implicated in these post-developmental processes and consequently also the ability to regenerate. Few arthropods and lizards are capable to form hydrated regenerative blastemas with a similar consistence of embryonic tissues. The present hypothesis submits that regeneration cannot be activated in the dry land environment and consequently was largely or completely abolished in terrestrial animals. After injury or organ loss, nematodes, most arthropods and terrestrial vertebrates can only form scars or a limited healing or regengrow in juveniles. This is a process where somatic growth is superimposed to wound healing so that the apparent regeneration derives from the combination from both processes. When full growth is terminated these terrestrial animals can only heal by scarring.

Comparative and integrative single cell analysis reveals new insights into the transcriptional immaturity of stem cell-derived ß cells.

Diabetes cell replacement therapy has the potential to be transformed by human pluripotent stem cell-derived ß cells (SC-ß cells). However, the precise identity of SC-ß cells in relationship to primary fetal and adult ß-cells remains unclear. Here, we used single-cell sequencing datasets to characterize the transcriptional identity of islets from in vitro differentiation, fetal islets, and adult islets. Our analysis revealed that SC-ß cells share a core ß-cell transcriptional identity with human adult and fetal ß-cells, however SC-ß cells possess a unique transcriptional profile characterized by the persistent expression and activation of progenitor and neural-biased gene networks. These networks are present in SC-ß cells, irrespective of the derivation protocol used. Notably, fetal ß-cells also exhibit this neural signature at the transcriptional level. Our findings offer insights into the transcriptional identity of SC-ß cells and underscore the need for further investigation of the role of neural transcriptional networks in their development.

Ovarian cancer classification and prognosis assessment model based on prognostic target genes in key microRNA-target gene networks.

BACKGROUND: The present study was designed to screen key microRNA (miRNA)-target gene networks for ovarian cancer (OC) and to classify and construct a risk assessment system for OC based on the target genes. METHODS: OC sample data of The Cancer Genome Atlas dataset and GSE26193, GSE30161, GSE63885 and GSE9891 datasets were retrospectively collected. Pearson correlation analysis and targeted analysis of miRNA and target gene were performed to screen key miRNA-target gene networks. Target genes associated with the prognosis of OC were screened from key miRNA-target gene networks for consensus clustering and least absolute shrinkage and selection operator-based regression machine learning analysis of OC samples. RESULTS: Twenty target genes of 2651 key miRNA-target gene pairs had significant prognostic correlation in each OC cohort, and OC was divided into three clusters. There were differences in prognostic outcome, biological pathways, immune cell abundance and susceptibility to immune checkpoint blockade (ICB) therapy and anti-tumor drugs among the three molecular clusters. S2 exhibited the least advantage in prognosis and immunotherapy response rate in the three molecular clusters, and the pathways regulating immunity, hypoxia, metabolism and promoting malignant progression of cancer, as well as infiltrating immune and stromal cell population abundance, were the highest in this cluster. An eight-target gene prognostic model was created, and the risk index obtained by using this model not only significantly distinguished the immune characteristics of the sample, but also predicted the response of the sample to ICB treatment, and helped to screen 36 potential anti-OC drugs. CONCLUSIONS: The present study provides a classification strategy for OC based on prognostic target genes in key miRNA-target gene networks, and creates a risk assessment system for predicting prognosis and response to ICB therapy in OC patients, providing molecular basis for prognosis and precise treatment of OC.

A developmental lineage-based gene co-expression network for mouse pancreatic ß-cells reveals a role for Zfp800 in pancreas development.

To gain a deeper understanding of pancreatic ß-cell development, we used iterative weighted gene correlation network analysis to calculate a gene co-expression network (GCN) from 11 temporally and genetically defined murine cell populations. The GCN, which contained 91 distinct modules, was then used to gain three new biological insights. First, we found that the clustered protocadherin genes are differentially expressed during pancreas development. Pcdhγ genes are preferentially expressed in pancreatic endoderm, Pcdhß genes in nascent islets, and Pcdhα genes in mature ß-cells. Second, after extracting sub-networks of transcriptional regulators for each developmental stage, we identified 81 zinc finger protein (ZFP) genes that are preferentially expressed during endocrine specification and ß-cell maturation. Third, we used the GCN to select three ZFPs for further analysis by CRISPR mutagenesis of mice. Zfp800 null mice exhibited early postnatal lethality, and at E18.5 their pancreata exhibited a reduced number of pancreatic endocrine cells, alterations in exocrine cell morphology, and marked changes in expression of genes involved in protein translation, hormone secretion and developmental pathways in the pancreas. Together, our results suggest that developmentally oriented GCNs have utility for gaining new insights into gene regulation during organogenesis.

Water deficit response in nodulated soybean roots: a comprehensive transcriptome and translatome network analysis.

BACKGROUND: Soybean establishes a mutualistic interaction with nitrogen-fixing rhizobacteria, acquiring most of its nitrogen requirements through symbiotic nitrogen fixation. This crop is susceptible to water deficit; evidence suggests that its nodulation status-whether it is nodulated or not-can influence how it responds to water deficit. The translational control step of gene expression has proven relevant in plants subjected to water deficit. RESULTS: Here, we analyzed soybean roots' differential responses to water deficit at transcriptional, translational, and mixed (transcriptional + translational) levels. Thus, the transcriptome and translatome of four combined-treated soybean roots were analyzed. We found hormone metabolism-related genes among the differentially expressed genes (DEGs) at the translatome level in nodulated and water-restricted plants. Also, weighted gene co-expression network analysis followed by differential expression analysis identified gene modules associated with nodulation and water deficit conditions. Protein-protein interaction network analysis was performed for subsets of mixed DEGs of the modules associated with the plant responses to nodulation, water deficit, or their combination. CONCLUSIONS: Our research reveals that the stand-out processes and pathways in the before-mentioned plant responses partially differ; terms related to glutathione metabolism and hormone signal transduction (2 C protein phosphatases) were associated with the response to water deficit, terms related to transmembrane transport, response to abscisic acid, pigment metabolic process were associated with the response to nodulation plus water deficit. Still, two processes were common: galactose metabolism and branched-chain amino acid catabolism. A comprehensive analysis of these processes could lead to identifying new sources of tolerance to drought in soybean.

Comparative transcriptomic analysis unveils candidate genes associated with sugarcane growth rate.

Sugarcane (Saccharum spp.) growth is regulated by intricate gene networks and hormone secretions, positively correlating with sugarcane yield. There is a rising interest in exploring how the candidate genes found in sugarcane respond to plant growth. In this study, we simulated a typical growth environment to obtain accurate phenotypic data and screened for potential genes associated with plant growth through transcriptomics. Compared to Saccharum GuiTang 42, the other variety Saccharum GuiTang 44 exhibited earlier germination, a higher emergence rate, thicker pseudostems, taller plants, and a more extensive root system. The middle buds formed the greatest number of roots, followed by the lower and upper buds. Indole-3-acetic acid (IAA) and jasmonic acid effectively promoted bud development, while abscisic acid and trans-zeatin exhibited negative correlations with sugarcane bud growth. Transcriptome data from the upper, middle, and lower buds revealed 24,158 differentially expressed genes in all three comparisons, with MAPK signaling emerging as a critical pathway. The photosynthesis-antenna protein pathway is vital for middle and lower bud development during root germination. Lastly, key gene modules related to differences in hormone content between the two varieties were defined through weighted correlation network analysis and identified. The module significantly associated with IAA was enriched in pathways such as Proteasome and Protein processing in the endoplasmic reticulum, and the upregulation of key genes involved in this gene module had a highly significant positive correlation with bud outgrowth combined with IAA secretion. In conclusion, we have elucidated the pathways of hormones during sugarcane growth and the interactions between IAA and critical genes. These in-depth findings may guide modern sugarcane breeding.

Providing biological context for GWAS results using eQTL regulatory and co-expression networks in Populus.

Our study utilized genome-wide association studies (GWAS) to link nucleotide variants to traits in Populus trichocarpa, a species with rapid linkage disequilibrium decay. The aim was to overcome the challenge of interpreting statistical associations at individual loci without sufficient biological context, which often leads to reliance solely on gene annotations from unrelated model organisms. We employed an integrative approach that included GWAS targeting multiple traits using three individual techniques for lignocellulose phenotyping, expression quantitative trait loci (eQTL) analysis to construct transcriptional regulatory networks around each candidate locus and co-expression analysis to provide biological context for these networks, using lignocellulose biosynthesis in Populus trichocarpa as a case study. The research identified three candidate genes potentially involved in lignocellulose formation, including one previously recognized gene (Potri.005G116800/VND1, a critical regulator of secondary cell wall formation) and two genes (Potri.012G130000/AtSAP9 and Potri.004G202900/BIC1) with newly identified putative roles in lignocellulose biosynthesis. Our integrative approach offers a framework for providing biological context to loci associated with trait variation, facilitating the discovery of new genes and regulatory networks.

Hybridization and gene expression: Beyond differentially expressed genes.

Gene expression has a key role in reproductive isolation, and studies of hybrid gene expression have identified mechanisms causing hybrid sterility. Here, we review the evidence for altered gene expression following hybridization and outline the mechanisms shown to contribute to altered gene expression in hybrids. Transgressive gene expression, transcending that of both parental species, is pervasive in early generation sterile hybrids, but also frequently observed in viable, fertile hybrids. We highlight studies showing that hybridization can result in transgressive gene expression, also in established hybrid lineages or species. Such extreme patterns of gene expression in stabilized hybrid taxa suggest that altered hybrid gene expression may result in hybridization-derived evolutionary novelty. We also conclude that while patterns of misexpression in hybrids are well documented, the understanding of the mechanisms causing misexpression is lagging. We argue that jointly assessing differences in cell composition and cell-specific changes in gene expression in hybrids, in addition to assessing changes in chromatin and methylation, will significantly advance our understanding of the basis of altered gene expression. Moreover, uncovering to what extent evolution of gene expression results in altered expression for individual genes, or entire networks of genes, will advance our understanding of how selection moulds gene expression. Finally, we argue that jointly studying the dual roles of altered hybrid gene expression, serving both as a mechanism for reproductive isolation and as a substrate for hybrid ecological adaptation, will lead to significant advances in our understanding of the evolution of gene expression.

Cuticle development and the underlying transcriptome-metabolome associations during early seedling establishment.

The plant cuticle is a complex extracellular lipid barrier that has multiple protective functions. This study investigated cuticle deposition by integrating metabolomics and transcriptomics data gathered from six different maize seedling organs of four genotypes, the inbred lines B73 and Mo17, and their reciprocal hybrids. These datasets captured the developmental transition of the seedling from heterotrophic skotomorphogenic growth to autotrophic photomorphogenic growth, a transition that is highly vulnerable to environmental stresses. Statistical interrogation of these data revealed that the predominant determinant of cuticle composition is seedling organ type, whereas the seedling genotype has a smaller effect on this phenotype. Gene-to-metabolite associations assessed by integrated statistical analyses identified three gene networks associated with the deposition of different elements of the cuticle: cuticular waxes; monomers of lipidized cell wall biopolymers, including cutin and suberin; and both of these elements. These gene networks reveal three metabolic programs that appear to support cuticle deposition, including processes of chloroplast biogenesis, lipid metabolism, and molecular regulation (e.g. transcription factors, post-translational regulators, and phytohormones). This study demonstrates the wider physiological metabolic context that can determine cuticle deposition and lays the groundwork for new targets for modulating the properties of this protective barrier.

Changes of Transcriptomic Activity in Rat Brain Cells under the Influence of Synthetic Adrenocorticotropic Hormone-Like Peptides.

Synthetic peptides have a wide range of clinical effects. Of particular interest are peptides based on adrenocorticotropic hormone (ACTH) both as already used and as potential drugs for preventing consequences of cerebral ischemia. However, it is necessary to study influence of the peptide on the brain cells under normal physiological conditions, including understanding the risks of their use. Here, we used high-throughput RNA sequencing (RNA-Seq) to identify differentially expressed genes (DEGs) in the brain frontal cortex of rat receiving intraperitoneal administration of ACTH-like peptides ACTH(4-7)PGP (Semax) and ACTH(6-9)PGP, or saline. We identified 258 and 228 DEGs, respectively, with the fold change > 1.5 and Padj < 0.05 at 22.5 h after the first administration of Semax and ACTH(6-9)PGP. Metabolic pathways, characterizing both common and specific effects of the peptides on the transcriptome were identified. Both peptides predominantly caused decrease in expression of the genes associated with the immune system. At the same time, when comparing the effects of ACTH(6-9)PGP relative to Semax, DEGs were identified that characterized the main differences in the effects of the peptides. These genes were mostly downregulated and associated with neurosignaling systems and regulation of ion channels, thus characterizing differences in the effects of the peptides. Our data show how differences in the structure of ACTH derivatives are associated with the changes in the brain cell transcriptome following exposure to these related peptides. Furthermore, our results demonstrate that when studying influence of regulatory peptides on transcriptome under pathological conditions, it is necessary to take into account their actions under normal physiological conditions.

Spatial biology of Ising-like synthetic genetic networks.

BACKGROUND: Understanding how spatial patterns of gene expression emerge from the interaction of individual gene networks is a fundamental challenge in biology. Developing a synthetic experimental system with a common theoretical framework that captures the emergence of short- and long-range spatial correlations (and anti-correlations) from interacting gene networks could serve to uncover generic scaling properties of these ubiquitous phenomena. RESULTS: Here, we combine synthetic biology, statistical mechanics models, and computational simulations to study the spatial behavior of synthetic gene networks (SGNs) in Escherichia coli quasi-2D colonies growing on hard agar surfaces. Guided by the combined mechanisms of the contact process lattice simulation and two-dimensional Ising model (CPIM), we describe the spatial behavior of bi-stable and chemically coupled SGNs that self-organize into patterns of long-range correlations with power-law scaling or short-range anti-correlations. These patterns, resembling ferromagnetic and anti-ferromagnetic configurations of the Ising model near critical points, maintain their scaling properties upon changes in growth rate and cell shape. CONCLUSIONS: Our findings shed light on the spatial biology of coupled and bistable gene networks in growing cell populations. This emergent spatial behavior could provide insights into the study and engineering of self-organizing gene patterns in eukaryotic tissues and bacterial consortia.

Resequencing of Brazilian locally adapted cattle breeds revealed variants in candidate genes and transcription factors for meat fatty acid profile.

The beef cattle industry has experienced a shift driven by a market demand for healthier meat, cost efficiency and environmental sustainability in recent years. Consequently, there has been a growing focus on the fatty acids content and functions of meat in cattle breeding programmes. Besides, a deeper understanding of the biological mechanisms influencing the expression of different phenotypes related to fatty acid profiles is crucial. In this study, we aimed to identify Single-Nucleotide Variants (SNV) and Insertion/Deletion (InDels) DNA variants in candidate genes related to fatty acid profiles described in genomic, transcriptomic and proteomic studies conducted in beef cattle breeds. Utilizing whole-genome re-sequencing data from Brazilian locally adapted bovine breeds, namely Caracu and Pantaneiro, we identified SNVs and InDels associated with 23,947 genes. From these, we identified 318 candidate genes related to fatty acid profiles that contain variants. Subsequently, we select only genes with SNVs and InDels in their promoter, 5' UTR and coding region. Through the gene-biological process network, approximately 19 genes were highlighted. Furthermore, considering the studied trait and a literature review, we selected the main transcription factors (TF). Functional analysis via gene-TF network allowed us to identify the 30 most likely candidate genes for meat fatty acid profile in cattle. LIPE, MFSD2A and SREBF1 genes were highlighted in networks due to their biological importance. Further dissection of these genes revealed 15 new variants found in promoter regions of Caracu and Pantaneiro sequences. The gene networks facilitated a better functional understanding of genes and TF, enabling the identification of variants potentially related to the expression of candidate genes for meat fatty acid profiles in cattle.