Evolution of regulatory networks associated with traits under selection in cichlids

Wiring Between Close Nodes in Molecular Networks Evolves More Quickly Than Between Distant Nodes

As species diverge, a wide range of evolutionary processes lead to changes in protein-protein interaction (PPI) networks and metabolic networks. The rate at which molecular networks evolve is an important question in evolutionary biology. Previous empirical work has focused on interactomes from model organisms to calculate rewiring rates, but this is limited by the relatively small number of species and sparse nature of network data across species. We present a proxy for variation in network topology: variation in drug-drug interactions (DDIs), obtained by studying drug combinations (DCs) across taxa. Here, we propose the rate at which DDIs change across species as an estimate of the rate at which the underlying molecular network changes as species diverge. We computed the evolutionary rates of DDIs using previously published data from a high-throughput study in gram-negative bacteria. Using phylogenetic comparative methods, we found that DDIs diverge rapidly over short evolutionary time periods, but that divergence saturates over longer time periods. In parallel, we mapped drugs with known targets in PPI and cofunctional networks. We found that the targets of synergistic DDIs are closer in these networks than other types of DCs and that synergistic interactions have a higher evolutionary rate, meaning that nodes that are closer evolve at a faster rate. Future studies of network evolution may use DC data to gain larger-scale perspectives on the details of network evolution within and between species.

Hierarchical network structure as the source of hierarchical dynamics (power-law frequency spectra) in living and non-living systems: How state-trait continua (body plans, personalities) emerge from first principles in biophysics

Living systems are hierarchical control systems that display a small world network structure. In such structures, many smaller clusters are nested within fewer larger ones, producing a fractal-like structure with a 'power-law' cluster size distribution (a mereology). Just like their structure, the dynamics of living systems shows fractal-like qualities: the timeseries of inner message passing and overt behavior contain high frequencies or 'states' (treble) that are nested within lower frequencies or 'traits' (bass), producing a power-law frequency spectrum that is known as a 'state-trait continuum' in the behavioral sciences. Here, we argue that the power-law dynamics of living systems results from their power-law network structure: organisms 'vertically encode' the deep spatiotemporal structure of their (anticipated) environments, to the effect that many small clusters near the base of the hierarchy produce high frequency signal changes and fewer larger clusters at its top produce ultra-low frequencies. Such ultra-low frequencies exert a tonic regulatory pressure that produces morphological as well as behavioral traits (i.e., body plans and personalities). Nested-modular structure causes higher frequencies to be embedded within lower frequencies, producing a power-law state-trait continuum. At the heart of such dynamics lies the need for efficient energy dissipation through networks of coupled oscillators, which also governs the dynamics of non-living systems (e.q., earthquakes, stock market fluctuations). Since hierarchical structure produces hierarchical dynamics, the development and collapse of hierarchical structure (e.g., during maturation and disease) should leave specific traces in system dynamics (shifts in lower frequencies, i.e. morphological and behavioral traits) that may serve as early warning signs to system failure. The applications of this idea range from (bio)physics and phylogenesis to ontogenesis and clinical medicine.

Chromatin accessibility in gill tissue identifies candidate genes and loci associated with aquaculture relevant traits in tilapia

The Nile tilapia (Oreochromis niloticus) accounts for ∼9% of global freshwater finfish production however, extreme cold weather and decreasing freshwater resources has created the need to develop resilient strains. By determining the genetic bases of aquaculture relevant traits, we can genotype and breed desirable traits into farmed strains. We generated ATAC-seq and gene expression data from O. niloticus gill tissues, and through the integration of SNPs from 27 tilapia species, identified 1168 highly expressed genes (4% of all Nile tilapia genes) with highly accessible promoter regions with functional variation at transcription factor binding sites (TFBSs). Regulatory variation at these TFBSs is likely driving gene expression differences associated with tilapia gill adaptations, and differentially segregate in freshwater and euryhaline tilapia species. The generation of novel integrative data revealed candidate genes e.g., prolactin receptor 1 and claudin-h, genetic relationships, and loci associated with aquaculture relevant traits like salinity and osmotic stress acclimation.

Molecular Mechanisms Underlying Vertebrate Adaptive Evolution: A Systematic Review

Adaptive evolution is a process in which variation that confers an evolutionary advantage in a specific environmental context arises and is propagated through a population. When investigating this process, researchers have mainly focused on describing advantageous phenotypes or putative advantageous genotypes. A recent increase in molecular data accessibility and technological advances has allowed researchers to go beyond description and to make inferences about the mechanisms underlying adaptive evolution. In this systematic review, we discuss articles from 2016 to 2022 that investigated or reviewed the molecular mechanisms underlying adaptive evolution in vertebrates in response to environmental variation. Regulatory elements within the genome and regulatory proteins involved in either gene expression or cellular pathways have been shown to play key roles in adaptive evolution in response to most of the discussed environmental factors. Gene losses were suggested to be associated with an adaptive response in some contexts. Future adaptive evolution research could benefit from more investigations focused on noncoding regions of the genome, gene regulation mechanisms, and gene losses potentially yielding advantageous phenotypes. Investigating how novel advantageous genotypes are conserved could also contribute to our knowledge of adaptive evolution.

Multilayered Networks of SalmoNet2 Enable Strain Comparisons of the Salmonella Genus on a Molecular Level

Serovars of the genus Salmonella primarily evolved as gastrointestinal pathogens in a wide range of hosts. Some serotypes later evolved further, adopting a more invasive lifestyle in a narrower host range associated with systemic infections. A system-level knowledge of these pathogens could identify the complex adaptations associated with the evolution of serovars with distinct pathogenicity, host range, and risk to human health. This promises to aid the design of interventions and serve as a knowledge base in the Salmonella research community. Here, we present SalmoNet2, a major update to SalmoNet1, the first multilayered interaction resource for Salmonella strains, containing protein-protein, transcriptional regulatory, and enzyme-enzyme interactions. The new version extends the number of Salmonella networks from 11 to 20. We now include a strain from the second species in the Salmonella genus, a strain from the Salmonella enterica subspecies arizonae and additional strains of importance from the subspecies enterica, including S. Typhimurium strain D23580, an epidemic multidrug-resistant strain associated with invasive nontyphoidal salmonellosis (iNTS). The database now uses strain specific metabolic models instead of a generalized model to highlight differences between strains. The update has increased the coverage of high-quality protein-protein interactions, and enhanced interoperability with other computational resources by adopting standardized formats. The resource website has been updated with tutorials to help researchers analyze their Salmonella data using molecular interaction networks from SalmoNet2. SalmoNet2 is accessible at http://salmonet.org/. IMPORTANCE Multilayered network databases collate interaction information from multiple sources, and are powerful both as a knowledge base and subject of analysis. Here, we present SalmoNet2, an integrated network resource containing protein-protein, transcriptional regulatory, and metabolic interactions for 20 Salmonella strains. Key improvements to the update include expanding the number of strains, strain-specific metabolic networks, an increase in high-quality protein-protein interactions, community standard computational formats to help interoperability, and online tutorials to help users analyze their data using SalmoNet2.

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.

High-performance single-cell gene regulatory network inference at scale: the Inferelator 3.0

MOTIVATION

Gene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above.

RESULTS

In this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type-specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data.

AVAILABILITY AND IMPLEMENTATION

The inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Common Themes and Future Challenges in Understanding Gene Regulatory Network Evolution

A major driving force behind the evolution of species-specific traits and novel structures is alterations in gene regulatory networks (GRNs). Comprehending evolution therefore requires an understanding of the nature of changes in GRN structure and the responsible mechanisms. Here, we review two insect pigmentation GRNs in order to examine common themes in GRN evolution and to reveal some of the challenges associated with investigating changes in GRNs across different evolutionary distances at the molecular level. The pigmentation GRN in Drosophila melanogaster and other drosophilids is a well-defined network for which studies from closely related species illuminate the different ways co-option of regulators can occur. The pigmentation GRN for butterflies of the Heliconius species group is less fully detailed but it is emerging as a useful model for exploring important questions about redundancy and modularity in cis-regulatory systems. Both GRNs serve to highlight the ways in which redeployment of trans-acting factors can lead to GRN rewiring and network co-option. To gain insight into GRN evolution, we discuss the importance of defining GRN architecture at multiple levels both within and between species and of utilizing a range of complementary approaches.

The Visual Opsin Gene Repertoires of Teleost Fishes: Evolution, Ecology, and Function

Visual opsin genes expressed in the rod and cone photoreceptor cells of the retina are core components of the visual sensory system of vertebrates. Here, we provide an overview of the dynamic evolution of visual opsin genes in the most species-rich group of vertebrates, teleost fishes. The examination of the rich genomic resources now available for this group reveals that fish genomes contain more copies of visual opsin genes than are present in the genomes of amphibians, reptiles, birds, and mammals. The expansion of opsin genes in fishes is due primarily to a combination of ancestral and lineage-specific gene duplications. Following their duplication, the visual opsin genes of fishes repeatedly diversified at the same key spectral-tuning sites, generating arrays of visual pigments sensitive from the ultraviolet to the red spectrum of the light. Species-specific opsin gene repertoires correlate strongly with underwater light habitats, ecology, and color-based sexual selection. Expected final online publication date for the Annual Review of Cell and Developmental Biology, Volume 37 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.