Snowshoe hares () maintain seasonal camouflage by molting to a white winter coat, but some hares remain brown during the winter in regions with low snow cover. We show that cis-regulatory variation controlling seasonal expression of the gene underlies this adaptive winter camouflage polymorphism. Genetic variation at clustered by winter coat color across multiple hare and jackrabbit species, revealing a history of recurrent interspecific gene flow. Brown winter coats in snowshoe hares likely originated from an introgressed black-tailed jackrabbit allele that has swept to high frequency in mild winter environments. These discoveries show that introgression of genetic variants that underlie key ecological traits can seed past and ongoing adaptation to rapidly changing environments.

Summary The yellow and red feather pigmentation of many bird species [1] plays pivotal roles in social signaling and mate choice [2, 3]. To produce red pigments, birds ingest yellow carotenoids and endogenously convert them into red ketocarotenoids via an oxidation reaction catalyzed by a previously unknown ketolase [4–6]. We investigated the genetic basis for red coloration in birds using whole-genome sequencing of red siskins (Spinus cucullata), common canaries (Serinus canaria), and 'red factor' canaries, which are the hybrid product of crossing red siskins with common canaries [7]. We identified two genomic regions introgressed from red siskins into red factor canaries that are required for red coloration. One of these regions contains a gene encoding a cytochrome P450 enzyme, CYP2J19. Transcriptome analysis demonstrates that CYP2J19 is significantly upregulated in the skin and liver of red factor canaries, strongly implicating CYP2J19 as the ketolase that mediates red coloration in birds. Interestingly, a second introgressed region required for red feathers resides within the epidermal differentiation complex, a cluster of genes involved in development of the integument. Lastly, we present evidence that CYP2J19 is involved in ketocarotenoid formation in the retina. The discovery of the carotenoid ketolase has important implications for understanding sensory function and signaling mediated by carotenoid pigmentation.

Mitochondrial DNA introgression from Lepus timidus into Lepus granatensis and Lepus europaeus was recently reported in Iberia, although L. timidus presumably retreated from this region at the end of the last ice age. Here we assess the extent of this ancient mtDNA introgression by RFLP analysis of 695 specimens representing the three hare species present in Iberia. The introgressed L. timidus lineage was found in 23 of the 37 populations sampled. It is almost fixed in L. europaeus across its Iberian range in the Pyrenean foothills, and in L. granatensis, which occupies the rest of the peninsula, it is predominant in the north and gradually disappears further south. We also found it in Lepus castroviejoi, a species endemic to Cantabria. Multiple hybridizations and, potentially, a selective advantage for the L. timidus lineage can explain the remarkable taxonomic and geographical range of this mitochondrial introgression.

In animals, the population genomic literature is dominated by two taxa, namely mammals and drosophilids, in which fully sequenced, well-annotated genomes have been available for years. Data from other metazoan phyla are scarce, probably because the vast majority of living species still lack a closely related reference genome. Here we achieve de novo, reference-free population genomic analysis from wild samples in five non-model animal species, based on next-generation sequencing transcriptome data. We introduce a pipe-line for cDNA assembly, read mapping, SNP/genotype calling, and data cleaning, with specific focus on the issue of hidden paralogy detection. In two species for which a reference genome is available, similar results were obtained whether the reference was used or not, demonstrating the robustness of our de novo inferences. The population genomic profile of a hare, a turtle, an oyster, a tunicate, and a termite were found to be intermediate between those of human and Drosophila, indicating that the discordant genomic diversity patterns that have been reported between these two species do not reflect a generalized vertebrate versus invertebrate gap. The genomic average diversity was generally higher in invertebrates than in vertebrates (with the notable exception of termite), in agreement with the notion that population size tends to be larger in the former than in the latter. The non-synonymous to synonymous ratio, however, did not differ significantly between vertebrates and invertebrates, even though it was negatively correlated with genetic diversity within each of the two groups. This study opens promising perspective regarding genome-wide population analyses of non-model organisms and the influence of population size on non-synonymous versus synonymous diversity.