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.

Effect of changes in hook pattern and size on catch rate, hooking location, injury and bleeding for a number of tropical reef fish species

The Queensland Great Barrier Reef line fishery in Australia is regulated via a range of input and output controls including minimum size limits, daily catch limits and commercial catch quotas. As a result of these measures a substantial proportion of the catch is released or discarded. The fate of these released fish is uncertain, but hook-related mortality can potentially be decreased by using hooks that reduce the rates of injury, bleeding and deep hooking. There is also the potential to reduce the capture of non-target species though gear selectivity. A total of 1053 individual fish representing five target species and three non-target species were caught using six hook types including three hook patterns (non-offset circle, J and offset circle), each in two sizes (small 4/0 or 5/0 and large 8/0). Catch rates for each of the hook patterns and sizes varied between species with no consistent results for target or non-target species. When data for all of the fish species were aggregated there was a trend for larger hooks, J hooks and offset circle hooks to cause a greater number of injuries. Using larger hooks was more likely to result in bleeding, although this trend was not statistically significant. Larger hooks were also more likely to foul-hook fish or hook fish in the eye. There was a reduction in the rates of injuries and bleeding for both target and non-target species when using the smaller hook sizes. For a number of species included in our study the incidence of deep hooking decreased when using non-offset circle hooks, however, these results were not consistent for all species. Our results highlight the variability in hook performance across a range of tropical demersal finfish species. The most obvious conservation benefits for both target and non-target species arise from using smaller sized hooks and non-offset circle hooks. Fishers should be encouraged to use these hook configurations to reduce the potential for post-release mortality of released fish.

Space allowances for confined livestock and their determination from allometric principles

The amount of space provided to animals governs important elements of their behaviour and, hence, is critical for their health and welfare. We review the use of allometric principles and equations to estimate the static space requirements of animals when standing and lying, and the space required for animals to feed, drink, stand-up and lie-down. We use the research literature relating to transportation and intensive housing of sheep and cattle to assess the validity of allometric equations for estimating space allowances. We investigated these areas because transportation and intensive housing provide points along a continuum in terms of the duration of confinement, (from hours to months) and spatial requirements are likely to increase with increasing duration of confinement, as animals will need to perform a greater behavioural repertoire for long-term survival, health and welfare. We find that, although there are theoretical reasons why allometric relationships to space allowances may vary slightly for different classes of stock, space allowances that have been demonstrated to have adverse effects on animal welfare during transportation correlated well with an inability to accommodate standing animals, as estimated from allometry. For intensive housing, we were able to detect a space allowance below which there were adverse effects on welfare. For short duration transportation during which animals remain standing, a space allowance per animal described by the allometric equation: area (m^2) = 0.020W^0.66, where W = liveweight (kg), would appear to be appropriate. Where it is desirable for all animals to lie simultaneously, then a minimum space allowance per animal described by the allometric equation: area (m^2) = 0.027W^0.66 appears to permit this, given that animals in a group time-share space. However, there are insufficient data to determine whether this allowance onboard a vehicle/vessel would enable animals to move and access food and water with ease. In intensive housing systems, a minimum space allowance per animal described by the allometric equation: area (m^2) = 0.033W^0.66 appears to be the threshold below which there are adverse effects on welfare. These suggested space allowances require verification with a range of species under different thermal conditions and, for transportation, under different conditions of vehicular/vessel stability. The minimum length of trough per animal (L in m) required for feeding and drinking can be determined from L = 0.064W^0.33, with the number of animals required to feed/drink simultaneously taken into account, together with any requirement to minimise competition. This also requires verification with a range of species. We conclude that allometric relationships are an appropriate basis for the formulation of space allowances for livestock.

Refined Global Analysis of Bemisia tabaci (Hemiptera: Sternorrhyncha: Aleyrodoidea: Aleyrodidae) Mitochondrial Cytochrome Oxidase 1 to Identify Species Level Genetic Boundaries.

Identifying species boundaries within morphologically indistinguishable cryptic species complexes is often contentious. For the whitefly Bemisia tabaci (Gennadius) (Hemiptera: Sternorrhyncha: Aleyrodoidea: Aleyrodidae), the lack of a clear understanding about the genetic limits of the numerous genetic groups and biotypes so far identified has resulted in a lack of consistency in the application of the terms, the approaches use to apply them and in our understanding of what genetic structure within B. tabaci means. Our response has been to use mitochondrial gene cytochrome oxidase one to consider how to clearly and consistently define genetic separation. Using Bayesian phylogenetic analysis and analysis of sequence pairwise divergence we found a considerably higher to number of genetic groups than had been previously determined with two breaks in the distribution, one at 11% and another at 3.5%. At >11% divergence, 11 distinct groups were resolved, whereas at >3.5% divergence 24 groups were identified. Consensus sequences for each of these groups were determined and were shown to be useful in the correct assignment of sequences of unknown origin. The 3.5% divergence bound is consistent with species level separations in other insect taxa and Suggests that B. tabaci is it cryptic species composed of at least 24 distinct species. We further show that the placement of Bemesia atriplex (Froggatt) within the B. tabaci in, group adds further weight to the argument for species level separation within B. tabaci. This new analysis, which constructs consensus sequences and uses these its a standard against which unknown sequences call be compared, provides for the first time it consistent means of identifying the genetic hounds of each species with it high degree of certainty.

QTL for root angle and number in a population developed from bread wheats (Triticum aestivum) with contrasting adaptation to water-limited environments

Root architecture traits in wheat are important in deep soil moisture acquisition and may be used to improve adaptation to water-limited environments. The genetic architecture of two root traits, seminal root angle and seminal root number, were investigated using a doubled haploid population derived from SeriM82 and Hartog. Multiple novel quantitative trait loci (QTL) were identified, each one having a modest effect. For seminal root angle, four QTL (-log10(P) >3) were identified on 2A, 3D, 6A and 6B, and two suggestive QTL (-log10(P) >2) on 5D and 6B. For root number, two QTL were identified on 4A and 6A with four suggestive QTL on 1B, 3A, 3B and 4A. QTL for root angle and root number did not co-locate. Transgressive segregation was found for both traits. Known major height and phenology loci appear to have little effect on root angle and number. Presence or absence of the T1BL.1RS translocation did not significantly influence root angle. Broad sense heritability (h 2) was estimated as 50 % for root angle and 31 % for root number. Root angle QTL were found to be segregating between wheat cultivars adapted to the target production region indicating potential to select for root angle in breeding programs. © 2013 Springer-Verlag Berlin Heidelberg.

Circulating concentrations of leptin, ovarian follicle number, and oocyte lipid content and active mitochondria, in Zebu crossbred cows maintained on standard or improved nutrition

Zebu (Bos indicus) crossbred beef cows (Droughtmaster) were maintained long-term (16 months) on standard nutrition (SN) or improved nutrition (IN). Cows on IN had better body condition and greater (P<0.05) circulating concentrations of leptin than cows on SN (0.7±0.1n/ml and 1.7±0.1n/ml, respectively). There were no outstanding differences between SN and IN cows in basal number of ovarian follicles (≤4mm, 5-8mm, and≥9mm) and there were also no differences in number of oocytes recovered by oocyte pick-up. Cows on IN had a greater (P<0.05) number of total follicles after stimulation with FSH than cows on SN. Oocytes from cows on IN had greater (P<0.05) lipid content than cows on SN (-0.23±0.16 and 0.20±0.18 arbitrary units, respectively) and oocytes of the former cows also tended to have more active mitochondria, although this was not significant. Cows on IN showed a positive relationship (R2=0.31, P<0.05) between plasma leptin and oocyte lipid content. Lipids are utilized by oocytes during high energy consumptive processes including fertilization and early cleavage. The greater lipid content of oocytes from IN cows could therefore confer a reproductive advantage. The present study has shown relationships between nutrition, body condition, circulating leptin, and oocyte lipid content, but a clear cause-and-effect requires further investigation in the cow. © 2013 Elsevier B.V.

Variable-number tandem repeats genotyping used to aid and inform management strategies for a bovine Johne's disease incursion in tropical and subtropical Australia

The application of variable-number tandem repeats (VNTR) genotyping of Mycobacterium avium subsp. paratuberculosis isolates to assist in investigating incidents of bovine Johne’s disease in a low-prevalence region of Australia is described in the current study. Isolates from a response to detection of bovine Johne’s disease in Queensland were compared with strains from national and international sources. The tandem application of mycobacterial interspersed repetitive unit (MIRU) and multilocus short sequence repeats (MLSSR) genotyping identified 2 strains, 1 that infected cattle on multiple properties with trace-forward histories from a common infected property, and 1 genotypically different strain recovered from a single property. The former strain showed an identical genotype to an isolate from India. Neither strain showed a genotypic link to regions of Australia with a higher prevalence of the disease. Genotyping has indicated incursions from 2 independent sources. This intelligence has informed investigations into potential routes of entry and the soundness of ongoing control measures, and supported strategy and policy decisions regarding management of Mycobacterium avium subsp. paratuberculosis incursions for Queensland.

A universal equation to predict methane production of forage-fed cattle in Australia

The methods for estimating methane emissions from cattle as used in the Australian national inventory are based on older data that have now been superseded by a large amount of more recent data. Recent data suggested that the current inventory emissions estimates can be improved. To address this issue, a total of 1034 individual animal records of daily methane production (MP) was used to reassess the relationship between MP and each of dry matter intake (DMI) and gross energy intake (GEI). Data were restricted to trials conducted in the past 10 years using open-circuit respiration chambers, with cattle fed forage-based diets (forage >70%). Results from diets considered to inhibit methanogenesis were omitted from the dataset. Records were obtained from dairy cattle fed temperate forages (220 records), beef cattle fed temperate forages (680 records) and beef cattle fed tropical forages (133 records). Relationships were very similar for all three production categories and single relationships for MP on a DMI or GEI basis were proposed for national inventory purposes. These relationships were MP (g/day) = 20.7 (±0.28) × DMI (kg/day) (R2 = 0.92, P < 0.001) and MP (MJ/day) = 0.063 (±0.008) × GEI (MJ/day) (R2 = 0.93, P < 0.001). If the revised MP (g/day) approach is used to calculate Australia’s national inventory, it will reduce estimates of emissions of forage-fed cattle by 24%. Assuming a global warming potential of 25 for methane, this represents a 12.6 Mt CO2-e reduction in calculated annual emissions from Australian cattle.