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.

Fine-scale spatial genetic correlation analyses reveal strong female philopatry within a brush-tailed rock-wallaby colony in southeast Queensland

We combine spatial data on home ranges of individuals and microsatellite markers to examine patterns of fine-scale spatial genetic structure and dispersal within a brush-tailed rock-wallaby (Petrogale penicillata) colony at Hurdle Creek Valley, Queensland. Brush-tailed rock-wallabies were once abundant and widespread throughout the rocky terrain of southeastern Australia; however, populations are nearly extinct in the south of their range and in decline elsewhere. We use pairwise relatedness measures and a recent multilocus spatial autocorrelation analysis to test the hypotheses that in this species, within-colony dispersal is male-biased and that female philopatry results in spatial clusters of related females within the colony. We provide clear evidence for strong female philopatry and male-biased dispersal within this rock-wallaby colony. There was a strong, significant negative correlation between pairwise relatedness and geographical distance of individual females along only 800 m of cliff line. Spatial genetic autocorrelation analyses showed significant positive correlation for females in close proximity to each other and revealed a genetic neighbourhood size of only 600 m for females. Our study is the first to report on the fine-scale spatial genetic structure within a rock-wallaby colony and we provide the first robust evidence for strong female philopatry and spatial clustering of related females within this taxon. We discuss the ecological and conservation implications of our findings for rock-wallabies, as well as the importance of fine-scale spatial genetic patterns in studies of dispersal behaviour.

A preliminary phylogenetic analysis of the Capsalidae (Platyhelminthes : Monogenea : Monopisthocotylea) inferred from large subunit rDNA sequences

Phylogenetic relationships within the Capsalidae (Monogenea) were examined Using large subunit ribosomal DNA sequences from 17 capsalid species (representing 7 genera, 5 subfamilies), 2 outgroup taxa (Monocotylidae) plus Udonella caligorum (Udonellidae). Trees were constructed using maximum likelihood, minimum evolution and maximum parsimony algorithms. An initial tree, generated from sequences 315 bases long, Suggests that Capsalinae, Encotyllabinae, Entobdellinae and Trochopodinae are monophyletic, but that Benedeniinae is paraphyletic. Analyses indicate that Neobenedenia, currently in the Benedeniinae, should perhaps be placed in 2 separate subfamily. An additional analysis was made which omitted 3 capsalid taxa (for which only short sequences were available) and all outgroup taxa because of alignment difficulties. Sequence length increased to 693 bases and good branch support was achieved. The Benedeniinae was again paraphyletic. Higher-level classification of the Capsalidae, evolution of the Entobdellinae and issues of species identity in Neobenedenia are discussed.

Mutation scanning analysis of Marteilia sydneyi populations from different geographical locations in eastern Australia

Marteilia sydneyi (Paramyxea) is the causative agent of QX disease in oysters. In spite of the economic impact of this disease, its origin and the precise reason(s) for its apparent spread in Australian waters are not yet known. Given such knowledge gaps, investigating the population genetic structure(s) of M. sydneyi populations could provide insights into the epidemiology and ecology of the parasite and could assist in its prevention and control. In this study, single strand conformation polymorphism (SSCP)-based analysis of a region (195 bp) of the first internal transcribed spacer (ITS-1) of ribosomal DNA was employed to investigate genetic variation within and among five populations of M. sydneyi from oysters from five different locations in eastern Australia. The analysis showed the existence of a genetic variant of M. sydneyi common to the Great Sandy Strait, and the Richmond and Georges Rivers, as distinct from variants at the Pimpama and Clarence Rivers. Together with historical and other information relating to the QX disease outbreaks in eastern Australia, the molecular findings support the proposal that the parasite originated in the Great Sandy Strait and/or Richmond River and then extended southward along the coast. From a technical perspective, the study demonstrated the usefulness of SSCP as a tool to study the population genetics and epidemiology of M. sydneyi. (C) 2003 Elsevier Ltd. All rights reserved.

Inherited copper toxicosis with emphasis on copper toxicosis in Bedlington terriers

Canine copper toxicosis is an important inherited disease in Bedlington terriers, because of its high prevalence rate and similarity to human copper storage disease. It can lead to chronic liver disease and occasional haemolytic anaemia due to impaired copper excretion. The responsible gene for copper toxicosis in Bedlington terriers has been recently identified and was found not to be related to human Wilson's disease gene ATP7B. Although our understanding of copper metabolism in mammals has improved through genetic molecular technology, the diversity of gene mutation related to copper metabolism in animals will help identify the responsible genes for non-Wilsonian copper toxicoses in human. This review paper discusses our knowledge of normal copper metabolism and the pathogenesis, molecular genetics and current research into copper toxicosis in Bedlington terriers, other animals and humans. (C) 2004 Elsevier GmbH. All rights reserved.

Phenotypic divergence along lines of genetic variance

Natural populations inhabiting the same environment often independently evolve the same phenotype. Is this replicated evolution a result of genetic constraints imposed by patterns of genetic covariation? We looked for associations between directions of morphological divergence and the orientation of the genetic variance-covariance matrix (G) by using an experimental system of morphological evolution in two allopatric nonsister species of rainbow fish. Replicate populations of both Melanotaenia eachamensis and Melanotaenia duboulayi have independently adapted to lake versus stream hydrodynamic environments. The major axis of divergence (z) among all eight study populations was closely associated with the direction of greatest genetic variance (g(max)), suggesting directional genetic constraint on evolution. However, the direction of hydrodynamic adaptation was strongly associated with vectors of G describing relatively small proportions of the total genetic variance, and was only weakly associated with g(max). In contrast, divergence between replicate populations within each habitat was approximately proportional to the level of genetic variance, a result consistent with theoretical predictions for neutral phenotypic divergence. Divergence between the two species was also primarily along major eigenvectors of G. Our results therefore suggest that hydrodynamic adaptation in rainbow fish was not directionally constrained by the dominant eigenvector of G. Without partitioning divergence as a consequence of the adaptation of interest (here, hydrodynamic adaptation) from divergence due to other processes, empirical studies are likely to overestimate the potential for the major eigenvectors of G to directionally constrain adaptive evolution.

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.

Gene-flow between populations of cotton bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) is highly variable between years

Both large and small scale migrations of Helicoverpa armigera Hübner in Australia were investigated using AMOVA analysis and genetic assignment tests. Five microsatellite loci were screened across 3142 individuals from 16 localities in eight major cotton and grain growing regions within Australia, over a 38-month period (November 1999 to January 2003). From November 1999 to March 2001 relatively low levels of migration were characterized between growing regions. Substantially higher than average gene-flow rates and limited differentiation between cropping regions characterized the period from April 2001 to March 2002. A reduced migration rate in the year from April 2002 to March 2003 resulted in significant genetic structuring between cropping regions. This differentiation was established within two or three generations. Genetic drift alone is unlikely to drive genetic differentiation over such a small number of generations, unless it is accompanied by extreme bottlenecks and/or selection. Helicoverpa armigera in Australia demonstrated isolation by distance, so immigration into cropping regions is more likely to come from nearby regions than from afar. This effect was most pronounced in years with limited migration. However, there is evidence of long distance dispersal events in periods of high migration (April 2001-March 2002). The implications of highly variable migration patterns for resistance management are considered.

Monitoring genetic diversity in tropical trees with multilocus dominant markers

Since no universal codominant markers are currently available, dominant genetic markers, such as amplified fragment length polymorphism (AFLP), are valuable tools for assessing genetic diversity in tropical trees. However, the measurement of genetic diversity (H) with dominant markers depends on the frequency of null homozygotes (Q) and the fixation index (F) of populations. While Q can be estimated for AFLP loci, F is less accessible. Through a modelling approach, we show that the monolocus estimation of genetic diversity is strongly dependent on the value of F, but that the multilocus diversity estimate is surprisingly robust to variations in F. The robustness of the estimate is due to a mechanistic effect of compensation between negative and positive biases of H by different AFLP loci exhibiting contrasting frequency profiles of Q. The robustness was tested across contrasting theoretical frequency profiles of Q and verified for 10 neotropical species. Practical recommendations for the implementation of this analytical method are given for genetic surveys in tropical trees, where such markers are widely applied.

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.

Sampling tissue for DNA analysis of trees: Trunk cambium as an alternative to canopy leaves

The number of studies of tropical tree species that use molecular tools is increasing, most of which collect leaf tissue for genomic DNA extraction. In tropical trees the canopy is not only frequently inaccessible, but also, once reached, the leaf tissue is often heavily defended against herbivory by high concentrations of anti-predation compounds, which may inhibit downstream applications, particularly PCR. Cambium tissue, accessed directly from the tree trunk at ground level, offers a readily accessible resource that is less hampered by the presence of defensive chemicals than leaf tissue. Here we describe a simple method for obtaining tissue from the cambial zone for DNA extraction and test the applicability of the method in a range of tropical tree species. The method was used successfully to extract DNA from 11 species in nine families. A subset of the DNA extracts was tested in more detail and proved to be highly suitable for AFLP analysis.

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.

Genetic structure of a reef-building coral from thermally distinct environments on the Great Barrier Reef

Adaptation to localised thermal regimes is facilitated by restricted gene flow, ultimately leading to genetic divergence among populations and differences in their physiological tolerances. Allozyme analysis of six polymorphic loci was used to assess genetic differentiation between nine populations of the reef-building coral Acropora millepora over a latitudinal temperature gradient on the inshore regions of the Great Barrier Reef (GBR). Small but significant genetic differentiation indicative of moderate levels of gene flow (pairwise F-ST 0.023 to 0.077) was found between southern populations of A. millepora in cooler regions of the GBR and the warmer, central or northern GBR populations. Patterns of genetic differentiation at these putatively neutral allozyme loci broadly matched experimental variation in thermal tolerance and were consistent with local thermal regimes (warmest monthly-averages) for the A. millepora populations examined. It is therefore hypothesized that natural selection has influenced the thermal tolerance of the A. millepora populations examined and greater genetic divergence is likely to be revealed by examination of genetic markers under the direct effects of natural selection.