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

Temperature is a cryptic factor shaping the geographical pattern of genetic variation in Ceratophyllum demersum across a subtropical freshwater lake

Macrophyte habitats exhibit remarkable heterogeneity, encompassing the spatial variation of abiotic and biotic components such as changes in water conditions and weather as well as anthropogenic stressors. Environmental factors are thought to be important drivers shaping the genetic and epigenetic variation of aquatic plants. However, the links among genetic diversity, epigenetic variation, and environmental variables remain largely unclear, especially for clonal aquatic plants. Here, we performed population genetic and epigenetic analyses in conjunction with habitat discrimination to elucidate the environmental factors driving intraspecies genetic and epigenetic variation in hornwort (Ceratophyllum demersum) in a subtropical lake. Environmental factors were highly correlated with the genetic and epigenetic variation of C. demersum, with temperature being a key driver of the genetic variation. Lower temperature was detected to be correlated with greater genetic and epigenetic variation. Genetic and epigenetic variation were positively driven by water temperature, but were negatively affected by ambient air temperature. These findings indicate that the genetic and epigenetic variation of this clonal aquatic herb is not related to the geographic feature but is instead driven by environmental conditions, and demonstrate the effects of temperature on local genetic and epigenetic variation in aquatic systems.

Individual-based landscape genomics for conservation: An analysis pipeline

Landscape genomics can harness environmental and genetic data to inform conservation decisions by providing essential insights into how landscapes shape biodiversity. The massive increase in genetic data afforded by the genomic era provides exceptional resolution for answering critical conservation genetics questions. The accessibility of genomic data for non-model systems has also enabled a shift away from population-based sampling to individual-based sampling, which now provides accurate and robust estimates of genetic variation that can be used to examine the spatial structure of genomic diversity, population connectivity and the nature of environmental adaptation. Nevertheless, the adoption of individual-based sampling in conservation genetics has been slowed due, in large part, to concerns over how to apply methods developed for population-based sampling to individual-based sampling schemes. Here, we discuss the benefits of individual-based sampling for conservation and describe how landscape genomic methods, paired with individual-based sampling, can answer fundamental conservation questions. We have curated key landscape genomic methods into a user-friendly, open-source workflow, which we provide as a new R package, A Landscape Genomics Analysis Toolkit in R (algatr). The algatr package includes novel added functionality for all of the included methods and extensive vignettes designed with the primary goal of making landscape genomic approaches more accessible and explicitly applicable to conservation biology.

Late Pleistocene Altitudinal Segregation and Demography Define Future Climate Change Distribution of the Peromyscus mexicanus Species Group: Conservation Implications

Mountains harbor a significant number of the World's biodiversity, both on tropical and temperate regions. Notably, one crucial gap in conservation is the consideration of historical and contemporary patterns influencing differential distribution in small mammal mountain species and how climate change will affect their distribution and survival. The mice Peromyscus mexicanus species group is distributed across mountains in Guatemala-Chiapas and Central America, which experienced significant effects of glacial and interglacial cycles. We determined phylogeographic and demographic patterns of lowlands and highlands mountain lineages, revealing that the radiation of modern P. mexicanus lineages occurred during the Pleistocene (ca. 2.6 mya) along Nuclear Central America. In concert with climatic cycles and the distribution of habitats, lowland and highland lineages showed recent population size increase and decrease, respectively. We also estimated the current and future distribution ranges for six lineages, finding marked area size increase for two lineages for which vegetation type and distribution would facilitate migrating towards higher elevations. Contrastingly, three lineages showed range size decrease; their ecological requirements make them highly susceptible to future habitat loss. Our findings are clear evidence of the negative impacts of future climate change, while our ability to manage and conserve these vulnerable ecosystems and mountain species is contingent on our understanding of the implications of climate change on the distribution, ecology, and genetics of wildlife populations.

A wide hybrid zone mediated by precipitation contributed to confused geographical structure of Scutiger boulengeri

Confused geographical structure of a population and mitonuclear discordance are shaped by a combination of rapid changes in population demographics and shifts in ecology. In this study, we generated a time-calibrated phylogeny of Scutiger boulengeri, an endemic Xizang alpine toad occurring in mountain streams on the Qinghai-Xizang (Tibet) Plateau (QTP). Based on three mitochondrial DNA (mtDNA) genes, eight clades were assigned to three deeply divergent lineages. Analysis of nuclear DNA (nuDNA) genes revealed three distinct clusters without geographic structure, indicating significantly high rates of gene flow. Coalescent theory framework analysis (approximate Bayesian computation model DIYABC and Migrate-N) suggested that divergence of the main intraspecific clusters was the result of hybridization after secondary contact in the Holocene around 0.59 million years ago (Ma). The ratio of mtDNA F _ST (fixation index) to nuDNA F _ST was 2.3, thus failing to show male-biased dispersal. Geographic cline analysis showed that a wide hybrid zone was initially established in southwestern China, without significant reproductive isolation but with strong introgression in S. boulengeri, suggesting high hybrid fitness. Furthermore, mtDNA genes exhibited isolation by distance (IBD) while nuDNA genes exhibited significant isolation by environment (IBE). Results suggested that mitonuclear discordance may have initially been caused by geographic isolation, followed by precipitation-mediated hybridization, producing a wide hybrid zone and geographic structure confusion of nuDNA genes in S. boulengeri. This study indicated that complicated historical processes may have led to specific genetic patterns, with a specific climate factor facilitating gene flow in the system.

Sexual imprinting leads to speciation in locally adapted populations

Sexual imprinting is widespread in birds and other species but its existence requires explanation. Our results suggest that sexual imprinting leads to speciation in locally-adapted populations if a neutral mating cue-e.g., novel plumage coloration-arises through mutation. Importantly, the mating cue locus is not linked to adaptation loci. Local adaptation is a necessary precursor to speciation and occurs when evolution results in stable genetic polymorphisms with one allele predominating in some areas while others predominate elsewhere. Here we use a deterministic two-niche population genetic model to map the set of migration and selection rates for which polymorphic evolutionary outcomes, i.e., local adaptations, can occur. Approximate equations for the boundaries of the set of polymorphic evolutionary outcomes were derived by Bulmer (American Naturalist, 106, 254, 1972), but our results, obtained by deterministic simulation of the evolutionary process, show that one of Bulmer's equations is inaccurate except when the level of dominance is 0.5, and fails if one of the alleles is dominant. Having an accurate map of the set of migration and selection rates for which polymorphic evolutionary outcomes can occur, we then show using the model of Sibly et al. (Ecology and Evolution, 9, 13506, 2019) that local adaptation in all analyzed cases leads to speciation if a new neutral mating cue arises by mutation. We finish by considering how genome sequencing makes possible testing our model and its predictions.

Axes of multivariate sexual signal divergence among incipient species: Concordance with selection, genetic variation and phenotypic plasticity

Sexual signalling traits are often observed to diverge rapidly among populations, thereby playing a potentially key early role in the evolution of reproductive isolation. While often assumed to reflect divergent sexual selection among populations, patterns of sexual trait diversification might sometimes be biased along axes of standing additive genetic variation and covariation among trait components. Additionally, theory predicts that environmentally induced phenotypic variation might facilitate rapid trait evolution, suggesting that patterns of divergence between populations should mirror phenotypic plasticity within populations. Here, we evaluate the concordance between observed axes of multivariate sexual trait divergence and predicted divergence based on (1) interpopulation variation in sexual selection, (2) additive genetic variances and (3) temperature-related phenotypic plasticity in male courtship song among geographically isolated populations of the Hawaiian swordtail cricket, Laupala cerasina, which exhibit sexual isolation due acoustic signalling traits. The major axis of multivariate divergence, d_max , accounted for 76% of variation among population male song trait means and was moderately correlated with interpopulation differences in directional sexual selection based on female preferences. However, the majority of additive genetic variance was largely oriented away from the direction of divergence, suggesting that standing genetic variation may not play a dominant role in the patterning of signal divergence. In contrast, the axis of phenotypic plasticity strongly mirrored patterns of interpopulation phenotypic divergence, which is consistent with a role for temperature-related plasticity in facilitating instead of inhibiting male song evolution and sexual isolation in these incipient species. We propose potential mechanisms by which sexual selection might interact with phenotypic plasticity to facilitate the rapid acoustic diversification observed in this species and clade.