Negative effects of the spatial clumping of thermal resources on lizard thermoregulation in a fragmented habitat.

In ecosystems threatened by the expansion of croplands, habitat fragmentation and climate change, two of the main extinction drivers, may have thermoregulation-mediated interacting effects on demographic trends of terrestrial ectotherms. We studied the thermal biology of a metapopulation of the widespread Mediterranean lacertid Psammodromus algirus in ten fragments of evergreen or deciduous oak forests interspersed among cereal fields. We obtained thermoregulation statistics (selected temperature range, body and operative temperatures, thermal quality of the habitat, and precision, accuracy, and effectiveness of thermoregulation) that could be compared among fragments and with conspecific populations living in unfragmented habitat. We also measured the selection (use vs. availability) and spatial distribution of sunlit and shaded patches used for behavioral thermoregulation in fragments, and we estimated operative temperatures and thermal habitat quality in the agricultural matrix surrounding the fragments. Variation of the thermal environment was much larger within fragments than among them, and thermoregulation was accurate, precise, and efficient throughout the fragmented landscape; its effectiveness was similar to that of previously studied unfragmented populations. The average distance between sunlit and shaded patches was shorter in deciduous than in evergreen fragments, producing a more clumped distribution of the mosaic of thermal resources. Consequently, in evergreen habitat the cost of thermoregulation was higher, because lizards were more selective in their choice of sunlit sites (i.e. they used sunlit patches closer to shade and refuge than expected at random, and the extent of such selection was larger than at deciduous habitat). Temperatures available in croplands were too high to allow lizard dispersal, at least in the post-breeding season. This result confirms the role of croplands as a thermal barrier that promotes inbreeding and associated fitness losses in isolated fragments, and it forecasts a dark future for populations of forest lizards in agricultural landscapes under the combined effects of habitat fragmentation and global warming.

Environmental association modelling with loci under divergent selection predicts the distribution range of a lizard.

During the historical building of a species range, individual colonizers have to confront different ecological challenges, and the capacity of the species to broaden its range may depend on the total amount of adaptive genetic variation supplied by evolution. We set out to increase our understanding of what defines a range and the role of underlying genetics by trying to predict an entire species' range from the geographical distribution of its genetic diversity under selection. We sampled five populations of the western Mediterranean lizard Psammodromus algirus that inhabit a noticeable environmental gradient of temperature and precipitation. We correlated the genotypes of 95 individuals (18-20 individuals per population) for 21 SNPs putatively under selection with environmental scores on a bioclimatic gradient, using 1 × 1 km2  grid cells as sampling units. By extrapolating the resulting model to all possible combinations of alleles, we inferred all the geographic cells that were theoretically suitable for a given amount of genetic variance under selection. The inferred distribution range overlapped to a large extent with the realized range of the species (77.46% of overlap), including an accurate prediction of internal gaps and range borders. Our results suggest an adaptability threshold determined by the amount of genetic variation available that would be required to warrant adaptation beyond a certain limit of environmental variation. These results support the idea that the expansion of a species' range can be ultimately linked to the arising of new variants under selection (either newly selected variants from standing genetic variation or innovative mutations under selection).

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species.

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

Low genome-wide divergence between two lizard populations with high adaptive phenotypic differentiation.

Usually, adaptive phenotypic differentiation is paralleled by genetic divergence between locally adapted populations. However, adaptation can also happen in a scenario of nonsignificant genetic divergence due to intense gene flow and/or recent differentiation. While this phenomenon is rarely published, findings on incipient ecologically driven divergence or isolation by adaptation are relatively common, which could confound our understanding about the frequency at which they actually occur in nature. Here, we explore genome-wide traces of divergence between two populations of the lacertid lizard Psammodromus algirus separated by a 600 m elevational gradient. These populations seem to be differentially adapted to their environments despite showing low levels of genetic differentiation (according to previously studies of mtDNA and microsatellite data). We performed a search for outliers (i.e., loci subject to selection) trying to identify specific loci with FST statistics significantly higher than those expected on the basis of overall, genome-wide estimates of genetic divergence. We find that local phenotypic adaptation (in terms of a wide diversity of characters) was not accompanied by genome-wide differentiation, even when we maximized the chances of unveiling such differentiation at particular loci with FST-based outlier detection tests. Instead, our analyses confirmed the lack of genome-wide differentiation on the basis of more than 70,000 SNPs, which is concordant with a scenario of local adaptation without isolation by environment. Our results add evidence to previous studies in which local adaptation does not lead to any kind of isolation (or early stages of ecological speciation), but maintains phenotypic divergence despite the lack of a differentiated genomic background.

Comparative Metagenomics of Palearctic and Neotropical Avian Cloacal Viromes Reveal Geographic Bias in Virus Discovery.

Our understanding about viruses carried by wild animals is still scarce. The viral diversity of wildlife may be best described with discovery-driven approaches to the study of viral diversity that broaden research efforts towards non-canonical hosts and remote geographic regions. Birds have been key organisms in the transmission of viruses causing important diseases, and wild birds are threatened by viral spillovers associated with human activities. However, our knowledge of the avian virome may be biased towards poultry and highly pathogenic diseases. We describe and compare the fecal virome of two passerine-dominated bird assemblages sampled in a remote Neotropical rainforest in French Guiana (Nouragues Natural Reserve) and a Mediterranean forest in central Spain (La Herrería). We used metagenomic data to quantify the degree of functional and genetic novelty of viruses recovered by examining if the similarity of the contigs we obtained to reference sequences differed between both locations. In general, contigs from Nouragues were significantly less similar to viruses in databases than contigs from La Herrería using Blastn but not for Blastx, suggesting that pristine regions harbor a yet unknown viral diversity with genetically more singular viruses than more studied areas. Additionally, we describe putative novel viruses of the families Picornaviridae, Reoviridae and Hepeviridae. These results highlight the importance of wild animals and remote regions as sources of novel viruses that substantially broaden the current knowledge of the global diversity of viruses.

The combined use of raw and phylogenetically independent methods of outlier detection uncovers genome-wide dynamics of local adaptation in a lizard.

Local adaptation is a dynamic process by which different allele combinations are selected in different populations at different times, and whose genetic signature can be inferred by genome-wide outlier analyses. We combined gene flow estimates with two methods of outlier detection, one of them independent of population coancestry (CIOA) and the other one not (ROA), to identify genetic variants favored when ecology promotes phenotypic convergence. We analyzed genotyping-by-sequencing data from five populations of a lizard distributed over an environmentally heterogeneous range that has been changing since the split of eastern and western lineages ca. 3 mya. Overall, western lizards inhabit forest habitat and are unstriped, whereas eastern ones inhabit shrublands and are striped. However, one population (Lerma) has unstriped phenotype despite its eastern ancestry. The analysis of 73,291 SNPs confirmed the east-west division and identified nonoverlapping sets of outliers (12 identified by ROA and 9 by CIOA). ROA revealed ancestral adaptive variation in the uncovered outliers that were subject to divergent selection and differently fixed for eastern and western populations at the extremes of the environmental gradient. Interestingly, such variation was maintained in Lerma, where we found high levels of heterozygosity for ROA outliers, whereas CIOA uncovered innovative variants that were selected only there. Overall, it seems that both the maintenance of ancestral variation and asymmetric migration have counterbalanced adaptive lineage splitting in our model species. This scenario, which is likely promoted by a changing and heterogeneous environment, could hamper ecological speciation of locally adapted populations despite strong genetic structure between lineages.