Bee Trypanosomatids: First Steps in the Analysis of the Genetic Variation and Population Structure of Lotmaria passim, Crithidia bombi and Crithidia mellificae

Trypanosomatids are among the most prevalent parasites in bees but, despite the fact that their impact on the colonies can be quite important and that their infectivity may potentially depend on their genotypes, little is known about the population diversity of these pathogens. Here we cloned and sequenced three non-repetitive single copy loci (DNA topoisomerase II, glyceraldehyde-3-phosphate dehydrogenase and RNA polymerase II large subunit, RPB1) to produce new genetic data from Crithidia bombi, C. mellificae and Lotmaria passim isolated from honeybees and bumblebees. These were analysed by applying population genetic tools in order to quantify and compare their variability within and between species, and to obtain information on their demography and population structure. The general pattern for the three species was that (1) they were subject to the action of purifying selection on nonsynonymous variants, (2) the levels of within species diversity were similar irrespective of the host, (3) there was evidence of recombination among haplotypes and (4) they showed no haplotype structuring according to the host. C. bombi exhibited the lowest levels of synonymous variation (π_S= 0.06 ± 0.04 %) - and a mutation frequency distribution compatible with a population expansion after a bottleneck - that contrasted with the extensive polymorphism displayed by C. mellificae (π_S= 2.24 ± 1.00 %), which likely has a more ancient origin. L. passim showed intermediate values (π_S= 0.40 ± 0.28 %) and an excess of variants a low frequencies probably linked to the spread of this species to new geographical areas.

An overview of current population genomics methods for the analysis of whole-genome resequencing data in eukaryotes

Characterizing the population history of a species and identifying loci underlying local adaptation is crucial in functional ecology, evolutionary biology, conservation and agronomy. The constant improvement of high-throughput sequencing techniques has facilitated the production of whole genome data in a wide range of species. Population genomics now provides tools to better integrate selection into a historical framework, and take into account selection when reconstructing demographic history. However, this improvement has come with a profusion of analytical tools that can confuse and discourage users. Such confusion limits the amount of information effectively retrieved from complex genomic data sets, and impairs the diffusion of the most recent analytical tools into fields such as conservation biology. It may also lead to redundancy among methods. To address these isssues, we propose an overview of more than 100 state-of-the-art methods that can deal with whole genome data. We summarize the strategies they use to infer demographic history and selection, and discuss some of their limitations. A website listing these methods is available at www.methodspopgen.com.

The differential genetic signatures related to climatic landscapes for jaguars and pumas on a continental scale

Modern and paleoclimate changes may have altered species dynamics by shifting species' niche suitability over space and time. We analyze whether the current genetic structure and isolation of the two large American felids, jaguar (Panthera onca) and puma (Puma concolor), are mediated by changes in climatic suitability and connection routes over modern and paleoclimatic landscapes. We estimate species distribution under 5 climatic landscapes (modern, Holocene, last maximum glaciations [LMG], average suitability, and climatic instability) and correlate them with individuals' genetic isolation through causal modeling on a resemblance matrix. Both species exhibit genetic isolation patterns correlated with LMG climatic suitability, suggesting that these areas may have worked as "allele refuges." However, the jaguar showed higher vulnerability to climate changes, responding to modern climatic suitability and connection routes, whereas the puma showed a continuous and gradual transition of genetic variation. Despite differential responsiveness to climate change, both species are subjected to the climatic effects on genetic configuration, which may make them susceptible to future climatic changes, since these are progressing faster and with higher intensity than changes in the paleoclimate. Thus, the effects of climatic changes should be considered in the design of conservation strategies to ensure evolutionary and demographic processes mediated by gene flow for both species.

Integrating Multiple Lines of Evidence to Explore Intraspecific Variability in a Rare Endemic Alpine Plant and Implications for Its Conservation

We studied the ecology, distribution, and phylogeography of Tephroseris balbisiana, a rare plant whose range is centered to the South-Western Alps. Our aim was to assess the extent of intraspecific variability within the nominal species and the conservation status of isolated populations. We studied genetic diversity across the whole species range. We analyzed leaf traits, which are distinctive morphological characters within the Tephroseris genus. A clear pattern of genetic variation was found among populations of T. balbisiana, which clustered according to their geographic position. On the contrary, there was a strong overlap in the morphological space of individuals across the species' range, with few peripheral populations diverging in their leaf morphology. Studying habitat suitability by means of species distribution models, we observed that T. balbisiana range is primarily explained by solar radiation and precipitation seasonality. Environmental requirements could explain the genetic and morphological uniformity of T. balbisiana in its core distribution area and justify genetic, morphological, and ecological divergences found among the isolated populations of the Apennines. Our findings emphasize the need to account for the whole diversity of a species, comprising peripheral populations, in order to better estimate its status and to prioritize areas for its conservation.

Morphological and Genetic Variation in Monocultures, Forestry Systems and Wild Populations of Agave maximiliana of Western Mexico: Implications for Its Conservation

Forestry systems in Mexico are examples of traditional management of land and biodiversity that integrates the use, conservation and restoration of forest elements. Current in situ management practices of Agave maximiliana in western Mexico include the tolerance of many forest elements, reintroduction of young Agave plants and germination of seeds. More intense forms of management include monocultures, which are agroindustrialized systems developed in more recent times and characterized by the establishment of high densities of A. maximiliana plants in deforested areas and abandoned agricultural lands. We compared monocultures, forestry systems and wild populations (i.e., non/slightly-exploited forests) in order to evaluate whether these practices have had an effect on intraspecific morphological and genetic variation and divergence. We also tested whether divergence has a positive relationship with environmental and geographic distance. We analyzed 16 phenotypic traits in 17 populations of A. maximiliana, and 14 populations were further examined by amplifying 9 SSR loci. We employed multivariate methods and analyses of variance in phenotypic and genetic traits to test whether clusters and the percentage of variation contained in the managed and wild categories can be identified. Tests of isolation by environment (IBE) and distance (IBD) were performed to detect the magnitude of divergence explained by climatic and geographic variables. We found that forestry systems are effective as reservoirs of morphological and genetic diversity, since they maintain levels similar to those of wild populations. Moreover, the monocultures showed similar levels, reflecting their recent emergence. While the species showed high morphological diversity (IMD = 0.638, SE ± 0.07), it had low to intermediate genetic diversity (A = 2.37, H _E = 0.418). Similar morphological and genetic divergences were found among populations, but these were not correlated with each other in population pairs. Non-significant morphological differentiation was found among categories. Only IBE was significant in the genetic structure (β = 0.32, p = 0.007), while neither IBE nor IBD was detected in the morphological differentiation. We discuss the implications of these results in the context of the weaknesses and strengths of A. maximiliana in the face of the socio-ecological changes predicted for the study area in the short term.

Effects of Sample Size and Full Sibs on Genetic Diversity Characterization: A Case Study of Three Syntopic Iberian Pond-Breeding Amphibians

Accurate characterization of genetic diversity is essential for understanding population demography, predicting future trends and implementing efficient conservation policies. For that purpose, molecular markers are routinely developed for nonmodel species, but key questions regarding sampling design, such as calculation of minimum sample sizes or the effect of relatives in the sample, are often neglected. We used accumulation curves and sibship analyses to explore how these 2 factors affect marker performance in the characterization of genetic diversity. We illustrate this approach with the analysis of an empirical dataset including newly optimized microsatellite sets for 3 Iberian amphibian species: Hyla molleri, Epidalea calamita, and Pelophylax perezi. We studied 17-21 populations per species (total n = 547, 652, and 516 individuals, respectively), including a reference locality in which the effect of sample size was explored using larger samples (77-96 individuals). As expected, FIS and tests for Hardy-Weinberg equilibrium and linkage disequilibrium were affected by the presence of full sibs, and most initially inferred disequilibria were no longer statistically significant when full siblings were removed from the sample. We estimated that to obtain reliable estimates, the minimum sample size (potentially including full sibs) was close to 20 for expected heterozygosity, and between 50 and 80 for allelic richness. Our pilot study based on a reference population provided a rigorous assessment of marker properties and the effects of sample size and presence of full sibs in the sample. These examples illustrate the advantages of this approach to produce robust and reliable results for downstream analyses.

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub‐legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west‐to‐east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.

Population Genetic Structure of the Endangered Kaiser's Mountain Newt, Neurergus kaiseri (Amphibia: Salamandridae)

Species often exhibit different levels of genetic structuring correlated to their environment. However, understanding how environmental heterogeneity influences genetic variation is difficult because the effects of gene flow, drift and selection are confounded. We investigated the genetic variation and its ecological correlates in an endemic and critically endangered stream breeding mountain newt, Neurergus kaiseri, within its entire range in southwestern Iran. We identified two geographic regions based on phylogenetic relationships using Bayesian inference and maximum likelihood of 779 bp mtDNA (D-loop) in 111 individuals from ten of twelve known breeding populations. This analysis revealed a clear divergence between northern populations, located in more humid habitats at higher elevation, and southern populations, from drier habitats at lower elevations regions. From seven haplotypes found in these populations none was shared between the two regions. Analysis of molecular variance (AMOVA) of N. kaiseri indicates that 94.03% of sequence variation is distributed among newt populations and 5.97% within them. Moreover, a high degree of genetic subdivision, mainly attributable to the existence of significant variance among the two regions is shown (θ_CT = 0.94, P = 0.002). The positive and significant correlation between geographic and genetic distances (r = 0.61, P = 0.002) following controlling for environmental distance suggests an important influence of geographic divergence of the sites in shaping the genetic variation and may provide tools for a possible conservation based prioritization policy for the endangered species.

Designing conservation strategies to preserve the genetic diversity of Astragalus edulis Bunge, an endangered species from western Mediterranean region

Astragalus edulis (Fabaceae) is an endangered annual species from the western Mediterranean region that colonized the SE Iberian Peninsula, NE and SW Morocco, and the easternmost Macaronesian islands (Lanzarote and Fuerteventura). Although in Spain some conservation measures have been adopted, it is still necessary to develop an appropriate management plan to preserve genetic diversity across the entire distribution area of the species. Our main objective was to use population genetics as well as ecological and phylogeographic data to select Relevant Genetic Units for Conservation (RGUCs) as the first step in designing conservation plans for A. edulis. We identified six RGUCs for in situ conservation, based on estimations of population genetic structure and probabilities of loss of rare alleles. Additionally, further population parameters, i.e. occupation area, population size, vulnerability, legal status of the population areas, and the historical haplotype distribution, were considered in order to establish which populations deserve conservation priority. Three populations from the Iberian Peninsula, two from Morocco, and one from the Canary Islands represent the total genetic diversity of the species and the rarest allelic variation. Ex situ conservation is recommended to complement the preservation of A. edulis, given that effective in situ population protection is not feasible in all cases. The consideration of complementary phylogeographic and ecological data is useful for management efforts to preserve the evolutionary potential of the species.