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.

Development of a DNA based aging technique for use in fisheries assessments. Final report, Australian Fisheries Research and Development Corporation Project 2007/033.

The productivity of a fisheries resource can be quantified from estimates of recruitment, individual growth and natural and fisheries-related mortality, assuming the spatial extent of the resource has been quantified and there is minimal immigration or emigration. The sustainability of a fisheries resource is facilitated by management controls such as minimum and maximum size limits and total allowable catch. Minimum size limits are often set to allow individuals the opportunity to reproduce at least once before the chance of capture. Total allowable catches are a proportion of the population biomass, which is estimated based on known reproduction, recruitment, mortality and growth rates. In some fisheries, however, management actions are put in place without quantification of the resource through the stock assessment process. This occurs because species-specific information, for example individual growth, may not be available. In these circumstances, management actions need to be precautionary to protect against future resource collapse, but this often means that the resource is lightly exploited. Consequently, the productivity of the resource is not fully realised. Australia’s most valuable fisheries are invertebrate fisheries (Australian Department of Agriculture Fisheries and Forestry, 2008). For example, Australian fisheries (i.e. excluding aquaculture) production of crustaceans (largely prawns, rock lobster and crab) was 41,000 tonnes in 2006/7, worth $778 million. Production from mollusc (largely abalone, scallops, oysters and squid) fisheries was 39,000 tonnes, worth $502 million. Together, in 2006/7 crustacean and mollusc fisheries represented 58% of the total value of Australian wild fisheries production. Sustainable management of Australia’s invertebrate fisheries is frustrated by the lack of data on species-specific growth rates. This project investigated a new method to estimate age, and hence individual growth rates, in invertebrate fisheries species. The principle behind the new aging method was that telomeres (i.e. DNA end-caps of chromosomes) get shorter as an individual gets older. We studied commercial crustacean and molluscan species. A vertebrate fish species (silver perch, Bidyanus bidyanus) was used as a control to standardise our work against the literature. We found a clear relationship between telomere length and shell size for temperate abalone (Haliotis rubra). Further research is recommended before the method can be implemented to assist management of wildharvested abalone populations. Age needs to be substituted for shell size in the relationship and it needs to be studied for abalone from several regions. This project showed that telomere length declined with increasing age in Sydney rock oysters (Saccostrea glomerata) and was affected by regional variation. A relationship was not apparent between telomere length and age (or size as a surrogate for age) for crustacean species (school prawns, Metapenaeus macleayi; eastern rock lobster, Sagmariasus verreauxi; southern rock lobster, Jasus edwardsii; and spanner crabs, Ranina ranina). For school prawns, there was no difference between telomere length in males and females. Further research is recommended, however, as telomeric DNA from crustaceans was difficult to analyse using the terminal restriction fragment (TRF) assay. Telomere lengths of spanner crabs and lobsters were at the upper limit of resolution of the assay used and results were affected by degradation and possible contamination of telomeric DNA. It is possible that telomere length is an indicator of remaining lifespan in molluscan and crustacean individuals, as suggested for some vertebrate species (e.g. Monaghan, 2010). Among abalone of similar shell size and among lobster pueruli, there was evidence of individuals having significantly longer or shorter telomeres than the group average. At a population level, this may be a surrogate for estimates of future natural mortality, which may have usefulness in the management of those populations. The method used to assay telomere length (terminal restriction fragment assay) performed adequately for most species, but it was too expensive and time-consuming to be considered a useful tool for gathering information for fisheries management. Research on alternative methods is strongly recommended.