Molecular evolution constraints in the floral organ specification gene regulatory network module across 18 angiosperm genomes

Comprehensive genome-wide identification and functional characterization of mapk gene family in northern snakeheads (Channa argus)

The mitogen-activated protein kinase (MAPK) signaling cascade, integral to cellular regulation, orchestrates cell growth, differentiation, stress response, and inflammatory reactions to adapt to challenging environments. The northern snakeheads (Channa argus), a valuable freshwater species known for its hypoxia tolerance, rapid growth, and high nutritional value, lacks comprehensive research on its mapk gene family. In this study, we identified 16 mapk genes in northern snakeheads, among which mapk8, mapk12 and mapk14 have duplicate copies. Phylogenetic analysis confirmed the evolutionary conservation of this gene family. Structural and motif analyses further underscored the conserved nature of these genes. Expression pattern analysis under abiotic and biotic stress conditions showed significant differences expression of mapks in the gills and suprabranchial organ (SBO) after air exposure, as well as in the brain following cold stress, highlighting the extensive role of mapks in stress regulation. It was worth noting that the significant expression differences of mapks were also observed in the spleen after N. seriolae infection, implicating that these genes may be involved in the regulation of innate immune responses. Additionally, analysis of protein-protein interaction (PPI) networks suggested that the co-activation of multiple MAPK signaling pathways may play a key role in regulating an organism's response to biotic and abiotic stresses. This study provides a detailed description of the mapk gene family in the northern snakeheads and elucidates its biological functions under various stress conditions, offering valuable insights into the regulatory mechanisms of the mapk gene family.

The Evolution of Hemocyanin Genes in Caenogastropoda: Gene Duplications and Intron Accumulation in Highly Diverse Gastropods

Hemocyanin is the oxygen transport protein of most molluscs and represents an important physiological factor that has to be well-adapted to their environments because of the strong influences of abiotic factors on its oxygen affinity. Multiple independent gene duplications and intron gains have been reported for hemocyanin genes of Tectipleura (Heterobranchia) and the caenogastropod species Pomacea canaliculata, which contrast with the uniform gene architectures of hemocyanins in Vetigastropoda. The goal of this study was to analyze hemocyanin gene evolution within the diverse group of Caenogastropoda in more detail. Our findings reveal multiple gene duplications and intron gains and imply that these represent general features of Apogastropoda hemocyanins. Whereas hemocyanin exon–intron structures are identical within different Tectipleura lineages, they differ strongly within Caenogastropoda among phylogenetic groups as well as between paralogous hemocyanin genes of the same species. Thus, intron accumulation took place more gradually within Caenogastropoda but finally led to a similar consequence, namely, a multitude of introns. Since both phenomena occurred independently within Heterobranchia and Caenogastropoda, the results support the hypothesis that introns may contribute to adaptive radiation by offering new opportunities for genetic variability (multiple paralogs that may evolve differently) and regulation (multiple introns). Our study indicates that adaptation of hemocyanin genes may be one of several factors that contributed to the evolution of the large diversity of Apogastropoda. While questions remain, this hypothesis is presented as a starting point for the further study of hemocyanin genes and possible correlations between hemocyanin diversity and adaptive radiation.

Variation in Pleiotropic Hub Gene Expression Is Associated with Interspecific Differences in Head Shape and Eye Size in Drosophila

Revealing the mechanisms underlying the breathtaking morphological diversity observed in nature is a major challenge in Biology. It has been established that recurrent mutations in hotspot genes cause the repeated evolution of morphological traits, such as body pigmentation or the gain and loss of structures. To date, however, it remains elusive whether hotspot genes contribute to natural variation in the size and shape of organs. As natural variation in head morphology is pervasive in Drosophila, we studied the molecular and developmental basis of differences in compound eye size and head shape in two closely related Drosophila species. We show differences in the progression of retinal differentiation between species and we applied comparative transcriptomics and chromatin accessibility data to identify the GATA transcription factor Pannier (Pnr) as central factor associated with these differences. Although the genetic manipulation of Pnr affected multiple aspects of dorsal head development, the effect of natural variation is restricted to a subset of the phenotypic space. We present data suggesting that this developmental constraint is caused by the coevolution of expression of pnr and its cofactor u-shaped (ush). We propose that natural variation in expression or function of highly connected developmental regulators with pleiotropic functions is a major driver for morphological evolution and we discuss implications on gene regulatory network evolution. In comparison to previous findings, our data strongly suggest that evolutionary hotspots are not the only contributors to the repeated evolution of eye size and head shape in Drosophila.

Transcriptome Analysis Reveals Putative Target Genes of APETALA3-3 During Early Floral Development in Nigella damascena L

Even though petals are homoplastic structures, their identity consistently involves genes of the APETALA3 (AP3) lineage. However, the extent to which the networks downstream of AP3 are conserved in species with petals of different evolutionary origins is unknown. In Ranunculaceae, the specificity of the AP3-III lineage offers a great opportunity to identify the petal gene regulatory network in a comparative framework. Using a transcriptomic approach, we investigated putative target genes of the AP3-III ortholog NdAP3-3 in Nigella damascena at early developmental stages when petal identity is determined, and we compared our data with that from selected eudicot species. We generated a de novo reference transcriptome to carry out a differential gene expression analysis between the wild-type and mutant NdAP3-3 genotypes differing by the presence vs. absence of petals at early stages of floral development. Among the 1,620 genes that were significantly differentially expressed between the two genotypes, functional annotation suggested a large involvement of nuclear activities, including regulation of transcription, and enrichment in processes linked to cell proliferation. Comparing with Arabidopsis data, we found that highly conserved genes between the two species are enriched in homologs of direct targets of the AtAP3 protein. Integrating AP3-3 binding site data from another Ranunculaceae species, Aquilegia coerulea, allowed us to identify a set of 18 putative target genes that were conserved between the three species. Our results suggest that, despite the independent evolutionary origin of petals in core eudicots and Ranunculaceae, a small conserved set of genes determines petal identity and early development in these taxa.

Comparative analysis of MAPK and MKK gene families reveals differential evolutionary patterns in Brachypodium distachyon inbred lines

BACKGROUND

Mitogen-activated protein kinase (MAPK) cascades are involved with signal transduction in almost every aspect of plant growth and development, as well as biotic and abiotic stress responses. The evolutionary analysis of MAPKs and MKKs in individual or entire plant species has been reported, but the evolutionary patterns in the diverse inbred lines of Brachypodium distachyon are still unclear.

RESULTS

We conducted the systematical molecular evolutionary analysis of B. distachyon. A total of 799 MAPKs and 618 MKKs were identified from 53 B. distachyon inbred lines. Remarkably, only three inbred lines had 16 MPKs and most of those inbred lines lacked MPK7-2 members, whereas 12 MKKs existed in almost all B. distachyon inbred lines. Phylogenetic analysis indicated that MAPKs and MKKs were divided into four groups as previously reported, grouping them in the same branch as corresponding members. MPK21-2 was the exception and fell into two groups, which may be due to their exon-intron patterns, especially the untranslated regions (UTRs). We also found that differential evolution patterns of MKK10 paralogues from ancient tandem duplicates may have undergone functional divergence. Expression analyses suggested that MAPKs and MKKs likely played different roles in different genetic contexts within various tissues and with abiotic stresses.

CONCLUSION

Our study revealed that UTRs affected the structure and evolution of MPK21-2 genes and the differential evolution of MKK10 paralogues with ancient tandem duplication might have functional divergences. Our findings provide new insights into the functional evolution of genes in closely inbred lines.

Unraveling the Developmental and Genetic Mechanisms Underpinning Floral Architecture in Proteaceae

Proteaceae are a basal eudicot family with a highly conserved floral groundplan but which displays considerable variation in other aspects of floral and inflorescence morphology. Their morphological diversity and phylogenetic position make them good candidates for understanding the evolution of floral architecture, in particular the question of the homology of the undifferentiated perianth with the differentiated perianth of core eudicots, and the mechanisms underlying the repeated evolution of zygomorphy. In this paper, we combine a morphological approach to explore floral ontogenesis and a transcriptomic approach to access the genes involved in floral organ identity and development, focusing on Grevillea juniperina, a species from subfamily Grevilleoideae. We present developmental data for Grevillea juniperina and three additional species that differ in their floral symmetry using stereomicroscopy, SEM and High Resolution X-Ray Computed Tomography. We find that the adnation of stamens to tepals takes place at early developmental stages, and that the establishment of bilateral symmetry coincides with the asymmetrical growth of the single carpel. To set a framework for understanding the genetic basis of floral development in Proteaceae, we generated and annotated de novo a reference leaf/flower transcriptome from Grevillea juniperina. We found Grevillea homologs of all lineages of MADS-box genes involved in floral organ identity. Using Arabidopsis thaliana gene expression data as a reference, we found homologs of other genes involved in floral development in the transcriptome of G. juniperina. We also found at least 21 class I and class II TCP genes, a gene family involved in the regulation of growth processes, including floral symmetry. The expression patterns of a set of floral genes obtained from the transcriptome were characterized during floral development to assess their organ specificity and asymmetry of expression.

Comparative analysis of whole flower transcriptomes in the Zingiberales

The advancement of next generation sequencing technologies (NGS) has revolutionized our ability to generate large quantities of data at a genomic scale. Despite great challenges, these new sequencing technologies have empowered scientists to explore various relevant biological questions on non-model organisms, even in the absence of a complete sequenced reference genome. Here, we analyzed whole flower transcriptome libraries from exemplar species across the monocot order Zingiberales, using a comparative approach in order to gain insight into the evolution of the molecular mechanisms underlying flower development in the group. We identified 4,153 coding genes shared by all floral transcriptomes analyzed, and 1,748 genes that are only retrieved in the Zingiberales. We also identified 666 genes that are unique to the ginger lineage, and 2,001 that are only found in the banana group, while in the outgroup species Dichorisandra thyrsiflora J.C. Mikan (Commelinaceae) we retrieved 2,686 unique genes. It is possible that some of these genes underlie lineage-specific molecular mechanisms of floral diversification. We further discuss the nature of these lineage-specific datasets, emphasizing conserved and unique molecular processes with special emphasis in the Zingiberales. We also briefly discuss the strengths and shortcomings of de novo assembly for the study of developmental processes across divergent taxa from a particular order. Although this comparison is based exclusively on coding genes, with particular emphasis in transcription factors, we believe that the careful study of other regulatory mechanisms, such as non-coding RNAs, might reveal new levels of complexity, which were not explored in this work.

Functional divergence and intron variability during evolution of angiosperm TERMINAL FLOWER1 (TFL1) genes

The protein encoded by the TERMINAL FLOWER1 (TFL1) gene maintains indeterminacy in inflorescence meristem to repress flowering, and has undergone multiple duplications. However, basal angiosperms have one copy of a TFL1-like gene, which clusters with eudicot TFL1/CEN paralogs. Functional conservation has been reported in the paralogs CENTRORADIALIS (CEN) in eudicots, and ROOTS CURL IN NPA (RCNs) genes in monocots. In this study, long-term functional conservation and selective constraints were found between angiosperms, while the relaxation of selective constraints led to subfunctionalisation between paralogs. Long intron lengths of magnoliid TFL1-like gene contain more conserved motifs that potentially regulate TFL1/CEN/RCNs expression. These might be relevant to the functional flexibility of the non-duplicate TFL1-like gene in the basal angiosperms in comparison with the short, lower frequency intron lengths in eudicot and monocot TFL1/CEN/RCNs paralogs. The functionally conserved duplicates of eudicots and monocots evolved according to the duplication-degeneration-complementation model, avoiding redundancy by relaxation of selective constraints on exon 1 and exon 4. These data suggest that strong purifying selection has maintained the relevant functions of TFL1/CEN/RCNs paralogs on flowering regulation throughout the evolution of angiosperms, and the shorter introns with radical amino acid changes are important for the retention of paralogous duplicates.

Identification of Two New Mechanisms That Regulate Fruit Growth by Cell Expansion in Tomato

Key mechanisms controlling fruit weight and shape at the levels of meristem, ovary or very young fruit have already been identified using natural tomato diversity. We reasoned that new developmental modules prominent at later stages of fruit growth could be discovered by using new genetic and phenotypic diversity generated by saturated mutagenesis. Twelve fruit weight and tissue morphology mutants likely affected in late fruit growth were selected among thousands of fruit size and shape EMS mutants available in our tomato EMS mutant collection. Their thorough characterization at organ, tissue and cellular levels revealed two major clusters controlling fruit growth and tissue morphogenesis either through (i) the growth of all fruit tissues through isotropic cell expansion or (ii) only the growth of the pericarp through anisotropic cell expansion. These likely correspond to new cell expansion modules controlling fruit growth and tissue morphogenesis in tomato. Our study therefore opens the way for the identification of new gene regulatory networks controlling tomato fruit growth and morphology.

Change of Fate and Staminodial Laminarity as Potential Agents of Floral Diversification in the Zingiberales

The evolution of floral morphology in the monocot order Zingiberales shows a trend in which androecial whorl organs are progressively modified into variously conspicuous "petaloid" structures with differing degrees of fertility. Petaloidy of androecial members results from extensive laminarization of an otherwise radially symmetric structure. The genetic basis of the laminarization of androecial members has been addressed through recent candidate gene studies focused on understanding the spatiotemporal expression patterns of genes known to be necessary to floral organ formation. Here, we explore the correlation between gene duplication events and floral and inflorescence morphological diversification across the Zingiberales by inferring ancestral character states and gene copy number using the most widely accepted phylogenetic hypotheses. Our results suggest that the duplication and differential loss of GLOBOSA (GLO) copies is correlated with a change in the degree of the laminarization of androecial members. We also find an association with increased diversification in most families. We hypothesize that retention of paralogs in flower development genes could have led to a developmental shift affecting androecial organs with potential adaptive consequences, thus favoring diversification in some lineages but not others.

A Century of Sex Determination in Flowering Plants

Plants have evolved a diverse array of strategies for sexual reproduction, particularly through the modification of male and female organs at distinct points in development. The immense variation in sexual systems across the land plants provides a unique opportunity to study the genetic, epigenetic, phylogenetic, and ecological underpinnings of sex determination. Here, we reflect on more than a century of research into flowering plant sex determination, placing a particular focus on the foundational genetic and cytogenetic observations, experiments, and hypotheses. Building on the seminal work on the genetics of plant sex, modern comparative genomic analyses now allow us to address longstanding questions about sex determination and the origins of sex chromosomes.

The MADS-box XAANTAL1 increases proliferation at the Arabidopsis root stem-cell niche and participates in transition to differentiation by regulating cell-cycle components

Background Morphogenesis depends on the concerted modulation of cell proliferation and differentiation. Such modulation is dynamically adjusted in response to various external and internal signals via complex transcriptional regulatory networks that mediate between such signals and regulation of cell-cycle and cellular responses (proliferation, growth, differentiation). In plants, which are sessile, the proliferation/differentiation balance is plastically adjusted during their life cycle and transcriptional networks are important in this process. MADS-box genes are key developmental regulators in eukaryotes, but their role in cell proliferation and differentiation modulation in plants remains poorly studied. Methods We characterize the XAL1 loss-of-function xal1-2 allele and overexpression lines using quantitative cellular and cytometry analyses to explore its role in cell cycle, proliferation, stem-cell patterning and transition to differentiation. We used quantitative PCR and cellular markers to explore if XAL1 regulates cell-cycle components and PLETHORA1 (PLT1) gene expression, as well as confocal microscopy to analyse stem-cell niche organization. Key Results We previously showed that XAANTAL1 (XAL1/AGL12) is necessary for Arabidopsis root development as a promoter of cell proliferation in the root apical meristem. Here, we demonstrate that XAL1 positively regulates the expression of PLT1 and important components of the cell cycle: CYCD3;1, CYCA2;3, CYCB1;1, CDKB1;1 and CDT1a. In addition, we show that xal1-2 mutant plants have a premature transition to differentiation with root hairs appearing closer to the root tip, while endoreplication in these plants is partially compromised. Coincidently, the final size of cortex cells in the mutant is shorter than wild-type cells. Finally, XAL1 overexpression-lines corroborate that this transcription factor is able to promote cell proliferation at the stem-cell niche. Conclusion XAL1 seems to be an important component of the networks that modulate cell proliferation/differentiation transition and stem-cell proliferation during Arabidopsis root development; it also regulates several cell-cycle components.

Evolutionary Developmental Biology (Evo-Devo) Research in Latin America

Famous for its blind cavefish and Darwin's finches, Latin America is home to some of the richest biodiversity hotspots of our planet. The Latin American fauna and flora inspired and captivated naturalists from the nineteenth and twentieth centuries, including such notable pioneers such as Fritz Müller, Florentino Ameghino, and Léon Croizat who made a significant contribution to the study of embryology and evolutionary thinking. But, what are the historical and present contributions of the Latin American scientific community to Evo-Devo? Here, we provide the first comprehensive overview of the Evo-Devo laboratories based in Latin America and describe current lines of research based on endemic species, focusing on body plans and patterning, systematics, physiology, computational modeling approaches, ecology, and domestication. Literature searches reveal that Evo-Devo in Latin America is still in its early days; while showing encouraging indicators of productivity, it has not stabilized yet, because it relies on few and sparsely distributed laboratories. Coping with the rapid changes in national scientific policies and contributing to solve social and health issues specific to each region are among the main challenges faced by Latin American researchers. The 2015 inaugural meeting of the Pan-American Society for Evolutionary Developmental Biology played a pivotal role in bringing together Latin American researchers eager to initiate and consolidate regional and worldwide collaborative networks. Such networks will undoubtedly advance research on the extremely high genetic and phenotypic biodiversity of Latin America, bound to be an almost infinite source of amazement and fascinating findings for the Evo-Devo community.

Dynamic network modelling to understand flowering transition and floral patterning

Differentiation and morphogenetic processes during plant development are particularly robust. At the cellular level, however, plants also show great plasticity in response to environmental conditions, and can even reverse apparently terminal differentiated states with remarkable ease. Can we understand and predict both robust and plastic systemic responses as a general consequence of the non-trivial interplay between intracellular regulatory networks, extrinsic environmental signalling, and tissue-level mechanical constraints? Flower development has become an ideal model system to study these general questions of developmental biology, which are especially relevant to understanding stem cell patterning in plants, animals, and human disease. Decades of detailed study of molecular developmental genetics, as well as novel experimental techniques for in vivo assays in both wild-type and mutant plants, enable the postulation and testing of experimentally grounded mathematical and computational network dynamical models. Research in our group aims to explain the emergence of robust transitions that occur at the shoot apical meristem, as well as flower development, as the result of the collective action of key molecular components in regulatory networks subjected to intra-organismal signalling and extracellular constraints. Here we present a brief overview of recent work from our group, and that of others, focusing on the use of simple dynamical models to address cell-fate specification and cell-state stochastic dynamics during flowering transition and cell-state transitions at the shoot apical meristem of Arabidopsis thaliana. We also focus on how our work fits within the general field of plant developmental modelling, which is being developed by many others.

The impact of Polycomb group (PcG) and Trithorax group (TrxG) epigenetic factors in plant plasticity

Current advances indicate that epigenetic mechanisms play important roles in the regulatory networks involved in plant developmental responses to environmental conditions. Hence, understanding the role of such components becomes crucial to understanding the mechanisms underlying the plasticity and variability of plant traits, and thus the ecology and evolution of plant development. We now know that important components of phenotypic variation may result from heritable and reversible epigenetic mechanisms without genetic alterations. The epigenetic factors Polycomb group (PcG) and Trithorax group (TrxG) are involved in developmental processes that respond to environmental signals, playing important roles in plant plasticity. In this review, we discuss current knowledge of TrxG and PcG functions in different developmental processes in response to internal and environmental cues and we also integrate the emerging evidence concerning their function in plant plasticity. Many such plastic responses rely on meristematic cell behavior, including stem cell niche maintenance, cellular reprogramming, flowering and dormancy as well as stress memory. This information will help to determine how to integrate the role of epigenetic regulation into models of gene regulatory networks, which have mostly included transcriptional interactions underlying various aspects of plant development and its plastic response to environmental conditions.

Molecular Evolution of the Neural Crest Regulatory Network in Ray-Finned Fish

Gene regulatory networks (GRN) are central to developmental processes. They are composed of transcription factors and signaling molecules orchestrating gene expression modules that tightly regulate the development of organisms. The neural crest (NC) is a multipotent cell population that is considered a key innovation of vertebrates. Its derivatives contribute to shaping the astounding morphological diversity of jaws, teeth, head skeleton, or pigmentation. Here, we study the molecular evolution of the NC GRN by analyzing patterns of molecular divergence for a total of 36 genes in 16 species of bony fishes. Analyses of nonsynonymous to synonymous substitution rate ratios (dN/dS) support patterns of variable selective pressures among genes deployed at different stages of NC development, consistent with the developmental hourglass model. Model-based clustering techniques of sequence features support the notion of extreme conservation of NC-genes across the entire network. Our data show that most genes are under strong purifying selection that is maintained throughout ray-finned fish evolution. Late NC development genes reveal a pattern of increased constraints in more recent lineages. Additionally, seven of the NC-genes showed signs of relaxation of purifying selection in the famously species-rich lineage of cichlid fishes. This suggests that NC genes might have played a role in the adaptive radiation of cichlids by granting flexibility in the development of NC-derived traits-suggesting an important role for NC network architecture during the diversification in vertebrates.

Reshaping the epigenetic landscape during early flower development: induction of attractor transitions by relative differences in gene decay rates

BACKGROUND

Gene regulatory network (GRN) dynamical models are standard systems biology tools for the mechanistic understanding of developmental processes and are enabling the formalization of the epigenetic landscape (EL) model.

METHODS

In this work we propose a modeling framework which integrates standard mathematical analyses to extend the simple GRN Boolean model in order to address questions regarding the impact of gene specific perturbations in cell-fate decisions during development.

RESULTS

We systematically tested the propensity of individual genes to produce qualitative changes to the EL induced by modification of gene characteristic decay rates reflecting the temporal dynamics of differentiation stimuli. By applying this approach to the flower specification GRN (FOS-GRN) we uncovered differences in the functional (dynamical) role of their genes. The observed dynamical behavior correlates with biological observables. We found a relationship between the propensity of undergoing attractor transitions between attraction basins in the EL and the direction of differentiation during early flower development - being less likely to induce up-stream attractor transitions as the course of development progresses. Our model also uncovered a potential mechanism at play during the transition from EL basins defining inflorescence meristem to those associated to flower organs meristem. Additionally, our analysis provided a mechanistic interpretation of the homeotic property of the ABC genes, being more likely to produce both an induced inter-attractor transition and to specify a novel attractor. Finally, we found that there is a close relationship between a gene's topological features and its propensity to produce attractor transitions.

CONCLUSIONS

The study of how the state-space associated with a dynamical model of a GRN can be restructured by modulation of genes' characteristic expression times is an important aid for understanding underlying mechanisms occurring during development. Our contribution offers a simple framework to approach such problem, as exemplified here by the case of flower development. Different GRN models and the effect of diverse inductive signals can be explored within the same framework. We speculate that the dynamical role of specific genes within a GRN, as uncovered here, might give information about which genes are more likely to link a module to other regulatory circuits and signaling transduction pathways.

Role of transcriptional regulation in the evolution of plant phenotype: A dynamic systems approach

A growing body of evidence suggests that alterations in transcriptional regulation of genes involved in modulating development are an important part of phenotypic evolution, and this can be documented among species and within populations. While the effects of differential transcriptional regulation in organismal development have been preferentially studied in animal systems, this phenomenon has also been addressed in plants. In this review, we summarize evidence for cis-regulatory mutations, trans-regulatory changes and epigenetic modifications as molecular events underlying important phenotypic alterations, and thus shaping the evolution of plant development. We postulate that a mechanistic understanding of why such molecular alterations have a key role in development, morphology and evolution will have to rely on dynamic models of complex regulatory networks that consider the concerted action of genetic and nongenetic components, and that also incorporate the restrictions underlying the genotype to phenotype mapping process. Developmental Dynamics 244:1074-1095, 2015. © 2015 Wiley Periodicals, Inc.

Pleiotropy constrains the evolution of protein but not regulatory sequences in a transcription regulatory network influencing complex social behaviors

It is increasingly apparent that genes and networks that influence complex behavior are evolutionary conserved, which is paradoxical considering that behavior is labile over evolutionary timescales. How does adaptive change in behavior arise if behavior is controlled by conserved, pleiotropic, and likely evolutionary constrained genes? Pleiotropy and connectedness are known to constrain the general rate of protein evolution, prompting some to suggest that the evolution of complex traits, including behavior, is fuelled by regulatory sequence evolution. However, we seldom have data on the strength of selection on mutations in coding and regulatory sequences, and this hinders our ability to study how pleiotropy influences coding and regulatory sequence evolution. Here we use population genomics to estimate the strength of selection on coding and regulatory mutations for a transcriptional regulatory network that influences complex behavior of honey bees. We found that replacement mutations in highly connected transcription factors and target genes experience significantly stronger negative selection relative to weakly connected transcription factors and targets. Adaptively evolving proteins were significantly more likely to reside at the periphery of the regulatory network, while proteins with signs of negative selection were near the core of the network. Interestingly, connectedness and network structure had minimal influence on the strength of selection on putative regulatory sequences for both transcription factors and their targets. Our study indicates that adaptive evolution of complex behavior can arise because of positive selection on protein-coding mutations in peripheral genes, and on regulatory sequence mutations in both transcription factors and their targets throughout the network.