Exploring a Pool‐seq‐only approach for gaining population genomic insights in nonmodel species

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool‐seq data to generate a de novo genome assembly for mining exons, upon which Pool‐seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual‐based single nucleotide polymorphisms [SNPs]) and from mapping the Pool‐seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F_ST) between the two introduced populations exceeds that of the naturally sympatric populations (F_ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( ≈ 0.002 and  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high‐quality reference assembly from a divergent species. We conclude that the Pool‐seq‐only approach can be suitable for detecting and quantifying genome‐wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.     

Whole‐genome resequencing reveals the pleistocene temporal dynamics of Branchiostoma belcheri and Branchiostoma floridae populations

Global climatic fluctuations governed the ancestral demographic histories of species and contributed to place the current population status into a more extensive ecological and evolutionary context. Genetic variations will leave unambiguous signatures in the patterns of intraspecific genetic variation in extant species since the genome of each individual is an imperfect mosaic of the ancestral genomes. Here, we report the genome sequences of 20 Branchiostoma individuals by whole‐genome resequencing strategy. We detected over 140 million genomic variations for each Branchiostoma individual. In particular, we applied the pairwise sequentially Markovian coalescent (PSMC) method to estimate the trajectories of changes in the effective population size (N_e) of Branchiostoma population during the Pleistocene. We evaluated the threshold of sequencing depth for proper inference of demographic histories using PSMC was ≥25×. The PSMC results highlight the role of historical global climatic fluctuations in the long‐term population dynamics of Branchiostoma. The inferred ancestral N_e of the Branchiostoma belcheri populations from Zhanjiang and Xiamen (China) seawaters was different in amplitude before the first (mutation rate = 3 × 10^?9) or third glaciation (mutation rate = 9 × 10^?9) of the Pleistocene, indicating that the two populations most probably started to evolve in isolation in their respective seas after the first or third glaciation of the Pleistocene. A pronounced population bottleneck coinciding with the last glacial maximum was observed in all Branchiostoma individuals, followed by a population expansion occurred during the late Pleistocene. Species that have experienced long‐term declines may be especially vulnerable to recent anthropogenic activities. Recently, the industrial pollution and the exploitation of sea sand have destroyed the harmonious living environment of amphioxus species. In the future, we need to protect the habitat of Branchiostoma and make full use of these detected genetic variations to facilitate the functional study of Branchiostoma for adaptation to local environments.     

PSMC analysis of effective population sizes in molecular ecology and its application to black‐and‐white Ficedula flycatchers

Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have previously been hindered by a lack of means and genetic data to model changes in effective population size (N_e) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in N_e over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black‐and‐white flycatchers. N_e curves of Atlas, collared, pied and semicollared flycatcher converged 1–2 million years ago (Ma) at an N_e of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate N_e trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100–200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in N_e differed among species, and for pied flycatcher, the temporal dynamics of N_e differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per‐site filter of ≥10 reads and no more than 25% of missing data.     

Locus‐specific ancestry to detect recent response to selection in admixed Swiss Fleckvieh cattle

Identification of selection signatures is one of the current endeavors of evolutionary genetics. Admixed populations may be used to infer post‐admixture selection. We calculated local ancestry for Swiss Fleckvieh, a composite of Simmental (SI) and Red Holstein Friesian (RHF), to infer such signals. Illumina Bovine SNP50 BeadChip data for 300 admixed, 88 SI and 97 RHF bulls were used. The average RHF ancestry across the whole genome was 0.70. To identify regions with high deviation from average, we considered two significance thresholds, based on a permutation test and extreme deviation from normal distribution. Regions on chromosomes 13 (46.3–47.3 Mb) and 18 (18.7–25.9 Mb) passed both thresholds in the direction of increased SI. Extended haplotype homozygosity within (iHS) and between (Rsb) populations was calculated to explore additional patterns of pre‐ and post‐admixture selection signals. The Rsb score of admixed and SI was significant in a wide region of chromosome 18 (6.6–24.6 Mb) overlapped with one area of strong local ancestry deviation. FTO, with pleiotropic effect on milk and fertility, NOD2 on dairy and NKD1 and SALL1 on fertility traits are located there. Genetic differentiation of RHF and SI (F_st), an alternative indicator of pre‐admixture selection in pure populations, was calculated. No considerable overlap of peaks of local ancestry deviations and F_st was observed. We found two regions with significant signatures of post‐admixture selection in this very young composite, applying comparatively stringent significance thresholds. The signals cover relatively large genomic areas and did not allow pinpointing of the gene(s) responsible for the apparent shift in ancestry proportions.     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

Mapping the genomic architecture of ecological speciation in the wild: does linkage disequilibrium hold the key?

The hunt for the genes underlying ecological speciation has now closed in on a number of candidates, but making the link from genotype to phenotype continues to pose a significant challenge. This is partly because genetic studies in many systems remain impeded by long generation times or an inability to perform controlled crosses. Now, in this issue of Molecular Ecology, Malek et al. (2012) demonstrate the utility of a novel admixture mapping approach that can be used to identify genomic regions contributing to adaptive trait divergence between natural populations. Remarkably, they validate their approach by mapping traits associated with mate choice in a wild limnetic and benthic threespine stickleback (Gasterosteus aculeatus) species pair, finding several loci associated with male nuptial coloration and shape. While this study benefited from tried‐and‐true microsatellites in a well‐characterized species with a detailed genetic map (and genome sequence), the field is quickly moving towards the use of next‐generation sequencing, especially for nonmodel systems. The ability to characterize molecular polymorphisms for any system suggests that molecular ecologists working on virtually any species may benefit from applying Malek et al.'s approach, if naturally admixed populations are available.     

Do plant populations on distinct inselbergs talk to each other? A case study of genetic connectivity of a bromeliad species in an Ocbil landscape

Here, we explore the historical and contemporaneous patterns of connectivity among Encholirium horridum populations located on granitic inselbergs in an Ocbil landscape within the Brazilian Atlantic Forest, using both nuclear and chloroplast microsatellite markers. Beyond to assess the E. horridum population genetic structure, we built species distribution models across four periods (current conditions, mid‐Holocene, Last Glacial Maximum [LGM], and Last Interglacial) and inferred putative dispersal corridors using a least‐cost path analysis to elucidate biogeographic patterns. Overall, high and significant genetic divergence was estimated among populations for both nuclear and plastid DNA (Φ_ST(n) = 0.463 and Φ_ST(plastid) = 0.961, respectively, p < .001). For nuclear genome, almost total absence of genetic admixture among populations and very low migration rates were evident, corroborating with the very low estimates of immigration and emigration rates observed among E. horridum populations. Based on the cpDNA results, putative dispersal routes in Sugar Loaf Land across cycles of climatic fluctuations in the Quaternary period revealed that the populations’ connectivity changed little during those events. Genetic analyses highlighted the low genetic connectivity and long‐term persistence of populations, and the founder effect and genetic drift seemed to have been very important processes that shaped the current diversity and genetic structure observed in both genomes. The genetic singularity of each population clearly shows the need for in situ conservation of all of them.     

Demographic expansion and genetic load of the halophyte model plant Eutrema salsugineum

The halophyte model plant Eutrema salsugineum (Brassicaceae) disjunctly occurs in temperate to subarctic Asia and North America. This vast, yet extremely discontinuous distribution constitutes an ideal system to examine long‐distance dispersal and the ensuing accumulation of deleterious mutations as expected in expanding populations of selfing plants. In this study, we resequenced individuals from 23 populations across the range of E. salsugineum. Our population genomic data indicate that E. salsugineum migrated “out of the Altai region” at least three times to colonize northern China, northeast Russia and western China. It then expanded its distribution into North America independently from northeast Russia and northern China, respectively. The species colonized northern China around 33.7 thousand years ago (kya) and underwent a considerable expansion in range size approximately 7–8 kya. The western China lineage is likely a hybrid derivative of the northern China and Altai lineages, originating approximately 25–30 kya. Deleterious alleles accumulated in a stepwise manner from (a) Altai to northern China and North America and (b) Altai to northeast Russia and North America. In summary, E. salsugineum dispersed from Asia to North America and deleterious mutations accumulated in a stepwise manner during the expansion of the species’ distribution.     

The population demography of Betula maximowicziana,a cool‐temperate tree species in Japan,in relation to the last glacial period: its admixture‐like genetic structure is the result of simple population splitting not admixing

Conservation of the local genetic variation and evolutionary integrity of economically and ecologically important trees is a key aspect of studies involving forest genetics, and a population demographic history of the target species provides valuable information for this purpose. Here, the genetic structure of 48 populations of Betula maximowicziana was assessed using 12 expressed sequence tag–simple sequence repeat (EST‐SSR) markers. Genetic diversity was lower in northern populations than southern ones and structure analysis revealed three groups: northern and southern clusters and an admixed group. Eleven more genomic‐SSR loci were added and the demographic history of these three groups was inferred by approximate Bayesian computation (ABC). The ABC revealed that a simple split scenario was much more likely than isolation with admixture, suggesting that the admixture‐like structure detected in this species was due to ancestral polymorphisms. The ABC analysis suggested that the population growth and divergence of the three groups occurred 96 800 (95% CI, 20 500–599 000) and 28 300 (95% CI, 8700–98 400) years ago, respectively. We need to be aware of several sources of uncertainty in the inference such as assumptions about the generation time, overlapping of generations, confidence intervals of the estimated parameters and the assumed model in the ABC. However, the results of the ABC together with the model‐based maps of reconstructed past species distribution and palaeoecological data suggested that the modern genetic structure of B. maximowicziana originated prior to the last glacial maximum (LGM) and that some populations survived in the northern range even during the LGM.     

ABC inference of multi‐population divergence with admixture from unphased population genomic data

Rapidly developing sequencing technologies and declining costs have made it possible to collect genome‐scale data from population‐level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three‐population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation‐based cross‐validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.     

Tracking the elusive history of diversification in plant–herbivorous insect–parasitoid food webs: insights from figs and fig wasps

The food webs consisting of plants, herbivorous insects and their insect parasitoids are a major component of terrestrial biodiversity. They play a central role in the functioning of all terrestrial ecosystems, and the number of species involved is mind‐blowing (Nyman et al. 2015 ). Nevertheless, our understanding of the evolutionary and ecological determinants of their diversity is still in its infancy. In this issue of Molecular Ecology, Sutton et al. ( 2016 ) open a window into the comparative analysis of spatial genetic structuring in a set of comparable multitrophic models, involving highly species‐specific interactions: figs and fig wasps. This is the first study to compare genetic structure using population genetics tools in a fig‐pollinating wasp (Pleistodontes imperialis sp1) and its main parasitoid (Sycoscapter sp.A). The fig‐pollinating wasp has a discontinuous spatial distribution that correlates with genetic differentiation, while the parasitoid bridges the discontinuity by parasitizing other pollinator species on the same host fig tree and presents basically no spatial genetic structure. The full implications of these results for our general understanding of plant–herbivorous insect–insect parasitoids diversification become apparent when envisioned within the framework of recent advances in fig and fig wasp biology.     

Deep divergence among mitochondrial lineages in African jackals

Recently, molecular analyses revealed that African and Eurasian golden jackals are distinct species. This finding suggests re‐investigation of the phylogenetic relationships and taxonomy of other African members of the Canidae. Here, we provide a study on the phylogenetic relationship between populations of African jackals Lupulella mesomelas and L. adusta inferred from 962 bp of the mitochondrial cytochrome b (cytb) gene. As expected from its disjunct distribution, with one population in eastern Africa and the other one in southern Africa, we found two mitochondrial lineages within L. mesomelas, which diverged about 2.5 million years ago (Ma). In contrast, in L. adusta with its more continuous distribution in sub‐Saharan Africa, we found only a shallower genetic diversification, with the exception of the West African population, which diverged around 1.4 Ma from the Central and East African populations. Both divergence ages are older than, for example the 1.1–0.9 million years between the grey wolf Canis lupus and the African golden wolf C. lupaster. One taxonomic implication of our findings might be that the two L. mesomelas populations warrant species status. However, genome‐wide data with adequate geographical sampling are needed to substantiate our results.     

Demography and genome divergence of lake and stream populations of an East African cichlid fish

Disentangling the processes and mechanisms underlying adaptive diversification is facilitated by the comparative study of replicate population pairs that have diverged along a similar environmental gradient. Such a setting is realized in a cichlid fish from southern Lake Tanganyika, Astatotilapia burtoni, which occurs within the lake proper as well as in various affluent rivers. Previously, we demonstrated that independent lake and stream populations show similar adaptations to the two habitat regimes. However, little is known about the evolutionary and demographic history of the A. burtoni populations in question and the patterns of genome divergence among them. Here, we apply restriction site‐associated DNA sequencing (RADseq) to examine the evolutionary history, the population structure and genomic differentiation of lake and stream populations in A. burtoni. A phylogenetic reconstruction based on genome‐wide molecular data largely resolved the evolutionary relationships among populations, allowing us to re‐evaluate the independence of replicate lake–stream population clusters. Further, we detected a strong pattern of isolation by distance, with baseline genomic divergence increasing with geographic distance and decreasing with the level of gene flow between lake and stream populations. Genome divergence patterns were heterogeneous and inconsistent among lake‐stream population clusters, which is explained by differences in divergence times, levels of gene flow and local selection regimes. In line with the latter, we only detected consistent outlier loci when the most divergent lake–stream population pair was excluded. Several of the thus identified candidate genes have inferred functions in immune and neuronal systems and show differences in gene expression between lake and stream populations.     

Local adaptation through genetic differentiation in highly fragmented Tilia cordata populations

We assessed the level of geographic differentiation of Tilia cordata in Denmark based on tests of 91 trees selected from 12 isolated populations. We used quantitative analysis of spring phenology and population genetic analysis based on SSR markers to infer the likely historical genetic processes within and among populations. High genetic variation within and among populations was observed in spring phenology, which correlated with spring temperatures at the origin of the tested T. cordata trees. The population genetic analysis revealed significant differentiation among the populations, but with no clear sign of isolation by distance. We infer the findings as indications of ongoing fine scale selection in favor of local growth conditions made possible by limited gene flow among the small and fragmented populations. This hypothesis fits well with reports of limited fruiting in the investigated Danish T. cordata populations, while the species is known for its ability to propagate vegetatively by root suckers. Our results suggest that both divergent selection and genetic drift may have played important roles in forming the genetic patterns of T. cordata at its northern distribution limit. However, we also speculate that epigenetic mechanism arising from the original population environment could have created similar patterns in regulating the spring phenology.     

Sequencing of Australian wild rice genomes reveals ancestral relationships with domesticated rice

The related A genome species of the Oryza genus are the effective gene pool for rice. Here, we report draft genomes for two Australian wild A genome taxa: O. rufipogon‐like population, referred to as Taxon A, and O. meridionalis‐like population, referred to as Taxon B. These two taxa were sequenced and assembled by integration of short‐ and long‐read next‐generation sequencing (NGS) data to create a genomic platform for a wider rice gene pool. Here, we report that, despite the distinct chloroplast genome, the nuclear genome of the Australian Taxon A has a sequence that is much closer to that of domesticated rice (O. sativa) than to the other Australian wild populations. Analysis of 4643 genes in the A genome clade showed that the Australian annual, O. meridionalis, and related perennial taxa have the most divergent (around 3 million years) genome sequences relative to domesticated rice. A test for admixture showed possible introgression into the Australian Taxon A (diverged around 1.6 million years ago) especially from the wild indica/O. nivara clade in Asia. These results demonstrate that northern Australia may be the centre of diversity of the A genome Oryza and suggest the possibility that this might also be the centre of origin of this group and represent an important resource for rice improvement.     

Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species

High‐throughput DNA sequencing facilitates the analysis of large portions of the genome in nonmodel organisms, ensuring high accuracy of population genetic parameters. However, empirical studies evaluating the appropriate sample size for these kinds of studies are still scarce. In this study, we use double‐digest restriction‐associated DNA sequencing (ddRADseq) to recover thousands of single nucleotide polymorphisms (SNPs) for two physically isolated populations of Amphirrhox longifolia (Violaceae), a nonmodel plant species for which no reference genome is available. We used resampling techniques to construct simulated populations with a random subset of individuals and SNPs to determine how many individuals and biallelic markers should be sampled for accurate estimates of intra‐ and interpopulation genetic diversity. We identified 3646 and 4900 polymorphic SNPs for the two populations of A. longifolia, respectively. Our simulations show that, overall, a sample size greater than eight individuals has little impact on estimates of genetic diversity within A. longifolia populations, when 1000 SNPs or higher are used. Our results also show that even at a very small sample size (i.e. two individuals), accurate estimates of F_ST can be obtained with a large number of SNPs (≥1500). These results highlight the potential of high‐throughput genomic sequencing approaches to address questions related to evolutionary biology in nonmodel organisms. Furthermore, our findings also provide insights into the optimization of sampling strategies in the era of population genomics.     

Speciation processes in putative island endemic sister bat species: false impressions from mitochondrial DNA and microsatellite data

Cases of geographically restricted co‐occurring sister taxa are rare and may point to potential divergence with gene flow. The two bat species Murina gracilis and Murina recondita are both endemic to Taiwan and are putative sister species. To test for nonallopatric divergence and gene flow in these taxa, we generated sequences using Sanger and next‐generation sequencing, and combined these with microsatellite data for coalescent‐based analyses. MtDNA phylogenies supported the reciprocally monophyletic sister relationship between M. gracilis and M. recondita; however, clustering of microsatellite genotypes revealed several cases of species admixture suggesting possible introgression. Sequencing of microsatellite flanking regions revealed that admixture signatures stemmed from microsatellite allele homoplasy rather than recent introgressive hybridization, and also uncovered an unexpected sister relationship between M. recondita and the continental species Murina eleryi, to the exclusion of M. gracilis. To dissect the basis of these conflicts between ncDNA and mtDNA, we analysed sequences from 10 anonymous ncDNA loci with *beast and isolation‐with‐migration and found two distinct clades of M. eleryi, one of which was sister to M. recondita. We conclude that Taiwan was colonized by the ancestor of M. gracilis first, followed by the ancestor of M. recondita after a period of allopatric divergence. After colonization, the mitochondrial genome of M. recondita was replaced by that of the resident M. gracilis. This study illustrates how apparent signatures of sympatric divergence can arise from complex histories of allopatric divergence, colonization and hybridization, thus highlighting the need for rigorous analyses to distinguish between such scenarios.     

Genetic isolation between two recently diverged populations of a symbiotic fungus

Fungi are an omnipresent and highly diverse group of organisms, making up a significant part of eukaryotic diversity. Little is currently known about the drivers of fungal population differentiation and subsequent divergence of species, particularly in symbiotic, mycorrhizal fungi. Here, we investigate the population structure and environmental adaptation in Suillus brevipes (Peck) Kuntze, a wind‐dispersed soil fungus that is symbiotic with pine trees. We assembled and annotated the reference genome for Su. brevipes and resequenced the whole genomes of 28 individuals from coastal and montane sites in California. We detected two clearly delineated coast and mountain populations with very low divergence. Genomic divergence was restricted to few regions, including a region of extreme divergence containing a gene encoding for a membrane Na⁺/H⁺ exchanger known for enhancing salt tolerance in plants and yeast. Our results are consistent with a very recent split between the montane and coastal Su. brevipes populations, with few small genomic regions under positive selection and a pattern of dispersal and/or establishment limitation. Furthermore, we identify a putatively adaptive gene that motivates further functional analyses to link genotypes and phenotypes and shed light on the genetic basis of adaptive traits.     

Genetic structure and demographic history of the endangered tree species Dysoxylum malabaricum (Meliaceae) in Western Ghats,India: implications for conservation in a biodiversity hotspot

The impact of fragmentation by human activities on genetic diversity of forest trees is an important concern in forest conservation, especially in tropical forests. Dysoxylum malabaricum (white cedar) is an economically important tree species, endemic to the Western Ghats, India, one of the world's eight most important biodiversity hotspots. As D. malabaricum is under pressure of disturbance and fragmentation together with overharvesting, conservation efforts are required in this species. In this study, range‐wide genetic structure of twelve D. malabaricum populations was evaluated to assess the impact of human activities on genetic diversity and infer the species’ evolutionary history, using both nuclear and chloroplast (cp) DNA simple sequence repeats (SSR). As genetic diversity and population structure did not differ among seedling, juvenile and adult age classes, reproductive success among the old‐growth trees and long distance seed dispersal by hornbills were suggested to contribute to maintain genetic diversity. The fixation index (F_IS) was significantly correlated with latitude, with a higher level of inbreeding in the northern populations, possibly reflecting a more severe ecosystem disturbance in those populations. Both nuclear and cpSSRs revealed northern and southern genetic groups with some discordance of their distributions; however, they did not correlate with any of the two geographic gaps known as genetic barriers to animals. Approximate Bayesian computation‐based inference from nuclear SSRs suggested that population divergence occurred before the last glacial maximum. Finally we discussed the implications of these results, in particular the presence of a clear pattern of historical genetic subdivision, on conservation policies.