PERSPECTIVE: HIGHLY VARIABLE LOCI AND THEIR INTERPRETATION IN EVOLUTION AND CONSERVATION

Testing for Hardy-Weinberg proportions: have we lost the plot?

Testing for Hardy-Weinberg proportions (HWP) is routine in almost all genetic studies of natural populations, but many researchers do not demonstrate a full understanding of the purposes of these tests or how to interpret the results. Common problems include a lack of understanding of statistical power and the difference between statistical significance and biological significance, how to interpret results of multiple tests, and how to distinguish between various factors that can cause statistically significant departures. In this perspective, which focuses on analysis of genetic data for nonmodel species, I 1) review factors that can cause departures from HWP at individual loci and linkage disequilibrium (LD) at pairs of loci; 2) discuss commonly used tests for HWP and LD, with an emphasis on multiple-testing issues; 3) show how to distinguish among possible causes of departures from HWP; and 4) outline some simple steps to follow when significant test results are found. Finally, I 5) identify some issues that merit particular attention as we move into an era in which analysis of genomics-scale datasets for nonmodel species is commonplace.

Comparative riverscape genetics reveals reservoirs of genetic diversity for conservation and restoration of Great Plains fishes

We used comparative landscape genetics to examine the relative roles of historical events, intrinsic traits and landscape factors in determining the distribution of genetic diversity of river fishes across the North American Great Plains. Spatial patterns of diversity were overlaid on a patch-based graphical model and then compared within and among three species that co-occurred across five Great Plains watersheds. Species differing in reproductive strategy (benthic vs. pelagic-spawning) were hypothesized to have different patterns of genetic diversity, but the overriding factor shaping contemporary patterns of diversity was the signature of past climates and geological history. Allelic diversity was significantly higher at southern latitudes for Cyprinella lutrensis and Hybognathus placitus, consistent with northward expansion from southern Pleistocene refugia. Within the historical context, all species exhibited lowered occupancy and abundance in heavily fragmented and drier upstream reaches, particularly H. placitus; a pelagic-spawning species, suggesting rates of extirpation have outpaced losses of genetic diversity in this species. Within most tributary basins, genetically diverse populations of each species persisted. Hence, reconnecting genetically diverse populations with those characterized by reduced diversity (regardless of their position within the riverine network) would provide populations with greater genetic and demographic resilience. We discuss cases where cross-basin transfer may be appropriate to enhance genetic diversity and mitigate negative effects of climate change. Overall, striking similarities in genetic patterns and in response to fragmentation and dewatering suggest a common strategy for genetic resource management in this unique riverine fish assemblage.

Challenges in analysis and interpretation of microsatellite data for population genetic studies

Advancing technologies have facilitated the ever-widening application of genetic markers such as microsatellites into new systems and research questions in biology. In light of the data and experience accumulated from several years of using microsatellites, we present here a literature review that synthesizes the limitations of microsatellites in population genetic studies. With a focus on population structure, we review the widely used fixation (F ST) statistics and Bayesian clustering algorithms and find that the former can be confusing and problematic for microsatellites and that the latter may be confounded by complex population models and lack power in certain cases. Clustering, multivariate analyses, and diversity-based statistics are increasingly being applied to infer population structure, but in some instances these methods lack formalization with microsatellites. Migration-specific methods perform well only under narrow constraints. We also examine the use of microsatellites for inferring effective population size, changes in population size, and deeper demographic history, and find that these methods are untested and/or highly context-dependent. Overall, each method possesses important weaknesses for use with microsatellites, and there are significant constraints on inferences commonly made using microsatellite markers in the areas of population structure, admixture, and effective population size. To ameliorate and better understand these constraints, researchers are encouraged to analyze simulated datasets both prior to and following data collection and analysis, the latter of which is formalized within the approximate Bayesian computation framework. We also examine trends in the literature and show that microsatellites continue to be widely used, especially in non-human subject areas. This review assists with study design and molecular marker selection, facilitates sound interpretation of microsatellite data while fostering respect for their practical limitations, and identifies lessons that could be applied toward emerging markers and high-throughput technologies in population genetics.

A nuclear DNA perspective on delineating evolutionarily significant lineages in polyploids: the case of the endangered shortnose sturgeon (Acipenser brevirostrum)

The shortnose sturgeon, Acipenser brevirostrum, oft considered a phylogenetic relic, is listed as an “endangered species threatened with extinction” in the US and “Vulnerable” on the IUCN Red List. Effective conservation of A. brevirostrum depends on understanding its diversity and evolutionary processes, yet challenges associated with the polyploid nature of its nuclear genome have heretofore limited population genetic analysis to maternally inherited haploid characters. We developed a suite of polysomic microsatellite DNA markers and characterized a sample of 561 shortnose sturgeon collected from major extant populations along the North American Atlantic coast. The 181 alleles observed at 11 loci were scored as binary loci and the data were subjected to multivariate ordination, Bayesian clustering, hierarchical partitioning of variance, and among-population distance metric tests. The methods uncovered moderately high levels of gene diversity suggesting population structuring across and within three metapopulations (Northeast, Mid-Atlantic, and Southeast) that encompass seven demographically discrete and evolutionarily distinct lineages. The predicted groups are consistent with previously described behavioral patterns, especially dispersal and migration, supporting the interpretation that A. brevirostrum exhibit adaptive differences based on watershed. Combined with results of prior genetic (mitochondrial DNA) and behavioral studies, the current work suggests that dispersal is an important factor in maintaining genetic diversity in A. brevirostrum and that the basic unit for conservation management is arguably the local population.

Upper bounds on FST in terms of the frequency of the most frequent allele and total homozygosity: the case of a specified number of alleles

FST is one of the most frequently-used indices of genetic differentiation among groups. Though FST takes values between 0 and 1, authors going back to Wright have noted that under many circumstances, FST is constrained to be less than 1. Recently, we showed that at a genetic locus with an unspecified number of alleles, FST for two subpopulations is strictly bounded from above by functions of both the frequency of the most frequent allele (M) and the homozygosity of the total population (HT). In the two-subpopulation case, FST can equal one only when the frequency of the most frequent allele and the total homozygosity are 1/2. Here, we extend this work by deriving strict bounds on FST for two subpopulations when the number of alleles at the locus is specified to be I. We show that restricting to I alleles produces the same upper bound on FST over much of the allowable domain for M and HT, and we derive more restrictive bounds in the windows M∈[1/I,1/(I-1)) and HT∈[1/I,I/(I(2)-1)). These results extend our understanding of the behavior of FST in relation to other population-genetic statistics.

A field guide to whole-genome sequencing, assembly and annotation

Genome sequencing projects were long confined to biomedical model organisms and required the concerted effort of large consortia. Rapid progress in high-throughput sequencing technology and the simultaneous development of bioinformatic tools have democratized the field. It is now within reach for individual research groups in the eco-evolutionary and conservation community to generate de novo draft genome sequences for any organism of choice. Because of the cost and considerable effort involved in such an endeavour, the important first step is to thoroughly consider whether a genome sequence is necessary for addressing the biological question at hand. Once this decision is taken, a genome project requires careful planning with respect to the organism involved and the intended quality of the genome draft. Here, we briefly review the state of the art within this field and provide a step-by-step introduction to the workflow involved in genome sequencing, assembly and annotation with particular reference to large and complex genomes. This tutorial is targeted at scientists with a background in conservation genetics, but more generally, provides useful practical guidance for researchers engaging in whole-genome sequencing projects.

The impact of translocations on neutral and functional genetic diversity within and among populations of the Seychelles warbler

Translocations are an increasingly common tool in conservation. The maintenance of genetic diversity through translocation is critical for both the short- and long-term persistence of populations and species. However, the relative spatio-temporal impacts of translocations on neutral and functional genetic diversity, and how this affects genetic structure among the conserved populations overall, have received little investigation. We compared the impact of translocating different numbers of founders on both microsatellite and major histocompatibility complex (MHC) class I diversity over a 23-year period in the Seychelles warbler (Acrocephalus sechellensis). We found low and stable microsatellite and MHC diversity in the source population and evidence for only a limited loss of either type of diversity in the four new populations. However, we found evidence of significant, but low to moderate, genetic differentiation between populations, with those populations established with fewer founders clustering separately. Stochastic genetic capture (as opposed to subsequent drift) was the main determinant of translocated population diversity. Furthermore, a strong correlation between microsatellite and MHC differentiation suggested that neutral processes outweighed selection in shaping MHC diversity in the new populations. These data provide important insights into how to optimize the use of translocation as a conservation tool.

Adaptive differentiation of traits related to resource use in a desert annual along a resource gradient

• Plant resource-use traits are generally hypothesized to be adaptively differentiated for populations distributed along resource gradients. Although nutrient limitations are expected to select for resource-conservative strategies, water limitations may select for either resource-conservative or -acquisitive strategies. We test whether population differentiation reflects local adaptation for traits associated with resource-use strategies in a desert annual (Helianthus anomalus) distributed along a gradient of positively covarying water and nutrient availability. • We compared quantitative trait variation (Q(ST)) with neutral genetic differentiation (F(ST)), in a common garden glasshouse study, for leaf economics spectrum (LES) and related traits: photosynthesis (A(mass), A(area)), leaf nitrogen (N(mass), N(area)), leaf lifetime (LL), leaf mass per area (LMA), leaf water content (LWC), water-use efficiency (WUE, estimated as δ(13)C) and days to first flower (DFF). • Q(ST)-F(ST) differences support adaptive differentiation for Amass , N(mass), N(area), LWC and DFF. The trait combinations associated with drier and lower fertility sites represent correlated trait evolution consistent with the more resource-acquisitive end of the LES. There was no evidence for adaptive differentiation for A(area), LMA and WUE. • These results demonstrate that hot dry environments can selectively favor correlated evolution of traits contributing to a resource-acquisitive and earlier reproduction 'escape' strategy, despite lower fertility.

DifferInt: compositional differentiation among populations at three levels of genetic integration

In many fields of study, it is important to know how different populations are genetically. Commonly used measures such as FST and its derivatives based on gene identity probabilities do not reliably reflect difference, as they can be maximal when almost all populations are identical and very small when populations are completely distinct. Compositional differentiation, in contrast, is maximal only for completely distinct populations. Moreover, underlying forces of fragmentation that act on single- and multilocus genotypes may be missed if only the allelic level is viewed. For these reasons, descriptive measures of compositional differentiation based on untransformed distributions of genetic types at different levels of genetic integration (alleles, single- and multilocus genotypes) were derived. Here, two measures of complementary differentiation and two new measures of dispersive differentiation are described. One of each considers genic differences between individuals, ensuring their consistency across integration levels. An increase from one level to the next indicates differences among the populations in their forms of gene association. The computer program DifferInt calculates these measures for one or more codominantly expressed gene loci at the gene-pool, single-locus and multilocus levels. Effects of gene pools and gene association on differentiation are compared by permutation analysis. Snail diagrams depict the contribution of each population to differentiation. Results are demonstrated using a numerical example. The applicability of the measures calculated by this program in conservation, phylogeography, mating system analysis and adaptation studies is discussed.

Present-day genetic structure of Atlantic salmon (Salmo salar) in Icelandic rivers and ice-cap retreat models

Due to an improved understanding of past climatological conditions, it has now become possible to study the potential concordance between former climatological models and present-day genetic structure. Genetic variability was assessed in 26 samples from different rivers of Atlantic salmon in Iceland (total of 2,352 individuals), using 15 microsatellite loci. F-statistics revealed significant differences between the majority of the populations that were sampled. Bayesian cluster analyses using both prior information and no prior information on sampling location revealed the presence of two distinguishable genetic pools - namely, the Northern (Group 1) and Southern (Group 2) regions of Iceland. Furthermore, the random permutation of different allele sizes among allelic states revealed a significant mutational component to the genetic differentiation at four microsatellite loci (SsaD144, Ssa171, SSsp2201 and SsaF3), and supported the proposition of a historical origin behind the observed variation. The estimated time of divergence, using two different ABC methods, suggested that the observed genetic pattern originated from between the Last Glacial Maximum to the Younger Dryas, which serves as additional evidence of the relative immaturity of Icelandic fish populations, on account of the re-colonisation of this young environment following the Last Glacial Maximum. Additional analyses suggested the presence of several genetic entities which were likely to originate from the original groups detected.

Patterns of genetic diversity reveal multiple introductions and recurrent founder effects during range expansion in invasive populations of Geranium carolinianum (Geraniaceae)

Genetic diversity, and thus the adaptive potential of invasive populations, is largely based on three factors: patterns of genetic diversity in the species' native range, the number and location of introductions and the number of founding individuals per introduction. Specifically, reductions in genetic diversity ('founder effects') should be stronger for species with low within-population diversity in their native range and few introductions of few individuals to the invasive range. We test these predictions with Geranium carolinianum, a winter annual herb native to North America and invasive in China. We measure the extent of founder effects using allozymes and microsatellites, and ask whether this is consistent with its colonization history and patterns of diversity in the native range. In the native range, genetic diversity is higher and structure is lower than expected based on life history traits. In China, our results provide evidence for multiple introductions near Nanjing, Jiangsu province, with subsequent range expansion to the west and south. Patterns of genetic diversity across China reveal weak founder effects that are driven largely by low-diversity populations at the expansion front, away from the introduction location. This suggests that reduced diversity in China has resulted from successive founder events during range expansion, and that the loss of genetic diversity in the Nanjing area was mitigated by multiple introductions from diverse source populations. This has implications for the future of G. carolinianum in China, as continued gene flow among populations should eventually increase genetic diversity within the more recently founded populations.

Antagonistic selection factors induce a continuous population divergence in a polymorphism

Understanding the relative importance of selection and stochastic factors in population divergence of adaptive traits is a classical topic in evolutionary biology. However, it is difficult to separate these factors and detect the effects of selection when two or more contrasting selective factors are simultaneously acting on a single locus. In the damselfly Ischnura senegalensis, females exhibit color dimorphism and morph frequencies change geographically. We here evaluated the role of selection and stochastic factors in population divergence of morph frequencies by comparing the divergences in color locus and neutral loci. Comparisons between population pairwise FST for neutral loci and for the color locus did not detect any stochastic factors affecting color locus. Although comparison between population divergence in color and neutral loci using all populations detected only divergent selection, we detected two antagonistic selective factors acting on the color locus, that is, balancing and divergent selection, when considering geographical distance between populations. Our results suggest that a combination of two antagonistic selective factors, rather than stochastic factors, establishes the geographic cline in morph frequency in this system.

Significant variance in genetic diversity among populations of Schistosoma haematobium detected using microsatellite DNA loci from a genome-wide database

BACKGROUND

Urogenital schistosomiasis caused by Schistosoma haematobium is widely distributed across Africa and is increasingly being targeted for control. Genome sequences and population genetic parameters can give insight into the potential for population- or species-level drug resistance. Microsatellite DNA loci are genetic markers in wide use by Schistosoma researchers, but there are few primers available for S. haematobium.

METHODS

We sequenced 1,058,114 random DNA fragments from clonal cercariae collected from a snail infected with a single Schistosoma haematobium miracidium. We assembled and aligned the S. haematobium sequences to the genomes of S. mansoni and S. japonicum, identifying microsatellite DNA loci across all three species and designing primers to amplify the loci in S. haematobium. To validate our primers, we screened 32 randomly selected primer pairs with population samples of S. haematobium.

RESULTS

We designed >13,790 primer pairs to amplify unique microsatellite loci in S. haematobium, (available at http://www.cebio.org/projetos/schistosoma-haematobium-genome). The three Schistosoma genomes contained similar overall frequencies of microsatellites, but the frequency and length distributions of specific motifs differed among species. We identified 15 primer pairs that amplified consistently and were easily scored. We genotyped these 15 loci in S. haematobium individuals from six locations: Zanzibar had the highest levels of diversity; Malawi, Mauritius, Nigeria, and Senegal were nearly as diverse; but the sample from South Africa was much less diverse.

CONCLUSIONS

About half of the primers in the database of Schistosoma haematobium microsatellite DNA loci should yield amplifiable and easily scored polymorphic markers, thus providing thousands of potential markers. Sequence conservation among S. haematobium, S. japonicum, and S. mansoni is relatively high, thus it should now be possible to identify markers that are universal among Schistosoma species (i.e., using DNA sequences conserved among species), as well as other markers that are specific to species or species-groups (i.e., using DNA sequences that differ among species). Full genome-sequencing of additional species and specimens of S. haematobium, S. japonicum, and S. mansoni is desirable to better characterize differences within and among these species, to develop additional genetic markers, and to examine genes as well as conserved non-coding elements associated with drug resistance.

Extensive mitochondrial introgression in North American Great Black-backed Gulls (Larus marinus) from the American Herring Gull (Larus smithsonianus) with little nuclear DNA impact

Recent genetic studies have shown that introgression rates among loci may greatly vary according to their location in the genome. In particular, several cases of mito-nuclear discordances have been reported for a wide range of organisms. In the present study, we examine the causes of discordance between mitochondrial (mtDNA) and nuclear DNA introgression detected in North American populations of the Great Black-backed Gull (Larus marinus), a Holarctic species, from the Nearctic North American Herring Gull (Larus smithsonianus). Our results show that extensive unidirectional mtDNA introgression from Larus smithsonianus into Larus marinus in North America cannot be explained by ancestral polymorphism but most likely results from ancient hybridization events occurring when Larus marinus invaded the North America. Conversely, our nuclear DNA results based on 12 microsatellites detected very little introgression from Larus smithsonianus into North American Larus marinus. We discuss these results in the framework of demographic and selective mechanisms that have been postulated to explain mito-nuclear discrepancies. We were unable to demonstrate selection as the main cause of mito-nuclear introgression discordance but cannot dismiss the possible role of selection in the observed pattern. Among demographic explanations, only drift in small populations and bias in mate choice in an invasive context may explain our results. As it is often difficult to demonstrate that selection may be the main factor driving the introgression of mitochondrial DNA in natural populations, we advocate that evaluating alternative demographic neutral hypotheses may help to indirectly support or reject hypotheses invoking selective processes.

Mitochondrial DNA diversity of mud crab Scylla olivacea (Portunidae) in Peninsular Malaysia: a preliminary assessment

A primary factor in population management and wildlife conservation is the delineation of population units derived from descriptions of population genetic structure. Yet, predicting factors that influence the patterns of gene flow in a population particularly at landscape scales remains a major challenge in evolutionary biology. Here we report a population genetic study of the mud crab Scylla olivacea examined based on a 542 bp segment of the mitochondrial DNA cytochrome c oxidase I gene among 91 individuals from six localities in the west and east coast of Peninsular Malaysia. In total 55 unique haplotypes were distinguished with 45 private haplotypes and a single common haplotype shared among all populations studied. The other ten haplotypes were shared among various populations. The sharing of this haplotype reflects the connection of the mangrove areas between east and west coast of Peninsular Malaysia. High haplotype diversity (h = 0.968 ± 0.021; mean ± SD) and low nucleotide diversity (π = 0.120 ± 0.015; mean ± SD) were displayed, which may be indicative of genetic bottleneck events. No significant phylogenetic lineages were recognized using neighbour-joining and maximum parsimony methods. Hierarchical AMOVA analysis indicated that 99.33 % of the genetic variation was contained within populations and 0.67 % occurred among populations, suggesting no geographical patterning among populations studied, supported by F st test. Mismatch distribution analysis showed that the observed distribution of the pairwise mutation differences among haplotypes was multimodal, which is not concordant with a sudden range expansion scenario. However, neutrality tests showed non-significant negative values suggesting that the populations studied may have experienced past population growth, but the expansion may have been restricted to separate local areas that resulted in the non-significant negative Fu's Fs and Tajima's D value. Overall, this present preliminary study was able to be a reference on the phylogenetic relationships and assessment of genetic structure of Scylla sp. in Malaysia.

Polarisation of major histocompatibility complex II host genotype with pathogenesis of European Brown Hare syndrome virus

A study was conducted in order to determine the occurrence of European Brown Hare Syndrome virus (EBHSV) in Denmark and possible relation between disease pathogenesis and Major Histocompatibility Complex (MHC) host genotype. Liver samples were examined from 170 brown hares (hunted, found sick or dead), collected between 2004 and 2009. Macroscopical and histopathological findings consistent with EBHS were detected in 24 (14.1%) hares; 35 (20.6%) had liver lesions not typical of the syndrome, 50 (29.4%) had lesions in other tissues and 61 (35.9%) had no lesions. Sixty five (38.2%) of 170 samples were found to be EBHSV-positive (RT-PCR, VP60 gene). In order to investigate associations between viral pathogenesis and host genotype, variation within the exon 2 DQA gene of MHC was assessed. DQA exon 2 analysis revealed the occurrence of seven different alleles in Denmark. Consistent with other populations examined so far in Europe, observed heterozygosity of DQA (H_o = 0.1180) was lower than expected (H_e = 0.5835). The overall variation for both nucleotide and amino acid differences (2.9% and 14.9%, respectively) were lower in Denmark than those assessed in other European countries (8.3% and 16.9%, respectively). Within the peptide binding region codons the number of nonsynonymous substitutions (dN) was much higher than synonymous substitutions (dS), which would be expected for MHC alleles under balancing selection. Allele frequencies did not significantly differ between EBHSV-positive and -negative hares. However, allele Leeu-DQA*30 was detected in significantly higher (P = 0.000006) frequency among the positive hares found dead with severe histopathological lesions than among those found sick or apparently healthy. In contrast, the latter group was characterized by a higher frequency of the allele Leeu-DQA*14 as well as the proportion of heterozygous individuals (P = 0.000006 and P = 0.027). These data reveal a polarisation between EBHSV pathogenesis and MHC class II genotype within the European brown hare in Denmark.

Genetic and morphological structure of a spruce hybrid (Picea sitchensis x P. glauca) zone along a climatic gradient

PREMISE OF THE STUDY

Historic colonization and contemporary evolutionary processes contribute to patterns of genetic variation and differentiation among populations. However, separating the respective influences of these processes remains a challenge, particularly for natural hybrid zones, where standing genetic variation may result from evolutionary processes both preceding and following contact, influencing the evolutionary trajectory of hybrid populations. Where adaptation to novel environments may be facilitated by interspecific hybridization, teasing apart these processes will have practical implications for forest management in changing environments.

METHODS

We evaluated the neutral genetic architecture of the Picea sitchensis (Sitka spruce) × P. glauca (white spruce) hybrid zone along the Nass and Skeena river valleys in northwestern British Columbia using chloroplast, mitochondrial, and nuclear microsatellite markers, in combination with cone morphological traits.

KEY RESULTS

Sitka spruce mitotype "capture", evidenced by this species dominating the maternal lineage, is consistent with earlier colonization of the region by Sitka spruce. This "capture" differs from the spatial distribution of chloroplast haplotypes, indicating pollen dispersal and its contribution to geographic structure. Genetic ancestry, based on nuclear markers, was strongly influenced by climate and geography. Highly parallel results for replicate transects along environmental gradients provide support for the bounded hybrid superiority model of hybrid zone maintenance. •

CONCLUSIONS

This broad-scale analysis of neutral genetic structure indicates the importance of historic and contemporary gene flow, environmental selection, and their interaction in shaping neutral genetic variation within this hybrid zone, informative to seed transfer development and reforestation for future climates.

Combined analyses of kinship and FST suggest potential drivers of chaotic genetic patchiness in high gene-flow populations

We combine kinship estimates with traditional F-statistics to explain contemporary drivers of population genetic differentiation despite high gene flow. We investigate range-wide population genetic structure of the California spiny (or red rock) lobster (Panulirus interruptus) and find slight, but significant global population differentiation in mtDNA (ΦST = 0.006, P = 0.001; D(est_Chao) = 0.025) and seven nuclear microsatellites (F(ST) = 0.004, P < 0.001; D(est_Chao) = 0.03), despite the species' 240- to 330-day pelagic larval duration. Significant population structure does not correlate with distance between sampling locations, and pairwise FST between adjacent sites often exceeds that among geographically distant locations. This result would typically be interpreted as unexplainable, chaotic genetic patchiness. However, kinship levels differ significantly among sites (pseudo-F(16,988) = 1.39, P = 0.001), and ten of 17 sample sites have significantly greater numbers of kin than expected by chance (P < 0.05). Moreover, a higher proportion of kin within sites strongly correlates with greater genetic differentiation among sites (D(est_Chao), R(2) = 0.66, P < 0.005). Sites with elevated mean kinship were geographically proximate to regions of high upwelling intensity (R(2) = 0.41, P = 0.0009). These results indicate that P. interruptus does not maintain a single homogenous population, despite extreme dispersal potential. Instead, these lobsters appear to either have substantial localized recruitment or maintain planktonic larval cohesiveness whereby siblings more likely settle together than disperse across sites. More broadly, our results contribute to a growing number of studies showing that low F(ST) and high family structure across populations can coexist, illuminating the foundations of cryptic genetic patterns and the nature of marine dispersal.

A Large Plasmodium vivax Reservoir and Little Population Structure in the South Pacific

Introduction

The importance of Plasmodium vivax in malaria elimination is increasingly being recognized, yet little is known about its population size and population genetic structure in the South Pacific, an area that is the focus of intensified malaria control.

Methods

We have genotyped 13 microsatellite markers in 295 P. vivax isolates from four geographically distinct sites in Papua New Guinea (PNG) and one site from Solomon Islands, representing different transmission intensities.

Results

Diversity was very high with expected heterozygosity values ranging from 0.62 to 0.98 for the different markers. Effective population size was high (12′872 to 19′533 per site). In PNG population structuring was limited with moderate levels of genetic differentiation. F_ST values (adjusted for high diversity of markers) were 0.14–0.15. Slightly higher levels were observed between PNG populations and Solomon Islands (F_ST = 0.16).

Conclusions

Low levels of population structure despite geographical barriers to transmission are in sharp contrast to results from regions of low P. vivax endemicity. Prior to intensification of malaria control programs in the study area, parasite diversity and effective population size remained high.

Testing the role of meander cutoff in promoting gene flow across a riverine barrier in ground skinks (Scincella lateralis)

Despite considerable attention, the long-term impact of rivers on species diversification remains uncertain. Meander loop cutoff (MLC) is one river phenomenon that may compromise a river's diversifying effects by passively transferring organisms from one side of the river to the other. However, the ability of MLC to promote gene flow across rivers has not been demonstrated empirically. Here, we test several predictions of MLC-mediated gene flow in populations of North American ground skinks (Scincella lateralis) separated by a well-established riverine barrier, the Mississippi River: 1) individuals collected from within meander cutoffs should be more closely related to individuals across the river than on the same side, 2) individuals within meander cutoffs should contain more immigrants than individuals away from meander cutoffs, 3) immigration rates estimated across the river should be highest in the direction of the cutoff event, and 4) the distribution of alleles native to one side of the river should be better predicted by the historical rather than current path of the river. To test these predictions we sampled 13 microsatellite loci and mitochondrial DNA from ground skinks collected near three ancient meander loops. These predictions were generally supported by genetic data, although support was stronger for mtDNA than for microsatellite data. Partial support for genetic divergence of samples within ancient meander loops also provides evidence for the MLC hypothesis. Although a role for MLC-mediated gene flow was supported here for ground skinks, the transient nature of river channels and morphologies may limit the long-term importance of MLC in stemming population divergence across major rivers.

Demographic history and the low genetic diversity in Dipteryx alata (Fabaceae) from Brazilian Neotropical savannas

Genetic effects of habitat fragmentation may be undetectable because they are generally a recent event in evolutionary time or because of confounding effects such as historical bottlenecks and historical changes in species' distribution. To assess the effects of demographic history on the genetic diversity and population structure in the Neotropical tree Dipteryx alata (Fabaceae), we used coalescence analyses coupled with ecological niche modeling to hindcast its distribution over the last 21 000 years. Twenty-five populations (644 individuals) were sampled and all individuals were genotyped using eight microsatellite loci. All populations presented low allelic richness and genetic diversity. The estimated effective population size was small in all populations and gene flow was negligible among most. We also found a significant signal of demographic reduction in most cases. Genetic differentiation among populations was significantly correlated with geographical distance. Allelic richness showed a spatial cline pattern in relation to the species' paleodistribution 21 kyr BP (thousand years before present), as expected under a range expansion model. Our results show strong evidences that genetic diversity in D. alata is the outcome of the historical changes in species distribution during the late Pleistocene. Because of this historically low effective population size and the low genetic diversity, recent fragmentation of the Cerrado biome may increase population differentiation, causing population decline and compromising long-term persistence.

The environment affects epistatic interactions to alter the topology of an empirical fitness landscape

The fitness effect of mutations can be influenced by their interactions with the environment, other mutations, or both. Previously, we constructed 32 ( = 2⁵) genotypes that comprise all possible combinations of the first five beneficial mutations to fix in a laboratory-evolved population of Escherichia coli. We found that (i) all five mutations were beneficial for the background on which they occurred; (ii) interactions between mutations drove a diminishing returns type epistasis, whereby epistasis became increasingly antagonistic as the expected fitness of a genotype increased; and (iii) the adaptive landscape revealed by the mutation combinations was smooth, having a single global fitness peak. Here we examine how the environment influences epistasis by determining the interactions between the same mutations in two alternative environments, selected from among 1,920 screened environments, that produced the largest increase or decrease in fitness of the most derived genotype. Some general features of the interactions were consistent: mutations tended to remain beneficial and the overall pattern of epistasis was of diminishing returns. Other features depended on the environment; in particular, several mutations were deleterious when added to specific genotypes, indicating the presence of antagonistic interactions that were absent in the original selection environment. Antagonism was not caused by consistent pleiotropic effects of individual mutations but rather by changing interactions between mutations. Our results demonstrate that understanding adaptation in changing environments will require consideration of the combined effect of epistasis and pleiotropy across environments.

The fitness effect of beneficial mutations can depend on how they interact with their genetic and external environment. The form of these interactions is important because it can alter adaptive outcomes, selecting for or against certain combinations of beneficial mutations. Here, we examine how interactions between beneficial mutations favored during adaptation of a lab strain of Escherichia coli to one simple environment are altered when the strain is grown in two novel environments. We found that fitness effects were greatly influenced by both the genetic and external environments. In several instances a change in environment reversed the effect of a mutation from beneficial to deleterious or caused combinations of beneficial mutations to become deleterious. Our results suggest that a complex or fluctuating environment may favor combinations of mutations whose interactions may be less sensitive to external conditions.

Analysis of Geographic and Pairwise Distances among Chinese Cashmere Goat Populations

This study investigated the geographic and pairwise distances of nine Chinese local Cashmere goat populations through the analysis of 20 microsatellite DNA markers. Fluorescence PCR was used to identify the markers, which were selected based on their significance as identified by the Food and Agriculture Organization of the United Nations (FAO) and the International Society for Animal Genetics (ISAG). In total, 206 alleles were detected; the average allele number was 10.30; the polymorphism information content of loci ranged from 0.5213 to 0.7582; the number of effective alleles ranged from 4.0484 to 4.6178; the observed heterozygosity was from 0.5023 to 0.5602 for the practical sample; the expected heterozygosity ranged from 0.5783 to 0.6464; and Allelic richness ranged from 4.7551 to 8.0693. These results indicated that Chinese Cashmere goat populations exhibited rich genetic diversity. Further, the Wright’s F-statistics of subpopulation within total (FST) was 0.1184; the genetic differentiation coefficient (GST) was 0.0940; and the average gene flow (Nm) was 2.0415. All pairwise FST values among the populations were highly significant (p<0.01 or p<0.001), suggesting that the populations studied should all be considered to be separate breeds. Finally, the clustering analysis divided the Chinese Cashmere goat populations into at least four clusters, with the Hexi and Yashan goat populations alone in one cluster. These results have provided useful, practical, and important information for the future of Chinese Cashmere goat breeding.

Phylogeographic patterns of mtDNA variation revealed multiple glacial refugia for the frog species Feirana taihangnica endemic to the Qinling Mountains

Diversification patterns and demography of montane species are affected by Pleistocene climate fluctuations. Empirical cases from the Qinling Mountains (QM) region, which is a major biogeographic divider of East Asia, are few. We used DNA sequence data of the complete mitochondrial ND2 gene to detect effects of the Pleistocene glaciations on phylogeographic profiles of a frog species, Feirana taihangnica, which is endemic to the QM. Four distinct lineages consisting of seven sublineages were revealed. The strongest signal of biogeographical structure (F(ct) = 0.971, P < 0.01) was found when populations were grouped according to these seven sublineages. One narrow secondary contact zone was detected in the middle QM between the lineage from middle QM and the lineage from eastern QM. Coalescent simulations indicated that this species colonized the QM region by a stepping-stone model. Divergences among lineages had likely been influenced by the uplift of the Tibetan Plateau during the late Miocene-to-late Pleistocene, as well as by the Pleistocene climatic cycles. Coalescent simulations also suggested that F. taihangnica populations have persisted through the Pleistocene glacial periods in multiple refugia across the QM region. Demographic analyses indicated that all lineages, except the lineage in the Funiu Mountains, have been experienced postglacial expansion of population size and distribution range. In conclusion, Pleistocene climate fluctuations and tectonic changes during the late Miocene-late Pleistocene have profoundly influenced the phylogeography and historical demography of F. taihangnica.

On the measurements of genetic differentiation among populations

FST, a measurement of the genetic differentiation among subpopulations, is a fundamental parameter in population genetics, with many valuable applications in molecular biology, evolutionary biology, conservation and forensics. One of its close relatives, GST, has been widely used to measure differentiation from highly polymorphic markers such as microsatellites. However, because of the high mutation rate of such markers, GST may underestimate the genomic differentiation due to demographic causes such as migration rate and subpopulation size. A new statistic proposed recently, Jost's D, was claimed to have better properties than GST and was advocated to replace GST as a measure of differentiation. This paper shows that D is not a proper measure of differentiation because it fails to meet some fundamental requirements as a differentiation statistic, and is hardly estimable without bias in practice. D is highly dependent on the gene diversity of a marker and on the unknown parameter of the number of subpopulations, is highly sensitive to how alleles and loci are defined and how data are analysed, does not increase monotonically with either divergence time or drift, and does not always have a maximal value of 1. The maximal D value can be zero or close to zero, depending on the number of alleles at a locus relative to the number of subpopulations. I suggest continuing the use of GST, with caution in its interpretation when highly polymorphic markers are used, before a better estimator of FST that explicitly accounts for mutations is developed.

Unraveling the effects of selection and demography on immune gene variation in free-ranging plains zebra (Equus quagga) populations

Demography, migration and natural selection are predominant processes affecting the distribution of genetic variation among natural populations. Many studies use neutral genetic markers to make inferences about population history. However, the investigation of functional coding loci, which directly reflect fitness, is critical to our understanding of species' ecology and evolution. Immune genes, such as those of the Major Histocompatibility Complex (MHC), play an important role in pathogen recognition and provide a potent model system for studying selection. We contrasted diversity patterns of neutral data with MHC loci, ELA-DRA and -DQA, in two southern African plains zebra (Equus quagga) populations: Etosha National Park, Namibia, and Kruger National Park, South Africa. Results from neutrality tests, along with observations of elevated diversity and low differentiation across populations, supported previous genus-level evidence for balancing selection at these loci. Despite being low, MHC divergence across populations was significant and may be attributed to drift effects typical of geographically separated populations experiencing little to no gene flow, or alternatively to shifting allele frequency distributions driven by spatially variable and fluctuating pathogen communities. At the DRA, zebra exhibited geographic differentiation concordant with microsatellites and reduced levels of diversity in Etosha due to highly skewed allele frequencies that could not be explained by demography, suggestive of spatially heterogeneous selection and local adaptation. This study highlights the complexity in which selection affects immune gene diversity and warrants the need for further research on the ecological mechanisms shaping patterns of adaptive variation among natural populations.

Purifying selection modulates the estimates of population differentiation and confounds genome-wide comparisons across single-nucleotide polymorphisms

An improved understanding of the biological and numerical properties of measures of population differentiation across loci is becoming increasingly more important because of their growing use in analyzing genome-wide polymorphism data for detecting population structures, inferring the rates of migration, and identifying local adaptations. In a genome-wide analysis, we discovered that the estimates of population differentiation (e.g., F(ST), θ, and Jost's D) calculated for human single-nucleotide polymorphisms (SNPs) are strongly and positively correlated to the position-specific evolutionary rates measured from multispecies alignments. That is, genomic positions (loci) experiencing higher purifying selection (lower evolutionary rates) produce lower values for the degree of population differentiation than those evolving with faster rates. We show that this pattern is completely mediated by the negative effects of purifying selection on the minor allele frequency (MAF) at individual loci. Our results suggest that inferences and methods relying on the comparison of population differentiation estimates (F(ST), θ, and Jost's D) based on SNPs across genomic positions should be restricted to loci with similar MAFs and/or the rates of evolution in genome scale surveys.

Human neutral genetic variation and forensic STR data

The forensic genetics field is generating extensive population data on polymorphism of short tandem repeats (STR) markers in globally distributed samples. In this study we explored and quantified the informative power of these datasets to address issues related to human evolution and diversity, by using two online resources: an allele frequency dataset representing 141 populations summing up to almost 26 thousand individuals; a genotype dataset consisting of 42 populations and more than 11 thousand individuals. We show that the genetic relationships between populations based on forensic STRs are best explained by geography, as observed when analysing other worldwide datasets generated specifically to study human diversity. However, the global level of genetic differentiation between populations (as measured by a fixation index) is about half the value estimated with those other datasets, which contain a much higher number of markers but much less individuals. We suggest that the main factor explaining this difference is an ascertainment bias in forensics data resulting from the choice of markers for individual identification. We show that this choice results in average low variance of heterozygosity across world regions, and hence in low differentiation among populations. Thus, the forensic genetic markers currently produced for the purpose of individual assignment and identification allow the detection of the patterns of neutral genetic structure that characterize the human population but they do underestimate the levels of this genetic structure compared to the datasets of STRs (or other kinds of markers) generated specifically to study the diversity of human populations.

The relationship between F(ST) and the frequency of the most frequent allele

F(ST) is frequently used as a summary of genetic differentiation among groups. It has been suggested that F(ST) depends on the allele frequencies at a locus, as it exhibits a variety of peculiar properties related to genetic diversity: higher values for biallelic single-nucleotide polymorphisms (SNPs) than for multiallelic microsatellites, low values among high-diversity populations viewed as substantially distinct, and low values for populations that differ primarily in their profiles of rare alleles. A full mathematical understanding of the dependence of F(ST) on allele frequencies, however, has been elusive. Here, we examine the relationship between F(ST) and the frequency of the most frequent allele, demonstrating that the range of values that F(ST) can take is restricted considerably by the allele-frequency distribution. For a two-population model, we derive strict bounds on F(ST) as a function of the frequency M of the allele with highest mean frequency between the pair of populations. Using these bounds, we show that for a value of M chosen uniformly between 0 and 1 at a multiallelic locus whose number of alleles is left unspecified, the mean maximum F(ST) is ∼0.3585. Further, F(ST) is restricted to values much less than 1 when M is low or high, and the contribution to the maximum F(ST) made by the most frequent allele is on average ∼0.4485. Using bounds on homozygosity that we have previously derived as functions of M, we describe strict bounds on F(ST) in terms of the homozygosity of the total population, finding that the mean maximum F(ST) given this homozygosity is 1 - ln 2 ≈ 0.3069. Our results provide a conceptual basis for understanding the dependence of F(ST) on allele frequencies and genetic diversity and for interpreting the roles of these quantities in computations of F(ST) from population-genetic data. Further, our analysis suggests that many unusual observations of F(ST), including the relatively low F(ST) values in high-diversity human populations from Africa and the relatively low estimates of F(ST) for microsatellites compared to SNPs, can be understood not as biological phenomena associated with different groups of populations or classes of markers but rather as consequences of the intrinsic mathematical dependence of F(ST) on the properties of allele-frequency distributions.

Birds in space and time: genetic changes accompanying anthropogenic habitat fragmentation in the endangered black-capped vireo (Vireo atricapilla)

Anthropogenic alterations in the natural environment can be a potent evolutionary force. For species that have specific habitat requirements, habitat loss can result in substantial genetic effects, potentially impeding future adaptability and evolution. The endangered black-capped vireo (Vireo atricapilla) suffered a substantial contraction of breeding habitat and population size during much of the 20th century. In a previous study, we reported significant differentiation between remnant populations, but failed to recover a strong genetic signal of bottlenecks. In this study, we used a combination of historical and contemporary sampling from Oklahoma and Texas to (i) determine whether population structure and genetic diversity have changed over time and (ii) evaluate alternate demographic hypotheses using approximate Bayesian computation (ABC). We found lower genetic diversity and increased differentiation in contemporary samples compared to historical samples, indicating nontrivial impacts of fragmentation. ABC analysis suggests a bottleneck having occurred in the early part of the 20th century, resulting in a magnitude decline in effective population size. Genetic monitoring with temporally spaced samples, such as used in this study, can be highly informative for assessing the genetic impacts of anthropogenic fragmentation on threatened or endangered species, as well as revealing the dynamics of small populations over time.

MMOD: an R library for the calculation of population differentiation statistics

MMOD is a library for the R programming language that allows the calculation of the population differentiation measures D(est), G″(ST) and φ'(ST). R provides a powerful environment in which to conduct and record population genetic analyses but, at present, no R libraries provide functions for the calculation of these statistics from standard population genetic files. In addition to the calculation of differentiation measures, mmod can produce parametric bootstrap and jackknife samples of data sets for further analysis. By integrating with and complimenting the existing libraries adegenet and pegas, mmod extends the power of R as a population genetic platform.

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.

The role of the Ord Arid Intrusion in the historical and contemporary genetic division of long-tailed finch subspecies in northern Australia

The effect of separation by biogeographic features followed by secondary contact can blur taxonomic boundaries and produce complex genetic signatures. We analyzed population structure and gene flow across the range of the long-tailed finch (Poephila acuticauda) in northern Australia (1) to test the hypothesis that Ord Arid Intrusion acted as the causative barrier that led to divergence of P. acuticauda subspecies, (2) to determine whether genetic data support the presence of a gradual cline across the range or a sudden shift, both of which have been suggested based on morphological data, and (3) to estimate levels of contemporary gene flow within this species complex. We collected samples from 302 individuals from 10 localities. Analyses of 12 microsatellite loci and sequence data from 333 base pairs of the mitochondrial control region were used to estimate population structure and gene flow, using analysis of molecular variance (AMOVA), haplotype network analysis, frequency statistics, and clustering methods. Mitochondrial sequence data indicated the presence of three genetic groups (regions) across the range of P. acuticauda. Genetic diversity was highest in the east and lowest in the west. The Ord Arid Intrusion appears to have functioned as a biogeographic barrier in the past, according to mtDNA evidence presented here and evidence from previous studies. The absence of isolation by distance between adjacent regions and the lack of population genetic structure of mtDNA within regions indicates that genetic changes across the range of P. acuticauda subspecies are characterized by discrete breaks between regions. While microsatellite data indicate a complete absence of genetic structure across this species' range, it appears unlikely that this results from high levels of gene flow. Mitochondrial data do not support the presence of contemporary gene flow across the range of this species.

Genetic consequences of climate change for northern plants

Climate change will lead to loss of range for many species, and thus to loss of genetic diversity crucial for their long-term persistence. We analysed range-wide genetic diversity (amplified fragment length polymorphisms) in 9581 samples from 1200 populations of 27 northern plant species, to assess genetic consequences of range reduction and potential association with species traits. We used species distribution modelling (SDM, eight techniques, two global circulation models and two emission scenarios) to predict loss of range and genetic diversity by 2080. Loss of genetic diversity varied considerably among species, and this variation could be explained by dispersal adaptation (up to 57%) and by genetic differentiation among populations (F(ST); up to 61%). Herbs lacking adaptations for long-distance dispersal were estimated to lose genetic diversity at higher rate than dwarf shrubs adapted to long-distance dispersal. The expected range reduction in these 27 northern species was larger than reported for temperate plants, and all were predicted to lose genetic diversity according to at least one scenario. SDM combined with F(ST) estimates and/or with species trait information thus allows the prediction of species' vulnerability to climate change, aiding rational prioritization of conservation efforts.

Identification of novel single nucleotide polymorphisms (SNPs) in deer (Odocoileus spp.) using the BovineSNP50 BeadChip

Single nucleotide polymorphisms (SNPs) are growing in popularity as a genetic marker for investigating evolutionary processes. A panel of SNPs is often developed by comparing large quantities of DNA sequence data across multiple individuals to identify polymorphic sites. For non-model species, this is particularly difficult, as performing the necessary large-scale genomic sequencing often exceeds the resources available for the project. In this study, we trial the Bovine SNP50 BeadChip developed in cattle (Bos taurus) for identifying polymorphic SNPs in cervids Odocoileus hemionus (mule deer and black-tailed deer) and O. virginianus (white-tailed deer) in the Pacific Northwest. We found that 38.7% of loci could be genotyped, of which 5% (n = 1068) were polymorphic. Of these 1068 polymorphic SNPs, a mixture of putatively neutral loci (n = 878) and loci under selection (n = 190) were identified with the F_ST-outlier method. A range of population genetic analyses were implemented using these SNPs and a panel of 10 microsatellite loci. The three types of deer could readily be distinguished with both the SNP and microsatellite datasets. This study demonstrates that commercially developed SNP chips are a viable means of SNP discovery for non-model organisms, even when used between very distantly related species (the Bovidae and Cervidae families diverged some 25.1−30.1 million years before present).

Genetic monitoring and complex population dynamics: insights from a 12-year study of the Rio Grande silvery minnow

The endangered Rio Grande silvery minnow persists as a remnant population in a highly fragmented and regulated arid-land river system. The species is subject to dramatic annual fluctuations in density. Since 2003, the wild population has been supplemented by hatchery-reared fish. We report on a 12-year (1999–2010) monitoring study of genetic diversity and effective population size (N_e) of wild and hatchery stocks. Our goals were to evaluate how genetic metrics responded to changes in wild fish density and whether they corresponded to the number and levels of diversity of hatchery-reared repatriates. Genetic diversity and all measures of N_e in the wild population did not correlate with wild fish density until hatchery supplementation began in earnest. Estimates of variance and inbreeding effective size were not correlated. Our results suggest source–sink dynamics where captive stocks form a genetically diverse source and the wild population behaves as a sink. Nevertheless, overall genetic diversity of silvery minnow has been maintained over the last decade, and we attribute this to a well-designed and executed propagation management plan. When multiple factors like environmental fluctuation and hatchery supplementation act simultaneously on a population, interpretation of genetic monitoring data may be equally complex and require considerable ecological data.

Genetic evidence for hybridization in red oaks (Quercus sect. Lobatae, Fagaceae)

PREMISE OF THE STUDY

Hybridization is pervasive in many plant taxa, with consequences for species taxonomy, local adaptation, and management. Oaks (Quercus spp.) are thought to hybridize readily yet retain distinct traits, drawing into question the biological species concept for such taxa, but the true extent of gene flow is controversial. Genetic data are beginning to shed new light on this issue, but red oaks (section Lobatae), an important component of North American forests, have largely been neglected. Moreover, gene flow estimates may be sensitive to the choice of life stage, marker type, or genetic structure statistic.

METHODS

We coupled genetic structure data with parentage analyses for two mixed-species stands in North Carolina. Genetic structure analyses of adults (including F(ST), R(ST), G'(ST), and structure) reflect long-term patterns of gene flow, while the percentage of seedlings with parents of two different species reflect current levels of gene flow.

KEY RESULTS

Genetic structure analyses revealed low differentiation in microsatellite allele frequencies between co-occurring species, suggesting past gene flow. However, methods differed in their sensitivity to differentiation, indicating a need for caution when drawing conclusions from a single method. Parentage analyses identified >20% of seedlings as potential hybrids. The species examined exhibit distinct morphologies, suggesting selection against intermediate phenotypes.

CONCLUSIONS

Our results suggest that hybridization between co-occurring red oaks occurs, but that selection may limit introgression, especially at functional loci. However, by providing a source of genetic variation, hybridization could influence the response of oaks and other hybridizing taxa to environmental change.

SNP Discovery through Next-Generation Sequencing and Its Applications

The decreasing cost along with rapid progress in next-generation sequencing and related bioinformatics computing resources has facilitated large-scale discovery of SNPs in various model and nonmodel plant species. Large numbers and genome-wide availability of SNPs make them the marker of choice in partially or completely sequenced genomes. Although excellent reviews have been published on next-generation sequencing, its associated bioinformatics challenges, and the applications of SNPs in genetic studies, a comprehensive review connecting these three intertwined research areas is needed. This paper touches upon various aspects of SNP discovery, highlighting key points in availability and selection of appropriate sequencing platforms, bioinformatics pipelines, SNP filtering criteria, and applications of SNPs in genetic analyses. The use of next-generation sequencing methodologies in many non-model crops leading to discovery and implementation of SNPs in various genetic studies is discussed. Development and improvement of bioinformatics software that are open source and freely available have accelerated the SNP discovery while reducing the associated cost. Key considerations for SNP filtering and associated pipelines are discussed in specific topics. A list of commonly used software and their sources is compiled for easy access and reference.

Fine-scale analysis reveals cryptic landscape genetic structure in desert tortoises

Characterizing the effects of landscape features on genetic variation is essential for understanding how landscapes shape patterns of gene flow and spatial genetic structure of populations. Most landscape genetics studies have focused on patterns of gene flow at a regional scale. However, the genetic structure of populations at a local scale may be influenced by a unique suite of landscape variables that have little bearing on connectivity patterns observed at broader spatial scales. We investigated fine-scale spatial patterns of genetic variation and gene flow in relation to features of the landscape in desert tortoise (Gopherus agassizii), using 859 tortoises genotyped at 16 microsatellite loci with associated data on geographic location, sex, elevation, slope, and soil type, and spatial relationship to putative barriers (power lines, roads). We used spatially explicit and non-explicit Bayesian clustering algorithms to partition the sample into discrete clusters, and characterize the relationships between genetic distance and ecological variables to identify factors with the greatest influence on gene flow at a local scale. Desert tortoises exhibit weak genetic structure at a local scale, and we identified two subpopulations across the study area. Although genetic differentiation between the subpopulations was low, our landscape genetic analysis identified both natural (slope) and anthropogenic (roads) landscape variables that have significantly influenced gene flow within this local population. We show that desert tortoise movements at a local scale are influenced by features of the landscape, and that these features are different than those that influence gene flow at larger scales. Our findings are important for desert tortoise conservation and management, particularly in light of recent translocation efforts in the region. More generally, our results indicate that recent landscape changes can affect gene flow at a local scale and that their effects can be detected almost immediately.

The colonization history of Juniperus brevifolia (Cupressaceae) in the Azores Islands

Background

A central aim of island biogeography is to understand the colonization history of insular species using current distributions, fossil records and genetic diversity. Here, we analyze five plastid DNA regions of the endangered Juniperus brevifolia, which is endemic to the Azores archipelago.

Methodology/Principal Findings

The phylogeny of the section Juniperus and the phylogeographic analyses of J. brevifolia based on the coalescence theory of allele (plastid) diversity suggest that: (1) a single introduction event likely occurred from Europe; (2) genetic diversification and inter-island dispersal postdated the emergence of the oldest island (Santa Maria, 8.12 Ma); (3) the genetic differentiation found in populations on the islands with higher age and smaller distance to the continent is significantly higher than that on the younger, more remote ones; (4) the high number of haplotypes observed (16), and the widespread distribution of the most frequent and ancestral ones across the archipelago, are indicating early diversification, demographic expansion, and recurrent dispersal. In contrast, restriction of six of the seven derived haplotypes to single islands is construed as reflecting significant isolation time prior to colonization.

Conclusions/Significance

Our phylogeographic reconstruction points to the sequence of island emergence as the key factor to explain the distribution of plastid DNA variation. The reproductive traits of this juniper species (anemophily, ornithochory, multi-seeded cones), together with its broad ecological range, appear to be largely responsible for recurrent inter-island colonization of ancestral haplotypes. In contrast, certain delay in colonization of new haplotypes may reflect intraspecific habitat competition on islands where this juniper was already present.

Prisoners in their habitat? Generalist dispersal by habitat specialists: a case study in southern water vole (Arvicola sapidus)

Habitat specialists inhabiting scarce and scattered habitat patches pose interesting questions related to dispersal such as how specialized terrestrial mammals do to colonize distant patches crossing hostile matrices. We assess dispersal patterns of the southern water vole (Arvicola sapidus), a habitat specialist whose habitat patches are distributed through less than 2% of the study area (overall 600 km²) and whose populations form a dynamic metapopulational network. We predict that individuals will require a high ability to move through the inhospitable matrix in order to avoid genetic and demographic isolations. Genotypes (N = 142) for 10 microsatellites and sequences of the whole mitochondrial Control Region (N = 47) from seven localities revealed a weak but significant genetic structure partially explained by geographic distance. None of the landscape models had a significant effect on genetic structure over that of the Euclidean distance alone and no evidence for efficient barriers to dispersal was found. Contemporary gene flow was not severely limited for A. sapidus as shown by high migration rates estimates (>10%) between non-neighbouring areas. Sex-biased dispersal tests did not support differences in dispersal rates, as shown by similar average axial parent-offspring distances, in close agreement with capture-mark-recapture estimates. As predicted, our results do not support any preferences of the species for specific landscape attributes on their dispersal pathways. Here, we combine field and molecular data to illustrate how a habitat specialist mammal might disperse like a habitat generalist, acquiring specific long-distance dispersal strategies as an adaptation to patchy, naturally fragmented, heterogeneous and unstable habitats.

Genetic structure and the North American postglacial expansion of the barnacle, Semibalanus balanoides

Population genetic characteristics are shaped by the life-history traits of organisms and the geologic history of their habitat. This study provides a neutral framework for understanding the population dynamics and opportunities for selection in Semibalanus balanoides, a species that figures prominently in ecological and evolutionary studies in the Atlantic intertidal. We used mitochondrial DNA (mtDNA) control region (N = 131) and microsatellite markers (∼40 individuals/site/locus) to survey populations of the broadly dispersing acorn barnacle from 8 sites spanning 800 km of North American coast and 1 site in Europe. Patterns of mtDNA sequence evolution were consistent with larger population sizes in Europe and population expansion at the conclusion of the last ice age, approximately 20 000 years ago, in North America. A significant portion of mitochondrial diversity was partitioned between the continents (ϕ(ST) = 0.281), but there was only weak structure observed from mtDNA within North America. Microsatellites showed significant structuring between the continents (F(ST) = 0.021) as well as within North America (F(ST) = 0.013). Isolation by distance in North America was largely driven by a split between populations south of Cape Cod and all others (P < 10(-4)). The glacial events responsible for generating allelic diversity at mtDNA and microsatellites may also be responsible for generating selectable variation at metabolic enzymes in S. balanoides.

Major Histocompatibility Complex (MHC) markers in conservation biology

Human impacts through habitat destruction, introduction of invasive species and climate change are increasing the number of species threatened with extinction. Decreases in population size simultaneously lead to reductions in genetic diversity, ultimately reducing the ability of populations to adapt to a changing environment. In this way, loss of genetic polymorphism is linked with extinction risk. Recent advances in sequencing technologies mean that obtaining measures of genetic diversity at functionally important genes is within reach for conservation programs. A key region of the genome that should be targeted for population genetic studies is the Major Histocompatibility Complex (MHC). MHC genes, found in all jawed vertebrates, are the most polymorphic genes in vertebrate genomes. They play key roles in immune function via immune-recognition and -surveillance and host-parasite interaction. Therefore, measuring levels of polymorphism at these genes can provide indirect measures of the immunological fitness of populations. The MHC has also been linked with mate-choice and pregnancy outcomes and has application for improving mating success in captive breeding programs. The recent discovery that genetic diversity at MHC genes may protect against the spread of contagious cancers provides an added impetus for managing and protecting MHC diversity in wild populations. Here we review the field and focus on the successful applications of MHC-typing for conservation management. We emphasize the importance of using MHC markers when planning and executing wildlife rescue and conservation programs but stress that this should not be done to the detriment of genome-wide diversity.

Demographic processes underlying subtle patterns of population structure in the scalloped hammerhead shark, Sphyrna lewini

Genetic diversity (θ), effective population size (N(e)), and contemporary levels of gene flow are important parameters to estimate for species of conservation concern, such as the globally endangered scalloped hammerhead shark, Sphyrna lewini. Therefore, we have reconstructed the demographic history of S. lewini across its Eastern Pacific (EP) range by applying classical and coalescent population genetic methods to a combination of 15 microsatellite loci and mtDNA control region sequences. In addition to significant population genetic structure and isolation-by-distance among seven coastal sites between central Mexico and Ecuador, the analyses revealed that all populations have experienced a bottleneck and that all current values of θ are at least an order of magnitude smaller than ancestral θ, indicating large decreases in N(e) (θ = 4N(e)μ), where μ is the mutation rate. Application of the isolation-with-migration (IM) model showed modest but significant genetic connectivity between most sampled sites (point estimates of Nm = 0.1-16.7), with divergence times (t) among all populations significantly greater than zero. Using a conservative (i.e., slow) fossil-based taxon-specific phylogenetic calibration for mtDNA mutation rates, posterior probability distributions (PPDs) for the onset of the decline in N(e) predate modern fishing in this region. The cause of decline over the last several thousand years is unknown but is highly atypical as a post-glacial demographic history. Regardless of the cause, our data and analyses suggest that S. lewini was far more abundant throughout the EP in the past than at present.