Find research datasets worth reusing

Foraging behaviour variations, for gene expression and transcriptomic divergence in parasitic wasp populations of Venturia canescens

Foraging behaviours encompass strategies to locate resources and to exploit them. In many taxa, these behaviours are driven by a major gene called for, but mechanisms vary between species. In the parasitoid wasp Venturia canescens, sexual and asexual populations coexist in sympatry but differ in life-history trait, physiology and behaviours, which could impact their foraging strategies. Here, we explored the molecular bases underpinning divergence in behaviours by testing two mutually nonexclusive hypotheses: first, the divergence in the for gene correlates with difference in foraging strategies, and second, the latter rely on a divergence in whole-genome expression. Using comparative genomics, we showed that the for gene was conserved across insects considering both sequence and gene model complexity. Polymorphism analysis did not support the occurrence of two allelic variants diverging across the two populations, yet the asexual population exhibited less polymorphism than the sexual population. Sexual and asexual transcriptomes sharply split, with 10.9% of differentially expressed genes, but these were not enriched in behavioural-related genes. We showed that the for gene was more highly expressed in asexual female heads than in sexual heads and that those differences correlate with divergence in foraging behaviours in our experiment since asexuals explored the environment more and exploited more host patches. Overall, these results suggested that fine tuning of for gene expression between populations may have led to distinct foraging behaviours. We hypothesized that reproductive polymorphism and coexistence in sympatry of sexual and asexual populations specialized to different ecological niches via divergent optima on phenotypic traits could imply adaptation through different expression patterns of the for gene and at many other loci throughout the genome.

Dataset for: Breeding populations of Marbled Godwits and Willets have high annual survival and strong site fidelity to managed wetlands

The Prairie Pothole Region of central Canada supports a diverse community of breeding waterbirds but many species have declining populations and the demographic mechanisms driving the declines remain unknown. We conducted a 7-year field study during 1995–2001 to investigate the demographic performance of Marbled Godwits (Limosa fedoa) and Willets (Tringa semipalmata) breeding in managed wetlands near Brooks, Alberta. Mark-recapture analyses based on Cormack-Jolly-Seber models revealed that the annual rates of apparent survival for Marbled Godwits (phi-hat = 0.934 ± 0.015SE) and Willets (phi-hat = 0.849 ± 0.018SE) are among the highest rates of survivorship reported for any breeding or nonbreeding population of large-bodied shorebirds. Our estimates of life expectancy for males were comparable to longevity records in godwits (E-hat = 17.3 years ± 5.8SE vs. 25-29+ years) and Willets (E- hat = 7.7 ± 1.5SE vs. 10+ years). The two species both showed strong breeding site fidelity but differed in rates of mate fidelity. Pairs that reunited and males that switched mates usually nested <300 m from their previous nests whereas females that switched mates usually moved longer distances >1.1 to 1.5 km. Returning pairs usually reunited in godwits (85%) but not in Willets (28%), possibly because of species differences in adult survival or patterns of migration. Baseline estimates of annual survival for banded-only birds will be useful for evaluating the potential effects of new tracking tags or the environmental changes that have occurred during the past 20 years. Conservation strategies for large-bodied shorebirds should be focused on reduction of exposure to anthropogenic mortality because low rates of natural mortality suggest that losses to collisions at breeding sites or harvest at nonbreeding areas are likely to cause additive mortality.

opencc-by-4.0Dec 2022View details →

zenodo36/100

Effects of the environmental conditions and seasonality on a population survey of the Andean Condor Vultur gryphus in the tropical Andes

Script, datasets, and other supplementary material in Márquez-Alvis et al. 2023 PeerJ. Background: Among the New World vultures, the Andean Condor is considered one of the most culturally and ecologically important species. However, their populations are declining over their entire distributional range. In response, conservation strategies have been implemented in many countries to reverse the increasing extinction risk of this species. The initiatives rely on extensive population surveys to gather basic information necessary to implement policies and to intervene efficiently. Still, there is a need to standardize the surveys based on seasonality and suitable environmental conditions throughout the species distribution.  Here, we provide the first assessment of how daily temperature, rainfall, and seasonality influence surveys of Andean Condors on a communal roost in the central Peruvian Andes. Methods: Using an autoregressive generalized linear model, we associated environmental variables with visual surveys of adult and young condors at three different times of the day and three times a week between June 2014 and March 2015. Results: We found that both adults and young Andean Condors showed a threefold reduction in the use of the communal roost after the beginning of the rainy season. Colder and drier days (dry season) are preferable for surveying, as we expect the total number of condors using communal roosts to reduce under rainy (rainfall = -0.53 ± 0.16) and warmer days (temperature= - 0.04 ± 0.02) days. Therefore, the significant variation in the use of roosts across seasons and hours should be carefully accounted for in national surveys, at the risk of undermining the full potential of the communal roost surveys. Moreover, we also found a strong bias towards immatures (about 76%) in the adult:immature ratio and a remarkable absence of Andean Condors during the wet season. These results suggest that the species might be using other unknown communal roosts hierarchically. Such results provide key information for selecting priority areas for conservation and selecting the best time to survey this species in the tropical Andes. Finally, it may open a fruitful avenue for further research on the protection of the Andean Condor.

Genetic population dynamics of the critically endangered scalloped hammerhead shark (Sphyrna lewini) in the Eastern Tropical Pacific

The scalloped hammerhead shark, Sphyrna lewini, is a Critically Endangered, migratory species known for its tendency to form iconic and visually spectacular large aggregations. Herein, we investigated the population genetic dynamics of the scalloped hammerhead across much of its distribution in the Eastern Tropical Pacific (ETP), ranging from Costa Rica to Ecuador, focusing on young‐of‐the‐year animals from putative coastal nursery areas and adult females from seasonal aggregations that form in the northern Galápagos Islands. Nuclear microsatellites and partial mitochondrial control region sequences showed little evidence of population structure, suggesting that scalloped hammerheads in this ETP region comprise a single genetic stock. Galápagos aggregations of adults were not comprised of related individuals, suggesting that kinship does not play a role in the formation of the repeated, annual gatherings at these remote offshore locations. Despite high levels of fisheries exploitation of this species in the ETP, the adult scalloped hammerheads here showed greater genetic diversity compared with adult conspecifics from other parts of the species' global distribution. A phylogeographic analysis of available, globally sourced, mitochondrial control region sequence data (n = 1,818 sequences) revealed that scalloped hammerheads comprise three distinct matrilines corresponding to the three major world ocean basins, highlighting the need for conservation of these evolutionarily unique lineages. This study provides the first view of the genetic properties of a scalloped hammerhead aggregation, and the largest sample size‐based investigation of population structure and phylogeography of this species in the ETP to date.

Data and MATLAB files for: Timescale analyses of fluctuations in coexisting populations of a native and invasive tree squirrel

1. Competition from invasive species is an increasing threat to biodiversity. In Southern California, the western gray squirrel (Sciurus griseus, WGS) is facing increasing competition from the fox squirrel (Sciurus niger, FS), an invasive congener. 2. We used spectral methods to analyze 140 consecutive monthly censuses of WGS and FS within a 11.3 ha section of the California Botanic Garden. Variation in the numbers for both species and their synchrony was distributed across long timescales (> 15 months). 3. After filtering out annual changes, concurrent mean monthly temperatures from nearby Ontario Airport (ONT) yielded a spectrum with a large semiannual peak and significant spectral power at long timescales (> 30 months). Squirrel-temperature cospectra showed significant negative covariation at long timescales (> 35 months) for WGS and smaller significant negative peaks at 6 months for both species. 4. Simulations from a Lotka-Volterra model of two competing species indicates that the risk of extinction for the weaker competitor increases quickly as environmental noise shifts from short to long timescales. 5. We analyzed the timescales of fluctuations in detrended mean annual temperatures for the time period 1915-2014 from 1218 locations across the continental USA. In the last two decades, significant shifts from short timescales to long timescales have occurred, changing from less than 3 years to 4-6 years. 6. Our results indicate that (i) population fluctuations in co-occurring native and invasive tree squirrels are synchronous, occur over long timescales, and may be driven by fluctuations in environmental conditions; (ii) long timescale population fluctuations increase the risk of extinction in competing species, especially for the inferior competitor; and (iii) the timescales of interannual environmental fluctuations may be increasing from recent historical values. These results have broad implications for the impact of climate change on the maintenance of biodiversity.

Methodological challenges in the genomic analysis of an endangered mammal population with low genetic diversity

Recently, populations of various species with very low genetic diversity have been discovered. Some of these persist in the long term, but others could face extinction due to accelerated loss of fitness. In this work, we characterize 45 individuals of one of these populations, belonging to the Iberian desman (Galemys pyrenaicus). For this, we used the ddRADseq technique, which generated 1,421 SNPs. The heterozygosity values of the analyzed individuals were among the lowest recorded for mammals, ranging from 26 to 91 SNPs/Mb. Furthermore, the individuals from one of the localities, highly isolated due to strong barriers, presented extremely high inbreeding coefficients, with values above 0.7. Under this scenario of low genetic diversity and elevated inbreeding levels, some individuals appeared to be almost genetically identical. We used different methods and simulations to determine if genetic identification and parentage analysis were possible in this population. Only one of the methods, which does not assume population homogeneity, was able to identify all individuals correctly. Therefore, genetically impoverished populations pose a great methodological challenge for their genetic study. However, these populations are of primary scientific and conservation interest, so it is essential to characterize them genetically and improve genomic methodologies for their research.

Flow cytometry YFP and CFP data and deep sequencing data of populations evolving in galactose

Copy-number and point mutations form the basis for most evolutionary novelty through the process of gene duplication and divergence. While a plethora of genomic sequence data reveals the long-term fate of diverging coding sequences and their cis-regulatory elements, little is known about the early dynamics around the duplication event itself. In microorganisms, selection for increased gene expression often drives the expansion of gene copy-number mutations, which serves as a crude adaptation, prior to divergence through refining point mutations. Using a simple synthetic genetic system that allows us to distinguish copy-number and point mutations, we study their early and transient adaptive dynamics in real-time in Escherichia coli. We find two qualitatively different routes of adaptation depending on the level of functional improvement selected for: In conditions of high gene expression demand, the two types of mutations occur as a combination. Under low gene expression demand, negative epistasis between the two types of mutations renders them mutually exclusive. Thus, owing to their higher frequency, adaptation is dominated by copy-number mutations. Ultimately, due to high rates of reversal and pleiotropic cost, copy-number mutations may not only serve as a crude and transient adaptation but also <a>constrain</a> sequence divergence over evolutionary time scales.

Divergent foraging strategies between populations of sympatric matrilineal killer whales

In gregarious species, collective behavior maximizes individual fitness benefits while minimizing costs. Despite the relevance of behavior to conservation, the link between the robustness of behavioral patterns across populations and population health is poorly understood. We studied the collective foraging behavior of two closed, sympatric populations of piscivorous killer whales, leveraging two contemporaneously-collected data sets from suction cup-attached bio-logging tags, to quantify patterns of fine-scale foraging movements and their relationships with demography. We reveal striking plasticity in collective foraging behavior between populations. Prey capture rate and foraging efficiency were greater for males in the endangered Southern Resident (SRKW) population, yet greater for females in the Northern Resident (NRKW) population. The presence of a calf (≤ 3 y) reduced the number of prey captured by adult females in both populations, with the greatest effect in SRKW, in which no mothers with calves captured prey while bearing tags. SRKW adult males with a living mother tended to capture more prey than those whose mother had died, whereas the opposite was true for NRKW adult males. Moreover, males generally tended to forage in areas with deeper bathymetry than females, and SRKW captured prey deeper than NRKW. These population-level differences in sex-specific foraging behavior challenge the existing paradigm that mothers are disproportionate foragers in gregarious killer whales, underscoring that an endangered population is employing a potentially unstable collective foraging strategy mismatched to recovery. Thus, our study provides a mechanistic link between fine-scale, collective foraging behavior and population health in an apex marine predator.

Data for: Species identification and population genetics of the Antarctic fish genera Lepidonotothen and Nototheniops (Perciformes, Notothenioidei)

Accurate species identification is essential to assess biodiversity and species richness in ecosystems threatened by rapid and recent environmental changes, such as warming in most Antarctic waters. The Lepidonotothen species complex comprises demersal notothenioid fishes which inhabit the shelf areas of the Antarctic Peninsula, the Scotia Arc and sub-Antarctic islands with a circum-Antarctic distribution. Species determination in this group has often been problematic. In particular, whether Lepidonotothen squamifrons and Lepidonotothen kempi are valid as separate species has been questioned. In this study, we analysed the genetic variation among four nominal southern polar species within this complex (L. kempi, L. squamifrons, Nototheniops larseni, Nototheniops nudifrons) by means of three different markers (ND2 and tRNA mitochondrial genes and a panel of 16 nuclear microsatellites). We tested whether individuals morphologically assigned to L. kempi showed genetic separation from L. squamifrons. Our analyses indicated a lack of differentiation between L. kempi and L. squamifrons. However, a genetically distinct population was found for L. squamifrons at the Shag Rocks islands near South Georgia. Antarctic and sub-Antarctic islands are known to be home to many cryptic species and further studies will elucidate if the genetically differentiated population we found potentially originated from this context and can be considered as an incipient species. Our analysis contributes to further characterize the species composition of the most abundant fish suborder in the Southern Ocean, which is amongst the regions most threatened by climate change.

TE invasion fuels molecular adaptation in laboratory populations of Drosophila melanogaster

Transposable elements are mobile genetic parasites that frequently invade new host genomes through horizontal transfer. Invading TEs often exhibit a burst of transposition, followed by reduced transposition rates as repression evolves in the host. We recreated the horizontal transfer of P-element DNA transposons into a D. melanogaster host and followed the expansion of TE copies and evolution of host repression in replicate laboratory populations reared at different temperatures. We observed that while populations maintained at high temperatures rapidly go extinct after TE invasion, those maintained at lower temperatures persist, allowing for TE spread and the evolution of host repression. We also surprisingly discovered that invaded populations experienced recurrent insertion of P-elements into a specific long non-coding RNA, lncRNA:CR43651, and that these insertion alleles are segregating at unusually high frequency in experimental populations, indicative of positive selection. We propose that, in addition to driving the evolution of repression, transpositional bursts of invading TEs can drive molecular adaptation.

Population admixtures in medaka inferred by multiple arbitrary amplicon sequencing

Cost-effective genotyping can be achieved by sequencing PCR amplicons. Short 3–10 base primers can arbitrarily amplify thousands of loci using only a few primers. To improve the sequencing efficiency of the multiple arbitrary amplicon sequencing (MAAS) approach, we designed new primers and examined their efficiency in sequencing and genotyping. To demonstrate the effectiveness of our method, we applied it to examine the population structure of the small freshwater fish, medaka (Oryzias latipes). We obtained 2,987 informative SNVs with no missing genotype calls for 67 individuals from 15 wild populations and three artificial strains. The estimated phylogenic and population genetic structures of the wild populations were consistent with previous studies, corroborating the accuracy of our genotyping method. We also attempted to reconstruct the genetic backgrounds of a commercial orange mutant strain, Himedaka, which has caused a genetic disturbance in wild populations. Our admixture analysis focusing on Himedaka showed that at least two wild populations had genetically contributed to the nuclear genome of this mutant strain. Our genotyping methods and results will be useful in quantitative assessments of genetic disturbance by this commercially available strain.

Genome-phenotype-environment associations identify signatures of selection in a panmictic population of threespine stickleback

Adaptive genetic divergence occurs when selection imposed by the environment causes the genomic component of the phenotype to differentiate. However, genomic signatures of natural selection are usually identified without information on which trait is responding to selection by which selective agent(s). Here we integrate whole-genome-sequencing with phenomics and measures of putative selective agents to assess the extent of adaptive divergence in threespine stickleback occupying the highly heterogeneous lake Mývatn, NE Iceland. We find negligible genome-wide divergence, yet multiple traits (body size, gill raker structure and defence traits) were divergent along known ecological gradients (temperature, predatory bird densities and water depth). SNP-based heritability of all measured traits was high (h2 = 0.42 – 0.65), indicating adaptive potential for all traits. Whilst environment-association analyses identified thousands of loci putatively involved in selection, related to genes linked to neuron development and protein phosphorylation, only loci linked to pelvic spine length were concurrently linked to environmental variation (water depth) – supporting the conclusion that divergence in pelvic spine length occurred in face of gene flow. Our results suggest that whilst there is substantial genetic variation in the traits measured, phenotypic divergence of Mývatn stickleback is mostly weakly associated with environmental gradients, potentially as a result of substantial gene flow. Our study illustrates the value of integrative studies that combine genomic assays of multivariate trait variation with landscape genomics.

Genetic assessment reveals inbreeding, possible hybridization, and low levels of genetic structure in a declining goose population

The population numbers of taiga bean goose (Anser fabalis fabalis) have halved during recent decades. Since this subspecies is hunted throughout most of its range, the decline is of management concern. Knowledge of the genetic population structure and diversity is important for guiding management and conservation efforts. Genetically unique subpopulations might be hunted to extinction if not managed separately, and any inbreeding depression or lack of genetic diversity may affect the ability to adapt to changing environments and increase extinction risk. We used microsatellite and mitochondrial DNA markers to study the genetic population structure and diversity among taiga bean geese breeding within the Central flyway management unit using non-invasively collected feathers. We found some genetic structuring with the maternally inherited mitochondrial DNA between four geographic regions (ɸST = 0.11-0.20) but none with the nuclear microsatellite markers (all pairwise FST-values 0.002- 0.005). These results could be explained by female natal philopatry and male-biased dispersal, which completely homogenizes the nuclear genome. Therefore, the population could be managed as a single unit. Genetic diversity was still at a moderate level (average HE = 0.69) and there were no signs of past population size reductions, although significantly positive inbreeding coefficients in all sampling sites (FIS = 0.05-0.10) and high relatedness values (r = 0.60-0.86) between some individuals could indicate inbreeding. In addition, there was evidence of either incomplete lineage sorting or introgression from the pink-footed goose (A. brachyrhynchus). The current population is not under threat by genetic impoverishment but monitoring in the future is desirable.

Genomic analyses point to a low evolutionary potential of prospective source populations for assisted migration in a forest herb

Climate change is increasingly impacting temperate forest ecosystems and many forest herbs might be unable to track the changing climate due to dispersal limitation. Forest herbs with a low adaptive capacity may therefore benefit from conservation strategies that mitigate dispersal limitation and evolutionary constraints, such as assisted migration. However, assisted migration strategies rarely consider evolutionary constraints of potential source populations that may jeopardize their success. In cases where climate adaptation is overshadowed by competing evolutionary processes, assisted migration is unlikely to support adaptation to future climates. Using a combination of population and landscape genomic analyses, we disentangled local adaptation drivers and quantified the adaptability and vulnerability to climate change of the self-incompatible deciduous forest herb Primula elatior. Southern populations displayed a sharp genetic turnover, and a considerable amount of local adaptation under diversifying selection was discovered. However, most of the outlier loci could not be linked to climate variables (71%) and were likely related to other local adaptation drivers, such as photoperiodism. Furthermore, specific adaptations to climate extremes, such as drought stress, could not be detected. This is in line with the typical occurrence of forest herbs in buffered climatic conditions, which can be expected to reduce selection pressures imposed by climate. Finally, populations in the south of the distribution area had increased sensitivity to climate change due to a reduced adaptive capacity and a moderate genetic offset, while central European populations were sensitive due to a high genetic offset. We conclude that assisted migration from southern source populations could bear significant risk due to nonclimatic maladaptation and a low adaptive capacity. Regional admixture and restoration of ecological connectivity to increase the adaptive capacity, and assisted range expansion to suitable habitat in the north might be more appropriate mitigation strategies.

Raising offspring increases ageing: Differences in senescence among three populations of a long-lived seabird, the Atlantic puffin

1. Actuarial senescence, the decline of survival with age, is well documented in the wild. Rates of senescence vary widely between taxa, to some extent also between sexes, with the fastest life histories showing the highest rates of senescence. Few studies have investigated differences in senescence among populations of the same species, although such variation is expected from population-level differences in environmental conditions, leading to differences in vital rates and thus life histories. 2. We predict that, within species, populations differing in productivity (suggesting different paces of life) should experience different rates of senescence, but with little or no sexual difference in senescence within populations of monogamous, monomorphic species where the sexes share breeding duties. 3. We compared rates of actuarial senescence among three contrasting populations of the Atlantic puffin Fratercula arctica. The data set comprised 31 years (1990–2020) of parallel capture-mark-recapture data from three breeding colonies, Isle of May (North Sea), Røst (Norwegian Sea) and Hornøya (Barents Sea), showing contrasting productivities (i.e. annual breeding success) and population trends. We used time elapsed since first capture (TFC) as a proxy for bird age, and productivity and the winter North Atlantic Oscillation Index (wNAO) as proxies for the environmental conditions experienced by the populations within and outside the breeding season, respectively. 4. In accordance with our predictions, we found that senescence rates differed among the study populations, with no evidence for sexual differences. There was no evidence for an effect of wNAO, but the population with the lowest productivity, Røst, showed the lowest rate of senescence. As a consequence, the negative effect of senescence on the population growth rate (λ) was up to 3–5 times smaller on Røst (Δλ = -0.009) than on the two other colonies. 5. Our findings suggest that environmentally induced differences in senescence rates among populations of a species should be accounted for when predicting effects of climate variation and change on species persistence. There is thus a need for more detailed information on how both actuarial and reproductive senescence influence vital rates of populations of the same species, calling for large-scale comparative studies.

Data from: Population connectivity patterns of genetic diversity, immune responses and exposure to infectious pneumonia in a metapopulation of desert bighorn sheep

Habitat fragmentation is an important driver of biodiversity loss and can be remediated through management actions aimed at maintenance of natural connectivity in metapopulations. Connectivity may protect populations from infectious diseases by preserving immunogenetic diversity and disease resistance. However, connectivity could exacerbate the risk of infectious disease spread across vulnerable populations. We tracked the spread of a novel strain of Mycoplasma ovipneumoniae in a metapopulation of desert bighorn sheep (Ovis canadensis nelsoni) in the Mojave desert to investigate how variation in connectivity among populations influenced disease outcomes. M. ovipneumoniae was detected throughout the metapopulation, indicating that the relative isolation of many of these populations did not protect them from pathogen invasion. However, we show that connectivity among bighorn sheep populations was correlated with higher immunogenetic diversity, a protective immune response and lower disease prevalence. Variation in protective immunity predicted infection risk in individual bighorn sheep and was associated with heterozygosity at genetic loci linked to adaptive and innate immune signalling. Together, these findings may indicate that population connectivity maintains immunogenetic diversity in bighorn sheep populations in this system and has direct effects on immune responses in individual bighorn sheep and their susceptibility to infection by a deadly pathogen. Our study suggests that the genetic benefits of population connectivity could outweigh the risk of infectious disease spread and supports conservation management that maintains natural connectivity in metapopulations.

The population genetics of speciation by cascade reinforcement

Species interactions drive diverse evolutionary outcomes. Speciation by cascade reinforcement represents one example of how species interactions can contribute to the proliferation of species. This process occurs when the divergence of mating traits in response to selection against interspecific hybridization incidentally leads to reproductive isolation among populations of the same species. Here, we investigated the population genetic outcomes of cascade reinforcement in North American chorus frogs (Hylidae: Pseudacris). Specifically, we estimated the frequency of hybridization among three taxa, assessed genetic structure within the focal species, P. feriarum, and ascertained the directionality of gene flow within P. feriarum across replicated contact zones via coalescent modeling. Through field observations and preliminary experimental crosses, we assessed whether hybridization is possible under natural and laboratory conditions. We found that hybridization occurs among P. feriarum and two conspecifics at a low rate in multiple contact zones and that gene flow within the former species is unidirectional from allopatry into sympatry with these other species in three of four contact zones studied. We found evidence of substantial genetic structuring within P. feriarum including a divergent western allopatric cluster, a behaviorally-distinct sympatric South Carolina cluster, and several genetically-overlapping clusters from the remainder of the distribution. Furthermore, we found sub-structuring between reinforced and non-reinforced populations in the two most intensely-sampled contact zones. Our literature review indicated that P. feriarum hybridizes with at least five heterospecifics at the periphery of its range providing a mechanism for further intraspecific diversification. This work strengthens the evidence for cascade reinforcement in this clade, revealing the geographic and genetic landscape upon which this process can contribute to the proliferation of species.

Data from: Female and male plumage brightness is positively correlated among populations of the dichromatic Variable Seedeater

Because studies of plumage coloration often focus on highly elaborate and conspicuous males, little is known about the factors driving the evolution of drab plumage in females of dichromatic species. We explored variation in the color properties of drab female plumage across subspecies of the Variable Seedeater (Sporophila corvina), and found a correlation between plumage brightness of male and female plumage patches. This correlation is consistent with the hypothesis that differences in plumage brightness among females result, at least partially, from a genetic correlation or shared developmental pathway with the more conspicuous male plumage traits that are shaped by sexual selection. Further, we modeled the perception of female plumage using the avian visual system and found that most of the variation in color among females is mostly undetectable by birds; therefore, it is unlikely for such variation to have evolved via direct sexual selection as a signal for recognition or assortative mating. Together, these results suggest that variation in S. corvina female color result from non-adaptive processes. However, we found a reduced variation on color of dorsal patches compared with that of ventral patches among females, which could be explained by natural selection operation on crypsis, reducing the variation on dorsal patches. Finally, we found little differences in plumage color between adult females and juveniles, which suggest selection for juveniles to resemble adult females possibly to reduce aggression from adult males. Additional within-population studies and behavioral experiments are needed to discern the evolutionary factors driving variation on S. corvina female color and its function.