Inferring population history from microsatellite and enzyme data in serially introduced cane toads, Bufo marinus

Much progress has been made on inferring population history from molecular data. However, complex demographic scenarios have been considered rarely or have proved intractable. The serial introduction of the South-Central American cane Load Bufo marinas in various Caribbean and Pacific islands involves four major phases: a possible genetic admixture during the first introduction, a bottleneck associated with founding, a transitory, population boom, and finally, a demographic stabilization. A large amount of historical and demographic information is available for those introductions and can be combined profitably with molecular data. We used a Bayesian approach to combine this information With microsatellite (10 loci) and enzyme (22 loci) data and used a rejection algorithm to simultaneously estimate the demographic parameters describing the four major phases of the introduction history,. The general historical trends supported by microsatellites and enzymes were similar. However, there was a stronger support for a larger bottleneck at introductions for microsatellites than enzymes and for a more balanced genetic admixture for enzymes than for microsatellites. Verb, little information was obtained from either marker about the transitory population boom observed after each introduction. Possible explanations for differences in resolution of demographic events and discrepancies between results obtained with microsatellites and enzymes were explored. Limits Of Our model and method for the analysis of nonequilibrium populations were discussed.

Isolation and characterization of nine microsatellite loci from the Hawaiian grouper Epinephelus quernus (Serranidae) for population genetic analyses

The availability of variable genetic markers for groupers (Serranidae) has generally been limited to mitochondrial DNA. For studies of population genetic structure, more loci are usually required; particularly useful are those that are nuclear in origin such as microsatellites. Here, we isolated and characterized 9 microsatellite loci from the endemic Hawaiian grouper Epinephelus quernus using a biotin-labeled oligonucleotide-streptavidin-coated magnetic bead approach. Of the 20 repeat-containing fragments isolated, 15 had sufficient flanking region in which to design primers. Among these, 9 produced consistent polymerase chain reaction product, and 6 were highly variable. These 6 loci were all composed of dinucleotide repeats, with the number of alleles ranging from 6 to 18, and heterozygosities from 33.3% to 91.7%. The high levels of variability observed should make these markers useful for population genetic studies of E. quernus, and potentially other epinephelines.

Oceanic interchange and nonequilibrium population structure in the estuarine dependent Indo-Pacific tasselfish, Polynemus sheridani

We assayed mtDNA haplotype [300 base pairs (bp) control region] geography and genealogy in the Indo-Pacific tasselfish, Polynemus sheridani from its contiguous estuarine distribution across northern Australia (n = 169). Eight estuaries were sampled from three oceanographic regions (Timor Sea, Gulf of Carpentaria and the Coral Sea) to assess the impact of Pleistocene sea level changes on the historical connectivity among P. sheridani populations. Specifically, we investigated the genetic consequences of disruption to Indian-Pacific Ocean connectivity brought about by the closure of the Torres Strait. Overall there was significant population subdivision among estuaries (F-ST = 0.161, (Phi(ST) = 0.187). Despite a linear distribution, P. sheridani did not show isolation by distance over the entire sampled range because of genetic similarity of estuaries greater than 3000 km apart. However, significant isolation by distance was detected between estuaries separated by less than 3000 km of coastline. Unlike many genetic studies of Indo-Pacific marine species, there was no evidence for an historical division between eastern and western populations. Instead, phylogeographical patterns were dominated by a starlike intraspecific phylogeny coupled with evidence for population expansion in both the Gulf of Carpentaria and the Coral Sea but not the Timor Sea. This was interpreted as evidence for recent west to east recolonization across of northern Australia following the last postglacial marine advance. We argue that although sufficient time has elapsed postcolonization for populations to approach gene flow-drift equilibrium over smaller spatial scales (< 3000 km), the signal of historical colonization persists to obscure the expected equilibrium pattern of isolation by distance over large spatial scales (> 3000 km).

Genetic diversity and population structure of Tasmanian devils, the largest marsupial carnivore

Genetic diversity and population structure were investigated across the core range of Tasmanian devils (Sarcophilus laniarius; Dasyuridae), a wide-ranging marsupial carnivore restricted to the island of Tasmania. Heterozygosity (0.386-0.467) and allelic diversity (2.7-3.3) were low in all subpopulations and allelic size ranges were small and almost continuous, consistent with a founder effect. Island effects and repeated periods of low population density may also have contributed to the low variation. Within continuous habitat, gene flow appears extensive up to 50 km (high assignment rates to source or close neighbour populations; nonsignificant values of pairwise F-ST), in agreement with movement data. At larger scales (150-250 km), gene flow is reduced (significant pairwise F-ST) but there is no evidence for isolation by distance. The most substantial genetic structuring was observed for comparisons spanning unsuitable habitat, implying limited dispersal of devils between the well-connected, eastern populations and a smaller northwestern population. The genetic distinctiveness of the northwestern population was reflected in all analyses: unique alleles; multivariate analyses of gene frequency (multidimensional scaling, minimum spanning tree, nearest neighbour); high self-assignment (95%); two distinct populations for Tasmania were detected in isolation by distance and in Bayesian model-based clustering analyses. Marsupial carnivores appear to have stronger population subdivisions than their placental counterparts.

An investigation into effects of long-distance seed dispersal on organelle population genetic structure and colonization rate: a model analysis

A simulation-based modelling approach is used to examine the effects of stratified seed dispersal (representing the distribution of the majority of dispersal around the maternal parent and also rare long-distance dispersal) on the genetic structure of maternally inherited genomes and the colonization rate of expanding plant populations. The model is parameterized to approximate postglacial oak colonization in the UK, but is relevant to plant populations that exhibit stratified seed dispersal. The modelling approach considers the colonization of individual plants over a large area (three 500 km x 10 km rolled transects are used to approximate a 500 km x 300 km area). Our approach shows how the interaction of plant population dynamics with stratified dispersal can result in a spatially patchy haplotype structure. We show that while both colonization speeds and the resulting genetic structure are influenced by the characteristics of the dispersal kernel, they are robust to changes in the periodicity of long-distance events, provided the average number of long-distance dispersal events remains constant. We also consider the effects of additional physical and environmental mechanisms on plant colonization. Results show significant changes in genetic structure when the initial colonization of different haplotypes is staggered over time and when a barrier to colonization is introduced. Environmental influences on survivorship and fecundity affect both the genetic structure and the speed of colonization. The importance of these mechanisms in relation to the postglacial spread and genetic structure of oak in the UK is discussed.

Assessing moth migration and population structuring in Helicoverpa annigera (Lepidoptera : Noctuidae) at the regional scale: Example from the Darling Downs, Australia

Analysis of gene flow and migration of Helicoverpa armigera (Hubner) in a major cropping region of Australia identified substantial genetic structuring, migration events, and significant population genotype changes over the 38-mo sample period from November 1999 to January 2003. Five highly variable microsatellite markers were used to analyze 916 individuals from 77 collections across 10 localities in the Darling Downs. The molecular data indicate that in some years (e.g., April 2002-March 2003), low levels of H. armigera migration and high differentiation between populations occurred, whereas in other years (e.g., April 2001-March 2002), there were higher levels of adult moth movement resulting in little local structuring of populations. Analysis of populations in other Australian cropping regions provided insight into the quantity and direction of immigration of H. armigera adults into the Darling Downs growing region of Australia. These data provide evidence adult moth movement differs from season to season, highlighting the importance of studies in groups such as the Lepidoptera extending over consecutive years, because short-term sampling may be misleading when population dynamics and migration change so significantly. This research demonstrates the importance of maintaining a coordinated insecticide resistance management strategy, because in some years H. armigera populations may be independent within a region and thus significantly influenced by local management practices; however, periods with high migration will occur and resistance may rapidly spread.

Regional and population-scale influences on genetic diversity partitioning within costa rican populations of the pioneer tree Vochysia ferruginea mart.

The neotropical pioneer species Vochysia ferruginea is locally important for timber and is being increasingly exploited. The sustainable utilisation of this species would benefit from an understanding of the level and partitioning of genetic diversity within remnant and secondary regrowth populations. We used data from total genome (amplified fragment length polymorphism, AFLP) and chloroplast genome markers to assay diversity levels within seven Costa Rican populations. Significant chloroplast differentiation between Atlantic and Pacific watersheds was observed, suggesting divergent historical origins for these populations. Contemporary gene flow, though extensive, is geographically constrained and a clear pattern of isolation by distance was detectable when an inter-population distance representing gene flow around the central Costa Rican mountain range was used. Overall population differentiation was low (F-ST = 0.15) and within-population diversity high, though variable (H-s=0.16-0.32), which fits with the overall pattern of population genetic structure expected for a widespread, outcrossed tropical tree. However genetic diversity was significantly lower and differentiation higher for recently colonised and disturbed populations compared to that at more established sites. Such a pattern seems indicative of a pioneer species undergoing repeated cycles of colonisation and succession.

Population dynamics and gene flow of Helicoverpa armigera (Lepidoptera : Noctuidae) on cotton and grain crops in the Murrumbidgee Valley, Australia

The population dynamics of Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae) in the Murrumbidgee Valley, Australia, has been characterized using five highly variable microsatellite loci. In the 2001-2002 growing season, there were very high levels of migration into the Murrumbidgee Valley with no detectable genetic structuring, consistent with previous analyses on a national scale. By contrast, there was significant genetic structuring over the 2002-2003 growing season, with three distinct genetic types detected. The first type corresponded to the first two generations and was derived from local individuals emerging from diapause and their progeny. The second genetic type corresponded to generation 3 and resulted from substantial immigration into the region. There was another genetic shift in generation 4, which accounts for the third genetic type of the season. This genetic shift occurred despite low levels of immigration. During the third generation of the 2002-2003 growing season, different population dynamics was characterized for H. armigera on maize, Zea mays L., and cotton Gossipium hirsutum L. Populations on cotton tended to cycle independently with very little immigration from outside the region or from maize within the region. Maize acted as a major sink for immigrants from cotton and from outside the region. If resistance were to develop on cotton under these circumstances, susceptible individuals from maize or from other regions would not dilute this resistance. In addition, resistance is likely to be transferred to maize and be perpetuated until diapause, from where it may reemerge next season. If low levels of immigration were to occur on transgenic cotton, this may undermine the effectiveness of refugia, especially noncotton refugia.

Significant patterns of population genetic structure and limited gene flow in a threatened macropodid marsupial despite continuous habitat in southeast Queensland, Australia

Many endangered species worldwide are found in remnant populations, often within fragmented landscapes. However, when possible, an understanding of the natural extent of population structure and dispersal behaviour of threatened species would assist in their conservation and management. The brush-tailed rock-wallaby (Petrogale penicillata), a once abundant and widespread rock-wallaby species across southeastern Australia, has become nearly extinct across much of the southern part of its range. However, the northern part of the species' range still sustains many small colonies closely distributed across suitable habitat, providing a rare opportunity to investigate the natural population dynamics of a listed threatened species. We used 12 microsatellite markers to investigate genetic diversity, population structure and gene flow among brush-tailed rock-wallaby colonies within and among two valley regions with continuous habitat in southeast Queensland. We documented high and signifcant levels of population genetic structure between rock-wallaby colonies embedded in continuous escarpment habitat and forest. We found a strong and significant pattern of isolation-by-distance among colonies indicating restricted gene flow over a small geographic scale (< 10 km) and conclude that gene flow is more likely limited by intrinsic factors rather than environmental factors. In addition, we provide evidence that genetic diversity was significantly lower in colonies located in a more isolated valley region compared to colonies located in a valley region surrounded by continuous habitat. These findings shed light on the processes that have resulted in the endangered status of rock-wallaby species in Australia and they have strong implications for the conservation and management of both the remaining 'connected' brush-tailed rock-wallaby colonies in the northern parts of the species' range and the remnant endangered populations in the south.

Carapace dentition patterns, morphometrics and allozyme differentiation amongst two toothed freshwater crab species (Potamonautes warreni and P-unispinus) (Decapoda : Brachyura : Potamonautidae) from river systems in South Africa

The taxonomic relationship between two toothed South African river crabs, Potamonautes warreni and P. unispinus, is unclear. The problem stems from the widespread variation in carapace dentition patterns amongst P. warreni individuals over its biogeographic range, where single toothed individuals may appear similar in carapace morphology to P. unispinus. Ten populations of P. warreni and 18 populations of P. unispinus were collected and the morphometric and genetic differentiation between the two taxa quantified. Patterns of morphometric and genetic variation were examined using multivariate statistics and protein gel electrophoresis, respectively. Principal component analyses of carapace characters showed that the two species are morphologically indistinguishable. However, discriminate functions analyses and additional statistical results corroborate the morphological distinction between the two taxa. Allozyme electrophoresis of 17 protein coding loci, indicated a close genetic similarity between the two species (I = 0.92). A fixed allelic difference at one locus (LT-2) and extensive genetic variability at another locus (PGM-1) indicate that two gene pools are present and that the two taxa are genetically isolated. Intraspecific genetic I values for both species were > 0.97 and indicated no apparent genetic structuring on a micro or macro-geographic scale. The variation in carapace dentition among P. warreni populations possesses no genetic basis and may possibly toe the product of ecogenesis. The value of dentition patterns in the systematics of river crabs is discussed. Dentition patterns among river crab species appear to be conserved and reliable as species specific diagnostic markers, but should ideally be used in combination with other morphological data sets and genetic evidence.

Reconciling paleodistribution models and comparative phylogeography in the Wet Tropics rainforest land snail Gnarosophia bellendenkerensis (Brazier 1875)

Comparative phylogeography has proved useful for investigating biological responses to past climate change and is strongest when combined with extrinsic hypotheses derived from the fossil record or geology. However, the rarity of species with sufficient, spatially explicit fossil evidence restricts the application of this method. Here, we develop an alternative approach in which spatial models of predicted species distributions under serial paleoclimates are compared with a molecular phylogeography, in this case for a snail endemic to the rainforests of North Queensland, Australia. We also compare the phylogeography of the snail to those from several endemic vertebrates and use consilience across all of these approaches to enhance biogeographical inference for this rainforest fauna. The snail mtDNA phylogeography is consistent with predictions from paleoclimate modeling in relation to the location and size of climatic refugia through the late Pleistocene-Holocene and broad patterns of extinction and recolonization. There is general agreement between quantitative estimates of population expansion from sequence data (using likelihood and coalescent methods) vs. distributional modeling. The snail phylogeography represents a composite of both common and idiosyncratic patterns seen among vertebrates, reflecting the geographically finer scale of persistence and subdivision in the snail. In general, this multifaceted approach, combining spatially explicit paleoclimatological models and comparative phylogeography, provides a powerful approach to locating historical refugia and understanding species' responses to them.