Geographic variation in resource use by specialist versus generalist butterflyfishes

Localised patterns of resource use can be constrained by multiple factors. Comparison of resource use at multiple locations with differing resource availability can allow fundamental specialists to be distinguished from species that simply feed predominantly on prey types that are locally abundant. This study investigates geographic variation in the feeding ecology of coral‐feeding butterflyfishes to examine whether patterns of resource use and levels of dietary specialisation vary among distinct locations, corresponding with changes in resource availability. Our specific aims were to investigate whether the dietary niche breadth of four butterflyfishes varies among five geographically separated locations and assess whether each species utilises similar resources in each location. Resource availability and dietary composition of four butterflyfishes were quantified at three sites across each of five geographic locations throughout the Pacific. Niche breadth, niche overlap, and resource selection functions were calculated for each species at each site and compared among locations. Availability of dietary resources varied significantly among locations and sites. Chaetodon vagabundus, C. citrinellus and C. lunulatus had low levels of dietary specialisation and used different resources in each location. Chaetodon trifascialis had high levels of dietary specialisation and used the same few resources in each location. Our results indicate that relative levels of dietary specialisation among different butterflyfishes do hold at larger spatial scales, however, geographical variation in the dietary composition of all butterflyfishes indicates that prey availability has a fundamental influence on dietary composition. Highly specialised species such as C. trifascialis will be highly vulnerable to coral loss as they appear to be largely inflexible in their dietary composition. However, the increased feeding plasticity observed here for C. trifascialis suggests this species may have a greater capacity to respond to coral loss than previously assumed.     

Interspecific variation in distributions and diets of coral reef butterflyfishes (Teleostei: Chaetodontidae)

This study examined within‐reef distributions for 19 species of butterflyfishes (Chaetodontidae) at Lizard Island, northern Great Barrier Reef, Australia, and compared spatial patterns of abundance among species with contrasting diets. Spatial variation in abundance of butterflyfishes was most prominent among physiognomic reef zones mainly due to significant zonation of eight species, including four obligate hard‐coral feeders (Chaetodon trifascialis, Chaetodon baronessa, Chaetodon plebeius and Chaetodon lunulatus) and four generalist species (Chaetodon auriga, Chaetodon citrinellus, Chaetodon kleinii and Chaetodon rafflesi). Distributions of obligate hard‐coral feeders were closely associated with spatial variation in percentage cover of scleractinian corals, but no more restricted compared with facultative hard‐coral feeders or non‐coral feeders. Species with highest dietary specialization (C. trifascialis and C. baronessa), however, exhibited the most pronounced zonation patterns and were restricted to habitats with greatest abundance of their preferred prey. While there are conspicuous links between dietary specialization v. spatial patterns in abundance of butterflyfishes, it remains unclear whether dietary specialization is the cause or consequence of more restricted distributions.     

Declines in the abundance of Chaetodon butterflyfishes following extensive coral depletion

This study documented temporal variation in the abundance of butterflyfishes (Chaetodontidae) at Trunk Reef, on the central Great Barrier Reef, Australia, from May 2000 to March 2005. During this period, live coral cover declined by >90%, mostly due to severe coral bleaching. There were no short-term changes (within 4 months) in the abundance of butterflyfishes following initial declines in live coral cover. Surveys conducted in 2005, however, revealed significant declines in the abundance of Chaetodon baronessa , Chaetodon lunulatus , Chaetodon trifascialis, Chaetodon plebeius and Chaetodon rainfordi , all of which are obligate hard-coral feeders. In contrast, there was no significant change in the abundance of Chaetodon auriga , Chaetodon aureofasciatus , Chaetodon citrinellus , Chaetodon melannotus or Chaetodon vagabundus , which are much less reliant on scleractinian coral for food. Clearly, extensive coral depletion, such as that caused by severe coral bleaching, can have a major effect on the abundance of butterflyfishes. Specific responses of butterflyfishes varied according to their reliance on hard corals for food and their ability to utilize alternate prey types.     

Regional population dynamics define the local genetic structure in Sorbus torminalis

Recent changes in sylvicultural practices in Central Europe have created forests with closed canopies, and tree species preferring open and sunny forests have declined in area and abundance. This led to increased isolation of populations of many rare insect-pollinated, fleshy-fruited species with a naturally scattered distribution. To gain insight into the regional population dynamics of such species, we investigated the consequences of spatial isolation, population size and density on the genetic structure of Sorbus torminalis and simultaneously considered the relationship between fecundity and habitat quality. Genotype data for biparentally (ISSRs) and maternally inherited (cpDNA PCR-RFLPs) molecular markers were generated for 26 Swiss populations of S. torminalis. We applied analyses of molecular variance (amova) to both marker types and separated the relative contributions of pollen and seed dispersal to historical gene flow. amova detected significant differentiation among populations (Phi(ST ISSR) = 0.107; Phi(ST cpDNA) = 0.370) in both marker types. The relative rate of pollen to seed gene flow was low (r = 2.919) and significantly different from equality. Isolation by distance was weak within Eastern and Western Switzerland, although populations were substantially differentiated. Within-population molecular variance was not explained by population size, whereas habitat quality (openness) positively influenced the percentage of fruiting trees and the degree of fruiting per tree, indicating that more open forests enhance sexual reproduction. Our findings of significant genetic differentiation in the absence of clear geographical structuring can be explained by the distinct ecology of S. torminalis and nondirectional colonization events in metapopulation-like dynamics.     

Microsatellite analysis reveals interpopulation differentiation and gene flow in the endangered tree Changiostyrax dolichocarpa (Styracaceae) with fragmented distribution in central China

Polymorphic simple sequence repeat (SSR) markers were used to investigate the impact of habitat fragmentation on the population structure and gene flow of Changiostyrax dolichocarpa, a critically endangered tree in central China. Intrapopulation genetic diversity, population structure and gene flow in the five extant populations of this species were analysed by eight SSR markers. Intrapopulation genetic diversity results suggest that C. dolichocarpa remnants maintained a relatively high degree of genetic diversity despite severe fragmentation. Low genetic differentiation among populations was found based on Wright's F(ST) and amova analysis. Both the F(ST)-based estimate and private allele method revealed high historical gene flow among the remnant populations. Recent immigrants, detected by assignment tests, tend to decrease from the grandparent generation to the current generation. The potentially highly restricted current gene flow among fragments may render the fragmented populations of C. dolichocarpa at a higher risk of local extinction several generations after fragmentation. Both in situ and ex situ conservation management for the remnant populations of C. dolichocarpa are therefore urgently needed to rescue remaining genetic diversity.     

Phylogeography of the Atlanto-Mediterranean sea cucumber Holothuria (Holothuria) mammata: the combined effects of historical processes and current oceanographical pattern

We assessed the genetic structure of populations of the widely distributed sea cucumber Holothuria (Holothuria) mammata Grube, 1840, and investigated the effects of marine barriers to gene flow and historical processes. Several potential genetic breaks were considered, which would separate the Atlantic and Mediterranean basins, the isolated Macaronesian Islands from the other locations analysed, and the Western Mediterranean and Aegean Sea (Eastern Mediterranean). We analysed mitochondrial 16S and COI gene sequences from 177 individuals from four Atlantic locations and four Mediterranean locations. Haplotype diversity was high (H=0.9307 for 16S and 0.9203 for COI), and the haplotypes were closely related (π=0.0058 for 16S and 0.0071 for COI). The lowest genetic diversities were found in the Aegean Sea population. Our results showed that the COI gene was more variable and more useful for the detection of population structure than the 16S gene. The distribution of mtDNA haplotypes, the pairwise F(ST) values and the results of exact tests and amova revealed: (i) a significant genetic break between the population in the Aegean Sea and those in the other locations, as supported by both mitochondrial genes, and (ii) weak differentiation of the Canary and Azores Islands from the other populations; however, the populations from the Macaronesian Islands, Algarve and West Mediterranean could be considered to be a panmictic metapopulation. Isolation by distance was not identified in H. (H.) mammata. Historical events behind the observed findings, together with the current oceanographic patterns, were proposed and discussed as the main factors that determine the population structure and genetic signature of H. (H.) mammata.     

Genetic structure and gene flow in a metapopulation of an endangered plant species,Silene tatarica

We investigated the distribution of genetic variation within and between seven subpopulations in a riparian population of Silene tatarica in northern Finland by using amplified fragment length polymorphism (AFLP) markers. A Bayesian approach-based clustering program indicated that the marker data contained not only one panmictic population, but consisted of seven clusters, and that each original sample site seems to consist of a distinct subpopulation. A coalescent-based simulation approach shows recurrent gene flow between subpopulations. Relative high FST values indicated a clear subpopulation differentiation. However, amova analysis and UPGMA-dendrogram did not suggest any hierarchical regional structuring among the subpopulations. There was no correlation between geographical and genetic distances among the subpopulations, nor any correlation between the subpopulation census size and amount of genetic variation. Estimates of gene flow suggested a low level of gene flow between the subpopulations, and the assignment tests proposed a few long-distance bidirectional dispersal events between the subpopulations. No apparent difference was found in within-subpopulation genetic diversity among upper, middle and lower regions along the river. Relative high amounts of linkage disequilibrium at subpopulation level indicated recent population bottlenecks or admixture, and at metapopulation levels a high subpopulation turnover rate. The overall pattern of genetic variation within and between subpopulations also suggested a 'classical' metapopulation structure of the species suggested by the ecological surveys.     

Population genetic structure and long-distance dispersal among seabird populations: implications for colony persistence

Dramatic local population decline brought about by anthropogenic-driven change is an increasingly common threat to biodiversity. Seabird life history traits make them particularly vulnerable to such change; therefore, understanding population connectivity and dispersal dynamics is vital for successful management. Our study used a 357-base pair mitochondrial control region locus sequenced for 103 individuals and 18 nuclear microsatellite loci genotyped for 245 individuals to investigate population structure in the Atlantic and Pacific populations of the pelagic seabird, Leach's storm-petrel Oceanodroma leucorhoa leucorhoa. This species is under intense predation pressure at one regionally important colony on St Kilda, Scotland, where a disparity between population decline and predation rates hints at immigration from other large colonies. AMOVA, F(ST), Φ(ST) and Bayesian cluster analyses revealed no genetic structure among Atlantic colonies (Global Φ(ST) = -0.02 P > 0.05, Global F(ST) = 0.003, P > 0.05, STRUCTURE K = 1), consistent with either contemporary gene flow or strong historical association within the ocean basin. The Pacific and Atlantic populations are genetically distinct (Global Φ(ST) = 0.32 P < 0.0001, Global F(ST) = 0.04, P < 0.0001, STRUCTURE K = 2), but evidence for interocean exchange was found with individual exclusion/assignment and population coalescent analyses. These findings highlight the importance of conserving multiple colonies at a number of different sites and suggest that management of this seabird may be best viewed at an oceanic scale. Moreover, our study provides an illustration of how long-distance movement may ameliorate the potentially deleterious impacts of localized environmental change, although direct measures of dispersal are still required to better understand this process.     

Molecular and quantitative genetic differentiation across Europe in yellow dung flies

Relating geographic variation in quantitative traits to underlying population structure is crucial for understanding processes driving population differentiation, isolation and ultimately speciation. Our study represents a comprehensive population genetic survey of the yellow dung fly Scathophaga stercoraria, an important model organism for evolutionary and ecological studies, over a broad geographic scale across Europe (10 populations from the Swiss Alps to Iceland). We simultaneously assessed differentiation in five quantitative traits (body size, development time, growth rate, proportion of diapausing individuals and duration of diapause), to compare differentiation in neutral marker loci (F(ST)) to that of quantitative traits (Q(ST)). Despite long distances and uninhabitable areas between sampled populations, population structuring was very low but significant (F(ST) = 0.007, 13 microsatellite markers; F(ST) = 0.012, three allozyme markers; F(ST) = 0.007, markers combined). However, only two populations (Iceland and Sweden) showed significant allelic differentiation to all other populations. We estimated high levels of gene flow [effective number of migrants (Nm) = 6.2], there was no isolation by distance, and no indication of past genetic bottlenecks (i.e. founder events) and associated loss of genetic diversity in any northern or island population. In contrast to the low population structure, quantitative traits were strongly genetically differentiated among populations, following latitudinal clines, suggesting that selection is responsible for life history differentiation in yellow dung flies across Europe.     

Contrasting microsatellite variation between subalpine and western larch, two closely related species with different distribution patterns

Subalpine larch (Larix lyallii Parl.) and western larch (Larix occidentalis Nutt.) represent two closely related species with contrasting abundance and distribution patterns in Western North America. Genetic diversity at seven informative microsatellite loci was determined for 19 populations of subalpine larch and nine populations of western larch. Contrasting genetic diversity and patterns of population differentiation were observed between the two species. The overall within-population genetic diversity parameters were lower in subalpine larch (A = 3.2; A(P) = 3.6; H(E) = 0.418) than in western larch (A(P) = 5.51; H(E) = 0.580), a pattern that is likely related to historical or demographic factors. No evidence of interspecific hybridization was observed. Significantly more population differentiation (theta = 0.15; R(ST) = 0.07), consistent with more restricted gene flow, was observed for subalpine larch as compared to western larch (theta = 0.05; R(ST) = 0.04). Under the assumption of an infinite allele mutation model, 12 of the 19 subalpine larch populations showed signs of deviation from the mutation-drift equilibrium, which suggests Holocene population bottlenecks and fluctuations in effective population size for this species. None of the western larch populations deviated significantly from the mutation-drift equilibrium. For both species, Mantel's test revealed a significant positive relationship between geographical and genetic distances indicative of isolation by distance. A similar geographical structure was detected in both species, suggesting at least two genetically distinct glacial populations in each species. The various implications for gene conservation are discussed.     

Genetic population structure of grey mackerel Scomberomorus semifasciatus in northern Australia

This study used mtDNA sequence and microsatellite markers to elucidate the population structure of Scomberomorus semifasciatus collected from 12 widespread sampling locations in Australia. Samples (n = 544) were genotyped with nine microsatellite loci, and 353 were sequenced for the control (384 bp) and ATPase (800 bp) mtDNA gene regions. Combined interpretation of microsatellite and mtDNA data identified four genetic stocks of S. semifasciatus: Western Australia, north-west coast of the Northern Territory, Gulf of Carpentaria and the eastern coast of Queensland. Connectivity among stocks across northern Australia from the Northern Territory to the eastern coast of Queensland was high (mean F(ST) = 0·003 for the microsatellite data and Φ(ST) = 0·033 and 0·009 for control region and ATPase, respectively) leading to some uncertainty about stock boundaries. In contrast, there was a clear genetic break between the stock in Western Australia compared to the rest of northern Australia (mean F(ST) = 0·132 for the microsatellite data and Φ(ST) = 0·135 and 0·188 for control region and ATPase, respectively). This indicates a restriction to gene flow possibly associated with suboptimal habitat along the Kimberley coast (north Western Australia). The appropriate scale of management for this species corresponds to the jurisdictions of the three Australian states, except that authorities in Queensland and Northern Territory should co-ordinate the management of the Gulf of Carpentaria stock.     

Influence of dietary specialization and resource availability on geographical variation in abundance of butterflyfish

Empirical evidence indicates that both niche breadth and resource availability are key drivers of a species' local abundance patterns. However, most studies have considered the influence of either niche breath or resource availability in isolation, while it is the interactive effects that are likely to influence local abundance. We examined geographic variation in the feeding ecology and distribution of coral-feeding butterflyfish to determine the influence of dietary specialization and dietary resource availability on their local abundance. Dietary composition and abundance of five butterflyfish and coral dietary resource availability were determined at 45 sites across five locations (Lizard Island and Heron Island, Great Barrier Reef; Kimbe Bay, Papua New Guinea; Noumea, New Caledonia; and Moorea, French Polynesia). Multiple regression models using variables representative of total dietary resource availability, availability of specific dietary resources, and interspecific competition were used to determine the best predictors of local abundance across all sites and locations for each species. Factors influencing local abundance varied between butterflyfish with specialized and generalized diets. Dietary resource availability had the strongest influence on the abundance of Chaetodon trifascialis-the most specialized species. Local abundance of C. trifascialis was best predicted by availability of the Acropora corals that it preferentially feeds on. In contrast, abundance of generalist butterflyfish was poorly described by variation in availability of specific resources. Rather, indices of total dietary resource availability best predicted their abundance. Overall, multiple regression models only explained a small proportion of the variation in local abundance for all five species. Despite their relatively specialized diets, dietary resource availability has limited influence on the local abundance of butterflyfish. Only the most specialized species appear to be consistently limited by prey availability. Local and total abundance of species are influenced by a wide range of different factors and there is definite need to conduct independent species assessments.     

Genetic connectivity among swarming sites in the wide ranging and recently declining little brown bat (Myotis lucifugus)

Characterizing movement dynamics and spatial aspects of gene flow within a species permits inference on population structuring. As patterns of structuring are products of historical and current demographics and gene flow, assessment of structure through time can yield an understanding of evolutionary dynamics acting on populations that are necessary to inform management. Recent dramatic population declines in hibernating bats in eastern North America from white‐nose syndrome have prompted the need for information on movement dynamics for multiple bat species. We characterized population genetic structure of the little brown bat, Myotis lucifugus, at swarming sites in southeastern Canada using 9 nuclear microsatellites and a 292‐bp region of the mitochondrial genome. Analyses of F_ST, Φ_ST, and Bayesian clustering (STRUCTURE) found weak levels of genetic structure among swarming sites for the nuclear and mitochondrial genome (Global F_ST = 0.001, P < 0.05, Global Φ_ST = 0.045, P < 0.01, STRUCTURE K = 1) suggesting high contemporary gene flow. Hierarchical AMOVA also suggests little structuring at a regional (provincial) level. Metrics of nuclear genetic structure were not found to differ between males and females suggesting weak asymmetries in gene flow between the sexes. However, a greater degree of mitochondrial structuring does support male‐biased dispersal long term. Demographic analyses were consistent with past population growth and suggest a population expansion occurred from approximately 1250 to 12,500 BP, following Pleistocene deglaciation in the region. Our study suggests high gene flow and thus a high degree of connectivity among bats that visit swarming sites whereby mainland areas of the region may be best considered as one large gene pool for management and conservation.     

Hierarchical components of genetic variation at a species boundary: population structure in two sympatric varieties of Lupinus microcarpus (Leguminosae)

Lupinus microcarpus is a self-compatible annual plant that forms a species complex of morphologically variable but indeterminate varieties. In order to examine the hypothesis that varieties of L. microcarpus comprise genetically differentiated and reproductively isolated species, populations of L. microcarpus var. horizontalis and var. densiflorus were sampled from an area of sympatry in central California and genotyped using six microsatellite loci. Bayesian clustering divided the total sample into two groups corresponding to the named varieties with extremely low levels of inferred coancestry. Similarly, maximum likelihood and distance methods for genetic assignment placed individuals in two nonoverlapping groups. Evidence for isolation by distance (IBD) within each variety was found at shorter distance classes, but varieties remained differentiated in sympatry. Furthermore, coalescent estimates of divergence time indicate separation within the past 950-5050 generations, with minimal gene flow after divergence. A four-level hierarchical analysis of molecular variance (amova) found significant levels of genetic differentiation among varieties (theta(P) = 0.292), populations within varieties (theta(S) = 0.449), subpopulations within populations (theta(SS) = 0.623), and individuals within subpopulations (f = 0.421); but the greatest degree of differentiation was at the subpopulation level. Although it is sometimes assumed that the magnitude of genetic differences (e.g. F(ST)) should be greater between species than among populations or subpopulations of the same species, shared ancestral polymorphism may lead to relatively low levels of differentiation at the species level, even as the stochastic effects of genetic drift generate higher levels of differentiation at lower hierarchical levels. These results suggest that L. microcarpus var. horizontalis and var. densiflorus are recently diverged yet reproductively isolated species, with high levels of inbreeding resulting from the combined effects of limited gene flow, demographic bottlenecks, and partial selfing in finite, geographically structured populations.     

Population structure and phylogeography of the short-tailed stingray, Dasyatis brevicaudata (Hutton 1875), in the Southern Hemisphere

There is accumulating evidence that the degree of vagility explains little of the extent of population subdivision found within elasmobranch species. Instead, patterns of gene flow in elasmobranchs appear more closely linked to the presence of dispersal barriers, either physical or biological. Here, we investigate the potential role of some of these isolating mechanisms in shaping the population structure of a member of the stingray family Dasyatidae (Dasyatis brevicaudata) at various scales (southern hemisphere vs. coastal New Zealand). Analyses of the mitochondrial DNA control region from 176 individuals revealed significant genetic structure between South Africa, Australia, and New Zealand populations (analysis of molecular variance [AMOVA], overall Ф(ST) = 0.67, P < 0.001), although New Zealand and Australia shared some haplotypes. Surprisingly, significant population differentiation was found among several coastal New Zealand locations (AMOVA, overall Ф(ST) = 0.05, P < 0.05). However, data did not support the genetic differentiation between individuals from an offshore breeding area and mainland individuals. Comparisons suggest that these stingrays exhibit similar levels of population differentiation as other coastal elasmobranchs, with high divergence across oceanic basins and lower differentiation along continuous coastal habitats. Differences in coastal population structuring in elasmobranch species studied to date may be attributed to species-specific preferences for coastal habitats, which may be linked to life history functions (e.g., feeding and pupping).     

Cryptic dietary components reduce dietary overlap among sympatric butterflyfishes (Chaetodontidae)

This study used three different methods to ascertain dietary composition for 21 Chaetodontidae species co-occurring on a single fringing reef in Derawan Island, East Kalimantan, Indonesia. The combination of in situ feeding observations, examination of gut contents and stable-isotope analyses was expected to identify previously unresolved prey items that contribute to dietary separation. In situ feeding observations identified five species that feed predominantly on anthozoans ( Chaetodon baronessa , Chaetodon bennetti , Chaetodon lunulatus , Chaetodon punctatofasciatus and Chaetodon speculum ). Stable-isotope ratios for these species, as well as for Chaetodon ornatissimus (for which no feeding observations were completed), were very similar and consistent with diets comprising mostly anthozoans. Feeding observations, however, showed that they mostly fed on different coral species, while the identifiable portion of their gut contents showed clear separation based on cryptic dietary components. For example, C. baronessa and C. bennetti appeared to ingest annelid worms during the course of coral feeding, whereas gut contents of C. punctatofasciatus and C. speculum were dominated by crustaceans. In situ feeding observations further identified the following groups: coral–bottom feeders, bottom feeders, sponge feeders and pelagic feeders, feeding on a wide variety of prey items such as Annelida, Crustacea, Cnidaria, Mollusca and macroalgae. Overall, many chaetodontid species had similar reliance on major prey items ( e.g . anthozoans or polychaetes) but differed greatly in the minor prey items that they utilized. Partitioning of minor prey items may be important in reducing interspecific competition and facilitating coexistence of chaetodontids on coral reefs.     

Historical and contemporary population genetic connectivity of the European short-snouted seahorse Hippocampus hippocampus and implications for management

This first genetic study of Hippocampus hippocampus covers the species' entire geographic range and employs two mtDNA markers (control region and cytochrome b) to establish patterns of population structuring. A total of 255 specimens from 21 locations were used to obtain 89 concatenated haplotypes. The common haplotype was present in all but one population, however, most haplotypes were unique. The haplotype network had a star-like construction, suggesting expansion from a bottleneck event. F(ST) and AMOVA revealed population subdivision into three geographic regions (English Channel + Bay of Biscay, Mediterranean Sea + Atlantic Ocean Iberian coast + Macaronesian Islands, and West Africa) with barriers to gene flow indentified at Cape Finisterre and the Cape Verde frontal zone. Neutrality tests and nested clade analysis suggest a complex demographic history, with both historic events and contemporary processes shaping patterns of genetic differentiation. The genetic population subdivision detected in this study indicates that H. hippocampus should be managed as three separate units. This is especially pertinent as H. hippocampus populations within the West African region are the only ones known to be specifically targeted for exploitation.     

Comparisons between QST and FST—how wrong have we been?

The comparison between quantitative genetic divergence (Q(ST) ) and neutral genetic divergence (F(ST) ) among populations has become the standard test for historical signatures of selection on quantitative traits. However, when the mutation rate of neutral markers is relatively high in comparison with gene flow, estimates of F(ST) will decrease, resulting in upwardly biased comparisons of Q(ST) vs. F(ST) . Reviewing empirical studies, the difference between Q(ST) and F(ST) is positively related to marker heterozygosity. After refuting alternative explanations for this pattern, we conclude that marker mutation rate indeed has had a biasing effect on published Q(ST) -F(ST) comparisons. Hence, it is no longer clear that populations have commonly diverged in response to divergent selection. We present and discuss potential solutions to this bias. Comparing Q(ST) with recent indices of neutral divergence that statistically correct for marker heterozygosity (Hedrick's G'st and Jost's D) is not advised, because these indices are not theoretically equivalent to Q(ST) . One valid solution is to estimate F(ST) from neutral markers with mutation rates comparable to those of the loci underlying quantitative traits (e.g. SNPs). Q(ST) can also be compared to Φ(ST) (Phi(ST) ) of amova, as long as the genetic distance among allelic variants used to estimate Φ(ST) reflects evolutionary history: in that case, neutral divergence is independent of mutation rate. In contrast to their common usage in comparisons of Q(ST) and F(ST) , microsatellites typically have high mutation rates and do not evolve according to a simple evolutionary model, so are best avoided in Q(ST) -F(ST) comparisons.     

Phylogeography of Varicorhinus barbatulus (Cyprinidae) in Taiwan based on nucleotide variation of mtDNA and allozymes

The phylogeographical patterns and population genetic structures of Varicorhinus barbatulus in Taiwan were investigated based on genetic diversity of 34 allozyme loci and nucleotide sequences of 3' end of the cytochrome b gene, tRNA genes, D-loop control region, and the 5' end of the 12S rRNA of mtDNA. Allozyme and mtDNA analyses revealing evident geographical structuring suggest limited gene flow between populations (F(ST)=0.511 and 0.791, respectively). Low genetic variability within populations (P=5.56%; He=0.018) based on allozymes and significantly negative Tajima's D statistics based on mtDNA suggest that most populations in Taiwan may have originated from a small number of founders followed by demographic expansion. The gene genealogy of mtDNA identified six lineages corresponding to major drainages that were separated by the geological barriers due to vicariant events. A minimum spanning network based on nucleotide substitutions reflects divergence from populations of the Miao-li Plateau to northern and southern regions of the island. In contrast to a previous hypothesis that suggests an early invasion to eastern part of Taiwan prior to the lifting of central mountain range some one million years ago, the mtDNA genealogy and molecular dating reveal very recent colonization of the eastern population. Nested clade analyses revealing significant associations between genetic structure and geographical division identify past fragmentation and range expansion as major phylogeographical events that shaped the geographical distribution of this species in Taiwan.     

The use of specialisation indices to predict vulnerability of coral‐feeding butterflyfishes to environmental change

In the absence of detailed assessments of extinction risk, ecological specialisation is often used as a proxy of vulnerability to environmental disturbances and extinction risk. Numerous indices can be used to estimate specialisation; however, the utility of these different indices to predict vulnerability to future environmental change is unknown. Here we compare the performance of specialisation indices using coral‐feeding butterflyfishes as a model group. Our aims were to 1) quantify the dietary preferences of three butterflyfish species across habitats with differing levels of resource availability; 2) investigate how estimates of dietary specialisation vary with the use of different specialisation indices; 3) determine which specialisation indices best inform predictions of vulnerability to environmental change; and 4) assess the utility of resource selection functions to inform predictions of vulnerability to environmental change. The relative level of dietary specialisation estimated for all three species varied when different specialisation indices were used, indicating that the choice of index can have a considerable impact upon estimates of specialisation. Specialisation indices that do not consider resource abundance may fail to distinguish species that primarily use common resources from species that actively target resources disproportionately more than they are available. Resource selection functions provided the greatest insights into the potential response of species to changes in resource availability. Examination of resource selection functions, in addition to specialisation indices, indicated that Chaetodon trifascialis was the most specialised feeder, with highly conserved dietary preferences across all sites, suggesting that this species is highly vulnerable to the impacts of climate‐induced coral loss on reefs. Our results indicate that vulnerability assessments based on some specialisation indices may be misleading and the best estimates of dietary specialisation will be provided by indices which incorporate resource availability measures, as well as assessing responses of species to changes in resource availability.