Individual variation in preference behavior in sailfin fish refines the neurotranscriptomic pathway for mate preference.

Social interactions can drive distinct gene expression profiles which may vary by social context. Here we use female sailfin molly fish (Poecilia latipinna) to identify genomic profiles associated with preference behavior in distinct social contexts: male interactions (mate choice) versus female interactions (shoaling partner preference). We measured the behavior of 15 females interacting in a non-contact environment with either two males or two females for 30 min followed by whole-brain transcriptomic profiling by RNA sequencing. We profiled females that exhibited high levels of social affiliation and great variation in preference behavior to identify an order of magnitude more differentially expressed genes associated with behavioral variation than by differences in social context. Using a linear model (limma), we took advantage of the individual variation in preference behavior to identify unique gene sets that exhibited distinct correlational patterns of expression with preference behavior in each social context. By combining limma and weighted gene co-expression network analyses (WGCNA) approaches we identified a refined set of 401 genes robustly associated with mate preference that is independent of shoaling partner preference or general social affiliation. While our refined gene set confirmed neural plasticity pathways involvement in moderating female preference behavior, we also identified a significant proportion of discovered that our preference-associated genes were enriched for 'immune system' gene ontology categories. We hypothesize that the association between mate preference and transcriptomic immune function is driven by the less well-known role of these genes in neural plasticity which is likely involved in higher-order learning and processing during mate choice decisions.

Distinct immune and transcriptomic profiles in dominant versus subordinate males in mouse social hierarchies.

Social status is a critical factor determining health outcomes in human and nonhuman social species. In social hierarchies with reproductive skew, individuals compete to monopolize resources and increase mating opportunities. This can come at a significant energetic cost leading to trade-offs between different physiological systems. In particular, changes in energetic investment in the immune system can have significant short and long-term effects on fitness and health. We have previously found that dominant alpha male mice living in social hierarchies have increased metabolic demands related to territorial defense. In this study, we tested the hypothesis that high-ranking male mice favor adaptive immunity, while subordinate mice show higher investment in innate immunity. We housed 12 groups of 10 outbred CD-1 male mice in a social housing system. All formed linear social hierarchies and subordinate mice had higher concentrations of plasma corticosterone (CORT) than alpha males. This difference was heightened in highly despotic hierarchies. Using flow cytometry, we found that dominant status was associated with a significant shift in immunophenotypes towards favoring adaptive versus innate immunity. Using Tag-Seq to profile hepatic and splenic transcriptomes of alpha and subordinate males, we identified genes that regulate metabolic and immune defense pathways that are associated with status and/or CORT concentration. In the liver, dominant animals showed a relatively higher expression of specific genes involved in major urinary production and catabolic processes, whereas subordinate animals showed relatively higher expression of genes promoting biosynthetic processes, wound healing, and proinflammatory responses. In spleen, subordinate mice showed relatively higher expression of genes facilitating oxidative phosphorylation and DNA repair and CORT was negatively associated with genes involved in lymphocyte proliferation and activation. Together, our findings suggest that dominant and subordinate animals adaptively shift immune profiles and peripheral gene expression to match their contextual needs.

Comparative Transcriptomics Reveals Distinct Patterns of Gene Expression Conservation through Vertebrate Embryogenesis.

Despite life's diversity, studies of variation often remind us of our shared evolutionary past. Abundant genome sequencing and analyses of gene regulatory networks illustrate that genes and entire pathways are conserved, reused, and elaborated in the evolution of diversity. Predating these discoveries, 19th-century embryologists observed that though morphology at birth varies tremendously, certain stages of vertebrate embryogenesis appear remarkably similar across vertebrates. In the mid to late 20th century, anatomical variability of early and late-stage embryos and conservation of mid-stages embryos (the "phylotypic" stage) was named the hourglass model of diversification. This model has found mixed support in recent analyses comparing gene expression across species possibly owing to differences in species, embryonic stages, and gene sets compared. We compare 186 microarray and RNA-seq data sets covering embryogenesis in six vertebrate species. We use an unbiased clustering approach to group stages of embryogenesis by transcriptomic similarity and ask whether gene expression similarity of clustered embryonic stages deviates from a null expectation. We characterize expression conservation patterns of each gene at each evolutionary node after correcting for phylogenetic nonindependence. We find significant enrichment of genes exhibiting early conservation, hourglass, late conservation patterns in both microarray and RNA-seq data sets. Enrichment of genes showing patterned conservation through embryogenesis indicates diversification of embryogenesis may be temporally constrained. However, the circumstances under which each pattern emerges remain unknown and require both broad evolutionary sampling and systematic examination of embryogenesis across species.

EDCs Reorganize Brain-Behavior Phenotypic Relationships in Rats.

All species, including humans, are exposed to endocrine-disrupting chemicals (EDCs). Previous experiments have shown behavioral deficits caused by EDCs that have implications for social competence and sexual selection. The neuromolecular mechanisms for these behavioral changes induced by EDCs have not been thoroughly explored. Here, we tested the hypothesis that EDCs administered to rats during a critical period of embryonic brain development would lead to the disruption of normal social preference behavior, and that this involves a network of underlying gene pathways in brain regions that regulate these behaviors. Rats were exposed prenatally to human-relevant concentrations of EDCs (polychlorinated biphenyls [PCBs], vinclozolin [VIN]), or vehicle. In adulthood, a sociosexual preference test was administered. We profiled gene expression of in preoptic area, medial amygdala, and ventromedial nucleus. Prenatal PCBs impaired sociosexual preference in both sexes, and VIN disrupted this behavior in males. Each brain region had unique sets of genes altered in a sex- and EDC-specific manner. The effects of EDCs on individual traits were typically small, but robust; EDC exposure changed the relationships between gene expression and behavior, a pattern we refer to as dis-integration and reconstitution. These findings underscore the effects that developmental exposure to EDCs can have on adult social behavior, highlight sex-specific and individual variation in responses, and provide a foundation for further work on the disruption of genes and behavior after prenatal exposure to EDCs.

Brain transcriptomics of agonistic behaviour in the weakly electric fish Gymnotus omarorum, a wild teleost model of non-breeding aggression.

Differences in social status are often mediated by agonistic encounters between competitors. Robust literature has examined social status-dependent brain gene expression profiles across vertebrates, yet social status and reproductive state are often confounded. It has therefore been challenging to identify the neuromolecular mechanisms underlying social status independent of reproductive state. Weakly electric fish, Gymnotus omarorum, display territorial aggression and social dominance independent of reproductive state. We use wild-derived G. omarorum males to conduct a transcriptomic analysis of non-breeding social dominance relationships. After allowing paired rivals to establish a dominance hierarchy, we profiled the transcriptomes of brain sections containing the preoptic area (region involved in regulating aggressive behaviour) in dominant and subordinate individuals. We identified 16 differentially expressed genes (FDR < 0.05) and numerous genes that co-varied with behavioural traits. We also compared our results with previous reports of differential gene expression in other teleost species. Overall, our study establishes G. omarorum as a powerful model system for understanding the neuromolecular bases of social status independent of reproductive state.

Conserved transcriptomic profiles underpin monogamy across vertebrates.

Social monogamy, typically characterized by the formation of a pair bond, increased territorial defense, and often biparental care, has independently evolved multiple times in animals. Despite the independent evolutionary origins of monogamous mating systems, several homologous brain regions and neuropeptides and their receptors have been shown to play a conserved role in regulating social affiliation and parental care, but little is known about the neuromolecular mechanisms underlying monogamy on a genomic scale. Here, we compare neural transcriptomes of reproductive males in monogamous and nonmonogamous species pairs of Peromyscus mice, Microtus voles, parid songbirds, dendrobatid frogs, and Xenotilapia species of cichlid fishes. We find that, while evolutionary divergence time between species or clades did not explain gene expression similarity, characteristics of the mating system correlated with neural gene expression patterns, and neural gene expression varied concordantly across vertebrates when species transition to monogamy. Our study provides evidence of a universal transcriptomic mechanism underlying the evolution of monogamy in vertebrates.

Linking conceptual mechanisms and transcriptomic evidence of plasticity-driven diversification.

The East African cichlid fishes provide text book examples of adaptive radiation. Diversification and speciation of cichlids associate with variation in diet and trophic morphologies among other ecological, behavioural and morphological phenotypes (Kocher 2004). Numerous case studies in cichlids reveal a role of developmental plasticity in generating jaw ecomorphs in response to variation in feeding ecology that can facilitate niche exploitation and subsequent diversification (e.g. Meyer 1987). Specifically, genetic divergence among such environmentally induced morphs can occur via reproductive isolation due to divergence in habitat and resource use in combination with genetic assimilation of environmentally induced phenotypes (West-Eberhard 2003; Pfennig et al. 2010). Expansion of this conceptual model has been hampered in part by the limited knowledge of the molecular mechanisms of plasticity in nonstandard model systems and the associated lack of evidence linking the molecular mechanisms of plasticity to those that generate phenotypic divergence among populations and taxa. In this issue of Molecular Ecology, Gunter et al. (2013) identify the transcriptional mechanisms of diet-induced lower pharyngeal jaw (LPJ) plasticity in the cichlid fish Astatoreochromis alluaudi. Natural populations of A. alluaudi exhibit variation in jaw morphology in relation to diet hardness. Among the plastic responses to diet are adjustments to the LPJ ranging from a robust molariform morph in response to a hard diet to a more gracile papilliform morph in response to a soft diet (Fig. 1). Gunter and colleagues induced developmental plasticity of the A. alluaudi jaw using diet manipulations and compared LPJ transcriptomic profiles of the resulting morphs. In this foundational work, the authors identify 187 differentially expressed genes that underlie the development and maintenance of diet-induced LPJ morphologies. This list includes a wide range of genes spanning from broad-acting transcription factors to signalling molecules and structural genes. Here, I examine the ontogeny of the molecular response to mechanical strain imposed by diet hardness and discuss the role of the stages of this response in the evolution of plasticity and plasticity-driven diversification.

Uterine gene expression in the live-bearing lizard, Chalcides ocellatus, reveals convergence of squamate reptile and mammalian pregnancy mechanisms.

Although the morphological and physiological changes involved in pregnancy in live-bearing reptiles are well studied, the genetic mechanisms that underlie these changes are not known. We used the viviparous African Ocellated Skink, Chalcides ocellatus, as a model to identify a near complete gene expression profile associated with pregnancy using RNA-Seq analyses of uterine transcriptomes. Pregnancy in C. ocellatus is associated with upregulation of uterine genes involved with metabolism, cell proliferation and death, and cellular transport. Moreover, there are clear parallels between the genetic processes associated with pregnancy in mammals and Chalcides in expression of genes related to tissue remodeling, angiogenesis, immune system regulation, and nutrient provisioning to the embryo. In particular, the pregnant uterine transcriptome is dominated by expression of proteolytic enzymes that we speculate are involved both with remodeling the chorioallantoic placenta and histotrophy in the omphaloplacenta. Elements of the maternal innate immune system are downregulated in the pregnant uterus, indicating a potential mechanism to avoid rejection of the embryo. We found a downregulation of major histocompatability complex loci and estrogen and progesterone receptors in the pregnant uterus. This pattern is similar to mammals but cannot be explained by the mammalian model. The latter finding provides evidence that pregnancy is controlled by different endocrinological mechanisms in mammals and reptiles. Finally, 88% of the identified genes are expressed in both the pregnant and the nonpregnant uterus, and thus, morphological and physiological changes associated with C. ocellatus pregnancy are likely a result of regulation of genes continually expressed in the uterus rather than the initiation of expression of unique genes.

Transcriptomic analysis of avian digits reveals conserved and derived digit identities in birds.

Morphological characters are the result of developmental gene expression. The identity of a character is ultimately grounded in the gene regulatory network directing development and thus whole-genome gene expression data can provide evidence about character identity. This approach has been successfully used to assess cell-type identity. Here we use transcriptomic data to address a long-standing uncertainty in evolutionary biology, the identity of avian wing digits. Embryological evidence clearly identifies the three wing digits as developing from digit positions 2, 3 and 4 (ref. 6), whereas palaeontological data suggest that they are digits I, II and III. We compare the transcriptomes of the wing and foot digits and find a strong signal that unites the first wing digit with the first foot digit, even though the first wing digit develops from embryological position 2. Interestingly, our transcriptomic data of the posterior digits show a higher degree of differentiation among forelimb digits compared with hindlimb digits. These data show that in the stem lineage of birds the first digit underwent a translocation from digit position 1 to position 2, and further indicate that the posterior wing digits have unique identities contrary to any model of avian digit identity proposed so far.

Why ontogenetic homology criteria can be misleading: lessons from digit identity transformations.

As the basis for comparative biology, correctly assigning character homology is critical. Yet, identifying homologous characters in practice is often challenging. Among the major roadblocks is that the mechanistic bases of character homology remain in question. Thus, investigators must rely on several independent lines of evidence (e.g., character anatomy, phylogenetic distribution, or embryological position); however, these distinct sources of evidence often lead to conflicting diagnoses of character homology. What is more, there is no consensus regarding the relative importance of distinct lines of evidence for determining character homology. Here, we review the difficulties that have hindered the search for the mechanistic bases of character identity, and relate these issues to a recently proposed mechanistic hypothesis of character identity--the Character Identity Network Hypothesis. Next, using two well-studied cases of homology conflict (i.e., avian and skink digit identity), we assess the utility of different lines of evidence in diagnosing homology. We conclude that, when comparing adult structures, because anatomical characters more closely reflect the actions of the developmental genetic mechanisms of character individuation they are more reliable than embryological homology criteria.

Identity of the avian wing digits: problems resolved and unsolved.

Controversy over bird wing digit identity has been a touchstone for various ideas in the phylogeny of birds, homology, and developmental evolution. This review summarizes the past 10 years of progress toward understanding avian digit identity. We conclude that the sum of evidence supports the Frame Shift Hypothesis, indicating that the avian wing digits have changed anatomical location. Briefly, the derivation of birds from theropod dinosaurs and the positional identities of the avian wing digits as 2, 3, and 4¹ are no longer in question. Additionally, increasing evidence indicates that the developmental programs for identity of the wing digits are of digits I, II, and III. Therefore, the attention moves from whether the digit identity frame shift occurred, to what the mechanisms of the frame shift were, and when in evolution it happened. There is considerable uncertainty about these issues and we identify exciting new research directions to resolve them.

Digital gene expression tag profiling of bat digits provides robust candidates contributing to wing formation.

BACKGROUND: As the only truly flying mammals, bats use their unique wing - consisting of four elongated digits (digits II-V) connected by membranes - to power their flight. In addition to the elongated digits II-V, the forelimb contains one shorter digit (digit I) that is morphologically similar to the hindlimb digits. Here, we capitalized on the morphological variation among the bat forelimb digits to investigate the molecular mechanisms underlying digit elongation and wing formation. Using next generation sequencing technology, we performed digital gene expression tag profiling (DGE-tag profiling) of developing digits in a pooled sample of two Myotis ricketti and validated our sequencing results using real-time quantitative PCR (RT-qPCR) of gene expression in the developing digits of two Hipposideros armiger. RESULTS: Among hundreds of genes exhibiting significant differences in expression between the short and long digits, we highlight 14 genes most related to digit elongation. These genes include two Tbx genes (Tbx3 and Tbx15), five BMP pathway genes (Bmp3, RGMB, Smad1, Smad4 and Nog), four Homeobox genes (Hoxd8, Hoxd9, Hoxa1 and Satb1), and three other genes (Twist1, Tmeff2 and Enpp2) related to digit malformations or cell proliferation. In addition, our results suggest that Tbx4 and Pitx2 contribute to the morphological similarity and five genes (Acta1, Tnnc2, Atp2a1, Hrc and Myoz1) contribute to the functional similarity between the thumb and hindlimb digits. CONCLUSIONS: Results of this study not only implicate many developmental genes as robust candidates underlying digit elongation and wing formation in bats, but also provide a better understanding of the genes involved in autopodial development in general.

Developmental plasticity links local adaptation and evolutionary diversification in foraging morphology.

Developmental plasticity is thought to reconcile the constraining role of natural selection in maintaining local adaptation with evolutionary diversification under novel conditions, but empirical documentations are rare. In vertebrates, growth and development of bones is partially guided by contractions of attached musculature and such muscle activity changes progressively through embryonic development from sporadic motility to direct functional effects. In species with short generation times, delayed skull maturation extends the guiding effects of muscle activity on formation of foraging morphology into adulthood, providing an opportunity to directly examine the links between plasticity of bone development, ecological adaptations, and evolutionary diversification in foraging morphology. In this case, the morphological consequences of inputs due to local functional requirements should be evident in adaptive divergence across taxa. Here we provide evidence that epigenetic regulation of bone growth in Soricid shrews may enable both development of local adaptations and evolutionary divergence in mandibular morphology. We contrast the effects of muscle stimulation on early- vs. late-maturing components of, foraging apparatus to show that the morphology of late-maturing components is more affected by functional requirements than are early-ossifying traits. Further, the divergence in foraging morphology across shrew species occurs along the directions delineated by inductive effects of muscle loading and bite force on bone formation in late-maturing but not early-maturing mandible components within species. These results support the hypothesis that developmental plasticity can link maintenance of local adaptations with evolutionary diversification in morphology.

Evolution of digit identity in the three-toed Italian skink Chalcides chalcides: a new case of digit identity frame shift.

Digit identity in the avian wing is a classical example of conflicting anatomical and embryological evidence regarding digit homology. Anatomical in conjunction with phylogenetic evidence supports the hypothesis that the three remaining digits in the bird wing are digits 1, 2, and 3. At the same time, various lines of embryological evidence support the notion that these digits develop in positions that normally produce digits 2, 3, and 4. In recent years, gene expression as well as experimental evidence was published that supports the hypothesis that this discrepancy arose from a digit identity shift in the evolution of the bird wing. A similar but less well-known controversy has been ongoing since the late 19th century regarding the identity of the digits of the three-toed Italian skink, Chalcides chalcides. Comparative anatomy identifies these digits as 1, 2, and 3, while embryological evidence suggests their derivation from embryological positions 2, 3, and 4. Here we re-examine this evidence and add gene expression data to determine the identity of the three digits of C. chalcides. The data confirm that the adult and the embryological evidence for digit identity are in conflict, and the expression of Hoxd11 suggests that digits 1, 2, and 3 develop in positions 2, 3, and 4. We conclude that in C. chalcides, and likely in its close relatives, a digit identity frame shift has occurred, similar to the one in avian evolution. This result suggests that changes in of digit identity might be a more frequent consequence of digit reduction than previously assumed.

Evolution on a local scale: developmental, functional, and genetic bases of divergence in bill form and associated changes in song structure between adjacent habitats.

Divergent selection on traits involved in both local adaptation and the production of mating signals can strongly facilitate population differentiation. Because of its links to foraging morphologies and cultural inheritance song of birds can contribute particularly strongly to maintenance of local adaptations. In two adjacent habitats--native Sonoran desert and urban areas--house finches (Carpodacus mexicanus) forage on seeds that are highly distinct in size and shell hardness and require different bite forces and bill morphologies. Here, we first document strong and habitat-specific natural selection on bill traits linked to bite force and find adaptive modifications of bite force and bill morphology and associated divergence in courtship song between the two habitats. Second, we investigate the developmental basis of this divergence and find that early ontogenetic tissue transformation in bill, but not skeletal traits, is accelerated in the urban population and that the mandibular primordia of the large-beaked urban finches express bone morphogenetic proteins (BMP) earlier and at higher level than those of the desert finches. Further, we show that despite being geographically adjacent, urban and desert populations are nevertheless genetically distinct corroborating findings of early developmental divergence between them. Taken together, these results suggest that divergent selection on function and development of traits involved in production of mating signals, in combination with localized learning of such signals, can be very effective at maintaining local adaptations, even at small spatial scales and in highly mobile animals.

Functional equivalence of morphologies enables morphological and ecological diversity.

Diversity in organismal forms among taxa is thought to reflect distinct selection pressures across environments. The central assumption underlying this expectation is that taxa experiencing similar selection have similar response to that selection. However, because selection acts on trait function, taxa similarity in selection response depends crucially on the relationship between function and morphology. Further, when a trait consists of multiple parts, changes in function in response to selection can result from modification of different parts, and adaptation to the same environment might result in functional but not morphological similarity. Here, we address the extent to which functional and morphological diversity in masticatory apparatus of soricid shrews reflects a shared ecological characteristic of their diet type. We examine the factors limiting morphological variation across shrew species by assessing the relative contribution of trait function (biomechanics of the jaw), ecology, and phylogeny to species similarity in mandibular traits. We found that species that shared diet type were functionally but not morphologically similar. The presence of multiple semi-independently varying traits enabled functional equivalence of composite foraging morphologies and resulted in variable response to selection exerted by similar diet. We show that functional equivalence of multiple morphologies enabled persistence of differences in habitat use (e.g., habitat moisture and coverage) among species that specialize on the same diet. We discuss the importance of developmental and functional integration among traits for evolutionary diversification of morphological structures that generate equivalent functions.

Evolution of ontogeny: linking epigenetic remodeling and genetic adaptation in skeletal structures.

Evolutionary diversifications are commonly attributed to the continued modifications of a conserved genetic toolkit of developmental pathways, such that complexity and convergence in organismal forms are assumed to be due to similarity in genetic mechanisms or environmental conditions. This approach, however, confounds the causes of organismal development with the causes of organismal differences and, as such, has only limited utility for addressing the cause of evolutionary change. Molecular mechanisms that are closely involved in both developmental response to environmental signals and major evolutionary innovations and diversifications are uniquely suited to bridge this gap by connecting explicitly the causes of within-generation variation with the causes of divergence of taxa. Developmental pathways of bone formation and a common role for bone morphogenetic proteins (BMPs) in both epigenetic bone remodeling and the evolution of major adaptive diversifications provide such opportunity. We show that variation in timing of ossification can result in similar phenotypic patterns through epigenetically induced changes in gene expression and propose that both genetic accommodation of environmentally induced developmental pathways and flexibility in development across environments evolve through heterochronic shifts in bone maturation relative to exposure to unpredictable environments. We suggest that such heterochronic shifts in ossification can not only buffer development under fluctuating environments while maintaining epigenetic sensitivity critical for normal skeletal formation, but also enable epigenetically induced gene expression to generate specialized morphological adaptations. We review studies of environmental sensitivity of BMP pathways and their regulation of formation, remodeling, and repair of cartilage and bone to examine the hypothesis that BMP-mediated skeletal adaptations are facilitated by evolved reactivity of BMPs to external signals. Surprisingly, no empirical study to date has identified the molecular mechanism behind developmental plasticity in skeletal traits. We outline a conceptual framework for future studies that focus on mediation of phenotypic plasticity in skeletal development by the patterns of BMP expression.

Evolutionary persistence of phenotypic integration: influence of developmental and functional relationships on complex trait evolution.

Examination of historical persistence of integration patterns provides an important insight into understanding the origin and evolution of complex traits. Specifically, the distinct effects of developmental and functional integration on the evolution of complex traits are often overlooked. Because patterns of functional integration are commonly shaped by selection exerted by the external environment, whereas patterns of developmental integration can be determined by relatively environment-independent selection for developmental homeostasis, examination of historical persistence of morphological integration patterns among species should reveal the relative importance of current selection in the evolution of complex traits. We compared historical persistence of integration patterns produced by current developmental versus ecological requirements by examining the evolution of complex mandibular structures in nine species of soricid shrews. We found that, irrespective of phylogenetic relatedness of species, patterns of developmental and functional integration were highly concordant, suggesting that strong selection for developmental homeostasis favors concordant channeling of both internal and external variation. Overall, our results suggest that divergence in mandible shape among species closely follows variation in functional demands and ecological requirements regardless of phylogenetic relatedness among species.

Evolution of morphological integration: developmental accommodation of stress-induced variation.

Extreme environmental change during growth often results in an increase in developmental abnormalities in the morphology of an organism. The evolutionary significance of such stress-induced variation depends on the recurrence of a stressor and on the degree to which developmental errors can be accommodated by an organism's ontogeny without significant loss of function. We subjected populations of four species of soricid shrews to an extreme environment during growth and measured changes in the patterns of integration and accommodation of stress-induced developmental errors in a complex of mandibular traits. Adults that grew under an extreme environment had lower integration of morphological variation among mandibular traits and highly elevated fluctuating asymmetry in these traits, compared to individuals that grew under the control conditions. However, traits differed strongly in the magnitude of response to a stressor--traits within attachments of the same muscle (functionally integrated traits) had lower response and changed their integration less than other traits. Cohesiveness in functionally integrated complexes of traits under stress was maintained by close covariation of their developmental variation. Such developmental accommodation of stress-induced variation might enable the individual's functioning and persistence under extreme environmental conditions and thus provides a link between individual adaptation to stress and the evolution of stress resistance.