Gene regulatory networks reused to build novel traits: co-option of an eye-related gene regulatory network in eye-like organs and red wing patches on insect wings is suggested by optix expression

Exploring gene regulation and biological processes in insects: Insights from omics data using gene regulatory network models

Gene regulatory network (GRN) comprises complicated yet intertwined gene-regulator relationships. Understanding the GRN dynamics will unravel the complexity behind the observed gene expressions. Insect gene regulation is often complicated due to their complex life cycles and diverse ecological adaptations. The main interest of this review is to have an update on the current mathematical modelling methods of GRNs to explain insect science. Several popular GRN architecture models are discussed, together with examples of applications in insect science. In the last part of this review, each model is compared from different aspects, including network scalability, computation complexity, robustness to noise and biological relevancy.

Gene regulatory network co-option is sufficient to induce a morphological novelty in Drosophila

Identifying the molecular origins by which new morphological structures evolve is one of the long standing problems in evolutionary biology. To date, vanishingly few examples provide a compelling account of how new morphologies were initially formed, thereby limiting our understanding of how diverse forms of life derived their complex features. Here, we provide evidence that the large projections on the Drosophila eugracilis phallus that are implicated in sexual conflict have evolved through co-option of the trichome genetic network. These unicellular apical projections on the phallus postgonal sheath are reminiscent of trichomes that cover the Drosophila body but are up to 20-fold larger in size. During their development, they express the transcription factor Shavenbaby, the master regulator of the trichome network. Consistent with the co-option of the Shavenbaby network during the evolution of the D. eugracilis projections, somatic mosaic CRISPR/Cas9 mutagenesis shows that shavenbaby is necessary for their proper length. Moreover, mis-expression of Shavenbaby in the sheath of D. melanogaster , a naïve species that lacks these extensions, is sufficient to induce small trichomes. These induced extensions rely on a genetic network that is shared to a large extent with the D. eugracilis projections, indicating its co-option but also some genetic rewiring. Thus, by leveraging a genetically tractable evolutionarily novelty, our work shows that the trichome-forming network is flexible enough that it can be co-opted in a new context, and subsequently refined to produce unique apical projections that are barely recognizable compared to their simpler ancestral beginnings.

Organogenetic transcriptomes of the Drosophila embryo at single cell resolution

To gain insight into the transcription programs activated during the formation of Drosophila larval structures, we carried out single cell RNA sequencing during two periods of Drosophila embryogenesis: stages 10-12, when most organs are first specified and initiate morphological and physiological specialization; and stages 13-16, when organs achieve their final mature architectures and begin to function. Our data confirm previous findings with regards to functional specialization of some organs - the salivary gland and trachea - and clarify the embryonic functions of another - the plasmatocytes. We also identify two early developmental trajectories in germ cells and uncover a potential role for proteolysis during germline stem cell specialization. We identify the likely cell type of origin for key components of the Drosophila matrisome and several commonly used Drosophila embryonic cell culture lines. Finally, we compare our findings with other recent related studies and with other modalities for identifying tissue-specific gene expression patterns. These data provide a useful community resource for identifying many new players in tissue-specific morphogenesis and functional specialization of developing organs.

Lateral thinking in syndromic congenital cardiovascular disease

Syndromic birth defects are rare diseases that can present with seemingly pleiotropic comorbidities. Prime examples are rare congenital heart and cardiovascular anomalies that can be accompanied by forelimb defects, kidney disorders and more. Whether such multi-organ defects share a developmental link remains a key question with relevance to the diagnosis, therapeutic intervention and long-term care of affected patients. The heart, endothelial and blood lineages develop together from the lateral plate mesoderm (LPM), which also harbors the progenitor cells for limb connective tissue, kidneys, mesothelia and smooth muscle. This developmental plasticity of the LPM, which founds on multi-lineage progenitor cells and shared transcription factor expression across different descendant lineages, has the potential to explain the seemingly disparate syndromic defects in rare congenital diseases. Combining patient genome-sequencing data with model organism studies has already provided a wealth of insights into complex LPM-associated birth defects, such as heart-hand syndromes. Here, we summarize developmental and known disease-causing mechanisms in early LPM patterning, address how defects in these processes drive multi-organ comorbidities, and outline how several cardiovascular and hematopoietic birth defects with complex comorbidities may be LPM-associated diseases. We also discuss strategies to integrate patient sequencing, data-aggregating resources and model organism studies to mechanistically decode congenital defects, including potentially LPM-associated orphan diseases. Eventually, linking complex congenital phenotypes to a common LPM origin provides a framework to discover developmental mechanisms and to anticipate comorbidities in congenital diseases affecting the cardiovascular system and beyond.

The genetic basis of wing spots in Pieris canidia butterflies

Spots in pierid butterflies and eyespots in nymphalid butterflies are likely non-homologous wing colour pattern elements, yet they share a few features in common. Both develop black scales that depend on the function of the gene spalt, and both might have central signalling cells. This suggests that both pattern elements may be sharing common genetic circuitry. Hundreds of genes have already been associated with the development of nymphalid butterfly eyespot patterns, but the genetic basis of the simpler spot patterns on the wings of pierid butterflies has not been investigated. To facilitate studies of pierid wing patterns, we report a high-quality draft genome assembly for Pieris canidia, the Indian cabbage white. We then conducted transcriptomic analyses of pupal wing tissues sampled from the spot and non-spot regions of P. canidia at 3-6 h post-pupation. A total of 1352 genes were differentially regulated between wing tissues with and without the black spot, including spalt, Krüppel-like factor 10, genes from the Toll, Notch, TGF-β, and FGFR signalling pathways, and several genes involved in the melanin biosynthetic pathway. We identified 14 genes that are up-regulated in both pierid spots and nymphalid eyespots and propose that spots and eyespots share regulatory modules despite their likely independent origins.

Evolution of gene regulatory network of C4 photosynthesis in the genus Flaveria reveals the evolutionary status of C3-C4 intermediate species

C₄ photosynthesis evolved from ancestral C₃ photosynthesis by recruiting pre-existing genes to fulfill new functions. The enzymes and transporters required for the C₄ metabolic pathway have been intensively studied and well documented; however, the transcription factors (TFs) that regulate these C₄ metabolic genes are not yet well understood. In particular, how the TF regulatory network of C₄ metabolic genes was rewired during the evolutionary process is unclear. Here, we constructed gene regulatory networks (GRNs) for four closely evolutionarily related species from the genus Flaveria, which represent four different evolutionary stages of C₄ photosynthesis: C₃ (F. robusta), type I C₃-C₄ (F. sonorensis), type II C₃-C₄ (F. ramosissima), and C₄ (F. trinervia). Our results show that more than half of the co-regulatory relationships between TFs and core C₄ metabolic genes are species specific. The counterparts of the C₄ genes in C₃ species were already co-regulated with photosynthesis-related genes, whereas the required TFs for C₄ photosynthesis were recruited later. The TFs involved in C₄ photosynthesis were widely recruited in the type I C₃-C₄ species; nevertheless, type II C₃-C₄ species showed a divergent GRN from C₄ species. In line with these findings, a ¹³CO₂ pulse-labeling experiment showed that the CO₂ initially fixed into C₄ acid was not directly released to the Calvin-Benson-Bassham cycle in the type II C₃-C₄ species. Therefore, our study uncovered dynamic changes in C₄ genes and TF co-regulation during the evolutionary process; furthermore, we showed that the metabolic pathway of the type II C₃-C₄ species F. ramosissima represents an alternative evolutionary solution to the ammonia imbalance in C₃-C₄ intermediate species.

Co-option of a carotenoid cleavage dioxygenase gene (CCD4a) into the floral symmetry gene regulatory network contributes to the polymorphic floral shape-color combinations in Chrysanthemum sensu lato

Morphological novelties, including formation of trait combinations, may result from de novo gene origination and/or co-option of existing genes into other developmental contexts. A variety of shape-color combinations of capitular florets occur in Chrysanthemum and its allies. We hypothesized that co-option of a carotenoid cleavage dioxygenase gene into the floral symmetry gene network would generate a white zygomorphic ray floret. We tested this hypothesis in an evolutionary context using species in Chrysanthemum sensu lato, a monophyletic group with diverse floral shape-color combinations, based on morphological investigation, interspecific crossing, molecular interaction and transgenic experiments. Our results showed that white color was significantly associated with floret zygomorphy. Specific expression of the carotenoid cleavage dioxygenase gene CCD4a in marginal florets resulted in white color. Crossing experiments between Chrysanthemum lavandulifolium and Ajania pacifica indicated that expression of CCD4a is trans-regulated. The floral symmetry regulator CYC2g can activate expression of CCD4a with a dependence on TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING (TCP) binding element 8 on the CCD4a promoter. Based on all experimental findings, we propose that gene co-option of carotenoid degradation into floral symmetry regulation, and the subsequent dysfunction or loss of either CYC2g or CCD4a, may have led to evolution of capitular shape-color patterning in Chrysanthemum sensu lato.

The Daphnia carapace and other novel structures evolved via the cryptic persistence of serial homologs

Understanding how novel structures arise is a central question in evolution. Novel structures are often defined as structures that are not derived from (homologous to) any structure in the ancestor.¹ The carapace of the crustacean Daphnia magna is a bivalved "cape" of exoskeleton. Shiga et al.² proposed that the carapace of crustaceans like Daphnia and many other plate-like outgrowths in arthropods are novel structures that arose through the repeated co-option of genes like vestigial that also pattern insect wings.^2-4 To determine whether the Daphnia carapace is a novel structure, we compare previous functional work² with the expression of genes known to pattern the proximal leg region (pannier, araucan, and vestigial)^5,6 between Daphnia, Parhyale, and Tribolium. Our results suggest that the Daphnia carapace did not arise by co-option but instead derived from an exite (lateral leg lobe) that emerges from an ancestral proximal leg segment that was incorporated into the Daphnia body wall. The Daphnia carapace, therefore, appears to be homologous to the Parhyale tergal plate and the insect wing.⁵ Remarkably, the vestigial-positive tissue that gives rise to the Daphnia carapace appears to be present in Parhyale⁷ and Tribolium as a small, inconspicuous protrusion. Thus, rather than a novel structure resulting from gene co-option, the Daphnia carapace appears to have arisen from a shared, ancestral tissue (morphogenetic field) that persists in a cryptic state in other arthropod lineages. Cryptic persistence of unrecognized serial homologs may thus be a general solution for the origin of novel structures.

Out from under the wing: reconceptualizing the insect wing gene regulatory network as a versatile, general module for body-wall lobes in arthropods

Body plan evolution often occurs through the differentiation of serially homologous body parts, particularly in the evolution of arthropod body plans. Recently, homeotic transformations resulting from experimental manipulation of gene expression, along with comparative data on the expression and function of genes in the wing regulatory network, have provided a new perspective on an old question in insect evolution: how did the insect wing evolve? We investigated the metamorphic roles of a suite of 10 wing- and body-wall-related genes in a hemimetabolous insect, Oncopeltus fasciatus. Our results indicate that genes involved in wing development in O. fasciatus play similar roles in the development of adult body-wall flattened cuticular evaginations. We found extensive functional similarity between the development of wings and other bilayered evaginations of the body wall. Overall, our results support the existence of a versatile development module for building bilayered cuticular epithelial structures that pre-dates the evolutionary origin of wings. We explore the consequences of reconceptualizing the canonical wing-patterning network as a bilayered body-wall patterning network, including consequences for long-standing debates about wing homology, the origin of wings and the origin of novel bilayered body-wall structures. We conclude by presenting three testable predictions that result from this reconceptualization.

Single-cell RNA sequencing of the Strongylocentrotus purpuratus larva reveals the blueprint of major cell types and nervous system of a non-chordate deuterostome

Identifying the molecular fingerprint of organismal cell types is key for understanding their function and evolution. Here, we use single-cell RNA sequencing (scRNA-seq) to survey the cell types of the sea urchin early pluteus larva, representing an important developmental transition from non-feeding to feeding larva. We identify 21 distinct cell clusters, representing cells of the digestive, skeletal, immune, and nervous systems. Further subclustering of these reveal a highly detailed portrait of cell diversity across the larva, including the identification of neuronal cell types. We then validate important gene regulatory networks driving sea urchin development and reveal new domains of activity within the larval body. Focusing on neurons that co-express Pdx-1 and Brn1/2/4, we identify an unprecedented number of genes shared by this population of neurons in sea urchin and vertebrate endocrine pancreatic cells. Using differential expression results from Pdx-1 knockdown experiments, we show that Pdx1 is necessary for the acquisition of the neuronal identity of these cells. We hypothesize that a network similar to the one orchestrated by Pdx1 in the sea urchin neurons was active in an ancestral cell type and then inherited by neuronal and pancreatic developmental lineages in sea urchins and vertebrates.

Gene regulation of adult skeletogenesis in starfish and modifications during gene network co-option

The larval skeleton of the echinoderm is believed to have been acquired through co-option of a pre-existing gene regulatory network (GRN); that is, the mechanism for adult skeleton formation in the echinoderm was deployed in early embryogenesis during echinoderm diversification. To explore the evolutionary changes that occurred during co-option, we examined the mechanism for adult skeletogenesis using the starfish Patiria pectinifera. Expression patterns of skeletogenesis-related genes (vegf, vegfr, ets1/2, erg, alx1, ca1, and clect) suggest that adult skeletogenic cells develop from the posterior coelom after the start of feeding. Treatment with inhibitors and gene knockout using transcription activator-like effector nucleases (TALENs) suggest that the feeding-nutrient sensing pathway activates Vegf signaling via target of rapamycin (TOR) activity, leading to the activation of skeletogenic regulatory genes in starfish. In the larval skeletogenesis of sea urchins, the homeobox gene pmar1 activates skeletogenic regulatory genes, but in starfish, localized expression of the pmar1-related genes phbA and phbB was not detected during the adult skeleton formation stage. Based on these data, we provide a model for the adult skeletogenic GRN in the echinoderm and propose that the upstream regulatory system changed from the feeding-TOR-Vegf pathway to a homeobox gene-system during co-option of the skeletogenic GRN.

Homology of process: developmental dynamics in comparative biology

Comparative biology builds up systematic knowledge of the diversity of life, across evolutionary lineages and levels of organization, starting with evidence from a sparse sample of model organisms. In developmental biology, a key obstacle to the growth of comparative approaches is that the concept of homology is not very well defined for levels of organization that are intermediate between individual genes and morphological characters. In this paper, we investigate what it means for ontogenetic processes to be homologous, focusing specifically on the examples of insect segmentation and vertebrate somitogenesis. These processes can be homologous without homology of the underlying genes or gene networks, since the latter can diverge over evolutionary time, while the dynamics of the process remain the same. Ontogenetic processes like these therefore constitute a dissociable level and distinctive unit of comparison requiring their own specific criteria of homology. In addition, such processes are typically complex and nonlinear, such that their rigorous description and comparison requires not only observation and experimentation, but also dynamical modelling. We propose six criteria of process homology, combining recognized indicators (sameness of parts, morphological outcome and topological position) with novel ones derived from dynamical systems modelling (sameness of dynamical properties, dynamical complexity and evidence for transitional forms). We show how these criteria apply to animal segmentation and other ontogenetic processes. We conclude by situating our proposed dynamical framework for homology of process in relation to similar research programmes, such as process structuralism and developmental approaches to morphological homology.

Visual mate preference evolution during butterfly speciation is linked to neural processing genes

Many animal species remain separate not because their individuals fail to produce viable hybrids but because they "choose" not to mate. However, we still know very little of the genetic mechanisms underlying changes in these mate preference behaviours. Heliconius butterflies display bright warning patterns, which they also use to recognize conspecifics. Here, we couple QTL for divergence in visual preference behaviours with population genomic and gene expression analyses of neural tissue (central brain, optic lobes and ommatidia) across development in two sympatric Heliconius species. Within a region containing 200 genes, we identify five genes that are strongly associated with divergent visual preferences. Three of these have previously been implicated in key components of neural signalling (specifically an ionotropic glutamate receptor and two regucalcins), and overall our candidates suggest shifts in behaviour involve changes in visual integration or processing. This would allow preference evolution without altering perception of the wider environment.

Tissue-specific developmental regulation and isoform usage underlie the role of doublesex in sex differentiation and mimicry in Papilio swallowtails

Adaptive phenotypes often arise by rewiring existing developmental networks. Co-option of transcription factors in novel contexts has facilitated the evolution of ecologically important adaptations. doublesex (dsx) governs fundamental sex differentiation during embryonic stages and has been co-opted to regulate diverse secondary sexual dimorphisms during pupal development of holometabolous insects. In Papilio polytes, dsx regulates female-limited mimetic polymorphism, resulting in mimetic and non-mimetic forms. To understand how a critical gene such as dsx regulates novel wing patterns while maintaining its basic function in sex differentiation, we traced its expression through metamorphosis in P. polytes using developmental transcriptome data. We found three key dsx expression peaks: (i) eggs in pre- and post-ovisposition stages; (ii) developing wing discs and body in final larval instar; and (iii) 3-day pupae. We identified potential dsx targets using co-expression and differential expression analysis, and found distinct, non-overlapping sets of genes-containing putative dsx-binding sites-in developing wings versus abdominal tissue and in mimetic versus non-mimetic individuals. This suggests that dsx regulates distinct downstream targets in different tissues and wing colour morphs and has perhaps acquired new, previously unknown targets, for regulating mimetic polymorphism. Additionally, we observed that the three female isoforms of dsx were differentially expressed across stages (from eggs to adults) and tissues and differed in their protein structure. This may promote differential protein-protein interactions for each isoform and facilitate sub-functionalization of dsx activity across its isoforms. Our findings suggest that dsx employs tissue-specific downstream effectors and partitions its functions across multiple isoforms to regulate primary and secondary sexual dimorphism through insect development.

Transcriptomic similarities and differences between the limb bud AER and unique carapacial ridge of turtle embryos

Evolutionary innovation may arise via major departures from an ancestral condition. Turtle shell morphogenesis depends on a unique structure known as the carapacial ridge (CR). This lateral tissue protrusion in turtle embryos exhibits similar properties as the apical ectodermal ridge (AER)-a well-known molecular signaling center involved in limb development. Still, how the CR influences shell morphogenesis is not entirely clear. The present study aimed to describe the CR transcriptome shortly before ribs were halted within its mesenchyme, as required for shell development. Analyses exposed that the mesenchymal marker VIM was one of the most highly co-expressed genes and numerous appendage formation genes were situated within the core of CR and AER co-expression networks. However, there were tissue-specific differences in the activity of these genes. For instance, WNT5A was most frequently assigned to appendage-related annotations of the CR network core, but not in the AER. Several homeobox transcription factors known to regulate limb bud patterning exhibited their highest expression levels in the AER, but were underexpressed in the CR. The results of this study corroborate that novel body plans often originate via alterations of pre-existing genetic networks. Altogether, this exploratory study enhances the groundwork for future experiments on the molecular underpinnings of turtle shell development and evolution.

Temporal flexibility of gene regulatory network underlies a novel wing pattern in flies

Developmental genes can be coopted to generate evolutionary novelties by changing their spatial regulation. However, developmental genes seldom act independently, but rather work in a gene regulatory network (GRN). How is it possible to recruit a single gene from a whole GRN? What are the properties that allow parallel cooptions of the same genes during evolution? Here, we show that a novel engrailed gene expression underlies a novel wing color pattern in flies. We show that cooption is facilitated 1) because of GRN flexibility over development and 2) because every single gene of the GRN has its own functional time window. We suggest these two temporal properties could explain why the same gene can be independently recruited several times during evolution.

Organisms have evolved endless morphological, physiological, and behavioral novel traits during the course of evolution. Novel traits were proposed to evolve mainly by orchestration of preexisting genes. Over the past two decades, biologists have shown that cooption of gene regulatory networks (GRNs) indeed underlies numerous evolutionary novelties. However, very little is known about the actual GRN properties that allow such redeployment. Here we have investigated the generation and evolution of the complex wing pattern of the fly Samoaia leonensis. We show that the transcription factor Engrailed is recruited independently from the other players of the anterior–posterior specification network to generate a new wing pattern. We argue that partial cooption is made possible because 1) the anterior–posterior specification GRN is flexible over time in the developing wing and 2) this flexibility results from the fact that every single gene of the GRN possesses its own functional time window. We propose that the temporal flexibility of a GRN is a general prerequisite for its possible cooption during the course of evolution.

On the specificity of gene regulatory networks: How does network co-option affect subsequent evolution?

The process of multicellular organismal development hinges upon the specificity of developmental programs: for different parts of the organism to form unique features, processes must exist to specify each part. This specificity is thought to be hardwired into gene regulatory networks, which activate cohorts of genes in particular tissues at particular times during development. However, the evolution of gene regulatory networks sometimes occurs by mechanisms that sacrifice specificity. One such mechanism is network co-option, in which existing gene networks are redeployed in new developmental contexts. While network co-option may offer an efficient mechanism for generating novel phenotypes, losses of tissue specificity at redeployed network genes could restrict the ability of the affected traits to evolve independently. At present, there has not been a detailed discussion regarding how tissue specificity of network genes might be altered due to gene network co-option at its initiation, as well as how trait independence can be retained or restored after network co-option. A lack of clarity about network co-option makes it more difficult to speculate on the long-term evolutionary implications of this mechanism. In this review, we will discuss the possible initial outcomes of network co-option, outline the mechanisms by which networks may retain or subsequently regain specificity after network co-option, and comment on some of the possible evolutionary consequences of network co-option. We place special emphasis on the need to consider selectively-neutral outcomes of network co-option to improve our understanding of the role of this mechanism in trait evolution.

The Invariant Nature of a Morphological Character and Character State: Insights from Gene Regulatory Networks

What constitutes a discrete morphological character versus character state has been long discussed in the systematics literature but the consensus on this issue is still missing. Different methods of classifying organismal features into characters and character states (CCSs) can dramatically affect the results of phylogenetic analyses. Here, I show that, in the framework of Markov models, the modular structure of the gene regulatory network (GRN) underlying trait development, and the hierarchical nature of GRN evolution, essentially remove the distinction between morphological CCS, thus endowing the CCS with an invariant property with respect to each other. This property allows the states of one character to be represented as several individual characters and vice versa. In practice, this means that a phenotype can be encoded using a set of characters or just one complex character with numerous states. The representation of a phenotype using one complex character can be implemented in Markov models of trait evolution by properly structuring transition rate matrix.

A transposable element insertion is associated with an alternative life history strategy

Tradeoffs affect resource allocation during development and result in fitness consequences that drive the evolution of life history strategies. Yet despite their importance, we know little about the mechanisms underlying life history tradeoffs. Many species of Colias butterflies exhibit an alternative life history strategy (ALHS) where females divert resources from wing pigment synthesis to reproductive and somatic development. Due to this reallocation, a wing color polymorphism is associated with the ALHS: either yellow/orange or white. Here we map the locus associated with this ALHS in Colias crocea to a transposable element insertion located downstream of the Colias homolog of BarH-1, a homeobox transcription factor. Using CRISPR/Cas9 gene editing, antibody staining, and electron microscopy we find white-specific expression of BarH-1 suppresses the formation of pigment granules in wing scales and gives rise to white wing color. Lipid and transcriptome analyses reveal physiological differences associated with the ALHS. Together, these findings characterize a mechanism for a female-limited ALHS.

Tradeoffs are central to life history theory and evolutionary biology, yet almost nothing is known about their mechanistic basis. Here the authors characterize one such mechanism and find a transposable element insertion is associated with the switch between alternative life history strategies.

Integration of Anatomy Ontologies and Evo-Devo Using Structured Markov Models Suggests a New Framework for Modeling Discrete Phenotypic Traits

Modeling discrete phenotypic traits for either ancestral character state reconstruction or morphology-based phylogenetic inference suffers from ambiguities of character coding, homology assessment, dependencies, and selection of adequate models. These drawbacks occur because trait evolution is driven by two key processes-hierarchical and hidden-which are not accommodated simultaneously by the available phylogenetic methods. The hierarchical process refers to the dependencies between anatomical body parts, while the hidden process refers to the evolution of gene regulatory networks (GRNs) underlying trait development. Herein, I demonstrate that these processes can be efficiently modeled using structured Markov models (SMM) equipped with hidden states, which resolves the majority of the problems associated with discrete traits. Integration of SMM with anatomy ontologies can adequately incorporate the hierarchical dependencies, while the use of the hidden states accommodates hidden evolution of GRNs and substitution rate heterogeneity. I assess the new models using simulations and theoretical synthesis. The new approach solves the long-standing "tail color problem," in which the trait is scored for species with tails of different colors or no tails. It also presents a previously unknown issue called the "two-scientist paradox," in which the nature of coding the trait and the hidden processes driving the trait's evolution are confounded; failing to account for the hidden process may result in a bias, which can be avoided by using hidden state models. All this provides a clear guideline for coding traits into characters. This article gives practical examples of using the new framework for phylogenetic inference and comparative analysis.

A homeotic shift late in development drives mimetic color variation in a bumble bee

Natural phenotypic radiations, with their high diversity and convergence, are well-suited for informing how genomic changes translate to natural phenotypic variation. New genomic tools enable discovery in such traditionally nonmodel systems. Here, we characterize the genomic basis of color pattern variation in bumble bees (Hymenoptera, Apidae, Bombus), a group that has undergone extensive convergence of setal color patterns as a result of Müllerian mimicry. In western North America, multiple species converge on local mimicry patterns through parallel shifts of midabdominal segments from red to black. Using genome-wide association, we establish that a cis-regulatory locus between the abdominal fate-determining Hox genes, abd-A and Abd-B, controls the red-black color switch in a western species, Bombus melanopygus Gene expression analysis reveals distinct shifts in Abd-B aligned with the duration of setal pigmentation at the pupal-adult transition. This results in atypical anterior Abd-B expression, a late developmental homeotic shift. Changing expression of Hox genes can have widespread effects, given their important role across segmental phenotypes; however, the late timing reduces this pleiotropy, making Hox genes suitable targets. Analysis of this locus across mimics and relatives reveals that other species follow independent genetic routes to obtain the same phenotypes.

Encounters across networks: Windows into principles of genomic regulation

Gene regulatory networks account for the ability of the genome to program development in complex multi-cellular organisms. Such networks are based on principles of gene regulation by combinations of transcription factors that bind to specific cis-regulatory DNA sites to activate transcription. These cis-regulatory regions mediate logic processing at each network node, enabling progressive increases in organismal complexity with development. Gene regulatory network explanations of development have been shown to account for patterning and cell type diversification in fly and sea urchin embryonic systems, where networks are characterized by fast coupling between transcriptional inputs and changes in target gene transcription rates, and crucial cis-regulatory elements are concentrated relatively close to the protein coding sequences of the target genes, thus facilitating their identification. Stem cell-based development in post-embryonic mammalian systems also depends on gene networks, but differs from the fly and sea urchin systems. First, the number of regulatory elements per gene and the distances between regulatory elements and the genes they control are considerably larger, forcing searches via genome-wide transcription factor binding surveys rather than functional assays. Second, the intrinsic timing of network state transitions can be slowed considerably by the need to undo stem-cell chromatin configurations, which presumably add stability to stem-cell states but retard responses to transcription factor changes during differentiation. The dispersed, partially redundant cis-regulatory systems controlling gene expression and the slow state transition kinetics in these systems already reveal new insights and opportunities to extend understanding of the repertoire of gene networks and regulatory system logic.

Gene networks and the evolution of plant morphology

Elaboration of morphology depends on the precise orchestration of gene expression by key regulatory genes. The hierarchy and relationship among the participating genes is commonly known as gene regulatory network (GRN). Therefore, the evolution of morphology ultimately occurs by the rewiring of gene network structures or by the co-option of gene networks to novel domains. The availability of high-resolution expression data combined with powerful statistical tools have opened up new avenues to formulate and test hypotheses on how diverse gene networks influence trait development and diversity. Here we summarize recent studies based on both big-data and genetics approaches to understand the evolution of plant form and physiology. We also discuss recent genome-wide investigations on how studying open-chromatin regions may help study the evolution of gene expression patterns.

Teleost Fish-Specific Preferential Retention of Pigmentation Gene-Containing Families After Whole Genome Duplications in Vertebrates

Vertebrate pigmentation is a highly diverse trait mainly determined by neural crest cell derivatives. It has been suggested that two rounds (1R/2R) of whole-genome duplications (WGDs) at the basis of vertebrates allowed changes in gene regulation associated with neural crest evolution. Subsequently, the teleost fish lineage experienced other WGDs, including the teleost-specific Ts3R before teleost radiation and the more recent Ss4R at the basis of salmonids. As the teleost lineage harbors the highest number of pigment cell types and pigmentation diversity in vertebrates, WGDs might have contributed to the evolution and diversification of the pigmentation gene repertoire in teleosts. We have compared the impact of the basal vertebrate 1R/2R duplications with that of the teleost-specific Ts3R and salmonid-specific Ss4R WGDs on 181 gene families containing genes involved in pigmentation. We show that pigmentation genes (PGs) have been globally more frequently retained as duplicates than other genes after Ts3R and Ss4R but not after the early 1R/2R. This is also true for non-pigmentary paralogs of PGs, suggesting that the function in pigmentation is not the sole key driver of gene retention after WGDs. On the long-term, specific categories of PGs have been repeatedly preferentially retained after ancient 1R/2R and Ts3R WGDs, possibly linked to the molecular nature of their proteins (e.g., DNA binding transcriptional regulators) and their central position in protein-protein interaction networks. Taken together, our results support a major role of WGDs in the diversification of the pigmentation gene repertoire in the teleost lineage, with a possible link with the diversity of pigment cell lineages observed in these animals compared to other vertebrates.

Taxon-restricted genes at the origin of a novel trait allowing access to a new environment

Taxon-restricted genes make up a considerable proportion of genomes, yet their contribution to phenotypic evolution is poorly understood. We combined gene expression with functional and behavioral assays to study the origin and adaptive value of an evolutionary innovation exclusive to the water strider genus Rhagovelia: the propelling fan. We discovered that two taxon-restricted genes, which we named geisha and mother-of-geisha, specifically control fan development. geisha originated through a duplication event at the base of the Rhagovelia lineage, and both duplicates acquired a novel expression in a specific cell population prefiguring fan development. These gene duplicates played a central role in Rhagovelia's adaptation to a new physical environment, demonstrating that the evolution of taxon-restricted genes can contribute directly to evolutionary novelties that allow access to unexploited ecological niches.

Waiting in the wings: what can we learn about gene co-option from the diversification of butterfly wing patterns?

A major challenge is to understand how conserved gene regulatory networks control the wonderful diversity of form that we see among animals and plants. Butterfly wing patterns are an excellent example of this diversity. Butterfly wings form as imaginal discs in the caterpillar and are constructed by a gene regulatory network, much of which is conserved across the holometabolous insects. Recent work in Heliconius butterflies takes advantage of genomic approaches and offers insights into how the diversification of wing patterns is overlaid onto this conserved network. WntA is a patterning morphogen that alters spatial information in the wing. Optix is a transcription factor that acts later in development to paint specific wing regions red. Both of these loci fit the paradigm of conserved protein-coding loci with diverse regulatory elements and developmental roles that have taken on novel derived functions in patterning wings. These discoveries offer insights into the 'Nymphalid Ground Plan', which offers a unifying hypothesis for pattern formation across nymphalid butterflies. These loci also represent 'hotspots' for morphological change that have been targeted repeatedly during evolution. Both convergent and divergent evolution of a great diversity of patterns is controlled by complex alleles at just a few genes. We suggest that evolutionary change has become focused on one or a few genetic loci for two reasons. First, pre-existing complex cis-regulatory loci that already interact with potentially relevant transcription factors are more likely to acquire novel functions in wing patterning. Second, the shape of wing regulatory networks may constrain evolutionary change to one or a few loci. Overall, genomic approaches that have identified wing patterning loci in these butterflies offer broad insight into how gene regulatory networks evolve to produce diversity.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Single master regulatory gene coordinates the evolution and development of butterfly color and iridescence

The optix gene has been implicated in butterfly wing pattern adaptation by genetic association, mapping, and expression studies. The actual developmental function of this gene has remained unclear, however. Here we used CRISPR/Cas9 genome editing to show that optix plays a fundamental role in nymphalid butterfly wing pattern development, where it is required for determination of all chromatic coloration. optix knockouts in four species show complete replacement of color pigments with melanins, with corresponding changes in pigment-related gene expression, resulting in black and gray butterflies. We also show that optix simultaneously acts as a switch gene for blue structural iridescence in some butterflies, demonstrating simple regulatory coordination of structural and pigmentary coloration. Remarkably, these optix knockouts phenocopy the recurring "black and blue" wing pattern archetype that has arisen on many independent occasions in butterflies. Here we demonstrate a simple genetic basis for structural coloration, and show that optix plays a deeply conserved role in butterfly wing pattern development.

Control of Hoxd gene transcription in the mammary bud by hijacking a preexisting regulatory landscape

Vertebrate Hox genes encode transcription factors operating during the development of multiple organs and structures. However, the evolutionary mechanism underlying this remarkable pleiotropy remains to be fully understood. Here, we show that Hoxd8 and Hoxd9, two genes of the HoxD complex, are transcribed during mammary bud (MB) development. However, unlike in other developmental contexts, their coexpression does not rely on the same regulatory mechanism. Hoxd8 is regulated by the combined activity of closely located sequences and the most distant telomeric gene desert. On the other hand, Hoxd9 is controlled by an enhancer-rich region that is also located within the telomeric gene desert but has no impact on Hoxd8 transcription, thus constituting an exception to the global regulatory logic systematically observed at this locus. The latter DNA region is also involved in Hoxd gene regulation in other contexts and strongly interacts with Hoxd9 in all tissues analyzed thus far, indicating that its regulatory activity was already operational before the appearance of mammary glands. Within this DNA region and neighboring a strong limb enhancer, we identified a short sequence conserved in therian mammals and capable of enhancer activity in the MBs. We propose that Hoxd gene regulation in embryonic MBs evolved by hijacking a preexisting regulatory landscape that was already at work before the emergence of mammals in structures such as the limbs or the intestinal tract.

Making sense of transcription networks

When transcription regulatory networks are compared among distantly related eukaryotes, a number of striking similarities are observed: a larger-than-expected number of genes, extensive overlapping connections, and an apparently high degree of functional redundancy. It is often assumed that the complexity of these networks represents optimized solutions, precisely sculpted by natural selection; their common features are often asserted to be adaptive. Here, we discuss support for an alternative hypothesis: the common structural features of transcription networks arise from evolutionary trajectories of "least resistance"--that is, the relative ease with which certain types of network structures are formed during their evolution.

Transcriptional control of sexual development in Cryptococcus neoformans

Developmental processes are essential for the normal life cycles of many pathogenic fungi, and they can facilitate survival in challenging environments, including the human host. Sexual development of the human fungal pathogen Cryptococcus neoformans not only produces infectious particles (spores) but has also enabled the evolution of new disease-related traits such as drug resistance. Transcription factor networks are essential to the development and pathogenesis of C. neoformans, and a variety of sequence-specific DNA-binding proteins control both key developmental transitions and virulence by regulating the expression of their target genes. In this review we discuss the roles of known transcription factors that harbor important connections to both development and virulence. Recent studies of these transcription factors have identified a common theme in which metabolic, stress, and other responses that are required for sexual development appear to have been co-opted for survival in the human host, thus facilitating pathogenesis. Future work elucidating the connection between development and pathogenesis will provide vital insights into the evolution of complex traits in eukaryotes as well as mechanisms that may be used to combat fungal pathogens.

Transcriptome analysis of the painted lady butterfly, Vanessa cardui during wing color pattern development

BACKGROUND

Butterfly wing color patterns are an important model system for understanding the evolution and development of morphological diversity and animal pigmentation. Wing color patterns develop from a complex network composed of highly conserved patterning genes and pigmentation pathways. Patterning genes are involved in regulating pigment synthesis however the temporal expression dynamics of these interacting networks is poorly understood. Here, we employ next generation sequencing to examine expression patterns of the gene network underlying wing development in the nymphalid butterfly, Vanessa cardui.

RESULTS

We identified 9, 376 differentially expressed transcripts during wing color pattern development, including genes involved in patterning, pigmentation and gene regulation. Differential expression of these genes was highest at the pre-ommochrome stage compared to early pupal and late melanin stages. Overall, an increasing number of genes were down-regulated during the progression of wing development. We observed dynamic expression patterns of a large number of pigment genes from the ommochrome, melanin and also pteridine pathways, including contrasting patterns of expression for paralogs of the yellow gene family. Surprisingly, many patterning genes previously associated with butterfly pattern elements were not significantly up-regulated at any time during pupation, although many other transcription factors were differentially expressed. Several genes involved in Notch signaling were significantly up-regulated during the pre-ommochrome stage including slow border cells, bunched and pebbles; the function of these genes in the development of butterfly wings is currently unknown. Many genes involved in ecdysone signaling were also significantly up-regulated during early pupal and late melanin stages and exhibited opposing patterns of expression relative to the ecdysone receptor. Finally, a comparison across four butterfly transcriptomes revealed 28 transcripts common to all four species that have no known homologs in other metazoans.

CONCLUSIONS

This study provides a comprehensive list of differentially expressed transcripts during wing development, revealing potential candidate genes that may be involved in regulating butterfly wing patterns. Some differentially expressed genes have no known homologs possibly representing genes unique to butterflies. Results from this study also indicate that development of nymphalid wing patterns may arise not only from melanin and ommochrome pigments but also the pteridine pigment pathway.

The functional basis of wing patterning in Heliconius butterflies: the molecules behind mimicry

Wing-pattern mimicry in butterflies has provided an important example of adaptation since Charles Darwin and Alfred Russell Wallace proposed evolution by natural selection >150 years ago. The neotropical butterfly genus Heliconius played a central role in the development of mimicry theory and has since been studied extensively in the context of ecology and population biology, behavior, and mimicry genetics. Heliconius species are notable for their diverse color patterns, and previous crossing experiments revealed that much of this variation is controlled by a small number of large-effect, Mendelian switch loci. Recent comparative analyses have shown that the same switch loci control wing-pattern diversity throughout the genus, and a number of these have now been positionally cloned. Using a combination of comparative genetic mapping, association tests, and gene expression analyses, variation in red wing patterning throughout Heliconius has been traced back to the action of the transcription factor optix. Similarly, the signaling ligand WntA has been shown to control variation in melanin patterning across Heliconius and other butterflies. Our understanding of the molecular basis of Heliconius mimicry is now providing important insights into a variety of additional evolutionary phenomena, including the origin of supergenes, the interplay between constraint and evolvability, the genetic basis of convergence, the potential for introgression to facilitate adaptation, the mechanisms of hybrid speciation in animals, and the process of ecological speciation.

How Do Genomes Create Novel Phenotypes? Insights from the Loss of the Worker Caste in Ant Social Parasites

A central goal of biology is to uncover the genetic basis for the origin of new phenotypes. A particularly effective approach is to examine the genomic architecture of species that have secondarily lost a phenotype with respect to their close relatives. In the eusocial Hymenoptera, queens and workers have divergent phenotypes that may be produced via either expression of alternative sets of caste-specific genes and pathways or differences in expression patterns of a shared set of multifunctional genes. To distinguish between these two hypotheses, we investigated how secondary loss of the worker phenotype in workerless ant social parasites impacted genome evolution across two independent origins of social parasitism in the ant genera Pogonomyrmex and Vollenhovia. We sequenced the genomes of three social parasites and their most-closely related eusocial host species and compared gene losses in social parasites with gene expression differences between host queens and workers. Virtually all annotated genes were expressed to some degree in both castes of the host, with most shifting in queen-worker bias across developmental stages. As a result, despite >1 My of divergence from the last common ancestor that had workers, the social parasites showed strikingly little evidence of gene loss, damaging mutations, or shifts in selection regime resulting from loss of the worker caste. This suggests that regulatory changes within a multifunctional genome, rather than sequence differences, have played a predominant role in the evolution of social parasitism, and perhaps also in the many gains and losses of phenotypes in the social insects.

Transcriptomics of two evolutionary novelties: how to make a sperm-transfer organ out of an anal fin and a sexually selected "sword" out of a caudal fin

Swords are exaggerated male ornaments of swordtail fishes that have been of great interest to evolutionary biologists ever since Darwin described them in the Descent of Man (1871). They are a novel sexually selected trait derived from modified ventral caudal fin rays and are only found in the genus Xiphophorus. Another phylogenetically more widespread and older male trait is the gonopodium, an intromittent organ found in all poeciliid fishes, that is derived from a modified anal fin. Despite many evolutionary and behavioral studies on both traits, little is known so far about the molecular mechanisms underlying their development. By investigating transcriptomic changes (utilizing a RNA-Seq approach) in response to testosterone treatment in the swordtail fish, Xiphophorus hellerii, we aimed to better understand the architecture of the gene regulatory networks underpinning the development of these two evolutionary novelties. Large numbers of genes with tissue-specific expression patterns were identified. Among the "sword genes" those involved in embryonic organ development, sexual character development and coloration were highly expressed, while in the gonopodium rather more morphogenesis-related genes were found. Interestingly, many genes and genetic pathways are shared between both developing novel traits derived from median fins: the sword and the gonopodium. Our analyses show that a larger set of gene networks was co-opted during the development and evolution of the "older" gonopodium than in the "younger," and morphologically less complex trait, the sword. We provide a catalog of candidate genes for future efforts to dissect the development of those sexually selected exaggerated male traits in swordtails.

Origin, development, and evolution of butterfly eyespots

This article reviews the latest developments in our understanding of the origin, development, and evolution of nymphalid butterfly eyespots. Recent contributions to this field include insights into the evolutionary and developmental origin of eyespots and their ancestral deployment on the wing, the evolution of eyespot number and eyespot sexual dimorphism, and the identification of genes affecting eyespot development and black pigmentation. I also compare features of old and more recently proposed models of eyespot development and propose a schematic for the genetic regulatory architecture of eyespots. Using this schematic I propose two hypotheses for why we observe limits to morphological diversity across these serially homologous traits.

Metamorphic labral axis patterning in the beetle Tribolium castaneum requires multiple upstream, but few downstream, genes in the appendage patterning network

The arthropod labrum is an anterior appendage-like structure that forms the dorsal side of the preoral cavity. Conflicting interpretations of fossil, nervous system, and developmental data have led to a proliferation of scenarios for labral evolution. The best supported hypothesis is that the labrum is a novel structure that shares development with appendages as a result of co-option. Here, we use RNA interference in the red flour beetle Tribolium castaneum to compare metamorphic patterning of the labrum to previously published data on ventral appendage patterning. As expected under the co-option hypothesis, depletion of several genes resulted in similar defects in the labrum and ventral appendages. These include proximal deletions and proximal-to-distal transformations resulting from depletion of the leg gap genes homothorax and extradenticle, large-scale deletions resulting from depletion of the leg gap gene Distal-less, and smaller distal deletions resulting from knockdown of the EGF ligand Keren. However, depletion of dachshund and many of the genes that function downstream of the leg gap genes in the ventral appendages had either subtle or no effects on labral axis patterning. This pattern of partial similarity suggests that upstream genes act through different downstream targets in the labrum. We also discovered that many appendage axis patterning genes have roles in patterning the epipharyngeal sensillum array, suggesting that they have become integrated into a novel regulatory network. These genes include Notch, Delta, and decapentaplegic, and the transcription factors abrupt, bric à brac, homothorax, extradenticle and the paralogs apterous a and apterous b.

Evolution of sexes from an ancestral mating-type specification pathway

Male and female sexes have evolved repeatedly in eukaryotes but the origins of dimorphic sexes and their relationship to mating types in unicellular species are not understood. Volvocine algae include isogamous species such as Chlamydomonas reinhardtii, with two equal-sized mating types, and oogamous multicellular species such as Volvox carteri with sperm-producing males and egg-producing females. Theoretical work predicts genetic linkage of a gamete cell-size regulatory gene(s) to an ancestral mating-type locus as a possible step in the evolution of dimorphic gametes, but this idea has not been tested. Here we show that, contrary to predictions, a single conserved mating locus (MT) gene in volvocine algae-MID, which encodes a RWP-RK domain transcription factor-evolved from its ancestral role in C. reinhardtii as a mating-type specifier, to become a determinant of sperm and egg development in V. carteri. Transgenic female V. carteri expressing male MID produced functional sperm packets during sexual development. Transgenic male V. carteri with RNA interference (RNAi)-mediated knockdowns of VcMID produced functional eggs, or self-fertile hermaphrodites. Post-transcriptional controls were found to regulate cell-type-limited expression and nuclear localization of VcMid protein that restricted its activity to nuclei of developing male germ cells and sperm. Crosses with sex-reversed strains uncoupled sex determination from sex chromosome identity and revealed gender-specific roles for male and female mating locus genes in sexual development, gamete fitness and reproductive success. Our data show genetic continuity between the mating-type specification and sex determination pathways of volvocine algae, and reveal evidence for gender-specific adaptations in the male and female mating locus haplotypes of Volvox. These findings will enable a deeper understanding of how a master regulator of mating-type determination in an ancestral unicellular species was reprogrammed to control sexually dimorphic gamete development in a multicellular descendant.

Multiple recent co-options of Optix associated with novel traits in adaptive butterfly wing radiations

Background

While the ecological factors that drive phenotypic radiations are often well understood, less is known about the generative mechanisms that cause the emergence and subsequent diversification of novel features. Heliconius butterflies display an extraordinary diversity of wing patterns due in part to mimicry and sexual selection. Identifying the genetic drivers of this crucible of evolution is now within reach, as it was recently shown that cis-regulatory variation of the optix transcription factor explains red pattern differences in the adaptive radiations of the Heliconius melpomene and Heliconius erato species groups.

Results

Here, we compare the developmental expression of the Optix protein across a large phylogenetic sample of butterflies and infer that its color patterning role originated at the base of the neotropical passion-vine butterfly clade (Lepidoptera, Nymphalidae, Tribe: Heliconiini), shortly predating multiple Optix-driven wing pattern radiations in the speciose Heliconius and Eueides genera. We also characterize novel Optix and Doublesex expression in the male-specific pheromone wing scales of the basal heliconiines Dryas and Agraulis, thus illustrating that within the Heliconinii lineage, Optix has been evolutionarily redeployed in multiple contexts in association with diverse wing features.

Conclusions

Our findings reveal that the repeated co-option of Optix into various aspects of wing scale specification was associated with multiple evolutionary novelties over a relatively short evolutionary time scale. In particular, the recruitment of Optix expression in colored scale cell precursors was a necessary condition to the explosive diversification of passion-vine butterfly wing patterns. The novel deployment of a gene followed by spatial modulation of its expression in a given cell type could be a common mode of developmental innovation for triggering phenotypic radiations.

Eye-specification genes in the bacterial light organ of the bobtail squid Euprymna scolopes, and their expression in response to symbiont cues

The squid Euprymna scolopes has evolved independent sets of tissues capable of light detection, including a complex eye and a photophore or 'light organ', which houses the luminous bacterial symbiont Vibrio fischeri. As the eye and light organ originate from different embryonic tissues, we examined whether the eye-specification genes, pax6, eya, six, and dac, are shared by these two organs, and if so, whether they are regulated in the light organ by symbiosis. We obtained sequences of the four genes with PCR, confirmed orthology with phylogenetic analysis, and determined that each was expressed in the eye and light organ. With in situ hybridization (ISH), we localized the gene transcripts in developing embryos, comparing the patterns of expression in the two organs. The four transcripts localized to similar tissues, including those associated with the visual system ∼1/4 into embryogenesis (Naef stage 18) and the light organ ∼3/4 into embryogenesis (Naef stage 26). We used ISH and quantitative real-time PCR to examine transcript expression and differential regulation in postembryonic light organs in response to the following colonization conditions: wild-type, luminescent V. fischeri; a mutant strain defective in light production; and as a control, no symbiont. In ISH experiments light organs showed down regulation of the pax6, eya, and six transcripts in response to wild-type V. fischeri. Mutant strains also induced down regulation of the pax6 and eya transcripts, but not of the six transcript. Thus, luminescence was required for down regulation of the six transcript. We discuss these results in the context of symbiont-induced light-organ development. Our study indicates that the eye-specification genes are expressed in light-interacting tissues independent of their embryonic origin and are capable of responding to bacterial cues. These results offer evidence for evolutionary tinkering or the recruitment of eye development genes for use in a light-sensing photophore.

A quantitative atlas of Even-skipped and Hunchback expression in Clogmia albipunctata (Diptera: Psychodidae) blastoderm embryos

Background

Comparative studies of developmental processes are one of the main approaches to evolutionary developmental biology (evo-devo). Over recent years, there has been a shift of focus from the comparative study of particular regulatory genes to the level of whole gene networks. Reverse-engineering methods can be used to computationally reconstitute and analyze the function and dynamics of such networks. These methods require quantitative spatio-temporal expression data for model fitting. Obtaining such data in non-model organisms remains a major technical challenge, impeding the wider application of data-driven mathematical modeling to evo-devo.

Results

We have raised antibodies against four segmentation gene products in the moth midge Clogmia albipunctata, a non-drosophilid dipteran species. We have used these antibodies to create a quantitative atlas of protein expression patterns for the gap gene hunchback (hb), and the pair-rule gene even-skipped (eve). Our data reveal differences in the dynamics of Hb boundary positioning and Eve stripe formation between C. albipunctata and Drosophila melanogaster. Despite these differences, the overall relative spatial arrangement of Hb and Eve domains is remarkably conserved between these two distantly related dipteran species.

Conclusions

We provide a proof of principle that it is possible to acquire quantitative gene expression data at high accuracy and spatio-temporal resolution in non-model organisms. Our quantitative data extend earlier qualitative studies of segmentation gene expression in C. albipunctata, and provide a starting point for comparative reverse-engineering studies of the evolutionary and developmental dynamics of the segmentation gene system.

Experimental approaches for gene regulatory network construction: the chick as a model system

Setting up the body plan during embryonic development requires the coordinated action of many signals and transcriptional regulators in a precise temporal sequence and spatial pattern. The last decades have seen an explosion of information describing the molecular control of many developmental processes. The next challenge is to integrate this information into logic "wiring diagrams" that visualize gene actions and outputs, have predictive power and point to key control nodes. Here, we provide an experimental workflow on how to construct gene regulatory networks using the chick as model system.

Unraveling the thread of nature's tapestry: the genetics of diversity and convergence in animal pigmentation

Animals display incredibly diverse color patterns yet little is known about the underlying genetic basis of these phenotypes. However, emerging results are reshaping our view of how the process of phenotypic evolution occurs. Here, we outline recent research from three particularly active areas of investigation: melanin pigmentation in Drosophila, wing patterning in butterflies, and pigment variation in lizards. For each system, we highlight (i) the function and evolution of color variation, (ii) various approaches that have been used to explore the genetic basis of pigment variation, and (iii) conclusions regarding the genetic basis of convergent evolution which have emerged from comparative analyses. Results from these studies indicate that natural variation in pigmentation is a particularly powerful tool to examine the molecular basis of evolution, especially with regard to convergent or parallel evolution. Comparison of these systems also reveals that the molecular basis of convergent evolution is heterogeneous, sometimes involving conserved mechanisms and sometimes not. In the near future, additional work in other emerging systems will substantially expand the scope of available comparisons.

A single origin for nymphalid butterfly eyespots followed by widespread loss of associated gene expression

Understanding how novel complex traits originate involves investigating the time of origin of the trait, as well as the origin of its underlying gene regulatory network in a broad comparative phylogenetic framework. The eyespot of nymphalid butterflies has served as an example of a novel complex trait, as multiple genes are expressed during eyespot development. Yet the origins of eyespots remain unknown. Using a dataset of more than 400 images of butterflies with a known phylogeny and gene expression data for five eyespot-associated genes from over twenty species, we tested origin hypotheses for both eyespots and eyespot-associated genes. We show that eyespots evolved once within the family Nymphalidae, approximately 90 million years ago, concurrent with expression of at least three genes associated with early eyespot development. We also show multiple losses of expression of most genes from this early three-gene cluster, without corresponding losses of eyespots. We propose that complex traits, such as eyespots, may have originated via co-option of a large pre-existing complex gene regulatory network that was subsequently streamlined of genes not required to fulfill its novel developmental function.

Butterfly eyespots play an essential role in natural and sexual selection, yet the evolutionary origins of eyespots and of their underlying gene regulatory network remain unknown. By scoring phenotypes and wing expression of five genes in 399 and 21 nymphalid species, respectively, we tested when eyespots and expression of their associated genes evolved. We found that the origin of eyespots was concurrent with the origin of the gene expression patterns, approximately 90 million years ago. Following this event, many genes expressed in eyespot development were lost in some lineages without a corresponding loss of eyespots, indicating substantial evolution in the cluster of genes associated with eyespots. This finding suggests that complex traits such as butterfly eyespots may initially evolve by re-deploying pre-existing gene regulatory networks, which are subsequently trimmed of genes that are unnecessary in the novel context.

Transcriptome analysis reveals novel patterning and pigmentation genes underlying Heliconius butterfly wing pattern variation

Background

Heliconius butterfly wing pattern diversity offers a unique opportunity to investigate how natural genetic variation can drive the evolution of complex adaptive phenotypes. Positional cloning and candidate gene studies have identified a handful of regulatory and pigmentation genes implicated in Heliconius wing pattern variation, but little is known about the greater developmental networks within which these genes interact to pattern a wing. Here we took a large-scale transcriptomic approach to identify the network of genes involved in Heliconius wing pattern development and variation. This included applying over 140 transcriptome microarrays to assay gene expression in dissected wing pattern elements across a range of developmental stages and wing pattern morphs of Heliconius erato.

Results

We identified a number of putative early prepattern genes with color-pattern related expression domains. We also identified 51 genes differentially expressed in association with natural color pattern variation. Of these, the previously identified color pattern “switch gene” optix was recovered as the first transcript to show color-specific differential expression. Most differentially expressed genes were transcribed late in pupal development and have roles in cuticle formation or pigment synthesis. These include previously undescribed transporter genes associated with ommochrome pigmentation. Furthermore, we observed upregulation of melanin-repressing genes such as ebony and Dat1 in non-melanic patterns.

Conclusions

This study identifies many new genes implicated in butterfly wing pattern development and provides a glimpse into the number and types of genes affected by variation in genes that drive color pattern evolution.