Comparative analysis of gene regulatory networks: from network reconstruction to evolution.

Regulation of gene expression is central to many biological processes. Although reconstruction of regulatory circuits from genomic data alone is therefore desirable, this remains a major computational challenge. Comparative approaches that examine the conservation and divergence of circuits and their components across strains and species can help reconstruct circuits as well as provide insights into the evolution of gene regulatory processes and their adaptive contribution. In recent years, advances in genomic and computational tools have led to a wealth of methods for such analysis at the sequence, expression, pathway, module, and entire network level. Here, we review computational methods developed to study transcriptional regulatory networks using comparative genomics, from sequence to functional data. We highlight how these methods use evolutionary conservation and divergence to reliably detect regulatory components as well as estimate the extent and rate of divergence. Finally, we discuss the promise and open challenges in linking regulatory divergence to phenotypic divergence and adaptation.

Dynamic regulatory module networks for inference of cell type-specific transcriptional networks.

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type-specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type-specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type-specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus.

Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials - a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

Human left ventricular cardiac fibroblasts undergo a dynamic shift in secretome and gene expression toward a cardiac myofibroblast phenotype during early passage in typical culture expansion conditions.

Cardiac fibroblasts (CFs) are critical components of the cardiac niche and primarily responsible for assembly and maintenance of the cardiac extracellular matrix (ECM). CFs are increasingly of interest for tissue engineering and drug development applications, as they provide synergistic support to cardiomyocytes through direct cell-to-cell signaling and cell-to-ECM interactions via soluble factors, including cytokines, growth factors and extracellular vesicles. CFs can be activated to a cardiac myofibroblast (CMF) phenotype upon injury or stimulation with transforming growth factor beta 1. Once activated, CMFs assemble collagen-rich ECM, which is vitally important to stabilize scar formation following myocardial infarction, for example. Although there is greater experience with culture expansion of CFs among non-human strains, very little is known about human CF-to-CMF transitions and expression patterns during culture expansion. In this study, we evaluated for shifts in inflammatory and angiogenic expression profiles of human CFs in typical culture expansion conditions. Understanding shifts in cellular expression patterns during CF culture expansion is critically important to establish quality benchmarks and optimize large-scale manufacturing for future clinical applications.

Leveraging epigenomes and three-dimensional genome organization for interpreting regulatory variation.

Understanding the impact of regulatory variants on complex phenotypes is a significant challenge because the genes and pathways that are targeted by such variants and the cell type context in which regulatory variants operate are typically unknown. Cell-type-specific long-range regulatory interactions that occur between a distal regulatory sequence and a gene offer a powerful framework for examining the impact of regulatory variants on complex phenotypes. However, high-resolution maps of such long-range interactions are available only for a handful of cell types. Furthermore, identifying specific gene subnetworks or pathways that are targeted by a set of variants is a significant challenge. We have developed L-HiC-Reg, a Random Forests regression method to predict high-resolution contact counts in new cell types, and a network-based framework to identify candidate cell-type-specific gene networks targeted by a set of variants from a genome-wide association study (GWAS). We applied our approach to predict interactions in 55 Roadmap Epigenomics Mapping Consortium cell types, which we used to interpret regulatory single nucleotide polymorphisms (SNPs) in the NHGRI-EBI GWAS catalogue. Using our approach, we performed an in-depth characterization of fifteen different phenotypes including schizophrenia, coronary artery disease (CAD) and Crohn's disease. We found differentially wired subnetworks consisting of known as well as novel gene targets of regulatory SNPs. Taken together, our compendium of interactions and the associated network-based analysis pipeline leverages long-range regulatory interactions to examine the context-specific impact of regulatory variation in complex phenotypes.

Transcriptome-wide transmission disequilibrium analysis identifies novel risk genes for autism spectrum disorder.

Recent advances in consortium-scale genome-wide association studies (GWAS) have highlighted the involvement of common genetic variants in autism spectrum disorder (ASD), but our understanding of their etiologic roles, especially the interplay with rare variants, is incomplete. In this work, we introduce an analytical framework to quantify the transmission disequilibrium of genetically regulated gene expression from parents to offspring. We applied this framework to conduct a transcriptome-wide association study (TWAS) on 7,805 ASD proband-parent trios, and replicated our findings using 35,740 independent samples. We identified 31 associations at the transcriptome-wide significance level. In particular, we identified POU3F2 (p = 2.1E-7), a transcription factor mainly expressed in developmental brain. Gene targets regulated by POU3F2 showed a 2.7-fold enrichment for known ASD genes (p = 2.0E-5) and a 2.7-fold enrichment for loss-of-function de novo mutations in ASD probands (p = 7.1E-5). These results provide a novel connection between rare and common variants, whereby ASD genes affected by very rare mutations are regulated by an unlinked transcription factor affected by common genetic variations.

Intracerebroventricular dosing of N-sulfoglucosamine sulfohydrolase in mucopolysaccharidosis IIIA mice reduces markers of brain lysosomal dysfunction.

Mucopolysaccharidosis type IIIA (MPS IIIA) is a lysosomal storage disorder caused by N-sulfoglucosamine sulfohydrolase (SGSH) deficiency. SGSH removes the sulfate from N-sulfoglucosamine residues on the nonreducing end of heparan sulfate (HS-NRE) within lysosomes. Enzyme deficiency results in accumulation of partially degraded HS within lysosomes throughout the body, leading to a progressive severe neurological disease. Enzyme replacement therapy has been proposed, but further evaluation of the treatment strategy is needed. Here, we used Chinese hamster ovary cells to produce a highly soluble and fully active recombinant human sulfamidase (rhSGSH). We discovered that rhSGSH utilizes both the CI-MPR and LRP1 receptors for uptake into patient fibroblasts. A single intracerebroventricular (ICV) injection of rhSGSH in MPS IIIA mice resulted in a tissue half-life of 9 days and widespread distribution throughout the brain. Following a single ICV dose, both total HS and the MPS IIIA disease-specific HS-NRE were dramatically reduced, reaching a nadir 2 weeks post dose. The durability of effect for reduction of both substrate and protein markers of lysosomal dysfunction and a neuroimmune response lasted through the 56 days tested. Furthermore, seven weekly 148 µg doses ICV reduced those markers to near normal and produced a 99.5% reduction in HS-NRE levels. A pilot study utilizing every other week dosing in two animals supports further evaluation of less frequent dosing. Finally, our dose-response study also suggests lower doses may be efficacious. Our findings show that rhSGSH can normalize lysosomal HS storage and markers of a neuroimmune response when delivered ICV.

Evolution of miRNA-Binding Sites and Regulatory Networks in Cichlids.

The divergence of regulatory regions and gene regulatory network (GRN) rewiring is a key driver of cichlid phenotypic diversity. However, the contribution of miRNA-binding site turnover has yet to be linked to GRN evolution across cichlids. Here, we extend our previous studies by analyzing the selective constraints driving evolution of miRNA and transcription factor (TF)-binding sites of target genes, to infer instances of cichlid GRN rewiring associated with regulatory binding site turnover. Comparative analyses identified increased species-specific networks that are functionally associated to traits of cichlid phenotypic diversity. The evolutionary rewiring is associated with differential models of miRNA- and TF-binding site turnover, driven by a high proportion of fast-evolving polymorphic sites in adaptive trait genes compared with subsets of random genes. Positive selection acting upon discrete mutations in these regulatory regions is likely to be an important mechanism in rewiring GRNs in rapidly radiating cichlids. Regulatory variants of functionally associated miRNA- and TF-binding sites of visual opsin genes differentially segregate according to phylogeny and ecology of Lake Malawi species, identifying both rewired, for example, clade-specific and conserved network motifs of adaptive trait associated GRNs. Our approach revealed several novel candidate regulators, regulatory regions, and three-node motifs across cichlid genomes with previously reported associations to known adaptive evolutionary traits.

Human iPSC Modeling Reveals Mutation-Specific Responses to Gene Therapy in a Genotypically Diverse Dominant Maculopathy.

Dominantly inherited disorders are not typically considered to be therapeutic candidates for gene augmentation. Here, we utilized induced pluripotent stem cell-derived retinal pigment epithelium (iPSC-RPE) to test the potential of gene augmentation to treat Best disease, a dominant macular dystrophy caused by over 200 missense mutations in BEST1. Gene augmentation in iPSC-RPE fully restored BEST1 calcium-activated chloride channel activity and improved rhodopsin degradation in an iPSC-RPE model of recessive bestrophinopathy as well as in two models of dominant Best disease caused by different mutations in regions encoding ion-binding domains. A third dominant Best disease iPSC-RPE model did not respond to gene augmentation, but showed normalization of BEST1 channel activity following CRISPR-Cas9 editing of the mutant allele. We then subjected all three dominant Best disease iPSC-RPE models to gene editing, which produced premature stop codons specifically within the mutant BEST1 alleles. Single-cell profiling demonstrated no adverse perturbation of retinal pigment epithelium (RPE) transcriptional programs in any model, although off-target analysis detected a silent genomic alteration in one model. These results suggest that gene augmentation is a viable first-line approach for some individuals with dominant Best disease and that non-responders are candidates for alternate approaches such as gene editing. However, testing gene editing strategies for on-target efficiency and off-target events using personalized iPSC-RPE model systems is warranted. In summary, personalized iPSC-RPE models can be used to select among a growing list of gene therapy options to maximize safety and efficacy while minimizing time and cost. Similar scenarios likely exist for other genotypically diverse channelopathies, expanding the therapeutic landscape for affected individuals.

Identification of FMR1-regulated molecular networks in human neurodevelopment.

RNA-binding proteins (RNA-BPs) play critical roles in development and disease to regulate gene expression. However, genome-wide identification of their targets in primary human cells has been challenging. Here, we applied a modified CLIP-seq strategy to identify genome-wide targets of the FMRP translational regulator 1 (FMR1), a brain-enriched RNA-BP, whose deficiency leads to Fragile X Syndrome (FXS), the most prevalent inherited intellectual disability. We identified FMR1 targets in human dorsal and ventral forebrain neural progenitors and excitatory and inhibitory neurons differentiated from human pluripotent stem cells. In parallel, we measured the transcriptomes of the same four cell types upon FMR1 gene deletion. We discovered that FMR1 preferentially binds long transcripts in human neural cells. FMR1 targets include genes unique to human neural cells and associated with clinical phenotypes of FXS and autism. Integrative network analysis using graph diffusion and multitask clustering of FMR1 CLIP-seq and transcriptional targets reveals critical pathways regulated by FMR1 in human neural development. Our results demonstrate that FMR1 regulates a common set of targets among different neural cell types but also operates in a cell type-specific manner targeting distinct sets of genes in human excitatory and inhibitory neural progenitors and neurons. By defining molecular subnetworks and validating specific high-priority genes, we identify novel components of the FMR1 regulation program. Our results provide new insights into gene regulation by a critical neuronal RNA-BP in human neurodevelopment.

An LCO-responsive homolog of NODULE INCEPTION positively regulates lateral root formation in Populus sp.

The transcription factor NODULE INCEPTION (NIN) has been studied extensively for its multiple roles in root nodule symbiosis within plants of the nitrogen-fixing clade (NFC) that associate with soil bacteria, such as rhizobia and Frankia. However, NIN homologs are present in plants outside the NFC, suggesting a role in other developmental processes. Here, we show that the biofuel crop Populus sp., which is not part of the NFC, contains eight copies of NIN with diversified protein sequence and expression patterns. Lipo-chitooligosaccharides (LCOs) are produced by rhizobia and a wide range of fungi, including mycorrhizal ones, and act as symbiotic signals that promote lateral root formation. RNAseq analysis of Populus sp. treated with purified LCO showed induction of the PtNIN2 subfamily. Moreover, the expression of PtNIN2b correlated with the formation of lateral roots and was suppressed by cytokinin treatment. Constitutive expression of PtNIN2b overcame the inhibition of lateral root development by cytokinin under high nitrate conditions. Lateral root induction in response to LCOs likely represents an ancestral function of NIN retained and repurposed in nodulating plants, as we demonstrate that the role of NIN in LCO-induced root branching is conserved in both Populus sp. and legumes. We further established a visual marker of LCO perception in Populus sp. roots, the putative sulfotransferase PtSS1 that can be used to study symbiotic interactions with the bacterial and fungal symbionts of Populus sp.

A network-based comparative framework to study conservation and divergence of proteomes in plant phylogenies.

Comparative functional genomics offers a powerful approach to study species evolution. To date, the majority of these studies have focused on the transcriptome in mammalian and yeast phylogenies. Here, we present a novel multi-species proteomic dataset and a computational pipeline to systematically compare the protein levels across multiple plant species. Globally we find that protein levels diverge according to phylogenetic distance but is more constrained than the mRNA level. Module-level comparative analysis of groups of proteins shows that proteins that are more highly expressed tend to be more conserved. To interpret the evolutionary patterns of conservation and divergence, we develop a novel network-based integrative analysis pipeline that combines publicly available transcriptomic datasets to define co-expression modules. Our analysis pipeline can be used to relate the changes in protein levels to different species-specific phenotypic traits. We present a case study with the rhizobia-legume symbiosis process that supports the role of autophagy in this symbiotic association.

Temporal change in chromatin accessibility predicts regulators of nodulation in Medicago truncatula.

BACKGROUND: Symbiotic associations between bacteria and leguminous plants lead to the formation of root nodules that fix nitrogen needed for sustainable agricultural systems. Symbiosis triggers extensive genome and transcriptome remodeling in the plant, yet an integrated understanding of the extent of chromatin changes and transcriptional networks that functionally regulate gene expression associated with symbiosis remains poorly understood. In particular, analyses of early temporal events driving this symbiosis have only captured correlative relationships between regulators and targets at mRNA level. Here, we characterize changes in transcriptome and chromatin accessibility in the model legume Medicago truncatula, in response to rhizobial signals that trigger the formation of root nodules. RESULTS: We profiled the temporal chromatin accessibility (ATAC-seq) and transcriptome (RNA-seq) dynamics of M. truncatula roots treated with bacterial small molecules called lipo-chitooligosaccharides that trigger host symbiotic pathways of nodule development. Using a novel approach, dynamic regulatory module networks, we integrated ATAC-seq and RNA-seq time courses to predict cis-regulatory elements and transcription factors that most significantly contribute to transcriptomic changes associated with symbiosis. Regulators involved in auxin (IAA4-5, SHY2), ethylene (EIN3, ERF1), and abscisic acid (ABI5) hormone response, as well as histone and DNA methylation (IBM1), emerged among those most predictive of transcriptome dynamics. RNAi-based knockdown of EIN3 and ERF1 reduced nodule number in M. truncatula validating the role of these predicted regulators in symbiosis between legumes and rhizobia. CONCLUSIONS: Our transcriptomic and chromatin accessibility datasets provide a valuable resource to understand the gene regulatory programs controlling the early stages of the dynamic process of symbiosis. The regulators identified provide potential targets for future experimental validation, and the engineering of nodulation in species is unable to establish that symbiosis naturally.

Functional and comparative genomics reveals conserved noncoding sequences in the nitrogen-fixing clade.

Nitrogen is one of the most inaccessible plant nutrients, but certain species have overcome this limitation by establishing symbiotic interactions with nitrogen-fixing bacteria in the root nodule. This root-nodule symbiosis (RNS) is restricted to species within a single clade of angiosperms, suggesting a critical, but undetermined, evolutionary event at the base of this clade. To identify putative regulatory sequences implicated in the evolution of RNS, we evaluated the genomes of 25 species capable of nodulation and identified 3091 conserved noncoding sequences (CNS) in the nitrogen-fixing clade (NFC). We show that the chromatin accessibility of 452 CNS correlates significantly with the regulation of genes responding to lipochitooligosaccharides in Medicago truncatula. These included 38 CNS in proximity to 19 known genes involved in RNS. Five such regions are upstream of MtCRE1, Cytokinin Response Element 1, required to activate a suite of downstream transcription factors necessary for nodulation in M. truncatula. Genetic complementation of an Mtcre1 mutant showed a significant decrease of nodulation in the absence of the five CNS, when they are driving the expression of a functional copy of MtCRE1. CNS identified in the NFC may harbor elements required for the regulation of genes controlling RNS in M. truncatula.

Cultured cardiac fibroblasts and myofibroblasts express Sushi Containing Domain 2 and assemble a unique fibronectin rich matrix.

Cardiac fibroblasts and myofibroblasts assemble and maintain extracellular matrix during normal development and following injury. Culture expansion of these cells yield a bioengineered matrix that could lead to intriguing therapeutic opportunities. For example, we reported that cultured rat cardiac fibroblasts form a matrix that can be used to delivery therapeutic stem cells. Furthermore, we reported that matrix derived from cultured human cardiac fibroblasts/myofibroblasts converted monocytes into macrophages that express interesting anti-inflammatory and pro-angiogenic properties. Expanding these matrix investigations require characterization of the source cells for quality control. In these efforts, we observed and herein report that Sushi Containing Domain 2 (SUSD2) is a novel and consistent marker for cultured human cardiac fibroblast and myofibroblasts.

Intracerebroventricular enzyme replacement therapy with ß-galactosidase reverses brain pathologies due to GM1 gangliosidosis in mice.

Autosomal recessive mutations in the galactosidase ß1 (GLB1) gene cause lysosomal ß-gal deficiency, resulting in accumulation of galactose-containing substrates and onset of the progressive and fatal neurodegenerative lysosomal storage disease, GM1 gangliosidosis. Here, an enzyme replacement therapy (ERT) approach in fibroblasts from GM1 gangliosidosis patients with recombinant human ß-gal (rhß-gal) produced in Chinese hamster ovary cells enabled direct and precise rhß-gal delivery to acidified lysosomes. A single, low dose (3 nm) of rhß-gal was sufficient for normalizing ß-gal activity and mediating substrate clearance for several weeks. We found that rhß-gal uptake by the fibroblasts is dose-dependent and saturable and can be competitively inhibited by mannose 6-phosphate, suggesting cation-independent, mannose 6-phosphate receptor-mediated endocytosis from the cell surface. A single intracerebroventricularly (ICV) administered dose of rhß-gal (100 µg) resulted in broad bilateral biodistribution of rhß-gal to critical regions of pathology in a mouse model of GM1 gangliosidosis. Weekly ICV dosing of rhß-gal for 8 weeks substantially reduced brain levels of ganglioside and oligosaccharide substrates and reversed well-established secondary neuropathology. Of note, unlike with the ERT approach, chronic lentivirus-mediated GLB1 overexpression in the GM1 gangliosidosis patient fibroblasts caused accumulation of a prelysosomal pool of ß-gal, resulting in activation of the unfolded protein response and endoplasmic reticulum stress. This outcome was unsurprising in light of our in vitro biophysical findings for rhß-gal, which include pH-dependent and concentration-dependent stability and dynamic self-association. Collectively, our results highlight that ICV-ERT is an effective therapeutic intervention for managing GM1 gangliosidosis potentially more safely than with gene therapy approaches.

Calorie restriction and SIRT3 trigger global reprogramming of the mitochondrial protein acetylome.

Calorie restriction (CR) extends life span in diverse species. Mitochondria play a key role in CR adaptation; however, the molecular details remain elusive. We developed and applied a quantitative mass spectrometry method to probe the liver mitochondrial acetyl-proteome during CR versus control diet in mice that were wild-type or lacked the protein deacetylase SIRT3. Quantification of 3,285 acetylation sites-2,193 from mitochondrial proteins-rendered a comprehensive atlas of the acetyl-proteome and enabled global site-specific, relative acetyl occupancy measurements between all four experimental conditions. Bioinformatic and biochemical analyses provided additional support for the effects of specific acetylation on mitochondrial protein function. Our results (1) reveal widespread reprogramming of mitochondrial protein acetylation in response to CR and SIRT3, (2) identify three biochemically distinct classes of acetylation sites, and (3) provide evidence that SIRT3 is a prominent regulator in CR adaptation by coordinately deacetylating proteins involved in diverse pathways of metabolism and mitochondrial maintenance.

Imputed gene associations identify replicable trans-acting genes enriched in transcription pathways and complex traits.

Regulation of gene expression is an important mechanism through which genetic variation can affect complex traits. A substantial portion of gene expression variation can be explained by both local (cis) and distal (trans) genetic variation. Much progress has been made in uncovering cis-acting expression quantitative trait loci (cis-eQTL), but trans-eQTL have been more difficult to identify and replicate. Here we take advantage of our ability to predict the cis component of gene expression coupled with gene mapping methods such as PrediXcan to identify high confidence candidate trans-acting genes and their targets. That is, we correlate the cis component of gene expression with observed expression of genes in different chromosomes. Leveraging the shared cis-acting regulation across tissues, we combine the evidence of association across all available Genotype-Tissue Expression Project tissues and find 2,356 trans-acting/target gene pairs with high mappability scores. Reassuringly, trans-acting genes are enriched in transcription and nucleic acid binding pathways and target genes are enriched in known transcription factor binding sites. Interestingly, trans-acting genes are more significantly associated with selected complex traits and diseases than target or background genes, consistent with percolating trans effects. Our scripts and summary statistics are publicly available for future studies of trans-acting gene regulation.

Chromatin module inference on cellular trajectories identifies key transition points and poised epigenetic states in diverse developmental processes.

Changes in chromatin state play important roles in cell fate transitions. Current computational approaches to analyze chromatin modifications across multiple cell types do not model how the cell types are related on a lineage or over time. To overcome this limitation, we developed a method called Chromatin Module INference on Trees (CMINT), a probabilistic clustering approach to systematically capture chromatin state dynamics across multiple cell types. Compared to existing approaches, CMINT can handle complex lineage topologies, capture higher quality clusters, and reliably detect chromatin transitions between cell types. We applied CMINT to gain novel insights in two complex processes: reprogramming to induced pluripotent stem cells (iPSCs) and hematopoiesis. In reprogramming, chromatin changes could occur without large gene expression changes, different combinations of activating marks were associated with specific reprogramming factors, there was an order of acquisition of chromatin marks at pluripotency loci, and multivalent states (comprising previously undetermined combinations of activating and repressive histone modifications) were enriched for CTCF. In the hematopoietic system, we defined critical decision points in the lineage tree, identified regulatory elements that were enriched in cell-type-specific regions, and found that the underlying chromatin state was achieved by specific erasure of preexisting chromatin marks in the precursor cell or by de novo assembly. Our method provides a systematic approach to model the dynamics of chromatin state to provide novel insights into the relationships among cell types in diverse cell-fate specification processes.