Influence of geographic isolation and the environment on gene flow among phenotypically diverse lizards

Lizards in the genus Anolis comprise hundreds of species that display a wide range of phenotypic variation closely related to their environment. One example is the Guadeloupean anole (Anolis marmoratus ssp.) that display extreme phenotypic variation, primarily in adult male color and pattern, with twelve described subspecies on the archipelago. Here we examine the relationship between phenotypic and genetic divergence among five subspecies on the two main islands and test the role of geographic isolation and the environment in reducing gene flow. We also examined two offshore island populations to assess the impact of complete geographic isolation on gene flow. We analyzed color phenotypes by measuring spectral reflectance and genomic diversity using SNPs. Genetic divergence was correlated with dorsolateral head and body color phenotypes, and slope and geographic distance were nearly equivalent at explaining this divergence. There was minimal genome-wide divergence at neutral loci among phenotypically disparate subspecies on the two main islands and their differentiation is consistent with a model of divergence with gene flow. Our spatial visualization of gene flow showed an impact of environmental features consistent with a hypothesis of ecologically driven divergence. Nonetheless, subspecies on the two main islands remain interconnected by substantial gene flow and their phenotypic variation is likely maintained at selection-gene flow equilibrium by divergent selection at loci associated with their color phenotypes. Greater isolation, such as inhabiting a remote island, may be required for reducing gene flow. Our findings highlight the role of the environment, adaptation, and geographic isolation on gene flow.

Environment influences the genetic structure and genetic differentiation of Sassafras tzumu (Lauraceae)

BACKGROUND

Sassafras tzumu, an elegant deciduous arboreal species, belongs to the esteemed genus Sassafras within the distinguished family Lauraceae. With its immense commercial value, escalating market demands and unforeseen human activities within its natural habitat have emerged as new threats to S. tzumu in recent decades, so it is necessary to study its genetic diversity and influencing factors, to propose correlative conservation strategies.

RESULTS

By utilizing genotyping-by-sequence (GBS) technology, we acquired a comprehensive database of single nucleotide polymorphisms (SNPs) from a cohort of 106 individuals sourced from 13 diverse Sassafras tzumu natural populations, scattered across various Chinese mountainous regions. Through our meticulous analysis, we aimed to unravel the intricate genetic diversity and structure within these S. tzumu populations, while simultaneously investigating the various factors that potentially shape genetic distance. Our preliminary findings unveiled a moderate level of genetic differentiation (FST = 0.103, p < 0.01), accompanied by a reasonably high genetic diversity among the S. tzumu populations. Encouragingly, our principal component analysis painted a vivid picture of two distinct genetic and geographical regions across China, where gene flow appeared to be somewhat restricted. Furthermore, employing the sophisticated multiple matrix regression with randomization (MMRR) analysis method, we successfully ascertained that environmental distance exerted a more pronounced impact on genetic distance when compared to geographical distance (βE = 0.46, p < 0.01; βD = 0.16, p < 0.01). This intriguing discovery underscores the potential significance of environmental factors in shaping the genetic landscape of S. tzumu populations.

CONCLUSIONS

The genetic variance among populations of S. tzumu in our investigation exhibited a moderate degree of differentiation, alongside a heightened level of genetic diversity. The environmental distance of S. tzumu had a greater impact on its genetic diversity than geographical distance. It is of utmost significance to formulate and implement meticulous management and conservation strategies to safeguard the invaluable genetic resources of S. tzumu.

Distribution and diversification of Adelphobates, emblematic poison frogs from Brazilian Amazonia

Adelphobates contains three species, and the inaccurate identification of A. quinquevittatus and the scarcity of records of A. castaneoticus complicate inference of their distributions; the latter species occurs in sympatry with A. galactonotus. Our objective was to revise the distributions of Adelphobates by compiling data and modeling habitat suitability, as range limits may be shaped by landscape features and biotic interactions. We initially analyzed the existence of operational taxonomic units within the nominal species and subsequently inferred the observed and potential distributions, taking into account the possible independent lineages for the three species, and we also generated a molecular timetree to understand the chronology of interspecific diversification events. Adelphobates quinquevittatus was found to have a more easterly distribution than previously described, and specimens with phenotypic variation were found to occur in areas inconsistent with the modeling, and A. castaneoticus was concentrated in the Tapajós-Xingu interfluve, surrounded by A. galactonotus. Models indicated that the right bank of the Xingu River is suitable for both species, indeed, both were found there. Despite Adelphobates species having their distributions delimited by major Amazonian rivers, estimated divergence times predate the formation of the modern river network, suggesting that other mechanisms were involved in their diversification.

Relationship between genetic diversity and morpho-functional characteristics of flight-related traits in Triatoma garciabesi (Hemiptera: Reduviidae)

BACKGROUND

Triatoma garciabesi, a potential vector of the parasitic protozoan Trypanosoma cruzi, which is the causative agent of Chagas disease, is common in peridomestic and wild environments and found throughout northwestern and central Argentina, western Paraguay and the Bolivian Chaco. Genetic differentiation of a species across its range can help to understand dispersal patterns and connectivity between habitats. Dispersal by flight is considered to be the main active dispersal strategy used by triatomines. In particular, the morphological structure of the hemelytra is associated with their function. The aim of this study was to understand how genetic diversity is structured, how morphological variation of dispersal-related traits varies with genetic diversity and how the morphological characteristics of dispersal-related traits may explain the current distribution of genetic lineages in this species.

METHODS

Males from 24 populations of T. garciabesi across its distribution range were examined. The cytochrome c oxidase I gene (coI) was used for genetic diversity analyses. A geometric morphometric method based on landmarks was used for morpho-functional analysis of the hemelytra. Centroid size (CS) and shape of the forewing, and contour of both parts of the forewing, the head and the pronotum were characterised. Length and area of the forewing were measured to estimate the aspect ratio.

RESULTS

The morphometric and phylogenetic analysis identified two distinct lineages, namely the Eastern and Western lineages, which coincide with different ecological regions. The Eastern lineage is found exclusively in the eastern region of Argentina (Chaco and Formosa provinces), whereas the Western lineage is prevalent in the rest of the geographical range of the species. CS, shape and aspect ratio of the hemelytra differed between lineages. The stiff portion of the forewing was more developed in the Eastern lineage. The shape of both portions of the hemelytra were significantly different between lineages, and the shape of the head and pronotum differed between lineages.

CONCLUSIONS

The results provide preliminary insights into the evolution and diversification of T. garciabesi. Variation in the forewing, pronotum and head is congruent with genetic divergence. Consistent with genetic divergence, morphometry variation was clustered according to lineages, with congruent variation in the size and shape of the forewing, pronotum and head.

Genetic structure, UV-vision, wing coloration and size coincide with colour polymorphism in Fabriciana adippe butterflies

Colour polymorphisms have long served as model systems in evolutionary studies and continue to inform about processes involved in the origin and dynamics of biodiversity. Modern sequencing tools allow for evaluating whether phenotypic differences between morphs reflect genetic differentiation rather than developmental plasticity, and for investigating whether polymorphisms represent intermediate stages of diversification towards speciation. We investigated phenotypic and genetic differentiation between two colour morphs of the butterfly Fabriciana adippe using a combination of ddRAD-sequencing and comparisons of body size, colour patterns and optical properties of bright wing spots. The silvery-spotted adippe form had larger and darker wings and reflected UV light, while the yellow cleodoxa form displayed more green scales and reflected very little UV, showcasing that they constitute distinct and alternative integrated phenotypes. Genomic analyses revealed genetic structuring according to source population, and to colour morph, suggesting that the phenotypic differentiation reflects evolutionary modifications. We report 17 outlier loci associated with colour morph, including ultraviolet-sensitive visual pigment (UVRh1), which is associated with intraspecific communication and mate choice in butterflies. Together with the demonstration that the wings of the adippe (but essentially not the cleodoxa) morph reflect UV light, that UV reflectance is higher in females than males and that morphs differ in wing size, this suggests that these colour morphs might represent genetically integrated phenotypes, possibly adapted to different microhabitats. We propose that non-random mating might contribute to the differentiation and maintenance of the polymorphism.

Landscape Heterogeneity Explains the Genetic Differentiation of a Forest Bird across the Sino-Himalayan Mountains

Mountains are the world's most important centers of biodiversity. The Sino-Himalayan Mountains are global biodiversity hotspot due to their extremely high species richness and endemicity. Ample research investigated the impact of the Qinghai-Tibet Plateau uplift and Quaternary glaciations in driving species diversification in plants and animals across the Sino-Himalayan Mountains. However, little is known about the role of landscape heterogeneity and other environmental features in driving diversification in this region. We utilized whole genomes and phenotypic data in combination with landscape genetic approaches to investigate population structure, demography, and genetic diversity in a forest songbird species native to the Sino-Himalayan Mountains, the red-billed leiothrix (Leiothrix lutea). We identified 5 phylogeographic clades, including 1 in the East of China, 1 in Yunnan, and 3 in Tibet, roughly consistent with differences in song and plumage coloration but incongruent with traditional subspecies boundaries. Isolation-by-resistance model best explained population differentiation within L. lutea, with extensive secondary contact after allopatric isolation leading to admixture among clades. Ecological niche modeling indicated relative stability in the extent of suitable distribution areas of the species across Quaternary glacial cycles. Our results underscore the importance of mountains in the diversification of this species, given that most of the distinct genetic clades are concentrated in a relatively small area in the Sino-Himalayan Mountain region, while a single shallow clade populates vast lower-lying areas to the east. This study highlights the crucial role of landscape heterogeneity in promoting differentiation and provides a deep genomic perspective on the mechanisms through which diversity hotspots form.

Anti-predator defences are linked with high levels of genetic differentiation in frogs

Predator-prey interactions have been suggested as drivers of diversity in different lineages, and the presence of anti-predator defences in some clades is linked to higher rates of diversification. Warning signals are some of the most widespread defences in the animal world, and there is evidence of higher diversification rates in aposematic lineages. The mechanisms behind such species richness, however, are still unclear. Here, we test whether lineages that use aposematism as anti-predator defence exhibit higher levels of genetic differentiation between populations, leading to increased opportunities for divergence. We collated from the literature more than 3000 pairwise genetic differentiation values across more than 700 populations from over 60 amphibian species. We find evidence that over short geographical distances, populations of species of aposematic lineages exhibit greater genetic divergence relative to species that are not aposematic. Our results support a scenario where the use of warning signals could restrict gene flow, and suggest that anti-predator defences could impact divergence between populations and potentially have effects at a macro-evolutionary scale.

Himalayan mountains imposing a barrier on gene flow of wheat yellow rust pathogen in the bordering regions of Pakistan and China

The wheat yellow rust pathogen has been shown to be diverse and potentially originated in the Himalayan region. Although Himalayan populations of Pakistan, Nepal and Bhutan have been previously compared, little is known about the relative divergence and diversity in Puccinia striiformis populations in the bordering regions of Pakistan and China. To assess the relative diversity and divergence in these regions of Pakistan (Gilgit-Baltistan, Hazara and Azad Jammu Kashmir) and China (Xinjiang, Qinghai, Tibet, Sichuan, Guizhou and Yunnan), a total of 1245 samples were genotyped using 17 microsatellite SSR markers. A clear divergence was observed between the bordering regions of Pakistan and China (F_ST = 0.28) without any resampling of genetic groups and multilocus genotypes across two sides of the Himalayan mountains. The closest subpopulations across the two countries were Xinjiang and Gilgit-Baltistan (Nei's distance = 0.147), which were close geographically. A very high diversity and recombinant population structure was observed in both populations, though slightly higher in China (Genotypic diversity = 0.970; r¯d = 0.000) than in Pakistan (Genotypic diversity = 0.902; r¯d = 0.065). The distribution of genetic groups and resampling of MLGs revealed more gene flow across Yunnan, Guizhou and Sichuan regions in China, while between Hazara and Azad-Jammu Kashmir in Pakistan. The lack of gene flow between Pakistan and China populations is due to geographical barriers and a large patch of land without wheat. The information on the relative diversity and divergence in different geographical zones of the pathogen center of diversity and neighboring region should be considered in resistant wheat deployment while considering the invasion potential of the pathogen at regional and global contexts.

The impact of Pleistocene glaciations and environmental gradients on the genetic structure of Embothrium coccineum

The South American temperate forests were subjected to drastic topographic and climatic changes during the Pliocene-Pleistocene as a consequence of the Andean orogeny and glacial cycles. Such changes are common drivers of genetic structure and adaptation. Embothrium coccineum (Proteaceae) is an emblematic tree of the South American temperate forest (around 20°S of latitude) that has strongly been affected by topographic and climatic events. Previous studies have shown a marked genetic structure in this species, and distinct ecotypes have been described. Yet, little is known about their adaptive genetic responses. The main goal of this study was to investigate the effects of historical and contemporary landscape features affecting the genetic diversity and connectivity of E. coccineum throughout its current natural distribution. Using over 2000 single nucleotide polymorphisms (SNPs), we identified two genetic groups (a Northern and a Central-Southern group) that diverged around 2.8 million years ago. The level of genetic structure was higher among populations within the Northern genetic group than within the Central-Southern group. We propose that these differences in genetic structure may be due to differences in the assemblages of pollinators and in the evolutionary histories of the two genetic groups. Moreover, the data displayed a strong pattern of isolation by the environment in E. coccineum, suggesting that selection could have led to adaptive divergence among localities. We propose that in the Chilean temperate forest, the patterns of genetic variation in E. coccineum reflect both a Quaternary phylogenetic imprint and signatures of selection as a consequence of a strong environmental gradient.

Host association, environment, and geography underlie genomic differentiation in a major forest pest

Diverse geographic, environmental, and ecological factors affect gene flow and adaptive genomic variation within species. With recent advances in landscape ecological modelling and high-throughput DNA sequencing, it is now possible to effectively quantify and partition their relative contributions. Here, we use landscape genomics to identify determinants of genomic differentiation in the forest tent caterpillar, Malacosoma disstria, a widespread and irruptive pest of numerous deciduous tree species in North America. We collected larvae from multiple populations across Eastern Canada, where the species experiences a diversity of environmental gradients and feeds on a number of different host tree species, including trembling aspen (Populus tremuloides), sugar maple (Acer saccharum), red oak (Quercus rubra), and white birch (Betula papyrifera). Using a combination of reciprocal causal modelling (RCM) and distance-based redundancy analyses (dbRDA), we show that differentiation of thousands of genome-wide single nucleotide polymorphisms (SNPs) among individuals is best explained by a combination of isolation by distance, isolation by environment (spatial variation in summer temperatures and length of the growing season), and differences in host association. Configuration of suitable habitat inferred from ecological niche models was not significantly related to genomic differentiation, suggesting that M. disstria dispersal is agnostic with respect to habitat quality. Although population structure was not discretely related to host association, our modelling framework provides the first molecular evidence of host-associated differentiation in M. disstria, congruent with previous documentation of reduced growth and survival of larvae moved between natal host species. We conclude that ecologically mediated selection is contributing to variation within M. disstria, and that divergent adaptation related to both environmental conditions and host association should be considered in ongoing research and management of this important forest pest.

Climate change will disproportionally affect the most genetically diverse lineages of a widespread African tree species

Global climate change is proceeding at an alarming rate with major ecological and genetic consequences for biodiversity, particularly in drylands. The response of species to climate change may differ between intraspecific genetic groups, with major implications for conservation. We used molecular data from 10 nuclear and two chloroplast genomes to identify phylogeographic groups within 746 individuals from 29 populations of Senegalia senegal, a savannah tree species in sub-Saharan Africa. Three phylogroups are identified corresponding to Sudano-Sahelian, Zambezian and Southern African biogeographic regions in West, East and Southern Africa. Genetic diversity was highest in Southern and Zambesian and lowest in the Sudano-Sahelian phylogroups. Using species distribution modeling, we infer highly divergent future distributions of the phylogroups under three climate change scenarios. Climate change will lead to severe reductions of distribution area of the genetically diverse Zambezian (- 41-- 54%) and Southern (- 63-- 82%) phylogroups, but to an increase for the genetically depauperate Sudano-Sahelian (+ 7- + 26%) phylogroups. This study improves our understanding of the impact of climate change on the future distribution of this species. This knowledge is particularly useful for biodiversity management as the conservation of genetic resources needs to be considered in complementary strategies of in-situ conservation and assisted migration.

Fine-scale genome-wide signature of Pleistocene glaciation in Thitarodes moths (Lepidoptera: Hepialidae), host of Ophiocordyceps fungus in the Hengduan Mountains

The Hengduan Mountains region is a biodiversity hotspot known for its topologically complex, deep valleys and high mountains. While landscape and glacial refugia have been evoked to explain patterns of inter-species divergence, the accumulation of intra-species (i.e. population level) genetic divergence across the mountain-valley landscape in this region has received less attention. We used genome-wide restriction site-associated DNA sequencing (RADseq) to reveal signatures of Pleistocene glaciation in populations of Thitarodes shambalaensis (Lepidoptera: Hepialidae), the host moth of parasitic Ophiocordyceps sinensis (Hypocreales: Ophiocordycipitaceae) or "caterpillar fungus" endemic to the glacier of eastern Mt. Gongga. We used moraine history along the glacier valleys to model the distribution and environmental barriers to gene flow across populations of T. shambalaensis. We found that moth populations separated by less than 10 km exhibited valley-based population genetic clustering and isolation-by-distance (IBD), while gene flow among populations was best explained by models using information about their distributions at the local last glacial maximum (LGM_L , 58 kya), not their contemporary distribution. Maximum likelihood lineage history among populations, and among subpopulations as little as 500 meters apart, recapitulated glaciation history across the landscape. We also found signals of isolated population expansion following the retreat of LGM_L glaciers. These results reveal the fine-scale, long-term historical influence of landscape and glaciation on the genetic structuring of populations of an endangered and economically important insect species. Similar mechanisms, given enough time and continued isolation, could explain the contribution of glacier refugia to the generation of species diversity among the Hengduan Mountains.

Genetic isolation by distance underlies colour pattern divergence in red-eyed treefrogs (Agalychnis callidryas)

Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of colour-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg colour patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple colour-pattern morphs. To explore possible causes of this variation, we generated a single nucleotide polymorphism data set to analyse population genetic structure, measure genetic diversity and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond to individual measures of colour pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that colour-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance indicates that colour-pattern variation was either retained through other processes such as ancestral colour polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on colour pattern at a small geographic scale.

Environmental differences explain subtle yet detectable genetic structure in a widespread pollinator

BACKGROUND

The environment is a strong driver of genetic structure in many natural populations, yet often neglected in population genetic studies. This may be a particular problem in vagile species, where subtle structure cannot be explained by limitations to dispersal. Consequently, these species might falsely be considered quasi-panmictic and hence potentially mismanaged. A species this might apply to, is the buff-tailed bumble bee (Bombus terrestris), an economically important and widespread pollinator, which is considered to be quasi-panmictic at mainland continental scales. Here we aimed to (i) quantify genetic structure in 21+ populations of the buff-tailed bumble bee, sampled throughout two Eastern European countries, and (ii) analyse the degree to which structure is explained by environmental differences, habitat permeability and geographic distance. Using 12 microsatellite loci, we characterised populations of this species with Fst analyses, complemented by discriminant analysis of principal components and Bayesian clustering approaches. We then applied generalized dissimilarity modelling to simultaneously assess the informativeness of geographic distance, habitat permeability and environmental differences among populations in explaining divergence.

RESULTS

Genetic structure of the buff-tailed bumble bee quantified by means of Fst was subtle and not detected by Bayesian clustering. Discriminant analysis of principal components suggested insignificant but still noticeable structure that slightly exceeded estimates obtained through Fst analyses. As expected, geographic distance and habitat permeability were not informative in explaining the spatial pattern of genetic divergence. Yet, environmental variables related to temperature, vegetation and topography were highly informative, explaining between 33 and 39% of the genetic variation observed.

CONCLUSIONS

In contrast to previous studies reporting quasi-panmixia in continental populations of this species, we demonstrated the presence of subtle population structure related to environmental heterogeneity. Environmental data proved to be highly useful in unravelling the drivers of genetic structure in this vagile and opportunistic species. We highlight the potential of including these data to obtain a better understanding of population structure and the processes driving it in species considered to be quasi-panmictic.

The impact of life stage and pigment source on the evolution of novel warning signal traits

Our understanding of how novel warning color traits evolve in natural populations is largely based on studies of reproductive stages and organisms with endogenously produced pigmentation. In these systems, genetic drift is often required for novel alleles to overcome strong purifying selection stemming from frequency‐dependent predation and positive assortative mating. Here, we integrate data from field surveys, predation experiments, population genomics, and phenotypic correlations to explain the origin and maintenance of geographic variation in a diet‐based larval pigmentation trait in the redheaded pine sawfly (Neodiprion lecontei), a pine‐feeding hymenopteran. Although our experiments confirm that N. lecontei larvae are indeed aposematic—and therefore likely to experience frequency‐dependent predation—our genomic data do not support a historical demographic scenario that would have facilitated the spread of an initially deleterious allele via drift. Additionally, significantly elevated differentiation at a known color locus suggests that geographic variation in larval color is currently maintained by selection. Together, these data suggest that the novel white morph likely spread via selection. However, white body color does not enhance aposematic displays, nor is it correlated with enhanced chemical defense or immune function. Instead, the derived white‐bodied morph is disproportionately abundant on a pine species with a reduced carotenoid content relative to other pine hosts, suggesting that bottom‐up selection via host plants may have driven divergence among populations. Overall, our results suggest that life stage and pigment source can have a substantial impact on the evolution of novel warning signals, highlighting the need to investigate diverse aposematic taxa to develop a comprehensive understanding of color variation in nature.

Otolith Microchemistry and Demographic History Provide New Insight into the Migratory Behavior and Heterogeneous Genetic Divergence of Coilia grayii in the Pearl River

Coilia grayii is the anadromous form of anchovy that is distributed in the East and South China Seas. It is a common fish species in the estuarine area of the Pearl River. Nevertheless, freshwater populations appear upstream in the Pearl River, but the migratory pathway has been mostly impeded by dam construction. Behavioral differences and constrained habitat within tributaries are suspected of promoting genetic divergence in these populations. In this study, we investigated the migratory behavior and genetic divergence of six populations of C. grayii fragmented by dams based on the otolith strontium/calcium (Sr/Ca) ratio, mitochondrial DNA, and microsatellite genotyping. All populations were in freshwater with low Sr/Ca ratios, except the estuarine population (Humen population) hatched in brackish water. Reduced nucleotide diversity corresponding to distance was observed. Populations from distant hydrological regions exhibited a decline in genetic diversity and a significant difference with the remaining populations after fitting the isolation by distance model. Pairwise fixation indices confirmed these results and moderate and significant differentiation was found between Hengxian site and downstream sites. Furthermore, STRUCTURE analyses revealed that all separated populations exhibited an admixed phylogenetic pattern except for individuals from the Hengxian locality. The upstream sites showed significantly increased resistance to gene flow from the estuarine population because of isolation by the dam. The results of the neutrality test and Bayesian skyline plots demonstrated complex demography—individuals’ experienced historical expansion and partial upper-dam populations had recently undergone a colonization, forming a new genetic structure. Accordingly, this study demonstrates differences in the migration pattern and genetic differentiation of C. grayii as a consequence of demographic history and current processes (habitat fragmentation and colonization).

Ecological divergence and synchronous Pleistocene diversification in the widespread South American butter frog complex

Phylogeographic studies primarily focus on the major role of landscape topography in driving lineage diversification. However, populational phylogeographic breaks may also occur as a result of either niche conservatism or divergence, in the absence of geographic barriers to gene flow. Furthermore, these two factors are not mutually exclusive and can act in concert, making it challenging to evaluate their relative importance on explaining genetic variation in nature. Herein, we use sequences of two mitochondrial and four nuclear genes to investigate the timing and diversification patterns of species pertaining to the Leptodactylus latrans complex, which harbors four morphologically cryptic species with broad distributions across environmental gradients in eastern South America. The origin of this species complex dates back to the late Miocene (ca. 5.5 Mya), but most diversification events occurred synchronically during the late Pleistocene likely as the result of ecological divergence driven by Quaternary climatic oscillations. Further, significant patterns of environmental niche divergences among species in the L. latrans complex imply that ecological isolation is the primary mode of genetic diversification, mostly because phylogenetic breaks are associated with environmental transitions rather than topographic barriers at both species and populational scale. We provided new insights about diversification patterns and processes within a species complex of broadly and continuously distributed group of frogs along South America.

Contrasting environmental drivers of genetic and phenotypic divergence in an Andean poison frog (Epipedobates anthonyi)

Phenotypic and genetic divergence are shaped by the homogenizing effects of gene flow and the differentiating processes of genetic drift and local adaptation. Herein, we examined the mechanisms that underlie phenotypic (size and color) and genetic divergence in 35 populations (535 individuals) of the poison frog Epipedobates anthonyi along four elevational gradients (0-1800 m asl) in the Ecuadorian Andes. We found phenotypic divergence in size and color despite relatively low genetic divergence at neutral microsatellite loci. Genetic and phenotypic divergence were both explained by landscape resistance between sites (isolation-by-resistance, IBR), likely due to a cold and dry mountain ridge between the northern and southern elevational transects that limits dispersal and separates two color morphs. Moreover, environmental differences among sites also explained genetic and phenotypic divergence, suggesting isolation-by-environment (IBE). When northern and southern transects were analyzed separately, genetic divergence was predicted either by distance (isolation-by-distance, IBD; northern) or environmental resistance between sites (IBR; southern). In contrast, phenotypic divergence was primarily explained by environmental differences among sites, supporting the IBE hypothesis. These results indicate that although distance and geographic barriers are important drivers of population divergence, environmental variation has a two-fold effect on population divergence. On the one hand, landscape resistance between sites reduces gene flow (IBR), while on the other hand, environmental differences among sites exert divergent selective pressures on phenotypic traits (IBE). Our work highlights the importance of studying both genetic and phenotypic divergence to better understand the processes of population divergence and speciation along ecological gradients.

Genetic structure and ecological niche segregation of Indian gray mongoose (Urva edwardsii) in Iran

Combining genetic data with ecological niche models is an effective approach for exploring climatic and nonclimatic environmental variables affecting spatial patterns of intraspecific genetic variation. Here, we adopted this combined approach to evaluate genetic structure and ecological niche of the Indian gray mongoose (Urva edwardsii) in Iran, as the most western part of the species range. Using mtDNA, we confirmed the presence of two highly differentiated clades. Then, we incorporated ensemble of small models (ESMs) using climatic and nonclimatic variables with genetic data to assess whether genetic differentiation among clades was coupled with their ecological niche. Climate niche divergence was also examined based on a principal component analysis on climatic factors only. The relative habitat suitability values predicted by the ESMs for both clades revealed their niche separation. Between-clade climate only niche comparison revealed that climate space occupied by clades is similar to some extent, but the niches that they utilize differ between the distribution ranges of clades. We found that in the absence of evidence for recent genetic exchanges, distribution models suggest the species occurs in different niches and that there are apparent areas of disconnection across the species range. The estimated divergence time between the two Iranian clades (4.9 Mya) coincides with the uplifting of the Zagros Mountains during the Early Pliocene. The Zagros mountain-building event seems to have prevented the distribution of U. edwardsii populations between the western and eastern parts of the mountains as a result of vicariance events. Our findings indicated that the two U. edwardsii genetic clades in Iran can be considered as two conservation units and can be utilized to develop habitat-specific and climate change-integrated management strategies.

A strategic sampling design revealed the local genetic structure of cold-water fluvial sculpin: a focus on groundwater-dependent water temperature heterogeneity

A key piece of information for ecosystem management is the relationship between the environment and population genetic structure. However, it is difficult to clearly quantify the effects of environmental factors on genetic differentiation because of spatial autocorrelation and analytical problems. In this study, we focused on stream ecosystems and the environmental heterogeneity caused by groundwater and constructed a sampling design in which geographic distance and environmental differences are not correlated. Using multiplexed ISSR genotyping by sequencing (MIG-seq) method, a fine-scale population genetics study was conducted in fluvial sculpin Cottus nozawae, for which summer water temperature is the determinant factor in distribution and survival. There was a clear genetic structure in the watershed. Although a significant isolation-by-distance pattern was detected in the watershed, there was no association between genetic differentiation and water temperature. Instead, asymmetric gene flow from relatively low-temperature streams to high-temperature streams was detected, indicating the importance of low-temperature streams and continuous habitats. The groundwater-focused sampling strategy yielded insightful results for conservation.

Weak spatial-genetic structure in a native invasive, the southern pine beetle (Dendroctonus frontalis), across the eastern United States

The southern pine beetle, Dendroctonus frontalis, is a native pest of pine trees that has recently expanded its range into the northeastern United States. Understanding its colonization, dispersal, and connectivity will be critical for mitigating negative economic and ecological impacts in the newly invaded areas. Characterization of spatial-genetic structure can contribute to this; however, previous studies have reached different conclusions about regional population genetic structure, with one study reporting a weak east-west pattern, and the most recent reporting an absence of structure. Here we systematically assessed several explanations for the absence of spatial-genetic structure. To do this, we developed nine new microsatellite markers and combined them with an existing 24-locus data matrix for the same individuals. We then reanalyzed this full dataset alongside datasets in which certain loci were omitted with the goal of creating more favorable signal to noise ratios. We also partitioned the data based on the sex of D. frontalis individuals, and then employed a broad suite of genotypic clustering and isolation-by-distance (IBD) analyses. We found that neither inadequate information content in the molecular marker set, nor unfavorable signal-to-noise ratio, nor insensitivity of the analytical approaches could explain the absence of structure. Regardless of dataset composition, there was little evidence for clusters (i.e., distinct geo-genetic groups) or clines (i.e., gradients of increasing allele frequency differences over larger geographic distances), with one exception: significant IBD was repeatedly detected using an individual-based measure of relatedness whenever datasets included males (but not for female-only datasets). This is strongly indicative of broad-scale female-biased dispersal, which has not previously been reported for D. frontalis, in part owing to logistical limitations of direct approaches (e.g., capture-mark-recapture). Weak spatial-genetic structure suggests long-distance connectivity and that gene flow is high, but additional research is needed to understand range expansion and outbreak dynamics in this species using alternate approaches.

Phenotypic and genomic diversification with isolation by environment along elevational gradients in a neotropical treefrog

Understanding how geographic and environmental heterogeneity drive local patterns of genetic variation is a major goal of ecological genomics and a key question in evolutionary biology. The tropical Andes and inter-Andean valleys are shaped by markedly heterogeneous landscapes, where species experience strong selective processes. We examined genome-wide SNP data together with behavioural and ecological traits (mating calls and body size) known to contribute to genetic isolation in anurans in the banana tree-dwelling frog, Boana platanera, distributed across an environmental gradient in Central Colombia (northern South America). Here, we analysed the relationships between environmentally (temperature and precipitation) associated genetic and phenotypic differentiation and the potential drivers of isolation by environment along an elevation gradient. We identified candidate SNPs associated with temperature and body size, which follow a clinal pattern of genome-wide differentiation tightly coupled with phenotypic variation: as elevation increases, B. platanera exhibits larger body size and longer call duration with more pulses but lower pulse rate and frequency. Thus, the environmental landscape has rendered a scenario where isolation by environment and candidate loci show concordance with phenotypic divergence in this tropical frog along an elevation gradient in the Colombian Andes. Our study sets the basis for evaluating the role of temperature in the genetic structure and local adaptation in tropical treefrogs and its putative effect on life cycle (embryos, tadpoles, adults) along elevation gradients.

Multiscale patterns of isolation by ecology and fine-scale population structure in Texas bobcats

Patterns of spatial genetic variation can be generated by a variety of ecological processes, including individual preferences based on habitat. These ecological processes act at multiple spatial and temporal scales, generating scale-dependent effects on gene flow. In this study, we focused on bobcats (Lynx rufus), a highly mobile, generalist felid that exhibits ecological and behavioral plasticity, high abundance, and broad connectivity across much of their range. However, bobcats also show genetic differentiation along habitat breaks, a pattern typically observed in cases of isolation-by-ecology (IBE). The IBE observed in bobcats is hypothesized to occur due to habitat-biased dispersal, but it is unknown if this occurs at other habitat breaks across their range or at what spatial scale IBE becomes most apparent. Thus, we used a multiscale approach to examine isolation by ecology (IBE) patterns in bobcats (Lynx rufus) at both fine and broad spatial scales in western Texas. We genotyped 102 individuals at nine microsatellite loci and used partial redundancy analysis (pRDA) to test if a suite of landscape variables influenced genetic variation in bobcats. Bobcats exhibited a latitudinal cline in population structure with a spatial signature of male-biased dispersal, and no clear barriers to gene flow. Our pRDA tests revealed high genetic similarity in similar habitats, and results differed by spatial scale. At the fine spatial scale, herbaceous rangeland was an important influence on gene flow whereas mixed rangeland and agriculture were significant at the broad spatial scale. Taken together, our results suggests that complex interactions between spatial-use behavior and landscape heterogeneity can create non-random gene flow in highly mobile species like bobcats. Furthermore, our results add to the growing body of data highlighting the importance of multiscale study designs when assessing spatial genetic structure.

Ecological adaptation drives wood frog population divergence in life history traits

Phenotypic variation among populations is thought to be generated from spatial heterogeneity in environments that exert selection pressures that overcome the effects of gene flow and genetic drift. Here, we tested for evidence of isolation by distance or by ecology (i.e., ecological adaptation) to generate variation in early life history traits and phenotypic plasticity among 13 wood frog populations spanning 1200 km and 7° latitude. We conducted a common garden experiment and related trait variation to an ecological gradient derived from an ecological niche model (ENM) validated to account for population density variation. Shorter larval periods, smaller body weight, and relative leg lengths were exhibited by populations with colder mean annual temperatures, greater precipitation, and less seasonality in precipitation and higher population density (high-suitability ENM values). After accounting for neutral genetic variation, the Q_ST-F_ST analysis supported ecological selection as the key process generating population divergence. Further, the relationship between ecology and traits was dependent upon larval density. Specifically, high-suitability/high-density populations in the northern part of the range were better at coping with greater conspecific competition, evidenced by greater postmetamorphic survival and no difference in body weight when reared under stressful conditions of high larval density. Our results support that both climate and competition selection pressures drive clinal variation in larval and metamorphic traits in this species. Range-wide studies like this one are essential for accurate predictions of population's responses to ongoing ecological change.

Genetic Distinctiveness but Low Diversity Characterizes Rear-Edge Thuja standishii (Gordon) Carr. (Cupressaceae) Populations in Southwest Japan

Rear-edge populations are of significant scientific interest because they can contain allelic variation not found in core-range populations. However, such populations can differ in their level of genetic diversity and divergence reflecting variation in life-history traits, demographic histories and human impacts. Using 13 EST-microsatellites, we investigated the genetic diversity and differentiation of rear-edge populations of the Japanese endemic conifer Thuja standishii (Gordon) Carr. in southwest Japan from the core-range in northeast Japan. Range-wide genetic differentiation was moderate (Fst = 0.087), with northeast populations weakly differentiated (Fst = 0.047), but harboring high genetic diversity (average population-level Ar = 4.76 and Ho = 0.59). In contrast, rear-edge populations were genetically diverged (Fst = 0.168), but contained few unique alleles with lower genetic diversity (Ar = 3.73, Ho = 0.49). The divergence between rear-edge populations exceeding levels observed in the core-range and results from ABC analysis and species distribution modelling suggest that these populations are most likely relicts of the Last Glacial Maximum. However, despite long term persistence, low effective population size, low migration between populations and genetic drift have worked to promote the genetic differentiation of southwest Japan populations of T. standishii without the accumulation of unique alleles.

Using environmental and geographic data to optimize ex situ collections and preserve evolutionary potential

Maintenance of biodiversity through seed banks and botanical gardens, where the wealth of species' genetic variation may be preserved ex situ, is a major goal of conservation. However, challenges can persist in optimizing ex situ collections if trade-offs exist among cost, effort, and conserving species evolutionary potential, particularly when genetic data are not available. We evaluated the genetic consequences of population preservation informed by geographic (isolation by distance [IBD]) and environmental (isolation by environment [IBE]) distance for ex situ collections for which population provenance is available. We used 19 genetic and genomic data sets from 15 plant species to assess the proportion of population genetic differentiation explained by geographic and environmental factors and to simulate ex situ collections prioritizing source populations based on pairwise geographic distance, environmental distance, or both. Specifically, we tested the impact prioritizing sampling based on these distances may have on the capture of neutral, functional, or putatively adaptive genetic diversity and differentiation. Individually, IBD and IBE explained limited population genetic differences across all 3 genetic marker classes (IBD, 10-16%; IBE, 1-5.5%). Together, they explained a substantial proportion of population genetic differences for functional (45%) and adaptive (71%) variation. Simulated ex situ collections revealed that inclusion of IBD, IBE, or both increased allelic diversity and genetic differentiation captured among populations, particularly for loci that may be important for adaptation. Thus, prioritizing population collections based on environmental and geographic distance data can optimize genetic variation captured ex situ. For the vast majority of plant species for which there is no genetic information, these data are invaluable to conservation because they can guide preservation of genetic variation needed to maintain evolutionary potential within collections.

Type and intensity of surrounding human land use, not local environment, shape genetic structure of a native grassland plant

Landscape heterogeneity can shape genetic structure and functional connectivity of populations. When this heterogeneity imposes variable costs of moving across the landscape, populations can be structured according to a pattern of "isolation by resistance" (IBR). At the same time, divergent local environmental filters can limit gene flow, creating an alternative pattern of "isolation by environment" (IBE). Here, we evaluate IBR and IBE in the insect-pollinated, biennial plant Sabatia angularis (L.) Pursh (Gentianaceae) across serpentine grasslands in the fragmented landscape of SE Pennsylvania, USA using ~4500 neutral SNP loci. Specifically, we test the extent to which radical alteration of the landscape matrix by humans has fundamentally altered the cost of movement, imprinting a pattern of IBR dictated by land use type and intensity, and the potential for IBE in relation to a gradient of heavy metal concentrations found in serpentine soil. We reveal a strong signal of IBR and a weak signal of IBE across sites, indicating the greater importance of the landscape matrix in shaping genetic structure of S. angularis populations in the study region. Based on Circuitscape and least cost path approaches, we find that both low- and high-intensity urbanization resist gene flow by orders of magnitude greater than "natural" habitats, although resistance to low-intensity urbanization weakens at larger spatial scales. While cropland presents a substantially lower barrier than urban development, cumulative human land use surrounding populations predicts within-population genetic diversity and inbreeding in S. angularis. Our results emphasize the role of forest buffers and corridors in facilitating gene flow between serpentine grassland patches and averting local extinction of plant populations.

Freshwater zooplankton metapopulations and metacommunities respond differently to environmental and spatial variation

Theory predicts that population genetic structure and metacommunity structure are linked by the common processes of drift and migration, but how population genetic structure and metacommunity structure are related in nature is still unknown. Deeper understanding of the processes influencing both genetic and community diversity is vital for better predicting how environmental change will impact biodiversity patterns. We examined how crustacean zooplankton and rotifer species' metapopulation genetic structure and metacommunities respond to environmental and spatial variation both within and across four regions of boreal Canada. Metapopulation and metacommunity variation partitioning results were compared within and across the four regions. Metapopulations and metacommunities responded differently to environmental variation and spatial structure both within and across regions, as metapopulations were influenced by different environmental variables compared to metacommunities. At larger spatial scales both metapopulations and metacommunities exhibited greater spatial and environmental structuring, again responding to a different subset of environmental variables. Our findings suggest that even though both genetic and species diversity are linked by the same processes, regional variation in environmental characteristics and spatial structure influence resulting biodiversity patterns differently. To date, no other empirical research has explored relationships between entire metapopulation and metacommunity assemblages at large regional spatial scales.

Niche conservatism promotes speciation in cycads: the case of Dioon merolae (Zamiaceae) in Mexico

Niche conservatism is the tendency of lineages to retain the same niche as their ancestors. It constrains biological groups and prevents ecological divergence. However, theory predicts that niche conservatism can hinder gene flow, strengthen drift and increase local adaptation: does it mean that it also can facilitate speciation? Why does this happen? We aim to answer these questions. We examined the variation of chloroplast DNA, genome-wide single nucleotide polymorphisms, morphological traits and environmental variables across the Dioon merolae cycad populations. We tested geographical structure, scenarios of demographic history, and niche conservatism between population groups. Lineage divergence is associated with the presence of a geographical barrier consisting of unsuitable habitats for cycads. There is a clear genetic and morphological distinction between the geographical groups, suggesting allopatric divergence. However, even in contrasting available environmental conditions, groups retain their ancestral niche, supporting niche conservatism. Niche conservatism is a process that can promote speciation. In D. merolae, lineage divergence occurred because unsuitable habitats represented a barrier against gene flow, incurring populations to experience isolated demographic histories and disparate environmental conditions. This study explains why cycads, despite their ancient lineage origin and biological stasis, have been able to diversify into modern ecosystems worldwide.

Gene flow and climate-associated genetic variation in a vagile habitat specialist

Previous work in landscape genetics suggests that geographic isolation is of greater importance to genetic divergence than variation in environmental conditions. This is intuitive when configurations of suitable habitat are a dominant factor limiting dispersal and gene flow, but has not been thoroughly examined for habitat specialists with strong dispersal capability. Here, we evaluate the effects of geographic and environmental isolation on genetic divergence for a vagile invertebrate with high habitat specificity and a discrete dispersal life stage: Dod's Old World swallowtail butterfly, Papilio machaon dodi. In Canada, P. m. dodi are generally restricted to eroding habitat along major river valleys where their larval host plant occurs. A series of causal and linear mixed effects models indicate that divergence of genome-wide single nucleotide polymorphisms is best explained by a combination of environmental isolation (variation in summer temperatures) and geographic isolation (Euclidean distance). Interestingly, least-cost path and circuit distances through a resistance surface parameterized as the inverse of habitat suitability were not supported. This suggests that, although habitat associations of many butterflies are specific due to reproductive requirements, habitat suitability and landscape permeability are not equivalent concepts due to considerable adult vagility. We infer that divergent selection related to variation in summer temperatures has produced two genetic clusters within P. m. dodi, differing in voltinism and diapause propensity. Within the next century, temperatures are predicted to rise by amounts greater than the present-day difference between regions of the genetic clusters, potentially affecting the persistence of the northern cluster under continued climate change.

Spatial predictors of genomic and phenotypic variation differ in a lowland Middle American bird (Icterus gularis)

Spatial patterns of intraspecific variation are shaped by geographical distance among populations, historical changes in gene flow and interactions with local environments. Although these factors are not mutually exclusive and operate on both genomic and phenotypic variation, it is unclear how they affect these two axes of variation. We address this question by exploring the predictors of genomic and phenotypic divergence in Icterus gularis, a broadly distributed Middle American bird that exhibits marked geographical variation in body size across its range. We combined a comprehensive single nucleotide polymorphism and phenotypic data set to test whether genome-wide genetic and phenotypic differentiation are best explained by (i) isolation by distance, (ii) isolation by history or (iii) isolation by environment. We find that the pronounced genetic and phenotypic variation in I. gularis are only partially correlated and differ regarding spatial predictors. Whereas genomic variation is largely explained by historical barriers to gene flow, phenotypic diversity can be best predicted by contemporary environmental heterogeneity. Our genomic analyses reveal strong phylogeographical structure coinciding with the Chivela Pass at the Isthmus of Tehuantepec that was formed during the Pleistocene, when populations were isolated in north-south refugia. In contrast, we found a strong association between body size and environmental variables, such as temperature and precipitation. The relationship between body size and local climate is consistent with a pattern produced by either natural selection or environmental plasticity. Overall, these results provide empirical evidence for why phenotypic and genomic data are often in conflict in taxonomic and phylogeographical studies.

Ecological, genetic and evolutionary drivers of regional genetic differentiation in Arabidopsis thaliana

BACKGROUND

Disentangling the drivers of genetic differentiation is one of the cornerstones in evolution. This is because genetic diversity, and the way in which it is partitioned within and among populations across space, is an important asset for the ability of populations to adapt and persist in changing environments. We tested three major hypotheses accounting for genetic differentiation-isolation-by-distance (IBD), isolation-by-environment (IBE) and isolation-by-resistance (IBR)-in the annual plant Arabidopsis thaliana across the Iberian Peninsula, the region with the largest genomic diversity. To that end, we sampled, genotyped with genome-wide SNPs, and analyzed 1772 individuals from 278 populations distributed across the Iberian Peninsula.

RESULTS

IBD, and to a lesser extent IBE, were the most important drivers of genetic differentiation in A. thaliana. In other words, dispersal limitation, genetic drift, and to a lesser extent local adaptation to environmental gradients, accounted for the within- and among-population distribution of genetic diversity. Analyses applied to the four Iberian genetic clusters, which represent the joint outcome of the long demographic and adaptive history of the species in the region, showed similar results except for one cluster, in which IBR (a function of landscape heterogeneity) was the most important driver of genetic differentiation. Using spatial hierarchical Bayesian models, we found that precipitation seasonality and topsoil pH chiefly accounted for the geographic distribution of genetic diversity in Iberian A. thaliana.

CONCLUSIONS

Overall, the interplay between the influence of precipitation seasonality on genetic diversity and the effect of restricted dispersal and genetic drift on genetic differentiation emerges as the major forces underlying the evolutionary trajectory of Iberian A. thaliana.

A ketocarotenoid-based colour polymorphism in the Sira poison frog Ranitomeya sirensis indicates novel gene interactions underlying aposematic signal variation

The accumulation of red ketocarotenoids is an important component of coloration in many organisms, but the underlying mechanisms are poorly understood. In some organisms, ketocarotenoids are sequestered from the diet and can accumulate when enzymes responsible for carotenoid breakdown are disrupted. In other organisms, ketocarotenoids are formed endogenously from dietary precursors via oxidation reactions carried out by carotenoid ketolase enzymes. Here, we study the genetic basis of carotenoid coloration in an amphibian. We demonstrate that a red/yellow polymorphism in the dendrobatid poison frog Ranitomeya sirensis is due to the presence/absence of ketocarotenoids. Using whole-transcriptome sequencing of skins and livers, we found that a transcript encoding a cytochrome P450 enzyme (CYP3A80) is expressed 3.4-fold higher in livers of red frogs versus yellow. As CYP3A enzymes are known carotenoid ketolases in other organisms, our results point to CYP3A80 as a strong candidate for a carotenoid ketolase in amphibians. Furthermore, in red frogs, the transcript encoding the carotenoid cleavage enzyme BCO2 is expressed at a low level or as a splice variant lacking key catalytic amino acids. This suggests that BCO2 function may be disrupted in red frogs, providing a mechanism whereby the accumulation of ketocarotenoids and their dietary precursors may be enhanced.

Environmental influences on evolvable robots

The field of Evolutionary Robotics addresses the challenge of automatically designing robotic systems. Furthermore, the field can also support biological investigations related to evolution. In this paper, we evolve (simulated) modular robots under diverse environmental conditions and analyze the influences that these conditions have on the evolved morphologies, controllers, and behavior. To this end, we introduce a set of morphological, controller, and behavioral descriptors that together span a multi-dimensional trait space. Using these descriptors, we demonstrate how changes in environmental conditions induce different levels of differentiation in this trait space. Our main goal is to gain deeper insights into the effect of the environment on a robotic evolutionary process.

Mechanisms for Color Convergence in a Mimetic Radiation of Poison Frogs

In animals, bright colors often evolve to mimic other species when a resemblance is selectively favored. Understanding the proximate mechanisms underlying such color mimicry can give insights into how mimicry evolves-for example, whether color convergence evolves from a shared set of mechanisms or through the evolution of novel color production mechanisms. We studied color production mechanisms in poison frogs (Dendrobatidae), focusing on the mimicry complex of Ranitomeya imitator. Using reflectance spectrometry, skin pigment analysis, electron microscopy, and color modeling, we found that the bright colors of these frogs, both within and outside the mimicry complex, are largely structural and produced by iridophores but that color production depends crucially on interactions with pigments. Color variation and mimicry are regulated predominantly by iridophore platelet thickness and, to a lesser extent, concentration of the red pteridine pigment drosopterin. Compared with each of the four morphs of model species that it resembles, R. imitator displays greater variation in both structural and pigmentary mechanisms, which may have facilitated phenotypic divergence in this species. Analyses of nonmimetic dendrobatids in other genera demonstrate that these mechanisms are widespread within the family and that poison frogs share a complex physiological "color palette" that can produce diverse and highly reflective colors.

The evolution of polymorphism in the warning coloration of the Amazonian poison frog Adelphobates galactonotus

While intraspecific variation in aposematic signals can be selected for by different predatory responses, their evolution is also contingent on other processes shaping genetic variation. We evaluate the relative contributions of selection, geographic isolation, and random genetic drift to the evolution of aposematic color polymorphism in the poison frog Adelphobates galactonotus, distributed throughout eastern Brazilian Amazonia. Dorsal coloration was measured for 111 individuals and genetic data were obtained from 220 individuals at two mitochondrial genes (mtDNA) and 7963 Single Nucleotide Polymorphisms (SNPs). Four color categories were described (brown, blue, yellow, orange) and our models of frog and bird visual systems indicated that each color was distinguishable for these taxa. Using outlier and correlative analyses we found no compelling genetic evidence for color being under divergent selection. A time-calibrated mtDNA tree suggests that the present distribution of dorsal coloration resulted from processes occurring during the Pleistocene. Separate phylogenies based on SNPs and mtDNA resolved the same well supported clades, each containing different colored populations. Ancestral character state analysis provided some evidence for evolutionary transitions in color type. Genetic structure was more strongly associated with geographic features, than color category, suggesting that the distribution of color is explained by localized processes. Evidence for geographic isolation together with estimates of low effective population size implicates drift as playing a key role in color diversification. Our results highlight the relevance of considering the neutral processes involved with the evolution of traits with important fitness consequences.

Fragmented and isolated: limited gene flow coupled with weak isolation by environment in the paleoendemic giant sequoia (Sequoiadendron giganteum)

PREMISE

Patterns of genetic structure across a species' range reflect the long-term interplay between genetic drift, gene flow, and selection. Given the importance of gene flow in preventing the loss of diversity through genetic drift among spatially isolated populations, understanding the dynamics of gene flow and the factors that influence connectivity across a species' range is a major goal for conservation of genetic diversity. Here we present a detailed look at gene flow dynamics of Sequoiadendron giganteum, a paleoendemic tree species that will likely face numerous threats due to climate change.

METHODS

We used microsatellite markers to examine nineteen populations of S. giganteum for patterns of genetic structure and to estimate admixture and rates of gene flow between eight population pairs. Also, we used Generalized Dissimilarity Models to elucidate landscape factors that shape genetic differentiation among populations.

RESULTS

We found minimal gene flow between adjacent groves in the northern disjunct range. In most of the southern portion of the range, groves showed a signal of connectivity which degrades to isolation in the extreme south. Geographic distance was the most important predictor of genetic dissimilarity across the range, with environmental conditions related to precipitation and temperature explaining a small, but significant, portion of the genetic variance.

CONCLUSIONS

Due to their isolation and unique genetic composition, northern populations of S. giganteum should be considered a high conservation priority. In this region, we suggest germplasm conservation as well as restoration planting to enhance genetic diversity.

Influence of environmental factors on the genetic variation of the aquatic macrophyte Ranunculus subrigidus on the Qinghai-Tibetan Plateau

BACKGROUND

Due to the environmental heterogeneity along elevation gradients, alpine ecosystems are ideal study objects for investigating how ecological variables shape the genetic patterns of natural species. The highest region in the world, the Qinghai-Tibetan Plateau, is a hotspot for the studies of evolutionary processes in plants. Many large rivers spring from the plateau, providing abundant habitats for aquatic and amphibious organisms. In the present study, we examined the genetic diversity of 13 Ranunculus subrigidus populations distributed throughout the plateau in order to elucidate the relative contribution of geographic distance and environmental dissimilarity to the spatial genetic pattern.

RESULTS

A relatively low level of genetic diversity within populations was found. No spatial genetic structure was suggested by the analyses of molecular variance, Bayesian clustering analysis and Mantel tests. Partial Mantel tests and multiple matrix regression analysis showed a significant influence of the environment on the genetic divergence of the species. Both climatic and water quality variables contribute to the habitat heterogeneity of R. subrigidus populations.

CONCLUSIONS

Our results suggest that historical processes involving long-distance dispersal and local adaptation may account for the genetic patterns of R. subrigidus and current environmental factors play an important role in the genetic differentiation and local adaptation of aquatic plants in alpine landscapes.

What drives genetic and phenotypic divergence in the Red-crowned Ant tanager (Habia rubica, Aves: Cardinalidae), a polytypic species?

AIM

The effects of geographic and environmental variables on patterns of genetic and phenotypic differentiation have been thoroughly studied. Ecological speciation involves reproductive isolation due to divergent natural selection that can result in a positive correlation between genetic divergence and adaptive phenotypic divergence (isolation by adaptation, IBA). If the phenotypic target of selection is unknown or not easily measured, environmental variation can be used as a proxy, expecting positive correlation between genetic and environmental distances, independent of geographic distances (isolation by environment, IBE). The null model is that the amount of gene flow between populations decreases as the geographic distance between them increases, and genetic divergence is due simply to the neutral effects of genetic drift (isolation by distance, IBD). However, since phenotypic differentiation in natural populations may be autocorrelated with geographic distance, it is often difficult to distinguish IBA from the neutral expectation of IBD. In this work, we test hypotheses of IBA, IBE, and IBD in the Red-crowned Ant tanager (Habia rubica).

LOCATION

Mesoamerica (Mexico-Central America) and South America.

TAXON

Habia rubica (Aves: Cardinalidae).

METHODS

We compiled genetic data, coloration, and morphometric data from specimens from collections in Mexico and the United States. We used the Multiple Matrix Regression with Randomization (MMRR) approach to evaluate the influence of geographic and environmental distances on genetic and phenotypic differentiation of H. rubica at both phylogroup and population levels.

RESULTS

Our results provide strong evidence that geographic distance is the main driver of genetic variation in H. rubica. We did not find evidence that climate variation is driving population differentiation in this species across a widespread geographic region.

MAIN CONCLUSIONS

Our data point to geographic isolation as the main factor structuring genetic variation within populations of H. rubica and suggest that climate is not playing a major role in genetic differentiation within this species.

Geographic isolation and climatic variability contribute to genetic differentiation in fragmented populations of the long-lived subalpine conifer Pinus cembra L. in the western Alps

Background

Genetic processes shape the modern-day distribution of genetic variation within and between populations and can provide important insights into the underlying mechanisms of evolution. The resulting genetic variation is often unequally partitioned within species’ distribution range and especially large differences can manifest at the range limit, where population fragmentation and isolation play a crucial role in species survival. Despite several molecular studies investigating the genetic diversity and differentiation of European Alpine mountain forests, the climatic and demographic constrains which influence the genetic processes are often unknown. Here, we apply non-coding microsatellite markers to evaluate the sporadic peripheral and continuous populations of cembra pine (Pinus cembra L.), a long-lived conifer species that inhabits the subalpine treeline ecotone in the western Alps to investigate how the genetic processes contribute to the modern-day spatial distribution. Moreover, we corroborate our findings with paleoecological records, micro and macro-remains, to infer the species’ possible glacial refugia and expansion scenarios.

Results

Four genetically distinct groups were identified, with Bayesian and F_ST based approaches, across the range of the species, situated in the northern, inner and south-western Alps. We found that genetic differentiation is substantially higher in marginal populations than at the center of the range, and marginal stands are characterized by geographic and genetic isolation due to spatial segregation and restricted gene flow. Moreover, multiple matrix regression approaches revealed effects of climatic heterogeneity in species’ spatial genetic pattern. Also, population stability tests indicated that all populations had experienced a severe historical bottleneck, no heterozygosity excess was detected, suggesting that more recently population sizes have remained relatively stable.

Conclusions

Our study demonstrated that cembra pine might have survived in multiple glacial refugia and subsequently recolonized the Alps by different routes. Modern-day marginal populations, at the edge of the species’ range, could maintain stable sizes over long periods without inbreeding depression and preserve high amounts of genetic variation. Moreover, our analyses indicate that climatic variability has played a major role in shaping differentiation, in addition to past historical events such as migration and demographic changes.

Invisible barriers: anthropogenic impacts on inter- and intra-specific interactions as drivers of landscape-independent fragmentation

Anthropogenically induced fragmentation constitutes a major threat to biodiversity. Presently, conservation research and actions focus predominantly on fragmentation caused directly by physical transformation of the landscape (e.g. deforestation, agriculture, urbanization, roads, etc.). While there is no doubt that landscape features play a key role in fragmenting populations or enhancing connectivity, fragmentation may also come about by processes other than the transformation of the landscape and which may not be readily visible. Such landscape-independent fragmentation (LIF) usually comes about when anthropogenic disturbance alters the inter- and intra-specific interactions among and within species. LIF and its drivers have received little attention in the scientific literature and in the management of wildlife populations. We discuss three major classes of LIF processes and their relevance for the conservation and management of species and habitats: (i) interspecific dispersal dependency, in which populations of species that rely on other species for transport and propagation become fragmented as the transporting species declines; (ii) interspecific avoidance induction, where species are excluded from habitats and corridors owing to interspecific interactions resulting from anthropogenically induced changes in community structure (e.g. exclusions by increased predation pressure); and (iii) intraspecific behavioural divergence, where populations become segregated owing to anthropogenically induced behavioural differentiation among them.

This article is part of the theme issue ‘Linking behaviour to dynamics of populations and communities: application of novel approaches in behavioural ecology to conservation'.

Divergent selection along elevational gradients promotes genetic and phenotypic disparities among small mammal populations

Species distributed along mountain slopes, facing contrasting habitats in short geographic scale, are of particular interest to test how ecologically based divergent selection promotes phenotypic and genetic disparities as well as to assess isolation-by-environment mechanisms. Here, we conduct the first broad comparative study of phenotypic variation along elevational gradients, integrating a large array of ecological predictors and disentangling population genetic driver processes. The skull form of nine ecologically distinct species distributed over a large altitudinal range (100-4200 m) was compared to assess whether phenotypic divergence is a common phenomenon in small mammals and whether it shows parallel patterns. We also investigated the relative contribution of biotic (competition and predation) and abiotic parameters on phenotypic divergence via mixed models. Finally, we assessed the population genetic structure of a rodent species (Niviventer confucianus) via analysis of molecular variance and FST along three mountain slopes and tested the isolation-by-environment hypothesis using Mantel test and redundancy analysis. We found a consistent phenotypic divergence and marked genetic structure along elevational gradients; however, the species showed mixed patterns of size and skull shape trends across mountain zones. Individuals living at lower altitudes differed greatly in both phenotype and genotype from those living at high elevations, while middle-elevation individuals showed more intermediate forms. The ecological parameters associated with phenotypic divergence along elevation gradients are partly related to species' ecological and evolutionary constraints. Fossorial and solitary animals are mainly affected by climatic factors, while terrestrial and more gregarious species are influenced by biotic and abiotic parameters. A novel finding of our study is that predator richness emerged as an important factor associated with the intraspecific diversification of the mammalian skull along elevational gradients, a previously overlooked parameter. Population genetic structure was mainly driven by environmental heterogeneity along mountain slopes, with no or a week spatial effect, fitting the isolation-by-environment scenario. Our study highlights the strong and multifaceted effects of heterogeneous steep habitats and ecologically based divergent selective forces in small mammal populations.