Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

High nuclear genetic diversity, high levels of outcrossing and low differentiation among remnant populations of Quercus petraea at the margin of its range in Ireland

Background and Aims Quercus petraea colonized Ireland after the last glaciation from refugia on mainland Europe. Deforestation. however. beginning in Neolithic times, has resulted in small, scattered forest fragments, now covering less than 12 000 ha. Methods Plastid (three fragments) and microsatellite variation (13 loci) were characterized in seven Irish populations sampled along a north-south gradient. Using Bayesian approaches and Wright's F-statistics, the effects of colonization and fragmentation on the genetic structure and mating patterns of extant oak populations were investigated. Key-Results All Populations possessed cytotypes common to the Iberian Peninsula. Despite the distance from the refugial core and the extensive deforestation in Ireland, nuclear genetic variation was high and comparable to mainland Europe. Low population differentiation was observed within Ireland and populations showed no evidence for isolation by distance. As expected of a marker with an effective Population size of one-quarter relative to the nuclear genome, plastid variation indicated higher differentiation. Individual inbreeding coefficients indicated high levels of outcrossing. Conclusions Consistent with a large effective Population size in the historical migrant gene pool and/or with high gene flow among populations, high within-population diversity and low population differentiation was observred within Ireland. It is proposed that native Q. petraea populations in Ireland share a common phylogeographic history and that the present genetic structure does not reflect founder effects. (C) 2004 Annals of Botany Company.

The genetic structure of Australian populations of Mycosphaerella musicola suggests restricted gene flow at the continental scale

Mycosphaerello musicolo causes Sigatoka disease of banana and is endemic to Australia. The population genetic structure of M. musicola in Australia was examined by applying single-copy restriction fragment length polymorphism probes to hierarchically sampled populations collected along the Australian cast coast. The 363 isolates studied were from 16 plantations at 12 sites in four different regions, and comprised 11 populations. These populations displayed moderate levels of gene diversity (H = 0.142 to 0.369) and similar levels of genotypic richness and evenness. Populations were dominated by unique genotypes, but isolates sharing the same genotype (putative clones) were detected. Genotype distribution was highly localized within each population, and the majority of putative clones were detected for isolates sampled from different sporodochia in the same lesion or different lesions on a plant. Multilocus gametic disequilibrium tests provided further evidence of a degree of clonality within the populations at the plant scale. A complex pattern of population differentiation was detected for M. musicola in Australia. Populations sampled from plantations outside the two major production areas were genetically very different to all other populations. Differentiation was much lower between populations of the two major production areas, despite their geographic separation of over 1,000 km. These results suggest low gene flow at the continental scale due to limited spore dispersal and the movement of infected plant material.

Optimal sampling strategy for estimation of spatial genetic structure in tree populations

Fine-scale spatial genetic structure (SGS) in natural tree populations is largely a result of restricted pollen and seed dispersal. Understanding the link between limitations to dispersal in gene vectors and SGS is of key interest to biologists and the availability of highly variable molecular markers has facilitated fine-scale analysis of populations. However, estimation of SGS may depend strongly on the type of genetic marker and sampling strategy (of both loci and individuals). To explore sampling limits, we created a model population with simulated distributions of dominant and codominant alleles, resulting from natural regeneration with restricted gene flow. SGS estimates from subsamples (simulating collection and analysis with amplified fragment length polymorphism (AFLP) and microsatellite markers) were correlated with the 'real' estimate (from the full model population). For both marker types, sampling ranges were evident, with lower limits below which estimation was poorly correlated and upper limits above which sampling became inefficient. Lower limits (correlation of 0.9) were 100 individuals, 10 loci for microsatellites and 150 individuals, 100 loci for AFLPs. Upper limits were 200 individuals, five loci for microsatellites and 200 individuals, 100 loci for AFLPs. The limits indicated by simulation were compared with data sets from real species. Instances where sampling effort had been either insufficient or inefficient were identified. The model results should form practical boundaries for studies aiming to detect SGS. However, greater sample sizes will be required in cases where SGS is weaker than for our simulated population, for example, in species with effective pollen/seed dispersal mechanisms.

Genetic differentiation in Scottish populations of the pine beauty moth, Panolis flammea (Lepidoptera : Noctuidae)

Pine beauty moth, Panolis flammea (Denis & Schiffermuller), is a recent but persistent pest of lodgepole pine plantations in Scotland, but exists naturally at low levels within remnants and plantations of Scots pine. To test whether separate host races occur in lodgepole and Scots pine stands and to examine colonization dynamics, allozyme, randomly amplified polymorphic DNA (RAPD) and mitochondrial variation were screened within a range of Scottish samples. RAPD analysis indicated limited long distance dispersal (F-ST=0.099), and significant isolation by distance (P < 0.05); but that colonization between more proximate populations was often variable, from extensive to limited exchange. When compared with material from Germany, Scottish samples were found to be more diverse and significantly differentiated for all markers. For mtDNA, two highly divergent groups of haplotypes were evident, one group contained both German and Scottish samples and the other was predominantly Scottish. No genetic differentiation was evident between P. flammea populations sampled from different hosts, and no diversity bottleneck was observed in the lodgepole group. Indeed, lodgepole stands appear to have been colonized on multiple occasions from Scots pine sources and neighbouring populations on different hosts are close to panmixia.

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.

Reduced genetic diversity and significant genetic differentiation after translocation: Comparison of the remnant and translocated populations of bridled nailtail wallabies (Onychogalea fraenata)

Loss of genetic diversity and increased population differentiation from source populations are common problems associated with translocation programmes established from captive-bred stock or a small number of founders. The bridled nailtail wallaby is one of the most endangered macropods in Australia, having been reduced to a single remnant population in the last 100 years. A translocated population of bridled nailtail wallabies was established using animals sourced directly from the remnant population (wild-released) as well as the progeny of animals collected for a captive breeding programme (captive-bred). The aims of this study were to compare genetic diversity among released animals and their wild-born progeny to genetic diversity observed in the remnant population, and to monitor changes in genetic diversity over time as more animals were released into the population. Heterozygosity did not differ between the translocated and remnant population; however, allelic diversity was significantly reduced across all released animals and their wild-born progeny. Animals bred in captivity and their wild-born progeny were also significantly differentiated from the source population after just four generations. Wild-released animals, however, were representative of the source population and several alleles were unique to this group. Both heterozygosity and allelic diversity among translocated animals decreased over time with the additional release of captive-bred animals, as no new genetic stock was added to the population. Captive breeding programmes can provide large numbers of animals for release, but this study highlights the importance of sourcing animals directly from remnant populations in order to maintain genetic diversity and minimise genetic drift.

Genetic variation of the scleractinian coral Stylophora pistillata, from western Pacific reefs

Most scleractinian coral species are widely distributed across the tropical and subtropical Indo-Pacific. However, the genetic connectivity between populations of corals separated by large distances (thousands of kilometers) is not well known. We analyzed variability in the nucleotide sequence of the internal transcribed spacer-1 (ITS-1) of the nuclear ribosomal gene unit in the ubiquitous coral Stylophora pistillata, across the western Pacific Ocean. Eight populations from Japan, Malaysia, and the northern and southern Great Barrier Reef (GBR) were studied. Phylogenetic analyses and analysis of molecular variance (AMOVA) clearly revealed that there is panmixia among these coral populations. AMOVA showed that ITS-1 sequence variability was greater within populations (78.37%) than among populations (12.06%). These patterns strongly suggest high levels of connectivity across the species' latitudinal distribution range in the western Pacific, as is seen in many marine invertebrates.

The genetic covariance among clinal environments after adaptation to an environmental gradient in Drosophila serrata

We examined the genetic basis of clinal adaptation by determining the evolutionary response of life-history traits to laboratory natural selection along a gradient of thermal stress in Drosophila serrata. A gradient of heat stress was created by exposing larvae to a heat stress of 36degrees for 4 hr for 0, 1, 2, 3, 4, or 5 days of larval development, with the remainder of development taking place at 25degrees. Replicated lines were exposed to each level of this stress every second generation for 30 generations. At the end of selection, we conducted a complete reciprocal transfer experiment where all populations were raised in all environments, to estimate the realized additive genetic covariance matrix among clinal environments in three life-history traits. Visualization of the genetic covariance functions of the life-history traits revealed that the genetic correlation between environments generally declined as environments became more different and even became negative between the most different environments in some cases. One exception to this general pattern was a life-history trait representing the classic trade-off between development time and body size, which responded to selection in a similar genetic fashion across all environments. Adaptation to clinal environments may involve a number of distinct genetic effects along the length of the cline, the complexity of which may not be fully revealed by focusing primarily on populations at the ends of the cline.

Factors supporting the non-persistence of fruit fly populations in South Australia

For purposes of interstate and international fruit trade, it is necessary to demonstrate that in areas in which fruit fly species have not previously established permanent populations, but which are subject to introductions of fruit flies from outside the area, the introduced population once detected, has not become established. In this paper, we apply methodology suggested mainly by Carey (1991, 1995) to introductions of Mediterranean fruit fly (Medfly), Ceratitis capitata Weid., and Queensland fruit fly (QFF) Bactrocera tryoni Froggatt (Diptera: Tephritidae) to South Australia, a state in which these species do not occur naturally and in which introductions, once detected, are actively treated. By analysing historical data associated with fruit fly outbreaks in South Australia, we demonstrate that: (i) fruit flies occur seasonally, as would occur in established populations, except there is no evidence of the critical spring generation of either species; (ii) there is no evidence of increasing frequency of outbreaks, trapped flies or larval occurrences over 29 years; (iii) there is no evidence of decreasing time between catches of adult flies as the years progress; (iv) there is no decrease in the mean number of years between outbreaks in the same locations; (v) there is no statistically significant recurrence of outbreaks in the same locations in successive years; (vi) there is no evidence of spread of outbreaks outwards from a central location; (vii) the likelihood of outbreaks in a city or town is related to the size of the human population; (viii) introduction pathways by road from Western Australia (for Medfly) and eastern Australia (for QFF) are shown to exist and to illegally or accidentally carry considerable amounts of fruit into South Australia; and (ix) there was no association between the numbers of either Queensland fruit fly or Medfly and the spatial pattern of either loquat or cumquat trees as sources of larval food in spring. This analysis supports the hypothesis that most fruit fly outbreaks in South Australia have been the result of separate introductions of infested fruit by vehicular traffic and that most of the resultant fly outbreaks were detected and died out within a few weeks of the application of eradication procedures. An alternative hypothesis, that populations of fruit flies are established in South Australia at below detectable levels, is impossible to disprove with conventional technology, but the likelihood of it being true is minimised by our analysis. Both hypotheses could be tested soon with newly developed genetic techniques.

Targeting genetic resource conservation in widespread species: a case study of Cedrela odorata L.

Conservation of genetic resources is a recognised necessity for the long term maintenance of evolutionary potential. Effective assessment and implementation Strategies are required to permit rapid evaluation and protection of resources. Here we use information from the chloroplast, total genome and quantitative characters assayed across wide-ranging populations to assess genetic resources in a Neotropical tree, Cedrela odorata. A major differentiation identified for organelle, total genomic and quantitative variation was found to coincide with an environmental gradient across Costa Rica. However, a major evolutionary divergence between the Yucatan region and Honduras/Nicaragua identified within the chloroplast genome was not differentiated using quantitative characters. Based on these and other results, a three-tiered conservation genetic prioritisation process is recommended. In order of importance, and where information is available, conservation units should be defined using quantitative (expressed genes), nuclear (genetic connectivity) and organellar (evolutionary) measures. Where possible, information from range wide and local scale studies should be combined and emphasis should be placed on coincidental disjunctions for two or more measures. However, if only rapid assessments of diversity are possible, then assessment of organelle variation provides the most cautious assessment of genetic resources, at least for C. odorata, and can be used to propose initial conservation units. When considering effective implementation of genetic resource management strategies a final tier should be considered, that of landuse/geopolitical divisions. (C) 2004 Elsevier B.V. All rights reserved.

Genetic effects of chronic habitat fragmentation on tree species: the case of Sorbus aucuparia in a deforested Scottish landscape

Sustainable forest restoration and management practices require a thorough understanding of the influence that habitat fragmentation has on the processes shaping genetic variation and its distribution in tree populations. We quantified genetic variation at isozyme markers and chloroplast DNA (cpDNA), analysed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) in severely fragmented populations of Sorbus aucuparia (Rosaceae) in a single catchment (Moffat) in southern Scotland. Remnants maintain surprisingly high levels of gene diversity (H-E) for isozymes (H-E = 0.195) and cpDNA markers (H-E = 0.490). Estimates are very similar to those from non-fragmented populations in continental Europe, even though the latter were sampled over a much larger spatial scale. Overall, no genetic bottleneck or departures from random mating were detected in the Moffat fragments. However, genetic differentiation among remnants was detected for both types of marker (isozymes Theta(n) = 0.043, cpDNA Theta(c) = 0.131; G-test, P-value < 0.001). In this self-incompatible, insect-pollinated, bird-dispersed tree species, the estimated ratio of pollen flow to seed flow between fragments is close to 1 (r = 1.36). Reduced pollen-mediated gene flow is a likely consequence of habitat fragmentation, but effective seed dispersal by birds is probably helping to maintain high levels of genetic diversity within remnants and reduce genetic differentiation between them.

Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power

Genetic assignment methods use genotype likelihoods to draw inference about where individuals were or were not born, potentially allowing direct, real-time estimates of dispersal. We used simulated data sets to test the power and accuracy of Monte Carlo resampling methods in generating statistical thresholds for identifying F-0 immigrants in populations with ongoing gene flow, and hence for providing direct, real-time estimates of migration rates. The identification of accurate critical values required that resampling methods preserved the linkage disequilibrium deriving from recent generations of immigrants and reflected the sampling variance present in the data set being analysed. A novel Monte Carlo resampling method taking into account these aspects was proposed and its efficiency was evaluated. Power and error were relatively insensitive to the frequency assumed for missing alleles. Power to identify F-0 immigrants was improved by using large sample size (up to about 50 individuals) and by sampling all populations from which migrants may have originated. A combination of plotting genotype likelihoods and calculating mean genotype likelihood ratios (D-LR) appeared to be an effective way to predict whether F-0 immigrants could be identified for a particular pair of populations using a given set of markers.