Isolation by distance and sharp discontinuities in gene frequencies: implications for the phylogeography of an alpine insect species, Carabus solieri

Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

Molecular and quantitative genetic variation within and between populations of the declining grassland species Saxifraga granulata

Formerly common plant species are expected to be particularly susceptible to recent habitat fragmentation. We studied the population genetics of 19 recently fragmented Saxifraga granulata populations (max. distance 61 km) in Luxembourg and neighboring Germany using RAPD markers and a common garden experiment. We assessed (1) the relationships between plant fitness, quantitative genetic variation, molecular genetic variation, and population size; and (2) the relative importance of genetic drift and selection in shaping genetic variation. Molecular genetic diversity was high but did not correlate with population size, habitat conditions, or plant performance. Genetic differentiation was low (F ST = 0.079 ± 0.135), and there was no isolation by distance. Longevity, clonality, and the long-lived seed bank of S. granulata may have prevented strong genetic erosion and genetic differentiation among populations. However, genetic distinctness increased with decreasing genetic diversity indicating that random genetic drift occurred in the studied populations. Quantitative and molecular genetic variations were correlated, and their differentiation (Q ST vs. F ST) among S. granulata populations was similar, suggesting that mainly random processes have shaped the quantitative genetic differentiation among populations. However, pairwise quantitative genetic distances increased with geographic and climatic distances, even when adjusted for molecular genetic distances, indicating diversifying selection. Our results indicate that long-lived clonal species may be buffered at least temporarily against the negative effects of fragmentation. The relationship between quantitative genetic and geographic distance may be a more sensitive indicator of selection than Q ST-F ST differences.

Fine-scale genetic structure reflects limited and coordinated dispersal in the colonial monk parakeet, Myiopsitta monachus

The genetic structure of animal populations has considerable behavioural, ecological and evolutionary implications and may arise from various demographic traits. Here, we use observational field data and molecular genetics to determine the genetic structure of an invasive population of monk parakeets, Myiopsitta monachus, at a range of spatial scales, and investigate the demographic processes that generate the observed structure. Monk parakeets construct large nests that can house several pairs occupying separate chambers; these nests are often aggregated within nesting trees. We determined patterns of relatedness within compound nests, within nesting trees and between trees. Spatial autocorrelation analyses of pairwise genetic relatedness revealed fine-scale genetic structure with relatives of both sexes spatially clustered within, but not beyond, nesting trees. In addition, males were more related to males sharing their compound nests than to other males occupying the same nesting tree. By contrast, males and females within compound nests were not significantly more closely related than elsewhere in the same tree, and we found no evidence for inbreeding. Adults showed high breeding site fidelity between years despite considerable disturbance of nest sites. Natal dispersal was female-biased, but dispersal distances were relatively short with some natal philopatry observed in both sexes. Sibling coalitions, typically of males, were observed amongst both philopatric and dispersing birds. Our results show significant clustering of kin within compound nests and nesting trees resulting from limited and coordinated natal dispersal, with subsequent breeding site fidelity. The resulting genetic structure has implications for social behaviour in this unusual parrot species.

Habitat fragmentation differentially shapes neutral and immune gene variation in a tropical bird species

Habitat fragmentation is a major cause of biodiversity loss, responsible for an alteration of intraspecific patterns of neutral genetic diversity and structure. Although neutral genetic variation can be informative for demographic inferences, it may be a poor predictor of adaptive genetic diversity and thus of the consequences of habitat fragmentation on selective evolutionary processes. In this context, we contrasted patterns of genetic diversity and structure of neutral loci (microsatellites) and immune genes (i.e., toll-like receptors) in an understorey bird species, the wedge-billed woodcreeper Glyphorynchus spirurus. The objectives were (1) to investigate forest fragmentation effects on population genetic diversity, (2) to disentangle the relative role of demography (genetic drift and migration) and selection, and (3) to assess whether immunogenetic patterns could be associated with variation of ectoparasite (i.e., ticks) pressures. Our results revealed an erosion of neutral genetic diversity and a substantial genetic differentiation among fragmented populations, resulting from a decrease in landscape connectivity and leading to the divergence of distinct genetic pools at a small spatial scale. Patterns of genetic diversity observed for TLR4 and TLR5 were concordant with neutral genetic patterns, whereas those observed for TLR3 and TLR21 were discordant. This result underlines that the dominant evolutionary force shaping immunogenetic diversity (genetic drift vs. selection) may be different depending on loci considered. Finally, tick prevalence was higher in fragmented environments. We discussed the hypothesis that pathogen selective pressures may contribute to maintain adaptive genetic diversity despite the negative demographic effect of habitat fragmentation on neutral genetic diversity.

Microgeographical adaptation corresponds to elevational distributions of congeneric montane grasshoppers

Local adaptation can occur at small spatial scales relative to the dispersal capacity of species. Alpine ecosystems have sharp environmental clines that offer an opportunity to investigate the effects of fine-scale shifts in species' niche breadth on adaptive genetic processes. Here we examine two grasshopper species endemic to the Australian Alps (Kosciuscola spp.) that differ in elevational niche breadth: one broader, K. usitatus (1400-2200 m), and one narrower, K. tristis (1600-2000 m). We examine signatures of selection with respect to environmental and morphological variables in two mountain regions using F_ST outlier tests and environmental association analyses (EAAs) applied to single nucleotide polymorphism (SNP) data (K. usitatus: 9017 SNPs, n = 130; K. tristis: 7363 SNPs, n = 135). Stronger genetic structure was found in the more narrowly distributed K. tristis, which showed almost twice the number of SNPs under putative selection (10.8%) compared with K. usitatus (5.3%). When examining SNPs in common across species (n = 3058), 260 SNPs (8.5%) were outliers shared across species, and these were mostly associated with elevation, a proxy for temperature, suggesting parallel adaptive processes in response to climatic drivers. Additive polygenic scores (an estimate of the cumulative signal of selection across all candidate loci) were nonlinearly and positively correlated with elevation in both species. However, a steeper correlation in K. tristis indicated a stronger signal of spatially varying selection towards higher elevations. Our study illustrates that the niche breadth of co-occurring and related species distributed along the same environmental cline is associated with differences in patterns of microgeographical adaptation.

Genealogical Position of Japanese Populations of the Globally Distributed Mayfly Cloeon dipterum and Related Species (Ephemeroptera, Baetidae): A Molecular Phylogeographic Analysis

In the present study, we add genetic data of the mayfly Cloeon dipterum collected from the Japanese Islands to the established molecular phylogenetic knowledge in the mitochondrial COI gene of Cloeon mayflies. Cloeon dipterum is a typical cosmopolitan species that includes six intraspecific haplotype groups. The present phylogenetic analysis revealed that haplotypes of the Japanese C. dipterum constitute a seventh group together with a haplotype from Korea. This East Asian group forms a sister group with previously known European and North American haplotype groups (i.e., the clade CT1 to CT3). The present phylogenetic analysis further revealed the occurrence of two described species (C. dipterum and Cloeon ryogokuensis) and possibly three species (Cloeon sp. 1 to 3) in Japan. Consideration is given to the degree of genetic differentiation, divergence time, and differentiation process among these seven genetic groups.

Habitat preference differentiates the Holocene range dynamics but not barrier effects on two sympatric, congeneric trees (Tristaniopsis, Myrtaceae)

Niche partitioning can lead to differences in the range dynamics of plant species through its impacts on habitat availability, dispersal, or selection for traits that affect colonization and persistence. We investigated whether niche partitioning into upland and riparian habitats differentiates the range dynamics of two closely related and sympatric eastern Australian trees: the mountain water gum (Tristaniopsis collina) and the water gum (T. laurina). Using genomic data from SNP genotyping of 480 samples, we assessed the impact of biogeographic barriers and tested for signals of range expansion. Circuit theory was used to model isolation-by-resistance across three palaeo-environment scenarios: the Last Glacial Maximum, the Holocene Climate Optimum and present-day (1950-2014). Both trees showed similar genetic structure across historically dry barriers, despite evidence of significant environmental niche differentiation and different post-glacial habitat shifts. Tristaniopsis collina exhibits the signature of serial founder effects consistent with recent or rapid range expansion, whilst T. laurina has genetic patterns consistent with long-term persistence in geographically isolated populations despite occupying a broader bioclimatic niche. We found the minor influence of isolation-by-resistance on both species, though other unknown factors appear to shape genetic variation. We postulate that specialized recruitment traits (adapted to flood-disturbance regimes) rather than habitat availability limited post-glacial range expansion in T. laurina. Our findings indicate that niche breadth does not always facilitate range expansion through colonization and migration across barriers, though it can promote long-term persistence in situ.

Mating barriers between genetically divergent strains of the parasitic nematode Haemonchus contortus suggest incipient speciation

Haemonchus contortus, in common with many nematode species, has extremely high levels of genetic variation within and between field populations derived from distant geographical locations. MHco10(CAVR), MHco3(ISE) and MHco4(WRS) are genetically divergent H. contortus strains, originally derived from Australia, Kenya and South Africa, respectively, that have been maintained by numerous rounds of in vivo experimental infection of sheep. In order to explore potential pre-zygotic competition or post-zygotic incompatibility between the strains, we have investigated the ability of MHco10(CAVR) to interbreed with either MHco3(ISE) or MHco4(WRS) during dual strain co-infections. Sheep were experimentally co-infected with 4000 infective larvae (L₃) per os of the MHco10(CAVR) strain and an equal number of either the MHco3(ISE) or the MHco4(WRS) strain L₃. The adult worm establishement rates and the proportions of F₁ progeny resulting from intra- and inter-strain mating events were determined by admixture analysis of microsatellite multi-locus genotypes. Although there was no difference in adult worm establishment rates, the proportions of F₁ progeny of both the MHco10(CAVR) × MHco3(ISE) and MHco10(CAVR) × MHco4(WRS) dual strain co-infections departed from Mendelian expectations. The proportions of inter-strain hybrid F₁ progeny were lower than the expected 50%, suggesting either pre-zygotic competition or post-zygotic incompatibility between the co-infecting strains. To investigate this further, both eggs and hatched L₁ of broods from single adult female worms recovered from each dual co-infection were genotyped. Unhatched eggs from the broods revealed no inter-strain hybrid genotype deficit, suggesting there is no pre-zygotic competition between the strains. In contrast, there was a deficit in L₁ inter-strain hybrid genotypes in the broods derived from MHco3(ISE) or MHco4(WRS) maternal parents, but not from MHco10(CAVR) maternal parents. This suggests that hybrid progeny of MHco10(CAVR) paternal parents have reduced post-zygotic development and/or viability consistent with incipient speciation of the MHco10(CAVR) strain. The presence of mating barriers between allopatric H. contortus strains has important implications for parasite ecology, including the ability of newly introduced anthelmintic-resistant parasite populations to compete and interbreed with populations already established in a region.

Next-generation phylogeography of the cockle Cerastoderma glaucum: Highly heterogeneous genetic differentiation in a lagoon species

AIM

Coastal lagoons form an intriguing example of fragmented marine habitats. Restricted gene flow among isolated populations of lagoon species may promote their genetic divergence and may thus provide a first step toward speciation. In the present study, the population genetic structure of the lagoon cockle Cerastoderma glaucum has been investigated to clarify the complex phylogeographic pattern found in previous studies, to localize major genetic breaks, and to discuss their origin and maintenance.

LOCATION

The Atlantic and Mediterranean coasts, including the Baltic, North Sea, and Black Sea.

METHODS

A total of 204 C. glaucum individuals from 14 populations were genotyped using restriction site-associated DNA sequencing (RADseq). The genetic diversity, divergence, and structure were analyzed using genome-wide single nucleotide polymorphisms (SNPs). Phylogenetic relationships were inferred under a coalescent model using svdquartets.

RESULTS

The RADseq approach allowed inferring phylogeographic relationships with an unprecedented resolution. Three deeply divergent lineages were identified within C. glaucum that are separated by many genetic barriers: one lineage in the Aegean-Black Sea region, one in the Ionian Sea, and the last one widely distributed from the Western Mediterranean to the Baltic Sea. The nested branching pattern displayed on the species tree largely agrees with the likely scenario of C. glaucum postglacial expansion from the Mediterranean to the Baltic Sea.

MAIN CONCLUSION

The genetic differentiations between geographically separated lagoons proved to be strong, highlighting the evolutionary influence of these naturally fragmented habitats. The postglacial expansion created complex patterns of spatial segregation of genetic diversity with allele frequency gradients in many outlier loci, but also discrepancies between the nuclear and mitochondrial genetic markers that probably arose from genetic surfing of mitochondrial variation.

Genetic divergence, population differentiation and phylogeography of the cicada Subpsaltria yangi based on molecular and acoustic data: an example of the early stage of speciation?

Background

Geographical isolation combined with historical climatic fluctuations have been identified as two major factors that contribute to the formation of new species. On the other hand, biotic factors such as competition and predation are also able to drive the evolution and diversification of organisms. To determine whether geographical barriers contributed to population divergence or speciation in the rare endemic cicada Subpsaltria yangi the population differentiation, genetic structure and phylogeography of the species were investigated in the Loess Plateau and adjacent areas of northwestern China by analysing mitochondrial and nuclear DNA and comparing the calling song structure of 161 male individuals.

Results

The results reveal a low level of genetic differentiation and relatively simple phylogeographic structure for this species, but two independent clades corresponding to geographically isolated populations were recognised. Genetic and geographical distances were significantly correlated among lineages. Results of divergence-time estimation are consistent with a scenario of isolation due to glacial refugia and interglacial climate oscillation in northwestern China. Significant genetic divergence was found between the population occurring in the Helan Mountains and other populations, and recent population expansion has occurred in the Helan Mountains and/or adjacent areas. This population is also significantly different in calling song structure from other populations.

Conclusions

Geographical barriers (i.e., the deserts and semi-deserts surrounding the Helan Mountains), possibly coupled with related ecological differences, may have driven population divergence and allopatric speciation. This provides a possible example of incipient speciation in Cicadidae, improves understanding of population differentiation, acoustic signal diversification and phylogeographic relationships of this rare cicada species of conservation concern, and informs future studies on population differentiation, speciation and phylogeography of other insects with a high degree of endemism in the Helan Mountains and adjacent areas.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1317-8) contains supplementary material, which is available to authorized users.

Partitioning drivers of spatial genetic variation for a continuously distributed population of boreal caribou: Implications for management unit delineation

Isolation by distance (IBD) is a natural pattern not readily incorporated into theoretical models nor traditional metrics for differentiating populations, although clinal genetic differentiation can be characteristic of many wildlife species. Landscape features can also drive population structure additive to baseline IBD resulting in differentiation through isolation-by-resistance (IBR). We assessed the population genetic structure of boreal caribou across western Canada using nonspatial (STRUCTURE) and spatial (MEMGENE) clustering methods and investigated the relative contribution of IBD and IBR on genetic variation of 1,221 boreal caribou multilocus genotypes across western Canada. We further introduced a novel approach to compare the partitioning of individuals into management units (MU) and assessed levels of genetic connectivity under different MU scenarios. STRUCTURE delineated five genetic clusters while MEMGENE identified finer-scale differentiation across the study area. IBD was significant and did not differ for males and females both across and among detected genetic clusters. MEMGENE landscape analysis further quantified the proportion of genetic variation contributed by IBD and IBR patterns, allowing for the relative importance of spatial drivers, including roads, water bodies, and wildfires, to be assessed and incorporated into the characterization of population structure for the delineation of MUs. Local population units, as currently delineated in the boreal caribou recovery strategy, do not capture the genetic variation and connectivity of the ecotype across the study area. Here, we provide the tools to assess fine-scale spatial patterns of genetic variation, partition drivers of genetic variation, and evaluate the best management options for maintaining genetic connectivity. Our approach is highly relevant to vagile wildlife species that are of management and conservation concern and demonstrate varying degrees of IBD and IBR with clinal spatial genetic structure that challenges the delineation of discrete population boundaries.

Level of genetic differentiation affects relative performances of expressed sequence tag and genomic SSRs

Microsatellites, also called simple sequence repeats (SSRs), are markers of choice to estimate relevant parameters for conservation genetics, such as migration rates, effective population size and kinship. Cross-amplification of SSRs is the simplest way to obtain sets of markers, and highly conserved SSRs have recently been developed from expressed sequence tags (EST) to improve SSR cross-species utility. As EST-SSRs are located in coding regions, the higher stability of their flanking regions reduces the frequency of null alleles and improves cross-species amplification. However, EST-SSRs have generally less allelic variability than genomic SSRs, potentially leading to differences in estimates of population genetic parameters such as genetic differentiation. To assess the potential of EST-SSRs in studies of within-species genetic diversity, we compared the relative performance of EST- and genomic SSRs following a multispecies approach on passerine birds. We tested whether patterns and levels of genetic diversity within and between populations assessed from EST- and from genomic SSRs are congruent, and we investigated how the relative efficiency of EST- and genomic SSRs is influenced by levels of differentiation. EST- and genomic SSRs ensured comparable inferences of population genetic structure in cases of strong genetic differentiation, and genomic SSRs performed slightly better than EST-SSRs when differentiation is moderate. However and interestingly, EST-SSRs had a higher power to detect weak genetic structure compared to genomic SSRs. Our study attests that EST-SSRs may be valuable molecular markers for conservation genetic studies in taxa such as birds, where the development of genomic SSRs is impeded by their low frequency.

Population subdivision of hydrothermal vent polychaete Alvinella pompejana across equatorial and Easter Microplate boundaries

Background

The Equator and Easter Microplate regions of the eastern Pacific Ocean exhibit geomorphological and hydrological features that create barriers to dispersal for a number of animals associated with deep-sea hydrothermal vent habitats. This study examined effects of these boundaries on geographical subdivision of the vent polychaete Alvinella pompejana. DNA sequences from one mitochondrial and eleven nuclear genes were examined in samples collected from ten vent localities that comprise the species’ known range from 23°N latitude on the East Pacific Rise to 38°S latitude on the Pacific Antarctic Ridge.

Results

Multi-locus genotypes inferred from these sequences clustered the individual worms into three metapopulation segments — the northern East Pacific Rise (NEPR), southern East Pacific Rise (SEPR), and northeastern Pacific Antarctic Ridge (PAR) — separated by the Equator and Easter Microplate boundaries. Genetic diversity estimators were negatively correlated with tectonic spreading rates. Application of the isolation-with-migration (IMa2) model provided information about divergence times and demographic parameters. The PAR and NEPR metapopulation segments were estimated to have split roughly 4.20 million years ago (Mya) (2.42–33.42 Mya, 95 % highest posterior density, (HPD)), followed by splitting of the SEPR and NEPR segments about 0.79 Mya (0.07–6.67 Mya, 95 % HPD). Estimates of gene flow between the neighboring regions were mostly low (2 Nm < 1). Estimates of effective population size decreased with southern latitudes: NEPR > SEPR > PAR.

Conclusions

Highly effective dispersal capabilities allow A. pompejana to overcome the temporal instability and intermittent distribution of active hydrothermal vents in the eastern Pacific Ocean. Consequently, the species exhibits very high levels of genetic diversity compared with many co-distributed vent annelids and mollusks. Nonetheless, its levels of genetic diversity in partially isolated populations are inversely correlated with tectonic spreading rates. As for many other vent taxa, this pioneering colonizer is similarly affected by local rates of habitat turnover and by major dispersal filters associated with the Equator and the Easter Microplate region.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0807-9) contains supplementary material, which is available to authorized users.

Genetic differentiation in red-bellied piranha populations (Pygocentrus nattereri, Kner, 1858) from the Solimões-Amazonas River

Red-bellied piranhas (Pygocentrus nattereri) are widely caught with different intensities throughout the region of Solimões-Amazonas River by local fishermen. Thus, the management of this resource is performed in the absence of any information on its genetic stock. P. nattereri is a voracious predator and widely distributed in the Neotropical region, and it is found in other regions of American continent. However, information about genetic variability and structure of wild populations of red-bellied piranha is unavailable. Here, we describe the levels of genetic diversity and genetic structure of red-bellied piranha populations collected at different locations of Solimões-Amazonas River system. We collected 234 red-bellied piranhas and analyzed throughout eight microsatellite markers. We identified high genetic diversity within populations, although the populations of lakes ANA, ARA, and MAR have shown some decrease in their genetic variability, indicating overfishing at these communities. Was identified the existence of two biological populations when the analysis was taken altogether at the lakes of Solimões-Amazonas River system, with significant genetic differentiation between them. The red-bellied piranha populations presented limited gene flow between two groups of populations, which were explained by geographical distance between these lakes. However, high level of gene flow was observed between the lakes within of the biological populations. We have identified high divergence between the Catalão subpopulation and all other subpopulations. We suggest the creation of sustainable reserve for lakes near the city of Manaus to better manage and protect this species, whose populations suffer from both extractive and sport fishing.

Habitat specialization predicts genetic response to fragmentation in tropical birds

Habitat fragmentation is one of the most severe threats to biodiversity as it may lead to changes in population genetic structure, with ultimate modifications of species evolutionary potential and local extinctions. Nonetheless, fragmentation does not equally affect all species and identifying which ecological traits are related to species sensitivity to habitat fragmentation could help prioritization of conservation efforts. Despite the theoretical link between species ecology and extinction proneness, comparative studies explicitly testing the hypothesis that particular ecological traits underlies species-specific population structure are rare. Here, we used a comparative approach on eight bird species, co-occurring across the same fragmented landscape. For each species, we quantified relative levels of forest specialization and genetic differentiation among populations. To test the link between forest specialization and susceptibility to forest fragmentation, we assessed species responses to fragmentation by comparing levels of genetic differentiation between continuous and fragmented forest landscapes. Our results revealed a significant and substantial population structure at a very small spatial scale for mobile organisms such as birds. More importantly, we found that specialist species are more affected by forest fragmentation than generalist ones. Finally, our results suggest that even a simple habitat specialization index can be a satisfying predictor of genetic and demographic consequences of habitat fragmentation, providing a reliable practical and quantitative tool for conservation biology.

Contrasted patterns of genetic differentiation across eight bird species in the Lesser Antilles

Archipelagoes are considered as "natural laboratories" for studying processes that shape the distribution of diversity. The Lesser Antilles provide a favorable geographical context for divergence to occur. However, although morphological subspecies have been described across this archipelago in numerous avian species, the potential for the Lesser Antilles in driving intra-specific genetic divergence in highly mobile organisms such as birds remains understudied. Here, we assessed level of intra-specific genetic diversity and differentiation between three islands of the Lesser Antilles (Guadeloupe, Dominica and Martinique) using a multi-species approach on eight bird species. For each species, we built a set of microsatellite markers from cross-species amplifications. Significant patterns of inter-island and/or within-island genetic differentiation were detected in all species. However, levels of intra-specific genetic differentiation among the eight bird species were not always consistent with the boundaries of subspecies previously described in the sampled islands. These results suggest different histories of colonization/expansion and/or different species-specific ecological traits affecting gene flow, advocating for multi-species studies of historical and contemporary factors shaping the distribution of diversity on islands.

Spatial structuring of the population genetics of a European subterranean termite species

In population genetics studies, detecting and quantifying the distribution of genetic variation can help elucidate ecological and evolutionary processes. In social insects, the distribution of population-level genetic variability is generally linked to colony-level genetic structure. It is thus especially crucial to conduct complementary analyses on such organisms to examine how spatial and social constraints interact to shape patterns of intraspecific diversity. In this study, we sequenced the mitochondrial COII gene for 52 colonies of the subterranean termite Reticulitermes grassei (Isoptera: Rhinotermitidae), sampled from a population in southwestern France. Three haplotypes were detected, one of which was found exclusively in the southern part of the study area (near the Pyrenees). After genotyping 6 microsatellite loci for 512 individual termites, we detected a significant degree of isolation by distance among individuals over the entire range; however, the cline of genetic differentiation was not continuous, suggesting the existence of differentiated populations. A spatial principal component analysis based on allele frequency data revealed significant spatial autocorrelation among genotypes: the northern and southern groups were strongly differentiated. This finding was corroborated by clustering analyses; depending on the randomized data set, two or three clusters, exhibiting significant degrees of differentiation, were identified. An examination of colony breeding systems showed that colonies containing related neotenic reproductives were prevalent, suggesting that inbreeding may contribute to the high level of homozygosity observed and thus enhance genetic contrasts among colonies. We discuss the effect of evolutionary and environmental factors as well as reproductive and dispersal modes on population genetic structure.

Multilocus microsatellite typing of Leishmania infantum isolates in monitored Leishmania/HIV coinfected patients

Background

Leishmania infantum is the main etiological agent of both visceral and cutaneous clinical forms of leishmaniasis in the Mediterranean area. Leishmania/HIV coinfection in this area is characterized by a chronic course and frequent recurrences of clinical episodes. The present study using Multilocus Microsatellite Typing (MLMT) analysis, a highly discriminative tool, aimed to genetically characterize L. infantum isolates taken from monitored Leishmania/HIV coinfected patients presenting successive clinical episodes.

Methods

In this study, by the analysis of 20 microsatellite loci, we studied the MLMT profiles of 25 L. infantum isolates from 8 Leishmania/HIV coinfected patients who had experienced several clinical episodes. Two to seven isolates per patient were taken before and after treatment, during clinical and non-clinical episodes, with time intervals of 6 days to 29 months. Genetic diversity, clustering and phenetic analyses were performed.

Results

MLMT enabled us to study the genetic characteristics of the 25 L. infantum isolates, differentiating 18 genotypes, corresponding to a genotypic diversity of 0.72. Fifteen genotypes were unique in the total sample set and only 3 were repeated, 2 of which were detected in different patients. Both clustering and phylogenetic analyses provided insights into the genetic links between the isolates; in five patients isolates showed clear genetic links: either the genotype was exactly the same or only slightly different. In contrast, the isolates of the other three patients were dispersed in different clusters and some could be the result of mixing between populations.

Conclusions

Our data indicated a great MLMT variability between isolates from coinfected patients and no predominant genotype was observed. Despite this, almost all clinical episodes could be interpreted as a relapse rather than a reinfection. The results showed that diverse factors like an intrapatient evolution over time or culture bias could influence the parasite population detected in the patient, making it difficult to differentiate between relapse and reinfection.

Persistence, isolation and diversification of a naturally fragmented species in local refugia: the case of Hydromantes strinatii

The study of the European plethodontid salamander Hydromantes strinatii using allozyme and mitochondrial markers showed a strong geographical genetic structure. This was likely the outcome of different evolutionary mechanisms leaving their signature despite the effects of the genetic drift due to the low population size typical of this species. Two highly divergent clades were identified in the eastern and central-western part of the range, with further geographic sub-structure. Nuclear and mitochondrial markers substantially recovered the same population groups but were conflicting in reconstructing their relationships. This apparent incongruence highlighted the action of different mechanisms such as secondary contacts and incomplete lineage sorting in originating the observed genetic variation. The troglophilic habit of this species provided the opportunity to show the importance of caves as local refugia in maintaining the genetic diversity through the persistence of local populations. Accordingly, high nucleotide and haplotype diversity, strong geographic genetic structuring and lack of expansion were evidenced. This signature was found in the populations from the Ligurian and Maritime Alps, in agreement with the complex orography and paleoclimatic history of this Mediterranean hotspot.

The microgeographical patterns of morphological and molecular variation of a mixed ploidy population in the species complex Actinidia chinensis

Polyploidy and hybridization are thought to have significant impacts on both the evolution and diversification of the genus Actinidia, but the structure and patterns of morphology and molecular diversity relating to ploidy variation of wild Actinidia plants remain much less understood. Here, we examine the distribution of morphological variation and ploidy levels along geographic and environmental variables of a large mixed-ploidy population of the A. chinensis species complex. We then characterize the extent of both genetic and epigenetic diversity and differentiation exhibited between individuals of different ploidy levels. Our results showed that while there are three ploidy levels in this population, hexaploids were constituted the majority (70.3%). Individuals with different ploidy levels were microgeographically structured in relation to elevation and extent of niche disturbance. The morphological characters examined revealed clear difference between diploids and hexaploids, however tetraploids exhibited intermediate forms. Both genetic and epigenetic diversity were high but the differentiation among cytotypes was weak, suggesting extensive gene flow and/or shared ancestral variation occurred in this population even across ploidy levels. Epigenetic variation was clearly correlated with changes in altitudes, a trend of continuous genetic variation and gradual increase of epigenomic heterogeneities of individuals was also observed. Our results show that complex interactions between the locally microgeographical environment, ploidy and gene flow impact A. chinensis genetic and epigenetic variation. We posit that an increase in ploidy does not broaden the species habitat range, but rather permits A. chinensis adaptation to specific niches.

Strong genetic structure corresponds to small-scale geographic breaks in the Australian alpine grasshopper Kosciuscola tristis

Background

Mountain landscapes are topographically complex, creating discontinuous ‘islands’ of alpine and sub-alpine habitat with a dynamic history. Changing climatic conditions drive their expansion and contraction, leaving signatures on the genetic structure of their flora and fauna. Australia’s high country covers a small, highly fragmented area. Although the area is thought to have experienced periods of relative continuity during Pleistocene glacial periods, small-scale studies suggest deep lineage divergence across low-elevation gaps. Using both DNA sequence data and microsatellite markers, we tested the hypothesis that genetic partitioning reflects observable geographic structuring across Australia’s mainland high country, in the widespread alpine grasshopper Kosciuscola tristis (Sjösted).

Results

We found broadly congruent patterns of regional structure between the DNA sequence and microsatellite datasets, corresponding to strong divergence among isolated mountain regions. Small and isolated mountains in the south of the range were particularly distinct, with well-supported divergence corresponding to climate cycles during the late Pliocene and Pleistocene. We found mixed support, however, for divergence among other mountain regions. Interestingly, within areas of largely contiguous alpine and sub-alpine habitat around Mt Kosciuszko, microsatellite data suggested significant population structure, accompanied by a strong signature of isolation-by-distance.

Conclusions

Consistent patterns of strong lineage divergence among different molecular datasets indicate genetic breaks between populations inhabiting geographically distinct mountain regions. Three primary phylogeographic groups were evident in the highly fragmented Victorian high country, while within-region structure detected with microsatellites may reflect more recent population isolation. Despite the small area of Australia’s alpine and sub-alpine habitats, their low topographic relief and lack of extensive glaciation, divergence among populations was on the same scale as that detected in much more extensive Northern hemisphere mountain systems. The processes driving divergence in the Australian mountains might therefore differ from their Northern hemisphere counterparts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0204-1) contains supplementary material, which is available to authorized users.

Contrasting patterns of genetic structuring in natural populations of Arabidopsis lyrata Subsp. petraea across different regions in northern Europe

Level and partitioning of genetic diversity is expected to vary between contrasting habitats, reflecting differences in strength of ecological and evolutionary processes. Therefore, it is necessary to consider processes acting on different time scales when trying to explain diversity patterns in different parts of species' distributions. To explore how historical and contemporary factors jointly may influence patterns of genetic diversity and population differentiation, we compared genetic composition in the perennial herb Arabidopsis lyrata ssp. petraea from the northernmost parts of its distribution range on Iceland to that previously documented in Scandinavia. Leaf tissue and soil were sampled from ten Icelandic populations of A. lyrata. Seedlings were grown from soil samples, and tissue from above-ground and seed bank individuals were genotyped with 21 microsatellite markers. Seed bank density in Icelandic populations was low but not significantly different from that observed in Norwegian populations. While within-population genetic diversity was relatively high on Iceland (H(E) = 0.35), among-population differentiation was low (F(ST) = 0.10) compared to Norwegian and Swedish populations. Population differentiation was positively associated with geographical distance in both Iceland and Scandinavia, but the strength of this relationship varied between regions. Although topography and a larger distribution range may explain the higher differentiation between mountainous Norwegian relative to lowland populations in Sweden, these factors cannot explain the lower differentiation in Icelandic compared to Swedish populations. We propose that low genetic differentiation among Icelandic populations is not caused by differences in connectivity, but is rather due to large historical effective population sizes. Thus, rather than contemporary processes, historical factors such as survival of Icelandic lineages in northern refugia during the last glacial period may have contributed to the observed pattern.

Identification of hemiclonal reproduction in three species of Hexagrammos marine reef fishes

Natural hybrids between the boreal species Hexagrammos octogrammus and two temperate species Hexagrammos agrammus and Hexagrammos otakii were observed frequently in southern Hokkaido, Japan. Previous studies revealed that H. octogrammus is a maternal ancestor of both hybrids; the hybrids are all fertile females and they frequently breed with paternal species. Although such rampant hybridization occurs, species boundaries have been maintained in the hybrid zone. Possible explanations for the absence of introgressions, despite the frequent backcrossing, might include clonal reproduction: parthenogenesis, gynogenesis and hybridogenesis. The natural hybrids produced haploid eggs that contained only the H. octogrammus genome (maternal ancestor) with discarded paternal genome and generated F1 -hybrid type offspring by fertilization with the haploid sperm of H. agrammus or H. otakii (paternal ancestor). This reproductive mode was found in an artificial backcross hybrid between the natural hybrid and a male of the paternal ancestor. These findings indicate that the natural hybrids adopt hybridogenesis with high possibility and produce successive generations through hybridogenesis by backcrossing with the paternal ancestor. These hybrids of Hexagrammos represent the first hybridogenetic system found from marine fishes that widely inhabit the North Pacific Ocean. In contrast with other hybridogenetic systems, these Hexagrammos hybrids coexist with all three ancestral species in the hybrid zone. The coexistence mechanism is also discussed.

Morphology and genetics reveal an intriguing pattern of differentiation at a very small geographic scale in a bird species, the forest thrush Turdus lherminieri

Mobile organisms are expected to show population differentiation only over fairly large geographical distances. However, there is growing evidence of discrepancy between dispersal potential and realized gene flow. Here we report an intriguing pattern of differentiation at a very small spatial scale in the forest thrush (Turdus lherminieri), a bird species endemic to the Lesser Antilles. Analysis of 331 individuals from 17 sampling sites distributed over three islands revealed a clear morphological and genetic differentiation between these islands isolated by 40-50 km. More surprisingly, we found that the phenotypic divergence between the two geographic zones of the island of Guadeloupe was associated with a very strong genetic differentiation (Fst from 0.073-0.153), making this pattern a remarkable case in birds given the very small spatial scale considered. Molecular data (mitochondrial control region sequences and microsatellite genotypes) suggest that this strong differentiation could have occurred in situ, although alternative hypotheses cannot be fully discarded. This study suggests that the ongoing habitat fragmentation, especially in tropical forests, may have a deeper impact than previously thought on avian populations.

Detecting hybridization between Iranian wild wolf (Canis lupus pallipes) and free-ranging domestic dog (Canis familiaris) by analysis of microsatellite markers

The genetic threat due to hybridization with free-ranging dogs is one major concern in wolf conservation. The identification of hybrids and extent of hybridization is important in the conservation and management of wolf populations. Genetic variation was analyzed at 15 unlinked loci in 28 dogs, 28 wolves, four known hybrids, two black wolves, and one dog with abnormal traits in Iran. Pritchard's model, multivariate ordination by principal component analysis and neighbor joining clustering were used for population clustering and individual assignment. Analysis of genetic variation showed that genetic variability is high in both wolf and dog populations in Iran. Values of H(E) in dog and wolf samples ranged from 0.75-0.92 and 0.77-0.92, respectively. The results of AMOVA showed that the two groups of dog and wolf were significantly different (F(ST) = 0.05 and R(ST) = 0.36; P < 0.001). In each of the three methods, wolf and dog samples were separated into two distinct clusters. Two dark wolves were assigned to the wolf cluster. Also these models detected D32 (dog with abnormal traits) and some other samples, which were assigned to more than one cluster and could be a hybrid. This study is the beginning of a genetic study in wolf populations in Iran, and our results reveal that as in other countries, hybridization between wolves and dogs is sporadic in Iran and can be a threat to wolf populations if human perturbations increase.

Parasitic genotypes appear to differ in leishmaniasis patients compared with asymptomatic related carriers

For numerous infectious diseases affecting humans, clinical manifestations range from asymptomatic forms to severe pathologies. The originality of this study was its focus on asymptomatic carriers of Leishmania infantum in southern France. The fundamental interest in these asymptomatic carriers is that they can be a reservoir of potentially pathogenic microorganisms. It remains to be established whether the parasitic genomes from asymptomatic carriers differ from those of patients. Multilocus microsatellite typing was used to investigate the genetic variation among 36 French strains of L. infantum. Nine Leishmania strains isolated from blood donors (asymptomatic carriers) were compared with 27 strains of L. infantum belonging to zymodemes, MON-1, -33 and -183. These strains were isolated from HIV positive or negative patients with visceral leishmaniasis, cutaneous leishmaniasis, from canine leishmaniasis or from phlebotomine sandflies. Multilocus microsatellite typing data generated using 33 loci were analyzed by a Bayesian model-based clustering algorithm and construction of a phylogenetic tree based on genetic distances. Both analyses structured the MON-1 sample into two main clusters. Furthermore, genetic analysis demonstrated that these nine asymptomatic carrier strains are divided into two clusters grouped with the MON-1 strains. One cluster with seven strains is related to, but different from, human symptomatic strains from the Alpes-Maritimes region whereas the other cluster has the two remaining strains together with canine leishmaniasis strains as well as one strain from a visceral leishmaniasis patient. Genetic diversity among asymptomatic carrier was very weak since the nine Leishmania strains belong to only two genotypes. Genetic differentiations were evidenced between asymptomatic carrier strains and non-asymptomatic carrier strains and especially between asymptomatic carrier and HIV+ populations, although these findings require confirmation with a larger sample size. We believe that our data explore for the first time, the genetic diversity among L. infantum from asymptomatic human carriers and reveal a weak polymorphism compared with Leishmania parasites isolated from human patients.

The impact of multiple biogeographic barriers and hybridization on species-level differentiation

PREMISE OF THE STUDY

The glacial cycles of the Quaternary did not impact Australia in the same way as Europe and North America. Here we investigate the history of population isolation, species differentiation, and hybridization in the southeastern Australian landscape, using five species of Lomatia (Proteaceae). We use a chloroplast DNA phylogeography to assess chloroplast haplotype (chlorotype) sharing among these species and whether species with shared distributions have been affected by shared biogeographic barriers. •

METHODS

We used six chloroplast DNA simple sequence repeats (cpSSR) across five species of Lomatia, sampled across their entire distributional range in southeastern Australia. Resulting size data were combined, presented as a network, and visualized on a map. Biogeographical barriers were tested using AMOVA. To explore hypotheses of chlorotype origin, we converted the network into a cladogram and reconciled with all possible species trees using parsimony-based tree mapping. •

KEY RESULTS

Some chlorotypes were shared across multiple species of Lomatia in the study, including between morphologically differentiated species. Chlorotypes were either widespread in distribution or geographically restricted to specific regions. Biogeographical structure was identified across the range of Lomatia. The most parsimonious reconciled tree incorporated horizontal transfer of chlorotypes. •

CONCLUSIONS

Lomatia shows evidence of both incomplete lineage sorting and extensive hybridization between co-occurring species. Although the species in the study appear to have responded to a number of biogeographic barriers to varying degrees, our findings identified the Hunter River Valley as the most important long-term biogeographic barrier for the genus in southeastern Australia.

Contemporary and historical evolutionary processes interact to shape patterns of within-lake phenotypic divergences in polyphenic pumpkinseed sunfish, Lepomis gibbosus

Historical and contemporary evolutionary processes can both contribute to patterns of phenotypic variation among populations of a species. Recent studies are revealing how interactions between historical and contemporary processes better explain observed patterns of phenotypic divergence than either process alone. Here, we investigate the roles of evolutionary history and adaptation to current environmental conditions in structuring phenotypic variation among polyphenic populations of sunfish inhabiting 12 postglacial lakes in eastern North America. The pumpkinseed sunfish polyphenism includes sympatric ecomorphs specialized for littoral or pelagic lake habitats. First, we use population genetic methods to test the evolutionary independence of within-lake phenotypic divergences of ecomorphs and to describe patterns of genetic structure among lake populations that clustered into three geographical groupings. We then used multivariate analysis of covariance (MANCOVA) to partition body shape variation (quantified with geometric morphometrics) among the effects of evolutionary history (reflecting phenotypic variation among genetic clusters), the shared phenotypic response of all populations to alternate habitats within lakes (reflecting adaptation to contemporary conditions), and unique phenotypic responses to habitats within lakes nested within genetic clusters. All effects had a significant influence on body form, but the effects of history and the interaction between history and contemporary habitat were larger than contemporary processes in structuring phenotypic variation. This highlights how divergence can be better understood against a known backdrop of evolutionary history.

Spatial genetic structure of the mountain pine beetle (Dendroctonus ponderosae) outbreak in western Canada: historical patterns and contemporary dispersal

Environmental change has a wide range of ecological consequences, including species extinction and range expansion. Many studies have shown that insect species respond rapidly to climatic change. A mountain pine beetle epidemic of record size in North America has led to unprecedented mortality of lodgepole pine, and a significant range expansion to the northeast of its historic range. Our goal was to determine the spatial genetic variation found among outbreak population from which genetic structure, and dispersal patterns may be inferred. Beetles from 49 sampling locations throughout the outbreak area in western Canada were analysed at 13 microsatellite loci. We found significant north-south population structure as evidenced by: (i) Bayesian-based analyses, (ii) north-south genetic relationships and diversity gradients; and (iii) a lack of isolation-by-distance in the northernmost cluster. The north-south structure is proposed to have arisen from the processes of postglacial colonization as well as recent climate-driven changes in population dynamics. Our data support the hypothesis of multiple sources of origin for the outbreak and point to the need for population specific information to improve our understanding and management of outbreaks. The recent range expansion across the Rocky Mountains into the jack/lodgepole hybrid and pure jack pine zones of northern Alberta is consistent with a northern British Columbia origin. We detected no loss of genetic variability in these populations, indicating that the evolutionary potential of mountain pine beetle to adapt has not been reduced by founder events. This study illustrates a rapid range-wide response to the removal of climatic constraints, and the potential for range expansion of a regional population.

Polyploid evolution and Pleistocene glacial cycles: A case study from the alpine primrose Primula marginata (Primulaceae)

Background

Recent studies highlighted the role of Pleistocene climatic cycles in polyploid speciation and of southern Alpine refugia as reservoirs of diversity during glacial maxima. The polyploid Primula marginata, endemic to the southwestern Alps, includes both hexaploid and dodecaploid cytotypes that show no ecological or morphological differences. We used flow cytometry to determine variation and geographic distribution of cytotypes within and between populations and analyses of chloroplast (cp) and nuclear ribosomal (nr) DNA sequences from the Internal Transcribed Spacer (ITS) region to infer the evolutionary history of the two cytotypes and the auto- vs. allopolyploid origin of dodecaploid populations.

Results

We did not detect any intermediate cytotypes or variation of ploidy levels within populations. Hexaploids occur in the western and dodecaploids in the eastern part of the distributional range, respectively. The cpDNA and nrDNA topologies are in conflict, for the former supports shared ancestry between P. marginata and P. latifolia, while the latter implies common origins between at least some ITS clones of P. marginata and P. allionii.

Conclusions

Our results suggest an initial episode of chloroplast capture involving ancestral lineages of P. latifolia and P. marginata, followed by polyploidization between P. marginata-like and P. allionii-like lineages in a southern refugium of the Maritime Alps. The higher proportion of ITS polymorphisms in dodecaploid than in hexaploid accessions of P. marginata and higher total nucleotide diversity of ITS clones in dodecaploid vs. hexaploid individuals sequences are congruent with the allopolyploid hypothesis of dodecaploid origin.

Comparing the genetic structure of codling moth Cydia pomonella (L.) from Greece and France: long distance gene-flow in a sedentary pest species

Codling moth Cydia pomonella L. (Lepidoptera: Tortricidae) is the most important insect pest of apple production in Europe. Despite the economic importance of this pest, there is not information about the genetic structure of its population in Greece and the patterns of gene-flow which might affect the success of control programs. In this study, we analysed nine samples from apple, pear and walnut from various regions of mainland Greece using 11 microsatellite loci. Six samples from the aforementioned hosts from southern France were also examined for comparison. Bayesian clustering and genetic distance analyses separated the codling moth samples in two genetic clusters. The first cluster consisted mainly of the individuals from Greece, and the second of those from France, although admixture and miss-classified individuals were also observed. The low genetic differentiation among samples within each country was also revealed by F(ST) statistics (0.009 among Greek samples and 0.0150 among French samples compared to 0.050 global value among all samples and 0.032 the mean of the pair-wise values between the two countries). These F(ST) values suggest little structuring at large geographical scales in agreement with previous published studies. The host species and local factors (climatic conditions, topography, pest control programs) did not affect the genetic structure of codling moth populations within each country. The results are discussed in relation to human-made activities that promote gene-flow even at large geographic distances. Possible factors for the genetic differentiation between the two genetic clusters are also discussed.

From global to local genetic structuring in the red gorgonian Paramuricea clavata: the interplay between oceanographic conditions and limited larval dispersal

Defining the scale of connectivity among marine populations and identifying the barriers to gene flow are tasks of fundamental importance for understanding the genetic structure of populations and for the design of marine reserves. Here, we investigated the population genetic structure at three spatial scales of the red gorgonian Paramuricea clavata (Cnidaria, Octocorallia), a key species dwelling in the coralligenous assemblages of the Mediterranean Sea. Colonies of P. clavata were collected from 39 locations across the Mediterranean Sea from Morocco to Turkey and analysed using microsatellite loci. Within three regions (Medes, Marseille and North Corsica), sampling was obtained from multiple locations and at different depths. Three different approaches (measures of genetic differentiation, Bayesian clustering and spatially explicit maximum-difference algorithm) were used to determine the pattern of genetic structure. We identified genetic breaks in the spatial distribution of genetic diversity, which were concordant with oceanographic conditions in the Mediterranean Sea. We revealed a high level of genetic differentiation among populations and a pattern of isolation by distance across the studied area and within the three regions, underlining short effective larval dispersal in this species. We observed genetic differentiation among populations in the same locality dwelling at different depths, which may be explained by local oceanographic conditions and which may allow a process of local adaptation of the populations to their environment. We discuss the implications of our results for the conservation of the species, which is exposed to various threats.

Human impact on genetic diversity of Toxoplasma gondii: example of the anthropized environment from French Guiana

In French Guiana, severe cases of toxoplasmosis in immunocompetent patients are associated with atypical strains of Toxoplasma gondii linked to a wild neotropical rainforest cycle and a higher genetic diversity than usually observed for T. gondii isolates from anthropized environment. This raises the question of the impact of anthropization of the natural environment, on genetic diversity and on the population structure of T. gondii. However, few data are available on strains circulating in the anthropized areas from French Guiana. Seropositive animals originating mainly from anthropized sub-urban areas and punctually from wild environment in French Guiana were analyzed for T. gondii isolation and genotyping. Thirty-three strains were obtained by bioassay in mice and compared with 18 previously reported isolates chiefly originating from the Amazon rainforest. The genotyping analysis performed with 15 microsatellite markers located on 12 different chromosomes revealed a lower genetic diversity in the anthropized environment. Results were analyzed in terms of population structure by clustering methods, Neighbor-joining trees reconstruction based on genetic distances, F(ST,) Mantel's tests and linkage disequilibrium. They clearly showed a genetic differentiation between strains associated to the anthropized environment and those associated to the wild, but with some inbreeding between them. The majority of strains from the anthropized environment were clustered into additional lineages of T. gondii that are common in the Caribbean. In conclusion the two environmental populations "wild" and "anthropized" were genetically well differentiated. The anthropization of the environment seems to be accompanied with a decreased diversity of T. gondii associated with a greater structure of the populations. We detected potential interpenetration and genetic exchanges between these two environmental populations. As a higher pathogenicity in human of "wild" genotypes has been described, the interpenetration of both environments leads to hybridization between strains that may be at risk for human health.