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.

Spatial model of site index based on Y-ray spectrometry and a digital elevation model for two Pinus species in Tuan Toolara State Forest, Queensland, Australia

Site index prediction models are an important aid for forest management and planning activities. This paper introduces a multiple regression model for spatially mapping and comparing site indices for two Pinus species (Pinus elliottii Engelm. and Queensland hybrid, a P. elliottii x Pinus caribaea Morelet hybrid) based on independent variables derived from two major sources: g-ray spectrometry (potassium (K), thorium (Th), and uranium (U)) and a digital elevation model (elevation, slope, curvature, hillshade, flow accumulation, and distance to streams). In addition, interpolated rainfall was tested. Species were coded as a dichotomous dummy variable; interaction effects between species and the g-ray spectrometric and geomorphologic variables were considered. The model explained up to 60% of the variance of site index and the standard error of estimate was 1.9 m. Uranium, elevation, distance to streams, thorium, and flow accumulation significantly correlate to the spatial variation of the site index of both species, and hillshade, curvature, elevation and slope accounted for the extra variability of one species over the other. The predicted site indices varied between 20.0 and 27.3 m for P. elliottii, and between 23.1 and 33.1 m for Queensland hybrid; the advantage of Queensland hybrid over P. elliottii ranged from 1.8 to 6.8 m, with the mean at 4.0 m. This compartment-based prediction and comparison study provides not only an overview of forest productivity of the whole plantation area studied but also a management tool at compartment scale.

A digital imaging technique to quantify the incidence of internal browning in potato (Solanum tuberosum) tubers

Internal browning disorders, including brown fleck (BF), in potato (Solanum tuberosum) tubers greatly reduce tuber quality, but the causes are not well understood. This is due, in part, to the highly variable data provided by visual value-based rating systems. A digital imaging technique was developed to quantify accurately the incidence of internal browning in potato tubers. Images of tuber sections were scanned using a flatbed scanner and digitally enhanced to highlight tuber BF lesions, and the area of affected tissue calculated using pixel quantification software. Digital imaging allowed for the determination of previously unused indices of the incidence and severity of internal browning in potato tubers. Statistical analysis of the comparison between digitally derived and visual-rating BF data from a glasshouse experiment showed that digital data greatly improved the delineation of treatment effects. The F-test probability was further improved through square root or logarithmic data transformations of the digital data, but not of the visual-rating data. Data from a field experiment showed that the area of tuber affected by BF and the number of small BF lesions increased with time and with increase in tuber size. The results from this study indicate that digital imaging of internal browning disorders of potato tubers holds much promise in determining their causes that heretofore have proved elusive.

Detecting the attributes of a wheat crop using digital imagery acquired from a low-altitude platform

A low-altitude platform utilising a 1.8-m diameter tethered helium balloon was used to position a multispectral sensor, consisting of two digital cameras, above a fertiliser trial plot where wheat (Triticum spp.) was being grown. Located in Cecil Plains, Queensland, Australia, the plot was a long-term fertiliser trial being conducted by a fertiliser company to monitor the response of crops to various levels of nutrition. The different levels of nutrition were achieved by varying nitrogen application rates between 0 and 120 units of N at 40 unit increments. Each plot had received the same application rate for 10 years. Colour and near-infrared images were acquired that captured the whole 2 ha plot. These images were examined and relationships sought between the captured digital information and the crop parameters imaged at anthesis and the at-harvest quality and quantity parameters. The statistical analysis techniques used were correlation analysis, discriminant analysis and partial least squares regression. A high correlation was found between the image and yield (R2 = 0.91) and a moderate correlation between the image and grain protein content (R2 = 0.66). The utility of the system could be extended by choosing a more mobile platform. This would increase the potential for the system to be used to diagnose the causes of the variability and allow remediation, and/or to segregate the crop at harvest to meet certain quality parameters.

Issues of Digital Ecosystems Approach to Biological Ecosystem Simulation Model Fit Optimization

This study examines the application of digital ecosystems concepts to a biological ecosystem simulation problem. The problem involves the use of a digital ecosystem agent to optimize the accuracy of a second digital ecosystem agent, the biological ecosystem simulation. The study also incorporates social ecosystems, with a technological solution design subsystem communicating with a science subsystem and simulation software developer subsystem to determine key characteristics of the biological ecosystem simulation. The findings show similarities between the issues involved in digital ecosystem collaboration and those occurring when digital ecosystems interact with biological ecosystems. The results also suggest that even precise semantic descriptions and comprehensive ontologies may be insufficient to describe agents in enough detail for use within digital ecosystems, and a number of solutions to this problem are proposed.

More fertile florets and grains per spike can be achieved at higher temperature in wheat lines with high spike biomass and sugar content at booting

An understanding of processes regulating wheat floret and grain number at higher temperatures is required to better exploit genetic variation. In this study we tested the hypothesis that at higher temperatures, a reduction in floret fertility is associated with a decrease in soluble sugars and this response is exacerbated in genotypes low in water soluble carbohydrates (WSC). Four recombinant inbred lines contrasting for stem WSC were grown at 20/10 degrees C and 11 h photoperiod until terminal spikelet, and then continued in a factorial combination of 20/10 degrees C or 28/14 degrees C with 11 h or 16 h photoperiod until anthesis. Across environments, High WSC lines had more grains per spike associated with more florets per spike. The number of fertile florets was associated with spike biomass at booting and, by extension, with glucose amount, both higher in High WSC lines. At booting, High WSC lines had higher fixed C-13 and higher levels of expression of genes involved in photosynthesis and sucrose transport and lower in sucrose degradation compared with Low WSC lines. At higher temperature, the intrinsic rate of floret development rate before booting was slower in High WSC lines. Grain set declined with the intrinsic rate of floret development before booting, with an advantage for High WSC lines at 28/14 degrees C and 16 h. Genotypic and environmental action on floret fertility and grain set was summarised in a model.

Comparison of the pregnancy rates and costs per calf born after fixed-time artificial insemination or artificial insemination after estrus detection in Bos indicus heifers

This study compared pregnancy rates (PRs) and costs per calf born after fixed-time artificial insemination (FTAI) or AI after estrus detection (i.e., estrus detection and AI, EDAI), before and after a single PGF2α treatment in Bos indicus (Brahman-cross) heifers. On Day 0, the body weight, body condition score, and presence of a CL (46% of heifers) were determined. The heifers were then alternately allocated to one of two FTAI groups (FTAI-1, n = 139) and (FTAI-2, n = 141) and an EDAI group (n = 273). Heifers in the FTAI groups received an intravaginal progesterone-releasing device (IPRD; 0.78 g of progesterone) and 1 mg of estradiol benzoate intramuscularly (im) on Day 0. Eight days later, the IPRD was removed and heifers received 500 μg of PGF2α and 300 IU of eCG im; 24 hours later, they received 1 mg estradiol benzoate im and were submitted to FTAI 30 to 34 hours later (54 and 58 hours after IPRD removal). Heifers in the FTAI-2 group started treatment 8 days after those in the FTAI-1 group. Heifers in the EDAI group were inseminated approximately 12 hours after the detection of estrus between Days 4 and 9 at which time the heifers that had not been detected in estrus received 500 μg of PGF2α im and EDAI continued until Day 13. Heifers in the FTAI groups had a higher overall PR (proportion pregnant as per the entire group) than the EDAI group (34.6% vs. 23.2%; P = 0.003), however, conception rate (PR of heifers submitted for AI) tended to favor the estrus detection group (34.6% vs. 44.1%; P = 0.059). The cost per AI calf born was estimated to be $267.67 and $291.37 for the FTAI and EDAI groups, respectively. It was concluded that in Brahman heifers typical of those annually mated in northern Australia FTAI compared with EDAI increases the number of heifers pregnant and reduces the cost per calf born.