Rapid genome wide mapping of phosphine resistance loci by a simple regional averaging analysis in the red flour beetle, Tribolium castaneum

Background Next-generation sequencing technology is an important tool for the rapid, genome-wide identification of genetic variations. However, it is difficult to resolve the ‘signal’ of variations of interest and the ‘noise’ of stochastic sequencing and bioinformatic errors in the large datasets that are generated. We report a simple approach to identify regional linkage to a trait that requires only two pools of DNA to be sequenced from progeny of a defined genetic cross (i.e. bulk segregant analysis) at low coverage (<10×) and without parentage assignment of individual SNPs. The analysis relies on regional averaging of pooled SNP frequencies to rapidly scan polymorphisms across the genome for differential regional homozygosity, which is then displayed graphically. Results Progeny from defined genetic crosses of Tribolium castaneum (F4 and F19) segregating for the phosphine resistance trait were exposed to phosphine to select for the resistance trait while the remainders were left unexposed. Next generation sequencing was then carried out on the genomic DNA from each pool of selected and unselected insects from each generation. The reads were mapped against the annotated T. castaneum genome from NCBI (v3.0) and analysed for SNP variations. Since it is difficult to accurately call individual SNP frequencies when the depth of sequence coverage is low, variant frequencies were averaged across larger regions. Results from regional SNP frequency averaging identified two loci, tc_rph1 on chromosome 8 and tc_rph2 on chromosome 9, which together are responsible for high level resistance. Identification of the two loci was possible with only 5-7× average coverage of the genome per dataset. These loci were subsequently confirmed by direct SNP marker analysis and fine-scale mapping. Individually, homozygosity of tc_rph1 or tc_rph2 results in only weak resistance to phosphine (estimated at up to 1.5-2.5× and 3-5× respectively), whereas in combination they interact synergistically to provide a high-level resistance >200×. The tc_rph2 resistance allele resulted in a significant fitness cost relative to the wild type allele in unselected beetles over eighteen generations. Conclusion We have validated the technique of linkage mapping by low-coverage sequencing of progeny from a simple genetic cross. The approach relied on regional averaging of SNP frequencies and was used to successfully identify candidate gene loci for phosphine resistance in T. castaneum. This is a relatively simple and rapid approach to identifying genomic regions associated with traits in defined genetic crosses that does not require any specialised statistical analysis.

Salmonella and on-farm risk factors in healthy slaughter-age cattle and sheep in eastern Australia

To examine healthy slaughter-age cattle and sheep on-farm for the excretion of Salmonella serovars in faeces and to identify possible risk factors using a questionnaire. The study involved 215 herds and flocks in the four eastern states of Australia, 56 with prior history of salmonellosis. Production systems examined included pasture beef cattle, feedlot beef cattle, dairy cattle, prime lambs and mutton sheep and animals were all at slaughter age. From each herd or flock, 25 animals were sampled and the samples pooled for Salmonella culture. All Salmonella isolated were serotyped and any Salmonella Typhimurium isolates were phage typed. Questionnaires on each production system, prepared in Epi Info 6.04, were designed to identify risk factors associated with Salmonella spp excretion, with separate questionnaires designed for each production system. Salmonellae were identified in all production systems and were more commonly isolated from dairies and beef feedlots than other systems. Statistical analysis revealed that dairy cattle were significantly more likely to shed Salmonella in faeces than pasture beef cattle, mutton sheep and prime lambs (P < 0.05). A wide diversity of Salmonella serovars, all of which have been isolated from humans in Australia, was identified in both cattle and sheep. Analysis of the questionnaires showed access to new arrivals was a significant risk factor for Salmonella excretion on dairy properties. For beef feedlots, the presence of large numbers of flies in the feedlot pens or around stored manure were significant risk factors for Salmonella excretion. Dairy cattle pose the highest risk of all the slaughter-age animals tested. Some of the identified risk factors can be overcome by improved management practices, especially in relation to hygiene.

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.

Australian Bat Lyssavirus

In Chapter 1, the literature relating to rabies virus and the rabies like lyssaviruses is reviewed. In Chapter 2, data are presented from 1170 diagnostic submissions for ABLV testing by fluorescent antibody test (Centocor FAT). All 27 non-bat submissions were ABLV-negative. Of 1143 bat accessions 74 (16%) were ABLV-positive, including 69 of 974 (7.1%) flying foxes (Pteropus spp.), 5 of 7 (71.4%) Saccolaimus flaviventris (Yellow-bellied sheathtail bats), none of 151 other microchiropteran bats, and none of 11 unidentified bats. Statistical analysis of data from 868 wild Black, Grey-headed, Little Red and Spectacled flying foxes (Pteropus alecto, P. poliocephalus, P. scapulatus, and P. conspicillatus) indicated that three factors; species, health status and age were associated with significant (p< 0.001) differences in the proportion of ABLV-positive bats. Other factors including sex, whether the bat bit a person or animal, region, year, and season submitted, were not associated with ABLV. Case data for 74 ABLV-positive bats, including the circumstances in which they were found and clinical signs, is presented. In Chapter 3, the aetiological diagnosis was investigated for 100 consecutive flying fox submissions with neurological signs. ABLV (32%), spinal and head injuries (29%), and neuro-angiostrongylosis (18%) accounted for most neurological syndromes in flying foxes. No evidence of lead poisoning was found in unwell (n=16) or healthy flying foxes (n=50). No diagnosis was reached for 16 cases, all of which were negative for ABLV by TaqMan PCR. The molecular diversity of ABLV was examined in Chapter 4 by sequencing 36 bases of the leader sequence, the entire N gene, and start of the P gene of 28 isolates from pteropid bats and 3 isolates from Yellow-bellied sheathtail (YBST) bats. Phylogenetic analysis indicated all ABLV isolates clustered together as a discrete group within the Lyssavirus genera closely related to rabies virus and European bat lyssavirus-2 isolates. The ABLV lineage consisted of two variants; one (ybst-ABLV) consisted of isolates only from YBST bats, the other (pteropid-ABLV) was common to Black, Grey-headed and Little Red flying foxes. No associations were found between the sequences and either the geographical location or year found, or individual flying fox species. In Chapter 5, 15 inocula prepared from the brains or salivary glands of naturally-infected bats were evaluated by intracerebral (IC) and footpad (FP) inoculation of Quackenbush mice in order to select and characterize a highly virulent inoculum for further use in bats (Inoculum 5). In Chapter 6, nine Grey-headed flying foxes were inoculated with 105.2 to 105.5 MICED50 of Inoculum 5 divided into four sites, left footpad, pectoral muscle, temporal muscle and muzzle. Another bat was inoculated with half this dose divided into the footpad and pectoral muscle only. Seven of 10 bats developed clinical disease of 1 to 4 days duration between PI-days 10 and 19 and were shown to be ABL-positive by FAT, HAM immunoperoxidase staining, virus isolation in mice, and TaqMan PCR. Five of the seven bats displayed overt aggression, one died during a seizure, and one showed intractable agitation, pacing, tremors, and ataxia. Viral antigen was demonstrated throughout the central and peripheral nervous systems and in the epithelial cells of the submandibular salivary glands (n=4). All affected bats had mild to moderate non-suppurative meningoencephalitis and severe ganglioneuritis. No ABLV was detected in three bats that remained well until the end of the experiment on day 82. One survivor developed a strong but transient antibody response. In Chapter 7, the relative virulence of inocula prepared from the brains and salivary glands of experimentally infected flying foxes was evaluated in mice by IC and FP inoculation and TaqMan assay. The effects in mice were correlated to the TaqMan CT value and indicated a crude association between virulence and CT value that has potential application in the selection of inocula. In Chapter 8, 36 Black and Grey-headed flying foxes were vaccinated with one (day 0) or two (+ day 28) doses of Nobivac rabies vaccine and co-vaccinated with keyhole limpet haemocyanin (KLH). All bats responded to the Nobivac vaccine with a rabies-RFFIT titer > 0.5 IU/mL that is nominally indicative of protective immunity. Plasma from bats with rabies titres >2 IU/mL had cross-neutralising ABLV titres >1:154. A specifically developed ELISA detected a strong but transient response to KLH.

Environment characterization as an aid to wheat improvement: interpreting genotype-environment interactions by modelling water-deficit patterns in North-Eastern Australia.

Genotype-environment interactions (GEI) limit genetic gain for complex traits such as tolerance to drought. Characterization of the crop environment is an important step in understanding GEI. A modelling approach is proposed here to characterize broadly (large geographic area, long-term period) and locally (field experiment) drought-related environmental stresses, which enables breeders to analyse their experimental trials with regard to the broad population of environments that they target. Water-deficit patterns experienced by wheat crops were determined for drought-prone north-eastern Australia, using the APSIM crop model to account for the interactions of crops with their environment (e.g. feedback of plant growth on water depletion). Simulations based on more than 100 years of historical climate data were conducted for representative locations, soils, and management systems, for a check cultivar, Hartog. The three main environment types identified differed in their patterns of simulated water stress around flowering and during grain-filling. Over the entire region, the terminal drought-stress pattern was most common (50% of production environments) followed by a flowering stress (24%), although the frequencies of occurrence of the three types varied greatly across regions, years, and management. This environment classification was applied to 16 trials relevant to late stages testing of a breeding programme. The incorporation of the independently-determined environment types in a statistical analysis assisted interpretation of the GEI for yield among the 18 representative genotypes by reducing the relative effect of GEI compared with genotypic variance, and helped to identify opportunities to improve breeding and germplasm-testing strategies for this region.

Stock assessment of the Queensland east coast common coral trout (Plectropomus leopardus) fishery

Common coral trout Plectropomus leopardus is an iconic fish of the Great Barrier Reef (GBR) and is the most important fish for the commercial fishery there. Most of the catch is exported live to Asia. This stock assessment was undertaken in response to falls in catch sizes and catch rates in recent years, in order to gauge the status of the stock. It is the first stock assessment ever conducted of coral trout on the GBR, and brings together a multitude of different data sources for the first time. The GBR is very large and was divided into a regional structure based on the Bioregions defined by expert committees appointed by the Great Barrier Reef Marine Park Authority (GBRMPA) as part of the 2004 rezoning of the GBR. The regional structure consists of six Regions, from the Far Northern Region in the north to the Swains and Capricorn–Bunker Regions in the south. Regions also closely follow the boundaries between Bioregions. Two of the northern Regions are split into Subregions on the basis of potential changes in fishing intensity between the Subregions; there are nine Subregions altogether, which include four Regions that are not split. Bioregions are split into Subbioregions along the Subregion boundaries. Finally, each Subbioregion is split into a “blue” population which is open to fishing and a “green” population which is closed to fishing. The fishery is unusual in that catch rates as an indicator of abundance of coral trout are heavily influenced by tropical cyclones. After a major cyclone, catch rates fall for two to three years, and rebound after that. This effect is well correlated with the times of occurrence of cyclones, and usually occurs in the same month that the cyclone strikes. However, statistical analyses correlating catch rates with cyclone wind energy did not provide significantly different catch rate trends. Alternative indicators of cyclone strength may explain more of the catch rate decline, and future work should investigate this. Another feature of catch rates is the phenomenon of social learning in coral trout populations, whereby when a population of coral trout is fished, individuals quickly learn not to take bait. Then the catch rate falls sharply even when the population size is still high. The social learning may take place by fish directly observing their fellows being hooked, or perhaps heeding a chemo-sensory cue emitted by fish that are hooked. As part of the assessment, analysis of data from replenishment closures of Boult Reef in the Capricorn–Bunker Region (closed 1983–86) and Bramble Reef in the Townsville Subregion (closed 1992–95) estimated a strong social learning effect. A major data source for the stock assessment was the large collection of underwater visual survey (UVS) data collected by divers who counted the coral trout that they sighted. This allowed estimation of the density of coral trout in the different Bioregions (expressed as a number of fish per hectare). Combined with mapping data of all the 3000 or so reefs making up the GBR, the UVS results provided direct estimates of the population size in each Subbioregion. A regional population dynamic model was developed to account for the intricacies of coral trout population dynamics and catch rates. Because the statistical analysis of catch rates did not attribute much of the decline to tropical cyclones, (and thereby implied “real” declines in biomass), and because in contrast the UVS data indicate relatively stable population sizes, model outputs were unduly influenced by the unlikely hypothesis that falling catch rates are real. The alternative hypothesis that UVS data are closer to the mark and declining catch rates are an artefact of spurious (e.g., cyclone impact) effects is much more probable. Judging by the population size estimates provided by the UVS data, there is no biological problem with the status of coral trout stocks. The estimate of the total number of Plectropomus leopardus on blue zones on the GBR in the mid-1980s (the time of the major UVS series) was 5.34 million legal-sized fish, or about 8400 t exploitable biomass, with an 2 additional 3350 t in green zones (using the current zoning which was introduced on 1 July 2004). For the offshore regions favoured by commercial fishers, the figure was about 4.90 million legal-sized fish in blue zones, or about 7700 t exploitable biomass. There is, however, an economic problem, as indicated by relatively low catch rates and anecdotal information provided by commercial fishers. The costs of fishing the GBR by hook and line (the only method compatible with the GBR’s high conservation status) are high, and commercial fishers are unable to operate profitably when catch rates are depressed (e.g., from a tropical cyclone). The economic problem is compounded by the effect of social learning in coral trout, whereby catch rates fall rapidly if fishers keep returning to the same fishing locations. In response, commercial fishers tend to spread out over the GBR, including the Far Northern and Swains Regions which are far from port and incur higher travel costs. The economic problem provides some logic to a reduction in the TACC. Such a reduction during good times, such as when the fishery is rebounding after a major tropical cyclone, could provide a net benefit to the fishery, as it would provide a margin of stock safety and make the fishery more economically robust by providing higher catch rates during subsequent periods of depressed catches. During hard times when catch rates are low (e.g., shortly after a major tropical cyclone), a change to the TACC would have little effect as even a reduced TACC would not come close to being filled. Quota adjustments based on catch rates should take account of long-term trends in order to mitigate variability and cyclone effects in data.

Probiotic Bacillus amyloliquefaciens Strain H57 Improves the Performance of Pregnant and Lactating Ewes Fed a Diet Based on Palm Kernel Meal

Probiotic supplements are single or mixed strain cultures of live microorganisms that benefit the host by improving the properties of the indigenous microflora (Seo et al 2010). In a pilot study at the University of Queensland, Norton et al (2008) found that Bacillus amyloliquefaciens Strain H57 (H57), primarily investigated as an inoculum to make high-quality hay, improved feed intake and nitrogen utilisation over several weeks in pregnant ewes. The purpose of the following study was to further challenge the potential of H57 -to show it survives the steam-pelleting process, and that it improves the performance of ewes fed pellets based on an agro-industrial by-product with a reputation for poor palatability, palm kernel meal (PKM), (McNeill 2013). Thirty-two first-parity White Dorper ewes (day 37 of pregnancy, mean liveweight = 47.3 kg, mean age = 15 months) were inducted into individual pens in the animal house at the University of Queensland, Gatton. They were adjusted onto PKM-based pellets (g/kg drymatter (DM): PKM, 408; sorghum, 430; chick pea hulls, 103; minerals and vitamins; Crude protein, 128; ME: 11.1MJ/kg DM) until day 89 of pregnancy and thereafter fed a predominately pelleted diet incorporating with or without H57 spores (10 9 colony forming units (cfu)/kg pellet, as fed), plus 100g/ewe/day oaten chaff, until day 7 of lactation. From day 7 to 20 of lactation the pelleted component of the diet was steadily reduced to be replaced by a 50:50 mix of lucerne: oaten chaff, fed ad libitum, plus 100g/ewe/day of ground sorghum grain with or without H57 (10 9 cfu/ewe/day). The period of adjustment in pregnancy (day 37-89) extended beyond expectations due to some evidence of mild ruminal acidosis after some initially high intakes that were followed by low intakes. During that time the diet was modified, in an attempt to improve palatability, by the addition of oaten chaff and the removal of an acidifying agent (NH4Cl) that was added initially to reduce the risk of urinary calculi. Eight ewes were removed due to inappetence, leaving 24 ewes to start the trial at day 90 of pregnancy. From day 90 of pregnancy until day 63 of lactation, liveweights of the ewes and their lambs were determined weekly and at parturition. Feed intakes of the ewes were determined weekly. Once lambing began, 1 ewe was removed as it gave birth to twin lambs (whereas the rest gave birth to a single lamb), 4 due to the loss of their lambs (2 to dystocia), and 1 due to copper toxicity. The PKM pellets were suspected to be the cause of the copper toxicity and so were removed in early lactation. Hence, the final statistical analysis using STATISTICA 8 (Repeated measures ANOVA for feed intake, One-way ANOVA for liveweight change and birth weight) was completed on 23 ewes for the pregnancy period (n = 11 fed H57; n = 12 control), and 18 ewes or lambs for the lactation period (n = 8 fed H57; n = 10 control). From day 90 of pregnancy until parturition the H57 supplemented ewes ate 17 more DM (g/day: 1041 vs 889, sed = 42.4, P = 0.04) and gained more liveweight (g/day: 193 vs 24.0, sed = 25.4, P = 0.0002), but produced lambs with a similar birthweight (kg: 4.18 vs 3.99, sed = 0.19, P = 0.54). Over the 63 days of lactation the H57 ewes ate similar amounts of DM but grew slower than the control ewes (g/day: 1.5 vs 97.0, sed = 21.7, P = 0.012). The lambs of the H57 ewes grew faster than those of the control ewes for the first 21 days of lactation (g/day: 356 vs 265, sed = 16.5, P = 0.006). These data support the findings of Norton et al (2008) and Kritas et al (2006) that certain Bacillus spp. supplements can improve the performance of pregnant and lactating ewes. In the current study we particularly highlighted the capacity of H57 to stimulate immature ewes to continue to grow maternal tissue through pregnancy, possibly through an enhanced appetite, which appeared then to stimulate a greater capacity to partition nutrients to their lambs through milk, at least for the first few weeks of lactation, a critical time for optimising lamb survival. To conclude, H57 can survive the steam pelleting process to improve feed intake and maternal liveweight gain in late pregnancy, and performance in early lactation, of first-parity ewes fed a diet based on PKM.

On-farm Campylobacter and Escherichia coli in commercial broiler chickens: Re-used bedding does not influence Campylobacter emergence and levels across sequential farming cycles

Limitations in quality bedding material have resulted in the growing need to re-use litter during broiler farming in some countries, which can be of concern from a food-safety perspective. The aim of this study was to compare the Campylobacter levels in ceca and litter across three litter treatments under commercial farming conditions. The litter treatments were (a) the use of new litter after each farming cycle; (b) an Australian partial litter re-use practice; and (c) a full litter re-use practice. The study was carried out on two farms over two years (Farm 1, from 2009–2010 and Farm 2, from 2010–2011), across three sheds (35,000 to 40,000 chickens/shed) on each farm, adopting three different litter treatments across six commercial cycles. A random sampling design was adopted to test litter and ceca for Campylobacter and Escherichia coli, prior to commercial first thin-out and final pick-up. Campylobacter levels varied little across litter practices and farming cycles on each farm and were in the range of log 8.0–9.0 CFU/g in ceca and log 4.0–6.0 MPN/g for litter. Similarly the E. coli in ceca were ∼log 7.0 CFU/g. At first thin-out and final pick-up, the statistical analysis for both litter and ceca showed that the three-way interaction (treatments by farms by times) was highly significant (P < 0.01), indicating that the patterns of Campylobacter emergence/presence across time vary between the farms, cycles and pickups. The emergence and levels of both organisms were not influenced by litter treatments across the six farming cycles on both farms. Either C. jejuni or C. coli could be the dominant species across litter and ceca, and this phenomenon could not be attributed to specific litter treatments. Irrespective of the litter treatments in place, cycle 2 on Farm 2 remained campylobacter-free. These outcomes suggest that litter treatments did not directly influence the time of emergence and levels of Campylobacter and E. coli during commercial farming.