Islands within islands: two montane palaeo‐endemic birds impacted by recent anthropogenic fragmentation

Anthropogenic habitat fragmentation of species that live in naturally patchy metapopulations such as mountaintops or sky islands experiences two levels of patchiness. Effects of such multilevel patchiness on species have rarely been examined. Metapopulation theory suggests that patchy habitats could have varied impacts on persistence, dependent on differential migration. It is not known whether montane endemic species, evolutionarily adapted to natural patchiness, are able to disperse between anthropogenic fragments at similar spatial scales as natural patches. We investigated historic and contemporary gene flow between natural and anthropogenic patches across the distribution range of a Western Ghats sky‐island‐endemic bird species complex. Data from 14 microsatellites for 218 individuals detected major genetic structuring by deep valleys, including one hitherto undescribed barrier. As expected, we found strong effects of historic genetic differentiation across natural patches, but not across anthropogenic fragments. Contrastingly, contemporary differentiation (D_PS) was higher relative to historic differentiation (F_ST) in anthropogenic fragments, despite the species’ ability to historically traverse shallow valleys. Simulations of recent isolation resulted in high D_PS/F_ST values, confirming recent isolation in Western Ghats anthropogenic fragments and also suggesting that this ratio can be used to identifying recent fragmentation in the context of historic connectedness. We suggest that in this landscape, in addition to natural patchiness affecting population connectivity, anthropogenic fragmentation additionally impacts connectivity, making anthropogenic fragments akin to islands within natural islands of montane habitat, a pattern that may be recovered in other sky‐island systems.     

Relative performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation

Traditional methods for characterizing genetic differentiation among populations rely on a priori grouping of individuals. Bayesian clustering methods avoid this limitation by using linkage and Hardy–Weinberg disequilibrium to decompose a sample of individuals into genetically distinct groups. There are several software programs available for Bayesian clustering analyses, all of which describe a decrease in the ability to detect distinct clusters as levels of genetic differentiation among populations decrease. However, no study has yet compared the performance of such methods at low levels of population differentiation, which may be common in species where populations have experienced recent separation or high levels of gene flow. We used simulated data to evaluate the performance of three Bayesian clustering software programs, PARTITION, STRUCTURE, and BAPS, at levels of population differentiation below F _ST=0.1. PARTITION was unable to correctly identify the number of subpopulations until levels of F _ST reached around 0.09. Both STRUCTURE and BAPS performed very well at low levels of population differentiation, and were able to correctly identify the number of subpopulations at F _ST around 0.03. The average proportion of an individual’s genome assigned to its true population of origin increased with increasing F _ST for both programs, reaching over 92% at an F _ST of 0.05. The average number of misassignments (assignments to the incorrect subpopulation) continued to decrease as F _ST increased, and when F _ST was 0.05, fewer than 3% of individuals were misassigned using either program. Both STRUCTURE and BAPS worked extremely well for inferring the number of clusters when clusters were not well-differentiated (F _ST=0.02–0.03), but our results suggest that F _ST must be at least 0.05 to reach an assignment accuracy of greater than 97%.     

Comparison of genetic differentiation at marker loci and quantitative traits

The comparison of the degree of differentiation in neutral marker loci and genes coding quantitative traits with standardized and equivalent measures of genetic differentiation (F_ST and Q_ST, respectively) can provide insights into two important but seldom explored questions in evolutionary genetics: (i) what is the relative importance of random genetic drift and directional natural selection as causes of population differentiation in quantitative traits, and (ii) does the degree of divergence in neutral marker loci predict the degree of divergence in genes coding quantitative traits? Examination of data from 18 independent studies of plants and animals using both standard statistical and meta‐analytical methods revealed a number of interesting points. First, the degree of differentiation in quantitative traits (Q_ST) typically exceeds that observed in neutral marker genes (F_ST), suggesting a prominent role for natural selection in accounting for patterns of quantitative trait differentiation among contemporary populations. Second, the F_ST – Q_ST difference is more pronounced for allozyme markers and morphological traits, than for other kinds of molecular markers and life‐history traits. Third, very few studies reveal situations were Q_ST < F_ST, suggesting that selection pressures, and hence optimal phenotypes, in different populations of the same species are unlikely to be often similar. Fourth, there is a strong correlation between Q_ST and F_ST indices across the different studies for allozyme (r=0.81), microsatellite (r=0.87) and combined (r=0.75) marker data, suggesting that the degree of genetic differentiation in neutral marker loci is closely predictive of the degree of differentiation in loci coding quantitative traits. However, these interpretations are subject to a number of assumptions about the data and methods used to derive the estimates of population differentiation in the two sets of traits.     

GENETIC DIFFERENTIATION AMONG CHEWING LOUSE POPULATIONS (MALLOPHAGA: TRICHODECTIDAE) IN A POCKET GOPHER CONTACT ZONE (RODENTIA: GEOMYIDAE)

Genetic variation among populations of chewing lice (Geomydoecus actuosi) was examined in relation to chromosomal and electrophoretic variation among populations of their hosts (Thomomys bottae) at a contact zone. Louse demes were characterized by low levels of genetic heterozygosity (H? = 0.039) that may result from founder effects during primary infestation of hosts, compounded by seasonal reductions in louse population size. Louse populations sampled from different hosts showed high levels of genetic structuring both within and among host localities. Microgeographic differentiation of louse populations is high (mean F_ST = 0.092) suggesting that properties of this host–parasite system promote differentiation of louse populations living on different individual hosts. Among-population differentiation in lice (F_ST = 0.240) was similar to that measured among host populations (F_ST = 0.236), suggesting a close association between gene flow in pocket gophers and gene flow in their lice.     

Estimation,variance and optimal sampling of gene diversity II. Diploid locus

Nei's analysis of diversity at a diploid locus is extended to a population subdivided into a large number of subpopulations. The diversities and the heterozygotes frequency are defined with respect to the total population and unbiasedly estimated in a two-stage random cluster sampling. The fixation indices F _IS, F _IT andF _ST are derived, then inter- and intra-population variances of the estimated parameters are studied. We show that there is a unique sample size per population which yields the best accuracy in estimatingF _ST and F _IS, respectively, at a given locus. These results are illustrated with an analysis of DNA diversity in a forest tree and compared to those obtained under the Hardy-Weinberg assumption.     

A new method for the partition of allelic diversity within and between subpopulations

A method is proposed for the analysis of allelic diversity in the context of subdivided populations. The definition of an allelic distance between subpopulations allows for the partition of total allelic diversity into within- and between-subpopulation components, in a way analogous to the classical partition of gene diversity. A new definition of allelic differentiation, A _ST , between subpopulations results from this partition, and is contrasted with the concept of allelic richness differentiation. The partition of allelic diversity makes it possible to establish the relative contribution of each subpopulation to within and between-subpopulation components of diversity with implications in priorisation for conservation. A comparison between this partition and that corresponding to allelic richness is illustrated with an example. Computer simulations are used to investigate the behaviour of the new statistic A _ST in comparison with F _ST for a finite island model under a range of mutation and migration rates. A _ST has less dependence on migration rate than F _ST for large values of migration rate, but the opposite occurs for low migration rates. In addition, the variance in the estimates of A _ST is higher than that of F _ST for low mutation rates, but the opposite for high mutation rates.     

PHASE III OF WRIGHT'S SHIFTING BALANCE PROCESS AND THE VARIANCE AMONG DEMES IN MIGRATION RATE

Interdemic selection by the differential migration of individuals out from demes of high fitness and into demes of low fitness (Phase III) is one of the most controversial aspects of Wright's Shifting Balance Theory. I derive a relationship between Phase III migration and the interdemic selection differential, S, and show its potential effect on F_ST. The relationship reveals a diversifying effect of interdemic selection by Phase III migration on the genetic structure of a metapopulation. Using experimental metapopulations, I explored the effect of Phase III migration on F_ST by comparing the genetic variance among demes for two different patterns of migration: (1) island model migration and (2) Wright's Phase III migration. Although mean migration rates were the same, I found that the variance among demes in migration rate was significantly higher with Phase III than with island model migration. As a result, F_ST for the frequency of a neutral marker locus was higher with Phase III than it was with island model migration. By increasing F_ST, Phase III enhanced the genetic differentiation among demes for traits not subject to interdemic selection. This feature makes Wright's process different from individual selection which, by reducing effective population size, decreases the genetic variance within demes for all other traits. I discussed this finding in relation to the efficacy of Phase III and random migration for effecting peak shifts, and the contribution of genes with indirect effects to among‐deme variation.     

gesp: A computer program for modelling genetic effective population size,inbreeding and divergence in substructured populations

The genetically effective population size (N_e) is of key importance for quantifying rates of inbreeding and genetic drift and is often used in conservation management to set targets for genetic viability. The concept was developed for single, isolated populations and the mathematical means for analysing the expected N_e in complex, subdivided populations have previously not been available. We recently developed such analytical theory and central parts of that work have now been incorporated into a freely available software tool presented here. gesp (Genetic Effective population size, inbreeding and divergence in Substructured Populations) is R‐based and designed to model short‐ and long‐term patterns of genetic differentiation and effective population size of subdivided populations. The algorithms performed by gesp allow exact computation of global and local inbreeding and eigenvalue effective population size, predictions of genetic divergence among populations (G_ST) as well as departures from random mating (F_IS, F_IT) while varying (i) subpopulation census and effective size, separately or including trend of the global population size, (ii) rate and direction of migration between all pairs of subpopulations, (iii) degree of relatedness and divergence among subpopulations, (iv) ploidy (haploid or diploid) and (v) degree of selfing. Here, we describe gesp and exemplify its use in conservation genetics modelling.     

Genetic hitchhiking in a subdivided population of <Emphasis Type="Italic">Mytilus edulis</Emphasis>

Background  

Few models of genetic hitchhiking in subdivided populations have been developed and the rarity of empirical examples is even more striking. We here provide evidences of genetic hitchhiking in a subdivided population of the marine mussel Mytilus edulis. In the Bay of Biscay (France), a patch of M. edulis populations happens to be separated from its North Sea conspecifics by a wide region occupied only by the sister species M. galloprovincialis. Although genetic differentiation between the two M. edulis regions is largely non-significant at ten marker loci (average F_ST~0.007), a strong genetic differentiation is observed at a single locus (F_ST = 0.25). We validated the outlier status of this locus, and analysed DNA sequence polymorphism in order to identify the nature of the selection responsible for the unusual differentiation.     

Behavioral responses of Ceratitis capitata flies to bait spray droplets and natural food

We usedF-statistics to quantify the population structure of two sympatric species of leaf beetles,Oreina cacaliae andO. speciosissima (Chrysomelinae, Coleoptera), which share the same microhabitat since they feed on the same herbaceous host plants. We measured genetic differentiation at six allozyme loci 1) among populations separated by relatively small distances (40 to 250 kilometers), 2) within each population, and 3) between sexes within populations. For both species, the populations studied are not panmictic. For each population, heterozygosities are relatively high, but the observed heterozygosities are generally lower than the expected values. Overall, within-population differentiation is high and similar for both species (F _is=0.326 forO. cacaliae and 0.332 forO. speciosissima). Additionally, populations of both species are highly differentiated (F _st=0.234 versus 0.051 forO. speciosissima). ForO. cacaliae,F _is andF _st are greater among females than among males, while forO. speciosissima,F _st is sustantially greater among the males whileF _is is slightly greater among males. Differences in gene frequency among the sexes were statistically tested using a modifiedF _st with sex as the defining category, and the sexes differed significantly with the exception of one population inO. cacaliae. Possible explanations for this difference are discussed.     

Estimation,variance and optimal sampling of gene diversity

An extension of Nei's analysis of diversity in a subdivided population is proposed for a haploid locus. The differentiation G _STbecomes a natural extension of Wright's F _STand generalizes Weir and Cockerham's parameter of co-ancestry by relaxing the assumption of identical correlation for all the alleles. Inter- and intrapopulation variances of the estimated diversities and differentiation are derived. Finally, the optimal sampling strategy for measuring G _STwhen a fixed number of individuals can be analysed is considered. It is shown that, at a given locus, there is a unique sample size per population which yields the smallest variance of G _ST,regardless of the number of populations studied. These theoretical developments are illustrated with an analysis of chloroplast DNA diversity in a forest tree. The results emphasize the necessity of sampling many populations, rather than many individuals per population, for an accurate measurement of the subdivision of gene diversity at a single locus.     

Contrasting patterns of morphological and neutral genetic divergence among geographically proximate populations of sockeye salmon Oncorhynchus nerka in Lake Aleknagik,Alaska

Levels of differentiation in morphological traits (age at maturity, body length at age, egg mass and body depth) and spawning time were examined in sockeye salmon Oncorhynchus nerka from three geographically proximate but physically distinct creeks in Lake Aleknagik, Alaska. Happy Creek fish had significantly greater values for most measured morphological traits, and Eagle Creek fish spawned significantly later than fish in the other creeks. Phenotypic differentiation between creeks, measured as P_ST, was then compared with microsatellite marker differentiation between creeks, measured as F_ST. No correlations were apparent between P_ST and F_ST values, and P_ST values were generally significantly larger than zero (P_ST= 0·0018–0·31) whereas F_ST values were not (F_ST=?0·0004 to 0·0016). The insignificant pair‐wise F_ST values between creek samples indicated that gene flow occurs between creeks, assuming the creek populations have reached migration–drift equilibrium. However, the strong homing behaviour of sockeye salmon precludes a scenario in which fish from the three creeks constitute a single population that segregates by body size. Rather, significant phenotypic differentiation suggests that strong divergent selection occurs on the phenotypic traits despite the homogenizing effects of gene flow.     

Generalization of the QST framework in hierarchically structured populations: Impacts of inbreeding and dominance

Q_ST is a differentiation parameter based on the decomposition of the genetic variance of a trait. In the case of additive inheritance and absence of selection, it is analogous to the genic differentiation measured on individual loci, F_ST. Thus, Q_ST?F_ST comparison is used to infer selection: selective divergence when Q_ST > F_ST, or convergence when Q_ST < F_ST. The definition of Q‐statistics was extended to two‐level hierarchical population structures with Hardy–Weinberg equilibrium. Here, we generalize the Q‐statistics framework to any hierarchical population structure. First, we developed the analytical definition of hierarchical Q‐statistics for populations not at Hardy–Weinberg equilibrium. We show that the Q‐statistics values obtained with the Hardy–Weinberg definition are lower than their corresponding F‐statistics when F_IS > 0 (higher when F_IS < 0). Then, we used an island model simulation approach to investigate the impact of inbreeding and dominance on the Q_ST?F_ST framework in a hierarchical population structure. We show that, while differentiation at the lower hierarchical level (Q_SR) is a monotonic function of migration, differentiation at the upper level (Q_RT) is not. In the case of additive inheritance, we show that inbreeding inflates the variance of Q_RT, which can increase the frequency of Q_RT > F_RT cases. We also show that dominance drastically reduces Q‐statistics below F‐statistics for any level of the hierarchy. Therefore, high values of Q‐statistics are good indicators of selection, but low values are not in the case of dominance.     

Mitochondrial genetic diversity and gene flow of common carp from main river drainages in China

1. The common carp (Cyprinus carpio L.) is one of the most widely distributed and important freshwater fishes in the world. In China, the common carp has been recognised as three subspecies: C. carpio haematopterus, C. carpio rubrofuscus and C. carpio carpio. The Nanling Mountains have been suggested as providing the dividing line between C. carpio haematopterus and C. carpio rubrofuscus. However, the demographic history and gene flow of the common carp in China is not clear. 2. We collected mitochondrial COII and D‐loop sequences (1494 bp) from 241 individuals distributed in eight major river drainages across China. The objective was to provide the first investigation into population genetic structure, demographic history and migration patterns of the common carp from these river drainages, and to assess the validity of the three subspecies. 3. Phylogenetic analysis did not result in three major monophyletic lineages corresponding to the three subspecies. The Nanling Mountains do not form a border separating C. carpio haematopterus and C. carpio rubrofuscus. amova showed low population differentiation, with 11.60% of the molecular variance found among river drainages. Pairwise F_ST values between river drainages were moderate (0.0331–0.2617). Substantial gene flow detected in coalescent analysis between drainages showed that the Yangtze drainage was the centre from which migrants moved northward, southward and north‐westward. Human‐mediated translocation has confounded our ability to identify subspecies of common carp in China.     

A Beta-mixture model for assessing genetic population structure

Summary An important fraction of recently generated molecular data is dominant markers. They contain substantial information about genetic variation but dominance makes it impossible to apply standard techniques to calculate measures of genetic differentiation, such as F‐statistics. In this article, we propose a new Bayesian beta‐mixture model that more accurately describes the genetic structure from dominant markers and estimates multiple F_ST s from the sample. The model also has important application for codominant markers and single‐nucleotide polymorphism (SNP) data. The number of F_ST is assumed unknown beforehand and follows a random distribution. The reversible jump algorithm is used to estimate the unknown number of multiple F_ST s. We evaluate the performance of three split proposals and the overall performance of the proposed model based on simulated dominant marker data. The model could reliably identify and estimate a spectrum of degrees of genetic differentiation present in multiple loci. The estimates of F_ST s also incorporate uncertainty about the magnitude of within‐population inbreeding coefficient. We illustrate the method with two examples, one using dominant marker data from a rare orchid and the other using codominant marker data from human populations.     

ESTIMATES OF INTRAGAMETOPHYTIC SELFING AND INTERPOPULATIONAL GENE FLOW IN HOMOSPOROUS FERNS

Relatively little information is available on mating systems and interpopulational gene flow in species of homosporous pteridophytes. Because of the proximity of antheridia and archegonia on the same thallus, it has long been maintained that intragametophytic selling is the predominant mode of reproduction in natural populations of homosporous ferns and other homosporous plants. Furthermore, quantitative estimates of interpopulational gene flow via spore dispersal are lacking. In this paper, we examine five species of homosporous ferns (Botrychium virginianum, Polystichum munitum, P. imbricans, Blechnum spicant, and Dryopteris expansa) and present estimates of 1) rates of intragametophytic selling, 2) levels of interpopulational gene flow, and 3) interpopulational genetic differentiation (F_ST). Our data demonstrate that mating systems vary among species of ferns, just as they do among species of seed plants. The data also suggest that levels of interpopulational gene flow are generally high. The F_ST values indicate little genetic divergence among populations for all species except Dryopteris expansa, which exhibits significant levels of interpopulational genetic differentiation. Patterns of genetic diversity in the five species examined are related to the mating system and rate of interpopulational gene flow in each species. The F_ST values for all species except Botrychium virginianum are in close agreement with those predicted for an island model of population structure.     

VOCAL DIALECTS AND THE STRUCTURE OF GEOGRAPHIC VARIATION IN MORPHOLOGICAL AND ALLOZYMIC CHARACTERS IN THE RUFOUS-COLLARED SPARROW,ZONOTRICHIA CAPENSIS

The geographical patterns of variation shown at 20 allozyme and non-enzymatic protein-coding loci, in 8 external, and in 12 skeletal morphological characters in the rufous-collared sparrow, Zonotrichia capensis, were analyzed in order to test the local (genetic) adaptation hypothesis regarding the origin and maintenance of vocal dialects in birds. Approximately 20 males were collected from each of four sites within each of six different dialect zones. There was significant variability in both external and skeletal morphology among all 24 sites and among dialect groups. Average Wright's corrected fixation coefficient (F_ST) was 0.118, indicating significant genetic differentiation among all sites, regardless of dialect. Hierarchical F statistics indicated that only 50% of among site variability was due to a dialect effect. Puna dialect sites were highly differentiated from all other sites with respect to both morphology (external and skeletal measures) and allozyme frequencies. Heterogeneity at the PGM-1 locus among puna scrub sites was the major cause of the high average F_ST across all sites, and within the puna scrub dialect. Average genetic differentiation among non-puna sites (F_ST = 0.018) was similar to differentiation among sites within each of the five non-puna dialect groups (mean F_ST = 0.0132 ± 0.0069). Hierarchical F statistics indicated that none of the among-site differentiation in this subset of samples was due to a dialect effect. These observations are not consistent with the local adaptation hypothesis. All significant genetic heterogeneity occurred among sites in mountainous habitats, and we suggest that topography and patchiness of habitat may have been major factors involved in population differentiation, rather than vocal dialects.     

Background selection and FST: Consequences for detecting local adaptation

Background selection is a process whereby recurrent deleterious mutations cause a decrease in the effective population size and genetic diversity at linked loci. Several authors have suggested that variation in the intensity of background selection could cause variation in F_ST across the genome, which could confound signals of local adaptation in genome scans. We performed realistic simulations of DNA sequences, using recombination maps from humans and sticklebacks, to investigate how variation in the intensity of background selection affects F_ST and other statistics of population differentiation in sexual, outcrossing species. We show that, in populations connected by gene flow, Weir and Cockerham's (1984; Evolution, 38 , 1358) estimator of F_ST is largely insensitive to locus‐to‐locus variation in the intensity of background selection. Unlike F_ST, however, d_XY is negatively correlated with background selection. Moreover, background selection does not greatly affect the false‐positive rate in F_ST outlier studies in populations connected by gene flow. Overall, our study indicates that background selection will not greatly interfere with finding the variants responsible for local adaptation.     

Landscape genetics of alpine Sierra Nevada salamanders reveal extreme population subdivision in space and time

Quantifying the influence of the landscape on the genetic structure of natural populations remains an important empirical challenge, particularly for poorly studied, ecologically cryptic species. We conducted an extensive microsatellite analysis to examine the population genetics of the southern long‐toed salamander (Ambystoma macrodactylum sigillatum) in a naturally complex landscape. Using spatially explicit modelling, we investigated the influence of the Sierra Nevada topography on potential dispersal corridors between sampled populations. Our results indicate very high‐genetic divergence among populations, high within‐deme relatedness, and little evidence of recent migration or population admixture. We also discovered unexpectedly high between‐year genetic differentiation (F_ST) for breeding sites, suggesting that breeding groups vary over localized space and time. While environmental factors associated with high‐elevation montane habitats apparently play an important role in shaping population differentiation, additional, species‐specific biological processes must also be operating to account for observed deviations from temporal, among‐year panmixia. Our study emphasizes the population‐level insights that can be gained from high‐density sampling in space and time, and the highly substructured population biology that may characterize amphibians in extreme montane habitats.     

Does GST underestimate genetic differentiation from marker data?

The widely applied genetic differentiation statistics F_ST and G_ST have recently been criticized for underestimating differentiation when applied to highly polymorphic markers such as microsatellites. New statistics claimed to be unaffected by marker polymorphisms have been proposed and advocated to replace the traditional F_ST and G_ST. This study shows that G_ST gives accurate estimates and underestimates of differentiation when demographic factors are more and less important than mutations, respectively. In the former case, all markers, regardless of diversity (H_S), have the same G_ST value in expectation and thus give replicated estimates of differentiation. In the latter case, markers of higher H_S have lower G_ST values, resulting in a negative, roughly linear correlation between G_ST and H_S across loci. I propose that the correlation coefficient between G_ST and H_S across loci, r_GH, can be used to distinguish the two cases and to detect mutational effects on G_ST. A highly negative and significant r_GH, when coupled with highly variable G_ST values among loci, would reveal that marker G_ST values are affected substantially by mutations and marker diversity, underestimate population differentiation, and are not comparable among studies, species and markers. Simulated and empirical data sets are used to check the power and statistical behaviour, and to demonstrate the usefulness of the correlation analysis.