Impacts of grazing management options on pasture and animal productivity in a Heteropogon contortus (black speargrass) pasture in central Queensland. 2. Population dynamics of Heteropogon contortus and Stylosanthes scabra cv. Seca.

The dynamics of Heteropogon contortus and Stylosanthes scabra cv. Seca populations were studied in a subset of treatments in an extensive grazing study conducted in central Queensland between 1988 and 2001. These treatments were 4 stocking rates in native pasture and 2 of these stocking rates in legume oversown and supplement/spring burning treatments. For the 1999-2000 summer, population data for H. contortus in 5 of these native pasture and supplement/burning treatments were compared with those for an additional burnt treatment. Seasonal rainfall throughout this study was below the long-term mean and mean annual pasture utilisation ranged from 24 to 61%. Increasing stocking rate from 5 to 2 ha/steer in native pasture reduced H. contortus plant density. Increasing stocking rate reduced seedling recruitment as a result of its effect on soil seedbanks. Seedling recruitment was the major determinant of change in plant density, although some individual H. contortus plants did survive throughout the study. Burning in spring 1999, particularly at light stocking rate, promoted seedling recruitment above that in both unburnt native and legume oversown pasture and resulted in increased H. contortus plant density. In the legume oversown treatments, S. scabra cv. Seca density increased rapidly from 15 plants/m2 in 1988 to 140 plants/m2 in 2001 following a lag phase between 1988 and 1993. This increased S. scabra density was associated with an eventual decline in H. contortus plant density through reduced seedling recruitment. It was concluded that H. contortus population density is sustainable at stocking rates of 4 and 5 ha/steer (30% pasture utilisation) and that spring burning at light stocking rate can promote H. contortus populations. Increasing densities of S. scabra need to be managed to prevent its dominance.

Raw data for "Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson)."

Abstract of Macbeth, G. M., Broderick, D., Buckworth, R. & Ovenden, J. R. (In press, Feb 2013). Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson). G3: Genes, Genomes and Genetics. Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. But, pitfalls are predicted when the effective size is large, as estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes around 5000. This was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest standard deviation of Ne estimates occurred when low frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or non-breeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were above 24,000.

NEESTIMATOR v2: re-implementation of software for the estimation of contemporary effective population size (N-e) from genetic data

NeEstimator v2 is a completely revised and updated implementation of software that produces estimates of contemporary effective population size, using several different methods and a single input file. NeEstimator v2 includes three single-sample estimators (updated versions of the linkage disequilibrium and heterozygote-excess methods, and a new method based on molecular coancestry), as well as the two-sample (moment-based temporal) method. New features include the following: (i) an improved method for accounting for missing data; (ii) options for screening out rare alleles; (iii) confidence intervals for all methods; (iv) the ability to analyse data sets with large numbers of genetic markers (10000 or more); (v) options for batch processing large numbers of different data sets, which will facilitate cross-method comparisons using simulated data; and (vi) correction for temporal estimates when individuals sampled are not removed from the population (Plan I sampling). The user is given considerable control over input data and composition, and format of output files. The freely available software has a new JAVA interface and runs under MacOS, Linux and Windows.

A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

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.

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.

Concordant molecular and phenotypic data delineate new taxonomy and conservation priorities for the endangered mountain yellow-legged frog

The mountain yellow-legged frog Rana muscosa sensu lato, once abundant in the Sierra Nevada of California and Nevada, and the disjunct Transverse Ranges of southern California, has declined precipitously throughout its range, even though most of its habitat is protected. The species is now extinct in Nevada and reduced to tiny remnants in southern California, where as a distinct population segment, it is classified as Endangered. Introduced predators (trout), air pollution and an infectious disease (chytridiomycosis) threaten remaining populations. A Bayesian analysis of 1901 base pairs of mitochondrial DNA confirms the presence of two deeply divergent clades that come into near contact in the Sierra Nevada. Morphological studies of museum specimens and analysis of acoustic data show that the two major mtDNA clades are readily differentiated phenotypically. Accordingly, we recognize two species, Rana sierrae, in the northern and central Sierra Nevada, and R. muscosa, in the southern Sierra Nevada and southern California. Existing data indicate no range overlap. These results have important implications for the conservation of these two species as they illuminate a profound mismatch between the current delineation of the distinct population segments (southern California vs. Sierra Nevada) and actual species boundaries. For example, our study finds that remnant populations of R. muscosa exist in both the southern Sierra Nevada and the mountains of southern California, which may broaden options for management. In addition, despite the fact that only the southern California populations are listed as Endangered, surveys conducted since 1995 at 225 historic (1899-1994) localities from museum collections show that 93.3% (n=146) of R. sierrae populations and 95.2% (n=79) of R. muscosa populations are extinct. Evidence presented here underscores the need for revision of protected population status to include both species throughout their ranges.

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.

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.

Mapping quantitative trait loci for resistance to Pratylenchus thornei from synthetic hexaploid wheat in the International Triticeae Mapping Initiative (ITMI) population

Root-lesion nematode (Pratylenchus thornei) is a serious pathogen of wheat in many countries. The International Triticeae Mapping Initiative (ITMI) population of recombinant inbred lines (RILs) was assessed for resistance to P. thornei to determine the chromosome locations of the resistance genes. The ITMI population is derived from a cross between the resistant synthetic hexaploid wheat W-7984 and a susceptible bread wheat cultivar Opata 85. Two years of phenotypic data for resistance to P. thornei were obtained in replicated glasshouse trials. Quantitative trait locus (QTL) analysis was performed using available segregation and map data for 114 RILs. A QTL on chromosome 6DS showed consistent effects for reduced nematode numbers (partial resistance) across years and accounted for 11% and 23% of the phenotypic variation. A second QTL for P. thornei resistance on chromosome 2BS accounted for an additional 19% and 5%. Restriction fragment length polymorphism (RFLP) and simple sequence repeat (SSR) markers associated with the QTLs are physically located in regions rich in major genes at the distal ends of the short chromosome arms of 6D and 2B. SSR markers with potential for marker-assisted selection of P. thornei resistance effective in different genetic backgrounds have been identified.

Inter-seasonal population dynamics and pest status of Bemisia tabaci (Gennadius) biotype B in an Australian cropping system

Bemisia tabaci, biotype B, commonly known as the silverleaf whitefly (SLW) is an alien species that invaded Australia in the mid-90s. This paper reports on the invasion ecology of SLW and the factors that are likely to have contributed to the first outbreak of this major pest in an Australian cotton cropping system, population dynamics of SLW within whitefly-susceptible crop (cotton and cucurbit) and non-crop vegetation (sowthistle, Sonchus spp.) components of the cropping system were investigated over four consecutive growing seasons (September-June) 2001/02-2004/05 in the Emerald Irrigation Area (EIA) of Queensland, Australia. Based on fixed geo-referenced sampling sites, variation in spatial and temporal abundance of SLW within each system component was quantified to provide baseline data for the development of ecologically sustainable pest management strategies. Parasitism of large (3rd and 4th instars) SLW nymphs by native aphelinid wasps was quantified to determine the potential for natural control of SLW populations. Following the initial outbreak in 2001/02, SLW abundance declined and stabilised over the next three seasons. The population dynamics of SLW is characterised by inter-seasonal population cycling between the non-crop (weed) and cotton components of the EIA cropping system. Cotton was the largest sink for and source of SLW during the study period. Over-wintering populations dispersed from weed host plant sources to cotton in spring followed by a reverse dispersal in late summer and autumn to broad-leaved crops and weeds. A basic spatial source-sink analysis showed that SLW adult and nymph densities were higher in cotton fields that were closer to over-wintering weed sources throughout spring than in fields that were further away. Cucurbit fields were not significant sources of SLW and did not appear to contribute significantly to the regional population dynamics of the pest. Substantial parasitism of nymphal stages throughout the study period indicates that native parasitoid species and other natural enemies are important sources of SLW mortality in Australian cotton production systems. Weather conditions and use of broad-spectrum insecticides for pest control are implicated in the initial outbreak and on-going pest status of SLW in the region.

Demography of feral camels in central Australia and its relevance to population control

Since their release over 100 years ago, camels have spread across central Australia and increased in number. Increasingly, they are being seen as a pest, with observed impacts from overgrazing and damage to infrastructure such as fences. Irregular aerial surveys since 1983 and an interview-based survey in 1966 suggest that camels have been increasing at close to their maximum rate. A comparison of three models of population growth fitted to these, albeit limited, data suggests that the Northern Territory population has indeed been growing at an annual exponential rate of r = 0.074, or 8% per year, with little evidence of a density-dependent brake. A stage-structured model using life history data from a central Australian camel population suggests that this rate approximates the theoretical maximum. Elasticity analysis indicates that adult survival is by far the biggest influence on rate of increase and that a 9% reduction in survival from 96% is needed to stop the population growing. In contrast, at least 70% of mature females need to be sterilised to have a similar effect. In a benign environment, a population of large mammals such as camels is expected to grow exponentially until close to carrying capacity. This will frustrate control programs, because an ever-increasing number of animals will need to be removed for zero growth the longer that culling or harvesting effort is delayed. A population projection for 2008 suggests ~10 500 animals need to be harvested across the Northern Territory. Current harvests are well short of this. The ability of commercial harvesting to control camel populations in central Australia will depend on the value of animals, access to animals and the presence of alternative species to harvest when camels are at low density.

Milkshark (Rhizoprionodon acutus) microsatellite genotype data (six loci) from samples collected from five locations in eastern Australia and central Indonesia.

Biodiversity of sharks in the tropical Indo-Pacific is high, but species-specific information to assist sustainable resource exploitation is scarce. The null hypothesis of population genetic homogeneity was tested for scalloped hammerhead shark (Sphyrna lewini, n=244) and the milkshark (Rhizoprionodon acutus, n=209) from northern and eastern Australia, using nuclear (S. lewini, eight microsatellite loci; R. acutus, six loci) and mitochondrial gene markers (873 base pairs of NADH dehydrogenase subunit 4). We were unable to reject genetic homogeneity for S. lewini, which was as expected based on previous studies of this species. Less expected were similar results for R. acutus, which is more benthic and less vagile than S. lewini. These features are probably driving the genetic break found between Australian and central Indonesian R. acutus (F-statistics; mtDNA, 0.751 to 0.903; microsatellite loci, 0.038 to 0.047). Our results support the spatially-homogeneous management plan for shark species in Queensland, but caution is advised for species yet to be studied.

Evidence for extensive population structure in the white-spotted eagle ray within the Indo-Pacific inferred from mitochondrial gene sequences

The white-spotted eagle ray Aetobatus narinari is a species complex that occurs circumglobally throughout warm-temperate waters. Aetobatus narinari is semi-pelagic and large (up to 300 cm disc width), suggesting high dispersal capabilities and gene flow on a wide spatial scale. Sequence data from two mitochondrial genes, cytochrome b (cytb) and NADH dehydrogenase subunit 4 (ND4), were used to determine the genetic variability within and among 18 sampling locations in the central Indo-Pacific biogeographical region. Populations in the Indo-Pacific were highly genetically structured with c. 70% of the total genetic variation found among three geographical regions (East China Sea, Southeast Asia and Australia). FST was 0.64 for cytb and 0.53 for ND4, with φST values being even larger, that is, 0.78 for cytb and 0.65 for ND4. This high-level genetic partitioning provides strong evidence against extensive gene flow in A. narinari. The degree of genetic population structuring in the Indo-Pacific was similar to that found on a global scale. Global FST was 0.63 for cytb and 0.57 for ND4, and global φST values were 0.94 for cytb and 0.82 for ND4. This suggests that the A. narinari complex may be more speciose than the two or three species proposed to date. Further sampling and genetic analyses are likely to uncover the ‘evolutionarily significant’ and ‘management’ units that are critical to determine the susceptibilities of individual populations to regional fishing pressures and to provide advice on management options. Network analyses showed a close genetic relationship between haplotypes from the central Indo-Pacific and South Africa, providing support for a proposed dispersal pathway from the possible centre of origin of the A. narinari species complex in the Indo-Pacific into the Atlantic Ocean.