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.

Extent of dense native woodland and exotic weed infestation in the extensive grazing lands of the Upper Herbert and Upper Burdekin River Catchments of far north Queensland: results of a producer survey

The Mt Garnet Landcare Group commissioned a survey of landholders within the Upper Herbert and Upper Burdekin River Catchments to assess the density of native woodlands and to gauge the extent of exotic weed infestation. Twenty-four of 49 landholders responded, representing an area of nearly 500 000 ha or 47% of the total area. Dense native woodland covers 24% (>117 000 ha) of the area surveyed, while a further 30% (140 000 ha) supports moderately dense stands. The dense stands are largely confined to the highly fertile alluvial soils (26% dense woodland) and the lower fertility sandy-surfaced soils (33% or >96 000 ha). Moderate and dense infestations of exotic weeds, principally Lantana camara, occur on 54% (20 000 ha) of alluvial soils and on 13% of sandy-surfaced soils (39 000 ha), where praxelis (Praxelis clematidia) is the major weed. Basaltic soils have low levels of both dense woodland and exotic weed infestation. Some implications of the results are discussed.

New records of the River Shark Glyphis (Carcharhinidae) reported from Cape York Peninsula, northern Australia

The distribution of the river shark Glyphis in northern Australia is extended with new records of occurrence in the Gulf of Carpentaria and a reassessment of historical survey data from Cape York Peninsula. Nine new specimens of Glyphis sp. A were collected in 2005 from the Weipa region on the Queensland coast of the Gulf of Carpentaria. A re-examination of archival records from 1978-86 marine and estuarine fish surveys in the Gulf of Carpentaria and along the northern Queensland East Coast allowed a further nineteen Glyphis specimens to be identified. Combined this gives twenty-eight new records of Glyphis specimens from the coasts of Cape York Peninsula, Queensland. Common habitat characteristics for all captures were turbid, shallow, fast running tidal water in the upper reaches of coastal rivers. The substrate was generally muddy and the rivers lined with mangrove. In all surveys the representation of Glyphis was low, being less than 1% of the total shark captures historically and 0.002 sharks 50 m net hour-1 in Weipa 2005. The size range captured was 1000-1800 mm total length historically and 705-1200 mm total length from Weipa 2005, with none recorded as sexually mature. Diagnostic characteristics of the Weipa specimens, identified as Glyphis sp. A, were: lower jaw teeth protruding and "spear-like"; second dorsal fin greater than half the height of the first dorsal fin; the snout relatively short and fleshy in the lateral view; pectoral fin ventral surface black in colouration; the precaudal vertebral count between 118 and 123; and the total vertebral count between 204 and 209.

In vivo delivery of bovine viral diahorrea virus, E2 protein using hollow mesoporous silica nanoparticles

Our work focuses on the application of mesoporous silica nanoparticles as a combined delivery vehicle and adjuvant for vaccine applications. Here we present results using the viral protein, E2, from bovine viral diarrhoea virus (BVDV). BVDV infection occurs in the target species of cattle and sheep herds worldwide and is therefore of economic importance. E2 is a major immunogenic determinant of BVDV and is an ideal candidate for the development of a subunit based nanovaccine using mesoporous silica nanoparticles. Hollow type mesoporous silica nanoparticles with surface amino functionalisation (termed HMSA) were characterised and assessed for adsorption and desorption of E2. A codon-optimised version of the E2 protein (termed Opti-E2) was produced in Escherichia coli. HMSA (120 nm) had an adsorption capacity of 80 [small mu ]g Opti-E2 per mg HMSA and once bound E2 did not dissociate from the HMSA. Immunisation studies in mice with a 20 [small mu ]g dose of E2 adsorbed to 250 [small mu ]g HMSA was compared to immunisation with Opti-E2 (50 [small mu ]g) together with the traditional adjuvant Quillaja saponaria Molina tree saponins (QuilA, 10 [small mu ]g). The humoral responses with the Opti-E2/HMSA nanovaccine although slightly lower than those obtained for the Opti-E2 + QuilA group demonstrated that HMSA particles are an effective adjuvant that stimulated E2-specific antibody responses. Importantly the cell-mediated immune responses were consistently high in all mice immunised with Opti-E2/HMSA nanovaccine formulation. Therefore we have shown the Opti-E2/HMSA nanoformulation acts as an excellent adjuvant that gives both T-helper 1 and T-helper 2 mediated responses in a small animal model. This study has provided proof-of-concept towards the development of an E2 subunit nanoparticle based vaccine.

A novel method for the age estimation of Saddletail snapper (Lutjanus malabaricus) using Fourier Transform-near infrared (FT-NIR) spectroscopy

Near infrared (NIR) spectroscopy was investigated as a potential rapid method of estimating fish age from whole otoliths of Saddletail snapper (Lutjanus malabaricus). Whole otoliths from 209 Saddletail snapper were extracted and the NIR spectral characteristics were acquired over a spectral range of 800–2780 nm. Partial least-squares models (PLS) were developed from the diffuse reflectance spectra and reference-validated age estimates (based on traditional sectioned otolith increments) to predict age for independent otolith samples. Predictive models developed for a specific season and geographical location performed poorly against a different season and geographical location. However, overall PLS regression statistics for predicting a combined population incorporating both geographic location and season variables were: coefficient of determination (R2) = 0.94, root mean square error of prediction (RMSEP) = 1.54 for age estimation, indicating that Saddletail age could be predicted within 1.5 increment counts. This level of accuracy suggests the method warrants further development for Saddletail snapper and may have potential for other fish species. A rapid method of fish age estimation could have the potential to reduce greatly both costs of time and materials in the assessment and management of commercial fisheries.

Continuous combined microwave and hot air treatment of apples for fruit fly (Bactrocera tryoni and B. jarvisi) disinfestation

Apples at 24 ± 2 °C were heated in a pilot scale hot air assisted (40 °C) continuous pentagonal microwave system, to evaluate the effectiveness of this treatment on insect mortality (variety Mutsu) and fruit quality (variety Granny Smith). An average temperature of 53.4 ± 1.3 °C at core, bottom and flesh of the apple was recorded at the end of the treatment. One hundred percent mortality of the most tolerant stage of Queensland fruit fly (Bactrocera tryoni, Froggatt) and Jarvis's fruit fly (Bactrocera jarvisi, Tryon), were observed when the Mortality value (M52, equivalent time of isothermal treatment at 52 °C) at the slowest heating point applicable for each experiment was ≥ 50 min and ≥ 37 min, respectively. Results showed that microwave heat treatment is effective for insect disinfestation without any adverse impact on total soluble solids, flesh or peel firmness of the treated apples. The treated apples recorded a significantly higher pH and lower ion leakage than the untreated apples after 3 or 4 weeks. Therefore, the microwave heat treatment has the potential to be developed as an alternative chemical free quarantine treatment against economically significant insect pests. Industrial relevance Hot air assisted microwave heating of fruits and vegetables, is more cost effective compared to vapour heat treatment and ionising radiation for disinfestation of insects. Microwave treatment is environmentally friendly compared to fumigation and chemical treatments. Hot air assisted microwave disinfestation can be performed at farms or centralised pack houses since the capital cost would be comparatively lower than vapour heat or ionising radiation treatments.

Maximum likelihood estimation of natural mortality and quantification of temperature effects on catchability of brown tiger prawn (Penaeus esculentus) in Moreton Bay (Australia) using logbook data

It is common to model the dynamics of fisheries using natural and fishing mortality rates estimated independently using two separate analyses. Fishing mortality is routinely estimated from widely available logbook data, whereas natural mortality estimations have often required more specific, less frequently available, data. However, in the case of the fishery for brown tiger prawn (Penaeus esculentus) in Moreton Bay, both fishing and natural mortality rates have been estimated from logbook data. The present work extended the fishing mortality model to incorporate an eco-physiological response of tiger prawn to temperature, and allowed recruitment timing to vary from year to year. These ecological characteristics of the dynamics of this fishery were ignored in the separate model that estimated natural mortality. Therefore, we propose to estimate both natural and fishing mortality rates within a single model using a consistent set of hypotheses. This approach was applied to Moreton Bay brown tiger prawn data collected between 1990 and 2010. Natural mortality was estimated by maximum likelihood to be equal to 0.032 ± 0.002 week−1, approximately 30% lower than the fixed value used in previous models of this fishery (0.045 week−1).