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 stable isotope index to assess decadal-scale vegetation change and application to woodlands of the Burdekin catchment, Australia

Forty-four study sites were established in remnant woodland in the Burdekin River catchment in tropical north-east Queensland, Australia, to assess recent (decadal) vegetation change. The aim of this study was further to evaluate whether wide-scale vegetation 'thickening' (proliferation of woody plants in formerly more open woodlands) had occurred during the last century, coinciding with significant changes in land management. Soil samples from several depth intervals were size separated into different soil organic carbon (SOC) fractions, which differed from one another by chemical composition and turnover times. Tropical (C4) grasses dominate in the Burdekin catchment, and thus δ13C analyses of SOC fractions with different turnover times can be used to assess whether the relative proportion of trees (C3) and grasses (C4) had changed over time. However, a method was required to permit standardized assessment of the δ13C data for the individual sites within the 13 Mha catchment, which varied in soil and vegetation characteristics. Thus, an index was developed using data from three detailed study sites and global literature to standardize individual isotopic data from different soil depths and SOC fractions to reflect only the changed proportion of trees (C3) to grasses (C3) over decadal timescales. When applied to the 44 individual sites distributed throughout the Burdekin catchment, 64% of the sites were shown to have experienced decadal vegetation thickening, while 29% had remained stable and the remaining 7% had thinned. Thus, the development of this index enabled regional scale assessment and comparison of decadal vegetation patterns without having to rely on prior knowledge of vegetation changes or aerial photography.

A new culture method for lesser mealworm, Alphitobius diaperinus

A new culture method for lesser mealworm, Alphitobius diaperinus (Panzer), was developed to provide large numbers of adult lesser mealworms of approximately the same age for insecticide resistance testing. Culturing entailed allowing 100 adults to reproduce for 4 days in a wheat-based culture medium contained inside a plastic culture box, removing the adults from the medium, and then rearing their progeny to adulthood therein, in approximately 56 days at 32 degrees C and 55% RH. During their development, progeny were supplied water via apple slices at 0, 21 and 35 days, and a foam substrate in which to pupate, also at 35 days. During 2004-2005, adult lesser mealworms were collected from six broiler-house populations and then cultured with this method. Each population produced 4500 adults required to complete resistance testing with one insecticide within ten culture boxes, at an average of 798 adults per culture box.

Vermiculite improves early development and survival of low chill stone-fruit embryos rescued in vitro

There are many reports of efficient embryo germination and the method has been optimized to suit subtropical low chill genotypes. However the subsequent growth, vigor, and ability of germinated embryos to develop and survive acclimatization is rarely reported. Many germinated embryos do not survive acclimatization, develop slowly, or fail to develop normally. Methods to improve plant development from in vitro embryo cultures are needed to improve the number of plants that survive to be useful in breeding programs. This paper describes an improved method of embryo rescue that significantly increases embryo shoot and root development that leads to increased plant survival. Four treatments: Woody Plant Media (WPM) solidified with agar, vermiculite with liquid WPM, vermiculite with WPM plus agar, and conventional stratification, were evaluated for embryo growth and subsequent plantlet development and survival for two low-chill peach and one low-chill nectarine cultivar. Highly significant improvements were found for shoot and root development of seedlings germinated in vermiculite based media compared to embryos germinated in conventional agar-based media. Vermiculite with WPM and agar improved plantlet growth subsequent to in vitro culture and significantly increased survival of germinated embryos resulting in more plants reaching the field.

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.

Development of commercial application of an avocado fruit robustness test

Fruit quality is one of the major factors limiting growth in avocado retail sales. Avocado growers are often unaware of their end-use fruit quality since quality problems only manifest upon fruit ripening and growers receive limited feedback from the supply chain. If growers were aware of their expected fruit quality they would be equipped to make better marketing decisions and if necessary to take remedial actions to improve their fruit quality. Avotest is being developed as a quick and easy method of determining expected end-use fruit quality before the start of the commercial fruit harvest. The test aims at distinguishing between blocks with robust fruit and those with less robust fruit. The test could also be used to predict the resulting fruit quality after the implementation of new farming practices.

Objective `Safe' Grazing Capacities for South-West Queensland Australia: Development of a Model for Individual Properties

Few tools are available to assist graziers, land administrators and financiers in making objective grazing capacity decisions on Australian rangelands, despite existing knowledge regarding stocking rate theory and the impact of stocking rates on land condition. To address this issue a model for objectively estimating 'safe' grazing capacities on individual grazing properties in south-west Queensland was developed. The method is based on 'safe' levels of utilisation (15%-20%) by domestic livestock of average annual forage grown for each land system on a property. Average annual forage grown (kglha) was calculated as the product of the rainfall use efficiency (kglhdmm) and average annual rainfall (mm) for a land system. This estimate included the impact of tree and shrub cover on forage production. The 'safe' levels of forage utilisation for south- west Queensland pastures were derived from the combined experience of (1) re-analysis of the results of grazing trials, (2) reaching a consensus on local knowledge and (3) examination of existing grazing practice on 'benchmark' grazing properties. We recognise the problems in defining, determining and using grazing capacity values, but consider that the model offers decision makers a tool that can be used to assess the grazing capacity of individual properties.

Bovine cysticercosis-Development of a real-time PCR to enhance classification of suspect cysts identified at meat inspection

Laboratory confirmation methods are important in bovine cysticerosis diagnosis as other pathologies can result in morphologically similar lesions resulting in false identifications. We developed a probe-based real-time PCR assay to identify Taenia saginata in suspect cysts encountered at meat inspection and compared its use with the traditional method of identification, histology, as well as a published nested PCR. The assay simultaneously detects T. saginata DNA and a bovine internal control using the cytochrome c oxidase subunit 1 gene of each species and shows specificity against parasites causing lesions morphologically similar to those of T. saginata. The assay was sufficiently sensitive to detect 1 fg (Ct 35.09 +/- 0.95) of target DNA using serially-diluted plasmid DNA in reactions spiked with bovine DNA as well as in all viable and caseated positive control cysts. A loss in PCR sensitivity was observed with increasing cyst degeneration as seen in other molecular methods. In comparison to histology, the assay offered greater sensitivity and accuracy with 10/19 (53%) T. saginata positives detected by real-time PCR and none by histology. When the results were compared with the reference PCR, the assay was less sensitive but offered advantages of faster turnaround times and reduced contamination risk. Estimates of the assay's repeatability and reproducibility showed the assay is highly reliable with reliability coefficients greater than 0.94. Crown Copyright (C) 2013 Published by Elsevier B.V. All rights reserved.

Methods to determine resistance to anthelmintics when continuing larval development occurs

The in vivo faecal egg count reduction test (FECRT) is the most commonly used test to detect anthelmintic resistance (AR) in gastrointestinal nematodes (GIN) of ruminants in pasture based systems. However, there are several variations on the method, some more appropriate than others in specific circumstances. While in some cases labour and time can be saved by just collecting post-drench faecal worm egg counts (FEC) of treatment groups with controls, or pre- and post-drench FEC of a treatment group with no controls, there are circumstances when pre- and post-drench FEC of an untreated control group as well as from the treatment groups are necessary. Computer simulation techniques were used to determine the most appropriate of several methods for calculating AR when there is continuing larval development during the testing period, as often occurs when anthelmintic treatments against genera of GIN with high biotic potential or high re-infection rates, such as Haemonchus contortus of sheep and Cooperia punctata of cattle, are less than 100% efficacious. Three field FECRT experimental designs were investigated: (I) post-drench FEC of treatment and controls groups, (II) pre- and post-drench FEC of a treatment group only and (III) pre- and post-drench FEC of treatment and control groups. To investigate the performance of methods of indicating AR for each of these designs, simulated animal FEC were generated from negative binominal distributions with subsequent sampling from the binomial distributions to account for drench effect, with varying parameters for worm burden, larval development and drench resistance. Calculations of percent reductions and confidence limits were based on those of the Standing Committee for Agriculture (SCA) guidelines. For the two field methods with pre-drench FEC, confidence limits were also determined from cumulative inverse Beta distributions of FEC, for eggs per gram (epg) and the number of eggs counted at detection levels of 50 and 25. Two rules for determining AR: (1) %reduction (%R) < 95% and lower confidence limit <90%; and (2) upper confidence limit <95%, were also assessed. For each combination of worm burden, larval development and drench resistance parameters, 1000 simulations were run to determine the number of times the theoretical percent reduction fell within the estimated confidence limits and the number of times resistance would have been declared. When continuing larval development occurs during the testing period of the FECRT, the simulations showed AR should be calculated from pre- and post-drench worm egg counts of an untreated control group as well as from the treatment group. If the widely used resistance rule 1 is used to assess resistance, rule 2 should also be applied, especially when %R is in the range 90 to 95% and resistance is suspected.

Development of a Multiplexed Bead-Based Suspension Array for the Detection and Discrimination of Pospiviroid Plant Pathogens

Efficient and reliable diagnostic tools for the routine indexing and certification of clean propagating material are essential for the management of pospiviroid diseases in horticultural crops. This study describes the development of a true multiplexed diagnostic method for the detection and identification of all nine currently recognized pospiviroid species in one assay using Luminex bead-based suspension array technology. In addition, a new data-driven, statistical method is presented for establishing thresholds for positivity for individual assays within multiplexed arrays. When applied to the multiplexed array data generated in this study, the new method was shown to have better control of false positives and false negative results than two other commonly used approaches for setting thresholds. The 11-plex Luminex MagPlex-TAG pospiviroid array described here has a unique hierarchical assay design, incorporating a near-universal assay in addition to nine species-specific assays, and a co-amplified plant internal control assay for quality assurance purposes. All assays of the multiplexed array were shown to be 100% specific, sensitive and reproducible. The multiplexed array described herein is robust, easy to use, displays unambiguous results and has strong potential for use in routine pospiviroid indexing to improve disease management strategies.

Flaws in the current method for calculating methane emissions during dairy manure management in New Zealand

New Zealand's Greenhouse Gas Inventory (the NZ Inventory) currently estimates methane (CH4) emissions from anaerobic dairy effluent ponds by: (1) determining the total pond volume across New Zealand; (2) dividing this volume by depth to obtain the total pond surface area; and (3) multiplying this area by an observational average CH4 flux. Unfortunately, a mathematically erroneous determination of pond volume has led to an imbalanced equation and a geometry error was made when scaling-up the observational CH4 flux. Furthermore, even if these errors are corrected, the nationwide estimate still hinges on field data from a study that used a debatable method to measure pond CH4 emissions at a single site, as well as a potentially inaccurate estimation of the amount of organic waste anaerobically treated. The development of a new methodology is therefore critically needed.

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.

Development of a Rapid Multiplex PCR Assay To Genotype Pasteurella multocida Strains by Use of the Lipopolysaccharide Outer Core Biosynthesis Locus

Pasteurella multocida is a Gram-negative bacterial pathogen that is the causative agent of a wide range of diseases in many animal species, including humans. A widely used method for differentiation of P. multocida strains involves the Heddleston serotyping scheme. This scheme was developed in the early 1970s and classifies P. multocida strains into 16 somatic or lipopolysaccharide (LPS) serovars using an agar gel diffusion precipitin test. However, this gel diffusion assay is problematic, with difficulties reported in accuracy, reproducibility, and the sourcing of quality serovar-specific antisera. Using our knowledge of the genetics of LPS biosynthesis in P. multocida, we have developed a multiplex PCR (mPCR) that is able to differentiate strains based on the genetic organization of the LPS outer core biosynthesis loci. The accuracy of the LPS-mPCR was compared with classical Heddleston serotyping using LPS compositional data as the "gold standard." The LPS-mPCR correctly typed 57 of 58 isolates; Heddleston serotyping was able to correctly and unambiguously type only 20 of the 58 isolates. We conclude that our LPS-mPCR is a highly accurate LPS genotyping method that should replace the Heddleston serotyping scheme for the classification of P. multocida strains.