Adaptive tail-length evolution in deer mice is associated with differential Hoxd13 expression in early development.

Variation in the size and number of axial segments underlies much of the diversity in animal body plans. Here we investigate the evolutionary, genetic and developmental mechanisms driving tail-length differences between forest and prairie ecotypes of deer mice (Peromyscus maniculatus). We first show that long-tailed forest mice perform better in an arboreal locomotion assay, consistent with tails being important for balance during climbing. We then identify six genomic regions that contribute to differences in tail length, three of which associate with caudal vertebra length and the other three with vertebra number. For all six loci, the forest allele increases tail length, indicative of the cumulative effect of natural selection. Two of the genomic regions associated with variation in vertebra number contain Hox gene clusters. Of those, we find an allele-specific decrease in Hoxd13 expression in the embryonic tail bud of long-tailed forest mice, consistent with its role in axial elongation. Additionally, we find that forest embryos have more presomitic mesoderm than prairie embryos and that this correlates with an increase in the number of neuromesodermal progenitors, which are modulated by Hox13 paralogues. Together, these results suggest a role for Hoxd13 in the development of natural variation in adaptive morphology on a microevolutionary timescale.

Evolution of gene expression across brain regions in behaviourally divergent deer mice.

The evolution of innate behaviours is ultimately due to genetic variation likely acting in the nervous system. Gene regulation may be particularly important because it can evolve in a modular brain-region specific fashion through the concerted action of cis- and trans-regulatory changes. Here, to investigate transcriptional variation and its regulatory basis across the brain, we perform RNA sequencing (RNA-Seq) on ten brain subregions in two sister species of deer mice (Peromyscus maniculatus and P. polionotus)-which differ in a range of innate behaviours, including their social system-and their F1 hybrids. We find that most of the variation in gene expression distinguishes subregions, followed by species. Interspecific differential expression (DE) is pervasive (52-59% of expressed genes), whereas the number of DE genes between sexes is modest overall (~3%). Interestingly, the identity of DE genes varies considerably across brain regions. Much of this modularity is due to cis-regulatory divergence, and while 43% of genes were consistently assigned to the same gene regulatory class across subregions (e.g. conserved, cis- or trans-regulatory divergence), a similar number were assigned to two or more different gene regulatory classes. Together, these results highlight the modularity of gene expression differences and divergence in the brain, which may be key to explain how the evolution of brain gene expression can contribute to the astonishing diversity of animal behaviours.

Accounting for diverse transposable element landscapes is key to developing and evaluating accurate de novo annotation strategies.

Transposable elements (TEs) are important drivers of genome evolution. Nonetheless, TE annotation remains a complex and challenging task. As more genomes from phylogenetically diverse species are published, a comprehensive pipeline for accurate annotation of diverse TEs is increasingly important. Recently, (Ou et al. Genome Biol. 20:275, 2019) developed a new comprehensive pipeline, Extensive De novo Transposable element Annotator (EDTA), and benchmarked its performance on the genomes of three species: maize, wheat, and fruit fly. Because TE landscapes can vary tremendously across species, we tested EDTA's performance on four additional genomes with different TE landscapes: mouse, zebrafish, zebra finch, and chicken. Our analysis reveals that EDTA faces challenges with repeat classification in these genomes and underperforms overall relative to its benchmark dataset. Notably, EDTA consistently misclassifies nonLTR retrotransposons as DNA transposons, resulting in erroneous TE annotations for species with considerable repertoires of nonLTR retrotransposons. Overall, we set expectations for EDTA's performance on genomes spanning additional diversity, urge caution when using EDTA on genomes with divergent TE repertoires from the species on which it was initially benchmarked, and hope to motivate the development of methods that are robust to both the diversity of TEs and TE landscapes observed across species.

The neural basis of defensive behaviour evolution in Peromyscus mice.

Evading imminent predator threat is critical for survival. Effective defensive strategies can vary, even between closely related species. However, the neural basis of such species-specific behaviours is still poorly understood. Here we find that two sister species of deer mice (genus Peromyscus) show different responses to the same looming stimulus: P. maniculatus, which occupy densely vegetated habitats, predominantly dart to escape, while the open field specialist, P. polionotus, pause their movement. This difference arises from species-specific escape thresholds, is largely context-independent, and can be triggered by both visual and auditory threat stimuli. Using immunohistochemistry and electrophysiological recordings, we find that although visual threat activates the superior colliculus in both species, the role of the dorsal periaqueductal gray (dPAG) in driving behaviour differs. While dPAG activity scales with running speed and involves both excitatory and inhibitory neurons in P. maniculatus, the dPAG is largely silent in P. polionotus, even when darting is triggered. Moreover, optogenetic activation of excitatory dPAG neurons reliably elicits darting behaviour in P. maniculatus but not P. polionotus. Together, we trace the evolution of species-specific escape thresholds to a central circuit node, downstream of peripheral sensory neurons, localizing an ecologically relevant behavioural difference to a specific region of the complex mammalian brain.

A single genetic locus lengthens deer mouse burrows via motor pattern evolution.

The question of how evolution builds complex behaviors has long fascinated biologists. To address this question from a genetic perspective, we capitalize on variation in innate burrowing behavior between two sister species of Peromyscus mice: P. maniculatus that construct short, simple burrows and P. polionotus that uniquely construct long, elaborate burrows. We identify three regions of the genome associated with differences in burrow length and then narrow in on one large-effect 12-Mb locus on chromosome 4. By introgressing the P. polionotus allele into a P. maniculatus background, we demonstrate this locus, on its own, increases burrow length by 20%. Next, by recording mice digging in a transparent tube, we find this locus has specific effects on burrowing behavior. This locus does not affect time spent digging or latency to dig, but rather affects usage of only two of the primary digging behaviors that differ between the focal species: forelimb digging, which loosens substrate, and hindlimb kicking, which powerfully ejects substrate. This locus has an especially large effect on hindkicking, explaining 56% and 22% of interspecific differences in latency and proportion of hindkicks, respectively. Together, these data provide genetic support for the hierarchical organization of complex behaviors, offering evolution the opportunity to tinker with specific behavioral components.

Two pup vocalization types are genetically and functionally separable in deer mice.

Vocalization is a widespread social behavior in vertebrates that can affect fitness in the wild. Although many vocal behaviors are highly conserved, heritable features of specific vocalization types can vary both within and between species, raising the questions of why and how some vocal behaviors evolve. Here, using new computational tools to automatically detect and cluster vocalizations into distinct acoustic categories, we compare pup isolation calls across neonatal development in eight taxa of deer mice (genus Peromyscus) and compare them with laboratory mice (C57BL6/J strain) and free-living, wild house mice (Mus musculus domesticus). Whereas both Peromyscus and Mus pups produce ultrasonic vocalizations (USVs), Peromyscus pups also produce a second call type with acoustic features, temporal rhythms, and developmental trajectories that are distinct from those of USVs. In deer mice, these lower frequency "cries" are predominantly emitted in postnatal days one through nine, whereas USVs are primarily made after day 9. Using playback assays, we show that cries result in a more rapid approach by Peromyscus mothers than USVs, suggesting a role for cries in eliciting parental care early in neonatal development. Using a genetic cross between two sister species of deer mice exhibiting large, innate differences in the acoustic structure of cries and USVs, we find that variation in vocalization rate, duration, and pitch displays different degrees of genetic dominance and that cry and USV features can be uncoupled in second-generation hybrids. Taken together, this work shows that vocal behavior can evolve quickly between closely related rodent species in which vocalization types, likely serving distinct functions in communication, are controlled by distinct genetic loci.

Transposable Element Interactions Shape the Ecology of the Deer Mouse Genome.

The genomic landscape of transposable elements (TEs) varies dramatically across species, with some TEs demonstrating greater success in colonizing particular lineages than others. In mammals, long interspersed nuclear element (LINE) retrotransposons are typically more common than any other TE. Here, we report an unusual genomic landscape of TEs in the deer mouse, Peromyscus maniculatus. In contrast to other previously examined mammals, long terminal repeat elements occupy more of the deer mouse genome than LINEs (11% and 10%, respectively). This pattern reflects a combination of relatively low LINE activity and a massive invasion of lineage-specific endogenous retroviruses (ERVs). Deer mouse ERVs exhibit diverse origins spanning the retroviral phylogeny suggesting they have been host to a wide range of exogenous retroviruses. Notably, we trace the origin of one ERV lineage, which arose ∼5-18 million years ago, to a close relative of feline leukemia virus, revealing inter-ordinal horizontal transmission. Several lineage-specific ERV subfamilies have very high copy numbers, with the top five most abundant accounting for ∼2% of the genome. We also observe a massive amplification of Kruppel-associated box domain-containing zinc finger genes, which likely control ERV activity and whose expansion may have been facilitated by ectopic recombination between ERVs. Finally, we find evidence that ERVs directly impacted the evolutionary trajectory of LINEs by outcompeting them for genomic sites and frequently disrupting autonomous LINE copies. Together, our results illuminate the genomic ecology that shaped the unique deer mouse TE landscape, shedding light on the evolutionary processes that give rise to variation in mammalian genome structure.

Post-mating parental behavior trajectories differ across four species of deer mice.

Among species, parental behaviors vary in their magnitude, onset relative to reproduction, and sexual dimorphism. In deer mice (genus Peromyscus), while most species are promiscuous with low paternal care, monogamy and biparental care have evolved at least twice under different ecological conditions. Here, in a common laboratory setting, we monitored parental behaviors of males and females of two promiscuous (eastern deer mouse P. maniculatus and white-footed mouse P. leucopus) and two monogamous (oldfield mouse P. polionotus and California mouse P. californicus) species from before mating to after giving birth. In the promiscuous species, females showed parental behaviors largely after parturition, while males showed little parental care. In contrast, both sexes of monogamous species performed parental behaviors. However, while oldfield mice began to display parental behaviors before mating, California mice showed robust parental care behaviors only postpartum. These different parental-care trajectories in the two monogamous species align with their socioecology. Oldfield mice have overlapping home ranges with relatives, so infants they encounter, even if not their own, are likely to be closely related. By contrast, California mice disperse longer distances into exclusive territories with possibly unrelated neighbors, decreasing the inclusive fitness benefits of caring for unfamiliar pups before parenthood. Together, we find that patterns of parental behaviors in Peromyscus are consistent with predictions from inclusive fitness theory.

Chromosomal inversion polymorphisms shape the genomic landscape of deer mice.

Chromosomal inversions are an important form of structural variation that can affect recombination, chromosome structure and fitness. However, because inversions can be challenging to detect, the prevalence and hence the significance of inversions segregating within species remains largely unknown, especially in natural populations of mammals. Here, by combining population-genomic and long-read sequencing analyses in a single, widespread species of deer mouse (Peromyscus maniculatus), we identified 21 polymorphic inversions that are large (1.5-43.8 Mb) and cause near-complete suppression of recombination when heterozygous (0-0.03 cM Mb-1). We found that inversion breakpoints frequently occur in centromeric and telomeric regions and are often flanked by long inverted repeats (0.5-50 kb), suggesting that they probably arose via ectopic recombination. By genotyping inversions in populations across the species' range, we found that the inversions are often widespread and do not harbour deleterious mutational loads, and many are likely to be maintained as polymorphisms by divergent selection. Comparisons of forest and prairie ecotypes of deer mice revealed 13 inversions that contribute to differentiation between populations, of which five exhibit significant associations with traits implicated in local adaptation. Taken together, these results show that inversion polymorphisms have a significant impact on recombination, genome structure and genetic diversity in deer mice and likely facilitate local adaptation across the widespread range of this species.

Interspecific variation in cooperative burrowing behavior by Peromyscus mice.

Animals often adjust their behavior according to social context, but the capacity for such behavioral flexibility can vary among species. Here, we test for interspecific variation in behavioral flexibility by comparing burrowing behavior across three species of deer mice (genus Peromyscus) with divergent social systems, ranging from promiscuous (Peromyscus leucopus and Peromyscus maniculatus) to monogamous (Peromyscus polionotus). First, we compared the burrows built by individual mice to those built by pairs of mice in all three species. Although burrow length did not differ in P. leucopus or P. maniculatus, we found that P. polionotus pairs cooperatively constructed burrows that were nearly twice as long as those built by individuals and that opposite-sex pairs dug longer burrows than same-sex pairs. Second, to directly observe cooperative digging behavior in P. polionotus, we designed a burrowing assay in which we could video-record active digging in narrow, transparent enclosures. Using this novel assay, we found, unexpectedly, that neither males nor females spent more time digging with an opposite-sex partner. Rather, we demonstrate that opposite-sex pairs are more socially cohesive and thus more efficient digging partners than same-sex pairs. Together, our study demonstrates how social context can modulate innate behavior and offers insight into how differences in behavioral flexibility may evolve among closely related species.

An enhancer of Agouti contributes to parallel evolution of cryptically colored beach mice.

Identifying the genetic basis of repeatedly evolved traits provides a way to reconstruct their evolutionary history and ultimately investigate the predictability of evolution. Here, we focus on the oldfield mouse (Peromyscus polionotus), which occurs in the southeastern United States, where it exhibits considerable color variation. Dorsal coats range from dark brown in mainland mice to near white in mice inhabiting sandy beaches; this light pelage has evolved independently on Florida's Gulf and Atlantic coasts as camouflage from predators. To facilitate genomic analyses, we first generated a chromosome-level genome assembly of Peromyscus polionotus subgriseus. Next, in a uniquely variable mainland population (Peromyscus polionotus albifrons), we scored 23 pigment traits and performed targeted resequencing in 168 mice. We find that pigment variation is strongly associated with an ∼2-kb region ∼5 kb upstream of the Agouti signaling protein coding region. Using a reporter-gene assay, we demonstrate that this regulatory region contains an enhancer that drives expression in the dermis of mouse embryos during the establishment of pigment prepatterns. Moreover, extended tracts of homozygosity in this Agouti region indicate that the light allele experienced recent and strong positive selection. Notably, this same light allele appears fixed in both Gulf and Atlantic coast beach mice, despite these populations being separated by >1,000 km. Together, our results suggest that this identified Agouti enhancer allele has been maintained in mainland populations as standing genetic variation and from there, has spread to and been selected in two independent beach mouse lineages, thereby facilitating their rapid and parallel evolution.

Behavioral genetics and genomics: Mendel's peas, mice, and bees.

The question of the heritability of behavior has been of long fascination to scientists and the broader public. It is now widely accepted that most behavioral variation has a genetic component, although the degree of genetic influence differs widely across behaviors. Starting with Mendel's remarkable discovery of "inheritance factors," it has become increasingly clear that specific genetic variants that influence behavior can be identified. This goal is not without its challenges: Unlike pea morphology, most natural behavioral variation has a complex genetic architecture. However, we can now apply powerful genome-wide approaches to connect variation in DNA to variation in behavior as well as analyses of behaviorally related variation in brain gene expression, which together have provided insights into both the genetic mechanisms underlying behavior and the dynamic relationship between genes and behavior, respectively, in a wide range of species and for a diversity of behaviors. Here, we focus on two systems to illustrate both of these approaches: the genetic basis of burrowing in deer mice and transcriptomic analyses of division of labor in honey bees. Finally, we discuss the troubled relationship between the field of behavioral genetics and eugenics, which reminds us that we must be cautious about how we discuss and contextualize the connections between genes and behavior, especially in humans.

A chromosomal inversion contributes to divergence in multiple traits between deer mouse ecotypes.

How locally adapted ecotypes are established and maintained within a species is a long-standing question in evolutionary biology. Using forest and prairie ecotypes of deer mice (Peromyscus maniculatus), we characterized the genetic basis of variation in two defining traits-tail length and coat color-and discovered a 41-megabase chromosomal inversion linked to both. The inversion frequency is 90% in the dark, long-tailed forest ecotype; decreases across a habitat transition; and is absent from the light, short-tailed prairie ecotype. We implicate divergent selection in maintaining the inversion at frequencies observed in the wild, despite high levels of gene flow, and explore fitness benefits that arise from suppressed recombination within the inversion. We uncover a key role for a large, previously uncharacterized inversion in the evolution and maintenance of classic mammalian ecotypes.

cis-Regulatory changes in locomotor genes are associated with the evolution of burrowing behavior.

How evolution modifies complex, innate behaviors is largely unknown. Divergence in many morphological traits, and some behaviors, is linked to cis-regulatory changes in gene expression. Given this, we compare brain gene expression of two interfertile sister species of Peromyscus mice that show large and heritable differences in burrowing behavior. Species-level differential expression and allele-specific expression in F1 hybrids indicate a preponderance of cis-regulatory divergence, including many genes whose cis-regulation is affected by burrowing behavior. Genes related to locomotor coordination show the strongest signals of lineage-specific selection on burrowing-induced cis-regulatory changes. Furthermore, genetic markers closest to these candidate genes associate with variation in burrow shape in a genetic cross, suggesting an enrichment for loci affecting burrowing behavior near these candidate locomotor genes. Our results provide insight into how cis-regulated gene expression can depend on behavioral context and how this dynamic regulatory divergence between species may contribute to behavioral evolution.

Tail Length Evolution in Deer Mice: Linking Morphology, Behavior, and Function.

Determining how variation in morphology affects animal performance (and ultimately fitness) is key to understanding the complete process of evolutionary adaptation. Long tails have evolved many times in arboreal and semi-arboreal rodents; in deer mice, long tails have evolved repeatedly in populations occupying forested habitat even within a single species (Peromyscus maniculatus). Here, we use a combination of functional modeling, laboratory studies, and museum records to test hypotheses about the function of tail-length variation in deer mice. First, we use computational models, informed by museum records documenting natural variation in tail length, to test whether differences in tail morphology between forest and prairie subspecies can influence performance in behavioral contexts relevant for tail use. We find that the deer- mouse tail plays little role in statically adjusting center of mass or in correcting body pitch and yaw, but rather it can affect body roll during arboreal locomotion. In this context, we find that even intraspecific tail-length variation could result in substantial differences in how much body rotation results from equivalent tail motions (i.e., tail effectiveness), but the relationship between commonly-used metrics of tail-length variation and effectiveness is non-linear. We further test whether caudal vertebra length, number, and shape are associated with differences in how much the tail can bend to curve around narrow substrates (i.e., tail curvature) and find that, as predicted, the shape of the caudal vertebrae is associated with intervertebral bending angle across taxa. However, although forest and prairie mice typically differ in both the length and number of caudal vertebrae, we do not find evidence that this pattern is the result of a functional trade-off related to tail curvature. Together, these results highlight how even simple models can both generate and exclude hypotheses about the functional consequences of trait variation for organismal-level performance.

Expanding evolutionary neuroscience: insights from comparing variation in behavior.

Neuroscientists have long studied species with convenient biological features to discover how behavior emerges from conserved molecular, neural, and circuit level processes. With the advent of new tools, from viral vectors and gene editing to automated behavioral analyses, there has been a recent wave of interest in developing new, "nontraditional" model species. Here, we advocate for a complementary approach to model species development, that is, model clade development, as a way to integrate an evolutionary comparative approach with neurobiological and behavioral experiments. Capitalizing on natural behavioral variation in and investing in experimental tools for model clades will be a valuable strategy for the next generation of neuroscience discovery.

Fishing for the genetic basis of migratory behavior.

For many species, migrating at just the right time is essential for both survival and reproduction. A new study in salmon localizes a small genomic region associated with migration timing, which in turn affects other physiological traits, suggesting that a seemingly complex suite of migration traits is linked by one "simple" phenotype.

The Tug1 lncRNA locus is essential for male fertility.

BACKGROUND: Several long noncoding RNAs (lncRNAs) have been shown to function as components of molecular machines that play fundamental roles in biology. While the number of annotated lncRNAs in mammalian genomes has greatly expanded, studying lncRNA function has been a challenge due to their diverse biological roles and because lncRNA loci can contain multiple molecular modes that may exert function. RESULTS: We previously generated and characterized a cohort of 20 lncRNA loci knockout mice. Here, we extend this initial study and provide a more detailed analysis of the highly conserved lncRNA locus, taurine-upregulated gene 1 (Tug1). We report that Tug1-knockout male mice are sterile with underlying defects including a low number of sperm and abnormal sperm morphology. Because lncRNA loci can contain multiple modes of action, we wanted to determine which, if any, potential elements contained in the Tug1 genomic region have any activity. Using engineered mouse models and cell-based assays, we provide evidence that the Tug1 locus harbors two distinct noncoding regulatory activities, as a cis-DNA repressor that regulates neighboring genes and as a lncRNA that can regulate genes by a trans-based function. We also show that Tug1 contains an evolutionary conserved open reading frame that when overexpressed produces a stable protein which impacts mitochondrial membrane potential, suggesting a potential third coding function. CONCLUSIONS: Our results reveal an essential role for the Tug1 locus in male fertility and uncover evidence for distinct molecular modes in the Tug1 locus, thus highlighting the complexity present at lncRNA loci.

Diet-based assortative mating through sexual imprinting.

Speciation is facilitated by "magic traits," where divergent natural selection on such traits also results in assortative mating. In animal populations, diet has the potential to act as a magic trait if populations diverge in consumed food that incidentally affects mating and therefore sexual isolation. While diet-based assortative mating has been observed in the laboratory and in natural populations, the mechanisms causing positive diet-based assortment remain largely unknown. Here, we experimentally created divergent diets in a sexually imprinting species of mouse, Peromyscus gossypinus (the cotton mouse), to test the hypothesis that sexual imprinting on diet could be a mechanism that generates rapid and significant sexual isolation. We provided breeding pairs with novel garlic- or orange-flavored water and assessed whether their offspring, exposed to these flavors in utero and in the nest before weaning, later preferred mates that consumed the same flavored water as their parents. While males showed no preference, females preferred males of their parental diet, which is predicted to yield moderate sexual isolation. Thus, our experiment demonstrates the potential for sexual imprinting on dietary cues learned in utero and/or postnatally to facilitate reproductive isolation and potentially speciation. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.n1qq6v3.