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.

Construction of microsatellite linkage maps for Corymbia

The genus Corymbia is closely related to the genus Eucalyptus, and like Eucalyptus contains tree species that are important for sub-tropical forestry. Corymbia's close relationship with Eucalyptus suggests genetic studies in Corymbia should benefit from transfer of genetic information from its more intensively studied relatives. Here we report a genetic map for Corymbia spp. based on microsatellite markers identified de novo in Corymbia sp or transferred from Eucalyptus. A framework consensus map was generated from an outbred F 2 population (n = 90) created by crossing two unrelated Corymbia torelliana x C. citriodora subsp. variegata F1 trees. The map had a total length of 367 cM (Kosambi) and was composed of 46 microsatellite markers distributed across 13 linkage groups (LOD 3). A high proportion of Eucalyptus microsatellites (90%) transferred to Corymbia. Comparative analysis between the Corymbia map and a published Eucalyptus map identified eight homeologous linkage groups in Corymbia with 13 markers mapping on one or both maps. Further comparative analysis was limited by low power to detect linkage due to low genome coverage in Corymbia, however, there was no convincing evidence for chromosomal structural differences because instances of non-synteny were associated with large distances on the Eucalyptus map. Segregation distortion was primarily restricted to a single linkage group and due to a deficit of hybrid genotypes, suggesting that hybrid inviability was one factor shaping the genetic composition of the F2 population in this inter-subgeneric hybrid. The conservation of microsatellite loci and synteny between Corymbia and Eucalyptus suggests there will be substantial value in exchanging information between the two groups.

The original native pasture ecosystems of the eastern and western Darling Downs - Can they be restored?

The original pasture ecosystems of southern inland Queensland ranged from treeless grasslands on cracking clays through grassy woodlands of varying density on a great range of soil types to those competing at the dynamic edges of forests and scrubs. Fire, both wild and aboriginal-managed, was a major factor, along with rainfall extremes, in shaping the pastures and tree:grass balance. Seedling recruitment was driven by rainfall extremes, availability of germinable seed and growing space, with seed availability and space being linked to the timing and intensity of recent fires and rain. The impact of insects, diseases, severe wind and hailstorms on recruitment should not be underestimated. The more fertile soils had denser grass growth, greater fire frequency and thinner tree cover than infertile soils, except where trees were so dense that grass growth was almost eliminated. The pastures were dominated by perennial tussock grasses of mid-height but included a wide array of minor herbaceous species whose abundance was linked to soil type and recent seasonal conditions. Many were strongly perennial with Asteraceae, Fabaceae, Malvaceae, Cyperaceae and Goodeniaceae most common in an environment, which can experience effective rainfall at any time of year. The former grassland communities that are now productive farming lands are not easily returned to their original composition. However, conservation of remnant examples of original pasture types is very achievable provided tree density is controlled, prescribed burning and grazing are used and rigorous control of invasive, exotic species is undertaken. This should be done with a clear understanding that significant short-and medium-term fluctuations in botanical composition are normal.

Wild dogma: An examination of recent "evidence" for dingo regulation of invasive mesopredator release in Australia

There is growing interest in the role that apex predators play in shaping terrestrial ecosystems and maintaining trophic cascades. In line with the mesopredator release hypothesis, Australian dingoes (Canis lupus dingo and hybrids) are assumed by many to regulate the abundance of invasive mesopredators, such as red foxes Vulpes vulpes and feral cats Felis catus, thereby providing indirect benefits to various threatened vertebrates. Several recent papers have claimed to provide evidence for the biodiversity benefits of dingoes in this way. Nevertheless, in this paper we highlight several critical weaknesses in the methodological approaches used in many of these reports, including lack of consideration for seasonal and habitat differences in activity, the complication of simple track-based indices by incorporating difficult-to-meet assumptions, and a reduction in sensitivity for assessing populations by using binary measures rather than potentially continuous measures. Of the 20 studies reviewed, 15 of them (75%) contained serious methodological flaws, which may partly explain the inconclusive nature of the literature nvestigating interactions between invasive Australian predators. We therefore assert that most of the “growing body of evidence” for mesopredator release is merely an inconclusive growing body of literature only. We encourage those interested in studying the ecological roles of dingoes relative to invasive mesopredators and native prey species to account for the factors we identify, and caution the value of studies that have not done so.

Increased virulence of rabbit haemorrhagic disease virus associated with genetic resistance in wild Australian rabbits (Oryctolagus cuniculus)

The release of myxoma virus (MYXV) and Rabbit Haemorrhagic Disease Virus (RHDV) in Australia with the aim of controlling overabundant rabbits has provided a unique opportunity to study the initial spread and establishment of emerging pathogens, as well as their co-evolution with their mammalian hosts. In contrast to MYXV, which attenuated shortly after its introduction, rapid attenuation of RHDV has not been observed. By studying the change in virulence of recent field isolates at a single field site we show, for the first time, that RHDV virulence has increased through time, likely because of selection to overcome developing genetic resistance in Australian wild rabbits. High virulence also appears to be favoured as rabbit carcasses, rather than diseased animals, are the likely source of mechanical insect transmission. These findings not only help elucidate the co-evolutionary interaction between rabbits and RHDV, but reveal some of the key factors shaping virulence evolution.