Bayesian model comparison in cosmology with Population Monte Carlo

We use Bayesian model selection techniques to test extensions of the standard flat LambdaCDM paradigm. Dark-energy and curvature scenarios, and primordial perturbation models are considered. To that end, we calculate the Bayesian evidence in favour of each model using Population Monte Carlo (PMC), a new adaptive sampling technique which was recently applied in a cosmological context. The Bayesian evidence is immediately available from the PMC sample used for parameter estimation without further computational effort, and it comes with an associated error evaluation. Besides, it provides an unbiased estimator of the evidence after any fixed number of iterations and it is naturally parallelizable, in contrast with MCMC and nested sampling methods. By comparison with analytical predictions for simulated data, we show that our results obtained with PMC are reliable and robust. The variability in the evidence evaluation and the stability for various cases are estimated both from simulations and from data. For the cases we consider, the log-evidence is calculated with a precision of better than 0.08. Using a combined set of recent CMB, SNIa and BAO data, we find inconclusive evidence between flat LambdaCDM and simple dark-energy models. A curved Universe is moderately to strongly disfavoured with respect to a flat cosmology. Using physically well-motivated priors within the slow-roll approximation of inflation, we find a weak preference for a running spectral index. A Harrison-Zel'dovich spectrum is weakly disfavoured. With the current data, tensor modes are not detected; the large prior volume on the tensor-to-scalar ratio r results in moderate evidence in favour of r=0.

Fluoroquinolone-resistant extraintestinal pathogenic Escherichia coli, including O25b-ST131, isolated from faeces of hospitalized dogs in an Australian veterinary referral centre

To determine rates of carriage of fluoroquinolone-resistant Escherichia coli and extraintestinal pathogenic E. coli (ExPEC) among dogs in a specialist referral hospital and to examine the population structure of the isolates. Fluoroquinolone-resistant faecal E. coli isolates (n232, from 23 of 123 dogs) recovered from hospitalized dogs in a veterinary referral centre in Sydney, Australia, over 140 days in 2009 were characterized by phylogenetic grouping, virulence genotyping and random amplified polymorphic DNA (RAPD) analysis. The RAPD dendrogram for representative isolates showed one group B2-associated cluster and three group D-associated clusters; each contained isolates with closely related ExPEC-associated virulence profiles. All group B2 faecal isolates represented the O25b-ST131 clonal group and were closely related to recent canine extraintestinal ST131 clinical isolates from the east coast of Australia by RAPD analysis. Hospitalized dogs may carry fluoroquinolone-resistant ExPEC in their faeces, including those representing O25b-ST131.

Impacts of level of utilisation by grazing on an Astrebla (Mitchell grass) grassland in north-western Queensland between 1984 and 2010. 1. Herbage mass and population dynamics of Astrebla spp

Managing large variations in herbage production, resulting from highly variable seasonal rainfall, provides a major challenge for the sustainable management of Astrebla (Mitchell grass) grasslands in Australia. A grazing study with sheep was conducted between 1984 and 2010 on an Astrebla grassland in northern Queensland to describe the effects of a range of levels of utilisation of the herbage at the end of the summer growing season (April–May in northern Australia) on the sustainability of these grasslands. In unreplicated paddocks, sheep numbers were adjusted annually to achieve 0, 10, 20, 30, 50 and 80% utilisation of the herbage mass at the end of the summer over the ensuing 12 months. Higher levels of utilisation reduced both total and Astrebla spp. herbage mass because of the effects of higher utilisation on Astrebla spp. and this effect was accentuated by drought. The tussock density of Astrebla spp. varied widely among years but with few treatment differences until 2005 when density was reduced at the 50% level of utilisation. A major change in density resulted from a large recruitment of Astrebla spp. in 1989 that influenced its density for the remainder of the study. Basal area of the tussocks fluctuated among years, with increases due to rainfall and decreases during droughts. Seasonal rainfall was more influential than level of utilisation in changes to the basal area of perennial grasses. Drought resulted in the death of Astrebla spp. tussocks and this effect was accentuated at higher levels of utilisation. A series of three grazing exclosures were used to examine the recovery of the density and basal area of Astrebla spp. after it had been reduced by 80% utilisation over the preceding 9 years. This recovery study indicated that, although grazing exclusion was useful in the recovery of Astrebla spp., above-average rainfall was the major factor driving increases in the basal area of perennial grasses. Spring values of the Southern Oscillation Index and associated rainfall probabilities were considered to have potential for understanding the dynamics of Astrebla spp. It was concluded that Astrebla grassland remained sustainable after 26 years when grazed at up to 30% utilisation, while, at 50% utilisation, they became unsustainable after 20 years. Results from this study emphasised the need to maintain the population of Astrebla spp. tussocks.

The (Non)Effects of Lethal Population Control on the Diet of Australian Dingoes

Top-predators contribute to ecosystem resilience, yet individuals or populations are often subject to lethal control to protect livestock, managed game or humans from predation. Such management actions sometimes attract concern that lethal control might affect top-predator function in ways ultimately detrimental to biodiversity conservation. The primary function of a predator is predation, which is often investigated by assessing their diet. We therefore use data on prey remains found in 4,298 Australian dingo scats systematically collected from three arid sites over a four year period to experimentally assess the effects of repeated broad-scale poison-baiting programs on dingo diet. Indices of dingo dietary diversity and similarity were either identical or near-identical in baited and adjacent unbaited treatment areas in each case, demonstrating no control-induced change to dingo diets. Associated studies on dingoes' movement behaviour and interactions with sympatric mesopredators were similarly unaffected by poison-baiting. These results indicate that mid-sized top-predators with flexible and generalist diets (such as dingoes) may be resilient to ongoing and moderate levels of population control without substantial alteration of their diets and other related aspects of their ecological function.

Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes

Despite international protection of white sharks Carcharodon carcharias, important conservation parameters such as abundance, population structure and genetic diversity are largely unknown. The tissue of 97 predominately juvenile white sharks sampled from spatially distant eastern and southwestern Australian coastlines was sequenced for the mitochondrial DNA (mtDNA) control region and genotyped with 6 nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (F-ST = 0.142, p < 0.0001), implying female reproductive philopatry. This concurs with recent satellite and acoustic tracking findings which suggest the sustained presence of discrete east coast nursery areas. Furthermore, population subdivision was found between the same regions with biparentally inherited micro satellite markers (F-ST = 0.009, p < 0.05), suggesting that males may also exhibit some degree of reproductive philopatry; 5 sharks captured along the east coast had mtDNA haplotypes that resembled western Indian Ocean sharks more closely than Australian/New Zealand sharks, suggesting that transoceanic dispersal, or migration resulting in breeding, may occur sporadically. Our most robust estimate of contemporary genetic effective population size was low and close to thresholds at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least 1, possibly 2, orders of magnitude below our historical effective size estimates. Population decline could expose these genetically isolated populations to detrimental genetic effects. Regional Australian white shark conservation management units should be implemented until genetic population structure, size and diversity can be investigated in more detail.

Proportion of Jews to general population of the European provinces of Russia

Digital image

Population density of the European provinces of Russia

Digital image

Population dynamics of Thaumastocoris peregrinus in Eucalyptus plantations of South Africa

Thaumastocoris peregrinus is a sap-sucking insect that infests non-native Eucalyptus plantations in Africa, New Zealand, South America and parts of Southern Europe, in addition to street trees in parts of its native range of Australia. In South Africa, pronounced fluctuations in the population densities have been observed. To characterise spatiotemporal variability in T. peregrinus abundance and the factors that might influence it, we monitored adult population densities at six sites in the main eucalypt growing regions of South Africa. At each site, twenty yellow sticky traps were monitored weekly for 30 months, together with climatic data. We also characterised the influence of temperature on growth and survival experimentally and used this to model how temperature may influence population dynamics. T. peregrinus was present throughout the year at all sites, with annual site-specific peaks in abundance. Peaks occurred during autumn (February–April) for the Pretoria site, summer (November–January) for the Zululand site and spring (August–October) for the Tzaneen, Sabie and Piet Retief monitoring sites. Temperature (both experimental and field-collected), humidity and rainfall were mostly weakly, or not at all, associated with population fluctuations. It is clear that a complex interaction of these and other factors (e.g. host quality) influence population fluctuations in an annual, site specific cycle. The results obtained not only provide insights into the biology of T. peregrinus, but will also be important for future planning of monitoring and control programs using semiochemicals, chemical insecticides or biological control agents.

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.

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.

Human-associated fluoroquinolone-resistant Escherichia coli clonal lineages, including ST354, isolated from canine feces and extraintestinal infections in Australia

Phylogenetic group D extraintestinal pathogenic Escherichia coli (ExPEC), including O15:K52:H1 and clonal group A, have spread globally and become fluoroquinolone-resistant. Here we investigated the role of canine feces as a reservoir of these (and other) human-associated ExPEC and their potential as canine pathogens. We characterized and compared fluoroquinolone-resistant E. coli isolates originally identified as phylogenetic group D from either the feces of hospitalized dogs (n = 67; 14 dogs) or extraintestinal infections (n = 53; 33 dogs). Isolates underwent phylogenetic grouping, random amplified polymorphic DNA (RAPD) analysis, virulence genotyping, resistance genotyping, human-associated ExPEC O-typing, and multi-locus sequence typing. Five of seven human-associated sequence types (STs) exhibited ExPEC-associated O-types, and appeared in separate RAPD clusters. The largest subgroup (16 fecal, 26 clinical isolates) were ST354 (phylogroup F) isolates. ST420 (phylogroup B2); O1-ST38, O15:K52:H1-ST393, and O15:K1-ST130 (phylogroup D); and O7-ST457, and O1-ST648 (phylogroup F) were also identified. Three ST-specific RAPD sub-clusters (ST354, ST393, and ST457) contained closely related isolates from both fecal or clinical sources. Genes encoding CTX-M and AmpC β-lactamases were identified in isolates from five STs. Major human-associated fluoroquinolone-resistant ± extended-spectrum cephalosporin-resistant ExPEC of public health importance may be carried in dog feces and cause extraintestinal infections in some dogs.

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.