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.

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.

Invisible opportunities and irreversible trends

Small firms are always vulnerable to complex technological change that may render their existing business model obsolete. This paper emphasises the need to understand how the Internet's ubiquitous World Wide Web is impacting on their operating environments. Consideration of evolutionary theory and the absorptive capacity construct provides the foundation for discussion of how learning and discovery take place within individuals, firms and the environments that interact with. Small firms, we argue, face difficulties identifying what routines and competencies are best aligned with the seemingly invisible dominant designs that support pursuit of new enterprise in web-impacted environments. We argue that such difficulties largely relate to an inability to acquire external knowledge and the subsequent reliance on existing internal selection processes that may reinforce the known, at the expense of the unknown. The paper concludes with consideration as to how managers can overcome the expected difficulties through the development of internal routines that support the continual search, evaluation and acquisition of specific external knowledge.

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.

Frost trends and their estimated impact on yield in the Australian wheatbelt

Radiant spring frosts occurring during reproductive developmental stages can result in catastrophic yield loss for wheat producers. To better understand the spatial and temporal variability of frost, the occurrence and impact of frost events on rain-fed wheat production was estimated across the Australian wheatbelt for 1957–2013 using a 0.05 ° gridded weather data set. Simulated yield outcomes at 60 key locations were compared with those for virtual genotypes with different levels of frost tolerance. Over the last six decades, more frost events, later last frost day, and a significant increase in frost impact on yield were found in certain regions of the Australian wheatbelt, in particular in the South-East and West. Increasing trends in frost-related yield losses were simulated in regions where no significant trend of frost occurrence was observed, due to higher mean temperatures accelerating crop development and causing sensitive post-heading stages to occur earlier, during the frost risk period. Simulations indicated that with frost-tolerant lines the mean national yield could be improved by up to 20 through (i) reduced frost damage (~10 improvement) and (ii) the ability to use earlier sowing dates (adding a further 10 improvement). In the simulations, genotypes with an improved frost tolerance to temperatures 1 °C lower than the current 0 °C reference provided substantial benefit in most cropping regions, while greater tolerance (to 3 °C lower temperatures) brought further benefits in the East. The results indicate that breeding for improved reproductive frost tolerance should remain a priority for the Australian wheat industry, despite warming climates.

Fauna community trends during early restoration of alluvial open forest/woodland ecosystems on former agricultural land

Vertebrate fauna was studied over 10 years following revegetation of a Eucalyptus tereticornis ecosystem on former agricultural land. We compared four vegetation types: remnant forest, plantings of a mix of native tree species on cleared land, natural regeneration of partially cleared land after livestock removal, and cleared pasture land with scattered paddock trees managed for livestock production. Pasture differed significantly from remnant in both bird and nonbird fauna. Although 10 years of ecosystem restoration is relatively short term in the restoration process, in this time bird assemblages in plantings and natural regeneration had diverged significantly from pasture, but still differed significantly from remnant. After 10 years, 70 and 66% of the total vertebrate species found in remnant had been recorded in plantings and natural regeneration, respectively. Although the fauna assemblages within plantings and natural regeneration were tracking toward those of remnant, significant differences in fauna between plantings and natural regeneration indicated community development along different restoration pathways. Because natural regeneration contained more mature trees (dbh > 30 cm), native shrub species, and coarse woody debris than plantings from the beginning of the study, these features possibly encouraged different fauna to the revegetation areas from the outset. The ability of plantings and natural regeneration to transition to the remnant state will be governed by a number of factors that were significant in the analyses, including shrub cover, herbaceous biomass, tree hollows, time since fire, and landscape condition. Both active and passive restoration produced significant change from the cleared state in the short term.

Integrating Rapid Phenotyping and Speed Breeding to Improve Stay-Green and Root Adaptation of Wheat in Changing, Water-Limited, Australian Environments

Temperatures have increased and in-crop rainfall decreased over recent decades in many parts of the Australian wheat cropping region. With these trends set to continue or intensify, improving crop adaptation in the face of climate change is particularly urgent in this, already drought-prone, cropping region. Importantly, improved performance under water-limitation must be achieved while retaining yield potential during more favourable seasons. A multi-trait-based approach to improve wheat yield and yield stability in the face of water-limitation and heat has been instigated in northern Australia using novel phenotyping techniques and a nested association mapping (NAM) approach. An innovative laboratory technique allows rapid root trait screening of hundreds of lines. Using soil grown seedlings, the method offers significant advantages over many other lab-based techniques. Another recently developed method allows novel stay-green traits to be quantified objectively for hundreds of genotypes in standard field trial plots. Field trials in multiple locations and seasons allow evaluation of targeted trait values and identification of superior germplasm. Traits, including yield and yield components are measured for hundreds of NAM lines in rain fed environments under various levels of water-limitation. To rapidly generate lines of interest, the University of Queensland “speed breeding” method is being employed, allowing up to 7 plant generations per annum. A NAM population of over 1000 wheat recombinant inbred lines has been progressed to the F5 generation within 18 months. Genotyping the NAM lines with the genome-wide DArTseq molecular marker system provides up to 40,000 markers. They are now being used for association mapping to validate QTL previously identified in bi-parental populations and to identify novel QTL for stay-green and root traits. We believe that combining the latest techniques in physiology, phenotyping, genetics and breeding will increase genetic progress toward improved adaptation to water-limited environments.

The diet-related health of the Australian population and potential opportunities for the food retail sector

This report describes the diet-related health of the Australian population and identifies potential opportunities for the food retail sector.

Early post-operative mortality after major lower limb amputation: A systematic review of population and regional based studies

OBJECTIVE: Lower limb amputation is often associated with a high risk of early post-operative mortality. Mortality rates are also increasingly being put forward as a possible benchmark for surgical performance. The primary aim of this systematic review is to investigate early post-operative mortality following a major lower limb amputation in population/regional based studies, and reported factors that might influence these mortality outcomes. METHODS: Embase, PubMed, Cinahl and Psycinfo were searched for publications in any language on 30 day or in hospital mortality after major lower limb amputation in population/regional based studies. PRISMA guidelines were followed. A self developed checklist was used to assess quality and susceptibility to bias. Summary data were extracted for the percentage of the population who died; pooling of quantitative results was not possible because of methodological differences between studies. RESULTS: Of the 9,082 publications identified, results were included from 21. The percentage of the population undergoing amputation who died within 30 days ranged from 7% to 22%, the in hospital equivalent was 4-20%. Transfemoral amputation and older age were found to have a higher proportion of early post-operative mortality, compared with transtibial and younger age, respectively. Other patient factors or surgical treatment choices related to increased early post-operative mortality varied between studies. CONCLUSIONS: Early post-operative mortality rates vary from 4% to 22%. There are very limited data presented for patient related factors (age, comorbidities) that influence mortality. Even less is known about factors related to surgical treatment choices, being limited to amputation level. More information is needed to allow comparison across studies or for any benchmarking of acceptable mortality rates. Agreement is needed on key factors to be reported.

The Population Consequences of Sex and Conflict

Natural selection generally operates at the level of the individual, or more specifically at the level of the gene. As a result, individual selection does not always favour traits which benefit the population or species as a whole. The spread of an individual gene may even act to the detriment of the organism in which it finds. Thus selection at the level of the individual can affect processes at the level of the organism, group or even at the level of the species. As most behaviours ultimately affect births, deaths and the distribution of individuals, it seems inevitable that behavioural decisions will have an impact on population dynamics and population densities. Behavioural decisions can often involve costs through allocation of energy into behavioural strategies, such as the investment into armaments involved in fighting over resources or increased mortality due to injury or increased predation risk. Similarly, behaviour may act o to benefit the population, in terms of higher survival and increased fecundity. Examples include increased investment through parental care, choosing a mate based on the nuptial gifts they may supply and choosing territories in the face of competition. Investigating the impact of behaviour on population ecology may seem like a trivial task, but it is likely to have important consequences at different levels. For example, antagonistic behaviour may occasionally become so extreme that it increases the risk of extinction, and such extinction risk may have important implications for conservation. As a corollary, any such behaviour may also act as a macroevolutionary force, weeding out populations with traits which, whilst beneficial to the individuals in the short term, ultimately result in population extinction. In this thesis, I examine how behaviours, specifically conflict and competition over a resource and aspects of behaviour involved in sexual selection, can affect population densities, and what the implications are for the evolution and ecology of the populations in question. It is found that both behaviours related to individual conflict and mating strategies can have an effect at the level of the population, but that various factors, such as a feedback between selection and population densities or macroevolution caused by species extinctions, may act to limit the intensity of conflicts that we observe in nature.

Phylogeography and Population Genetics of the Moor frog (Rana arvalis Nilsson) in Northern Europe

The ongoing climate change along with increasing levels of pollutants, diseases, habitat loss and fragmentation constitute global threats to the persistence of many populations, species and ecosystems. However, for the long-term persistence of local populations, one of the biggest threats is the intrinsic loss of genetic variation. In order to adapt to changes in the environment, organisms must have a sufficient supply of heritable variation in traits important for their fitness. With a loss of genetic variation, the risk of extinction will increase. For conservational practices, one should therefore understand the processes that shape the genetic population structure and also the broader (historical) phylogenetic patterning of the species in focus. In this thesis, microsatellite markers were applied to study genetic diversity and population differentiation of the protected moor frog (Rana arvalis) in Fennoscandia from both historical (evolutionary) and applied (conservation) perspectives. The results demonstrate that R. arvalis populations are highly structured over rather short geographic distances. Moreover, the results suggest that R. arvalis recolonized Fennoscandia from two directions after the last ice age. This has had implications for the genetic structuring and population differentiation, especially in the northernmost parts where the two lineages have met. Compared to more southern populations, the genetic variation decreases and the interpopulation differentiation increases dramatically towards north. This could be an outcome of serial population bottlenecking along the recolonization route. Also, current isolation and small population sizes increase the effect of drift, thus reinforcing the observed pattern. The same pattern can also be seen in island populations. However, though R. arvalis on the island of Gotland has lost most of its neutral genetic variability, our results indicate that the levels of additive genetic variation have remained high. This conforms to the conjecture that though neutral markers are widely used in conservation purposes, they may be quite uninformative about the levels of genetic variation in ecologically important traits. Finally, the evolutionary impact of the typical amphibian mating behaviour on genetic diversity was investigated. Given the short time available for larval development, it is important that mating takes place as early as possible. The genetic data and earlier capture-recapture data suggest that R. arvalis gather at mating grounds they are familiar with. However, by forming leks in random to relatedness, and having multiple paternities in single clutches, the risk of inbreeding may be minimized in this otherwise highly philopatric species.