Conspicuous colors (e.g., red, yellow, blue) have evolved numerous times across animals. But the function of this coloration can differ radically among species. Many species use this coloration as a sexual signal to conspecifics, whereas others use it as a warning signal to predators. Why do different species evolve conspicuous coloration in association with one function as opposed to the other? We address this question in terrestrial vertebrates (tetrapods) using phylogenetic approaches, and test whether day-night activities of species help determine these patterns. Using phylogenetic logistic regression, we found that conspicuous, sexually dimorphic coloration is significantly associated with diurnal lineages (e.g., many birds and lizards). By contrast, the evolution of warning signals was significantly associated with large-scale clades that were ancestrally nocturnal (e.g., snakes, amphibians), regardless of the current diel activity of species. Overall, we show that the evolution of conspicuous coloration as warning signals or sexual signals is influenced by the ecology of species, both recently and in the ancient past.

Rensch's rule suggests that sexual size dimorphism (SSD) increases with species size when males are the larger sex, whereas it decreases when females are the larger sex. However, the process responsible for this pattern remains obscure. SSD can result from sexual selection, such as intrasexual competition for access to mates, or from natural selection, due to resource partitioning or fecundity selection. We studied SSD in Caribbean Eleutherodactylus frogs using phylogenetic comparative methods to investigate the influence of microhabitat, fecundity, and parental care. Our results show that in Caribbean Eleutherodactylus females tend to be larger and, contrary to Rensch's rule, dimorphism increases with species size. SSD was not related to microhabitat use. However, SSD was positively correlated with fecundity, mediated by a greater increase in female size. SSD was also influenced by parental care, suggesting that male care promotes larger male size and reduces the female bias in SSD. As suggested for other anurans, female-biased SSD in Caribbean Eleutherodactylus results from fecundity selection, although the magnitude is countered by increased male size in species with paternal care. Our results highlight the importance of considering various selective forces that may act in concert to influence the evolution of SSD.

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

The occurrence of similar phenotypes in multiple independent populations derived from common ancestral conditions (viz. parallel evolution) is a testimony of evolution by natural selection. Parallel evolution implies that populations share a common phenotypic response to a common selection pressure associated with habitat similarity. Examples of parallel evolution at genetic and phenotypic levels are fairly common, but the driving selective agents often remain elusive. Similarly, the role of phenotypic plasticity in facilitating early stages of parallel evolution is unclear. We investigated whether the relaxation of predation pressure associated with the colonization of freshwater ponds by nine-spined sticklebacks (Pungitius pungitius) likely explains the divergence in complex behaviors between marine and pond populations, and whether this divergence is parallel. Using laboratory-raised individuals exposed to different levels of perceived predation risk, we calculated vectors of phenotypic divergence for four behavioral traits between habitats and predation risk treatments. We found a significant correlation between the directions of evolutionary divergence and phenotypic plasticity, suggesting that divergence in behavior between habitats is aligned with the response to relaxation of predation pressure. Finally, we show alignment across multiple pairs of populations, and that relaxation of predation pressure has likely driven parallel evolution of behavior in this species.

Multiple ant lineages have evolved a bizarre system called social hybridogenesis, involving multiple co-occurring genetic lineages, in which mating between lineages produces workers but mating within a lineage produces daughter queens. A new study reveals that this system evolved multiple times within harvester ants, each time from interspecific hybridization. A third finding, that the system likely evolves in small or isolated populations, could be explained either by exploitation of heterospecific males for their sperm, or simply by failure to avoid interspecific mating.

The availability of genome sequences from large numbers of species offers the prospect of studying genotype-phenotype correlations across various phylogenetic scales using only available data. A new study illustrates the power of this approach, showing an association across primates between morphological sexual dimorphism and the prevalence of a class of DNA elements that stimulate gene expression in response to male androgens.

Genes of the major histocompatibility complex (MHC) code for immune proteins that are crucial for pathogen recognition in vertebrates. MHC research in nonmodel taxa has long been hampered by its genomic complexity that makes the locus-specific genotyping challenging. The recent progress in sequencing and genotyping methodologies allows an extensive phylogenetic coverage in studies of MHC evolution. Here, we analyzed the peptide-binding region of MHC class I (MHC-I)in 30 species of salamanders from six families representative of Urodela phylogeny. This extensive dataset revealed an extreme diversity of MHC-I in salamanders, both in terms of sequence diversity (about 3000 variants) and architecture (2-22 gene copies per species). The signal of positive selection was moderate and consistent between both peptide-binding domains, but varied greatly between genera. Positions of positively selected sites mostly coincided with human peptide-binding sites, suggesting similar structural properties of MHC-I molecules across distant vertebrate lineages. Finally, we provided evidence for the common intraexonic recombination at MHC-I and for the role of life history traits in the processes of MHC-I expansion/contraction. Our study revealed novel evolutionary trajectories of amphibian MHC and it contributes to the understanding of the mechanisms that generated extraordinary MHC diversity throughout vertebrate evolution.

The archipelagos of Wallacea extend between the Sunda and Sahul Shelves, serving as a semipermeable two-way filter influencing faunal exchange between Asia and Australo-Papua. Forest skinks (Genus Sphenomorphus) are widespread throughout southern Wallacea and exhibit complex clinal, ontogenetic, sexual, and seasonal morphological variation, rendering species delimitation difficult. We screened a mitochondrial marker for 245 Sphenomorphus specimens from this area to inform the selection of 104 samples from which we used targeted sequence capture to generate a dataset of 1154 nuclear genes (∼1.8 Mb) plus complete mitochondrial genomes. Phylogenomic analyses recovered many deeply divergent lineages, three pairs of which are now sympatric, that began to diversify in the late Miocene shortly after the oldest islands are thought to have become emergent. We infer a complex and nonstepping-stone pattern of island colonization, with the group having originated in the Sunda Arc islands before using Sumba as a springboard for colonization of the Banda Arcs. Estimates of population structure and gene flow across the region suggest total isolation except between two Pleistocene Aggregate Island Complexes that become episodically land-bridged during glacial maxima. These historical processes have resulted in at least 11 Sphenomorphus species in the region, nine of which require formal description. This fine-scale geographic partitioning of undescribed species highlights the importance of utilizing comprehensive genomic studies for defining biodiversity hotspots to be considered for conservation protection.

In neutral models of quantitative trait evolution, both genetic and phenotypic divergence scale as random walks, producing a correlation between the two measures. However, complexity in the genotype-phenotype map may alter the correlation between genotypic and phenotypic divergence, even when both are evolving neutrally or nearly so. Understanding this correlation between phenotypic and genetic variation is critical for accurately interpreting the fossil record. This study compares the geographic structure and scaling of morphological variation of the shape of the first lower molar of 77 individuals of the northern grasshopper mouse Onychomys leucogaster to genome-wide SNP variation in the same sample. We found strong genetic structure but weak or absent morphological structure indicating that the scaling of each type of variation is decoupled from one another. Low PST values relative to FST values are consistent with a lack of morphological divergence in contrast to genetic divergence between groups. This lack of phenotypic structure and the presence of notable within-sample phenotypic variance are consistent with uniform selection or constraints on molar shape across a wide geographic and environmental range. Over time, this kind of decoupling may result in patterns of phenotypic stasis masking underlying genetic patterns.