The biology of pygmy whitefish, Prosopium coulterii, in a closed sub-boreal lake: spatial distribution and diel movements

Using gillnets and trap nets, we examined the spatial distribution, diel movements, and environmental tolerances of pygmy whitefish, Prosopium coulterii, in a small boreal lake in north-central British Columbia. Most gillnets were set below the thermocline but we also fished a shore net in the littoral zone. During the ice-free season (May to November) there was a strong diel onshore–offshore movement: during the day pygmy whitefish were offshore and below the thermocline (water temperatures of 4–6°C) but at night they were inshore and above the thermocline (water temperatures of 12–18°C). This onshore–offshore movement occurred close to the bottom and, regardless of where they were caught, most fish were <4 m off the bottom. Oxygen concentrations in most of the hypolimnion dropped to <5.0 mg l⁻¹ in June and by late August to <1.0 mg l⁻¹; indicating pygmy whitefish can tolerate low oxygen conditions. The catch of pygmy whitefish in gillnets set below the thermocline was highly skewed: 53% of the nets were empty, 37% caught 18 or less fish, and 10% caught 70% of the total catch (742 fish). Trap nets produced similarly skewed results: most trap net sets caught no pygmy whitefish but one set caught over 2,000 individuals. Our catch data suggest that in Dina Lake #1 some pygmy whitefish aggregate.     

Phylogeography and the origins of range disjunctions in a north temperate fish,the pygmy whitefish (Prosopium coulterii), inferred from mitochondrial and nuclear DNA sequence analysis

Aim To investigate the degree of phylogeographical divergence within pygmy whitefish (Prosopium coulterii) and to test hypotheses concerning the origin of disjunct populations within North America. Location North America from western Alaska to Lake Superior. Methods Mitochondrial (ATPase subunit VI) and nuclear (ITS‐1, ITS‐2) DNA sequence variation was assessed across the species’ North American range to test for the existence of distinct phylogeographical groupings of pygmy whitefish associated with known glacial refugia. Coalescent simulations of the mitochondrial DNA (mtDNA) data were used to test hypotheses of population structure. Results This species is composed of two monophyletic mitochondrial clades across its North American range. The two mtDNA clades differed by an average 3.3% nucleotide sequence divergence. These clades were also distinguished by ITS‐2, but the relationships among lineages were not resolved by the ITS‐1 analysis. Coalescent analyses rejected the null hypothesis of the current disjunct distributions being a result of fragmentation of a single widespread ancestral lineage across a variety of effective population sizes and divergence times. Main conclusions The current range disjunctions of pygmy whitefish in North America probably resulted from isolation, genetic divergence, and selective dispersal from at least two major Pleistocene glacial refugia: Beringia and Cascadia. More recent isolation and dispersal from an upper Mississippi refugium is suggested by relationships within one of the clades and by distributional evidence from co‐distributed species. The Beringian and Cascadian refugia have played major roles in the zoogeography of Nearctic temperate aquatics, but the roles of smaller refugia appear more variable among other species.     

Real‐time PCR identification of lake whitefish Coregonus clupeaformis in the Laurentian Great Lakes

The purpose of this study was to develop a real‐time PCR assay to specifically identify lake whitefish Coregonus clupeaformis in larval fish assemblages based on a 122 bp amplicon from the mitochondrial genome. The efficiency of the reaction, as calculated from the standard curve, was 90·77% with the standard curve having an r² value of 0·998. Specificity of the assay provided single melt peak in a melt‐curve analysis and amplification of only the target species. The assay successfully identified target DNA in as low as 0·1% proportion of a DNA mixture. This assay was designed on the portable Smart Cycler II platform and can be used in both field and laboratory settings to successfully identify C. clupeaformis.     

Unique sympatric quartet of limnetic,benthic, profundal and piscivorous brown trout populations resolved by 3D sampling and focused molecular marker selection

相似文献   

Divergence,gene flow,and the origin of leapfrog geographic distributions: The history of colour pattern variation in Phyllobates poison‐dart frogs

The geographic distribution of phenotypic variation among closely related populations is a valuable source of information about the evolutionary processes that generate and maintain biodiversity. Leapfrog distributions, in which phenotypically similar populations are disjunctly distributed and separated by one or more phenotypically distinct populations, represent geographic replicates for the existence of a phenotype, and are therefore especially informative. These geographic patterns have mostly been studied from phylogenetic perspectives to understand how common ancestry and divergent evolution drive their formation. Other processes, such as gene flow between populations, have not received as much attention. Here, we investigate the roles of divergence and gene flow between populations in the origin and maintenance of a leapfrog distribution in Phyllobates poison frogs. We found evidence for high levels of gene flow between neighbouring populations but not over long distances, indicating that gene flow between populations exhibiting the central phenotype may have a homogenizing effect that maintains their similarity, and that introgression between ‘leapfroging’ taxa has not played a prominent role as a driver of phenotypic diversity in Phyllobates. Although phylogenetic analyses suggest that the leapfrog distribution was formed through independent evolution of the peripheral (i.e. leapfrogging) populations, the elevated levels of gene flow between geographically close populations poise alternative scenarios, such as the history of phenotypic change becoming decoupled from genome‐averaged patterns of divergence, which we cannot rule out. These results highlight the importance of incorporating gene flow between populations into the study of geographic variation in phenotypes, both as a driver of phenotypic diversity and as a confounding factor of phylogeographic inferences.     

Phylogeography of the brown hare (Lepus europaeus) in Europe: a legacy of south‐eastern Mediterranean refugia?

Aim We analysed the population genetics of the brown hare (Lepus europaeus) in order to test the hypothesis that this species migrated into central Europe from a number of late glacial refugia, including some in Asia Minor. Location Thirty‐three localities in Greece, Bulgaria, Italy, Croatia, Serbia, Poland, Switzerland, Austria, France, Germany, the Netherlands, Spain, the United Kingdom, Turkey and Israel. Methods In total, 926 brown hares were analysed for mitochondrial DNA (mtDNA) variation by restriction fragment length polymorphism (RFLP) performed on polymerase chain reaction‐amplified products spanning cytochrome b (cyt b)/control region (CR), cytochrome oxidase I (COI) and 12S–16S rRNA. In addition, sequence analysis of the mtDNA CR‐I region was performed on 69 individuals, and the data were compared with 137 mtDNA CR‐I sequences retrieved from GenBank. Results The 112 haplotypes detected were partitioned into five phylogeographically well‐defined major haplogroups, namely the ‘south‐eastern European type haplogroup’ (SEEh), ‘Anatolian/Middle Eastern type haplogroup’ (AMh), ‘European type haplogroup, subgroup A’ (EUh‐A), ‘European type haplogroup, subgroup B’ (EUh‐B) and ‘Intermediate haplogroup’ (INTERh). Sequence data retrieved from GenBank were consistent with the haplogroups determined in this study. In Bulgaria and north‐eastern Greece numerous haplotypes of all five haplogroups were present, forming a large overlap zone. Main conclusions The mtDNA results allow us to infer post‐glacial colonization of large parts of Europe from a late glacial/early Holocene source population in the central or south‐central Balkans. The presence of Anatolian/Middle Eastern haplotypes in the large overlap zone in Bulgaria and north‐eastern Greece reveals gene flow from Anatolia to Europe across the late Pleistocene Bosporus land‐bridge. Although various restocking operations could be partly responsible for the presence of unexpected haplotypes in certain areas, we nevertheless trace a strong phylogeographic signal throughout all regions under study. Throughout Europe, mtDNA results indicate that brown hares are not separated into discernable phyletic groups.     

Pan‐African phylogeography of a model organism,the African clawed frog ‘Xenopus laevis’

The African clawed frog Xenopus laevis has a large native distribution over much of sub‐Saharan Africa and is a model organism for research, a proposed disease vector, and an invasive species. Despite its prominent role in research and abundance in nature, surprisingly little is known about the phylogeography and evolutionary history of this group. Here, we report an analysis of molecular variation of this clade based on 17 loci (one mitochondrial, 16 nuclear) in up to 159 individuals sampled throughout its native distribution. Phylogenetic relationships among mitochondrial DNA haplotypes were incongruent with those among alleles of the putatively female‐specific sex‐determining gene DM‐W, in contrast to the expectation of strict matrilineal inheritance of both loci. Population structure and evolutionarily diverged lineages were evidenced by analyses of molecular variation in these data. These results further contextualize the chronology, and evolutionary relationships within this group, support the recognition of X. laevis sensu stricto, X. petersii, X. victorianus and herein revalidated X. poweri as separate species. We also propose that portions of the currently recognized distributions of X. laevis (north of the Congo Basin) and X. petersii (south of the Congo Basin) be reassigned to X. poweri.     

Substantial contribution of invertebrates to the diet of a winter seed‐eater,the reed bunting Emberiza schoeniclus,wintering in a sewage farm in south‐western Poland

We investigated whether the winter diet of a typical seed‐eating bird is hard‐wired in the context of evolutionary hypothesis for granivory. We examined the diet composition of ‘a small‐billed form’ of the reed bunting Emberiza schoeniclus wintering in a sewage farm in south‐western Poland (Central Europe), where unfrozen wastewater provides various groups of invertebrate prey. The analysis of droppings (N = 151) collected from four different feeding grounds located in reedbeds and grasslands inundated with waste‐water showed the substantial contribution of invertebrates in the diet of reed buntings. Across four sample areas, the frequency of invertebrates in faecal samples ranged between 37% to 80%. In total, we identified 194 animal prey, mainly spiders Araneae (53% of all identified invertebrate prey), and several taxa of Coleoptera (43%). Among plant food (N = 8357 identified items), the most numerous were shells of weed seeds, namely Amaranthus sp. (56.8%), Urtica dioica (22.6%), and Chenopodium sp. (19.6%). Our results showed that, during winter, the reed bunting is not an obligatory seed‐eater. This species may exploit both animal and plant food; hence, our results indicate that the reed buntings wintering in temperate Europe may feed more opportunistically than was previously assumed. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 108 , 429–433.     

In defence of model‐based inference in phylogeography

Recent papers have promoted the view that model‐based methods in general, and those based on Approximate Bayesian Computation (ABC) in particular, are flawed in a number of ways, and are therefore inappropriate for the analysis of phylogeographic data. These papers further argue that Nested Clade Phylogeographic Analysis (NCPA) offers the best approach in statistical phylogeography. In order to remove the confusion and misconceptions introduced by these papers, we justify and explain the reasoning behind model‐based inference. We argue that ABC is a statistically valid approach, alongside other computational statistical techniques that have been successfully used to infer parameters and compare models in population genetics. We also examine the NCPA method and highlight numerous deficiencies, either when used with single or multiple loci. We further show that the ages of clades are carelessly used to infer ages of demographic events, that these ages are estimated under a simple model of panmixia and population stationarity but are then used under different and unspecified models to test hypotheses, a usage the invalidates these testing procedures. We conclude by encouraging researchers to study and use model‐based inference in population genetics.     

Life‐history predicts past and present population connectivity in two sympatric sea stars

Life‐history traits, especially the mode and duration of larval development, are expected to strongly influence the population connectivity and phylogeography of marine species. Comparative analysis of sympatric, closely related species with differing life histories provides the opportunity to specifically investigate these mechanisms of evolution but have been equivocal in this regard. Here, we sample two sympatric sea stars across the same geographic range in temperate waters of Australia. Using a combination of mitochondrial DNA sequences, nuclear DNA sequences, and microsatellite genotypes, we show that the benthic‐developing sea star, Parvulastra exigua, has lower levels of within‐ and among‐population genetic diversity, more inferred genetic clusters, and higher levels of hierarchical and pairwise population structure than Meridiastra calcar, a species with planktonic development. While both species have populations that have diverged since the middle of the second glacial period of the Pleistocene, most P. exigua populations have origins after the last glacial maxima (LGM), whereas most M. calcar populations diverged long before the LGM. Our results indicate that phylogenetic patterns of these two species are consistent with predicted dispersal abilities; the benthic‐developing P. exigua shows a pattern of extirpation during the LGM with subsequent recolonization, whereas the planktonic‐developing M. calcar shows a pattern of persistence and isolation during the LGM with subsequent post‐Pleistocene introgression.     

Population structure of Squatina guggenheim (Squatiniformes,Squatinidae) from the south‐western Atlantic Ocean

Population genetic analyses based on both mitochondrial cytochrome b and the internal transcribed spacer 2 of recombinant (r)DNA genes were implemented to examine hypotheses of population differentiation in the angular angel shark Squatina guggenheim, one of the four most‐widespread endemic species inhabiting coastal ecosystems in the south‐western Atlantic Ocean. A total of 82 individuals of S. guggenheim from 10 sampling sites throughout the Río de la Plata mouth, its maritime front, the outer shelf at the subtropical confluence and the coastal areas of the south‐west Atlantic Ocean, were included. The analysis of molecular variance (AMOVA) based on the second internal transcribed spacer (its‐2) region supports that the samples from the outer shelf represent an isolated group from other sites. Historical gene flow in a coalescent‐based approach revealed significant immigration and emigration asymmetry between sampling sites. Based on the low level of genetic diversity, the existence of a long‐term population decline or a past recent population expansion following a population bottleneck could be proposed in S. guggenheim. This demographic differentiation suggests a degree of vulnerability to overexploitation in this endemic and endangered south‐west Atlantic Ocean shark, given its longevity and low reproductive potential.     

A comparison of DNA extraction methods for high‐throughput DNA analyses

The inclusion of next‐generation sequencing technologies in population genetic and phylogenetic studies has elevated the need to balance time and cost of DNA extraction without compromising DNA quality. We tested eight extraction methods – ranging from low‐ to high‐throughput techniques – and eight phyla: Annelida, Arthropoda, Cnidaria, Chordata, Echinodermata, Mollusca, Ochrophyta and Porifera. We assessed DNA yield, purity, efficacy and cost of each method. Extraction efficacy was quantified using the proportion of successful polymerase chain reaction (PCR) amplification of two molecular markers for metazoans (mitochondrial COI and nuclear histone 3) and one for Ochrophyta (mitochondrial nad6) at four time points – 0.5, 1, 2 and 3 years following extraction. DNA yield and purity were quantified using NanoDrop absorbance ratios. Cost was estimated in terms of time and material expense. Results show differences in DNA yield, purity and PCR success between extraction methods and that performance also varied by taxon. The traditional time‐intensive, low‐throughput CTAB phenol–chloroform extraction performed well across taxa, but other methods also performed well and provide the opportunity to reduce time spent at the bench and increase throughput.     

Trophic ecology of yellownose skate Zearaja chilensis,a top predator in the south‐western Atlantic Ocean

The diet and trophic level (T_L) of the yellownose skate Zearaja chilensis in the south‐western Atlantic Ocean (35°–54° S), and how these varied in relation to body size, sex, maturity stage, depth and region were determined by analysis of stomach contents. From 776 specimens analysed, 671 (86·5%) ranging from 180 to 1190 mm total length (L_T) had prey in their stomachs. The diet was dominated by fishes, mainly the notothenioid Patagonotothen ramsayi and the Argentine hake Merluccius hubbsi. The consumption of fishes and crabs increased with increasing predator size, and these preys were more important in the north than in the south. Isopods and other crustaceans were consumed more in the south and their consumption decreased as the size of Z. chilensis increased. The T_L of Z. chilensis increased with L_T from 4·29 to 4·59 (mean 4·53), confirming their ecological role as a top predator. The small and large size classes exhibited a low diet overlap and the highest spatial segregation, whereas medium and large specimens had higher co‐occurrence and dietary overlap indices. A clear distinction in tooth shape was noted between sexes in adult specimens, with males having longer cusps. This sexual heterodonty may be related to reproductive behaviour, increasing the grasping ability of males during courtship, because there were no differences in diet between the sexes.     

Morphofunctional analysis of Svobodaina species (Brachiopoda,Heterorthidae) from south‐western Europe

Geometric morphometric methods applied to the ventral muscle field outline of the various species of Svobodaina from south‐western Europe identify criteria for discriminating among this key cluster of brachiopod species. These data indicate the close relationship between the patterns of the ventral muscle field in the Svobodaina species and the environmental conditions where each inhabited; the fields are better developed in species related to high‐energy environments. For example, S. armoricana, with the smallest diductor scars, would have inhabited the quiet marine environments of the lower offshore or within protected lagoonal settings. S. feisti would have inhabited the upper offshore, a more energetic environment than S. armoricana. Finally, S. havliceki, with the largest diductor scars, would have thrived in the most energetic environments among all the south‐western European Svobodaina species, living just above the fair‐weather wave base in the lower shoreface. The palaeoecological results suggest a distribution of Svobodaina species during the Late Ordovician along an onshore–offshore transect across the shallow marine platforms of the Mediterranean margin of Gondwana. On the other hand, the occurrence in some localities of several species with overlapping ranges or within the same assemblage indicates that the biostratigraphical efficacy of the genus is restricted. Thus, the previously defined taxon‐range biozones characterized by Svobodaina species of the north Gondwanan margin are in need of reassessment. The morphology of Svobodaina may be a considerable aid to environmental analyses rather than to precise biostratigraphical correlations.     

The advantages of going large: genome‐wide SNPs clarify the complex population history and systematics of the threatened western pond turtle

As the field of phylogeography has matured, it has become clear that analyses of one or a few genes may reveal more about the history of those genes than the populations and species that are the targets of study. To alleviate these concerns, the discipline has moved towards larger analyses of more individuals and more genes, although little attention has been paid to the qualitative or quantitative gains that such increases in scale and scope may yield. Here, we increase the number of individuals and markers by an order of magnitude over previously published work to comprehensively assess the phylogeographical history of a well‐studied declining species, the western pond turtle (Emys marmorata). We present a new analysis of 89 independent nuclear SNP markers and one mitochondrial gene sequence scored for rangewide sampling of >900 individuals, and compare these to smaller‐scale, rangewide genetic and morphological analyses. Our enlarged SNP data fundamentally revise our understanding of evolutionary history for this lineage. Our results indicate that the gains from greatly increasing both the number of markers and individuals are substantial and worth the effort, particularly for species of high conservation concern such as the pond turtle, where accurate assessments of population history are a prerequisite for effective management.     

Comparative phylogeography of two co‐distributed but ecologically distinct rainbowfishes of far‐northern Australia

Aim

To test the influence of historical and contemporary environment in shaping the genetic diversity of freshwater fauna we contrast genetic structure in two co‐distributed, but ecologically distinct, rainbowfish; a habitat generalist (Melanotaenia splendida) and a habitat specialist (M. trifasciata).

Location

Fishes were sampled from far northern Australia (Queensland and Northern Territory).

Methods

We used sequence data from one mitochondrial gene and one nuclear gene to investigate patterns of genetic diversity in M. splendida and M. trifasciata to determine how differences in habitat preference and historical changes in drainage boundaries affected patterns of connectivity.

Results

Melanotaenia splendida showed high levels of genetic diversity and little population structure across its range. In contrast, M. trifasciata showed high levels of population structure. Whereas phylogeographic patterns differed, both species showed a strong relationship between geographical distance and genetic differentiation between populations. Melanotaenia splendida showed a shallower relationship with geographical distance, and genetic differentiation was best explained by stream length and a lower scaled ocean distance (11.98 times coast length). For M. trifasciata, genetic differentiation was best explained by overwater distance between catchments and ocean distance scaled at 1.16 × 10⁶ times coast length.

Main conclusions

Connectivity of freshwater populations inhabiting regions periodically interconnected during glacial periods appears to have been affected by ecological differences between species. Species‐specific differences are epitomized here by the contrast between co‐distributed congeners with different habitat requirements: for the habitat generalist, M. splendida, there was evidence for greater historical genetic connectivity with oceans as a weaker barrier to gene exchange in contrast with the habitat specialist, M. trifasciata.     

Investigating the extent of parallelism in morphological and genomic divergence among lake trout ecotypes in Lake Superior

Understanding the emergence of species through the process of ecological speciation is a central question in evolutionary biology which also has implications for conservation and management. Lake trout (Salvelinus namaycush) is renowned for the occurrence of different ecotypes linked to resource and habitat use throughout North America. We aimed to unravel the fine genetic structure of the four lake trout ecotypes in Lake Superior. A total of 486 individuals from four sites were genotyped at 6822 filtered SNPs using RADseq technology. Our results revealed different extent of morphological and genetic differentiation within the different sites. Overall, genetic differentiation was weak but significant and was on average three times higher between sites (mean F_ST = 0.016) than between ecotypes within sites (mean F_ST = 0.005) indicating higher level of gene flow or a more recent shared ancestor between ecotypes within each site than between populations of the same ecotype. Evidence of divergent selection was also found between ecotypes and/or in association with morphological variation. Outlier loci found in genes related to lipid metabolism and visual acuity were of particular interest in this context of ecotypic divergence. However, we did not find clear indication of parallelism at the genomic level, despite the presence of phenotypic parallelism among some ecotypes from different sampling sites. Overall, the occurrence of different levels of both genomic and phenotypic differentiation between ecotypes within each site with several differentiated loci linked to relevant biological functions supports the presence of a continuum of divergence in lake trout.     

Weak large‐scale population genetic structure in a philopatric seabird,the European Shag Phalacrocorax aristotelis

Quantifying population genetic structure is fundamental to testing hypotheses regarding gene flow, population divergence and dynamics across large spatial scales. In species with highly mobile life‐history stages, where it is unclear whether such movements translate into effective dispersal among discrete philopatric breeding populations, this approach can be particularly effective. We used seven nuclear microsatellite loci and mitochondrial DNA (ND2) markers to quantify population genetic structure and variation across 20 populations (447 individuals) of one such species, the European Shag, spanning a large geographical range. Despite high breeding philopatry, rare cross‐sea movements and recognized subspecies, population genetic structure was weak across both microsatellites and mitochondrial markers. Furthermore, although isolation‐by‐distance was detected, microsatellite variation provided no evidence that open sea formed a complete barrier to effective dispersal. These data suggest that occasional long‐distance, cross‐sea movements translate into gene flow across a large spatial scale. Historical factors may also have shaped contemporary genetic structure: cluster analyses of microsatellite data identified three groups, comprising colonies at southern, mid‐ and northern latitudes, and similar structure was observed at mitochondrial loci. Only one private mitochondrial haplotype was found among subspecies, suggesting that this current taxonomic subdivision may not be mirrored by genetic isolation.     

Mitochondrial DNA and nuclear microsatellites reveal high diversity and genetic structure in an avian top predator,the white‐tailed sea eagle,in central Europe

We analysed 123 white‐tailed sea eagles (Haliaeetus albicilla) from (primarily central) Europe with respect to variability and differentiation based on 499 bp of the mitochondrial control region and genotypes at seven unlinked nuclear microsatellites. Variability was high (overall expected heterozygosity, haplotype and nucleotide diversity being 0.70, 0.764 and 0.00698, respectively) and both marker systems showed a subdivision into two main genetic clusters (microsatellites) or haplogroups (mtDNA). In line with earlier analyses focusing on populations from northern and eastern Europe, as well as from Asia, we found a high level of admixture in Europe and no signs of a bottleneck – despite a severe decline of white‐tailed sea eagle populations during the 20^th century. Europe is thus a global stronghold for this species not only with respect to the number of breeding pairs but also regarding the proportion of species‐wide genetic diversity. Our dense sampling revealed a possibly clinal variation within central Europe from north‐west to south‐east that was reflected by the distribution of mtDNA haplotypes as well as the two microsatellite‐based clusters. This population differentiation in central Europe probably originated from a geographically structured postglacial colonization and was later enhanced by recent demographic fluctuations. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 727–737.     

Testing models of speciation from genome sequences: divergence and asymmetric admixture in Island South‐East Asian Sus species during the Plio‐Pleistocene climatic fluctuations

In many temperate regions, ice ages promoted range contractions into refugia resulting in divergence (and potentially speciation), while warmer periods led to range expansions and hybridization. However, the impact these climatic oscillations had in many parts of the tropics remains elusive. Here, we investigate this issue using genome sequences of three pig (Sus) species, two of which are found on islands of the Sunda‐shelf shallow seas in Island South‐East Asia (ISEA). A previous study revealed signatures of interspecific admixture between these Sus species (Genome biology, 14 , 2013, R107). However, the timing, directionality and extent of this admixture remain unknown. Here, we use a likelihood‐based model comparison to more finely resolve this admixture history and test whether it was mediated by humans or occurred naturally. Our analyses suggest that interspecific admixture between Sunda‐shelf species was most likely asymmetric and occurred long before the arrival of humans in the region. More precisely, we show that these species diverged during the late Pliocene but around 23% of their genomes have been affected by admixture during the later Pleistocene climatic transition. In addition, we show that our method provides a significant improvement over D‐statistics which are uninformative about the direction of admixture.