Managing glyphosate resistance in Australian cotton farming: modelling shows how to delay evolution and maintain long-term population control

Glyphosate resistance is a rapidly developing threat to profitability in Australian cotton farming. Resistance causes an immediate reduction in the effectiveness of in-crop weed control in glyphosate-resistant transgenic cotton and summer fallows. Although strategies for delaying glyphosate resistance and those for managing resistant populations are qualitatively similar, the longer resistance can be delayed, the longer cotton growers will have choice over which tactics to apply and when to apply them. Effective strategies to avoid, delay, and manage resistance are thus of substantial value. We used a model of glyphosate resistance dynamics to perform simulations of resistance evolution in Sonchus oleraceus (common sowthistle) and Echinochloa colona (awnless barnyard grass) under a range of resistance prevention, delaying, and management strategies. From these simulations, we identified several elements that could contribute to effective glyphosate resistance prevention and management strategies. (i) Controlling glyphosate survivors is the most robust approach to delaying or preventing resistance. High-efficacy, high-frequency survivor control almost doubled the useful lifespan of glyphosate from 13 to 25 years even with glyphosate alone used in summer fallows. (ii) Two non-glyphosate tactics in-crop plus two in-summer fallows is the minimum intervention required for long-term delays in resistance evolution. (iii) Pre-emergence herbicides are important, but should be backed up with non-glyphosate knockdowns and strategic tillage; replacing a late-season, pre-emergence herbicide with inter-row tillage was predicted to delay glyphosate resistance by 4 years in awnless barnyard grass. (iv) Weed species' ecological characteristics, particularly seed bank dynamics, have an impact on the effectiveness of resistance strategies; S. oleraceus, because of its propensity to emerge year-round, was less exposed to selection with glyphosate than E. colona, resulting in an extra 5 years of glyphosate usefulness (18 v. 13 years) even in the most rapid cases of resistance evolution. Delaying tactics are thus available that can provide some or many years of continued glyphosate efficacy. If glyphosate-resistant cotton cropping is to remain profitable in Australian farming systems in the long-term, however, growers must adapt to the probability that they will have to deal with summer weeds that are no longer susceptible to glyphosate. Robust resistance management systems will need to include a diversity of weed control options, used appropriately.

Effect of fumigation temperature on the efficacy of phosphine against strongly resistant psocids liposcelis bostrychophila (Pscodoptera: Liposcelididae)

Species of Liposcelis psocids have emerged as major pests of stored grain in Australia in recent years. Several populations have been detected with high resistance to phosphine, the major chemical treatment. Highest resistance has been detected in the cosmopolitan species Liposcelis bostrychophila. As part of a national resistance management strategy to maintain the viability of phosphine, we are developing minimum effective dosage regimes (concentration x time) required to control all life stages of resistant L. bostrychophila at a range of grain temperatures. Four concentrations of phosphine, 0.1, 0.17, 0.3 aid 1 mg/L, were evaluated for their effectiveness against strongly resistant L. bostrychophila at a series of fumigation temperatures: 20, 25, 30 and 35°C. Results were recorded as the least number of days taken to achieve population extinction. We found that, at any fixed concentration of phosphine, time to population extinction decreased as fumigation temperature increased from 20 to 30°C. For example, at 0.1 mg/L, it took more than 14 days at 20°C to completely control these insects, whereas at 30°C it took only seven days. Increase in fumigation temperature from 25OC to 30°C dramatically reduced the exposure period needed to achieve population extinction of resistant psocids. For example, a dose of 0.17 mg/L over six days at 30°C completely controlled strongly resistant L. bostrychophila populations that can survive at 1 mg/L and 25°C over the same exposure period. Findings from our study will be used to formulate recommendations for registered dosage rates and fumigation periods for use in Australia.

Pronounced genetic population structure in a potentially vagile fish species (Pristipomoides multidens, Teleostei; Perciformes; Lutjanidae) from the East Indies triangle

The East Indies triangle, bordered by the Phillipines, Malay Peninsula and New Guinea, has a high level of tropical marine species biodiversity. Pristipomoides multidens is a large, long-lived, fecund snapper species that is distributed throughout the East Indies and Indo-Pacific. Samples were analysed from central and eastern Indonesia and northern Australia to test for genetic discontinuities in population structure. Fish (n = 377) were collected from the Indonesian islands of Bali, Sumbawa, Flores, West Timor, Tanimbar and Tual along with 131 fish from two northern Australian locations (Arafura and Timor Seas) from a previous study. Genetic variation in the control region of the mitochondrial genome was assayed using restriction fragment length polymorphism and direct sequencing. Haplotype diversity was high (0.67-0.82), as was intraspecific sequence divergence (range 0-5.8%). FST between pairs of populations ranged from 0 to 0.2753. Genetic subdivision was apparent on a small spatial scale; FST was 0.16 over 191 km (Bali/Sumbawa) and 0.17 over 491 km (Bali/Flores). Constraints to dispersal that contribute to, and maintain, the observed degree of genetic subdivision are experienced presumably by all life history stages of this tropical marine finfish. The constraints may include (1) little or no movement of eggs or larvae, (2) little or no home range or migratory movement of adults and (3) loss of larval cohorts due to transport of larvae away from suitable habitat by prevailing currents

Low-density koala (Phascolarctos cinereus) populations in the mulgalands of south-west Queensland. IV. Abundance and conservation status

Urban encroachment on dense, coastal koala populations has ensured that their management has received increasing government and public attention. The recently developed National Koala Conservation Strategy calls for maintenance of viable populations in the wild. Yet the success of this, and other, conservation initiatives is hampered by lack of reliable and generally accepted national and regional population estimates. In this paper we address this problem in a potentially large, but poorly studied, regional population in the State that is likely to have the largest wild populations. We draw on findings from previous reports in this series and apply the faecal standing-crop method (FSCM) to derive a regional estimate of more than 59 000 individuals. Validation trials in riverine communities showed that estimates of animal density obtained from the FSCM and direct observation were in close agreement. Bootstrapping and Monte Carlo simulations were used to obtain variance estimates for our population estimates in different vegetation associations across the region. The most favoured habitat was riverine vegetation, which covered only 0.9% of the region but supported 45% of the koalas. We also estimated that between 1969 and 1995 -30% of the native vegetation associations that are considered as potential koala habitat were cleared, leading to a decline of perhaps 10% in koala numbers. Management of this large regional population has significant implications for the national conservation of the species: the continued viability of this population is critically dependent on the retention and management of riverine and residual vegetation communities, and future vegetation-management guidelines should be cognisant of the potential impacts of clearing even small areas of critical habitat. We also highlight eight management implications.

Pheromone-mediated mass trapping and population diversion as strategies for suppressing Carpophilus spp. (Coleoptera: Nitidulidae) in Australian stone fruit orchards

1 Five experiments were conducted during 1995-99 in stone fruit orchards on the Central Coast and in inland New South Wales, Australia, on the use of synthetic aggregation pheromones and a coattractant to suppress populations of the ripening fruit pests Carpophilus spp. (Coleoptera: Nitidulidae). 2 Perimeter-based suppression traps baited with pheromone and coattractant placed at 3m intervals around small fruit blocks, caught large numbers of Carpophilus spp. Very small populations of Carpophilus spp. occurred within blocks, and fruit damage was minimal. 3 Carpophilus spp. populations in stone fruit blocks 15-370m from suppression traps were also small and non-damaging, indicating a large zone of pheromone attractivity. 4 Pheromone/coattractant-baited suppression traps appeared to divert Carpophilus spp. from nearby (130 m) ripening stone fruit. Ten metal drums containing decomposing fruit, baited with pheromone and treated with insecticide, attracted Carpophilus spp. and appeared to reduce populations and damage to ripening fruit at distances of 200-500 m. Populations and damage were significantly greater within 200m of the drums and may have been caused by ineffective poisoning or poor quality/overcrowding of fruit resources in the drums. 5 Suppression of Carpophilus spp. populations using synthetic aggregation pheromones and a coattractant appears to be a realistic management option in stone fruit orchards. Pheromone-mediated diversion of beetle populations from ripening fruit may be more practical than perimeter trapping, but more research is needed on the effective range of Carpophilus pheromones and the relative merits of trapping compared to attraction to insecticide-treated areas.

Stock assessment of the Queensland-New South Wales tailor fishery (Pomatomus saltatrix)

Data on catch sizes, catch rates, length-frequency and age composition from the Australian east coast tailor fishery are analysed by three different population dynamic models: a surplus production model, an age-structured model, and a model in which the population is structured by both age and length. The population is found to be very heavily exploited, with its ability to reproduce dependent on the fishery’s incomplete selectivity of one-year-old fish. Estimates of recent harvest rates (proportion of fish available to the fishery that are actually caught in a single year) are over 80%. It is estimated that only 30–50% of one-year-old fish are available to the fishery. Results from the age-length-structured model indicate that both exploitable biomass (total mass of fish selected by the fishery) and egg production have fallen to about half the levels that prevailed in the 1970s, and about 40% of virgin levels. Two-year-old fish appear to have become smaller over the history of the fishery. This is assumed to be due to increased fishing pressure combined with non-selectivity of small one-year-old fish, whereby the one-year-old fish that survive fishing are small and grow into small two-year-old fish the following year. An alternative hypothesis is that the stock has undergone a genetic change towards smaller fish; the true explanation is unknown. The instantaneous natural mortality rate of tailor is hypothesised to be higher than previously thought, with values between 0.8 and 1.3 yr–1 consistent with the models. These values apply only to tailor up to about three years of age, and it is possible that a lower value applies to fish older than three. The analysis finds no evidence that fishing pressure has yet affected recruitment. If a recruitment downturn were to occur, however, under current management and fishing pressure there is a strong chance that the fishery would need a complete closure for several years to recover, and even then recovery would be uncertain. Therefore it is highly desirable to better protect the spawning stock. The major recommendations are • An increase in the minimum size limit from 30cm to 40cm in order to allow most one-year-old fish to spawn, and • An experiment on discard mortality to gauge the proportion of fish between 30cm and 40cm that are likely to survive being caught and released by recreational line fishers (the dominant component of the fishery, currently harvesting roughly 1000t p.a. versus about 200t p.a. from the commercial fishery).

Population ecology of Heteronyx piceus (Coleoptera: Scarabaeidae) in a peanut/maize cropping system

Large larval populations of the scarabaeid beetle Heteronyx piceus Blanchard that occur under peanuts, but not maize, in the South Burnett region of Australia are the result of a high rate and prolonged period of egg production by females feeding on peanut foliage. Heteronyx piceus is a relatively sedentary species and movement of females between adjacent fields is low. Populations of H. piceus varied markedly with landscape position. High larval populations are more likely (1 in 4 chance) to be encountered on the ‘scrub’ soils in the upper parts of the landscape than in the ‘forest’ soils in the lower half (1 in 20 chance), indicating that soil type/landscape position is a key risk factor in assessing the need for management intervention. The studies indicate that, because of the species' sedentary nature, the most meaningful population entity for management of H. piceus is the individual field, rather than the whole-farm or the region. The implications of this population ecology for management of the pest are discussed in relation to control strategies.

The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus)

This study compares estimates of the census size of the spawning population with genetic estimates of effective current and long-term population size for an abundant and commercially important marine invertebrate, the brown tiger prawn (Penaeus esculentus). Our aim was to focus on the relationship between genetic effective and census size that may provide a source of information for viability analyses of naturally occurring populations. Samples were taken in 2001, 2002 and 2003 from a population on the east coast of Australia and temporal allelic variation was measured at eight polymorphic microsatellite loci. Moments-based and maximum-likelihood estimates of current genetic effective population size ranged from 797 to 1304. The mean long-term genetic effective population size was 9968. Although small for a large population, the effective population size estimates were above the threshold where genetic diversity is lost at neutral alleles through drift or inbreeding. Simulation studies correctly predicted that under these experimental conditions the genetic estimates would have non-infinite upper confidence limits and revealed they might be overestimates of the true size. We also show that estimates of mortality and variance in family size may be derived from data on average fecundity, current genetic effective and census spawning population size, assuming effective population size is equivalent to the number of breeders. This work confirms that it is feasible to obtain accurate estimates of current genetic effective population size for abundant Type III species using existing genetic marker technology.

Fusarium wilt of cotton: Population diversity and implications for management

Fusarium wilt of cotton, caused by the fungus Fusarium oxysporum Schlechtend. f. sp. vasinfectum (Atk.) Snyd. & Hans, was first identified in 1892 in cotton growing in sandy acid soils in Alabama (8). Although the disease was soon discovered in other major cotton-producing areas, it did not become global until the end of the next century. After its original discovery, Fusarium wilt of cotton was reported in Egypt (1902) (30), India (1908) (60), Tanzania (1954) (110), California (1959) (33), Sudan (1960) (44), Israel (1970) (27), Brazil (1978) (5), China (1981) (17), and Australia (1993) (56). In addition to a worldwide distribution, Fusarium wilt occurs in all four of the domesticated cottons, Gossypium arboretum L., G. barbadense L., G. herbaceum L., and G. hirsutum L. (4,30). Disease losses in cotton are highly variable within a country or region. In severely infested fields planted with susceptible cultivars, yield losses can be high. In California, complete crop losses in individual fields have been observed (R. M. Davis, unpublished). Disease loss estimates prepared by the National Cotton Disease Council indicate losses of over 109,000 bales (227 kg or 500 lb) in the United States in 2004 (12).

Identification of QTL for sugar-related traits in a sweet x grain sorghum (Sorghum bicolor L. Moench) recombinant inbred population

QTL for stem sugar-related and other agronomic traits were identified in a converted sweet (R9188) × grain (R9403463-2-1) sorghum population. QTL analyses were conducted using phenotypic data for 11 traits measured in two field experiments and a genetic map comprising 228 SSR and AFLP markers grouped into 16 linkage groups, of which 11 could be assigned to the 10 sorghum chromosomes (SBI-01 to SBI-10). QTL were identified for all traits and were generally co-located to five locations (SBI-01, SBI-03, SBI-05, SBI-06 and SBI-10). QTL alleles from R9188 were detected for increased sucrose content and sugar content on SBI-01, SBI-05 and SBI-06. R9188 also contributed QTL alleles for increased Brix on SBI-05 and SBI-06, and increased sugar content on SBI-03. QTL alleles from R9403463-2-1 were found for increased sucrose content and sucrose yield on SBI-10, and increased glucose content on SBI-07. QTL alleles for increased height, later flowering and greater total dry matter yield were located on SBI-01 of R9403463-2-1, and SBI-06 of R9188. QTL alleles for increased grain yield from both R9403463-2-1 and R9188 were found on SBI-03. As an increase in stem sugars is an important objective in sweet sorghum breeding, the QTL identified in this study could be further investigated for use in marker-assisted selection of sweet sorghum.

Genetic population structure of red snappers (Lutjanus malabaricus Bloch & Schneider, 1801 and Lutjanus erythropterus Bloch, 1790) in central and eastern Indonesia and northern Australia

The genetic population structure of red snapper Lutjanus malabaricus and Lutjanus erythropterus in eastern Indonesia and northern Australia was investigated by allozyme electrophoresis and sequence variation in the control region of mtDNA. Samples were collected from eight sites in Indonesia and four sites in northern Australia for both species. A total of 13 allozyme loci were scored. More variable loci were observed in L. malabaricus than in L. erythropterus. Sequence variation in the control region (left domain) of the mitochondrial genome was assessed by RFLP and direct sequencing. MtDNA haplotype diversity was high (L. erythropterus, 0.95 and L. malabaricus, 0.97), as was intraspecific sequence divergence, (L. erythropterus, 0.0-12.5% and L. malabaricus, 0.0-9.5%). The pattern of mtDNA haplotype frequencies grouped both species into two broad fisheries stocks with a genetic boundary either between Kupang and Sape (L. malabaricus) or between Kupang and Australian Timor Sea (L. erythropertus). The allozyme analyses revealed similar boundaries for L. erythropterus. Seven allozymes stocks compared to two mtDNA stocks of L. malabaricus including Ambon, which was not sampled with mtDNA, however, were reported. Possible reasons for differences in discrimination between the methods include: i) increased power of multiple allozyme loci over the single mtDNA locus, ii) insufficient gene sampling in the mtDNA control region and iii) relative evolutionary dynamics of nuclear (allozyme loci) and mitochondrial DNA in these taxa. Allozyme and haplotype data did not distinguish separate stocks among the four Australian locations nor the central Indonesian (Bali and Sape locations) for both L. malabaricus and L. erythropterus.

Sampling and management of Bemisia tabaci (Genn.) biotype B in Australian cotton

Data on seasonal population abundance of Bemisia tabaci biotype B (silverleaf whitefly (SLW)) in Australian cotton fields collected over four consecutive growing seasons (2002/2003-2005/2006) were used to develop and validate a multiple-threshold-based management and sampling plan. Non-linear growth trajectories estimated from the field sampling data were used as benchmarks to classify adult SLW field populations into six density-based management zones with associated control recommendations in the context of peak flowering and open boll crop growth stages. Control options based on application of insect growth regulators (IGRs) are recommended for high-density populations (>2 adults/leaf) whereas conventional (non-IGR) products are recommended for the control of low to moderate population densities. A computerised re-sampling program was used to develop and test a binomial sampling plan. Binomial models with thresholds of T=1, 2 and 3 adults/leaf were tested using the field abundance data. A binomial plan based on a tally threshold of T=2 adults/leaf and a minimum sample of 20 leaves at nodes 3, 4 or 5 below the terminal is recommended as the most parsimonious and practical sampling protocol for Australian cotton fields. A decision support guide with management zone boundaries expressed as binomial counts and control options appropriate for various SLW density situations is presented. Appropriate use of chemical insecticides and tactics for successful field control of whiteflies are discussed.

Population genetic structure of the brown tiger prawn, Penaeus esculentus, in tropical northern Australia

Eight polymorphic microsatellite loci were analysed in six population samples from four locations of the Australian endemic brown tiger prawn, Penaeus esculentus. Tests of Hardy-Weinberg equilibrium were generally in accord with expectations, with only one locus, in two samples, showing significant deviations. Three samples were taken in different years from the Exmouth Gulf. These showed no significant heterogeneity, and it was concluded that they were from a single panmictic population. A sample from Shark Bay, also on the west coast of Australia, showed barely detectable differentiation from Exmouth Gulf (F (ST) = 0 to 0.0014). A northeast sample from the Gulf of Carpentaria showed low (F (ST) = 0.008) but significant differentiation from Moreton Bay, on the east coast. However, Exmouth Gulf/Shark Bay samples were well differentiated from the Gulf of Carpentaria/Moreton Bay (F (ST) = 0.047-0.063). The data do not fit a simple isolation by distance model. It is postulated that the east-west differentiation largely reflects the isolation of east and west coast populations that occurred at the last glacial maximum when there was a land bridge between north-eastern Australia and New Guinea.

Optimal management of agricultural systems

To remain competitive, many agricultural systems are now being run along business lines. Systems methodologies are being incorporated, and here evolutionary computation is a valuable tool for identifying more profitable or sustainable solutions. However, agricultural models typically pose some of the more challenging problems for optimisation. This chapter outlines these problems, and then presents a series of three case studies demonstrating how they can be overcome in practice. Firstly, increasingly complex models of Australian livestock enterprises show that evolutionary computation is the only viable optimisation method for these large and difficult problems. On-going research is taking a notably efficient and robust variant, differential evolution, out into real-world systems. Next, models of cropping systems in Australia demonstrate the challenge of dealing with competing objectives, namely maximising farm profit whilst minimising resource degradation. Pareto methods are used to illustrate this trade-off, and these results have proved to be most useful for farm managers in this industry. Finally, land-use planning in the Netherlands demonstrates the size and spatial complexity of real-world problems. Here, GIS-based optimisation techniques are integrated with Pareto methods, producing better solutions which were acceptable to the competing organizations. These three studies all show that evolutionary computation remains the only feasible method for the optimisation of large, complex agricultural problems. An extra benefit is that the resultant population of candidate solutions illustrates trade-offs, and this leads to more informed discussions and better education of the industry decision-makers.

Thermal requirements, field mortality and population phenology modelling of Paropsis atomaria Olivier, an emergent pest in subtropical hardwood plantations

Paropsis atomaria is a recently emerged pest of eucalypt plantations in subtropical Australia. Its broad host range of at least 20 eucalypt species and wide geographical distribution provides it the potential to become a serious forestry pest both within Australia and, if accidentally introduced, overseas. Although populations of P. atomaria are genetically similar throughout its range, population dynamics differ between regions. Here, we determine temperature-dependent developmental requirements using beetles sourced from temperate and subtropical zones by calculating lower temperature thresholds, temperature-induced mortality, and day-degree requirements. We combine these data with field mortality estimates of immature life stages to produce a cohort-based model, ParopSys, using DYMEX™ that accurately predicts the timing, duration, and relative abundance of life stages in the field and number of generations in a spring–autumn (September–May) field season. Voltinism was identified as a seasonally plastic trait dependent upon environmental conditions, with two generations observed and predicted in the Australian Capital Territory, and up to four in Queensland. Lower temperature thresholds for development ranged between 4 and 9 °C, and overall development rates did not differ according to beetle origin. Total immature development time (egg–adult) was approximately 769.2 ± S.E. 127.8 DD above a lower temperature threshold of 6.4 ± S.E. 2.6 °C. ParopSys provides a basic tool enabling forest managers to use the number of generations and seasonal fluctuations in abundance of damaging life stages to estimate the pest risk of P. atomaria prior to plantation establishment, and predict the occurrence and duration of damaging life stages in the field. Additionally, by using local climatic data the pest potential of P. atomaria can be estimated to predict the risk of it establishing if accidentally introduced overseas. Improvements to ParopSys’ capability and complexity can be made as more biological data become available.