The influence of locus number and information content on species delimitation: an empirical test case in an endangered Mexican salamander

Perhaps the most important recent advance in species delimitation has been the development of model‐based approaches to objectively diagnose species diversity from genetic data. Additionally, the growing accessibility of next‐generation sequence data sets provides powerful insights into genome‐wide patterns of divergence during speciation. However, applying complex models to large data sets is time‐consuming and computationally costly, requiring careful consideration of the influence of both individual and population sampling, as well as the number and informativeness of loci on species delimitation conclusions. Here, we investigated how locus number and information content affect species delimitation results for an endangered Mexican salamander species, Ambystoma ordinarium. We compared results for an eight‐locus, 137‐individual data set and an 89‐locus, seven‐individual data set. For both data sets, we used species discovery methods to define delimitation models and species validation methods to rigorously test these hypotheses. We also used integrated demographic model selection tools to choose among delimitation models, while accounting for gene flow. Our results indicate that while cryptic lineages may be delimited with relatively few loci, sampling larger numbers of loci may be required to ensure that enough informative loci are available to accurately identify and validate shallow‐scale divergences. These analyses highlight the importance of striking a balance between dense sampling of loci and individuals, particularly in shallowly diverged lineages. They also suggest the presence of a currently unrecognized, endangered species in the western part of A. ordinarium's range.     

Genetic signatures through space,time and multiple disturbances in a ubiquitous brooding coral

The predominance of self‐recruitment in many reef‐building corals has fundamental and complex consequences for their genetic diversity, population persistence and responses to climate change. Knowledge of genetic structure over local scales needs to be placed within a broad spatial context, and also integrated with genetic monitoring through time to disentangle these consequences. Here, we examined patterns of genetic diversity over multiple spatio‐temporal scales across tropical Australia in the ubiquitous brooding coral, Seriatopora hystrix. We also analysed complimentary environmental and demographic data to elucidate the seascape drivers of these patterns. Large genetic differences were detected between the east vs. west coasts of Australia. In northwest Australia, geographic differentiation dominated genetic structure over multiple scales. However, three sympatric lineages were detected at the largest offshore reef system (Scott Reef). Similar to the differences observed among putative species in eastern Australia, these lineages were associated with different levels of wave exposure. Local genetic structure within the Scott Reef system was relatively stable over 10 years, but temporal differences were observed that reflected small but important genetic changes over a few generations during recovery after severe bleaching. These results highlight the importance of self‐recruitment together with occasional longer distance connectivity for the persistence of a metapopulation across spatially and temporally variable environments. Our multidimensional research provides a foundation for further long‐term genetic monitoring to inform conservation strategies and highlights that sampling scales, ecological effects and cryptic diversity are important considerations to develop realistic understanding of the evolutionary resilience of corals.     

Effective population size as a driver for divergence of an antimicrobial peptide (Hymenoptaecin) in two common European bumblebee species

Social insects are the target of numerous pathogens. This is because the high density of closely‐related individuals frequently interacting with each other enhances the transmission and establishment of pathogens. This high selective pressure results in the rapid evolution of immune genes, which might be counteracted by a reduced effective population size (N_e) lowering the effectiveness of selection. We tested the effect of N_e on the evolutionary rate of an important immune gene for the antimicrobial peptide Hymenoptaecin in two common central European bumblebee species: Bombus terrestris and Bombus lapidarius. Both species are similar in their biology and are expected to be under similar selective pressures because pathogen prevalence does not differ between species. However, previous studies indicated a higher N_e in B. terrestris compared to B. lapidarius. We found high intraspecific variability in the coding sequence but low variability for silent polymorphisms in B. lapidarius. Estimates of long‐ and short‐term N_e were three‐ to four‐fold higher N_e in B. terrestris, although the species did not differ in census population sizes. The difference in N_e might result in less efficient selection and suboptimal adaptation of immune genes (e.g. hymenoptaecin) in B. lapidarius, and thus this species might become less resistant and more tolerant, turning into a superspreader of diseases.     

Contrasting responses to climate change at Himalayan treelines revealed by population demographics of two dominant species

1. Alpine treelines are expected to shift upward due to recent climate change. However, interpretation of changes in montane systems has been problematic because effects of climate change are frequently confounded with those of land use changes. The eastern Himalaya, particularly Langtang National Park, Central Nepal, has been relatively undisturbed for centuries and thus presents an opportunity for studying climate change impacts on alpine treeline uncontaminated by potential confounding factors.
2. We studied two dominant species, Abies spectabilis (AS) and Rhododendron campanulatum (RC), above and below the treeline on two mountains. We constructed 13 transects, each spanning up to 400 m in elevation, in which we recorded height and state (dead or alive) of all trees, as well as slope, aspect, canopy density, and measures of anthropogenic and animal disturbance.
3. All size classes of RC plants had lower mortality above treeline than below it, and young RC plants (<2 m tall) were at higher density above treeline than below. AS shows little evidence of a position change from the historic treeline, with a sudden extreme drop in density above treeline compared to below. Recruitment, as measured by size–class distribution, was greater above treeline than below for both species but AS is confined to ~25 m above treeline whereas RC is luxuriantly growing up to 200 m above treeline.
4. Synthesis. Evidence suggests that the elevational limits of RC have shifted upward both because (a) young plants above treeline benefited from facilitation of recruitment by surrounding vegetation, allowing upward expansion of recruitment, and (b) temperature amelioration to mature plants increased adult survival. We predict that the current pure stand of RC growing above treeline will be colonized by AS that will, in turn, outshade and eventually relegate RC to be a minor component of the community, as is the current situation below the treeline.

     

Genomic and phenotypic signatures of climate adaptation in an Anolis lizard

Integrated knowledge on phenotype, physiology, and genomic adaptations is required to understand the effects of climate on evolution. The functional genomic basis of organismal adaptation to changes in the abiotic environment, its phenotypic consequences, and its possible convergence across vertebrates are still understudied. In this study, we use a comparative approach to verify predicted gene functions for vertebrate thermal adaptation with observed functions underlying repeated genomic adaptations in response to elevation in the lizard Anolis cybotes. We establish a direct link between recurrently evolved phenotypes and functional genomics of altitude‐related climate adaptation in three highland and lowland populations in the Dominican Republic. We show that across vertebrates, genes contained in this interactome are expressed within the brain, the endocrine system, and during development. These results are relevant to elucidate the effect of global climate change across vertebrates and might aid in furthering insight into gene–environment relationships under disturbances to homeostasis.     

Evaluating the ability of Bayesian clustering methods to detect hybridization and introgression using an empirical red wolf data set

Bayesian clustering methods have emerged as a popular tool for assessing hybridization using genetic markers. Simulation studies have shown these methods perform well under certain conditions; however, these methods have not been evaluated using empirical data sets with individuals of known ancestry. We evaluated the performance of two clustering programs, baps and structure , with genetic data from a reintroduced red wolf (Canis rufus) population in North Carolina, USA. Red wolves hybridize with coyotes (C. latrans), and a single hybridization event resulted in introgression of coyote genes into the red wolf population. A detailed pedigree has been reconstructed for the wild red wolf population that includes individuals of 50–100% red wolf ancestry, providing an ideal case study for evaluating the ability of these methods to estimate admixture. Using 17 microsatellite loci, we tested the programs using different training set compositions and varying numbers of loci. structure was more likely than baps to detect an admixed genotype and correctly estimate an individual's true ancestry composition. However, structure was more likely to misclassify a pure individual as a hybrid. Both programs were outperformed by a maximum‐likelihood‐based test designed specifically for this system, which never misclassified a hybrid (50–75% red wolf) as a red wolf or vice versa. Training set composition and the number of loci both had an impact on accuracy but their relative importance varied depending on the program. Our findings demonstrate the importance of evaluating methods used for detecting admixture in the context of endangered species management.     

Identifying refugia from climate change using coupled ecological and genetic data in a transitional Mediterranean‐temperate tree species

Populations occurring in areas of overlap between the current and future distribution of a species are particularly important because they can represent “refugia from climate change”. We coupled ecological and range‐wide genetic variation data to detect such areas and to evaluate the impacts of habitat suitability changes on the genetic diversity of the transitional Mediterranean‐temperate tree Fraxinus angustifolia. We sampled and genotyped 38 natural populations comprising 1006 individuals from across Europe. We found the highest genetic diversity in western and northern Mediterranean populations, as well as a significant west to east decline in genetic diversity. Areas of potential refugia that correspond to approximately 70% of the suitable habitat may support the persistence of more than 90% of the total number of alleles in the future. Moreover, based on correlations between Bayesian genetic assignment and climate, climate change may favour the westward spread of the Black Sea gene pool in the long term. Overall, our results suggest that the northerly core areas of the current distribution contain the most important part of the genetic variation for this species and may serve as in situ macrorefugia from ongoing climate change. However, rear‐edge populations of the southern Mediterranean may be exposed to a potential loss of unique genetic diversity owing to habitat suitability changes unless populations can persist in microrefugia that have facilitated such persistence in the past.     

Calling behaviour under climate change: geographical and seasonal variation of calling temperatures in ectotherms

Calling behaviour is strongly temperature‐dependent and critical for sexual selection and reproduction in a variety of ectothermic taxa, including anuran amphibians, which are the most globally threatened vertebrates. However, few studies have explored how species respond to distinct thermal environments at time of displaying calling behaviour, and thus it is still unknown whether ongoing climate change might compromise the performance of calling activity in ectotherms. Here, we used new audio‐trapping techniques (automated sound recording and detection systems) between 2006 and 2009 to examine annual calling temperatures of five temperate anurans and their patterns of geographical and seasonal variation at the thermal extremes of species ranges, providing insights into the thermal breadths of calling activity of species, and the mechanisms that enable ectotherms to adjust to changing thermal environments. All species showed wide thermal breadths during calling behaviour (above 15 °C) and increases in calling temperatures in extremely warm populations and seasons. Thereby, calling temperatures differed both geographically and seasonally, both in terrestrial and aquatic species, and were 8–22 °C below the specific upper critical thermal limits (CT_max) and strongly associated with the potential temperatures of each thermal environment (operative temperatures during the potential period of breeding). This suggests that calling behaviour in ectotherms may take place at population‐specific thermal ranges, diverging when species are subjected to distinct thermal environments, and might imply plasticity of thermal adjustment mechanisms (seasonal and developmental acclimation) that supply species with means of coping with climate change. Furthermore, the thermal thresholds of calling at the onset of the breeding season were dissimilar between conspecific populations, suggesting that other factors besides temperature are needed to trigger the onset of reproduction. Our findings imply that global warming would not directly inhibit calling behaviour in the study species, although might affect other temperature‐dependent features of their acoustic communication system.     

A genomic assessment of population structure and gene flow in an aquatic salamander identifies the roles of spatial scale,barriers, and river architecture

相似文献   

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

Gene flow in the anemone Anthopleura elegantissima limits signatures of local adaptation across an extensive geographic range

Species inhabiting marine environments face a wide range of environmental conditions that vary spatially across several orders of magnitude. The selective pressures that these conditions impose on marine organisms, in combination with potentially high rates of gene flow between distant populations, make it difficult to predict the extent to which these populations can locally adapt. Here, I identify how selection and gene flow influence the population genetic structure of the anemone Anthopleura elegantissima along the Pacific coast of North America. Isolation by distance is the dominant pattern across the range of this species, with a genetic break near Pt. Conception, CA. Furthermore, demographic modelling suggests that this species was historically confined to southerly latitudes before expanding northward. Outlier analyses identify 24 loci under selection (out of ~1,100), but the same analysis on simulated genetic data generated using the most likely demographic model erroneously identified the same number of loci under selection, if not more. Taken together, these results suggest that demographic processes are the dominant force shaping population genetic patterns in A. elegantissima along the Pacific coast of North America. I discuss these patterns in terms of the evolutionary history of A. elegantissima, the potential for local adaptation, and their consequences with respect to interactions with the endosymbiont Breviolum muscatinei across their geographic range.     

Detecting an elusive invasive species: a diagnostic PCR to detect Burmese python in Florida waters and an assessment of persistence of environmental DNA

Recent studies have demonstrated that detection of environmental DNA (eDNA) from aquatic vertebrates in water bodies is possible. The Burmese python, Python bivittatus, is a semi‐aquatic, invasive species in Florida where its elusive nature and cryptic coloration make its detection difficult. Our goal was to develop a diagnostic PCR to detect P. bivittatus from water‐borne eDNA, which could assist managers in monitoring this invasive species. First, we used captive P. bivittatus to determine whether reptilian DNA could be isolated and amplified from water samples. We also evaluated the efficacy of two DNA isolation methods and two DNA extraction kits commonly used in eDNA preparation. A fragment of the mitochondrial cytochrome b gene from P. bivittatus was detected in all water samples isolated with the sodium acetate precipitate and the QIAamp DNA Micro Kit. Next, we designed P. bivittatus‐specific primers and assessed the degradation rate of eDNA in water. Our primers did not amplify DNA from closely related species, and we found that P. bivittatus DNA was consistently detectable up to 96 h. Finally, we sampled water from six field sites in south Florida. Samples from five sites, where P. bivittatus has been observed, tested positive for eDNA. The final site was negative and had no prior documented evidence of P. bivittatus. This study shows P. bivittatus eDNA can be isolated from water samples; thus, this method is a new and promising technique for the management of invasive reptiles.     

Increased genetic divergence between two closely related fir species in areas of range overlap

Because of introgressive hybridization, closely related species can be more similar to each other in areas of range overlap (parapatry or sympatry) than in areas where they are geographically isolated from each other (allopatry). Here, we report the reverse situation based on nuclear genetic divergence between two fir species, Abies chensiensis and Abies fargesii, in China, at sites where they are parapatric relative to where they are allopatric. We examined genetic divergence across 126 amplified fragment length polymorphism (AFLP) markers in a set of 172 individuals sampled from both allopatric and parapatric populations of the two species. Our analyses demonstrated that AFLP divergence was much greater between the species when comparisons were made between parapatric populations than between allopatric populations. We suggest that selection in parapatry may have largely contributed to this increased divergence.     

Environmental gradients shape the genetic structure of two medicinal Salvia species in Jordan

• Environmental gradients, and particularly climatic variables, exert a strong influence on plant distribution and, potentially, population genetic diversity and differentiation. Differences in water availability can cause among‐population variation in ecological processes and can thus interrupt populations’ connectivity and isolate them environmentally. The present study examines the effect of environmental heterogeneity on plant populations due to environmental isolation unrelated to geographic distance.
• Using AFLP markers, we analyzed genetic diversity and differentiation among 12 Salvia spinosa populations and 13 Salvia syriaca populations from three phytogeographical regions (Mediterranean, Irano‐Turanian and Saharo‐Arabian) representing the extent of the species’ geographic range in Jordan. Differences in geographic location and climate were considered in the analyses.
• For both species, flowering phenology varied among populations and regions. Irano‐Turanian and Saharo‐Arabian populations had higher genetic diversity than Mediterranean populations, and genetic diversity increased significantly with increasing temperature. Genetic diversity in Salvia syriaca was affected by population size, while genetic diversity responded to drought in S. spinosa. For both species, high levels of genetic differentiation were found as well as two well‐supported phytogeographical groups of populations, with Mediterranean populations clustering in one group and the Irano‐Turanian and Saharo‐Arabian populations in another. Genetic distance was significantly correlated to environmental distance, but not to geographic distance.
• Our data indicate that populations from moist vs. arid environments are environmentally isolated, where environmental gradients affect their flowering phenology, limit gene flow and shape their genetic structure. We conclude that environmental heterogeneity may act as driver for the observed variation in genetic diversity.

     

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80‐year period. Screening of >1000 gene‐associated single‐nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.     

Evolution of sociality in spiders leads to depleted genomic diversity at both population and species levels

Across several animal taxa, the evolution of sociality involves a suite of characteristics, a “social syndrome,” that includes cooperative breeding, reproductive skew, primary female‐biased sex ratio, and the transition from outcrossing to inbreeding mating system, factors that are expected to reduce effective population size (Ne). This social syndrome may be favoured by short‐term benefits but come with long‐term costs, because the reduction in Ne amplifies loss of genetic diversity by genetic drift, ultimately restricting the potential of populations to respond to environmental change. To investigate the consequences of this social life form on genetic diversity, we used a comparative RAD‐sequencing approach to estimate genomewide diversity in spider species that differ in level of sociality, reproductive skew and mating system. We analysed multiple populations of three independent sister‐species pairs of social inbreeding and subsocial outcrossing Stegodyphus spiders, and a subsocial outgroup. Heterozygosity and within‐population diversity were sixfold to 10‐fold lower in social compared to subsocial species, and demographic modelling revealed a tenfold reduction in Ne of social populations. Species‐wide genetic diversity depends on population divergence and the viability of genetic lineages. Population genomic patterns were consistent with high lineage turnover, which homogenizes the genetic structure that builds up between inbreeding populations, ultimately depleting genetic diversity at the species level. Indeed, species‐wide genetic diversity of social species was 5–8 times lower than that of subsocial species. The repeated evolution of species with this social syndrome is associated with severe loss of genomewide diversity, likely to limit their evolutionary potential.