Population structure, effective population size and adverse effects of stocking in the endangered Australian eastern freshwater cod Maccullochella ikei

Microsatellite markers were used to examine spatio-temporal genetic variation in the endangered eastern freshwater cod Maccullochella ikei in the Clarence River system, eastern Australia. High levels of population structure were detected. A model-based clustering analysis of multilocus genotypes identified four populations that were highly differentiated by F-statistics (FST = 0· 09 − 0· 49; P < 0· 05), suggesting fragmentation and restricted dispersal particularly among upstream sites. Hatchery breeding programmes were used to re-establish locally extirpated populations and to supplement remnant populations. Bayesian and frequency-based analyses of hatchery fingerling samples provided evidence for population admixture in the hatchery, with the majority of parental stock sourced from distinct upstream sites. Comparison between historical and contemporary wild-caught samples showed a significant loss of heterozygosity (21%) and allelic richness (24%) in the Mann and Nymboida Rivers since the commencement of stocking. Fragmentation may have been a causative factor; however, temporal shifts in allele frequencies suggest swamping with hatchery-produced M. ikei has contributed to the genetic decline in the largest wild population. This study demonstrates the importance of using information on genetic variation and population structure in the management of breeding and stocking programmes, particularly for threatened species.

Raw data for "Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson)."

Abstract of Macbeth, G. M., Broderick, D., Buckworth, R. & Ovenden, J. R. (In press, Feb 2013). Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson). G3: Genes, Genomes and Genetics. Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. But, pitfalls are predicted when the effective size is large, as estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes around 5000. This was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest standard deviation of Ne estimates occurred when low frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or non-breeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were above 24,000.

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.

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.

Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes

Despite international protection of white sharks Carcharodon carcharias, important conservation parameters such as abundance, population structure and genetic diversity are largely unknown. The tissue of 97 predominately juvenile white sharks sampled from spatially distant eastern and southwestern Australian coastlines was sequenced for the mitochondrial DNA (mtDNA) control region and genotyped with 6 nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (F-ST = 0.142, p < 0.0001), implying female reproductive philopatry. This concurs with recent satellite and acoustic tracking findings which suggest the sustained presence of discrete east coast nursery areas. Furthermore, population subdivision was found between the same regions with biparentally inherited micro satellite markers (F-ST = 0.009, p < 0.05), suggesting that males may also exhibit some degree of reproductive philopatry; 5 sharks captured along the east coast had mtDNA haplotypes that resembled western Indian Ocean sharks more closely than Australian/New Zealand sharks, suggesting that transoceanic dispersal, or migration resulting in breeding, may occur sporadically. Our most robust estimate of contemporary genetic effective population size was low and close to thresholds at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least 1, possibly 2, orders of magnitude below our historical effective size estimates. Population decline could expose these genetically isolated populations to detrimental genetic effects. Regional Australian white shark conservation management units should be implemented until genetic population structure, size and diversity can be investigated in more detail.

Extremely Low Microsatellite Diversity but Distinct Population Structure in a Long-Lived Threatened Species, the Australian Lungfish Neoceratodus forsteri (Dipnoi)

The Australian lungfish is a unique living representative of an ancient dipnoan lineage, listed as ‘vulnerable’ to extinction under Australia’s Environment Protection and Biodiversity Conservation Act 1999. Historical accounts indicate this species occurred naturally in two adjacent river systems in Australia, the Burnett and Mary. Current day populations in other rivers are thought to have arisen by translocation from these source populations. Early genetic work detected very little variation and so had limited power to answer questions relevant for management including how genetic variation is partitioned within and among sub-populations. In this study, we use newly developed microsatellite markers to examine samples from the Burnett and Mary Rivers, as well as from two populations thought to be of translocated origin, Brisbane and North Pine. We test whether there is significant genetic structure among and within river drainages; assign putatively translocated populations to potential source populations; and estimate effective population sizes. Eleven polymorphic microsatellite loci genotyped in 218 individuals gave an average within-population heterozygosity of 0.39 which is low relative to other threatened taxa and for freshwater fishes in general. Based on F_ST values (average over loci = 0.11) and STRUCTURE analyses, we identify three distinct populations in the natural range, one in the Burnett and two distinct populations in the Mary. These analyses also support the hypothesis that the Mary River is the likely source of translocated populations in the Brisbane and North Pine rivers, which agrees with historical published records of a translocation event giving rise to these populations. We were unable to obtain bounded estimates of effective population size, as we have too few genotype combinations, although point estimates were low, ranging from 29 - 129. We recommend that, in order to preserve any local adaptation in the three distinct populations that they be managed separately.

Phenology of Carpophilus spp. (Coleoptera: Nitidulidae)in stone fruit orchards as determined by pheromone trapping: Implications for prediction of crop damage

Traps baited with synthetic aggregation pheromone and fermenting bread dough were used to monitor seasonal incidence and abundance of the ripening fruit pests, Carpophilus hemipterus (L.), C. mutilatus Erichson and C. davidsoni Dobson in stone fruit orchards in the Leeton district of southern New South Wales during five seasons (1991-96). Adult beetles were trapped from September-May, but abundance varied considerably between years with the amount of rainfall in December-January having a major influence on population size and damage potential during the canning peach harvest (late February-March). Below average rainfall in December-January was associated with mean trap catches of < 10 beetles/trap/week in low dose pheromone traps during the harvest period in 1991/92 and 1993/94 and no reported damage to ripening fruit. Rainfall in December-January 1992/93 was more than double the average and mean trap catches ranged from 8-27 beetles/week during the harvest period with substantial damage to the peach crop. December-January rainfall was also above average in 1994/95 and 1995/96 and means of 50-300 beetles/trap/week were recorded in high dose pheromone traps during harvest periods. Carpophilus spp. caused economic damage to peach crops in both seasons. These data indicate that it may be possible to predict the likelihood of Carpophilus beetle damage to ripening stone fruit in inland areas of southern Australia, by routine pheromone-based monitoring of beetle populations and summer temperatures and rainfall.

A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

Frequency and precision of aerial surveys for kangaroo management

The appropriate frequency and precision for surveys of wildlife populations represent a trade-off between survey cost and the risk of making suboptimal management decisions because of poor survey data. The commercial harvest of kangaroos is primarily regulated through annual quotas set as proportions of absolute estimates of population size. Stochastic models were used to explore the effects of varying precision, survey frequency and harvest rate on the risk of quasiextinction for an arid-zone and a more mesic-zone kangaroo population. Quasiextinction probability increases in a sigmoidal fashion as survey frequency is reduced. The risk is greater in more arid regions and is highly sensitive to harvest rate. An appropriate management regime involves regular surveys in the major harvest areas where harvest rate can be set close to the maximum sustained yield. Outside these areas, survey frequency can be reduced in relatively mesic areas and reduced in arid regions when combined with lowered harvest rates. Relative to other factors, quasiextinction risk is only affected by survey precision (standard error/mean × 100) when it is >50%, partly reflecting the safety of the strategy of harvesting a proportion of a population estimate.

Seasonal incidence of Stenodiplosis sorghicola (Coquillett) (Diptera: Cecidomyiidae) and its parasitoids on Sorghum halepense (L.) Pers. in south-eastern Queensland, Australia

To quantify the role of Johnson grass, Sorghum halepense, in the population dynamics of the sorghum midge, Stenodiplosis sorghicola, patterns of flowering of Johnson grass and infestation by sorghum midge were studied in two different climatic environments in the Lockyer Valley and on the Darling Downs in south-eastern Queensland for 3 years. Parasitism levels of S. sorghicola were also recorded. In the Lockyer Valley, Johnson grass panicles were produced throughout the year but on the Darling Downs none were produced between June and September. In both areas, most panicle production occurred between November and March and infestation by S. sorghicola was the greatest during this period. The parasitism levels were between 20% and 50%. After emergence from winter diapause, one to two generations of S. sorghicola developed on S. halepense before commercial grain sorghum crops were available for infestation. Parasitoids recorded were: Aprostocetus diplosidis, Eupelmus australiensis and two species of Tetrastichus. Relationships between sorghum midge population growth rate and various environmental and population variables were investigated. Population size had a significant negative effect (P < 0.0001) on population growth rate. Mortality due to parasitism showed a significant positive density response (P < 0.0001). Temperature, rainfall, open pan evaporation, degree-days and host availability showed no significant effect on population growth rate. Given the phenology of sorghum production in south-eastern Queensland, Johnson grass provides an important bridging host, sustaining one to two generations of sorghum midge. Critical studies relating population change and build-up in sorghum to sorghum midge populations in Johnson grass are yet to be performed.

Recovery of the red-finned blue-eye: An endangered fish from springs of the Great Artesian Basin

The red-finned blue-eye (Scaturiginichthys vermeilipinnis) is endemic to a single complex of springs emanating from the Great Artesian Basin, Australia. The species has been recorded as naturally occurring in eight separate very shallow (generally <20 mm) springs, with a combined wetland area of ~0.3 ha. Since its discovery in 1990, five red-finned blue-eye (RFBE) populations have been lost and subsequent colonisation has occurred in two spring wetlands. Current population size is estimated at <3000 individuals. Artesian bores have reduced aquifer pressure, standing water levels and spring-flows in the district. There is evidence of spatial separation within the spring pools where RFBE and the introduced fish gambusia (Gambusia holbrooki) co-occur, although both species are forced together when seasonal extremes affect spring size and water temperature. Gambusia was present in four of the five springs where RFBE populations have been lost. Four out of the five remaining subpopulations of RFBE are Gambusia free. Circumstantial evidence suggests that gambusia is a major threat to red-finned blue-eyes. The impact of Gambusia is probably exacerbated by domestic stock (cattle and sheep), feral goats and pigs that utilise the springs and can negatively affect water quality and flow patterns. Three attempts to translocate RFBE to apparently suitable springs elsewhere within the complex have failed. Opportunities to mitigate threats are discussed, along with directions for future research to improve management of this extremely threatened fish and habitat.

Predicting dispersal and recruitment of Miconia calvescens (Melastomataceae) in Australian tropical rainforests

Miconia calvescens (Melastomataceae) is a serious invader in the tropical Pacific, including the Hawaiian and Tahitian Islands, and currently poses a major threat to native biodiversity in the Wet Tropics of Australia. The species is fleshy-fruited, small-seeded and shade tolerant, and thus has the potential to be dispersed widely and recruit in relatively intact rainforest habitats, displacing native species. Understanding and predicting the rate of spread is critical for the design and implementation of effective management actions. We used an individual-based model incorporating a dispersal function derived from dispersal curves for similar berry-fruited native species, and life-history parameters of fecundity and mortality to predict the spatial structure of a Miconia population after a 30 year time period. We compared the modelled population spatial structure to that of an actual infestation in the rainforests of north Queensland. Our goal was to assess how well the model predicts actual dispersion and to identify potential barriers and conduits to seed movement and seedling establishment. The model overpredicts overall population size and the spatial extent of the actual infestation, predicting individuals to occur at a maximum 1,750 m from the source compared with the maximum distance of any detected individual in the actual infestation of 1,191 m. We identify several characteristic features of managed invasive populations that make comparisons between modelled outcomes and actual infestations difficult. Our results suggest that the model’s ability to predict both spatial structure and spread of the population will be improved by incorporating a spatially explicit element, with dispersal and recruitment probabilities that reflect the relative suitability of different parts of the landscape for these processes.

New efforts at biological control of Mikania micrantha H.B.K. (Asteraceae) in Papua New Guinea and Fiji.

New efforts at biological control of Miconia calvescens (Melastomataceae) is a serious invader in the tropical Pacific, including the Hawaiian and Tahitian Islands, and currently poses a major threat to native biodiversity in the Wet Tropics of Australia. The species is fleshy-fruited, small-seeded and shade tolerant, and thus has the potential to be dispersed widely and recruit in relatively intact rainforest habitats, displacing native species. Understanding and predicting the rate of spread is critical for the design and implementation of effective management actions. We used an individual-based model incorporating a dispersal function derived from dispersal curves for similar berry-fruited native species, and life-history parameters of fecundity and mortality to predict the spatial structure of a Miconia population after a 30 year time period. We compared the modelled population spatial structure to that of an actual infestation in the rainforests of north Queensland. Our goal was to assess how well the model predicts actual dispersion and to identify potential barriers and conduits to seed movement and seedling establishment. The model overpredicts overall population size and the spatial extent of the actual infestation, predicting individuals to occur at a maximum 1,750 m from the source compared with the maximum distance of any detected individual in the actual infestation of 1,191 m. We identify several characteristic features of managed invasive populations that make comparisons between modelled outcomes and actual infestations difficult. Our results suggest that the model’s ability to predict both spatial structure and spread of the population will be improved by incorporating a spatially explicit element, with dispersal and recruitment probabilities that reflect the relative suitability of different parts of the landscape for these processes. Mikania micrantha H.B.K. (Asteraceae) in Papua New Guinea and Fiji.

Using harvest statistics to monitor temporal variation in kangaroo density and harvest rate

Management of the commercial harvest of kangaroos relies on quotas set annually as a proportion of regular estimates of population size. Surveys to generate these estimates are expensive and, in the larger states, logistically difficult; a cheaper alternative is desirable. Rainfall is a disappointingly poor predictor of kangaroo rate of increase in many areas, but harvest statistics (sex ratio, carcass weight, skin size and animals shot per unit time) potentially offer cost-effective indirect monitoring of population abundance (and therefore trend) and status (i.e. under-or overharvest). Furthermore, because harvest data are collected continuously and throughout the harvested areas, they offer the promise of more intensive and more representative coverage of harvest areas than aerial surveys do. To be useful, harvest statistics would need to have a close and known relationship with either population size or harvest rate. We assessed this using longterm (11-22 years) data for three kangaroo species (Macropus rufus, M. giganteus and M. fuliginosus) and common wallaroos (M. robustus) across South Australia, New South Wales and Queensland. Regional variation in kangaroo body size, population composition, shooter efficiency and selectivity required separate analyses in different regions. Two approaches were taken. First, monthly harvest statistics were modelled as a function of a number of explanatory variables, including kangaroo density, harvest rate and rainfall. Second, density and harvest rate were modelled as a function of harvest statistics. Both approaches incorporated a correlated error structure. Many but not all regions had relationships with sufficient precision to be useful for indirect monitoring. However, there was no single relationship that could be applied across an entire state or across species. Combined with rainfall-driven population models and applied at a regional level, these relationships could be used to reduce the frequency of aerial surveys without compromising decisions about harvest management.

Temporal genetic structure patterns in tropical maize populations under reciprocal recurrent selection.

In order to investigate the effect of long term recurrent selection on the pattern of gene diversity, thirty randomly-selected individuals from the progenitors (p) and four selection cycles (C0, C3, C6 and C11) were sampled for DNA analysis from the tropical maize (Zea mays L.) breeding populations, Atherton 1 (AT1) and Atherton 2 (AT2). Fifteen polymorphic Simple Sequence Repeat markers amplified a total of 284 and 257 alleles in AT1 and AT2 populations, respectively. Reductions in the number of alleles were observed at advanced selection cycles. About 11 and 12% of the alleles in AT1 and AT2 populations respectively, were near to fixation. However, a higher number of alleles (37% in AT1 and 33% in AT2) were close to extinction. Fisher's exact test and analysis of molecular variance (AMOVA) showed significant population differentiations. Gene diversity estimates and AMOVA revealed increased genetic differentiations at the expense of loss of heterozygosity. Population differentiations were mainly due to fixation of complementary alleles at a locus in the two breeding populations. The estimates of effective population at an advanced selection cycle were close to the population size predicted by the breeding method.