Regulatory gene networks that shape the development of adaptive phenotypic plasticity in a cichlid fish

Innervation and osteoclast distribution in the inferior pharyngeal jaw of the cichlid Nile tilapia (Oreochromis niloticus)

In addition to an oral jaw, cichlids have a pharyngeal jaw, which is used for crushing and processing captured prey. The teeth and morphology of the pharyngeal jaw bones adapt to changes in prey in response to changes in the growing environment. This study aimed to explore the possible involvement of the peripheral nervous system in remodeling the cichlid pharyngeal jaw by examining the innervation of the inferior pharyngeal jaw in the Nile tilapia, Oreochromis niloticus. Vagal innervation was identified in the Nile tilapia inferior pharyngeal jaw. Double staining with tartrate-resistant acid phosphatase and immunostaining with the neuronal markers, protein gene product 9.5, and acetylated tubulin, revealed that osteoclasts, which play an important role in remodeling, were distributed in the vicinity of the nerves and were in apposition with the nerve terminals. This contact between peripheral nerves and osteoclasts suggests that the peripheral nervous system may play a role in remodeling the inferior pharyngeal jaw in cichlids.

Foraging-induced craniofacial plasticity is associated with an early, robust and dynamic transcriptional response

Phenotypic plasticity is the ability of a single genotype to vary its phenotype in response to the environment. Plasticity of the skeletal system in response to mechanical input is widely studied, but the timing of its transcriptional regulation is not well understood. Here, we used the cichlid feeding apparatus to examine the transcriptional dynamics of skeletal plasticity over time. Using three closely related species that vary in their ability to remodel bone and a panel of 11 genes, including well-studied skeletal differentiation markers and newly characterized environmentally sensitive genes, we examined plasticity at one, two, four and eight weeks following the onset of alternate foraging challenges. We found that the plastic species exhibited environment-specific bursts in gene expression beginning at one week, followed by a sharp decline in levels, while the species with more limited plasticity exhibited consistently low levels of gene expression. This trend held across nearly all genes, suggesting that it is a hallmark of the larger plasticity regulatory network. We conclude that plasticity of the cichlid feeding apparatus is not the result of slowly accumulating gene expression difference over time, but rather is stimulated by early bursts of environment-specific gene expression followed by a return to homeostatic levels.

Genome-Wide Network Analysis of Above- and Below-Ground Co-growth in Populus euphratica

Tree growth is the consequence of developmental interactions between above- and below-ground compartments. However, a comprehensive view of the genetic architecture of growth as a cohesive whole is poorly understood. We propose a systems biology approach for mapping growth trajectories in genome-wide association studies viewing growth as a complex (phenotypic) system in which above- and below-ground components (or traits) interact with each other to mediate systems behavior. We further assume that trait-trait interactions are controlled by a genetic system composed of many different interactive genes and integrate the Lotka-Volterra predator-prey model to dissect phenotypic and genetic systems into pleiotropic and epistatic interaction components by which the detailed genetic mechanism of above- and below-ground co-growth can be charted. We apply the approach to analyze linkage mapping data of Populus euphratica, which is the only tree species that can grow in the desert, and characterize several loci that govern how above- and below-ground growth is cooperated or competed over development. We reconstruct multilayer and multiplex genetic interactome networks for the developmental trajectories of each trait and their developmental covariation. Many significant loci and epistatic effects detected can be annotated to candidate genes for growth and developmental processes. The results from our model may potentially be useful for marker-assisted selection and genetic editing in applied tree breeding programs. The model provides a general tool to characterize a complete picture of pleiotropic and epistatic genetic architecture in growth traits in forest trees and any other organisms.

Growth dynamics and molecular bases of evolutionary novel jaw extensions in halfbeaks and needlefishes (Beloniformes)

Evolutionary novelties-derived traits without clear homology found in the ancestors of a lineage-may promote ecological specialization and facilitate adaptive radiations. Examples for such novelties include the wings of bats, pharyngeal jaws of cichlids and flowers of angiosperms. Belonoid fishes (flying fishes, halfbeaks and needlefishes) feature an astonishing diversity of extremely elongated jaw phenotypes with undetermined evolutionary origins. We investigate the development of elongated jaws in a halfbeak (Dermogenys pusilla) and a needlefish (Xenentodon cancila) using morphometrics, transcriptomics and in situ hybridization. We confirm that these fishes' elongated jaws are composed of distinct base and novel 'extension' portions. These extensions are morphologically unique to belonoids, and we describe the growth dynamics of both bases and extensions throughout early development in both studied species. From transcriptomic profiling, we deduce that jaw extension outgrowth is guided by populations of multipotent cells originating from the anterior tip of the dentary. These cells are shielded from differentiation, but proliferate and migrate anteriorly during the extension's allometric growth phase. Cells left behind at the tip leave the shielded zone and undergo differentiation into osteoblast-like cells, which deposit extracellular matrix with both bone and cartilage characteristics that mineralizes and thereby provides rigidity. Such bone has characteristics akin to histological observations on the elongated 'kype' process on lower jaws of male salmon, which may hint at common conserved regulatory underpinnings. Future studies will evaluate the molecular pathways that govern the anterior migration and proliferation of these multipotent cells underlying the belonoids' evolutionary novel jaw extensions.

Speciation and repeated origins of hypertrophied lips in parallel adaptive radiations of cyprinid fish from East Africa

The evolution of convergent phenotypes is one of the most interesting phenomena of repeated adaptive radiations. Here, we examined the repeated patterns of thick-lipped or "rubberlip" phenotype of cyprinid fish of the genus Labeobarbus discovered in riverine environments of the Ethiopian Highlands, East Africa. To test the adaptive value of thickened lips, identify the ecological niche of the thick-lipped ecomorphs, and test whether these ecomorphs are the products of adaptive divergence, we studied six sympatric pairs of ecomorphs with hypertrophied lips and the normal lip structure from different riverine basins. Trophic morphology, diet, stable isotope (δ15N and δ13C) signatures, as well as mtDNA markers and genome-wide SNP variation, were analyzed. Our results show that thick-lipped ecomorphs partition trophic resources with generalized ecomorphs in only one-half of the examined sympatric pairs despite the pronounced divergence in lip structure. In these thick-lipped ecomorphs that were trophically diverged, the data on their diet along with the elevated 15N values suggest an invertivorous specialization different from the basal omnivorous-detritivouros feeding mode of the generalized ecomorphs. Genetic data confirmed an independent and parallel origin of all six lipped ecomorphs. Yet, only one of those six thick-lipped ecomorphs had a notable genetic divergence with sympatric non-lipped ecomorphs based on nuclear SNPs data (F ST = 0.21). Sympatric pairs can be sorted by combinations of phenotypic, ecological, and genetic divergence from an ecologically non-functional mouth polymorphism via ecologically functional polymorphism to a matured speciation stage via divergent evolution.

The transcriptional state and chromatin landscape of cichlid jaw shape variation across species and environments

Adaptive phenotypes are shaped by a combination of genetic and environmental forces, but how they interact remains poorly understood. Here, we utilize the cichlid oral jaw apparatus to better understand these gene-by-environment effects. First, we employed RNA-seq in bony and ligamentous tissues important for jaw opening to identify differentially expressed genes between species and across foraging environments. We used two Lake Malawi species adapted to different foraging habitats along the pelagic-benthic ecomorphological axis. Our foraging treatments were designed to force animals to employ either suction or biting/scraping, which broadly mimic pelagic or benthic modes of feeding. We found a large number of differentially expressed genes between species, and while we identified relatively few differences between environments, species differences were far more pronounced when they were challenged with a pelagic versus benthic foraging mode. Expression data carried the signature of genetic assimilation, and implicated cell cycle regulation in shaping the jaw across species and environments. Next, we repeated the foraging experiment and performed ATAC-seq procedures on nuclei harvested from the same tissues. Cross-referencing results from both analyses revealed subsets of genes that were both differentially expressed and differentially accessible. This reduced dataset implicated notable candidate genes including the Hedgehog effector, KIAA0586 and the ETS transcription factor, etv4, which connects environmental stress and craniofacial morphogenesis. Taken together, these data provide novel insights into the epigenetic, genetic and cellular bases of species- and environment-specific bone shapes.

East African cichlid fishes

Cichlid fishes are a very diverse and species-rich family of teleost fishes that inhabit lakes and rivers of India, Africa, and South and Central America. Research has largely focused on East African cichlids of the Rift Lakes Tanganyika, Malawi, and Victoria that constitute the biodiversity hotspots of cichlid fishes. Here, we give an overview of the study system, research questions, and methodologies. Research on cichlid fishes spans many disciplines including ecology, evolution, physiology, genetics, development, and behavioral biology. In this review, we focus on a range of organismal traits, including coloration phenotypes, trophic adaptations, appendages like fins and scales, sensory systems, sex, brains, and behaviors. Moreover, we discuss studies on cichlid phylogenies, plasticity, and general evolutionary patterns, ranging from convergence to speciation rates and the proximate and ultimate mechanisms underlying these processes. From a methodological viewpoint, the last decade has brought great advances in cichlid fish research, particularly through the advent of affordable deep sequencing and advances in genetic manipulations. The ability to integrate across traits and research disciplines, ranging from developmental biology to ecology and evolution, makes cichlid fishes a fascinating research system.

Computational modelling of plasticity-led evolution

Plasticity-led evolution is a form of evolution where a change in the environment induces novel traits via phenotypic plasticity, after which the novel traits are genetically accommodated over generations under the novel environment. This mode of evolution is expected to resolve the problem of gradualism (i.e., evolution by the slow accumulation of mutations that induce phenotypic variation) implied by the Modern Evolutionary Synthesis, in the face of a large environmental change. While experimental works are essential for validating that plasticity-led evolution indeed happened, we need computational models to gain insight into its underlying mechanisms and make qualitative predictions. Such computational models should include the developmental process and gene-environment interactions in addition to genetics and natural selection. We point out that gene regulatory network models can incorporate all the above notions. In this review, we highlight results from computational modelling of gene regulatory networks that consolidate the criteria of plasticity-led evolution. Since gene regulatory networks are mathematically equivalent to artificial recurrent neural networks, we also discuss their analogies and discrepancies, which may help further understand the mechanisms underlying plasticity-led evolution.

Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

Metabolic response provides insights into the mechanism of adaption to hypoxia in largemouth bass (Micropterus salmoides) under intermittent hypoxic conditions

In metabolism, molecular oxygen is a necessary substrate. Oxygen imbalances are linked to a variety of circumstances in the organism's homeostasis. Recently, the positive effects of hypoxia treatment in improving exercise ability and hypoxia tolerance have become a research focus. We explored the effects of intermittent hypoxia exposure (IHE, for one hour or three hours per day) on the hypoxia tolerance of largemouth bass in this study. The results showed that (1) IHE significantly reduced the LOEcrit (the critical O₂ tension for loss of equilibrium) value of largemouth bass, indicating that its hypoxia tolerance was enhanced. (2) The level of oxidative stress in the liver decreased in the HH3 group (exposed to a hypoxic condition for 3 h per day) compared to HH1 group (exposed to a hypoxic condition for 1 h per day). (3) IHE reduced the content of lactic acid and enhanced the process of gluconeogenesis in the liver. (4) Importantly, lipid mobilization and fatty acid oxidation in the liver of largemouth bass were significantly enhanced during IHE. In short, the results of this study indicate that IHE can improve hypoxia tolerance by regulating the energy metabolism of largemouth bass.

Evolved transcriptional responses and their trade-offs after long-term adaptation of Bemisia tabaci to a marginally-suitable host

Although generalist insect herbivores can migrate and rapidly adapt to a broad range of host plants, they can face significant difficulties when accidentally migrating to novel and marginally-suitable hosts. What happens, both in performance and gene expression regulation, if these marginally-suitable hosts must be used for multiple generations before migration to a suitable host can take place, largely remains unknown. In this study, we established multigenerational colonies of the whitefly Bemisia tabaci, a generalist phloem-feeding species, adapted to a marginally-suitable host (habanero pepper) or an optimal host (cotton). We used reciprocal host tests to estimate the differences in performance of the populations on both hosts under optimal (30 oC) and mild-stressful (24 oC) temperature conditions, and documented the associated transcriptomic changes. The habanero pepper-adapted population greatly improved its performance on habanero pepper but did not reach its performance level on cotton, the original host. It also showed reduced performance on cotton, relative to the non-adapted population, and an antagonistic effect of the lower-temperature stressor. The transcriptomic data revealed that most of the expression changes, associated with long-term adaptation to habanero pepper, can be categorized as "evolved" with no initial plastic response. Three molecular functions dominated: enhanced formation of cuticle structural constituents, enhanced activity of oxidation-reduction processes involved in neutralization of phytotoxins and reduced production of proteins from the cathepsin B family. Taken together, these findings indicate that generalist insects can adapt to novel host plants by modifying the expression of a relatively small set of specific molecular functions.

Phylogenomics of trophically diverse cichlids disentangles processes driving adaptive radiation and repeated trophic transitions

Cichlid fishes of the tribe Tropheini are a striking case of adaptive radiation, exemplifying multiple trophic transitions between herbivory and carnivory occurring in sympatry with other established cichlid lineages. Tropheini evolved highly specialized eco-morphologies to exploit similar trophic niches in different ways repeatedly and rapidly. To better understand the evolutionary history and trophic adaptations of this lineage, we generated a dataset of 532 targeted loci from 21 out of the 22 described Tropheini species. We resolved the Tropheini into seven monophyletic genera and discovered one to be polyphyletic. The polyphyletic genus, Petrochromis, represents three convergent origins of the algae grazing trophic specialization. This repeated evolution of grazing may have been facilitated by adaptive introgression as we found evidence for gene flow among algae grazing genera. We also found evidence of gene flow among algae browsing genera, but gene flow was restricted between herbivorous and carnivorous genera. Furthermore, we observed no evidence supporting a hybrid origin of this radiation. Our molecular evolutionary analyses suggest that opsin genes likely evolved in response to selection pressures associated with trophic ecology in the Tropheini. We found surprisingly little evidence of positive selection in coding regions of jaw-shaping genes in this trophically diverse lineage. This suggests low degrees of freedom for further change in these genes, and possibly a larger role for regulatory variation in driving jaw adaptations. Our study emphasizes Tropheini cichlids as an important model for studying the evolution of trophic specialization and its role in speciation.

My road to the ants: A model clade for eco-evo-devo

This chapter is the story of how I pioneered ants as a system for studying eco-evo-devo, a field that integrates developmental biology with ecology and evolutionary biology. One aim of eco-evo-devo is to understand how the interactions between genes and their environments during development facilitates the origin and evolution of novel phenotypes. In a series of six parts, I review some of the key discoveries from my lab on how novel worker caste systems in ants--soldiers and supersoldiers--originated and evolved. I also discuss some of the ideas that emerged from these discoveries, including the role that polyphenisms, hidden developmental potentials, and rudimentary organs play in facilitating developmental and evolutionary change. As superorganisms, I argue that ants are uniquely positioned to reveal types of variation that are often difficult to observe in nature. In doing so, they have the potential to transform our view of biology and provide new perspectives in medicine, agriculture, and biodiversity conservation. With my story I hope to inspire the next generation of biologists to continue exploring the unknown regions of phenotypic space to solve some of our most pressing societal challenges.

Gene by Environment Interactions reveal new regulatory aspects of signaling network plasticity

Phenotypes can change during exposure to different environments through the regulation of signaling pathways that operate in integrated networks. How signaling networks produce different phenotypes in different settings is not fully understood. Here, Gene by Environment Interactions (GEIs) were used to explore the regulatory network that controls filamentous/invasive growth in the yeast Saccharomyces cerevisiae. GEI analysis revealed that the regulation of invasive growth is decentralized and varies extensively across environments. Different regulatory pathways were critical or dispensable depending on the environment, microenvironment, or time point tested, and the pathway that made the strongest contribution changed depending on the environment. Some regulators even showed conditional role reversals. Ranking pathways' roles across environments revealed an under-appreciated pathway (OPI1) as the single strongest regulator among the major pathways tested (RAS, RIM101, and MAPK). One mechanism that may explain the high degree of regulatory plasticity observed was conditional pathway interactions, such as conditional redundancy and conditional cross-pathway regulation. Another mechanism was that different pathways conditionally and differentially regulated gene expression, such as target genes that control separate cell adhesion mechanisms (FLO11 and SFG1). An exception to decentralized regulation of invasive growth was that morphogenetic changes (cell elongation and budding pattern) were primarily regulated by one pathway (MAPK). GEI analysis also uncovered a round-cell invasion phenotype. Our work suggests that GEI analysis is a simple and powerful approach to define the regulatory basis of complex phenotypes and may be applicable to many systems.

Inferring multilayer interactome networks shaping phenotypic plasticity and evolution

Phenotypic plasticity represents a capacity by which the organism changes its phenotypes in response to environmental stimuli. Despite its pivotal role in adaptive evolution, how phenotypic plasticity is genetically controlled remains elusive. Here, we develop a unified framework for coalescing all single nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS) into a quantitative graph. This framework integrates functional genetic mapping, evolutionary game theory, and predator-prey theory to decompose the net genetic effect of each SNP into its independent and dependent components. The independent effect arises from the intrinsic capacity of a SNP, only expressed when it is in isolation, whereas the dependent effect results from the extrinsic influence of other SNPs. The dependent effect is conceptually beyond the traditional definition of epistasis by not only characterizing the strength of epistasis but also capturing the bi-causality of epistasis and the sign of the causality. We implement functional clustering and variable selection to infer multilayer, sparse, and multiplex interactome networks from any dimension of genetic data. We design and conduct two GWAS experiments using Staphylococcus aureus, aimed to test the genetic mechanisms underlying the phenotypic plasticity of this species to vancomycin exposure and Escherichia coli coexistence. We reconstruct the two most comprehensive genetic networks for abiotic and biotic phenotypic plasticity. Pathway analysis shows that SNP-SNP epistasis for phenotypic plasticity can be annotated to protein-protein interactions through coding genes. Our model can unveil the regulatory mechanisms of significant loci and excavate missing heritability from some insignificant loci. Our multilayer genetic networks provide a systems tool for dissecting environment-induced evolution.

Polyphenism of a Novel Trait Integrated Rapidly Evolving Genes into Ancestrally Plastic Networks

Developmental polyphenism, the ability to switch between phenotypes in response to environmental variation, involves the alternating activation of environmentally sensitive genes. Consequently, to understand how a polyphenic response evolves requires a comparative analysis of the components that make up environmentally sensitive networks. Here, we inferred coexpression networks for a morphological polyphenism, the feeding-structure dimorphism of the nematode Pristionchus pacificus. In this species, individuals produce alternative forms of a novel trait—moveable teeth, which in one morph enable predatory feeding—in response to environmental cues. To identify the origins of polyphenism network components, we independently inferred coexpression modules for more conserved transcriptional responses, including in an ancestrally nonpolyphenic nematode species. Further, through genome-wide analyses of these components across the nematode family (Diplogastridae) in which the polyphenism arose, we reconstructed how network components have changed. To achieve this, we assembled and resolved the phylogenetic context for five genomes of species representing the breadth of Diplogastridae and a hypothesized outgroup. We found that gene networks instructing alternative forms arose from ancestral plastic responses to environment, specifically starvation-induced metabolism and the formation of a conserved diapause (dauer) stage. Moreover, loci from rapidly evolving gene families were integrated into these networks with higher connectivity than throughout the rest of the P. pacificus transcriptome. In summary, we show that the modular regulatory outputs of a polyphenic response evolved through the integration of conserved plastic responses into networks with genes of high evolutionary turnover.

Convergent Evolution of Cichlid Fish Pharyngeal Jaw Dentitions in Mollusk-Crushing Predators: Comparative X-Ray Computed Tomography of Tooth Sizes, Numbers, and Replacement

Dental convergence is a hallmark of cichlid fish adaptive radiations. This type of repeated evolution characterizes both the oral jaws of these fishes as well as their pharyngeal jaws that are modified gill arches used to functionally process prey like hard-shelled mollusks. To test several hypotheses regarding the evolution of cichlid crushing pharyngeal dentitions, we used X-ray computed tomography scans to comparatively examine dental evolution in the pharyngeal jaw of a diversity of New World Heroine cichlid lineages. The substantial variation in erupted tooth sizes and numbers as well as replacement teeth found in these fishes showed several general patterns. Larger toothed species tended to have fewer teeth suggesting a potential role of spatial constraints in cichlid dental divergence. Species with larger numbers of erupted pharyngeal teeth also had larger numbers of replacement teeth. Replacement tooth size is almost exactly predicted (r = 0.99) from the size of erupted teeth across all of the species. Mollusk crushing was, therefore, highly associated with not only larger pharyngeal teeth, but also larger replacement teeth. Whether dental divergence arises as a result of environmental induced plasticity or originates via trophic polymorphism as found in the species Herichthys minckleyi, there appear to be general rules that structure interspecific divergence in cichlid pharyngeal erupted and replacement dentitions.

Phenotypic Plasticity in Vertebrate Dentitions

Vertebrates interact directly with food items through their dentition, and these interactions with trophic resources could often feedback to influence tooth structure. Although dentitions are often considered to be a fixed phenotype, there is the potential for environmentally induced phenotypic plasticity in teeth to extensively influence their diversity. Here, we review the literature concerning phenotypic plasticity of vertebrate teeth. Even though only a few taxonomically disparate studies have focused on phenotypic plasticity in teeth, there are a number of ways teeth can change their size, shape, or patterns of replacement as a response to the environment. Elucidating the underlying physiological, developmental, and genetic mechanisms that generate phenotypic plasticity can clarify its potential role in the evolution of dental phenotypes.

Gene coexpression networks reveal molecular interactions underlying cichlid jaw modularity

Background

The oral and pharyngeal jaw of cichlid fishes are a classic example of evolutionary modularity as their functional decoupling boosted trophic diversification and contributed to the success of cichlid adaptive radiations. Most studies until now have focused on the functional, morphological, or genetic aspects of cichlid jaw modularity. Here we extend this concept to include transcriptional modularity by sequencing whole transcriptomes of the two jaws and comparing their gene coexpression networks.

Results

We show that transcriptional decoupling of gene expression underlies the functional decoupling of cichlid oral and pharyngeal jaw apparatus and the two units are evolving independently in recently diverged cichlid species from Lake Tanganyika. Oral and pharyngeal jaw coexpression networks reflect the common origin of the jaw regulatory program as there is high preservation of gene coexpression modules between the two sets of jaws. However, there is substantial rewiring of genetic architecture within those modules. We define a global jaw coexpression network and highlight jaw-specific and species-specific modules within it. Furthermore, we annotate a comprehensive in silico gene regulatory network linking the Wnt and AHR signalling pathways to jaw morphogenesis and response to environmental cues, respectively. Components of these pathways are significantly differentially expressed between the oral and pharyngeal jaw apparatus.

Conclusion

This study describes the concerted expression of many genes in cichlid oral and pharyngeal jaw apparatus at the onset of the independent life of cichlid fishes. Our findings suggest that – on the basis of an ancestral gill arch network—transcriptional rewiring may have driven the modular evolution of the oral and pharyngeal jaws, highlighting the evolutionary significance of gene network reuse. The gene coexpression and in silico regulatory networks presented here are intended as resource for future studies on the genetics of vertebrate jaw morphogenesis and trophic adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01787-9.

A Genomic Cluster Containing Novel and Conserved Genes is Associated with Cichlid Fish Dental Developmental Convergence

The two toothed jaws of cichlid fishes provide textbook examples of convergent evolution. Tooth phenotypes such as enlarged molar-like teeth used to process hard-shelled mollusks have evolved numerous times independently during cichlid diversification. Although the ecological benefit of molar-like teeth to crush prey is known, it is unclear whether the same molecular mechanisms underlie these convergent traits. To identify genes involved in the evolution and development of enlarged cichlid teeth, we performed RNA-seq on the serially homologous-toothed oral and pharyngeal jaws as well as the fourth toothless gill arch of Astatoreochromis alluaudi. We identified 27 genes that are highly upregulated on both tooth-bearing jaws compared with the toothless gill arch. Most of these genes have never been reported to play a role in tooth formation. Two of these genes (unk, rpfA) are not found in other vertebrate genomes but are present in all cichlid genomes. They also cluster genomically with two other highly expressed tooth genes (odam, scpp5) that exhibit conserved expression during vertebrate odontogenesis. Unk and rpfA were confirmed via in situ hybridization to be expressed in developing teeth of Astatotilapia burtoni. We then examined expression of the cluster's four genes in six evolutionarily independent and phylogenetically disparate cichlid species pairs each with a large- and a small-toothed species. Odam and unk commonly and scpp5 and rpfA always showed higher expression in larger toothed cichlid jaws. Convergent trophic adaptations across cichlid diversity are associated with the repeated developmental deployment of this genomic cluster containing conserved and novel cichlid-specific genes.

Cross-species transcriptomics uncovers genes underlying genetic accommodation of developmental plasticity in spadefoot toads

That hardcoded genomes can manifest as plastic phenotypes responding to environmental perturbations is a fascinating feature of living organisms. How such developmental plasticity is regulated at the molecular level is beginning to be uncovered aided by the development of -omic techniques. Here, we compare the transcriptome-wide responses of two species of spadefoot toads with differing capacity for developmental acceleration of their larvae in the face of a shared environmental risk: pond drying. By comparing gene expression profiles over time and performing cross-species network analyses, we identified orthologues and functional gene pathways whose environmental sensitivity in expression have diverged between species. Genes related to lipid, cholesterol and steroid biosynthesis and metabolism make up most of a module of genes environmentally responsive in one species, but canalized in the other. The evolutionary changes in the regulation of the genes identified through these analyses may have been key in the genetic accommodation of developmental plasticity in this system.

Rearing environment affects the genetic architecture and plasticity of DNA methylation in Chinook salmon

Genetic architecture and phenotypic plasticity are important considerations when studying trait variation within and among populations. Since environmental change can induce shifts in the genetic architecture and plasticity of traits, it is important to consider both genetic and environmental sources of phenotypic variation. While there is overwhelming evidence for environmental effects on phenotype, the underlying mechanisms are less clear. Variation in DNA methylation is a potential mechanism mediating environmental effects on phenotype due to its sensitivity to environmental stimuli, transgenerational inheritance, and influences on transcription. To characterize the effect of environment on methylation, we created two 6 × 6 (North Carolina II) Chinook salmon breeding crosses and reared the offspring in two environments: uniform hatchery tanks and seminatural stream channels. We sampled the fish twice during development, at the alevin (larval) and fry (juvenile) stages. We measured DNA methylation at 13 genes using a PCR-based bisulfite sequencing protocol. The genetic architecture of DNA methylation differed between rearing environments, with greater additive and nonadditive genetic variance in hatchery fish and greater maternal effects in seminatural channel fish, though gene-specific variation was evident. We observed plasticity in methylation across all assayed genes, as well as gene-specific effects at two genes in alevin and six genes in fry, indicating developmental stage-specific effects of rearing environment on methylation. Characterizing genetic and environmental influences on methylation is critical for future studies on DNA methylation as a potential mechanism for acclimation and adaptation.

Climate explains population divergence in drought-induced plasticity of functional traits and gene expression in a South African Protea

Long-term environmental variation often drives local adaptation and leads to trait differentiation across populations. Additionally, when traits change in an environment-dependent way through phenotypic plasticity, the genetic variation underlying plasticity will also be under selection. These processes could create a landscape of differentiation across populations in traits and their plasticity. Here, we performed a dry-down experiment under controlled conditions to measure responses in seedlings of a shrub species from the Cape Floristic Region, the common sugarbush (Protea repens). We measured morphological and physiological traits, and sequenced whole transcriptomes of leaf tissues from eight populations that represent both the climatic and the geographical distribution of this species. We found that there is substantial variation in how populations respond to drought, but we also observed common patterns such as reduced leaf size and leaf thickness, and up-regulation of stress-related and down-regulation of growth-related gene groups. Both high environmental heterogeneity and milder source site climates were associated with higher plasticity in various traits and co-expression gene networks. Associations between traits, trait plasticity, gene networks and the source site climate suggest that temperature may play a greater role in shaping these patterns when compared to precipitation, in line with recent changes in the region due to climate change. We also found that traits respond to climatic variation in an environment-dependent manner: some associations between traits and climate were apparent only under certain growing conditions. Together, our results uncover common responses of P. repens populations to drought, and climatic drivers of population differentiation in functional traits, gene expression and their plasticity.

SEGN: Inferring real-time gene networks mediating phenotypic plasticity

The capacity of an organism to alter its phenotype in response to environmental perturbations changes over developmental time and is a process determined by multiple genes that are co-expressed in intricate but organized networks. Characterizing the spatiotemporal change of such gene networks can offer insight into the genomic signatures underlying organismic adaptation, but it represents a major methodological challenge. Here, we integrate the holistic view of systems biology and the interactive notion of evolutionary game theory to reconstruct so-called systems evolutionary game networks (SEGN) that can autonomously detect, track, and visualize environment-induced gene networks along the time axis. The SEGN overcomes the limitations of traditional approaches by inferring context-specific networks, encapsulating bidirectional, signed, and weighted gene-gene interactions into fully informative networks, and monitoring the process of how networks topologically alter across environmental and developmental cues. Based on the design principle of SEGN, we perform a transcriptional plasticity study by culturing Euphrates poplar, a tree that can grow in the saline desert, in saline-free and saline-stress conditions. SEGN characterize previously unknown gene co-regulation that modulates the time trajectories of the trees' response to salt stress. As a marriage of multiple disciplines, SEGN shows its potential to interpret gene interdependence, predict how transcriptional co-regulation responds to various regimes, and provides a hint for exploring the mass, energetic, or signal basis that drives various types of gene interactions.

[Phenotypic plasticity: a brief introduction]

Phenotypic plasticity describes the ability of a given genotype to produce different phenotypes in response to distinct environmental conditions. It has major implications in agronomy, animal husbandry and medicine and is also thought to facilitate evolution. Phenotypic plasticity is widely observed in the wild. It is only relatively recently that the mechanisms involved in phenotypic plasticity have been analysed. Thanks to laboratory experiments we understand better how environmental conditions are involved in phenotypic variations. This article introduces major concepts from the phenotypic plasticity field, presents briefly mechanisms involved in phenotypic plasticity and discusses the links between phenotypic plasticity and evolution.

Hedgehog signaling is necessary and sufficient to mediate craniofacial plasticity in teleosts

Phenotypic plasticity has emerged as an important concept in evolutionary biology. It is thought to contribute to an organism’s ability to adapt to environmental change within a single generation, which may facilitate survival and increase fitness. Furthermore, plasticity has the potential to bias the direction and/or speed of evolution by changing patterns of phenotypic variation and exposing new genetic variation to selection (i.e., flexible stem evolution). Our understanding of this important phenomenon is incomplete owing to limited knowledge of the molecular underpinnings of reaction norm evolution. Using the teleost feeding apparatus as a model, we explore this open question and show that the Hh signaling pathway underlies the ability of this structure to respond plastically to alternate feeding regimes.

Phenotypic plasticity, the ability of a single genotype to produce multiple phenotypes under different environmental conditions, is critical for the origins and maintenance of biodiversity; however, the genetic mechanisms underlying plasticity as well as how variation in those mechanisms can drive evolutionary change remain poorly understood. Here, we examine the cichlid feeding apparatus, an icon of both prodigious evolutionary divergence and adaptive phenotypic plasticity. We first provide a tissue-level mechanism for plasticity in craniofacial shape by measuring rates of bone deposition within functionally salient elements of the feeding apparatus in fishes forced to employ alternate foraging modes. We show that levels and patterns of phenotypic plasticity are distinct among closely related cichlid species, underscoring the evolutionary potential of this trait. Next, we demonstrate that hedgehog (Hh) signaling, which has been implicated in the evolutionary divergence of cichlid feeding architecture, is associated with environmentally induced rates of bone deposition. Finally, to demonstrate that Hh levels are the cause of the plastic response and not simply the consequence of producing more bone, we use transgenic zebrafish in which Hh levels could be experimentally manipulated under different foraging conditions. Notably, we find that the ability to modulate bone deposition rates in different environments is dampened when Hh levels are reduced, whereas the sensitivity of bone deposition to different mechanical demands increases with elevated Hh levels. These data advance a mechanistic understanding of phenotypic plasticity in the teleost feeding apparatus and in doing so contribute key insights into the origins of adaptive morphological radiations.

Appetite regulating genes may contribute to herbivory versus carnivory trophic divergence in haplochromine cichlids

Feeding is a complex behaviour comprised of satiety control, foraging, ingestion and subsequent digestion. Cichlids from the East African Great Lakes are renowned for their diverse trophic specializations, largely predicated on highly variable jaw morphologies. Thus, most research has focused on dissecting the genetic, morphological and regulatory basis of jaw and teeth development in these species. Here for the first time we explore another aspect of feeding, the regulation of appetite related genes that are expressed in the brain and control satiety in cichlid fishes. Using qPCR analysis, we first validate stably expressed reference genes in the brain of six haplochromine cichlid species at the end of larval development prior to foraging. We next evaluate the expression of 16 appetite related genes in herbivorous and carnivorous species from the parallel radiations of Lake Tanganyika, Malawi and Victoria. Interestingly, we find increased expression of two appetite-regulating genes (anorexigenic genes), cart and npy2r, in the brain of carnivorous species in all the three lakes. This supports the notion that appetite gene regulation might play a part in determining trophic niche specialization in divergent cichlid species, already prior to exposure to different diets. Our study contributes to the limited body of knowledge on the neurological circuitry that controls feeding transitions and adaptations in cichlids and other teleosts.

Divergent Allometric Trajectories in Gene Expression and Coexpression Produce Species Differences in Sympatrically Speciating Midas Cichlid Fish

The mechanisms of speciation without geographic isolation (i.e., sympatric speciation) remain debated. This is due in part to the fact that the genomic landscape that could promote or hinder species divergence in the presence of gene flow is still largely unknown. However, intensive research is now centered on understanding the genetic architecture of adaptive traits associated with this process as well as how gene expression might affect these traits. Here, using RNA-Seq data, we investigated gene expression of sympatrically speciating benthic and limnetic Neotropical cichlid fishes at two developmental stages. First, we identified groups of coexpressed genes (modules) at each stage. Although there are a few large and well-preserved modules, most of the other modules are not preserved across life stages. Second, we show that later in development more and larger coexpression modules are associated with divergence between benthic and limnetic fish compared with the earlier life stage. This divergence between benthic and limnetic fish in coexpression mirrors divergence in overall expression between benthic and limnetic fish, which is more pronounced later in life. Our results reveal that already at 1-day posthatch benthic and limnetic fish diverge in (co)expression, and that this divergence becomes more substantial when fish are free-swimming but still unlikely to have divergent swimming and feeding habits. More importantly, our study describes how the coexpression of several genes through development, as opposed to individual genes, is associated with benthic–limnetic species differences, and how two morphogenetic trajectories diverge as fish grow older.

Conservation and novelty in the microRNA genomic landscape of hyperdiverse cichlid fishes

MicroRNAs (miRNAs) play crucial roles in the post-transcriptional control of messenger RNA (mRNA). These miRNA-mRNA regulatory networks are present in nearly all organisms and contribute to development, phenotypic divergence, and speciation. To examine the miRNA landscape of cichlid fishes, one of the most species-rich families of vertebrates, we profiled the expression of both miRNA and mRNA in a diverse set of cichlid lineages. Among these, we found that conserved miRNAs differ from recently arisen miRNAs (i.e. lineage specific) in average expression levels, number of target sites, sequence variability, and physical clustering patterns in the genome. Furthermore, conserved miRNA target sites tend to be enriched at the 5' end of protein-coding gene 3' UTRs. Consistent with the presumed regulatory role of miRNAs, we detected more negative correlations between the expression of miRNA-mRNA functional pairs than in random pairings. Finally, we provide evidence that novel miRNA targets sites are enriched in genes involved in protein synthesis pathways. Our results show how conserved and evolutionarily novel miRNAs differ in their contribution to the genomic landscape and highlight their particular evolutionary roles in the adaptive diversification of cichlids.

Divergence in larval jaw gene expression reflects differential trophic adaptation in haplochromine cichlids prior to foraging

Background

Understanding how variation in gene expression contributes to morphological diversity is a major goal in evolutionary biology. Cichlid fishes from the East African Great lakes exhibit striking diversity in trophic adaptations predicated on the functional modularity of their two sets of jaws (oral and pharyngeal). However, the transcriptional basis of this modularity is not so well understood, as no studies thus far have directly compared the expression of genes in the oral and pharyngeal jaws. Nor is it well understood how gene expression may have contributed to the parallel evolution of trophic morphologies across the replicate cichlid adaptive radiations in Lake Tanganyika, Malawi and Victoria.

Results

We set out to investigate the role of gene expression divergence in cichlid fishes from these three lakes adapted to herbivorous and carnivorous trophic niches. We focused on the development stage prior to the onset of exogenous feeding that is critical for understanding patterns of gene expression after oral and pharyngeal jaw skeletogenesis, anticipating environmental cues. This framework permitted us for the first time to test for signatures of gene expression underlying jaw modularity in convergent eco-morphologies across three independent adaptive radiations. We validated a set of reference genes, with stable expression between the two jaw types and across species, which can be important for future studies of gene expression in cichlid jaws. Next we found evidence of modular and non-modular gene expression between the two jaws, across different trophic niches and lakes. For instance, prdm1a, a skeletogenic gene with modular anterior-posterior expression, displayed higher pharyngeal jaw expression and modular expression pattern only in carnivorous species. Furthermore, we found the expression of genes in cichlids jaws from the youngest Lake Victoria to exhibit low modularity compared to the older lakes.

Conclusion

Overall, our results provide cross-species transcriptional comparisons of modularly-regulated skeletogenic genes in the two jaw types, implicating expression differences which might contribute to the formation of divergent trophic morphologies at the stage of larval independence prior to foraging.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1483-3) contains supplementary material, which is available to authorized users.

Gene expression in the phenotypically plastic Arctic charr (Salvelinus alpinus): A focus on growth and ossification at early stages of development

Gene expression during development shapes the phenotypes of individuals. Although embryonic gene expression can have lasting effects on developmental trajectories, few studies consider the role of maternal effects, such as egg size, on gene expression. Using qPCR, we characterize relative expression of 14 growth and/or skeletal promoting genes across embryonic development in Arctic charr (Salvelinus alpinus). We test to what extent their relative expression is correlated with egg size and size at early life‐stages within the study population. We predict smaller individuals to have higher expression of growth and skeletal promoting genes, due to less maternal resources (i.e., yolk) and prioritization of energy toward ossification. We found expression levels to vary across developmental stages and only three genes (Mmp9, Star, and Sgk1) correlated with individual size at a given developmental stage. Contrary to our hypothesis, expression of Mmp9 and Star showed a non‐linear relationship with size (at post fertilization and hatching, respectively), whilst Sgk1 was higher in larger embryos at hatching. Interestingly, these genes are also associated with craniofacial divergence of Arctic charr morphs. Our results indicate that early life‐stage variation in gene expression, concomitant to maternal effects, can influence developmental plasticity and potentially the evolution of resource polymorphism in fishes.

Phylogenomics uncovers early hybridization and adaptive loci shaping the radiation of Lake Tanganyika cichlid fishes

Lake Tanganyika is the oldest and phenotypically most diverse of the three East African cichlid fish adaptive radiations. It is also the cradle for the younger parallel haplochromine cichlid radiations in Lakes Malawi and Victoria. Despite its evolutionary significance, the relationships among the main Lake Tanganyika lineages remained unresolved, as did the general timescale of cichlid evolution. Here, we disentangle the deep phylogenetic structure of the Lake Tanganyika radiation using anchored phylogenomics and uncover hybridization at its base, as well as early in the haplochromine radiation. This suggests that hybridization might have facilitated these speciation bursts. Time-calibrated trees support that the radiation of Tanganyika cichlids coincided with lake formation and that Gondwanan vicariance concurred with the earliest splits in the cichlid family tree. Genes linked to key innovations show signals of introgression or positive selection following colonization of lake habitats and species' dietary adaptations are revealed as major drivers of colour vision evolution. These findings shed light onto the processes shaping the evolution of adaptive radiations.

bric à brac (bab), a central player in the gene regulatory network that mediates thermal plasticity of pigmentation in Drosophila melanogaster

Drosophila body pigmentation has emerged as a major Evo-Devo model. Using two Drosophila melanogaster lines, Dark and Pale, selected from a natural population, we analyse here the interaction between genetic variation and environmental factors to produce this complex trait. Indeed, pigmentation varies with genotype in natural populations and is sensitive to temperature during development. We demonstrate that the bric à brac (bab) genes, that are differentially expressed between the two lines and whose expression levels vary with temperature, participate in the pigmentation difference between the Dark and Pale lines. The two lines differ in a bab regulatory sequence, the dimorphic element (called here bDE). Both bDE alleles are temperature-sensitive, but the activity of the bDE allele from the Dark line is lower than that of the bDE allele from the Pale line. Our results suggest that this difference could partly be due to differential regulation by AbdB. bab has been previously reported to be a repressor of abdominal pigmentation. We show here that one of its targets in this process is the pigmentation gene tan (t), regulated via the tan abdominal enhancer (t_MSE). Furthermore, t expression is strongly modulated by temperature in the two lines. Thus, temperature sensitivity of t expression is at least partly a consequence of bab thermal transcriptional plasticity. We therefore propose that a gene regulatory network integrating both genetic variation and temperature sensitivity modulates female abdominal pigmentation. Interestingly, both bDE and t_MSE were previously shown to have been recurrently involved in abdominal pigmentation evolution in drosophilids. We propose that the environmental sensitivity of these enhancers has turned them into evolutionary hotspots.

Complex traits such as size or disease susceptibility are typically modulated by both genetic variation and environmental conditions. Model organisms such as fruit flies (Drosophila) are particularly appropriate to analyse the interactions between genetic variation and environmental factors during the development of complex phenotypes. Natural populations carry high genetic variation and can be grown in controlled conditions in the laboratory. Here, we use Drosophila melanogaster female abdominal pigmentation, which is both genetically variable and modulated by the environment (temperature) to dissect this kind of interaction. We show that the pigmentation difference between two inbred fly lines is caused by genetic variation in an enhancer of the bab locus, which encodes two transcription factors controlling abdominal pigmentation. Indeed, this enhancer drives differential expression between the two lines. Interestingly, this enhancer is sensitive to temperature in both lines. We show that the effect of bab on pigmentation is mediated by the pigmentation gene tan (t) that is repressed by bab. Thus, the previously reported temperature-sensitive expression of t is a direct consequence of bab transcriptional plasticity.

Gene(s) and individual feeding behavior: Exploring eco-evolutionary dynamics underlying left-right asymmetry in the scale-eating cichlid fish Perissodus microlepis

The scale‐eating cichlid fish Perissodus microlepis is a textbook example of bilateral asymmetry due to its left or right‐bending heads and of negative frequency‐dependent selection, which is proposed to maintain this stable polymorphism. The mechanisms that underlie this asymmetry remain elusive. Several studies had initially postulated a simple genetic basis for this trait, but this explanation has been questioned, particularly by reports observing a unimodal distribution of mouth shapes. We hypothesize that this unimodal distribution might be due to a combination of genetic and phenotypically plastic components. Here, we expanded on previous work by investigating a formerly identified candidate SNP associated to mouth laterality, documenting inter‐individual variation in feeding preference using stable isotope analyses, and testing their association with mouth asymmetry. Our results suggest that this polymorphism is influenced by both a polygenic basis and inter‐individual non‐genetic variation, possibly due to feeding experience, individual specialization, and intraspecific competition. We introduce a hypothesis potentially explaining the simultaneous maintenance of left, right, asymmetric and symmetric mouth phenotypes due to the interaction between diverse eco‐evolutionary dynamics including niche construction and balancing selection. Future studies will have to further tease apart the relative contribution of genetic and environmental factors and their interactions in an integrated fashion.

Towards a gene regulatory network shaping the fins of the Princess cichlid

Variation in fin shape and size contributes to the outstanding morphological diversity of teleost fishes, but the regulation of fin growth has not yet been studied extensively outside the zebrafish model. A previous gene expression study addressing the ornamental elongations of unpaired fins in the African cichlid fish Neolamprologus brichardi identified three genes (cx43, mmp9 and sema3d) with strong and consistent expression differences between short and elongated fin regions. Remarkably, the expression patterns of these genes were not consistent with inferences on their regulatory interactions in zebrafish. Here, we identify a gene expression network (GRN) comprising cx43, mmp9, and possibly also sema3d by a stepwise approach of identifying co-expression modules and predicting their upstream regulators. Among the transcription factors (TFs) predicted as potential upstream regulators of 11 co-expressed genes, six TFs (foxc1, foxp1, foxd3, myc, egr2, irf8) showed expression patterns consistent with their cooperative transcriptional regulation of the gene network. Some of these TFs have already been implicated in teleost fish fin regeneration and formation. We particularly discuss the potential function of foxd3 as driver of the network and its role in the unexpected gene expression correlations observed in N. brichardi.

Why does the magnitude of genotype-by-environment interaction vary?

Genotype-by-environment interaction (G × E), that is, genetic variation in phenotypic plasticity, is a central concept in ecology and evolutionary biology. G×E has wide-ranging implications for trait development and for understanding how organisms will respond to environmental change. Although G × E has been extensively documented, its presence and magnitude vary dramatically across populations and traits. Despite this, we still know little about why G × E is so evident in some traits and populations, but minimal or absent in others. To encourage synthetic research in this area, we review diverse hypotheses for the underlying biological causes of variation in G × E. We extract common themes from these hypotheses to develop a more synthetic understanding of variation in G × E and suggest some important next steps.

Similarities in temperature-dependent gene expression plasticity across timescales in threespine stickleback (Gasterosteus aculeatus)

Phenotypic plasticity occurs at a variety of timescales, but little is known about the degree to which plastic responses at different timescales are associated with similar underlying molecular processes, which is critical for assessing the effects of plasticity on evolutionary trajectories. To address this issue, we identified differential gene expression in response to developmental temperature in the muscle transcriptome of adult threespine stickleback (Gasterosteus aculeatus) exposed to 12, 18 and 24°C until hatch and then held at 18°C for 9 months and compared these results to differential gene expression in response to adult thermal acclimation in stickleback developed at 18°C and then acclimated to 5 and 25°C as adults. Adult thermal acclimation affected the expression of 7,940 and 7,015 genes in response to cold and warm acclimation, respectively, and 4,851 of these genes responded in both treatments. In contrast, the expression of only 33 and 29 genes was affected by cold and warm development, respectively. The majority of the genes affected by developmental temperature were also affected by adult acclimation temperature. Many genes that were differentially expressed as a result of adult acclimation were associated with previously identified temperature-dependent effects on DNA methylation patterns, suggesting a role of epigenetic mechanisms in regulating gene expression plasticity during acclimation. Taken together, these results demonstrate similarities between the persistent effects of developmental plasticity on gene expression and the effects of adult thermal acclimation, emphasizing the potential for mechanistic links between plasticity acting at these different life stages.

Topological and statistical analyses of gene regulatory networks reveal unifying yet quantitatively different emergent properties

Understanding complexity in physical, biological, social and information systems is predicated on describing interactions amongst different components. Advances in genomics are facilitating the high-throughput identification of molecular interactions, and graphs are emerging as indispensable tools in explaining how the connections in the network drive organismal phenotypic plasticity. Here, we describe the architectural organization and associated emergent topological properties of gene regulatory networks (GRNs) that describe protein-DNA interactions (PDIs) in several model eukaryotes. By analyzing GRN connectivity, our results show that the anticipated scale-free network architectures are characterized by organism-specific power law scaling exponents. These exponents are independent of the fraction of the GRN experimentally sampled, enabling prediction of properties of the complete GRN for an organism. We further demonstrate that the exponents describe inequalities in transcription factor (TF)-target gene recognition across GRNs. These observations have the important biological implication that they predict the existence of an intrinsic organism-specific trans and/or cis regulatory landscape that constrains GRN topologies. Consequently, architectural GRN organization drives not only phenotypic plasticity within a species, but is also likely implicated in species-specific phenotype.

Advancing behavioural genomics by considering timescale

Animal behavioural traits often covary with gene expression, pointing towards a genomic constraint on organismal responses to environmental cues. This pattern highlights a gap in our understanding of the time course of environmentally responsive gene expression, and moreover, how these dynamics are regulated. Advances in behavioural genomics explore how gene expression dynamics are correlated with behavioural traits that range from stable to highly labile. We consider the idea that certain genomic regulatory mechanisms may predict the timescale of an environmental effect on behaviour. This temporally minded approach could inform both organismal and evolutionary questions ranging from the remediation of early life social trauma to understanding the evolution of trait plasticity.

A gene expression study of ornamental fin shape in Neolamprologus brichardi, an African cichlid species

The diversity of fin morphology within and across fish taxa offers great, but still largely unexplored, opportunities to investigate the proximate mechanisms underlying fin shape variation. Relying on available genetic knowledge brought forth mainly by the comprehensive study of the zebrafish caudal fin, we explored candidate molecular mechanisms for the maintenance and formation of the conspicuously elongated filaments adorning the unpaired fins of the East African "princess cichlid" Neolamprologus brichardi. Via qPCR assays, we detected expression differences of candidate genes between elongated and short regions of intact and regenerating fins. The identified genes include skeletogenic and growth factors (igf2b, fgf3, bmp2 and bmp4), components of the WNT pathway (lef1, wnt5b and wnt10) and a regulatory network determining fin ray segment size and junction (cx43, esco2 and sema3d), as well as other genes with different roles (mmp9, msxb and pea3). Interestingly, some of these genes showed fin specific expression differences which are often neglected in studies of model fish that focus on the caudal fin. Moreover, while the observed expression patterns were generally consistent with zebrafish results, we also detected deviating expression correlations and gene functions.

Molecular investigation of genetic assimilation during the rapid adaptive radiations of East African cichlid fishes

Adaptive radiations are characterized by adaptive diversification intertwined with rapid speciation within a lineage resulting in many ecologically specialized, phenotypically diverse species. It has been proposed that adaptive radiations can originate from ancestral lineages with pronounced phenotypic plasticity in adaptive traits, facilitating ecologically driven phenotypic diversification that is ultimately fixed through genetic assimilation of gene regulatory regions. This study aimed to investigate how phenotypic plasticity is reflected in gene expression patterns in the trophic apparatus of several lineages of East African cichlid fishes, and whether the observed patterns support genetic assimilation. This investigation used a split brood experimental design to compare adaptive plasticity in species from within and outside of adaptive radiations. The plastic response was induced in the crushing pharyngeal jaws through feeding individuals either a hard or soft diet. We find that nonradiating, basal lineages show higher levels of adaptive morphological plasticity than the derived, radiated lineages, suggesting that these differences have become partially genetically fixed during the formation of the adaptive radiations. Two candidate genes that may have undergone genetic assimilation, gif and alas1, were identified, in addition to alterations in the wiring of LPJ patterning networks. Taken together, our results suggest that genetic assimilation may have dampened the inducibility of plasticity related genes during the adaptive radiations of East African cichlids, flattening the reaction norms and canalizing their feeding phenotypes, driving adaptation to progressively more narrow ecological niches.

The Role of Alternative Splicing and Differential Gene Expression in Cichlid Adaptive Radiation

Species diverge eco-morphologically through the continuous action of natural selection on functionally important structures, producing alternative adaptive morphologies. In cichlid fishes, the oral and pharyngeal jaws are such key structures. Adaptive variation in jaw morphology contributes to trophic specialization, which is hypothesized to fuel their rapid speciation in the East African Great Lakes. Much is known about the genes involved in cichlid jaw and craniofacial development. However, it is still unclear what salient sources of variation gave rise to trophic-niche specialization, facilitating adaptive radiation. Here, we explore two sources of transcriptional variation that may underlie species-specific disparities in jaw morphology. Using whole transcriptome RNA-sequencing, we analyze differences in gene expression and alternative splicing, at the end of postlarval development, in fully functional jaws of six species of cichlids from the Lake Tanganyika tribe Tropheini. Our data reveal a surprisingly high degree of alternative splicing events compared with gene expression differences among species and trophic types. This suggests that differential trophic adaptation of the jaw apparatus may have been shaped by transcriptional rewiring of splicing as well as gene expression variation during the rapid radiation of the Tropheini. Specifically, genes undergoing splicing across most species were found to be enriched for pharyngeal jaw gene ontology terms. Overall, jaw transcriptional patterns at postlarval developmental stage were highly dynamic and species-specific. In conclusion, this work indicates that shifts in alternative splicing could have played a more important role in cichlid adaptive radiation, and possibly adaptive radiation in general, than currently recognized.

The potential role of genetic assimilation during maize domestication

Domestication research has largely focused on identification of morphological and genetic differences between extant populations of crops and their wild relatives. Little attention has been paid to the potential effects of environment despite substantial known changes in climate from the time of domestication to modern day. In recent research, the exposure of teosinte (i.e., wild maize) to environments similar to the time of domestication, resulted in a plastic induction of domesticated phenotypes in teosinte. These results suggest that early agriculturalists may have selected for genetic mechanisms that cemented domestication phenotypes initially induced by a plastic response of teosinte to environment, a process known as genetic assimilation. To better understand this phenomenon and the potential role of environment in maize domestication, we examined differential gene expression in maize (Zea mays ssp. mays) and teosinte (Zea mays ssp. parviglumis) between past and present conditions. We identified a gene set of over 2000 loci showing a change in expression across environmental conditions in teosinte and invariance in maize. In fact, overall we observed both greater plasticity in gene expression and more substantial changes in co-expressionnal networks in teosinte across environments when compared to maize. While these results suggest genetic assimilation played at least some role in domestication, genes showing expression patterns consistent with assimilation are not significantly enriched for previously identified domestication candidates, indicating assimilation did not have a genome-wide effect.

The combined effects of temporal autocorrelation and the costs of plasticity on the evolution of plasticity

Adaptive phenotypic plasticity is an important source of intraspecific variation, and for many plastic traits, the costs or factors limiting plasticity seem cryptic. However, there are several different factors that may constrain the evolution of plasticity, but few models have considered costs and limiting factors simultaneously. Here we use a simulation model to investigate how the optimal level of plasticity in a population depends on a fixed maintenance fitness cost for plasticity or an incremental fitness cost for producing a plastic response in combination with environmental unpredictability (environmental fluctuation speed) limiting plasticity. Our model identifies two mechanisms that act, almost separately, to constrain the evolution of plasticity: (i) the fitness cost of plasticity scaled by the nonplastic environmental tolerance, and (ii) the environmental fluctuation speed scaled by the rate of phenotypic change. That is, the evolution of plasticity is constrained by the high cost of plasticity in combination with high tolerance for environmental variation, or fast environmental changes in combination with slow plastic response. Qualitatively similar results are found when maintenance and incremental fitness costs of plasticity are incorporated, although a larger degree of plasticity is selected for with an incremental cost. Our model highlights that it is important to consider direct fitness costs and phenotypic limitations in relation to nonplastic environmental tolerance and environmental fluctuations, respectively, to understand what constrains the evolution of phenotypic plasticity.