Pronounced genetic population structure in a potentially vagile fish species (Pristipomoides multidens, Teleostei; Perciformes; Lutjanidae) from the East Indies triangle

The East Indies triangle, bordered by the Phillipines, Malay Peninsula and New Guinea, has a high level of tropical marine species biodiversity. Pristipomoides multidens is a large, long-lived, fecund snapper species that is distributed throughout the East Indies and Indo-Pacific. Samples were analysed from central and eastern Indonesia and northern Australia to test for genetic discontinuities in population structure. Fish (n = 377) were collected from the Indonesian islands of Bali, Sumbawa, Flores, West Timor, Tanimbar and Tual along with 131 fish from two northern Australian locations (Arafura and Timor Seas) from a previous study. Genetic variation in the control region of the mitochondrial genome was assayed using restriction fragment length polymorphism and direct sequencing. Haplotype diversity was high (0.67-0.82), as was intraspecific sequence divergence (range 0-5.8%). FST between pairs of populations ranged from 0 to 0.2753. Genetic subdivision was apparent on a small spatial scale; FST was 0.16 over 191 km (Bali/Sumbawa) and 0.17 over 491 km (Bali/Flores). Constraints to dispersal that contribute to, and maintain, the observed degree of genetic subdivision are experienced presumably by all life history stages of this tropical marine finfish. The constraints may include (1) little or no movement of eggs or larvae, (2) little or no home range or migratory movement of adults and (3) loss of larval cohorts due to transport of larvae away from suitable habitat by prevailing currents

Phylogeography of the green turtle, Chelonia mydas, in the Southwest Indian Ocean

Patterns of mitochondrial DNA (mtDNA) variation were used to analyse the population genetic structure of southwestern Indian Ocean green turtle (Chelonia mydas) populations. Analysis of sequence variation over 396 bp of the mtDNA control region revealed seven haplotypes among 288 individuals from 10 nesting sites in the Southwest Indian Ocean. This is the first time that Atlantic Ocean haplotypes have been recorded among any Indo-Pacific nesting populations. Previous studies indicated that the Cape of Good Hope was a major biogeographical barrier between the Atlantic and Indian Oceans because evidence for gene flow in the last 1.5 million years has yet to emerge. This study, by sampling localities adjacent to this barrier, demonstrates that recent gene flow has occurred from the Atlantic Ocean into the Indian Ocean via the Cape of Good Hope. We also found compelling genetic evidence that green turtles nesting at the rookeries of the South Mozambique Channel (SMC) and those nesting in the North Mozambique Channel (NMC) belong to separate genetic stocks. Furthermore, the SMC could be subdivided in two different genetic stocks, one in Europa and the other one in Juan de Nova. We suggest that this particular genetic pattern along the Mozambique Channel is attributable to a recent colonization from the Atlantic Ocean and is maintained by oceanic conditions in the northern and southern Mozambique Channel that influence early stages in the green turtle life cycle.

Characterization of 26 new microsatellite loci in the dugong (Dugong dugon)

Twenty-six microsatellite loci have been isolated from a dugong (Dugong dugon). The average heterozygosity was 0.52 with two to 10 alleles per locus surveyed from 50 individuals. The markers are suitable for genetic mark-recapture (PID = 5 × 10-16) in dugongs and they could also be used to quantify physical tag loss, estimate relatedness, assign paternity, elucidate population structure and identify migrants.

The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus)

This study compares estimates of the census size of the spawning population with genetic estimates of effective current and long-term population size for an abundant and commercially important marine invertebrate, the brown tiger prawn (Penaeus esculentus). Our aim was to focus on the relationship between genetic effective and census size that may provide a source of information for viability analyses of naturally occurring populations. Samples were taken in 2001, 2002 and 2003 from a population on the east coast of Australia and temporal allelic variation was measured at eight polymorphic microsatellite loci. Moments-based and maximum-likelihood estimates of current genetic effective population size ranged from 797 to 1304. The mean long-term genetic effective population size was 9968. Although small for a large population, the effective population size estimates were above the threshold where genetic diversity is lost at neutral alleles through drift or inbreeding. Simulation studies correctly predicted that under these experimental conditions the genetic estimates would have non-infinite upper confidence limits and revealed they might be overestimates of the true size. We also show that estimates of mortality and variance in family size may be derived from data on average fecundity, current genetic effective and census spawning population size, assuming effective population size is equivalent to the number of breeders. This work confirms that it is feasible to obtain accurate estimates of current genetic effective population size for abundant Type III species using existing genetic marker technology.

New microsatellite loci for Carcharhinid sharks (Carcharhinus tilstoni and C. sorrah) and their cross-amplification in other shark species

Twelve microsatellite DNA markers were isolated in the spot-tail shark (Carcharhinus sorrah) and nine were isolated in Australian black-tip shark (Carcharhinus tilstoni). These loci plus 18 others developed for sharks from the genera Negaprion, Ginglymostoma, Carcharodon and Isurus were tested for amplification success on four species of Carcharhinus (including C. sorrah and C. tilstoni) and four other species representing three diverse families. Cross-amplification was most common within families. Five loci were subsequently tested for polymorphism on 50 C. sorrah and 60 C. tilstoni. The number of alleles per locus was two to 24 and the average heterozygosity was 0.54 (range 0.16-0.87) for C. sorrah and 0.64 (range 0.44-0.78) for C. tilstoni. These loci may be useful tools for genetic analyses of the Carcharhinidae.

The genetic structure of Australasian green turtles (Chelonia mydas): Exploring the geographical scale of genetic exchange

Ecological and genetic studies of marine turtles generally support the hypothesis of natal homing, but leave open the question of the geographical scale of genetic exchange and the capacity of turtles to shift breeding sites. Here we combine analyses of mitochondrial DNA (mtDNA) variation and recapture data to assess the geographical scale of individual breeding populations and the distribution of such populations through Australasia. We conducted multiscale assessments of mtDNA variation among 714 samples from 27 green turtle rookeries and of adult female dispersal among nesting sites in eastern Australia. Many of these rookeries are on shelves that were flooded by rising sea levels less than 10 000 years (c. 450 generations) ago. Analyses of sequence variation among the mtDNA control region revealed 25 haplotypes, and their frequency distributions indicated 17 genetically distinct breeding stocks (Management Units) consisting either of individual rookeries or groups of rookeries in general that are separated by more than 500 km. The population structure inferred from mtDNA was consistent with the scale of movements observed in long-term mark-recapture studies of east Australian rookeries. Phylogenetic analysis of the haplotypes revealed five clades with significant partitioning of sequence diversity (Φ = 68.4) between Pacific Ocean and Southeast Asian/Indian Ocean rookeries. Isolation by distance was indicated for rookeries separated by up to 2000 km but explained only 12% of the genetic structure. The emerging general picture is one of dynamic population structure influenced by the capacity of females to relocate among proximal breeding sites, although this may be conditional on large population sizes as existed historically across this region.

Polymorphic microsatellite loci for the zebra shark Stegostoma fasciatum

We report on the development and characterization of 14 polymorphic microsatellite loci in the zebra shark (Stegostoma fasciatum). Five tetranucleotide and nine dinucleotide loci were polymorphic with heterozygosities ranging from 0.400 to 0.967 and from three to 22 alleles per locus. Cross-species amplification of these zebra shark primers on four other species of orectolobid sharks was not successful.

A PCR assay for gender assignment in dugong (Dugong dugon) and West Indian manatee (Trichechus manatus)

Gender assignment for some aquatic mammals in the field is difficult. Molecular sexing from tissue biopsies is possible as males are heterogametic. Here we describe a multiplex PCR assay that amplifies the male specific SRY gene and differentiates ZFX and ZFY gametologues in two sirenian species, dugong (Dugong dugon) and West Indian manatee (Trichechus manatus). The assay was validated with animals of known gender and proved accurate and robust to experimental failure.

ShadowBoxer and LocusEater: programs to optimise experimental design and multiplexing strategies for genetic mark-recapture

Genetic mark–recapture requires efficient methods of uniquely identifying individuals. 'Shadows' (individuals with the same genotype at the selected loci) become more likely with increasing sample size, and bias harvest rate estimates. Finding loci is costly, but better loci reduce analysis costs and improve power. Optimal microsatellite panels minimize shadows, but panel design is a complex optimization process. locuseater and shadowboxer permit power and cost analysis of this process and automate some aspects, by simulating the entire experiment from panel design to harvest rate estimation.

IUCN classification zones concord with, but underestimate, the population genetic structure of the zebra shark Stegostoma fasciatum in the Indo-West Pacific

The Indo-West Pacific (IWP), from South Africa in the western Indian Ocean to the western Pacific Ocean, contains some of the most biologically diverse marine habitats on earth, including the greatest biodiversity of chondrichthyan fishes. The region encompasses various densities of human habitation leading to contrasts in the levels of exploitation experienced by chondrichthyans, which are targeted for local consumption and export. The demersal chondrichthyan, the zebra shark, Stegostoma fasciatum, is endemic to the IWP and has two current regional International Union for the Conservation of Nature (IUCN) Red List classifications that reflect differing levels of exploitation: ‘Least Concern’ and ‘Vulnerable’. In this study, we employed mitochondrial ND4 sequence data and 13 microsatellite loci to investigate the population genetic structure of 180 zebra sharks from 13 locations throughout the IWP to test the concordance of IUCN zones with demographic units that have conservation value. Mitochondrial and microsatellite data sets from samples collected throughout northern Australia and Southeast Asia concord with the regional IUCN classifications. However, we found evidence of genetic subdivision within these regions, including subdivision between locations connected by habitat suitable for migration. Furthermore, parametric FST analyses and Bayesian clustering analyses indicated that the primary genetic break within the IWP is not represented by the IUCN classifications but rather is congruent with the Indonesian throughflow current. Our findings indicate that recruitment to areas of high exploitation from nearby healthy populations in zebra sharks is likely to be minimal, and that severe localized depletions are predicted to occur in zebra shark populations throughout the IWP region.

Characterization of highly informative cross-species microsatellite panels for the Australian dugong (Dugong dugon) and Florida manatee (Trichechus manatus latirostris) including five novel primers

The Australian dugong (Dugong dugon) and Florida manatee (Trichechus manatus latirostris) are threatened species of aquatic mammals in the order Sirenia. Sirenian conservation and management actions would benefit from a more complete understanding of genetic diversity and population structure. Generally, species-specific microsatellite markers are employed in conservation genetic studies; however, robust markers can be difficult and costly to isolate. To increase the number of available markers, dugong and manatee microsatellite primers were evaluated for cross-species amplification. Furthermore, one manatee and four dugong novel primers are reported. After polymerase chain reaction optimization, 23 (92%) manatee primers successfully amplified dugong DNA, of which 11 (48%) were polymorphic. Of the 32 dugong primers tested, 27 (84%) yielded product in the manatee, of which 17 (63%) were polymorphic. Dugong and manatee primers were compared and the most informative markers were selected to create robust and informative marker-panels for each species. These cross-species microsatellite marker-panels can be employed to assess other sirenian populations and can provide beneficial information for the protection and management of these unique mammals.

Fifteen microsatellite loci for the jungle perch, Kuhlia rupestris

We developed and optimized 15 polymorphic microsatellite loci in the jungle perch, Kuhlia rupestris. Loci were screened in a single population (n = 24) from Fraser Island, Queensland, Australia. Number of alleles per locus ranged from 3 to 19 and observed heterozygosity from 0.25 to 1. No significant linkage disequilibrium was detected between any pair of loci. Genotype proportions for these loci in the population sampled were in Hardy–Weinberg equilibrium.

The spatiotemporal dynamics of Tribolium castaneum (Herbst): adult flight and gene flow.

Tribolium castaneum (Herbst) has been used as a model organism to develop and test important ecological and evolutionary concepts and is also a major pest of grain and grain products globally. This beetle species is assumed to be a good colonizer of grain storages through anthropogenic movement of grain, and active dispersal by flight is considered unlikely. Studies using T. castaneum have therefore used confined or walking insects. We combine an ecological study of dispersal with an analysis of gene flow using microsatellites to investigate the spatiotemporal dynamics and adult flight of T. castaneum in an ecological landscape in eastern Australia. Flying beetles were caught in traps at grain storages and in fields at least 1 km from the nearest stored grain at regular intervals for an entire year. Significantly more beetles were trapped at storages than in fields, and almost no beetles were caught in native vegetation reserves many kilometres from the nearest stored grain. Genetic differentiation between beetles caught at storages and in fields was low, indicating that although T. castaneum is predominantly aggregated around grain storages, active dispersal takes place to the extent that significant gene flow occurs between them, mitigating founder effects and genetic drift. By combining ecological and molecular techniques, we reveal much higher levels of active dispersal through adult flight in T. castaneum than previously thought. We conclude that the implications of adult flight to previous and future studies on this model organism warrant consideration.

Rapid microsatellite marker development for African mahogany ( Khaya senegalensis, Meliaceae) using next-generation sequencing and assessment of its intra-specific genetic diversity.

Khaya senegalensis (African mahogany or dry-zone mahogany) is a high-value hardwood timber species with great potential for forest plantations in northern Australia. The species is distributed across the sub-Saharan belt from Senegal to Sudan and Uganda. Because of heavy exploitation and constraints on natural regeneration and sustainable planting, it is now classified as a vulnerable species. Here, we describe the development of microsatellite markers for K. senegalensis using next-generation sequencing to assess its intra-specific diversity across its natural range, which is a key for successful breeding programs and effective conservation management of the species. Next-generation sequencing yielded 93943 sequences with an average read length of 234bp. The assembled sequences contained 1030 simple sequence repeats, with primers designed for 522 microsatellite loci. Twenty-one microsatellite loci were tested with 11 showing reliable amplification and polymorphism in K. senegalensis. The 11 novel microsatellites, together with one previously published, were used to assess 73 accessions belonging to the Australian K. senegalensis domestication program, sampled from across the natural range of the species. STRUCTURE analysis shows two major clusters, one comprising mainly accessions from west Africa (Senegal to Benin) and the second based in the far eastern limits of the range in Sudan and Uganda. Higher levels of genetic diversity were found in material from western Africa. This suggests that new seed collections from this region may yield more diverse genotypes than those originating from Sudan and Uganda in eastern Africa.

NEESTIMATOR v2: re-implementation of software for the estimation of contemporary effective population size (N-e) from genetic data

NeEstimator v2 is a completely revised and updated implementation of software that produces estimates of contemporary effective population size, using several different methods and a single input file. NeEstimator v2 includes three single-sample estimators (updated versions of the linkage disequilibrium and heterozygote-excess methods, and a new method based on molecular coancestry), as well as the two-sample (moment-based temporal) method. New features include the following: (i) an improved method for accounting for missing data; (ii) options for screening out rare alleles; (iii) confidence intervals for all methods; (iv) the ability to analyse data sets with large numbers of genetic markers (10000 or more); (v) options for batch processing large numbers of different data sets, which will facilitate cross-method comparisons using simulated data; and (vi) correction for temporal estimates when individuals sampled are not removed from the population (Plan I sampling). The user is given considerable control over input data and composition, and format of output files. The freely available software has a new JAVA interface and runs under MacOS, Linux and Windows.