Pigmentation Genes Show Evidence of Repeated Divergence and Multiple Bouts of Introgression in Setophaga Warblers

Museum genomics provide insight into the extinction of a specialist North American warbler species

Museum genomics provide an opportunity to investigate population demographics of extinct species, especially valuable when research prior to extinction was minimal. The Bachman's warbler (Vermivora bachmanii) is hypothesized to have gone extinct due to loss of its specialized habitat. However, little is known about other potential contributing factors such as natural rarity or changes to connectivity following habitat fragmentation. We examined mitochondrial DNA (mtDNA) and genome-wide SNPs using specimens collected from breeding and migration sites across the range of the Bachman's warbler. We found no signals of strong population structuring across the breeding range of Bachman's warblers in both mtDNA and genome-wide SNPs. Thus, long-term population isolation did not appear to be a significant contributor to the extinction of the Bachman's warbler. Instead, our findings support the theory that Bachman's warblers underwent a rapid decline likely driven by habitat destruction, which may have been exacerbated by the natural rarity, habitat specificity and low genetic diversity of the species.

Genetic Basis and Evolutionary Forces of Sexually Dimorphic Color Variation in a Toad-Headed Agamid Lizard

In the animal kingdom, sexually dimorphic color variation is a widespread phenomenon that significantly influences survival and reproductive success. However, the genetic underpinnings of this variation remain inadequately understood. Our investigation into sexually dimorphic color variation in the desert-dwelling Guinan population of the toad-headed agamid lizard (Phrynocephalus putjatai) utilized a multidisciplinary approach, encompassing phenotypic, ultrastructural, biochemical, genomic analyses, and behavioral experiments. Our findings unveil the association between distinct skin colorations and varying levels of carotenoid and pteridine pigments. The red coloration in males is determined by a genomic region on chromosome 14, housing four pigmentation genes: BCO2 and three 6-pyruvoyltetrahydropterin synthases. A Guinan population-specific nonsynonymous single nucleotide polymorphism in BCO2 is predicted to alter the electrostatic potential within the binding domain of the BCO2-β-carotene complex, influencing their interaction. Additionally, the gene MAP7 on chromosome 2 emerges as a potential contributor to the blue coloration in subadults and adult females. Sex-specific expression patterns point to steroid hormone-associated genes (SULT2B1 and SRD5A2) as potential upstream regulators influencing sexually dimorphic coloration. Visual modeling and field experiments support the potential selective advantages of vibrant coloration in desert environments. This implies that natural selection, potentially coupled with assortative mating, might have played a role in fixing color alleles, contributing to prevalence in the local desert habitat. This study provides novel insights into the genetic basis of carotenoid and pteridine-based color variation, shedding light on the evolution of sexually dimorphic coloration in animals. Moreover, it advances our understanding of the driving forces behind such intricate coloration patterns.

Whole-Genome Resequencing Reveals Genetic Diversity and Wool Trait-Related Genes in Liangshan Semi-Fine-Wool Sheep

Understanding the genetic makeup of local sheep breeds is essential for their scientific conservation and sustainable utilization. The Liangshan semi-fine-wool sheep (LSS), a Chinese semi-fine-wool breed renowned for its soft wool, was analyzed using whole-genome sequencing data including 35 LSS, 84 sheep from other domestic breeds, and 20 Asiatic mouflons. We investigated the genetic composition of LSS by conducting analyses of the population structure, runs of homozygosity, genomic inbreeding coefficients, and selection signature. Our findings indicated that LSS shares greater genetic similarity with Border Leicester and Romney sheep than with Tibetan (TIB), Yunnan (YNS), and Chinese Merino sheep. Genomic analysis indicated low to moderate inbreeding coefficients, ranging from 0.014 to 0.154. In identifying selection signals across the LSS genome, we pinpointed 195 candidate regions housing 74 annotated genes (e.g., IRF2BP2, BVES, and ALOX5). We also found the overlaps between the candidate regions and several known quantitative trait loci related to wool traits, such as the wool staple length and wool fiber diameter. A selective sweep region, marked by the highest value of cross-population extended haplotype homozygosity, encompassed IRF2BP2-an influential candidate gene affecting fleece fiber traits. Furthermore, notable differences in genotype frequency at a mutation site (c.1051 + 46T > C, Chr25: 6,784,190 bp) within IRF2BP2 were observed between LSS and TIB and YNS sheep (Fisher's exact test, p < 2.2 × 10-16). Taken together, these findings offer insights crucial for the conservation and breeding enhancement of LSS.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Multi-character approach reveals a new mangrove population of the Yellow Warbler complex, Setophaga petechia, on Cozumel Island, Mexico

The Setophaga petechia complex includes 43 subspecies distributed within the new world, of which some are migratory and others are resident, with only two resident subspecies in the Mexican Caribbean: Setophaga petechia bryanti a mangrove subspecies belonging to the erithachorides group resident on the mainland of the Yucatan Peninsula and Setophaga petechia rufivertex endemic to Cozumel Island and belonging to the petechia group. Recently, a new population of individuals presenting intermediate phenotypic traits and living in mangrove ecosystems was discovered and reported for Cozumel Island. In this study, we used a multi-character approach including genetic (five ISSR genetic markers), morphometric (eight traits), phenotypic (four characteristics of males), and acoustic dataset (11 parameters) to understand the process of differentiation and the status of these new island individuals in relation to the two well-established subspecies using a total of 60 individuals (20 for each group). Through multivariate analyses based on different dataset used in our study, we show how the new population is related to the endemic island subspecies, S. p. rufivertex and to the mainland subspecies, S. p. bryanti while demonstrating finite differences. We conclude that the new population of S. petechia on Cozumel Island is a well-established population with high level of differentiation.

Genomes of the extinct Bachman's warbler show high divergence and no evidence of admixture with other extant Vermivora warblers

Bachman's warbler¹ (Vermivora bachmanii)-last sighted in 1988-is one of the only North American passerines to recently go extinct.^2,3,4 Given extensive ongoing hybridization of its two extant congeners-the blue-winged warbler (V. cyanoptera) and golden-winged warbler (V. chrysoptera)^5,6,7,8-and shared patterns of plumage variation between Bachman's warbler and hybrids between those extant species, it has been suggested that Bachman's warbler might have also had a component of hybrid ancestry. Here, we use historic DNA (hDNA) and whole genomes of Bachman's warblers collected at the turn of the 20^th century to address this. We combine these data with the two extant Vermivora species to examine patterns of population differentiation, inbreeding, and gene flow. In contrast to the admixture hypothesis, the genomic evidence is consistent with V. bachmanii having been a highly divergent, reproductively isolated species, with no evidence of introgression. We show that these three species have similar levels of runs of homozygosity (ROH), consistent with effects of a small long-term effective population size or population bottlenecks, with one V. bachmanii outlier showing numerous long ROH and a F_ROH greater than 5%. We also found-using population branch statistic estimates-previously undocumented evidence of lineage-specific evolution in V. chrysoptera near a pigmentation gene candidate, CORIN, which is a known modifier of ASIP, which is in turn involved in melanic throat and mask coloration in this family of birds. Together, these genomic results also highlight how natural history collections are such invaluable repositories of information about extant and extinct species.

Long divergent haplotypes introgressed from wild sheep are associated with distinct morphological and adaptive characteristics in domestic sheep

The worldwide sheep population comprises more than 1000 breeds. Together, these exhibit a considerable morphological diversity, which has not been extensively investigated at the molecular level. Here, we analyze whole-genome sequencing individuals of 1,098 domestic sheep from 154 breeds, and 69 wild sheep from seven Ovis species. On average, we detected 6.8%, 1.0% and 0.2% introgressed sequence in domestic sheep originating from Iranian mouflon, urial and argali, respectively, with rare introgressions from other wild species. Interestingly, several introgressed haplotypes contributed to the morphological differentiations across sheep breeds, such as a RXFP2 haplotype from Iranian mouflon conferring the spiral horn trait, a MSRB3 haplotype from argali strongly associated with ear morphology, and a VPS13B haplotype probably originating from urial and mouflon possibly associated with facial traits. Our results reveal that introgression events from wild Ovis species contributed to the high rate of morphological differentiation in sheep breeds, but also to individual variation within breeds. We propose that long divergent haplotypes are a ubiquitous source of phenotypic variation that allows adaptation to a variable environment, and that these remain intact in the receiving population probably due to reduced recombination.

Gut microbiome composition better reflects host phylogeny than diet diversity in breeding wood-warblers

Understanding the factors that shape microbiomes can provide insight into the importance of host-symbiont interactions and on co-evolutionary dynamics. Unlike for mammals, previous studies have found little or no support for an influence of host evolutionary history on avian gut microbiome diversity and instead have suggested a greater influence of the environment or diet due to fast gut turnover. Because effects of different factors may be conflated by captivity and sampling design, examining natural variation using large sample sizes is important. Our goal was to overcome these limitations by sampling wild birds to compare environmental, dietary and evolutionary influences on gut microbiome structure. We performed faecal metabarcoding to characterize both the gut microbiome and diet of 15 wood-warbler species across a 4-year period and from two geographical localities. We find host taxonomy generally explained ~10% of the variation between individuals, which is ~6-fold more variation of any other factor considered, including diet diversity. Further, gut microbiome similarity was more congruent with the host phylogeny than with host diet similarity and we found little association between diet diversity and microbiome diversity. Together, our results suggest evolutionary history is the strongest predictor of gut microbiome differentiation among wood-warblers. Although the phylogenetic signal of the warbler gut microbiome is not very strong, our data suggest that a stronger influence of diet (as measured by diet diversity) does not account for this pattern. The mechanism underlying this phylogenetic signal is not clear, but we argue host traits may filter colonization and maintenance of microbes.

Genetic confirmation of a hybrid between two highly divergent cardinalid species: A rose-breasted grosbeak (Pheucticus ludovicianus) and a scarlet tanager (Piranga olivacea)

Using low‐coverage whole‐genome sequencing, analysis of vocalizations, and inferences from natural history, we document a first‐generation hybrid between a rose‐breasted grosbeak (Pheucticus ludovicianus) and a scarlet tanager (Piranga olivacea). These two species occur sympatrically throughout much of eastern North America, although were not previously known to interbreed. Following the field identification of a putative hybrid, we use genetic and bioacoustic data to show that a rose‐breasted grosbeak was the maternal parent and a scarlet tanager was the paternal parent of the hybrid, whose song was similar to the latter species. These two species diverged >10 million years ago, and thus it is surprising to find a hybrid formed under natural conditions in the wild. Notably, the hybrid has an exceptionally heterozygous genome, with a conservative estimate of a heterozygous base every 100 bp. The observation that this hybrid of such highly divergent parental taxa has survived until adulthood serves as another example of the capacity for hybrid birds to survive with an exceptionally divergent genomic composition.

Estimating hybridization in the wild using citizen science data: a path forward

Genomic evidence of introgression in natural populations has reinvigorated the study of hybridization in recent years. Still, it is largely unknown how frequently individual organisms mate across species lines. Recently, Justyn et al. 2020 suggested that eBird, one of the world's largest citizen science databases, may supply adequate data for estimating hybridization rates. Here, we compare Justyn et al.'s estimates-and their conclusions that hybridization is rare-with estimates from museum and molecular data. We also estimate hybridization using eBird observations from areas and times when hybridization is possible, namely, in contact zones during the breeding season. These estimates are all considerably higher than those reported in Justyn et al., emphasizing that inferences from multiple datasets can differ radically. Finally, we demonstrate an approach for predicting the location of hybrid zones using eBird data, which can be done with high confidence and with unprecedented resolution. We show that citizen science data, far from settling the question of how frequently bird species hybridize, instead offer a promising step toward more focused study of hybrid zones. This article is protected by copyright. All rights reserved.

A multispecies BCO2 beak color polymorphism in the Darwin's finch radiation

Carotenoid-based polymorphisms are widespread in populations of birds, fish, and reptiles,¹ but generally little is known about the factors affecting their maintenance in populations.² We report a combined field and molecular-genetic investigation of a nestling beak color polymorphism in Darwin's finches. Beaks are pink or yellow, and yellow is recessive.³ Here we show that the polymorphism arose in the Galápagos half a million years ago through a mutation associated with regulatory change in the BCO2 gene and is shared by 14 descendant species. The polymorphism is probably a balanced polymorphism, maintained by ecological selection associated with survival and diet. In cactus finches, the frequency of the yellow genotype is correlated with cactus fruit abundance and greater hatching success and may be altered by introgressive hybridization. Polymorphisms that are hidden as adults, as here, may be far more common than is currently recognized, and contribute to diversification in ways that are yet to be discovered.

Signatures of mitonuclear coevolution in a warbler species complex

Divergent mitonuclear coadaptation could facilitate speciation. We investigate this possibility in two hybridizing species of warblers, Setophaga occidentalis and S. townsendi, in western North America. Inland S. townsendi harbor distinct mitochondrial DNA haplotypes from those of S. occidentalis. These populations also differ in several nuclear DNA regions. Coastal S. townsendi demonstrate mixed mitonuclear ancestry from S. occidentalis and inland S. townsendi. Of the few highly-differentiated chromosomal regions between inland S. townsendi and S. occidentalis, a 1.2 Mb gene block on chromosome 5 is also differentiated between coastal and inland S. townsendi. Genes in this block are associated with fatty acid oxidation and energy-related signaling transduction, thus linked to mitochondrial functions. Genetic variation within this candidate gene block covaries with mitochondrial DNA and shows signatures of divergent selection. Spatial variation in mitonuclear ancestries is correlated with climatic conditions. Together, these observations suggest divergent mitonuclear coadaptation underpins cryptic differentiation in this species complex.

Sexual selection and introgression in avian hybrid zones: Spotlight on Manacus

Hybrid zones offer a window into the processes and outcomes of evolution, from species formation or fusion to genomic underpinnings of specific traits and isolating mechanisms. Sexual selection is believed to be an important factor in speciation processes, and hybrid zones present special opportunities to probe its impact. The manakins (Aves, Pipridae) are a promising group in which to study the interplay of sexual selection and natural hybridization: they show substantial variation across the family in the strength of sexual selection they experience, they readily hybridize within and between genera, and they appear to have formed hybrid species, a rare event in birds. A hybrid zone between two manakins in the genus Manacus is unusual in that plumage and behavioral traits of one species have introgressed asymmetrically into populations of the second species through positive sexual selection, then apparently stalled at a river barrier. This is one of a handful of documented examples of asymmetric sexual trait introgression with a known selective mechanism. It offers opportunities to examine reproductive isolation, introgression, plumage color evolution, and natural factors enhancing or constraining the effects of sexual selection in real time. Here, we review previous work in this system, propose new hypotheses for observed patterns, and recommend approaches to test them.

Comparing divergence landscapes from reduced-representation and whole genome resequencing in the yellow-rumped warbler (Setophaga coronata) species complex

Researchers seeking to generate genomic data for non-model organisms are faced with a number of trade-offs when deciding which method to use. The selection of reduced representation approaches versus whole genome resequencing will ultimately affect the marker density, sequencing depth, and the number of individuals that can multiplexed. These factors can affect researchers' ability to accurately characterize certain genomic features, such as landscapes of divergence-how F_ST varies across the genomes. To provide insight into the effect of sequencing method on the estimation of divergence landscapes, we applied an identical bioinformatic pipeline to three generations of sequencing data (GBS, ddRAD, and WGS) produced for the same system, the yellow-rumped warbler species complex. We compare divergence landscapes generated using each method for the myrtle warbler (Setophaga coronata coronata) and the Audubon's warbler (S. c. auduboni), and for Audubon's warblers with deeply divergent mtDNA resulting from mitochondrial introgression. We found that most high-F_ST peaks were not detected in the ddRAD data set, and that while both GBS and WGS were able to identify the presence of large peaks, WGS was superior at a finer scale. Comparing Audubon's warblers with divergent mitochondrial haplotypes, only WGS allowed us to identify small (10-20 kb) regions of elevated differentiation, one of which contained the nuclear-encoded mitochondrial gene NDUFAF3. We calculated the cost per base pair for each method and found it was comparable between GBS and WGS, but significantly higher for ddRAD. These comparisons highlight the advantages of WGS over reduced representation methods when characterizing landscapes of divergence.

Dog colour patterns explained by modular promoters of ancient canid origin

Distinctive colour patterns in dogs are an integral component of canine diversity. Colour pattern differences are thought to have arisen from mutation and artificial selection during and after domestication from wolves but important gaps remain in understanding how these patterns evolved and are genetically controlled. In other mammals, variation at the ASIP gene controls both the temporal and spatial distribution of yellow and black pigments. Here, we identify independent regulatory modules for ventral and hair cycle ASIP expression, and we characterize their action and evolutionary origin. Structural variants define multiple alleles for each regulatory module and are combined in different ways to explain five distinctive dog colour patterns. Phylogenetic analysis reveals that the haplotype combination for one of these patterns is shared with Arctic white wolves and that its hair cycle-specific module probably originated from an extinct canid that diverged from grey wolves more than 2 million years ago. Natural selection for a lighter coat during the Pleistocene provided the genetic framework for widespread colour variation in dogs and wolves.