Breeding and fecundity of the endemic Australian gudgeon, sleepy cod Oxyeleotris lineolatus (Steindachner 1867) (Eleotridae)

Sleepy cod (Oxyeleotris lineolatus Steindachner) is a tropical species of eleotrid native to northern Australia. A related species, sand or marbled goby, is the highest priced freshwater fish in Asia, and a market for a similar fish exists in expatriate Chinese communities. Sleepy cod breed when minimum temperatures reach 24 °C for more than 3 days. During the breeding season the genital papilla is broad and flattened in females compared to the triangular papilla of males and juveniles. Spawning pairs were usually of approximately equal size. Females could spawn up to 10 times during one breeding season. Wet weather increased the frequency of spawning. Eggs were usually laid hanging from the underside of a surface. Most spawning occurred between 05:00 and 10:00 h. Females attended egg masses immediately after spawning, after which males cared for eggs until hatching, 3–5 days later. Agitation of the egg mass was essential for development. The mean number of eggs per spawning was 43 130. Larvae commenced feeding 2–5 days after hatching, on plankton from 100 to 250 m in size. A spawning trap used to collect egg masses is described. The breeding biology of sleepy cod is considered to be an adaptation to the monsoonal tropics.

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

Functionally associated molecular genetic marker map construction in perennial ryegrass (Lolium perenne L.)

A molecular marker-based map of perennial ryegrass (Lolium perenne L.) has been constructed through the use of polymorphisms associated with expressed sequence tags (ESTs). A pair-cross between genotypes from a North African ecotype and the cultivar Aurora was used to generate a two-way pseudo-testcross population. A selection of 157 cDNAs assigned to eight different functional categories associated with agronomically important biological processes was used to detect polymorphic EST–RFLP loci in the F1(NA6 × AU6) population. A comprehensive set of EST–SSR markers was developed from the analysis of 14,767 unigenes, with 310 primer pairs showing efficient amplification and detecting 113 polymorphic loci. Two parental genetic maps were produced: the NA6 genetic map contains 88 EST–RFLP and 71 EST–SSR loci with a total map length of 963 cM, while the AU6 genetic map contains 67 EST–RFLP and 58 EST–SSR loci with a total map length of 757 cM. Bridging loci permitted the alignment of homologous chromosomes between the parental maps, and a sub-set of genomic DNA-derived SSRs was used to relate linkage groups to the perennial ryegrass reference map. Regions of segregation distortion were identified, in some instances in common with other perennial ryegrass maps. The EST-derived marker-based map provides the basis for in silico comparative genetic mapping, as well as the evaluation of co-location between QTLs and functionally associated genetic loci.

Isolation and characterization of novel microsatellite markers and their application for diversity assessment in cultivated groundnut (Arachis hypogaea)

Background: Cultivated peanut or groundnut (Arachis hypogaea L.) is the fourth most important oilseed crop in the world, grown mainly in tropical, subtropical and warm temperate climates. Due to its origin through a single and recent polyploidization event, followed by successive selection during breeding efforts, cultivated groundnut has a limited genetic background. In such species, microsatellite or simple sequence repeat (SSR) markers are very informative and useful for breeding applications. The low level of polymorphism in cultivated germplasm, however, warrants a need of larger number of polymorphic microsatellite markers for cultivated groundnut. Results: A microsatellite- enriched library was constructed from the genotype TMV2. Sequencing of 720 putative SSR-positive clones from a total of 3,072 provided 490 SSRs. 71.2% of these SSRs were perfect type, 13.1% were imperfect and 15.7% were compound. Among these SSRs, the GT/CA repeat motifs were the most common (37.6%) followed by GA/CT repeat motifs (25.9%). The primer pairs could be designed for a total of 170 SSRs and were optimized initially on two genotypes. 104 (61.2%) primer pairs yielded scorable amplicon and 46 (44.2%) primers showed polymorphism among 32 cultivated groundnut genotypes. The polymorphic SSR markers detected 2 to 5 alleles with an average of 2.44 per locus. The polymorphic information content (PIC) value for these markers varied from 0.12 to 0.75 with an average of 0.46. Based on 112 alleles obtained by 46 markers, a phenogram was constructed to understand the relationships among the 32 genotypes. Majority of the genotypes representing subspecies hypogaea were grouped together in one cluster, while the genotypes belonging to subspecies fastigiata were grouped mainly under two clusters. Conclusion. Newly developed set of 104 markers extends the repertoire of SSR markers for cultivated groundnut. These markers showed a good level of PIC value in cultivated germplasm and therefore would be very useful for germplasm analysis, linkage mapping, diversity studies and phylogenetic relationships in cultivated groundnut as well as related Arachis species.

Concordant molecular and phenotypic data delineate new taxonomy and conservation priorities for the endangered mountain yellow-legged frog

The mountain yellow-legged frog Rana muscosa sensu lato, once abundant in the Sierra Nevada of California and Nevada, and the disjunct Transverse Ranges of southern California, has declined precipitously throughout its range, even though most of its habitat is protected. The species is now extinct in Nevada and reduced to tiny remnants in southern California, where as a distinct population segment, it is classified as Endangered. Introduced predators (trout), air pollution and an infectious disease (chytridiomycosis) threaten remaining populations. A Bayesian analysis of 1901 base pairs of mitochondrial DNA confirms the presence of two deeply divergent clades that come into near contact in the Sierra Nevada. Morphological studies of museum specimens and analysis of acoustic data show that the two major mtDNA clades are readily differentiated phenotypically. Accordingly, we recognize two species, Rana sierrae, in the northern and central Sierra Nevada, and R. muscosa, in the southern Sierra Nevada and southern California. Existing data indicate no range overlap. These results have important implications for the conservation of these two species as they illuminate a profound mismatch between the current delineation of the distinct population segments (southern California vs. Sierra Nevada) and actual species boundaries. For example, our study finds that remnant populations of R. muscosa exist in both the southern Sierra Nevada and the mountains of southern California, which may broaden options for management. In addition, despite the fact that only the southern California populations are listed as Endangered, surveys conducted since 1995 at 225 historic (1899-1994) localities from museum collections show that 93.3% (n=146) of R. sierrae populations and 95.2% (n=79) of R. muscosa populations are extinct. Evidence presented here underscores the need for revision of protected population status to include both species throughout their ranges.

Impacts of feeding system and season on milk composition and Cheddar cheese yield in a subtropical environment

Milk obtained from cows on 2 subtropical dairy feeding systems were compared for their suitability for Cheddar cheese manufacture. Cheeses were made in a small-scale cheesemaking plant capable of making 2 blocks ( about 2 kg each) of Cheddar cheese concurrently. Its repeatability was tested over 10 separate cheesemaking days with no significant differences being found between the 2 vats in cheesemaking parameters or cheese characteristics. In the feeding trial, 16 pairs of Holstein - Friesian cows were used in 2 feeding systems (M1, rain-grown tropical grass pastures and oats; and M5, a feedlot, based on maize/barley silage and lucerne hay) over 2 seasons ( spring and autumn corresponding to early and late lactation, respectively). Total dry matter, crude protein (kg/cow. day) and metabolisable energy (MJ/cow.day) intakes were 17, 2.7, and 187 for M1 and 24, 4, 260 for M5, respectively. M5 cows produced higher milk yields and milk with higher protein and casein levels than the M1 cows, but the total solids and fat levels were similar (P > 0.05) for both M1 and M5 cows. The yield and yield efficiency of cheese produced from the 2 feeding systems were also not significantly different. The results suggest that intensive tropical pasture systems can produce milk suitable for Cheddar cheese manufacture when cows are supplemented with a high energy concentrate. Season and stage of lactation had a much greater effect than feeding system on milk and cheesemaking characteristics with autumn ( late lactation) milk having higher protein and fat contents and producing higher cheese yields.

A one-tube fluorescent assay for the quarantine detection and identification of Tilletia indica and other grass bunts in wheat

A molecular assay with enhanced specificity and sensitivity has been developed to assist in the surveillance of Karnal bunt, a quarantineable disease with a significant impact on international trade. The protocol involves the release of DNA from spores, PCR amplification to enrich Tilletia-specific templates from released DNA and a five-plex, real-time PCR assay to detect, identify and distinguish T. indica and other Tilletia species (T. walkeri, T. ehrhartae, T. horrida and a group comprising T. caries, T. laevis, T. contraversa, T. bromi and T. fusca) in wheat grains. This fluorescent molecular tool has a detection sensitivity of one spore and thus bypasses the germination step, which in the current protocol is required for confirmation when only a few spores have been found in grain samples. The assay contains five dual-labelled, species-specific probes and associated species-specific primer pairs in a PCR mix in one tube. The different amplification products are detected simultaneously by five different fluorescence spectra. This specific and sensitive assay with reduced labour and reagent requirements makes it an effective and economically sustainable tool to be used in a Karnal bunt surveillance program. This protocol will also be valuable for the identification of some contaminant Tilletia sp. in wheat grains.

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

Background: Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers. Results: The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci ( 1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions. Conclusion: The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

Scalloped Hammerhead microsatellite genotypes (eight loci) representing five regions on Australia's east coast and one location in central Indonesia.

Biodiversity of sharks in the tropical Indo-Pacific is high, but species-specific information to assist sustainable resource exploitation is scarce. The null hypothesis of population genetic homogeneity was tested for scalloped hammerhead shark (Sphyrna lewini, n=244) and the milkshark (Rhizoprionodon acutus, n=209) from northern and eastern Australia, using nuclear (S. lewini, eight microsatellite loci; R. acutus, six loci) and mitochondrial gene markers (873 base pairs of NADH dehydrogenase subunit 4). We were unable to reject genetic homogeneity for S. lewini, which was as expected based on previous studies of this species. Less expected were similar results for R. acutus, which is more benthic and less vagile than S. lewini. These features are probably driving the genetic break found between Australian and central Indonesian R. acutus (F-statistics; mtDNA, 0.751 to 0.903; microsatellite loci, 0.038 to 0.047). Our results support the spatially-homogeneous management plan for shark species in Queensland, but caution is advised for species yet to be studied.

Milkshark (Rhizoprionodon acutus) microsatellite genotype data (six loci) from samples collected from five locations in eastern Australia and central Indonesia.

Biodiversity of sharks in the tropical Indo-Pacific is high, but species-specific information to assist sustainable resource exploitation is scarce. The null hypothesis of population genetic homogeneity was tested for scalloped hammerhead shark (Sphyrna lewini, n=244) and the milkshark (Rhizoprionodon acutus, n=209) from northern and eastern Australia, using nuclear (S. lewini, eight microsatellite loci; R. acutus, six loci) and mitochondrial gene markers (873 base pairs of NADH dehydrogenase subunit 4). We were unable to reject genetic homogeneity for S. lewini, which was as expected based on previous studies of this species. Less expected were similar results for R. acutus, which is more benthic and less vagile than S. lewini. These features are probably driving the genetic break found between Australian and central Indonesian R. acutus (F-statistics; mtDNA, 0.751 to 0.903; microsatellite loci, 0.038 to 0.047). Our results support the spatially-homogeneous management plan for shark species in Queensland, but caution is advised for species yet to be studied.

Quantitative and molecular genetics of juvenile wood traits in radiata and slash/Caribbean pines.

The Juvenile Wood Initiative (JWI) project has been running successfully since July 2003 under a Research Agreement with FWPA and Letters of Association with the consortium partners STBA (Southern Tree Breeding Association), ArborGen and FPQ (Forestry Plantations Queensland). Over the last five and half years, JWI scientists in CSIRO, FPQ, and STBA have completed all 12 major milestones and 28 component milestones according to the project schedule. We have made benchmark progress in understanding the genetic control of wood formation and interrelationships among wood traits. The project has made 15 primary scientific findings and several results have been adopted by industry as summarized below. This progress was detailed in 10 technical reports to funding organizations and industry clients. Team scientists produced 16 scientific manuscripts (8 published, 1 in press, 2 submitted, and several others in the process of submission) and 15 conference papers or presentations. Primary Scientific Findings. The 15 major scientific findings related to wood science, inheritance and the genetic basis of juvenile wood traits are: 1. An optimal method to predict stiffness of standing trees in slash/Caribbean pine is to combine gravimetric basic density from 12 mm increment cores with a standing tree prediction of MoE using a time of flight acoustic tool. This was the most accurate and cheapest way to rank trees for breeding selection for slash/Caribbean hybrid pine. This method was also recommended for radiata pine. 2. Wood density breeding values were predicted for the first time in the STBA breeding population using a large sample of 7,078 trees (increment cores) and it was estimated that selection of the best 250 trees for deployment will produce wood density gains of 12.4%. 3. Large genetic variation for a suite of wood quality traits including density, MFA, spiral grain, shrinkage, acoustic and non-acoustic stiffness (MoE) for clear wood and standing trees were observed. Genetic gains of between 8 and 49% were predicted for these wood quality traits with selection intensity between 1 to 10% for radiata pine. 4. Site had a major effect on juvenile-mature wood transition age and the effect of selective breeding for a shorter juvenile wood formation phase was only moderate (about 10% genetic gain with 10% selection intensity, equivalent to about 2 years reduction of juvenile wood). 5. The study found no usable site by genotype interactions for the wood quality traits of density, MFA and MoE for both radiata and slash/Caribbean pines, suggesting that assessment of wood properties on one or two sites will provide reliable estimates of the genetic worth of individuals for use in future breeding. 6. There were significant and sizable genotype by environment interactions between the mainland and Tasmanian regions and within Tasmania for DBH and branch size. 7. Strong genetic correlations between rings for density, MFA and MoE for both radiata and slash/Caribbean pines were observed. This suggests that selection for improved wood properties in the innermost rings would also result in improvement of wood properties in the subsequent rings, as well as improved average performance of the entire core. 8. Strong genetic correlations between pure species and hybrid performance for each of the wood quality traits were observed in the hybrid pines. Parental performance can be used to identify the hybrid families which are most likely to have superior juvenile wood properties of the slash/Caribbean F1 hybrid in southeast Queensland. 9. Large unfavourable genetic correlations between growth and wood quality traits were a prominent feature in radiata pine, indicating that overcoming this unfavourable genetic correlation will be a major technical issue in progressing radiata pine breeding. 10. The project created the first radiata pine 18 k cDNA microarray and generated 5,952 radiata pine xylogenesis expressed sequence tags (ESTs) which assembled into 3,304 unigenes. 11. A total of 348 genes were identified as preferentially expressed genes in earlywood or latewood while a total of 168 genes were identified as preferentially expressed genes in either juvenile or mature wood. 12. Juvenile earlywood has a distinct transcriptome relative to other stages of wood development. 13. Discovered rapid decay of linkage disequilibrium (LD) in radiata pine with LD decaying to approximately 50% within 1,700 base pairs (within a typical gene). A total of 913 SNPS from sequencing 177,380 base pairs were identified for association genetic studies. 14. 149 SNPs from 44 genes and 255 SNPs from a further 51 genes (total 95 genes) were selected for association analysis with 62 wood traits, and 30 SNPs were shortlisted for their significant association with variation of wood quality traits (density, MFA and MoE) with individual significant SNPs accounting for between 1.9 and 9.7% of the total genetic variation in traits. 15. Index selection using breeding objectives was the most profitable selection method for radiata pine, but in the long term it may not be the most effective in dealing with negative genetic correlations between wood volume and quality traits. A combination of economic and biological approaches may be needed to deal with the strong adverse correlation.

Raw data for project "Development of a DNA based aging technique for use in fisheries assessments" Fisheries and Development Corporation project 2007/033

The primary aim of this study was to determine the relationship between telomere length and age in a range of marine invertebrates including abalone (Haliotis spp) oysters (Saccostrea glomerata), spiny lobsters (Sagmariasus verreauxi formerly Jasus verreauxi and Jasus edwardsii) and school prawns (Metapenaeus macleayi). Additionally, this relationship was studied in a vertebrate organism using the freshwater fish Silver perch (Bidyanus bidyanus). Telomere length differences between tissues were also examined in some species such as Saccostrea glomerata, Sagmariasus verreauxi and Bidyanus bidyanus. In some cases cultured specimens of known age were used and this is quoted in the spreadsheets. For other wild-caught specimens where age was not known, size was used as a proxy for age. This may be a broad size class, or be determined by shell size or carapace length depending on the organism. Each spreadsheet contains raw data of telomere length estimates from Terminal Restriction Fragment Assays (TRF) for various individuals of each species including appropriate details such as age or size and tissue. Telomere length estimates are given in base pairs (bp). In most cases replicate experiments were conducted on groups of samples three times but on a small number of occasions only two replicate experiments were conducted. Further description of the samples can be found in final report of FRDC 2007/033. The arithmetic average for each individual (sample ID) across the two or three replicate experiments is also given. Bidyanus bidyanus (SilverPerch) Two sheets are contained within. a) Comparison of telomere length between different tissues (heart, liver and muscle) within the three year old age class - two replicate experiments were conducted. b) Comparison of telomere length between fish of different but known ages (0.25, 1, 2, and 3 years old) in each of three tissues, heart, liver and muscle – three replicate experiments were conducted per tissue. Haliotis spp (Abalone species) Three species were tested. H. asinina Telomere length was compared in two age classes-11 month and 18 month old abalone using muscle tissue from the foot. Within gel-variation was also estimated using a single sample run three times on one gel (replicate experiment). H. laevigata x H. rubra hybrids Telomere length was compared in three known age classes – two, three and four years old using muscle tissue from the foot. H. rubra Telomere length was compared in a range of different sized abalone using muscle tissue from the foot. Shell size is also given for each abalone Saccostrea glomerata Three sheets are contained within the file. a) Samples came from Moreton Bay Queensland in 2007. Telomere length was compared in two tissues (gill and mantle) of oysters in three age groups (1, 3 and 4 years) b) Samples came from Moreton Bay Queensland in 2009. Telomere length was compared in three age classes using DNA from gill tissue only c) Samples came from Wallis Lake, New South Wales. Telomere length was estimated from whole body minus the shell from 1 year old oysters, gill tissue of 3 age classes (1.5 years, 3 and 4 years), mantle tissue of two age classes (3 and 4 years). Sagmariasus verreauxi (formerly Jasus verreauxi) Telomere length was estimated from abdomen tissue of puerulus, gill and muscle tissue of 3 year old, large and very large size classes of lobsters. Jasus edwardsii Telomere length was measured in two size classes of lobsters- adults of varying sizes using muscle tissue and puerulus using tissues from the abdomen minus the exoskeleton. Metapenaeus macleayi Telomere length was measured in three size classes of school prawns adults. Muscle tissue was used, minus the exoskeleton.

Negligible evidence for regional genetic population structure for two shark species Rhizoprionodon acutus (Rüppell, 1837) and Sphyrna lewini (Griffith & Smith, 1834) with contrasting biology

Biodiversity of sharks in the tropical Indo-Pacific is high, but species-specific information to assist sustainable resource exploitation is scarce. The null hypothesis of population genetic homogeneity was tested for scalloped hammerhead shark (Sphyrna lewini, n = 237) and the milk shark (Rhizoprionodon acutus, n = 207) from northern and eastern Australia, using nuclear (S. lewini, eight microsatellite loci; R. acutus, six loci) and mitochondrial gene markers (873 base pairs of NADH dehydrogenase subunit 4). We were unable to reject genetic homogeneity for S. lewini, which was as expected based on previous studies of this species. Less expected were similar results for R. acutus, which is more benthic and less vagile than S. lewini. These features are probably driving the genetic break found between Australian and central Indonesian R. acutus (F-statistics; mtDNA, 0.751–0.903, respectively; microsatellite loci, 0.038–0.047 respectively). Our results support the spatially homogeneous monitoring and management plan for shark species in Queensland, Australia.

CattleTickBase: An integrated Internet-based bioinformatics resource for Rhipicephalus (Boophilus) microplus

The Rhipicephalus micro plus genome is large and complex in structure, making it difficult to assemble a genome sequence and costly to resource the required bioinformatics. In light of this, a consortium of international collaborators was formed to pool resources to begin sequencing this genome. We have acquired and assembled genomic DNA into contigs that represent over 1.8 Gigabase pairs of DNA from gene-enriched regions of the R. micro plus genome. We also have several datasets containing transcript sequences from a number of gene expression experiments conducted by the consortium. A web-based resource was developed to enable the scientific community to access our datasets and conduct analysis through a web-based bioinformatics environment called YABI. The collective bioinformatics resource is termed CattleTickBase. Our consortium has acquired genomic and transcriptomic sequence data at approximately 0.9X coverage of the gene-coding regions of the R. microplus genome. The YABI tool will facilitate access and manipulation of cattle tick genome sequence data as the genome sequencing of R. microplus proceeds. During this process the CattleTickBase resource will continue to be updated. Published by Elsevier Ltd. on behalf of Australian Society for Parasitology Inc.

Evidence for reproductive philopatry in the bull shark Carcharhinus leucas

Reproductive philopatry in bull sharks Carcharhinus leucas was investigated by comparing mitochondrial (NADH dehydrogenase subunit 4, 797 base pairs and control region genes 837 base pairs) and nuclear (three microsatellite loci) DNA of juveniles sampled from 13 river systems across northern Australia. High mitochondrial and low microsatellite genetic diversity among juveniles sampled from different rivers (mitochondrial fST = 0.0767, P < 0.05; microsatellite FST = -0.0022, P > 0.05) supported female reproductive philopatry. Genetic structure was not further influenced by geographic distance (P > 0.05) or long-shore barriers to movement (P > 0.05). Additionally, results suggest that C. leucas in northern Australia has a long-term effective population size of 11 000-13 000 females and has undergone population bottlenecks and expansions that coincide with the timing of the last ice-ages.