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 feeding sorghum ergot (Claviceps africana) to sows during mid-lactation on plasma prolactin and litter performance

Diets containing 3% sorghum ergot (16 mg alkaloids/kg, including 14 mg dihydroergosine/kg) were fed to 12 sows from 14 days post-farrowing until weaning 14 days later, and their performance was compared with that of 10 control sows. Ergot-fed sows displayed a smaller weight loss during lactation of 24 kg/head vs. 29 kg/head in control sows (p > 0.05) despite feed consumption being less (61 kg/head total feed intake vs. 73 kg/head by control sows; p < 0.05). Ergot-fed sows had poorer weight gain of litters over the 14-day period (16.6 kg/litter vs. 28.3 kg/litter for controls; p < 0.05) despite an increase in consumption of creep feed by the piglets from the ergot-fed sows (1.9 kg/litter compared with 1.1 kg/litter by the control; p > 0.05). Sow plasma prolactin was reduced with ergot feeding after 7 days to 4.8 μg/l compared with 15.1 μg/l in the control sows (p < 0.01) and then at weaning was 4.9 μg/l compared with 8.0 μg/l (p < 0.01) in the control sows. Two sows fed ergot ceased lactation early, and the above sow feed intakes, body weight losses with litter weight gains and creep consumption indirectly indicate an ergot effect on milk production.

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.

Genetic variation and responses in reproductive performance of sows in lines selected for growth rate under restricted feeding

The objective of this study was to examine genetic changes in reproduction traits in sows (total number born (TNB), number born alive (NBA), average piglet birth weight (ABW) and number of piglets weaned (NW), body weight prior to mating (MW), gestation length (GL) and daily food intake during lactation (DFI)) in lines of Large White pigs divergently selected over 4 years for high and low post-weaning growth rate on a restricted ration. Heritabilities and repeatabilities of the reproduction traits were also determined. The analyses were carried out on 913 litter records using average information-restricted maximum likelihood method applied to single trait animal models. Estimates of heritability for most traits were small, except for ABW (0·33) and MW (0·35). Estimates of repeatability were slightly higher than those of heritability for TNB, NBA and NW, but they were almost identical for ABW, MW, GL and DFI. After 4 years of selection, the high growth line sows had significantly heavier body weight prior to mating and produced significantly more piglets born alive with heavier average birth weight than the low line sows. There were, however, no statistical differences between the selected lines in TNB or NW. The lower food intake of high relative to low line sows during lactation was not significant, indicating that daily food intake differences found between grower pigs in the high and low lines (2·71 v. 2·76 kg/day, s.e.d. 0·024) on ad libitum feeding were not fully expressed in lactating sows. It is concluded that selection for growth rate on the restricted ration resulted in beneficial effects on important measures of reproductive performance of the sows.

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.

Multi-trait assessment of early-in-life female, male and genomic measures for use in genetic selection to improve female reproductive performance of Brahman cattle

Early-in-life female and male measures with potential to be practical genetic indicators were chosen from earlier analyses and examined together with genomic measures for multi-trait use to improve female reproduction of Brahman cattle. Combinations of measures were evaluated on the genetic gains expected from selection of sires and dams for each of age at puberty (AGECL, i.e. first observation of a corpus luteum), lactation anoestrous interval in 3-year-old cows (LAI), and lifetime annual weaning rate (LAWR, i.e. the weaning rate of cows based on the number of annual matings they experienced over six possible matings). Selection was on an index of comparable records for each combination. Selection intensities were less than theoretically possible but assumed a concerted selection effort was able to be made across the Brahman breed. The results suggested that substantial genetic gains could be possible but need to be confirmed in other data. The estimated increase in LAWR in 10 years, for combinations without or with genomic measures, ranged from 8 to 12 calves weaned per 100 cows from selection of sires, and from 12 to 15 calves weaned per 100 cows from selection of sires and dams. Corresponding reductions in LAI were 60-103 days or 94-136 days, and those for AGECL were 95-125 or 141-176 days, respectively. Coat score (a measure of the sleekness or wooliness of the coat) and hip height in females, and preputial eversion and liveweight in males, were measures that may warrant wider recording for Brahman female reproduction genetic evaluation. Pregnancy-test outcomes from Matings 1 and 2 also should be recorded. Percentage normal sperm may be important to record for reducing LAI and scrotal size and serum insulin-like growth factor-I concentration in heifers at 18 months for reducing AGECL. Use of a genomic estimated breeding value (EBV) in combination with other measures added to genetic gains, especially at genomic EBV accuracies of 40%. Accuracies of genomic EBVs needed to approach 60% for the genomic EBV to be the most important contributor to gains in the combinations of measures studied.

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.

Immediate and residual effects of heat stress and restricted intake on milk protein and casein composition and energy metabolism

The effects of heat stress on dairy production can be separated into 2 distinct causes: those effects that are mediated by the reduced voluntary feed intake associated with heat stress, and the direct physiological and metabolic effects of heat stress. To distinguish between these, and identify their effect on milk protein and casein concentration, mid-lactation Holstein-Friesian cows (n = 24) were housed in temperature-controlled chambers and either subjected to heat stress HS; temperature-humidity index (THI) ~78 or kept in a THI < 70 environment and pair-fed with heat-stressed cows (TN-R) for 7 d. A control group of cows was kept in a THI < 70 environment with ad libitum feeding (TN-AL). A subsequent recovery period (7 d), with THI < 70 and ad libitum feeding followed. Intake accounted for only part of the effects of heat stress. Heat stress reduced the milk protein concentration, casein number, and casein concentration and increased the urea concentration in milk beyond the effects of restriction of intake. Under HS, the proportion in total casein of αS1-casein increased and the proportion of αS2-casein decreased. Because no effect of HS on milk fat or lactose concentration was found, these effects appeared to be the result of specific downregulation of mammary protein synthesis, and not a general reduction in mammary activity. No residual effects were found of HS or TN-R on milk production or composition after THI < 70 and ad libitum intake were restored. Heat-stressed cows had elevated blood concentrations of urea and Ca, compared with TN-R and TN-AL. Cows in TN-R had higher serum nonesterified fatty acid concentrations than cows in HS. It was proposed that HS and TN-R cows may mobilize different tissues as endogenous sources of energy.