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.

A comparison of molecular markers for genetic analysis of macadamia

Various marker systems exist for genetic analysis of horticultural species. Isozymes were first applied to the woody perennial nut crop, macadamia, in the early 1990s. The advent of DNA markers saw the development, for macadamia, of STMS (sequence-tagged microsatellite site), RAPD (randomly amplified polymorphic DNA), and RAF (randomly amplified DNA fingerprinting). The RAF technique typically generates dominant markers, but within the dominant marker profiles, certain primers also amplify multi-allelic co-dominant markers that are suspected to be microsatellites. In this paper, we confirm this for one such marker, and describe how RAF primers can be chosen that amplify one or more putative microsatellites. This approach of genotyping anonymous microsatellite markers via RAF is designated RAMiFi (randomly amplified microsatellite fingerprinting). Several marker systems were compared for the type, amount, and cost-efficiency of the information generated, using data from published studies on macadamia. The markers were also compared for the way they clustered a common set of accessions. The RAMiFi approach was identified as the most efficient and economical. The availability of such a versatile tool offers many advantages for the genetic characterisation of horticultural species.

The importance of mating system in translocation programs: reproductive success of released male bridled nailtail wallabies

Translocation is an important tool for the conservation of species that have suffered severe range reductions. The success of a translocation should be measured not only by the survival of released animals, but by the reproductive output of individuals and hence the establishment of a self-sustaining population. The bridled nailtail wallaby is an endangered Australian macropod that suffered an extensive range contraction to a single remaining wild population. A translocated population was established and subsequently monitored over a four year period. The aim of this study was to measure the reproductive success of released males using genetic tools and to determine the factors that predicted reproductive success. Captive-bred and wild-caught animals were released and we found significant variation in male reproductive success among release groups. Variation in reproductive success was best explained by individual male weight, survival and release location rather than origin. Only 26% of candidate males were observed to sire an offspring during the study. The bridled nailtail wallaby is a sexually dimorphic, polygynous macropod and reproductive success is skewed toward large males. Males over 5800 g were six times more likely to sire an offspring than males below this weight. This study highlights the importance of considering mating system when choosing animals for translocation. Translocation programs for polygynous species should release a greater proportion of females, and only release males of high breeding potential. By maximizing the reproductive output of released animals, conservation managers will reduce the costs of translocation and increase the chance of successfully establishing a self-sustaining population. (C) 2004 Elsevier Ltd. All rights reserved.

Isolation and characterization of microsatellite loci from mother-of-millions, Bryophyllum delagoense (Crassulaceae), and its hybrid with Bryophyllum daigremontianum, 'Houghton's hybrid'

Nine microsatellite loci were developed, and are transferable, across the Madagascan succulents Bryophyllum daigremontianum, Bryophyllum delagoense (mother-of-millions) and their horticultural hybrid (Houghton's), from enriched libraries of the later two species. For B. delagoense, a tetraploid, three to 13 alleles per locus were found for native Madagascan (H-O = 0.4-1.0), and one to nine in invasive Australian (H-O = 0.0-1.0) samples. In addition for 91 Australian samples, only five multilocus genotypes were found (95% of individuals were of two genotypes), suggesting extensive clonality in its introduced range. These loci will be used to examine genetic diversity, hybrid origin and mating system in natural and introduced populations.

Historical and contemporary mating patterns in remnant populations of the forest tree Fraxinus excelsior L.

Genetic variation at microsatellite markers was used to quantify genetic structure and mating behavior in a severely fragmented population of the wind-pollinated, wind-dispersed temperate tree Fraxinus excelsior in a deforested catchment in Scotland. Remnants maintain high levels of genetic diversity, comparable with those reported for continuous populations in southeastern Europe, and show low interpopulation differentiation (Theta = 0.080), indicating that historical gene exchange has not been limited (Nm = 3.48). We estimated from seeds collected from all trees producing fruits in three of five remnants that F. excelsior is predominantly outcrossing (t(m). = 0.971 +/- 0.028). Use of a neighborhood model approach to describe the relative contribution of local and long-distance pollen dispersal indicates that pollen gene flow into each of the three remnants is extensive (46-95%) and pollen dispersal has two components. The first is very localized and restricted to tens of meters around the mother trees. The second is a long-distance component with dispersal occurring over several kilometers. Effective dispersal distances, accounting for the distance and directionality to mother trees of sampled pollen donors, average 328 m and are greater than values reported for a continuous population. These results suggest that the opening of the landscape facilitates airborne pollen movement and may alleviate the expected detrimental genetic effects of fragmentation.

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.

Lack of association between genetic markers on chromosome 16q22-Q24 and type 1 diabetes in Russian affected families

Aim To evaluate whether the T1D susceptibility locus on chromosome 16q contributes to the genetic susceptibility to T1D in Russian patients. Method Thirteen microsatellite markers, spanning a 47-centimorgan genomic region on 16q22-q24 were evaluated for linkage to T1D in 98 Russian multiplex families. Multipoint logarithm of odds (LOD) ratio (MLS) and nonparametric LOD (NPL) values were computed for each marker, using GENEHUNTER 2.1 software. Four microsatellites (D16S422, D16S504, D16S3037, and D16S3098) and 6 biallelic markers in 2 positional candidate genes, ICSBP1 and NQO1, were additionally tested for association with T1D in 114 simplex families, using transmission disequilibrium test (TDT). Results A peak of linkage (MLS = 1.35, NPL = 0.91) was shown for marker D16S750, but this was not significant (P = 0.18). The subsequent linkage analysis in the subset of 46 multiplex families carrying a common risk HLA-DR4 haplotype increased peak MLS and NPL values to 1.77 and 1.22, respectively, but showed no significant linkage (P = 0.11) to T1D in the 16q22-q24 genomic region. TDT analysis failed to find significant association between these markers and disease, even after the conditioning for the predisposing HLA-DR4 haplotype. Conclusion Our results did not support the evidence for the susceptibility locus to T1D on chromosome 16q22-24 in the Russian family data set. The lack of association could reflect genetic heterogeneity of type 1 diabetes in diverse ethnic groups.

The TAF5L gene on chromosome 1q42 is associated with type 1 diabetes in Russian affected patients

Type 1 diabetes (T1D) is a multifactorial autoimmune disease, with strong genetic component. Several susceptibility loci contribute to genetic predisposition to T1D. One of these loci have been mapped to chromosome 1q42 in UK and US joined affected family data sets but needs to be replicated in other populations. In this study, we evaluated sixteen microsatellites located on 1q42 for linkage with T1D in 97 Russian affected sibling pairs. A 2.7-cm region of suggestive linkage to T1D between markers D1S1644 and D1S225 was found by multipoint linkage analysis. The peak of linkage was shown for D1S2847 (P = 0.0005). Transmission disequilibrium test showed significant undertransmission of the 156-bp allele of D1S2847 from parents to diabetic children (28 transmissions vs. 68 nontransmissions, P = 0.043) in Russian affected families. A preferential transmission from parents to diabetic offspring was also shown for the T(-25) and T1362 alleles of the C/T(-25) and C/T1362 dimorphisms, both located at the TAF5L gene, which is situated 103 kb from D1S2847. Together with the A/C744 TAF5L SNP, these markers share common T(-25)/A744/T1362 and C(-25)/C744/T1362 haplotypes associated with higher and lower risk of diabetes (Odds Ratio = 2.15 and 0.62, respectively). Our results suggest that the TAF5L gene, encoding TAF5L-like RNA polymerase II p300/CBP associated factor (PCAF)-associated factor, could represent the susceptibility gene for T1D on chromosome 1q42 in Russian affected patients.

Low level microsatellite instability may be associated with reduced cancer specific survival in sporadic stage C colorectal carcinoma

Background: Colorectal cancers (CRCs) may be categorised according to the degree of microsatellite instability (MSI) exhibited, as MSI-high (MSI-H), MSI-low (MSI-L), or microsatellite stable (MSS). MSI-H status confers a survival advantage to patients with sporadic CRC. Aims: To determine if low levels of MSI are related to the clinicopathological features and prognosis of sporadic stage C CRC. Patients: A total of 255 patients who underwent resection for sporadic stage C CRC were studied. No patient received chemotherapy. Minimum follow up was five years. Methods: DNA extracted from archival malignant and non-malignant tissue was amplified by polymerase chain reaction using a panel of 11 microsatellites. MSI-H was defined as instability at greater than or equal to40% of markers, MSS as no instability, and MSI-L as instability at >0% but,40% of markers. Patients with MSI-H CRC were excluded from analysis as they have previously been shown to have better survival. Results: Thirty three MSI-L and 176 MSS CRCs were identified. There was no difference in biological characteristics or overall survival of MSI-L compared with MSS CRC but MSI-L was associated with poorer cancer specific survival (hazard ratio 2.0 (95% confidence interval 1.1-3.6)). Conclusions: Sporadic MSI-L and MSS CRCs have comparable clinicopathological features. Further studies are required to assess the impact of MSI-L on prognosis.

Application of DNA parentage analyses for determining relative growth rates of Penaeus japonicus families reared in commercial ponds

The ability to track large numbers of individuals and families is a key determinant of the power and precision of breeding programs, including the capacity to quantify interactions between genotypes and their environment. Until recently, most family based selective breeding programs for shrimp, and other highly fecund aquaculture species, have been restricted by the number of animals that can be physically tagged and individually selected. Advances in the development of molecular markers, such as microsatellite loci, are now providing the means to track large numbers of individuals and families in commercial production systems. In this study microsatellites, coupled with DNA parentage analyses, were used to determine the relative performance of 22 families of R japonicus reared in commercial production ponds. In the experimental design 6000 post-larvae from each of 22 families, whose maternal parents had been genotyped at 8 microsatellite loci, were stocked into each of four I ha ponds. After 6 months the ponds were harvested and a total of 6000 individuals were randomly weighed from each pond. Mean wet weight of the shrimp from one pond was significantly lower than that of the other three ponds demonstrating a possible pond effect on growth rate. The representation of families in the top 10% of each pond's weight distribution was then determined by randomly genotyping up to 300 individuals from this upper weight class. Parentage analyses based on individual genotypic data demonstrated that some families were over-represented in the top 10% in all ponds, while others were under-represented due to slower growth rates. The results also revealed some weak, but significant, male genotype x environment (G x E) interactions in the expression of shrimp growth for some families. This indicates that G x E effects may need to be factored into future R japonicus selective breeding programs. This study demonstrated the utility of DNA parentage analyses for tracking individual family performance in communally stocked shrimp pond populations and, its application to examining G x E effects on trait expression under commercial culture conditions. Crown Copyright (c) 2005 Published by Elsevier B.V. All rights reserved.

Microsatellite markers for the pest fruit fly, Bactrocera papayae (Diptera : Tephritidae) and other Bactrocera species

The dorsalis complex contains some of the most economically important fruit fly pests of the Asia-Pacific regions, including Bactrocera dorsalis, Bactrocera papayae and Bactrocera carambolae. These species are morphologically indistinct and genetically very similar. We describe the development of 12 microsatellite markers isolated from a representative of the dorsalis complex, B. papayae. We show the potential utility of the B. papayae microsatellites and a set of microsatellites isolated from Bactrocera tryoni as population and species markers for the dorsalis complex.

Restricted mating dispersal and strong breeding group structure in a mid-sized marsupial mammal (Petrogale penicillata)

Ecological genetic studies have demonstrated that spatial patterns of mating dispersal, the dispersal of gametes through mating behaviour, can facilitate inbreeding avoidance and strongly influence the structure of populations, particularly in highly philopatric species. Elements of breeding group dynamics, such as strong structuring and sex-biased dispersal among groups, can also minimize inbreeding and positively influence levels of genetic diversity within populations. Rock-wallabies are highly philopatric mid-sized mammals whose strong dependence on rocky terrain has resulted in series of discreet, small colonies in the landscape. Populations show no signs of inbreeding and maintain high levels of genetic diversity despite strong patterns of limited gene flow within and among colonies. We used this species to investigate the importance of mating dispersal and breeding group structure to inbreeding avoidance within a 'small' population. We examined the spatial patterns of mating dispersal, the extent of kinship within breeding groups, and the degree of relatedness among brush-tailed rock-wallaby breeding pairs within a colony in southeast Queensland. Parentage data revealed remarkably restricted mating dispersal and strong breeding group structuring for a mid-sized mammal. Breeding groups showed significant levels of female kinship with evidence of male dispersal among groups. We found no evidence for inbreeding avoidance through mate choice; however, anecdotal data suggest the importance of life history traits to inbreeding avoidance between first-degree relatives. We suggest that the restricted pattern of mating dispersal and strong breeding group structuring facilitates inbreeding avoidance within colonies. These results provide insight into the population structure and maintenance of genetic diversity within colonies of the threatened brush-tailed rock-wallaby.

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.