Genetic differentiation and reduced genetic diversity at the northern range edge of two species with different dispersal modes

Genotype of pioneer plant Miscanthus is not a key factor in the structure of rhizosphere bacterial community in heavy metal polluted sites

Miscanthus is a common pioneer plant with abundant genetic variation in abandoned mines in southern China. However, the extent to which genetic differentiation among species modulates rhizosphere bacterial communities remains unclear. Miscanthus samples were collected from 26 typical abandoned heavy-metal mines with different soil types in southern China, tested using 14 pairs of simple sequence repeats (SSR) primers, and classified into two genotypes based on Nei's genetic distance. The structure and diversity of rhizosphere bacterial communities were examined using 16 S rRNA sequencing. The results showed that among the factors affecting the rhizosphere bacterial community structure of Miscanthus samples, the role of genotype was not significant, and geographical conditions were the most important factors, followed by pH and total organic carbon (TOC). The process of rhizospheric community assembly varied among different genotypes; however, the recruited species and their abundances were similar. Collectively, we provided an approach based on genetic differentiation to quantify the relative contribution of genotypes to the rhizosphere bacterial community, demonstrating that genotypes contribute less than soil conditions. Our findings provide new insights into the role of host genetics in the ecological processes of plant rhizosphere bacterial communities in abandoned mines and provide theoretical support for microbe-assisted phytoremediation.

Multiple genetic sources facilitate the northward range expansion of an intertidal oyster along China's coast

Coastal artificial structures on the former mudflats provide available habitats for the rocky intertidal species which can establish new populations in these emerging habitats over their former distribution range limits. As a former southern species, the oyster Crassostrea sikamea has become a pioneer and rapidly invaded the artificial shorelines in northern China. We used a seascape genomics approach to investigate the population structure and genetic sources of C. sikamea on the coastal artificial structures, which is crucial for understanding the genetic mechanisms driving species distribution range expansion and invasion pathway of intertidal species. Five C. sikamea populations, including two artificial substrate populations (WGZ and ZAP), one oyster reef population (LS), and two natural rocky shore populations (ZS and XM), were measured using single nucleotide polymorphism (SNPs) obtained from double digest restriction-site associated DNA sequencing (ddRAD-Seq). Redundancy analyses (RDA) were implemented for investigating the relationship between local temperature variables and the temperature adaptability of C. sikamea. Genetic diversity, direction and strength of gene flow, and population structure all revealed that the LS and ZS populations were the genetic sources for the oyster populations on the emerging northern coastal artificial structures. Results of RDA showed that there were different adaptive potentials for northern and southern populations to local temperature variables and the oyster reef population which frequently suffers from heat stress owned high heat adaptability. The ZS population as a genetic source nearby the Yangtze River estuary provided mass larvae for the northern populations, and the other genetic source, the heat-tolerant LS population, in the oyster reef played an important role in the post-settlement success by providing preadapted genotypes. These results highlight the importance of multiple sources with divergent adaptative capabilities for biological invasion, and also emphasize the importance of the oyster reef in coastal biodiversity and conservation.

Population transcriptomics reveals the effect of gene flow on the evolution of range limits

One of the most important questions in evolutionary biology is how the spatial distribution of species is limited. Asymmetric gene flow from core populations is suggested to increase the number of poorly adapted immigrants in the populations at the range edge. Genetic load due to migration, i.e., migration load, should prevent adaptation to the local habitat, leading to decreases in distribution range via local extinction or the limiting range expansion. However, few experimental studies have examined the effects of immigration on fitness and natural selection within recipient populations. To investigate the influence of migration load on the evolution of distribution range, we performed field and laboratory observations as well as population transcriptomics for the common river snail, Semisulcospira reiniana. This species meets the conditions that migration from source populations can prevent local adaptation in a sink population because they inhabit the broader range of environments, including middle/upper reaches of a river and estuaries within a single river and they may be more vulnerable to being swept away by water currents due to lowered spontaneous (upward) locomotion activity. We found that river steepness was related to the lower distribution limit of S. reiniana, with a narrower distribution range in the steeper river. Population transcriptomic analysis showed that gene flow was heavily asymmetric from the upstream populations to downstream ones in the steep river, suggesting a greater migration load in the steep river. The number of genes putatively involved in adaptation to the local habitat was lower in the steep river than in the gentle river. Gene expression profiles suggested that individuals achieve better local adaptation in the gentle river. Laboratory experiments suggested that evolutionary differences in salinity tolerance among local populations were only found in the gentle river. Our results consistent with the hypothesis that migration load owing to asymmetric gene flow disturbs local adaptation and restricts the distribution range of river snails.

Life history, climate and biogeography interactively affect worldwide genetic diversity of plant and animal populations

Understanding how biological and environmental factors interactively shape the global distribution of plant and animal genetic diversity is fundamental to biodiversity conservation. Genetic diversity measured in local populations (GDP) is correspondingly assumed representative for population fitness and eco-evolutionary dynamics. For 8356 populations across the globe, we report that plants systematically display much lower GDP than animals, and that life history traits shape GDP patterns both directly (animal longevity and size), and indirectly by mediating core-periphery patterns (animal fecundity and plant dispersal). Particularly in some plant groups, peripheral populations can sustain similar GDP as core populations, emphasizing their potential conservation value. We further find surprisingly weak support for general latitudinal GDP trends. Finally, contemporary rather than past climate contributes to the spatial distribution of GDP, suggesting that contemporary environmental changes affect global patterns of GDP. Our findings generate new perspectives for the conservation of genetic resources at worldwide and taxonomic-wide scales.

Chaotic genetic structure and past demographic expansion of the invasive gastropod Tritia neritea in its native range, the Mediterranean Sea

To better predict population evolution of invasive species in introduced areas it is critical to identify and understand the mechanisms driving genetic diversity and structure in their native range. Here, we combined analyses of the mitochondrial COI gene and 11 microsatellite markers to investigate both past demographic history and contemporaneous genetic structure in the native area of the gastropod Tritia neritea, using Bayesian skyline plots (BSP), multivariate analyses and Bayesian clustering. The BSP framework revealed population expansions, dated after the last glacial maximum. The haplotype network revealed a strong geographic clustering. Multivariate analyses and Bayesian clustering highlighted the strong genetic structure at all scales, between the Black Sea and the Adriatic Sea, but also within basins. Within basins, a random pattern of genetic patchiness was observed, suggesting a superimposition of processes involving natural biological effects (no larval phase and thus limited larval dispersal) and putative anthropogenic transport of specimens. Contrary to the introduced area, no isolation-by-distance patterns were recovered in the Mediterranean or the Black Seas, highlighting different mechanisms at play on both native and introduced areas, triggering unknown consequences for species’ evolutionary trajectories. These results of Tritia neritea populations on its native range highlight a mixture of ancient and recent processes, with the effects of paleoclimates and life history traits likely tangled with the effects of human-mediated dispersal.

Genetic diversity, population structure, and effective population size in two yellow bat species in south Texas

There are increasing concerns regarding bat mortality at wind energy facilities, especially as installed capacity continues to grow. In North America, wind energy development has recently expanded into the Lower Rio Grande Valley in south Texas where bat species had not previously been exposed to wind turbines. Our study sought to characterize genetic diversity, population structure, and effective population size in Dasypterus ega and D. intermedius, two tree-roosting yellow bats native to this region and for which little is known about their population biology and seasonal movements. There was no evidence of population substructure in either species. Genetic diversity at mitochondrial and microsatellite loci was lower in these yellow bat taxa than in previously studied migratory tree bat species in North America, which may be due to the non-migratory nature of these species at our study site, the fact that our study site is located at a geographic range end for both taxa, and possibly weak ascertainment bias at microsatellite loci. Historical effective population size (NEF) was large for both species, while current estimates of Ne had upper 95% confidence limits that encompassed infinity. We found evidence of strong mitochondrial differentiation between the two putative subspecies of D. intermedius (D. i. floridanus and D. i. intermedius) which are sympatric in this region of Texas, yet little differentiation using microsatellite loci. We suggest this pattern is due to secondary contact and hybridization and possibly incomplete lineage sorting at microsatellite loci. We also found evidence of some hybridization between D. ega and D. intermedius in this region of Texas. We recommend that our data serve as a starting point for the long-term genetic monitoring of these species in order to better understand the impacts of wind-related mortality on these populations over time.

Signatures of local adaptation in the spatial genetic structure of the ascidian Pyura chilensis along the southeast Pacific coast

The highly heterogeneous Humboldt Current System (HCS) and the 30°S transition zone on the southeast Pacific coast, represent an ideal scenario to test the influence of the environment on the spatial genomic structure in marine near-shore benthic organisms. In this study, we used seascape genomic tools to evaluate the genetic structure of the commercially important ascidian Pyura chilensis, a species that exhibits a low larval transport potential but high anthropogenic dispersal. A recent study in this species recorded significant genetic differentiation across a transition zone around 30°S in putatively adaptive SNPs, but not in neutral ones, suggesting an important role of environmental heterogeneity in driving genetic structure. Here, we aim to understand genomic-oceanographic associations in P. chilensis along the Southeastern Pacific coast using two combined seascape genomic approaches. Using 149 individuals from five locations along the HCS, a total of 2,902 SNPs were obtained by Genotyping-By-Sequencing, of which 29–585 were putatively adaptive loci, depending on the method used for detection. In adaptive loci, spatial genetic structure was better correlated with environmental differences along the study area (mainly to Sea Surface Temperature, upwelling-associated variables and productivity) than to the geographic distance between sites. Additionally, results consistently showed the presence of two groups, located north and south of 30°S, which suggest that local adaptation processes seem to allow the maintenance of genomic differentiation and the spatial genomic structure of the species across the 30°S biogeographic transition zone of the Humboldt Current System, overriding the homogenizing effects of gene flow.

Rejection of the genetic implications of the "Abundant Centre Hypothesis" in marine mussels

The ‘Abundant-Centre Hypothesis’ is a well-established but controversial hypothesis stating that the abundance of a species is highest at the centre of its range and decreases towards the edges, where conditions are unfavourable. As genetic diversity depends on population size, edge populations are expected to show lower intra-population genetic diversity than core populations, while showing high inter-population genetic divergence. Here, the genetic implications of the Abundant-Centre Hypothesis were tested on two coastal mussels from South Africa that disperse by means of planktonic larvae, the native Perna perna and the invasive Mytilus galloprovincialis. Genetic structure was found within P. perna, which, together with evidence from Lagrangian particle simulations, points to significant reductions in gene flow between sites. Despite this, the expected diversity pattern between centre and edge populations was not found for either species. We conclude that the genetic predictions of the Abundant-Centre Hypothesis are unlikely to be met by high-dispersal species with large population sizes, and may only become evident in species with much lower levels of connectivity.

Chloroplast population genetics reveals low levels of genetic variation and conformation to the central-marginal hypothesis in Taxus wallichiana var. mairei, an endangered conifer endemic to China

The central-marginal hypothesis predicts that geographically peripheral populations should exhibit reduced genetic diversity and increased genetic differentiation than central populations due to smaller effective population size and stronger geographical isolation. We evaluated these predictions in the endangered conifer Taxus wallichiana var. mairei. Eight plastid simple sequence repeats (cpSSRs) were used to investigate plastid genetic variation in 22 populations of Taxus wallichiana var. mairei, encompassing nearly its entire distribution range. Low levels of plastid genetic variation and differentiation were detected in the populations, and the findings were attributed to low mutation rates, small population sizes, habitat fragmentation and isolation, and effective pollen or seed dispersal. Hunan and Hubei were identified as major refugia based on the number of private haplotypes and species distribution modeling. Trends in plastid genetic diversity and genetic differentiation from central to peripheral populations supported the predictions of the central-marginal hypothesis. In scenarios wherein the future climate becomes warmer, we predict that some peripheral populations will disappear and southern and southeastern regions will become significantly less habitable. Factors that include the levels of precipitation during the driest month, annual precipitation level, and annual temperature range will be decisive in shaping the future distribution of these populations. This study provides a theoretical basis for the conservation of T. wallichiana var. mairei.

A Review of the Effects of Climate Change on Chelonians

Climate change is occurring at an unprecedented rate and has begun to modify the distribution and phenology of organisms worldwide. Chelonians are expected to be particularly vulnerable due to limited dispersal capabilities as well as widespread temperature-dependent sex determination. The number of papers published about the effects of climate change on turtles has increased exponentially over the last decade; here, I review the data from peer-reviewed publications to assess the likely impacts of climate change on individuals, populations, and communities. Based upon these studies future research should focus on: (1) Individual responses to climate change, particularly with respect to thermal biology, phenology, and microhabitat selection; (2) improving species distribution models by incorporating fine-scale environmental variables as well as physiological processes; (3) identifying the consequences of skewed sex ratios; and (4) assessments of community resilience and the development of methods to mitigate climate change impacts. Although detailed management recommendations are not possible at this point, careful consideration should be given regarding how to manage low vagility species as habitats shift poleward. In the worst-case scenario, proactive management may be required in order to ensure that widespread losses do not occur.

Antibody Cross-Reactivity between Proteins of Chia Seed ( Salvia hispanica L.) and Other Food Allergens

Chia seeds are becoming increasingly common in Europe because of their functional and nutritional properties. Despite this, few studies have focused on the allergic potential and antibody cross-reactivity among storage proteins in chia seed and other plants. The aim of this study was to identify chia seed's immunoglobulin G (IgG) and immunoglobulin E (IgE) binding proteins ( Salvia hispanica L.) and to investigate the antibody cross-reactivity among its storage proteins and those of other seeds. Extracted chia seed proteins were separated by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). Immunodetection was performed with commercial antibodies against sesame seed, hazelnut, and peanut and sera from 33 patients with a hazelnut allergy and five with a sesame allergy. Cross-reactivity of certain antibodies with storage proteins of chia seed, sesame seed, and hazelnut was assessed using an enzyme-linked immunosorbent assay (ELISA) inhibition, blot inhibition, and surface plasmon resonance (SPR) spectroscopy. IgG binding proteins were identified at molecular weight (MW) 70, 49, 34, 23, and 20 kDa by applying commercial antibodies. Furthermore, the interaction of chia proteins with sera from sesame-allergic patients led to identify IgE binding proteins at MW 49, 45, 31, 20, and 12 kDa, while IgEs in sera from hazelnut-allergic patients reacted with proteins at MW 300, 140, 49, 45, 31, 20, and 6 kDa. The results of ELISA inhibition and blot inhibition indicated chia seed proteins are similar to sesame seed and hazelnut proteins in the primary structure. The antisesame antibodies' binding to sesame proteins was more strongly inhibited by the chia globulin fraction (GLO) than the antihazelnut antibodies' binding to hazelnut proteins. SPR results confirmed the presence of IgG binding proteins in GLO and the high similarity of epitopes on globulins of chia seed and sesame seed. Thus, chia seed consumption might lead to cross-sensitization in patients with a sesame allergy.

Development of EST-SSR markers and their application in an analysis of the genetic diversity of the endangered species Magnolia sinostellata

Magnolia sinostellata is an endemic species of Magnoliaceae that is narrowly distributed in the south of Zhejiang Province, China. To explore the genetic diversity and population structure of this endangered species, this study developed sequence tag-simple sequence repeat (EST-SSR) markers based on transcriptome data of M. sinostellata. In total, 25472 SSRs were identified among 110644 unique assembled sequences with a total of 90.83 Mb and an average frequency of 23.02%. The mononucleotide (33.53%) and dinucleotide (42.08%) motifs appeared to be the most abundant. In total, 150 potential loci were randomly selected to validate the quality of the developed SSR markers; an effective PCR rate of 32.00% and a polymorphism rate of 15.33% were obtained for these loci. After performing sequencing and cloning for validation, 23 pairs of SSR primers were retained and used to characterize the genetic diversity and population structure of M. sinostellata. Overall, 204 alleles were amplified. The results of Shannon's information index (I), heterozygosity (Ho), heterozygosity (He) and Nei's expected heterozygosity (H) indicated rich genetic diversity in M. sinostellata. However, the high inbreeding coefficient and differential coefficient suggest that serious genetic drift occurred within populations, and genetic differentiation is apparent among the populations. Consequently, although M. sinostellata has high genetic diversity among populations, it is still in a serious and dangerous condition. Habitat destruction caused by human activities is the main threat to this species, and enhancing the species abundance by adopting some conservation measures should be favourable for saving the species.

Discovering variation of secondary metabolite diversity and its relationship with disease resistance in Cornus florida L

Understanding intraspecific relationships between genetic and functional diversity is a major goal in the field of evolutionary biology and is important for conserving biodiversity. Linking intraspecific molecular patterns of plants to ecological pressures and trait variation remains difficult due to environment-driven plasticity. Next-generation sequencing, untargeted liquid chromatography-mass spectrometry (LC-MS) profiling, and interdisciplinary approaches integrating population genomics, metabolomics, and community ecology permit novel strategies to tackle this problem. We analyzed six natural populations of the disease-threatened Cornus florida L. from distinct ecological regions using genotype-by-sequencing markers and LC-MS-based untargeted metabolite profiling. We tested the hypothesis that higher genetic diversity in C. florida yielded higher chemical diversity and less disease susceptibility (screening hypothesis), and we also determined whether genetically similar subpopulations were similar in chemical composition. Most importantly, we identified metabolites that were associated with candidate loci or were predictive biomarkers of healthy or diseased plants after controlling for environment. Subpopulation clustering patterns based on genetic or chemical distances were largely congruent. While differences in genetic diversity were small among subpopulations, we did observe notable similarities in patterns between subpopulation averages of rarefied-allelic and chemical richness. More specifically, we found that the most abundant compound of a correlated group of putative terpenoid glycosides and derivatives was correlated with tree health when considering chemodiversity. Random forest biomarker and genomewide association tests suggested that this putative iridoid glucoside and other closely associated chemical features were correlated to SNPs under selection.

Is the sky the limit? On the expansion threshold of a species' range

More than 100 years after Grigg's influential analysis of species' borders, the causes of limits to species' ranges still represent a puzzle that has never been understood with clarity. The topic has become especially important recently as many scientists have become interested in the potential for species' ranges to shift in response to climate change-and yet nearly all of those studies fail to recognise or incorporate evolutionary genetics in a way that relates to theoretical developments. I show that range margins can be understood based on just two measurable parameters: (i) the fitness cost of dispersal-a measure of environmental heterogeneity-and (ii) the strength of genetic drift, which reduces genetic diversity. Together, these two parameters define an 'expansion threshold': adaptation fails when genetic drift reduces genetic diversity below that required for adaptation to a heterogeneous environment. When the key parameters drop below this expansion threshold locally, a sharp range margin forms. When they drop below this threshold throughout the species' range, adaptation collapses everywhere, resulting in either extinction or formation of a fragmented metapopulation. Because the effects of dispersal differ fundamentally with dimension, the second parameter-the strength of genetic drift-is qualitatively different compared to a linear habitat. In two-dimensional habitats, genetic drift becomes effectively independent of selection. It decreases with 'neighbourhood size'-the number of individuals accessible by dispersal within one generation. Moreover, in contrast to earlier predictions, which neglected evolution of genetic variance and/or stochasticity in two dimensions, dispersal into small marginal populations aids adaptation. This is because the reduction of both genetic and demographic stochasticity has a stronger effect than the cost of dispersal through increased maladaptation. The expansion threshold thus provides a novel, theoretically justified, and testable prediction for formation of the range margin and collapse of the species' range.

A multispecies approach reveals hot spots and cold spots of diversity and connectivity in invertebrate species with contrasting dispersal modes

Genetic diversity is crucial for species' maintenance and persistence, yet is often overlooked in conservation studies. Species diversity is more often reported due to practical constraints, but it is unknown if these measures of diversity are correlated. In marine invertebrates, adults are often sessile or sedentary and populations exchange genes via dispersal of gametes and larvae. Species with a larval period are expected to have more connected populations than those without larval dispersal. We assessed the relationship between measures of species and genetic diversity, and between dispersal ability and connectivity. We compiled data on genetic patterns and life history traits in nine species across five phyla. Sampling sites spanned 600 km in the northwest Mediterranean Sea and focused on a 50-km area near Marseilles, France. Comparative population genetic approaches yielded three main results. (i) Species without larvae showed higher levels of genetic structure than species with free-living larvae, but the role of larval type (lecithotrophic or planktotrophic) was negligible. (ii) A narrow area around Marseilles, subject to offshore advection, limited genetic connectivity in most species. (iii) We identified sites with significant positive contributions to overall genetic diversity across all species, corresponding with areas near low human population densities. In contrast, high levels of human activity corresponded with a negative contribution to overall genetic diversity. Genetic diversity within species was positively and significantly linearly related to local species diversity. Our study suggests that local contribution to overall genetic diversity should be taken into account for future conservation strategies.

Genome-wide survey of single-nucleotide polymorphisms reveals fine-scale population structure and signs of selection in the threatened Caribbean elkhorn coral, Acropora palmata

The advent of next-generation sequencing tools has made it possible to conduct fine-scale surveys of population differentiation and genome-wide scans for signatures of selection in non-model organisms. Such surveys are of particular importance in sharply declining coral species, since knowledge of population boundaries and signs of local adaptation can inform restoration and conservation efforts. Here, we use genome-wide surveys of single-nucleotide polymorphisms in the threatened Caribbean elkhorn coral, Acropora palmata, to reveal fine-scale population structure and infer the major barrier to gene flow that separates the eastern and western Caribbean populations between the Bahamas and Puerto Rico. The exact location of this break had been subject to discussion because two previous studies based on microsatellite data had come to differing conclusions. We investigate this contradiction by analyzing an extended set of 11 microsatellite markers including the five previously employed and discovered that one of the original microsatellite loci is apparently under selection. Exclusion of this locus reconciles the results from the SNP and the microsatellite datasets. Scans for outlier loci in the SNP data detected 13 candidate loci under positive selection, however there was no correlation between available environmental parameters and genetic distance. Together, these results suggest that reef restoration efforts should use local sources and utilize existing functional variation among geographic regions in ex situ crossing experiments to improve stress resistance of this species.

Mechanisms of peripheral phylogeographic divergence in the indo-Pacific: lessons from the spiny lobster Panulirus homarus

Background

There is increasing recognition of the concordance between marine biogeographic and phylogeographic boundaries. However, it is still unclear how population-level divergence translates into species-level divergence, and what are the principal factors that first initiate that divergence, and then maintain reproductive isolation. This study examines the likely forces driving population and lineage divergences in the broadly-distributed Indo-Pacific spiny lobster Panulirus homarus, which has peripheral divergent lineages in the west and east. The study focuses particularly on the West Indian Ocean, which is emerging as a region of unexpected diversity. Mitochondrial control region (mtCR) and COI sequences as well as genotypes of 9 microsatellite loci were examined in 410 individuals from 17 locations grouped into 7 regions from South Africa in the west, and eastward across to Taiwan and the Marquesas Islands. Phylogenetic and population-level analyses were used to test the significance and timing of divergences and describe the genetic relationships among populations.

Results

Analyses of the mtCR revealed high levels of divergence among the seven regions (Ф_ST = 0.594, P < 0.001). Microsatellite analyses also revealed significant divergence among regions, but at a much lower level (F_ST = 0.066, P < 0.001). The results reveal different patterns of mtCR v. nDNA divergence between the two distinct peripheral lineages: a subspecies in South Africa and Madagascar, and a phylogeographically diverged population in the Marquesas. The results also expose a number of other more fine-scale population divergences, particularly in the Indian Ocean.

Conclusions

The divergence of peripheral lineages in the west and east of the species’ range appear to have been initiated and maintained by very different processes. The pattern of mitochondrial and nuclear divergence of the western lineage, implicates processes of parapatric isolation, secondary contact and introgression, and suggests possible maintenance through adaptation and behavioural reproductive isolation. In contrast, the eastern lineage appears to have diverged through a rare colonisation event, maintained through long-term isolation, and matches expectations of the core-periphery hypothesis. The process of active peripheral speciation may be a common force in the Indo-Pacific that helps drive some of the regions’ recognized biogeographic boundaries.

Electronic supplementary material

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

Lineage divergence, local adaptation across a biogeographic break, and artificial transport, shape the genetic structure in the ascidian Pyura chilensis

Marine benthic organisms inhabit a heterogeneous environment in which connectivity between populations occurs mainly through dispersive larval stages, while local selective pressures acting on early life history stages lead to non-random mortality, shaping adaptive genetic structure. In order to test the influence of local adaptation and neutral processes in a marine benthic species with low dispersal, in this study we used Genotyping by Sequencing technology to compare the neutral and putatively selected signals (neutral and outlier loci, respectively) in SNPs scattered throughout the genome in six local populations of the commercially exploited ascidian Pyura chilensis along the southeast Pacific coast (24°–42°S). This species is sessile as an adult, has a short-lived larval stage, and may also be dispersed by artificial transport as biofouling. We found that the main signal in neutral loci was a highly divergent lineage present at 39°S, and a subjacent signal that indicated a separation at 30°S (north/south), widely reported in the area. North/south separation was the main signal in outlier loci, and the linage divergence at 39°S was subjacent. We conclude that the geographic structure of the genetic diversity of outlier and neutral loci was established by different strengths of environmental, historical and anthropogenic factors.

Solar Radiation-Associated Adaptive SNP Genetic Differentiation in Wild Emmer Wheat, Triticum dicoccoides

Whole-genome scans with large number of genetic markers provide the opportunity to investigate local adaptation in natural populations and identify candidate genes under positive selection. In the present study, adaptation genetic differentiation associated with solar radiation was investigated using 695 polymorphic SNP markers in wild emmer wheat originated in a micro-site at Yehudiyya, Israel. The test involved two solar radiation niches: (1) sun, in-between trees; and (2) shade, under tree canopy, separated apart by a distance of 2-4 m. Analysis of molecular variance showed a small (0.53%) but significant portion of overall variation between the sun and shade micro-niches, indicating a non-ignorable genetic differentiation between sun and shade habitats. Fifty SNP markers showed a medium (0.05 ≤ FST ≤ 0.15) or high genetic differentiation (FST > 0.15). A total of 21 outlier loci under positive selection were identified by using four different FST -outlier testing algorithms. The markers and genome locations under positive selection are consistent with the known patterns of selection. These results suggested that genetic differentiation between sun and shade habitats is substantial, radiation-associated, and therefore ecologically determined. Hence, the results of this study reflected effects of natural selection through solar radiation on EST-related SNP genetic diversity, resulting presumably in different adaptive complexes at a micro-scale divergence. The present work highlights the evolutionary theory and application significance of solar radiation-driven natural selection in wheat improvement.