Genetic and behavioural factors affecting interpopulation colour pattern variation in two congeneric chameleon species.

We conducted a study on interpopulation variation of colour patterns in two congeneric chameleon species, which have an analogous life history. Both species are able to rapidly change colour pattern, and their context-dependent colour patterns often vary across a wide geographical range. Specifically, we tested four hypotheses that can explain the observed interpopulation variation of colour patterns by a series of behavioural field trials where the colour patterns of individuals were recorded and later analysed by a deep neural network algorithm. We used redundancy analysis to relate genetic, spectral and behavioural predictors to interpopulation colour pattern distance. Our results showed that both isolation by distance (IBD) and alternative mating tactics were significant predictors for interpopulation colour pattern variation in Chamaeleo chamaeleon males. By contrast, in Chamaeleo dilepis, the interpopulation colour pattern variation was largely explained by IBD, and evidence for alternative mating tactics was absent. In both chameleon species, the environmental colours showed no evidence of influencing chameleon interpopulation colour pattern variation, regardless of sex or behavioural context. This contrasting finding suggests that interpopulation context-dependent colour pattern variations in each species are maintained under a different set of selective pressures or circumstances.

Monitoring of species' genetic diversity in Europe varies greatly and overlooks potential climate change impacts.

Genetic monitoring of populations currently attracts interest in the context of the Convention on Biological Diversity but needs long-term planning and investments. However, genetic diversity has been largely neglected in biodiversity monitoring, and when addressed, it is treated separately, detached from other conservation issues, such as habitat alteration due to climate change. We report an accounting of efforts to monitor population genetic diversity in Europe (genetic monitoring effort, GME), the evaluation of which can help guide future capacity building and collaboration towards areas most in need of expanded monitoring. Overlaying GME with areas where the ranges of selected species of conservation interest approach current and future climate niche limits helps identify whether GME coincides with anticipated climate change effects on biodiversity. Our analysis suggests that country area, financial resources and conservation policy influence GME, high values of which only partially match species' joint patterns of limits to suitable climatic conditions. Populations at trailing climatic niche margins probably hold genetic diversity that is important for adaptation to changing climate. Our results illuminate the need in Europe for expanded investment in genetic monitoring across climate gradients occupied by focal species, a need arguably greatest in southeastern European countries. This need could be met in part by expanding the European Union's Birds and Habitats Directives to fully address the conservation and monitoring of genetic diversity.

Large variance in inbreeding within the Iberian wolf population.

The gray wolf (Canis lupus) population on the Iberian Peninsula was the largest in western and central Europe during most of the 20th century, with its size apparently never under a few hundred individuals. After partial legal protection in the 1970s in Spain, the northwest Iberian population increased to about 300-350 packs and then stabilized. In contrast to many current European wolf populations, which have been connected through gene flow, the Iberian wolf population has been isolated for decades. Here we measured changes on genomic diversity and inbreeding through the last decades in a geographic context. We find that the level of genomic diversity in Iberian wolves is low compared to other Eurasian wolf populations. Despite population expansion in the last 50 years, some modern wolves had very high inbreeding, especially in the recently recolonized and historical edge areas. These individuals contrast with others with low inbreeding within the same population. The high variance in inbreeding despite population expansion seems associated with small-scale fragmentation of the range that is revealed by the genetic similarity between modern and historical samples from close localities despite being separated by decades, remaining differentiated from other individuals that are just over 100 km away, a small distance for a species with great dispersal capacity inhabiting a continuous range. This illustrates that, despite its demographically stable condition, the population would probably benefit from favoring connectivity within the population as well as genetic exchange with other European wolf populations to avoid excessive fragmentation and local inbreeding depression.

Evolutionary history of an island endemic, the Azorean common quail.

Oceanic islands are characterized by conditions that favour diversification into endemic lineages that can be very different from their mainland counterparts. This can be the result of fast phenotypic divergence due to drift or the result of slower adaptation to local conditions. This uniqueness can obscure their evolutionary history. Here we used morphological, stable isotope, genetic and genomic data to characterize common quails (Coturnix coturnix) in the Azores archipelago and assess the divergence from neighbouring common quail populations. Historical documents suggested that these quails could have a recent origin associated with the arrival of humans in the last centuries. Our results show that Azorean quails constitute a well-differentiated lineage with small size and dark throat pigmentation that has lost the migratory ability and that diverged from mainland quail lineages more than 0.8 mya, contrary to the notion of a recent human-mediated arrival. Even though some Azorean quails carry an inversion that affects 115 Mbp of chromosome 1 and that has been associated with the loss of the migratory behaviour in other common quail populations, half of the analysed individuals do not have that inversion and still do not migrate. The long coexistence and evolution in isolation in the Azores of two chromosomal variants (with and without the inversion) is best explained by balancing selection. Thus, a unique and long evolutionary history led to the island endemic that we know today, C. c. conturbans.

High Investment in Reproduction Is Associated with Reduced Life Span in Dogs.

AbstractProminent differences in aging among and within species present an evolutionary puzzle. The theories proposed to explain evolutionary differences in aging are based on the axiom that selection maximizes fitness, not necessarily life span. This implies trade-offs between investment in self-maintenance and investment in reproduction, where high investments in growth and current reproduction are associated with short life spans. Fast growth and large adult size are related to shorter life spans in the domestic dog, a bourgeoning model in aging research; however, whether reproduction influences life span in this system remains unknown. Here we test the relationship between reproduction and differences in life span among dog breeds, simultaneously controlling for shared ancestry and recent gene flow. We found that shared ancestry explains a higher proportion of the among-breed variation in life history traits, in comparison with recent gene flow. Our results also show that reproductive investment negatively impacts life span, and more strongly so in large breeds, an effect that is not merely a correlated response of adult size. These results suggest that basic life history trade-offs are apparent in a domestic animal whose diversity is the result of artificial selection and that among-breed differences in life span are due to a combination of size and reproduction.

Diversification dynamics in Caribbean rain frogs (Eleutherodactylus) are uncoupled from the anuran community and consistent with adaptive radiation.

Adaptive radiation is proposed to play a key role in generating differences in species richness among lineages and geographical regions. Due to the importance of ecological divergence in adaptive radiation, species richness is predicted to be influenced by equilibrium diversity dynamics, although the concept continues to generate much debate. An additional important question is whether radiating clades have intrinsic biological characteristics that make them particularly prone to diversify. We tackle these questions by analysing (i) the temporal patterns of diversification of Caribbean Eleutherodactylus frogs, and (ii) assembly of the complete native anuran community of the Caribbean archipelago (197 species), testing for the presence of equilibrium dynamics and whether diversification patterns of Eleutherodactylus differ from those of the rest of the Caribbean anurans. Diversification rates follow the predicted pattern of rapid diversification early in the radiation which gradually decreases towards the present. Eleutherodactylus diversification is significantly faster than that of the Caribbean anuran community, and although equilibrium dynamics influence richness of all Caribbean anurans, Eleutherodactylus shows higher carrying capacity. Our results indicate that ecological opportunity per se is not sufficient for adaptive radiation and that diverse lineages present intrinsic characteristics that enable them to make the most of available opportunity.

The Phylogenetic Limits to Diversity-Dependent Diversification.

While the theory of micro-evolution by natural selection assigns a crucial role to competition, its role in macroevolution is less clear. Phylogenetic evidence for a decelerating accumulation of lineages suggests a feedback of lineage diversity on diversification. However, does this feedback only occur between close relatives, or do distant relatives also influence each other's diversification? In other words: are there phylogenetic limits to this diversity-dependence? Islands form ideal systems to answer these questions because their boundedness facilitates an overview of all potential competitors. The DAISIE (Dynamic Assembly of Island biota through Speciation Immigration and Extinction) framework allows for testing the presence of diversity-dependence on islands given phylogenetic data on colonization and branching times. The current inference models in DAISIE assume that this diversity-dependence only applies within a colonizing clade, i.e., all mainland species can colonize and diversify independently from one another. We term this clade-specific (CS) diversity-dependence. Here we introduce a new DAISIE model that assumes that diversity-dependence applies to all island species of a taxonomic group regardless of their mainland ancestry, i.e., diversity-dependence applies both to species within the same clade and between different clades established by different mainland species. We call this island-wide (IW) diversity-dependence. We present a method to compute a likelihood for this model given phylogenetic data on colonization and branching events and use likelihood ratio bootstrapping to compare it to the likelihood of the CS model in order to overcome biases known for standard model selection. We apply it to the diversification of Eleutherodactylus frogs on Hispaniola. Across the Greater Antilles archipelago, this radiation shows repeated patterns of diversification in ecotypes that are similar across clades. This could be suggestive of overlapping niche space and hence between-clade interactions, i.e., IW diversity-dependence. But it could also be suggestive of only within-clade interactions because between-clade interactions would have blocked the same ecotype from re-appearing. We find that the CS model fits the data much better than the IW model, indicating that different colonizations while resulting in similar ecotypes, are sufficiently distinct to avoid interacting strongly. We argue that non-overlapping distributions between clades (both spatially and in terms of ecotypes) cannot be used as evidence of CS diversity-dependence, because this pattern may be a consequence of IW diversity-dependence. By contrast, by using phylogenetic data rather than distributional data our method does allow for inferring the phylogenetic limits to diversity-dependent diversification. We discuss possibilities for future extensions and applications of our modelling approach. [Adaptive radiation; birth-death model; Caribbean; diversity-dependence; Eleutherodactylus; island biogeography.].

Influence of microhabitat, fecundity, and parental care on the evolution of sexual size dimorphism in Caribbean Eleutherodactylus frogs.

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.

Massive genome inversion drives coexistence of divergent morphs in common quails.

The presence of population-specific phenotypes often reflects local adaptation or barriers to gene flow. The co-occurrence of phenotypic polymorphisms that are restricted within the range of a highly mobile species is more difficult to explain. An example of such polymorphisms is in the common quail Coturnix coturnix, a small migratory bird that moves widely during the breeding season in search of new mating opportunities, following ephemeral habitats,1,2 and whose females may lay successive clutches at different locations while migrating.3 In spite of this vagility, previous studies reported a higher frequency of heavier males with darker throat coloration in the southwest of the distribution (I. Jiménez-Blasco et al., 2015, Int. Union Game Biol., conference). We used population genomics and cytogenetics to explore the basis of this polymorphism and discovered a large inversion in the genome of the common quail. This inversion extends 115 Mbp in length and encompasses more than 7,000 genes (about 12% of the genome), producing two very different forms. Birds with the inversion are larger, have darker throat coloration and rounder wings, are inferred to have poorer flight efficiency, and are geographically restricted despite the high mobility of the species. Stable isotope analyses confirmed that birds carrying the inversion have shorter migratory distances or do not migrate. However, we found no evidence of pre- or post-zygotic isolation, indicating the two forms commonly interbreed and that the polymorphism remains locally restricted because of the effect on behavior. This illustrates a genomic mechanism underlying maintenance of geographically structured polymorphisms despite interbreeding with a lineage with high mobility.

Loss of Mitochondrial Genetic Diversity despite Population Growth: The Legacy of Past Wolf Population Declines.

Gray wolves (Canis lupus) in the Iberian Peninsula declined substantially in both range and population size in the last few centuries due to human persecution and habitat fragmentation. However, unlike many other western European populations, gray wolves never went extinct in Iberia. Since the minimum number was recorded around 1970, their numbers have significantly increased and then stabilized in recent decades. We analyzed mitochondrial genomes from 54 historical specimens of Iberian wolves from across their historical range using ancient DNA methods. We compared historical and current mitochondrial diversity in Iberian wolves at the 5' end of the control region (n = 17 and 27) and the whole mitochondrial genome excluding the control region (n = 19 and 29). Despite an increase in population size since the 1970s, genetic diversity declined. We identified 10 whole mitochondrial DNA haplotypes in 19 historical specimens, whereas only six of them were observed in 29 modern Iberian wolves. Moreover, a haplotype that was restricted to the southern part of the distribution has gone extinct. Our results illustrate a lag between demographic and genetic diversity changes, and show that after severe population declines, genetic diversity can continue to be lost in stable or even expanding populations. This suggests that such populations may be of conservation concern even after their demographic trajectory has been reversed.

Automated genotyping of microsatellite loci from feces with high throughput sequences.

Ecological and conservation genetic studies often use noninvasive sampling, especially with elusive or endangered species. Because microsatellites are generally short in length, they can be amplified from low quality samples such as feces. Microsatellites are highly polymorphic so few markers are enough for reliable individual identification, kinship determination, or population characterization. However, the genotyping process from feces is expensive and time consuming. Given next-generation sequencing (NGS) and recent software developments, automated microsatellite genotyping from NGS data may now be possible. These software packages infer the genotypes directly from sequence reads, increasing throughput. Here we evaluate the performance of four software packages to genotype microsatellite loci from Iberian wolf (Canis lupus) feces using NGS. We initially combined 46 markers in a single multiplex reaction for the first time, of which 19 were included in the final analyses. Megasat was the software that provided genotypes with fewer errors. Coverage over 100X provided little additional information, but a relatively high number of PCR replicates were necessary to obtain a high quality genotype from highly unoptimized, multiplexed reactions (10 replicates for 18 of the 19 loci analyzed here). This could be reduced through optimization. The use of new bioinformatic tools and next-generation sequencing data to genotype these highly informative markers may increase throughput at a reasonable cost and with a smaller amount of laboratory work. Thus, high throughput sequencing approaches could facilitate the use of microsatellites with fecal DNA to address ecological and conservation questions.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients.

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

Biased assessment of ongoing admixture using STRUCTURE in the absence of reference samples.

Detection of hybridization and introgression is important in ecological research as in conservation and evolutionary biology. STRUCTURE is one of the most popular software to study introgression and allows estimating what proportion of the genome of each individual belongs to each ancestral population, even in cases where no reference sample from the ancestral nonadmixed populations is previously identified. In spite of its frequent use, some studies have indicated that ancestry estimates may not always be reliable. We simulated population data under different conditions with regard to the genetic differentiation between ancestral populations, number of loci considered, number of alleles per marker and hybridization rate, and analysed data with STRUCTURE. When reference samples were not included, the comparison of the known degree of admixture for each simulated individual and the value estimated with STRUCTURE revealed a strong underestimation of the level of introgression, classifying many admixed individuals as nonadmixed. This derives from an inaccurate estimation of the ancestral allele frequencies. When samples from the nonadmixed ancestral population were included as reference in the analyses, the bias in the estimations was reduced. The most accurate estimates were obtained when potentially admixed samples were few in relation to reference samples. Thus, whenever possible, a very large proportion of nonadmixed reference samples should be included in admixture assessments and different approaches should be combined. The misestimate of the amount of introgression can impair our understanding of the evolutionary history of species and misguide conservation efforts.

Vastly underestimated species richness of Amazonian salamanders (Plethodontidae: Bolitoglossa) and implications about plethodontid diversification.

We present data showing that the number of salamander species in Amazonia is vastly underestimated. We used DNA sequences of up to five genes (3 mitochondrial and 2 nuclear) of 366 specimens, 189 corresponding to 89 non-Amazonian nominal species and 177 Amazonian specimens, including types or topotypes, of eight of the nine recognized species in the region. By including representatives of all known species of Amazonian Bolitoglossa, except for one, and 73% of the currently 132 recognized species of the genus, our dataset represents the broadest sample of Bolitoglossa species, specimens, and geographic localities studied to date. We performed phylogenetic analyses using parsimony with tree-alignment and maximum likelihood (ML) with similarity alignment, with indels as binary characters. Our optimal topologies were used to delimit lineages that we assigned to nominal species and candidate new species following criteria that maximize the consilience of the current species taxonomy, monophyly, gaps in branch lengths, genetic distances, and geographic distribution. We contrasted the results of our species-delimitation protocol with those of Automated Barcode Gap Discovery (ABGD) and multi-rate Poisson Tree Processes (mPTP). Finally, we inferred the historical biogeography of South American salamanders by dating the trees and using dispersal-vicariance analysis (DIVA). Our results revealed a clade including almost all Amazonian salamanders, with a topology incompatible with just the currently recognized nine species. Following our species-delimitation criteria, we identified 44 putative species in Amazonia. Both ABGD and mPTP inferred more species than currently recognized, but their numbers (23-49) and limits vary. Our biogeographic analysis suggested a stepping-stone colonization of the Amazonian lowlands from Central America through the Chocó and the Andes, with several late dispersals from Amazonia back into the Andes. These biogeographic events are temporally concordant with an early land bridge between Central and South America (~10-15 MYA) and major landscape changes in Amazonia during the late Miocene and Pliocene, such as the drainage of the Pebas system, the establishment of the Amazon River, and the major orogeny of the northern Andes.

From groups to communities in western lowland gorillas.

Social networks are the result of interactions between individuals at different temporal scales. Thus, sporadic intergroup encounters and individual forays play a central role in defining the dynamics of populations in social species. We assessed the rate of intergroup encounters for three western lowland gorilla ( Gorilla gorilla gorilla) groups with daily observations over 5 years, and non-invasively genotyped a larger population over four months. Both approaches revealed a social system much more dynamic than anticipated, with non-aggressive intergroup encounters that involved social play by immature individuals, exchanges of members between groups likely modulated by kinship, and absence of infanticide evidenced by infants not fathered by the silverback of the group where they were found. This resulted in a community composed of groups that interacted frequently and not-aggressively, contrasting with the more fragmented and aggressive mountain gorilla ( G. beringei beringei) societies. Such extended sociality can promote the sharing of behavioural and cultural traits, but might also increase the susceptibility of western lowland gorillas to infectious diseases that have decimated their populations in recent times.

Ecomorphological convergence in Eleutherodactylus frogs: a case of replicate radiations in the Caribbean.

Replicate radiations, the repeated multiplication of species associated with ecological divergence, have attracted much attention and generated as much debate. Due to the few well-studied cases, it remains unclear whether replicate radiations are an exceptional result of evolution or a relatively common example of the power of adaptation by natural selection. We examined the case of Eleutherodactylus frogs, which radiated in the Caribbean islands resulting in more than 160 species that occupy very diverse habitats. A time-calibrated phylogeny revealed that these frogs independently diversified on all larger islands producing species that occupy a broad range of microhabitats in different islands. Using phylogenetic comparative methods, we found an association between morphological traits and particular microhabitats, and for most microhabitats detected significant morphological convergence. Our results indicate Caribbean Eleutherodactylus are a novel example of replicate radiations, and highlight the predictability of evolutionary processes, as similar ecological opportunities can lead to similar outcomes.

The genomic basis of adaptation to high-altitude habitats in the eastern honey bee (Apis cerana).

The eastern honey bee (Apis cerana) is of central importance for agriculture in Asia. It has adapted to a wide variety of environmental conditions across its native range in southern and eastern Asia, which includes high-altitude regions. eastern honey bees inhabiting mountains differ morphologically from neighbouring lowland populations and may also exhibit differences in physiology and behaviour. We compared the genomes of 60 eastern honey bees collected from high and low altitudes in Yunnan and Gansu provinces, China, to infer their evolutionary history and to identify candidate genes that may underlie adaptation to high altitude. Using a combination of FST -based statistics, long-range haplotype tests and population branch statistics, we identified several regions of the genome that appear to have been under positive selection. These candidate regions were strongly enriched for coding sequences and had high haplotype homozygosity and increased divergence specifically in highland bee populations, suggesting they have been subjected to recent selection in high-altitude habitats. Candidate loci in these genomic regions included genes related to reproduction and feeding behaviour in honey bees. Functional investigation of these candidate loci is necessary to fully understand the mechanisms of adaptation to high-altitude habitats in the eastern honey bee.

On the path to extinction: Inbreeding and admixture in a declining grey wolf population.

Allee effects reduce the viability of small populations in many different ways, which act synergistically to lead populations towards extinction vortexes. The Sierra Morena wolf population, isolated in the south of the Iberian Peninsula and composed of just one or few packs for decades, represents a good example of how diverse threats act additively in very small populations. We sequenced the genome of one of the last wolves identified (and road-killed) in Sierra Morena and that of another wolf in the Iberian Wolf Captive Breeding Program and compared them with other wolf and dog genomes from around the world (including two previously published genome sequences from northern Iberian wolves). The results showed relatively low overall genetic diversity in Iberian wolves, but diverse population histories including past introgression of dog genes. The Sierra Morena wolf had an extraordinarily high level of inbreeding and long runs of homozygosity, resulting from the long isolation. In addition, about one-third of the genome was of dog origin. Despite the introgression of dog genes, heterozygosity remained low because of continued inbreeding after several hybridization events. The results thus illustrate the case of a small and isolated wolf population where the low population density may have favoured hybridization and introgression of dog alleles, but continued inbreeding may have resulted in large chromosomal fragments of wolf origin completely disappearing from the population, and being replaced by chromosomal fragments of dog origin. The latest population surveys suggest that this population may have gone extinct.