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

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

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.

Reassessing spatial and temporal dynamics of kangaroo populations

Long-running datasets from aerial surveys of kangaroos (Macropus giganteus, Macropus [uliginosus, Macropus robustus and Macropus rufus) across Queensland, New South Wales and South Australia have been analysed, seeking better predictors of rates of increase which would allow aerial surveys to be undertaken less frequently than annually. Early models of changes in kangaroo numbers in response to rainfall had shown great promise, but much variability. We used normalised difference vegetation index (NDVI) instead, reasoning that changes in pasture condition would provide a better predictor than rainfall. However, except at a fine scale, NDVI proved no better; although two linked periods of rainfall proved useful predictors of rates of increase, this was only in some areas for some species. The good correlations reported in earlier studies were a consequence of data dominated by large droughtinduced adult mortality, whereas over a longer time frame and where changes between years are less dramatic, juvenile survival has the strongest influence on dynamics. Further, harvesting, density dependence and competition with domestic stock are additional and important influences and it is now clear that kangaroo movement has a greater influence on population dynamics than had been assumed. Accordingly, previous conclusions about kangaroo populations as simple systems driven by rainfall need to be reassessed. Examination of this large dataset has permitted descriptions of shifts in distribution of three species across eastern Australia, changes in dispersion in response to rainfall, and an evaluation of using harvest statistics as an index of density and harvest rate. These results have been combined into a risk assessment and decision theory framework to identify optimal monitoring strategies.

Scallop fishery in south east Australia

Establishing fine-scale industry based spatial management and harvest strategies for the Commercial Scallop fishery in south east Australia.

White shark (Carcharodon carcharias) microsatellite genotypes (six loci) from regions along Australia's southern and eastern coastline.

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 six nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (FST = 0.142, p < 0.001), implying female natal philopatry. This concords 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 microsatellite markers (FST = 0.009, p <0.05), suggesting that males may also exhibit some degree of reproductive philopatry. Five 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 below the threshold at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least one, possibly two orders of magnitude below our historical effective size estimates. Further 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. Reference: Blower, D. C., Pandolfi, J. M., Gomez-Cabrera, M. del C., Bruce, B. D. & Ovenden, J. R. (In press - April 2012). Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes. Marine Ecology Progress Series.

The short-term effects of a routine poisoning campaign on the movements and detectability of a social top-predator

Top-predators can be important components of resilient ecosystems, but they are still controlled in many places to mitigate a variety of economic, environmental and/or social impacts. Lethal control is often achieved through the broad-scale application of poisoned baits. Understanding the direct and indirect effects of such lethal control on subsequent movements and behaviour of survivors is an important pre-requisite for interpreting the efficacy and ecological outcomes of top-predator control. In this study, we use GPS tracking collars to investigate the fine-scale and short-term movements of dingoes (Canis lupus dingo and other wild dogs) in response to a routine poison-baiting program as an example of how a common, social top-predator can respond (behaviourally) to moderate levels of population reduction. We found no consistent control-induced differences in home range size or location, daily distance travelled, speed of travel, temporal activity patterns or road/trail usage for the seven surviving dingoes we monitored immediately before and after a typical lethal control event. These data suggest that the spatial behaviour of surviving dingoes was not altered in ways likely to affect their detectability, and if control-induced changes in dingoes' ecological function did occur, these may not be related to altered spatial behaviour or movement patterns.

Population Expansion and Genetic Structure in Carcharhinus brevipinna in the Southern Indo-Pacific

Background:Quantifying genetic diversity and metapopulation structure provides insights into the evolutionary history of a species and helps develop appropriate management strategies. We provide the first assessment of genetic structure in spinner sharks (Carcharhinus brevipinna), a large cosmopolitan carcharhinid, sampled from eastern and northern Australia and South Africa. Methods and Findings:Sequencing of the mitochondrial DNA NADH dehydrogenase subunit 4 gene for 430 individuals revealed 37 haplotypes and moderately high haplotype diversity (h = 0.6770 ±0.025). While two metrics of genetic divergence (ΦST and FST) revealed somewhat different results, subdivision was detected between South Africa and all Australian locations (pairwise ΦST, range 0.02717–0.03508, p values ≤ 0.0013; pairwise FST South Africa vs New South Wales = 0.04056, p = 0.0008). Evidence for fine-scale genetic structuring was also detected along Australia’s east coast (pairwise ΦST = 0.01328, p < 0.015), and between south-eastern and northern locations (pairwise ΦST = 0.00669, p < 0.04).Conclusions: The Indian Ocean represents a robust barrier to contemporary gene flow in C. brevipinna between Australia and South Africa. Gene flow also appears restricted along a continuous continental margin in this species, with data tentatively suggesting the delineation of two management units within Australian waters. Further sampling, however, is required for a more robust evaluation of the latter finding. Evidence indicates that all sampled populations were shaped by a substantial demographic expansion event, with the resultant high genetic diversity being cause for optimism when considering conservation of this commercially-targeted species in the southern Indo-Pacific.

Understanding the barriers to the implementation of precision agriculture in the central region

There is an increasing requirement for more astute land resource management through efficiencies in agricultural inputs in a sugar cane production system. A precision agriculture (PA) approach can provide a pathway for a sustainable sugarcane production system. One of the impediments to the adoption of PA practices is access to paddock-scale mapping layers displaying variability in soil properties, crop growth and surface drainage. Variable rate application (VRA) of nutrients is an important component of PA. However, agronomic expertise within PA systems has fallen well behind significant advances in PA technologies. Generally, advisers in the sugar industry have a poor comprehension of the complex interaction of variables that contribute to within-paddock variations in crop growth. This is regarded as a significant impediment to the progression of PA in sugarcane and is one of the reasons for the poor adoption of VRA of nutrients in a PA approach to improved sugar cane production. This project therefore has established a number of key objectives which will contribute to the adoption of PA and the staged progression of VRA supported by relevant and practical agronomic expertise. These objectives include provision of base soils attribute mapping that can be determined using Veris 3100 Electrical Conductivity (EC) and digital elevation datasets using GPS mapping technology for a large sector of the central cane growing region using analysis of archived satellite imagery to determine the location and stability of yield patterns over time and in varying seasonal conditions on selected project study sites. They also include the stablishment of experiments to determine appropriate VRA nitrogen rates on various soil types subjected to extended anaerobic conditions, and the establishment of trials to determine nitrogen rates applicable to a declining yield potential associated with the aging of ratoons in the crop cycle. Preliminary analysis of archived yield estimation data indicates that yield patterns remain relatively stable overtime. Results also indicate the where there is considerable variability in EC values there is also significant variation in yield.

Guardian dogs protect sheep by guarding sheep, not by establishing territories and excluding predators

Guardian animals have been a common non-lethal method for reducing predator impacts on livestock for centuries in Europe. But elsewhere, livestock producers sometimes doubt whether such methods work or are compatible with modern livestock husbandry practices in extensive grazing systems. In this study we evaluate the hypothesis that guardian dogs primarily ‘work’ by establishing and defending territories from which canid predators are excluded. Eight maremmas and six free-ranging wild dogs of different sexes were fitted with GPS collars and monitored for 7 months on a large sheep property in north Queensland, Australia. Wild dog incursions into the territories of adjacent wild dogs and maremmas were recorded. Wild dog territories never overlapped and their home ranges infrequently overlapped. In contrast, 713 hourly locations from 120 wild dog incursions into maremma territories were recorded, mostly from three wild dogs. These three wild dogs spent a mean of 2.5–5.9 h inside maremma territories during incursions. At this location, maremmas worked by guarding sheep and prohibiting fine-scale interaction between wild dogs and sheep, not by establishing a territory respected by wild dogs. We conclude that shepherding behaviour and boisterous vocalisations of guardian dogs combined with the flocking behaviour of sheep circumvents attacks on sheep but does not prevent nor discourage wild dogs from foraging in close proximity. Certain husbandry practices and the behaviour of sheep at parturition may incur greater predation risk.

Fine mapping of qGW1, a major QTL for grain weight in sorghum

Key message We detected seven QTLs for 100-grain weight in sorghum using an F 2 population, and delimited qGW1 to a 101-kb region on the short arm of chromosome 1, which contained 13 putative genes. Abstract Sorghum is one of the most important cereal crops. Breeding high-yielding sorghum varieties will have a profound impact on global food security. Grain weight is an important component of grain yield. It is a quantitative trait controlled by multiple quantitative trait loci (QTLs); however, the genetic basis of grain weight in sorghum is not well understood. In the present study, using an F2 population derived from a cross between the grain sorghum variety SA2313 (Sorghum bicolor) and the Sudan-grass variety Hiro-1 (S. bicolor), we detected seven QTLs for 100-grain weight. One of them, qGW1, was detected consistently over 2 years and contributed between 20 and 40 % of the phenotypic variation across multiple genetic backgrounds. Using extreme recombinants from a fine-mapping F3 population, we delimited qGW1 to a 101-kb region on the short arm of chromosome 1, containing 13 predicted gene models, one of which was found to be under purifying selection during domestication. However, none of the grain size candidate genes shared sequence similarity with previously cloned grain weight-related genes from rice. This study will facilitate isolation of the gene underlying qGW1 and advance our understanding of the regulatory mechanisms of grain weight. SSR markers linked to the qGW1 locus can be used for improving sorghum grain yield through marker-assisted selection.

Fine Mapping of qDor7, a Major QTL Affecting Seed Dormancy in Sorghum (Sorghum bicolor (L.) Moench)

Seed dormancy is a key domestication trait for major crops, which is acquired in long-term systems development processes and enables the survival of plants in adverse natural conditions. It is a complex trait under polygenic control and is affected by endogenous and environmental factors. In the present study, a major seed dormancy QTL in sorghum (Sorghum bicolor (L.) Moench), qDor7, detected previously, was fine mapped using a large, multi-generational population. The qDor7 locus was delimited to a 96-kb region which contains 16 predicted gene models. These results lay a solid foundation for cloning qDor7. In addition, the functional markers tightly linked to the seed dormancy QTL may be used in marker-assisted selection for seed dormancy in sorghum.

Predictive modelling to aid the regional-scale management of a vertebrate pest

Extensive resources are allocated to managing vertebrate pests, yet spatial understanding of pest threats, and how they respond to management, is limited at the regional scale where much decision-making is undertaken. We provide regional-scale spatial models and management guidance for European rabbits (Oryctolagus cuniculus) in a 260,791 km(2) region in Australia by determining habitat suitability, habitat susceptibility and the effects of the primary rabbit management options (barrier fence, shooting and baiting and warren ripping) or changing predation or disease control levels. A participatory modelling approach was used to develop a Bayesian network which captured the main drivers of suitability and spread, which in turn was linked spatially to develop high resolution risk maps. Policy-makers, rabbit managers and technical experts were responsible for defining the questions the model needed to address, and for subsequently developing and parameterising the model. Habitat suitability was determined by conditions required for warren-building and by above-ground requirements, such as food and harbour, and habitat susceptibility by the distance from current distributions, habitat suitability, and the costs of traversing habitats of different quality. At least one-third of the region had a high probability of being highly suitable (support high rabbit densities), with the model supported by validation. Habitat susceptibility was largely restricted by the current known rabbit distribution. Warren ripping was the most effective control option as warrens were considered essential for rabbit persistence. The anticipated increase in disease resistance was predicted to increase the probability of moderately suitable habitat becoming highly suitable, but not increase the at-risk area. We demonstrate that it is possible to build spatial models to guide regional-level management of vertebrate pests which use the best available knowledge and capture fine spatial-scale processes.