Determining changes in the distribution and abundance of a Rhyzopertha dominica phosphine resistance allele in farm grain storages using a DNA marker

BACKGROUND: The lesser grain borer, Rhyzopertha dominica (F.), is a highly destructive pest of stored grain that is strongly resistant to the fumigant phosphine (PH3). Phosphine resistance is due to genetic variants at the rph2 locus that alter the function of the dihydrolipoamide dehydrogenase (DLD) gene. This discovery now enables direct detection of resistance variants at the rph2 locus in field populations. RESULTS: A genotype assay was developed for direct detection of changes in distribution and frequency of a phosphine resistance allele in field populations of R. dominica. Beetles were collected from ten farms in south-east Queensland in 2006 and resampled in 2011. Resistance allele frequency increased in the period from 2006 to 2011 on organic farms with no history of phosphine use, implying that migration of phosphine-resistant R. dominica had occurred from nearby storages. CONCLUSION: Increasing resistance allele frequencies on organic farms suggest local movement of beetles and dispersal of insects from areas where phosphine has been used. This research also highlighted for the first time the utility of a genetic DNA marker in accurate and rapid determination of the distribution of phosphine-resistant insects in the grain value chain. Extending this research over larger landscapes would help in identifying resistance problems and enable timely pest management decisions. © 2013 Society of Chemical Industry © 2013 Society of Chemical Industry 69 6 June 2013 10.1002/ps.3514 Rapid Report Rapid Report © 2013 Society of Chemical Industry.

Molecular Biomonitoring during Rhizoremediation of Oil-Contaminated Soil

Rhizoremediation is the use of microbial populations present in the rhizosphere of plants for environmental cleanup. The idea of this work was that bacteria living in the rhizosphere of a nitrogen-fixing leguminous plant, goat's rue (Galega orientalis), could take part in the degradation of harmful monoaromatic hydrocarbons, such as benzene, toluene and xylene (BTEX), from oil-contaminated soils. In addition to chemical (e.g. pollutant concentration) and physical (e.g. soil structure) information, the knowledge of biological aspects (e.g. bacteria and their catabolic genes) is essential when developing the rhizoremediation into controlled and effective bioremediation practice. Therefore, the need for reliable biomonitoring methods is obvious. The main aims of this thesis were to evaluate the symbiotic G. orientalis - Rhizobium galegae system for rhizoremediation of oil-contaminated soils, to develop molecular methods for biomonitoring, and to apply these methods for studying the microbiology of rhizoremediation. In vitro, Galega plants and rhizobia remained viable in m-toluate concentrations up to 3000 mg/l. Plant growth and nodulation were inhibited in 500 mg/l m-toluate, but were restored when plants were transferred to clean medium. In the greenhouse, Galega showed good growth, nodulation and nitrogen fixation, and developed a strong rhizosphere in soils contaminated with oil or spiked with 2000 mg/l m-toluate. The high aromatic tolerance of R. galegae and the viability of Galega plants in oil-polluted soils proved this legume system to be a promising method for the rhizoremediation of oil-contaminated soils. Molecular biomonitoring methods were designed and/or developed further for bacteria and their degradation genes. A combination of genomic fingerprinting ((GTG)5-PCR), taxonomic ribotyping of 16S rRNA genes and partial 16S rRNA gene sequencing were chosen for molecular grouping of culturable, heterogeneous rhizosphere bacteria. PCR primers specific for the xylE gene were designed for TOL plasmid detection. Amplified enzyme-coding DNA restriction analysis (AEDRA) with AluI was used to profile both TOL plasmids (xylE primers) and, in general, aromatics-degrading plasmids (C230 primers). The sensitivity of the direct monitoring of TOL plasmids in soil was enhanced by nested C23O-xylE-PCR. Rhizosphere bacteria were isolated from the greenhouse and field lysimeter experiments. High genetic diversity was observed among the 50 isolated, m-toluate tolerating rhizosphere bacteria in the form of five major lineages of the Bacteria domain. Gram-positive Rhodococcus, Bacillus and Arthrobacter and gram-negative Pseudomonas were the most abundant genera. The inoculum Pseudomonas putida PaW85/pWW0 was not found in the rhizosphere samples. Even if there were no ecological niches available for the bioaugmentation bacterium itself, its conjugative catabolic plasmid might have had some additional value for other bacterial species and thus, for rhizoremediation. Only 10 to 20% of the isolated, m-toluate tolerating bacterial strains were also able to degrade m-toluate. TOL plasmids were a major group of catabolic plasmids among these bacteria. The ability to degrade m-toluate by using enzymes encoded by a TOL plasmid was detected only in species of the genus Pseudomonas, and the best m-toluate degraders were these Pseudomonas species. Strain-specific differences in degradation abilities were found for P.oryzihabitans and P. migulae: some of these strains harbored a TOL plasmid - a new finding observed in this work, indicating putative horizontal plasmid transfer in the rhizosphere. One P. oryzihabitans strain harbored the pWW0 plasmid that had probably conjugated from the bioaugmentation Pseudomonas. Some P. migulae and P. oryzihabitans strains seemed to harbor both the pWW0- and the pDK1-type TOL plasmid. Alternatively, they might have harbored a TOL plasmid with both the pWW0- and the pDK1-type xylE gene. The breakdown of m-toluate by gram-negative bacteria was not restricted to the TOL pathway. Also some gram-positive Rhodococcus erythropolis and Arthrobacter aurescens strains were able to degrade m-toluate in the absence of a TOL plasmid. Three aspects of the rhizosphere effect of G. orientalis were manifested in oil-contaminated soil in the field: 1) G. orientalis and Pseudomonas bioaugmentation increased the amount of rhizosphere bacteria. G. orientalis especially together with Pseudomonas bioaugmentation increased the numbers of m-toluate utilizing and catechol positive bacteria indicating an increase in degradation potential. 2) Also the bacterial diversity, when measured as the amount of ribotypes, was increased in the Galega rhizosphere with or without Pseudomonas bioaugmentation. However, the diversity of m-toluate utilizing bacteria did not significantly increase. At the community level, by using the 16S rRNA gene PCR-DGGE method, the highest diversity of species was also observed in vegetated soils compared with non-vegetated soils. Diversified communities may best guarantee the overall success in rhizoremediation by offering various genetic machineries for catabolic processes. 3) At the end of the experiment, no TOL plasmid could be detected by direct DNA analysis in soil treated with both G. orientalis and Pseudomonas. The detection limit for TOL plasmids was encountered indicating decreased amount of degradation plasmids and thus, the success of rhizoremediation. The use of G. orientalis for rhizoremediation is unique. In this thesis new information was obtained about the rhizosphere effect of Galega orientalis in BTEX contaminated soils. The molecular biomonitoring methods can be applied for several purposes within environmental biotechnology, such as for evaluating the intrinsic biodegradation potential, monitoring the enhanced bioremediation, and estimating the success of bioremediation. Environmental protection by using nature's own resources and thus, acting according to the principle of sustainable development, would be both economically and environmentally beneficial for society. Keywords: molecular biomonitoring, genetic fingerprinting, soil bacteria, bacterial diversity, TOL plasmid, catabolic genes, horizontal gene transfer, rhizoremediation, rhizosphere effect, Galega orientalis, aerobic biodegradation, petroleum hydrocarbons, BTEX

Implications of methoprene resistance for managing Rhyzopertha dominica (F.) in stored grain

Commercial formulations of methoprene have been used in a number of countries such as Australia and the USA to provide long-term protection to grain from a range of storage pests. The level of resistance in Rhyzopertha dominica (F.), the lesser grain borer, was investigated in laboratory experiments by direct exposure of adults on treated wheat. Adults of a reference homozygous resistant strain of R.dominica were exposed to treatments of 0, 2, 3, 5, 10, 20 and 40kg-1 of s-methoprene for 7 days. Mortality and progeny production were dose dependent with 98.7% mortality and complete progeny suppression at 40mgkg-1, which is 67 times the registered rate at which s-methoprene is applied as a grain protectant in Australia (0.6mgkg-1) and eight times the rate which has been used in the USA (5mgkg-1). This strain was also tested by adding adults to wheat treated at 0, 1, 3, 10 and 30mgkg-1 and determining the number of adults (progeny plus original parental adults) after 6, 8, 10 or 12 weeks of continuous exposure. The effect of s-methoprene was consistent regardless of the number of weeks of continuous exposure and average population suppression was 99.5% at 30mgkg-1. Screening of 162 field samples collected from southeast Australia in 2009 showed that 93% of samples reproduced when exposed to wheat treated with the Australian registered rate of 0.6mgkg-1. When four unselected resistant field samples from this region were tested, progeny production but not mortality was dose dependent, and all four samples produced live progeny at the highest dose of 30mgkg-1. Our results show that methoprene resistance is a serious threat to the management of R.dominica, and that strategies need to be developed to minimise the further development and spread of resistance. © 2013.

Estimating an Aggregate Economic Genotype to Facilitate Breeding and Selection of Strawberry Genotypes in Southeast Queensland

The Queensland strawberry (Fragaria ×ananassa) breeding program in subtropical Australia aims to improve sustainable profitability for the producer. Selection must account for the relative economic importance of each trait and the genetic architecture underlying these traits in the breeding population. Our study used estimates of the influence of a trait on production costs and profitability to develop a profitability index (PI) and an economic weight (i.e., change in PI for a unit change in level of trait) for each trait. The economic weights were then combined with the breeding values for 12 plant and fruit traits on over 3000 genotypes that were represented in either the current breeding population or as progenitors in the pedigree of these individuals. The resulting linear combination (i.e., sum of economic weight × breeding value for all 12 traits) estimated the overall economic worth of each genotype as H, the aggregate economic genotype. H values were validated by comparisons among commercial cultivars and were also compared with the estimated gross margins. When the H value of ‘Festival’ was set as zero, the H values of genotypes in the pedigree ranged from –0.36 to +0.28. H was highly correlated (R2 = 0.77) with the year of selection (1945–98). The gross margins were highly linearly related (R2 > 0.98) to H values when the genotype was planted on less than 50% of available area, but the relationship was non-linear [quadratic with a maximum (R2 > 0.96)] when the planted area exceeded 50%. Additionally, with H values above zero, the variation in gross margin increased with increasing H values as the percentage of area planted to a genotype increased. High correlations among some traits allowed the omission of any one of three of the 12 traits with little or no effect on ranking (Spearman’s rank correlation 0.98 or greater). Thus, these traits may be dropped from the aggregate economic genotype, leading to either cost reductions in the breeding program or increased selection intensities for the same resources. H was efficient in identifying economically superior genotypes for breeding and deployment, but because of the non-linear relationship with gross margin, calculation of a gross margin for genotypes with high H is also necessary when cultivars are deployed across more than 50% of the available area.

Modelling interactions between farm-level structural adjustment and a regional economy: A case of the Australian rice industry

Climate change and on-going water policy reforms will likely contribute to on-farm and regional structural adjustment in Australia. This paper gathers empirical evidence of farm-level structural adjustments and integrates these with a regional equilibrium model to investigate sectoral and regional impacts of climate change and recent water use policy on rice industry. We find strong evidence of adjustments to the farming system, enabled by existing diversity in on-farm production. A further loss of water with additional pressures to adopt less intensive and larger-scale farming, will however reduce the net number of farm businesses, which may affect regional rice production. The results from a regional CGE model show impacts on the regional economy over and above the direct cost of the environmental water, although a net reduction in real economic output and real income is partially offset by gains in rest of the Australia through the reallocation or resources. There is some interest within the industry and from potential new corporate entrants in the relocation of some rice production to the north. However, strong government support would be crucial to implement such relocation.

Use of pooled samples to assess human exposure to parabens, benzophenone-3 and triclosan in Queensland, Australia

Parabens, benzophenone-3 and triclosan are common ingredients used as preservatives, ultraviolet radiation filters and antimicrobial agents, respectively. Human exposure occurs through consumption of processed food and use of cosmetics and consumer products. The aim of this study was to provide a preliminary characterisation of exposure to selected personal care product chemicals in the general Australian population. De-identified urine specimens stratified by age and sex were obtained from a community-based pathology laboratory and pooled (n= 24 pools of 100). Concentrations of free and total (sum of free plus conjugated) species of methyl, ethyl, propyl and butyl paraben, benzophenone-3 and triclosan were quantified using isotope dilution tandem mass spectrometry; with geometric means 232, 33.5, 60.6, 4.32, 61.5 and 87.7. ng/mL, respectively. Age was inversely associated with paraben concentration, and females had concentrations approximately two times higher than males. Total paraben and benzophenone-3 concentrations are significantly higher than reported worldwide, and the average triclosan concentration was more than one order of magnitude higher than in many other populations. This study provides the first data on exposure of the general Australian population to a range of common personal care product chemical ingredients, which appears to be prevalent and warrants further investigation.

Synthesis of Isoflavone Conjugates

During this study different approaches were studied to obtain isoflavone sulphates, glucuronides and sulphoglucuronides. Three isoflavone disulphates (daidzein-di-O-sulphate, genistein-di-O-sulphate and glycitein-di-O-sulphate) and three isoflavonoid disulphates (dihydrodaidzein-di-O-sulphate, dihydrogenistein-di-O-sulphate and equol-di-O-sulphate) were synthesised in moderate yields by using in situ prepared pyridine sulphur trioxide complex, made from chlorosulphonic acid and pyridine. These disulphated compounds can be used to develop analytical procedures and study the biological activity of disulphated products. As the use of the HPLC-MS methods in the field of isoflavones has increased its popularity, deuterated isoflavone disulphates were synthesised. A new microwave assisted deuteration method, using CF3COOD, was developed for this purpose. Three polydeuterated isoflavone disulphates (daidzein-d6-di-O-sulphate, genistein-d4-di-O-sulphate and glycitein-d6-di-O-sulphate) were obtained in moderate yields with high isotopic purity. A synthetic method was developed for daidzein sulphoglucuronide (daidzein-7-O-b-D-glucuronide-4´-O-sulphate), which is a major metabolite in rat bile. By using protection/deprotection steps, the desired product was finally obtained in moderate yield. The method developed can be used in further studies of synthesis of isoflavonoid mixed conjugates. As a part of this study, the structure of naturally occurring daidzein-4´-O-b-glucoside was verified. Different glycosidation methods are reviewed and possible factors affecting the stereoselectivity are discussed. The study of the selective chlorination of isoflavones was a consequence of the observed unexpected chlorination during the synthesis of isoflavone acid chlorides by thionyl chloride. This fascinating phenomenon was investigated further with various isoflavones and as a result a method for producing isoflavone chlorides (8-chlorogenistein, 6,8-dichlorogenistein and 6,8-dichlorobiochanin A) was developed. Protecting groups played a great role during this study, which led to an intensive study on them. A regioselective protection method was developed by using direct introduction of the protecting group (Benzyl and Benzoyl) to positions 7-O or 4´-O in daidzein, genistein and glycitein with t-BuOK as a base in DMF in moderate yields. The possibility of exploiting the transesterification was also investigated. It was observed that by using K2CO3 as a base in DMF, daidzein, genistein and glycitein could be benzoylated at position 4´-O selectively, in the presence of the more acidic 7 hydroxy group. Transesterification also proved to be useful in the glycosidation of isoflavones at position 7-O, starting from 7-O-benzoylated isoflavones. Different carboxylic acid derivatives were synthesised for use either in the development of radioimmunoassay (7-O-carboxymethylglycitein and 4´-O-carboxymethylglycitein) or synthesis of daunorubicin isoflavone derivative for biological testing (7-O-carboxypropylbiochanin A and 7-O-carboxypropylgenistein).

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.

Spatial variation of climate and the impact of disturbances on local climate and forest recovery in northern Finland

The structure and function of northern ecosystems are strongly influenced by climate change and variability and by human-induced disturbances. The projected global change is likely to have a pronounced effect on the distribution and productivity of different species, generating large changes in the equilibrium at the tree-line. In turn, movement of the tree-line and the redistribution of species produce feedback to both the local and the regional climate. This research was initiated with the objective of examining the influence of natural conditions on the small-scale spatial variation of climate in Finnish Lapland, and to study the interaction and feedback mechanisms in the climate-disturbances-vegetation system near the climatological border of boreal forest. The high (1 km) resolution spatial variation of climate parameters over northern Finland was determined by applying the Kriging interpolation method that takes into account the effect of external forcing variables, i.e., geographical coordinates, elevation, sea and lake coverage. Of all the natural factors shaping the climate, the geographical position, local topography and altitude proved to be the determining ones. Spatial analyses of temperature- and precipitation-derived parameters based on a 30-year dataset (1971-2000) provide a detailed description of the local climate. Maps of the mean, maximum and minimum temperatures, the frost-free period and the growing season indicate that the most favourable thermal conditions exist in the south-western part of Lapland, around large water bodies and in the Kemijoki basin, while the coldest regions are in highland and fell Lapland. The distribution of precipitation is predominantly longitudinally dependent but with the definite influence of local features. The impact of human-induced disturbances, i.e., forest fires, on local climate and its implication for forest recovery near the northern timberline was evaluated in the Tuntsa area of eastern Lapland, damaged by a widespread forest fire in 1960 and suffering repeatedly-failed vegetation recovery since that. Direct measurements of the local climate and simulated heat and water fluxes indicated the development of a more severe climate and physical conditions on the fire-disturbed site. Removal of the original, predominantly Norway spruce and downy birch vegetation and its substitution by tundra vegetation has generated increased wind velocity and reduced snow accumulation, associated with a large variation in soil temperature and moisture and deep soil frost. The changed structural parameters of the canopy have determined changes in energy fluxes by reducing the latter over the tundra vegetation. The altered surface and soil conditions, as well as the evolved severe local climate, have negatively affected seedling growth and survival, leading to more unfavourable conditions for the reproduction of boreal vegetation and thereby causing deviations in the regional position of the timberline. However it should be noted that other factors, such as an inadequate seed source or seedbed, the poor quality of the soil and the intensive logging of damaged trees could also exacerbate the poor tree regeneration. In spite of the failed forest recovery at Tunsta, the position and composition of the timberline and tree-line in Finnish Lapland may also benefit from present and future changes in climate. The already-observed and the projected increase in temperature, the prolonged growing season, as well as changes in the precipitation regime foster tree growth and new regeneration, resulting in an advance of the timberline and tree-line northward and upward. This shift in the distribution of vegetation might be decelerated or even halted by local topoclimatic conditions and by the expected increase in the frequency of disturbances.

Overcoming barriers to the successful implementation of a classical biological control strategy for the exotic invasive weed Mikania micrantha in the Asia-Pacific region

Mikania micrantha (Asteraceae) commonly known as mikania, is a major invasive alien plant (IAP) in the tropical humid agricultural and forest zones of the Asia-Pacific region. This fast-growing Neotropical vine is able to smother plants in agricultural ecosystems, agroforestry and natural habitats, reducing productivity and biodiversity. Fungal pathogens were first investigated for the classical biological control of this weed in 1996. This resulted in the selection and screening of the highly host-specific and damaging rust pathogen, Puccinia spegazzinii (Pucciniales). It was first released in India and China in 2005/6, although it is not believed to have established. Since then, it has been released successfully in Taiwan, Papua New Guinea (PNG), Fiji and most recently Vanuatu. The rust has established and is spreading rapidly after applying lessons learned from the first releases on the best rust pathotype and release strategy. In PNG, direct monitoring of vegetation change has demonstrated that the rust is having a significant impact on M. micrantha, with no unpredicted non-target impacts. Despite this, the authorities in many countries where mikania is a problem remain cautious about releasing the rust. In Western Samoa, introduction of the rust was not pursued because of a conflict of interest, and the perception that mikania suppresses even worse weeds. For some, ‘pathophobia’ is still a major obstacle. In Indonesia, where insects for weed CBC have been introduced, pathogens will currently not be considered. In other countries such as Bhutan and Myanmar, there are no baseline data on the presence and impact of IAPs and, with no history of CBC, no institutional framework for implementing this approach. Malaysia has a well-developed framework, but capacity needs to be built in the country. Overall, it remains critical to have champions at decision making levels. Hence, even with an effective ‘off-the-shelf’ agent available, implementation of mikania CBC still requires significant inputs tailored to the countries’ specific needs.

Integral Transformations of Volterra Gaussian Processes

We study integral representations of Gaussian processes with a pre-specified law in terms of other Gaussian processes. The dissertation consists of an introduction and of four research articles. In the introduction, we provide an overview about Volterra Gaussian processes in general, and fractional Brownian motion in particular. In the first article, we derive a finite interval integral transformation, which changes fractional Brownian motion with a given Hurst index into fractional Brownian motion with an other Hurst index. Based on this transformation, we construct a prelimit which formally converges to an analogous, infinite interval integral transformation. In the second article, we prove this convergence rigorously and show that the infinite interval transformation is a direct consequence of the finite interval transformation. In the third article, we consider general Volterra Gaussian processes. We derive measure-preserving transformations of these processes and their inherently related bridges. Also, as a related result, we obtain a Fourier-Laguerre series expansion for the first Wiener chaos of a Gaussian martingale. In the fourth article, we derive a class of ergodic transformations of self-similar Volterra Gaussian processes.

Water addition, evaporation and water holding capacity of poultry litter

Litter moisture content has been related to ammonia, dust and odour emissions as well as bird health and welfare. Improved understanding of the water holding properties of poultry litter as well as water additions to litter and evaporation from litter will contribute to improved litter moisture management during the meat chicken grow-out. The purpose of this paper is to demonstrate how management and environmental conditions over the course of a grow-out affect the volume of water A) applied to litter, B) able to be stored in litter, and C) evaporated from litter on a daily basis. The same unit of measurement has been used to enable direct comparison—litres of water per square metre of poultry shed floor area, L/m2, assuming a litter depth of 5 cm. An equation was developed to estimate the amount of water added to litter from bird excretion and drinking spillage, which are sources of regular water application to the litter. Using this equation showed that water applied to litter from these sources changes over the course of a grow-out, and can be as much as 3.2 L/m2/day. Over a 56 day grow-out, the total quantity of water added to the litter was estimated to be 104 L/m2. Litter porosity, water holding capacity and water evaporation rates from litter were measured experimentally. Litter porosity decreased and water holding capacity increased over the course of a grow-out due to manure addition. Water evaporation rates at 25 °C and 50% relative humidity ranged from 0.5 to 10 L/m2/day. Evaporation rates increased with litter moisture content and air speed. Maintaining dry litter at the peak of a grow-out is likely to be challenging because evaporation rates from dry litter may be insufficient to remove the quantity of water added to the litter on a daily basis.

Fate of key food-borne pathogens associated with intensive pig and poultry farming environments

Intensive pig and poultry farming in Australia can be a source of pathogens with implications for food-safety and/or human illness. Seven studies were undertaken with the following objectives: · Assess the types of zoonotic pathogens in waste · Assess the transfer of pathogens during re-use both within the shed and externally in the environment · The potential for movement of pathogens via aerosols In the first and second studies the extent of zoonotic pathogens was evaluated in both piggery effluent and chicken litter and Salmonella and Campylobacter were detected in both wastes. In the third study the dynamics of Salmonella during litter re-use was examined and results showed a trend for lower Salmonella levels and serovar diversity in re-used litter compared to new litter. Thus, re-use within the poultry farming system posed no increased risk. The fourth study addressed the direct risks of pathogens to farm workers due to reuse of piggery effluent within the pig shed. Based on air-borne Escherichia coli (E. coli) levels, re-using effluent did not pose a risk. In the fifth study high levels of Arcobacter spp. were detected in effluent ponds and freshly irrigated soils with potential food-safety risks during the irrigation of food-crops and pasture. The sixth and seventh studies addressed the risks from aerosols from mechanically ventilated sheds. Staphylococci were shown to have potential as markers, with airborne levels gradually dropping and reaching background levels at 400 m distance. Salmonella was detected (at low levels) both inside and outside the shed (at 10 m). Campylobacter was detected only once inside the shed during the 3-year period (at low levels). Results showed there was minimal risk to humans living adjacent to poultry farms This is the first comprehensive analysis studying key food-safety pathogens and potential public health risks associated with intensively farmed pigs and poultry in Australia.

Computational Techniques for Haplotype Inference and for Local Alignment Significance

This thesis which consists of an introduction and four peer-reviewed original publications studies the problems of haplotype inference (haplotyping) and local alignment significance. The problems studied here belong to the broad area of bioinformatics and computational biology. The presented solutions are computationally fast and accurate, which makes them practical in high-throughput sequence data analysis. Haplotype inference is a computational problem where the goal is to estimate haplotypes from a sample of genotypes as accurately as possible. This problem is important as the direct measurement of haplotypes is difficult, whereas the genotypes are easier to quantify. Haplotypes are the key-players when studying for example the genetic causes of diseases. In this thesis, three methods are presented for the haplotype inference problem referred to as HaploParser, HIT, and BACH. HaploParser is based on a combinatorial mosaic model and hierarchical parsing that together mimic recombinations and point-mutations in a biologically plausible way. In this mosaic model, the current population is assumed to be evolved from a small founder population. Thus, the haplotypes of the current population are recombinations of the (implicit) founder haplotypes with some point--mutations. HIT (Haplotype Inference Technique) uses a hidden Markov model for haplotypes and efficient algorithms are presented to learn this model from genotype data. The model structure of HIT is analogous to the mosaic model of HaploParser with founder haplotypes. Therefore, it can be seen as a probabilistic model of recombinations and point-mutations. BACH (Bayesian Context-based Haplotyping) utilizes a context tree weighting algorithm to efficiently sum over all variable-length Markov chains to evaluate the posterior probability of a haplotype configuration. Algorithms are presented that find haplotype configurations with high posterior probability. BACH is the most accurate method presented in this thesis and has comparable performance to the best available software for haplotype inference. Local alignment significance is a computational problem where one is interested in whether the local similarities in two sequences are due to the fact that the sequences are related or just by chance. Similarity of sequences is measured by their best local alignment score and from that, a p-value is computed. This p-value is the probability of picking two sequences from the null model that have as good or better best local alignment score. Local alignment significance is used routinely for example in homology searches. In this thesis, a general framework is sketched that allows one to compute a tight upper bound for the p-value of a local pairwise alignment score. Unlike the previous methods, the presented framework is not affeced by so-called edge-effects and can handle gaps (deletions and insertions) without troublesome sampling and curve fitting.

Algorithms for 13C metabolic flux analysis

The metabolism of an organism consists of a network of biochemical reactions that transform small molecules, or metabolites, into others in order to produce energy and building blocks for essential macromolecules. The goal of metabolic flux analysis is to uncover the rates, or the fluxes, of those biochemical reactions. In a steady state, the sum of the fluxes that produce an internal metabolite is equal to the sum of the fluxes that consume the same molecule. Thus the steady state imposes linear balance constraints to the fluxes. In general, the balance constraints imposed by the steady state are not sufficient to uncover all the fluxes of a metabolic network. The fluxes through cycles and alternative pathways between the same source and target metabolites remain unknown. More information about the fluxes can be obtained from isotopic labelling experiments, where a cell population is fed with labelled nutrients, such as glucose that contains 13C atoms. Labels are then transferred by biochemical reactions to other metabolites. The relative abundances of different labelling patterns in internal metabolites depend on the fluxes of pathways producing them. Thus, the relative abundances of different labelling patterns contain information about the fluxes that cannot be uncovered from the balance constraints derived from the steady state. The field of research that estimates the fluxes utilizing the measured constraints to the relative abundances of different labelling patterns induced by 13C labelled nutrients is called 13C metabolic flux analysis. There exist two approaches of 13C metabolic flux analysis. In the optimization approach, a non-linear optimization task, where candidate fluxes are iteratively generated until they fit to the measured abundances of different labelling patterns, is constructed. In the direct approach, linear balance constraints given by the steady state are augmented with linear constraints derived from the abundances of different labelling patterns of metabolites. Thus, mathematically involved non-linear optimization methods that can get stuck to the local optima can be avoided. On the other hand, the direct approach may require more measurement data than the optimization approach to obtain the same flux information. Furthermore, the optimization framework can easily be applied regardless of the labelling measurement technology and with all network topologies. In this thesis we present a formal computational framework for direct 13C metabolic flux analysis. The aim of our study is to construct as many linear constraints to the fluxes from the 13C labelling measurements using only computational methods that avoid non-linear techniques and are independent from the type of measurement data, the labelling of external nutrients and the topology of the metabolic network. The presented framework is the first representative of the direct approach for 13C metabolic flux analysis that is free from restricting assumptions made about these parameters.In our framework, measurement data is first propagated from the measured metabolites to other metabolites. The propagation is facilitated by the flow analysis of metabolite fragments in the network. Then new linear constraints to the fluxes are derived from the propagated data by applying the techniques of linear algebra.Based on the results of the fragment flow analysis, we also present an experiment planning method that selects sets of metabolites whose relative abundances of different labelling patterns are most useful for 13C metabolic flux analysis. Furthermore, we give computational tools to process raw 13C labelling data produced by tandem mass spectrometry to a form suitable for 13C metabolic flux analysis.