Effect of water quality and season on the population dynamics of Cabomba caroliniana in subtropical Queensland, Australia

Cabomba caroliniana is a submersed aquatic macrophyte that originates from the Americas and is currently invading temperate, subtropical, and tropical freshwater habitats around the world. Despite being a nuisance in many countries, little is known about its ecology. We monitored C. caroliniana populations in three reservoirs in subtropical Queensland, Australia, over 5.5 years. Although biomass, stem length, and plant density of the C. caroliniana stands fluctuated over time, they did not exhibit clear seasonal patterns. Water depth was the most important environmental factor explaining C. caroliniana abundance. Plant biomass was greatest at depths from 2–4 m and rooted plants were not found beyond 5 m. Plant density was greatest in shallow water and decreased with depth, most likely as a function of decreasing light and increasing physical stress. We tested the effect of a range of water physico-chemical parameters. The concentration of phosphorus in the water column was the variable that explained most of the variation in C. caroliniana population parameters. We found that in subtropical Australia, C. caroliniana abundance does not appear to be affected by seasonal conditions but is influenced by other environmental variables such as water depth and nutrient loading. Therefore, further spread will more likely be governed by local habitat rather than climatic conditions.

Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

Extremely Low Microsatellite Diversity but Distinct Population Structure in a Long-Lived Threatened Species, the Australian Lungfish Neoceratodus forsteri (Dipnoi)

The Australian lungfish is a unique living representative of an ancient dipnoan lineage, listed as ‘vulnerable’ to extinction under Australia’s Environment Protection and Biodiversity Conservation Act 1999. Historical accounts indicate this species occurred naturally in two adjacent river systems in Australia, the Burnett and Mary. Current day populations in other rivers are thought to have arisen by translocation from these source populations. Early genetic work detected very little variation and so had limited power to answer questions relevant for management including how genetic variation is partitioned within and among sub-populations. In this study, we use newly developed microsatellite markers to examine samples from the Burnett and Mary Rivers, as well as from two populations thought to be of translocated origin, Brisbane and North Pine. We test whether there is significant genetic structure among and within river drainages; assign putatively translocated populations to potential source populations; and estimate effective population sizes. Eleven polymorphic microsatellite loci genotyped in 218 individuals gave an average within-population heterozygosity of 0.39 which is low relative to other threatened taxa and for freshwater fishes in general. Based on F_ST values (average over loci = 0.11) and STRUCTURE analyses, we identify three distinct populations in the natural range, one in the Burnett and two distinct populations in the Mary. These analyses also support the hypothesis that the Mary River is the likely source of translocated populations in the Brisbane and North Pine rivers, which agrees with historical published records of a translocation event giving rise to these populations. We were unable to obtain bounded estimates of effective population size, as we have too few genotype combinations, although point estimates were low, ranging from 29 - 129. We recommend that, in order to preserve any local adaptation in the three distinct populations that they be managed separately.

Joint Regression and Association Models for Repeated Categorical Responses

The focus of this study is on statistical analysis of categorical responses, where the response values are dependent of each other. The most typical example of this kind of dependence is when repeated responses have been obtained from the same study unit. For example, in Paper I, the response of interest is the pneumococcal nasopharengyal carriage (yes/no) on 329 children. For each child, the carriage is measured nine times during the first 18 months of life, and thus repeated respones on each child cannot be assumed independent of each other. In the case of the above example, the interest typically lies in the carriage prevalence, and whether different risk factors affect the prevalence. Regression analysis is the established method for studying the effects of risk factors. In order to make correct inferences from the regression model, the associations between repeated responses need to be taken into account. The analysis of repeated categorical responses typically focus on regression modelling. However, further insights can also be gained by investigating the structure of the association. The central theme in this study is on the development of joint regression and association models. The analysis of repeated, or otherwise clustered, categorical responses is computationally difficult. Likelihood-based inference is often feasible only when the number of repeated responses for each study unit is small. In Paper IV, an algorithm is presented, which substantially facilitates maximum likelihood fitting, especially when the number of repeated responses increase. In addition, a notable result arising from this work is the freely available software for likelihood-based estimation of clustered categorical responses.

Optimising initial population density, growth time and nitrogen nutrition for assessing resistance of wheat cultivars to root-lesion nematode (Pratylenchus thornei)

Pratylenchus thornei is a major pathogen of wheat in Australia. Two glasshouse experiments with four wheat cultivars that had different final populations (Pf) of P. thornei in the field were used to optimise conditions for assessing resistance. With different initial populations (Pi) ranging up to 5250 P. thornei/kg soil, Pf of P. thornei increased to 16 weeks after sowing, and then decreased at 20 weeks in some cultivar x Pi combinations. The population dynamics of P. thornei up to 16 weeks were best described by a modified exponential equation P f (t) = aP i e kt where P f (t) is the final population density at time t, P i is the initial population density, a is the proportion of P i that initiates population development, and k is the intrinsic rate of increase of the population. The cultivar GS50a had very low k values at Pi of 5250 and 1050 indicating its resistance, Suneca and Potam had high k values indicating susceptibility, whereas intolerant Gatcher had a low value at the higher Pi and a high value at the lower Pi. Nitrate fertiliser increased plant growth and Pf values of susceptible cultivars, but in unplanted soil it decreased Pf. Nematicide (aldicarb 5 mg/kg soil) killed P. thornei more effectively in planted than in unplanted soil and increased plant growth particularly in the presence of N fertiliser. In both experiments, the wheat cultivars Suneca and Potam were more susceptible than the cultivar GS50a reflecting field results. The method chosen to discriminate wheat cultivars was to assess Pf after growth for 16 weeks in soil with Pi ~1050–5250 P. thornei/kg soil and fertilised with 200 mg NO3–N/kg soil.

Integrating pest population models with biophysical crop models to better represent the farming system

Farming systems frameworks such as the Agricultural Production Systems simulator (APSIM) represent fluxes through the soil, plant and atmosphere of the system well, but do not generally consider the biotic constraints that function within the system. We designed a method that allowed population models built in DYMEX to interact with APSIM. The simulator engine component of the DYMEX population-modelling platform was wrapped within an APSIM module allowing it to get and set variable values in other APSIM models running in the simulation. A rust model developed in DYMEX is used to demonstrate how the developing rust population reduces the crop's green leaf area. The success of the linking process is seen in the interaction of the two models and how changes in rust population on the crop's leaves feedback to the APSIM crop modifying the growth and development of the crop's leaf area. This linking of population models to simulate pest populations and biophysical models to simulate crop growth and development increases the complexity of the simulation, but provides a tool to investigate biotic constraints within farming systems and further moves APSIM towards being an agro-ecological framework.

Performance and physiological responses of combined t-bar and PIT tagged giant mud crabs (Scylla serrata)

Mud crabs (Scylla spp.) are intensively caught throughout South-East Asia and support a very substantial commercial, recreational fishing and aquaculture industry. Identification of individual animals is important to improve understanding and management of this species. However, tagging of crustaceans is difficult as they frequently molt and internal tags can pose a hazard to consumers. In this pilot study we tested a new method combining passive integrated transponder tags and t-bar tags externally. 45 giant mud crabs (Scylla serrata) were captured from the wild and kept in tanks for a maximum of 10 months. We inserted tags into the abdomen of 35 giant mud crabs and tested a modified method where the combined t-bar/PIT-tag was inserted into the muscle tissue of the rear leg between the dorsal carapace plate and the top of the abdominal flap. Tagged crabs with the modified method showed 85% tag retention for molting crabs. We tested the same method in the field where 852 individuals were tagged with combined t-bar/PIT-tags of which 82 were recaptured showing 100% tag retention but without any evidence of molting having occurred. The tested method of combined t-bar/PIT-tags in giant mud crabs can further improve monitoring for wild and aquaculture populations and can be deployed widely with low cost.

Performance and physiological responses of combined t-bar and PIT tagged giant mud crabs (Scylla serrata)

Mud crabs (Scylla spp.) are intensively caught throughout South-East Asia and support a very substantial commercial, recreational fishing and aquaculture industry. Identification of individual animals is important to improve understanding and management of this species. However, tagging of crustaceans is difficult as they frequently molt and internal tags can pose a hazard to consumers. In this pilot study we tested a new method combining passive integrated transponder tags and t-bar tags externally. 45 giant mud crabs (Scylla serrata) were captured from the wild and kept in tanks for a maximum of 10 months. We inserted tags into the abdomen of 35 giant mud crabs and tested a modified method where the combined t-bar/PIT-tag was inserted into the muscle tissue of the rear leg between the dorsal carapace plate and the top of the abdominal flap. Tagged crabs with the modified method showed 85% tag retention for molting crabs. We tested the same method in the field where 852 individuals were tagged with combined t-bar/PIT-tags of which 82 were recaptured showing 100% tag retention but without any evidence of molting having occurred. The tested method of combined t-bar/PIT-tags in giant mud crabs can further improve monitoring for wild and aquaculture populations and can be deployed widely with low cost.

The study of genetic diversity and population structure of the Korean fleshy shrimp, Fenneropenaeus chinensis, using newly developed microsatellite markers

The fleshy shrimp, Fenneropenaeus chinensis, is the family of Penaeidae and one of the most economically important marine culture species in Korea. However, its genetic characteristics have never been studied. In this study, a total of 240 wild F. chinensis individuals were collected from four locations as follows: Narodo (NRD, n = 60), Beopseongpo (BSP, n = 60), Chaesukpo (CSP, n = 60), and Cheonsuman (CSM, n = 60). Genetic variability and the relationships among four wild F. chinensis populations were analyzed using 13 newly developed microsatellite loci. Relatively high levels of genetic variability (mean allelic richness = 16.87; mean heterozygosity = 0.845) were found among localities. Among the 52 population loci, 13 showed significant deviation from the Hardy–Weinberg equilibrium. Neighbor-joining, principal coordinate, and molecular variance analyses revealed the presence of three subpopulations (NRD, CSM, BSP and CSP), which was consistent with clustering based on genetic distance. The mean observed heterozygosity values of the NRD, CSM, BSP, and CSP populations were 0.724, 0.821, 0.814, and 0.785 over all loci, respectively. These genetic variability and differentiation results of the four wild populations can be applied for future genetic improvement using selective breeding and to design suitable management guidelines for Korean F. chinensis culture.

Preliminary investigation of some physiological responses of Bos indicus heifers to surgical spaying

Objective To determine the value of peripheral blood concentrations of cortisol, creatine phosphokinase (CPK), aspartate aminotransferase (AST), non-esterified fatty acids (NEFAs) and haptoglobin as indicators of welfare in Brahman heifers spayed by either the Willis dropped ovary technique (WDOT) or the flank laparotomy method. Design A total of 24, 2-year-old Brahman heifers were allocated to: crush (head-bail) restraint alone (Control, n = 5); crush restraint and ear-punch (Ear-punch, n = 5); crush restraint, WDOT spay and ear-punch (WDOT, n = 9); or crush restraint, elecrtoimmobilisation, flank spay and ear-punch (Flank; n = 5). Cattle were blood sampled frequently to 8 h, and then daily to day 4 and were monitored to 42 days post-procedure. Peripheral blood concentrations of bound and unbound cortisol, CPK, AST, NEFAs and haptoglobin were determined. Results Concentrations of plasma bound cortisol peaked in the spayed heifers 3-4 h post-procedure; values in the Flank (1603 nmol/L) and WDOT (1290 nmol/L) groups were similar and significantly greater (P < 0.05) than in the Controls (519 nmol/L). Flank heifers had elevated plasma haptoglobin levels to day 4 postprocedure. Liveweights were significantly lower in the spayed compared with the Control heifers at 21 and 42 days post-procedure, with liveweight gains also significantly reduced at day 21. Conclusions Bound cortisol responses in spayed heifers were elevated to 6 h post-procedure and similar in WDOT- and flank-spayed animals, indicating comparable levels of pain and stress. An inflammatory response, indicated by haptoglobin concentrations, was sustained for longer in Flank than in WDOT spayed heifers, suggesting longer-lasting adverse effects on welfare from flank spaying than WDOT spaying. © 2011 The State of Queensland (Department of Employment, Economic Development and Innovation). Australian Veterinary Journal © 2011 Australian Veterinary Association.

Genetic characterization of field-evolved resistance to phosphine in the rusty grain beetle, Cryptolestes ferrugineus (Laemophloeidae: Coleoptera)

Inheritance of resistance to phosphine fumigant was investigated in three field-collected strains of rusty grain beetle, Cryptolestes ferrugineus, Susceptible (S-strain), Weakly Resistant (Weak-R) and Strongly Resistant (Strong-R). The strains were purified for susceptibility, weak resistance and strong resistance to phosphine, respectively, to ensure homozygosity of resistance genotype. Crosses were established between S-strain × Weak-R, S-strain × Strong-R and Weak-R × Strong-R, and the dose mortality responses to phosphine of these strains and their F1, F2 and F1-backcross progeny were obtained. The fumigations were undertaken at 25 °C and 55% RH for 72 h. Weak-R and Strong-R showed resistance factors of 6.3 × and 505 × compared with S-strain at the LC50. Both weak and strong resistances were expressed as incompletely recessive with degrees of dominance of − 0.48 and − 0.43 at the LC50, respectively. Responses of F2 and F1-backcross progeny indicated the existence of one major gene in Weak-R, and at least two major genes in Strong-R, one of which was allelic with the major factor in Weak-R. Phenotypic variance analyses also estimated that the number of independently segregating genes conferring weak resistance was 1 (nE = 0.89) whereas there were two genes controlling strong resistance (nE = 1.2). The second gene, unique to Strong-R, interacted synergistically with the first gene to confer a very high level of resistance (~ 80 ×). Neither of the two major resistance genes was sex linked. Despite the similarity of the genetics of resistance to that previously observed in other pest species, a significant proportion (~ 15 to 30%) of F1 individuals survived at phosphine concentrations higher than predicted. Thus it is likely that additional dominant heritable factors, present in some individuals in the population, also influenced the resistance phenotype. Our results will help in understanding the process of selection for phosphine resistance in the field which will inform resistance management strategies. In addition, this information will provide a basis for the identification of the resistance genes.

Genetic characterization of field-evolved resistance to phosphine in the rusty grain beetle, Cryptolestes ferrugineus (Laemophloeidae: Coleoptera)

Inheritance of resistance to phosphine fumigant was investigated in three field-collected strains of rusty grain beetle, Cryptolestes ferrugineus, Susceptible (S-strain), Weakly Resistant (Weak-R) and Strongly Resistant (Strong-R). The strains were purified for susceptibility, weak resistance and strong resistance to phosphine, respectively, to ensure homozygosity of resistance genotype. Crosses were established between S-strain × Weak-R, S-strain × Strong-R and Weak-R × Strong-R, and the dose mortality responses to phosphine of these strains and their F1, F2 and F1-backcross progeny were obtained. The fumigations were undertaken at 25 °C and 55% RH for 72 h. Weak-R and Strong-R showed resistance factors of 6.3 × and 505 × compared with S-strain at the LC50. Both weak and strong resistances were expressed as incompletely recessive with degrees of dominance of − 0.48 and − 0.43 at the LC50, respectively. Responses of F2 and F1-backcross progeny indicated the existence of one major gene in Weak-R, and at least two major genes in Strong-R, one of which was allelic with the major factor in Weak-R. Phenotypic variance analyses also estimated that the number of independently segregating genes conferring weak resistance was 1 (nE = 0.89) whereas there were two genes controlling strong resistance (nE = 1.2). The second gene, unique to Strong-R, interacted synergistically with the first gene to confer a very high level of resistance (~ 80 ×). Neither of the two major resistance genes was sex linked. Despite the similarity of the genetics of resistance to that previously observed in other pest species, a significant proportion (~ 15 to 30%) of F1 individuals survived at phosphine concentrations higher than predicted. Thus it is likely that additional dominant heritable factors, present in some individuals in the population, also influenced the resistance phenotype. Our results will help in understanding the process of selection for phosphine resistance in the field which will inform resistance management strategies. In addition, this information will provide a basis for the identification of the resistance genes.

Responses of fruit flies (Tephritidae: Dacinae) to novel male attractants in north Queensland, Australia, and improved lures for some pest species

Male fruit fly attractants, cue-lure and methyl eugenol (ME), have been successfully used for the last 50 years in the monitoring and control of Dacini fruit flies (Bactrocera and Dacus species). However, over 50% of Dacini are non-responsive to either lure, including some pest species. A new lure, zingerone, has been found to weakly attract cue- and ME-responsive species in Malaysia. In Australia it attracted a weakly cue-responsive minor pest Bactrocera jarvisi (Tryon) and three non-responsive' species. Similar compounds were tested in Queensland and attracted cue- and ME-responsive species and two non-responsive' species. In this study, 14 novel compounds, including raspberry ketone formate (RKF) (Melolure) and zingerone, were field tested in comparison with cue-lure and ME at 17 sites in north Queensland. The most attractive novel lures were isoeugenol, methyl-isoeugenol, dihydroeugenol and zingerone. Several non-responsive' species responded to the new lures: Bactrocera halfordiae (Tryon), a species of some market access concern, was most attracted to isoeugenol; B.barringtoniae (Tryon), B.bidentata (May) and B.murrayi (Perkins) responded to isoeugenol, methyl-isoeugenol and dihydroeugenol; two new species of Dacus responded to zingerone. Bactrocera kraussi (Hardy), a cue-responsive minor pest in north Queensland, was significantly more attracted to isoeugenol than cue-lure. The cue-responsive D.absonifacies (May) and D.secamoneaeDrew were significantly more attracted to zingerone than cue-lure. Bactrocera yorkensisDrew & Hancock, a ME-responsive species was significantly more attracted to isoeugenol, methyl-isoeugenol and dihydroeugenol than ME. The preferential response to RKF or cue-lure was species specific. Six species were significantly more attracted to RKF, including the pests B.tryoni (Froggatt), B.frauenfeldi (Schiner) and minor pest B.bryoniae (Tryon); eight species were significantly more attracted to cue-lure including the pest B.neohumeralis (Hardy). These findings have significance in the search for optimal male lures for pest species elsewhere in the world.

Silica Vesicle Nanovaccine Formulations Stimulate Long-Term Immune Responses to the Bovine Viral Diarrhoea Virus E2 Protein

Bovine Viral Diarrhoea Virus (BVDV) is one of the most serious pathogen, which causes tremendous economic loss to the cattle industry worldwide, meriting the development of improved subunit vaccines. Structural glycoprotein E2 is reported to be a major immunogenic determinant of BVDV virion. We have developed a novel hollow silica vesicles (SV) based platform to administer BVDV-1 Escherichia coli-expressed optimised E2 (oE2) antigen as a nanovaccine formulation. The SV-140 vesicles (diameter 50 nm, wall thickness 6 nm, perforated by pores of entrance size 16 nm and total pore volume of 0.934 cm(3)g(-1)) have proven to be ideal candidates to load oE2 antigen and generate immune response. The current study for the first time demonstrates the ability of freeze-dried (FD) as well as non-FD oE2/SV140 nanovaccine formulation to induce long-term balanced antibody and cell mediated memory responses for at least 6 months with a shortened dosing regimen of two doses in small animal model. The in vivo ability of oE2 (100 mu g)/SV-140 (500 mu g) and FD oE2 (100 mu g)/SV-140 (500 mu g) to induce long-term immunity was compared to immunisation with oE2 (100 mu g) together with the conventional adjuvant Quil-A from the Quillaja saponira (10 mu g) in mice. The oE2/SV-140 as well as the FD oE2/SV-140 nanovaccine generated oE2-specific antibody and cell mediated responses for up to six months post the final second immunisation. Significantly, the cell-mediated responses were consistently high in mice immunised with oE2/SV-140 (1,500 SFU/million cells) at the six-month time point. Histopathology studies showed no morphological changes at the site of injection or in the different organs harvested from the mice immunised with 500 mu g SV-140 nanovaccine compared to the unimmunised control. The platform has the potential for developing single dose vaccines without the requirement of cold chain storage for veterinary and human applications.

Comparison of the impact of six heat-load management strategies on thermal responses and milk production of feed-pad and pasture fed dairy cows in a subtropical environment

Exposure to hot environments affects milk yield (MY) and milk composition of pasture and feed-pad fed dairy cows in subtropical regions. This study was undertaken during summer to compare MY and physiology of cows exposed to six heat-load management treatments. Seventy-eight Holstein-Friesian cows were blocked by season of calving, parity, milk yield, BW, and milk protein (%) and milk fat (%) measured in 2 weeks prior to the start of the study. Within blocks, cows were randomly allocated to one of the following treatments: open-sided iron roofed day pen adjacent to dairy (CID) + sprinklers (SP); CID only; non-shaded pen adjacent to dairy + SP (NSD + SP); open-sided shade cloth roofed day pen adjacent to dairy (SCD); NSD + sprinkler (sprinkler on for 45 min at 1100 h if mean respiration rate >80 breaths per minute (NSD + WSP)); open-sided shade cloth roofed structure over feed bunk in paddock + 1 km walk to and from the dairy (SCP + WLK). Sprinklers for CID + SP and NSD + SP cycled 2 min on, 12 min off when ambient temperature >26°C. The highest milk yields were in the CID + SP and CID treatments (23.9 L cow−1 day−1), intermediate for NSD + SP, SCD and SCP + WLK (22.4 L cow−1 day−1), and lowest for NSD + WSP (21.3 L cow−1 day−1) (P < 0.05). The highest (P < 0.05) feed intakes occurred in the CID + SP and CID treatments while intake was lowest (P < 0.05) for NSD + WSP and SCP + WLK. Weather data were collected on site at 10-min intervals, and from these, THI was calculated. Nonlinear regression modelling of MY × THI and heat-load management treatment demonstrated that cows in CID + SP showed no decline in MY out to a THI break point value of 83.2, whereas the pooled MY of the other treatments declined when THI >80.7. A combination of iron roof shade plus water sprinkling throughout the day provided the most effective control of heat load.