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.

Dingoes at the Doorstep: Home Range Sizes and Activity Patterns of Dingoes and Other Wild Dogs around Urban Areas of North-Eastern Australia

Top-predators around the world are becoming increasingly intertwined with humans, sometimes causing conflict and increasing safety risks in urban areas. In Australia, dingoes and dingo � domestic dog hybrids are common in many urban areas, and pose a variety of human health and safety risks. However, data on urban dingo ecology is scant. We GPS-collared 37 dingoes in north-eastern Australia and continuously monitored them each 30 min for 11–394 days. Most dingoes were nocturnal, with an overall mean home range size of 17.47 km2. Overall mean daily distance travelled was 6.86 km/day. At all times dingoes were within 1000 m of houses and buildings. Home ranges appeared to be constrained to patches of suitable vegetation fragments within and around human habitation. These data can be used to reallocate dingo management effort towards mitigating actual conflicts between humans and dingoes in urban areas.

Temporal and Spatial Patterns in Recruitment of Three Penaeid Prawns in Moreton Bay, Queensland, Australia

The size at recruitment, temporal and spatial distribution, and abiotic factors influencing abundance of three commercially important species of penaeid prawns in the sublittoral trawl grounds of Moreton Bay (Queensland, Australia) were compared. Metapenaeus bennettae and Penaeus plebejus recruit to the trawl grounds at sizes which are relatively small (14-15 mm carapace length, CL) and below that at which prawns are selected for, and retained, in the fleet's cod-ends. In contrast, Penaeus esculenlus recruit at the relatively large size of 27 mm CL from February to May, well above the size ranges selected for. Recruitment of M. bennettae extends over several months, September-October and February March, and was thus likely to be bi-annual, while the recruitment period of P. plebejus was distinct, peaking in October-November each year. Size classes of M . bennettae were the most spatially stratified of the three species. Catch rates of recruits were negatively correlated with depth for all three species, and were also negatively correlated with salinity for M. bennettae.

Pesticide residues in food - acute dietary exposure

Consumer risk assessment is a crucial step in the regulatory approval of pesticide use on food crops. Recently, an additional hurdle has been added to the formal consumer risk assessment process with the introduction of short-term intake or exposure assessment and a comparable short-term toxicity reference, the acute reference dose. Exposure to residues during one meal or over one day is important for short-term or acute intake. Exposure in the short term can be substantially higher than average because the consumption of a food on a single occasion can be very large compared with typical long-term or mean consumption and the food may have a much larger residue than average. Furthermore, the residue level in a single unit of a fruit or vegetable may be higher by a factor (defined as the variability factor, which we have shown to be typically ×3 for the 97.5th percentile unit) than the average residue in the lot. Available marketplace data and supervised residue trial data are examined in an investigation of the variability of residues in units of fruit and vegetables. A method is described for estimating the 97.5th percentile value from sets of unit residue data. Variability appears to be generally independent of the pesticide, the crop, crop unit size and the residue level. The deposition of pesticide on the individual unit during application is probably the most significant factor. The diets used in the calculations ideally come from individual and household surveys with enough consumers of each specific food to determine large portion sizes. The diets should distinguish the different forms of a food consumed, eg canned, frozen or fresh, because the residue levels associated with the different forms may be quite different. Dietary intakes may be calculated by a deterministic method or a probabilistic method. In the deterministic method the intake is estimated with the assumptions of large portion consumption of a ‘high residue’ food (high residue in the sense that the pesticide was used at the highest recommended label rate, the crop was harvested at the smallest interval after treatment and the residue in the edible portion was the highest found in any of the supervised trials in line with these use conditions). The deterministic calculation also includes a variability factor for those foods consumed as units (eg apples, carrots) to allow for the elevated residue in some single units which may not be seen in composited samples. In the probabilistic method the distribution of dietary consumption and the distribution of possible residues are combined in repeated probabilistic calculations to yield a distribution of possible residue intakes. Additional information such as percentage commodity treated and combination of residues from multiple commodities may be incorporated into probabilistic calculations. The IUPAC Advisory Committee on Crop Protection Chemistry has made 11 recommendations relating to acute dietary exposure.

Variability in Voluntary Intake of Loose Mineral Mix Supplements by Grazing Heifers

Loose mineral mix (LMM) supplements are often fed to ruminants in extensive grazing situations to provide minerals and nitrogen likely to be deficient in pasture. However a large proportion of animals offered such supplements may not consume any supplement, while among consumer animals the variability in supplement intake may be high (Wheeler et al., 1980; Dixon et al., 1996). Two experiments examined the distribution of intake of LMM supplements offered to heifers grazing in mob and paddock sizes representative of commercial cattle properties.

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).

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 structure of Australasian green turtles (Chelonia mydas): Exploring the geographical scale of genetic exchange

Ecological and genetic studies of marine turtles generally support the hypothesis of natal homing, but leave open the question of the geographical scale of genetic exchange and the capacity of turtles to shift breeding sites. Here we combine analyses of mitochondrial DNA (mtDNA) variation and recapture data to assess the geographical scale of individual breeding populations and the distribution of such populations through Australasia. We conducted multiscale assessments of mtDNA variation among 714 samples from 27 green turtle rookeries and of adult female dispersal among nesting sites in eastern Australia. Many of these rookeries are on shelves that were flooded by rising sea levels less than 10 000 years (c. 450 generations) ago. Analyses of sequence variation among the mtDNA control region revealed 25 haplotypes, and their frequency distributions indicated 17 genetically distinct breeding stocks (Management Units) consisting either of individual rookeries or groups of rookeries in general that are separated by more than 500 km. The population structure inferred from mtDNA was consistent with the scale of movements observed in long-term mark-recapture studies of east Australian rookeries. Phylogenetic analysis of the haplotypes revealed five clades with significant partitioning of sequence diversity (Φ = 68.4) between Pacific Ocean and Southeast Asian/Indian Ocean rookeries. Isolation by distance was indicated for rookeries separated by up to 2000 km but explained only 12% of the genetic structure. The emerging general picture is one of dynamic population structure influenced by the capacity of females to relocate among proximal breeding sites, although this may be conditional on large population sizes as existed historically across this region.

A high-density consensus map of barley linking DArT markers to SSR, RFLP and STS loci and agricultural traits

Background: Molecular marker technologies are undergoing a transition from largely serial assays measuring DNA fragment sizes to hybridization-based technologies with high multiplexing levels. Diversity Arrays Technology (DArT) is a hybridization-based technology that is increasingly being adopted by barley researchers. There is a need to integrate the information generated by DArT with previous data produced with gel-based marker technologies. The goal of this study was to build a high-density consensus linkage map from the combined datasets of ten populations, most of which were simultaneously typed with DArT and Simple Sequence Repeat (SSR), Restriction Enzyme Fragment Polymorphism (RFLP) and/or Sequence Tagged Site (STS) markers. Results: The consensus map, built using a combination of JoinMap 3.0 software and several purpose-built perl scripts, comprised 2,935 loci (2,085 DArT, 850 other loci) and spanned 1,161 cM. It contained a total of 1,629 'bins' (unique loci), with an average inter-bin distance of 0.7 ± 1.0 cM (median = 0.3 cM). More than 98% of the map could be covered with a single DArT assay. The arrangement of loci was very similar to, and almost as optimal as, the arrangement of loci in component maps built for individual populations. The locus order of a synthetic map derived from merging the component maps without considering the segregation data was only slightly inferior. The distribution of loci along chromosomes indicated centromeric suppression of recombination in all chromosomes except 5H. DArT markers appeared to have a moderate tendency toward hypomethylated, gene-rich regions in distal chromosome areas. On the average, 14 ± 9 DArT loci were identified within 5 cM on either side of SSR, RFLP or STS loci previously identified as linked to agricultural traits. Conclusion: Our barley consensus map provides a framework for transferring genetic information between different marker systems and for deploying DArT markers in molecular breeding schemes. The study also highlights the need for improved software for building consensus maps from high-density segregation data of multiple populations.

A comparison of plotless density estimators using Monte Carlo simulation on totally enumerated field data sets

Background: Plotless density estimators are those that are based on distance measures rather than counts per unit area (quadrats or plots) to estimate the density of some usually stationary event, e.g. burrow openings, damage to plant stems, etc. These estimators typically use distance measures between events and from random points to events to derive an estimate of density. The error and bias of these estimators for the various spatial patterns found in nature have been examined using simulated populations only. In this study we investigated eight plotless density estimators to determine which were robust across a wide range of data sets from fully mapped field sites. They covered a wide range of situations including animal damage to rice and corn, nest locations, active rodent burrows and distribution of plants. Monte Carlo simulations were applied to sample the data sets, and in all cases the error of the estimate (measured as relative root mean square error) was reduced with increasing sample size. The method of calculation and ease of use in the field were also used to judge the usefulness of the estimator. Estimators were evaluated in their original published forms, although the variable area transect (VAT) and ordered distance methods have been the subjects of optimization studies. Results: An estimator that was a compound of three basic distance estimators was found to be robust across all spatial patterns for sample sizes of 25 or greater. The same field methodology can be used either with the basic distance formula or the formula used with the Kendall-Moran estimator in which case a reduction in error may be gained for sample sizes less than 25, however, there is no improvement for larger sample sizes. The variable area transect (VAT) method performed moderately well, is easy to use in the field, and its calculations easy to undertake. Conclusion: Plotless density estimators can provide an estimate of density in situations where it would not be practical to layout a plot or quadrat and can in many cases reduce the workload in the field.

Effect of changes in hook pattern and size on catch rate, hooking location, injury and bleeding for a number of tropical reef fish species

The Queensland Great Barrier Reef line fishery in Australia is regulated via a range of input and output controls including minimum size limits, daily catch limits and commercial catch quotas. As a result of these measures a substantial proportion of the catch is released or discarded. The fate of these released fish is uncertain, but hook-related mortality can potentially be decreased by using hooks that reduce the rates of injury, bleeding and deep hooking. There is also the potential to reduce the capture of non-target species though gear selectivity. A total of 1053 individual fish representing five target species and three non-target species were caught using six hook types including three hook patterns (non-offset circle, J and offset circle), each in two sizes (small 4/0 or 5/0 and large 8/0). Catch rates for each of the hook patterns and sizes varied between species with no consistent results for target or non-target species. When data for all of the fish species were aggregated there was a trend for larger hooks, J hooks and offset circle hooks to cause a greater number of injuries. Using larger hooks was more likely to result in bleeding, although this trend was not statistically significant. Larger hooks were also more likely to foul-hook fish or hook fish in the eye. There was a reduction in the rates of injuries and bleeding for both target and non-target species when using the smaller hook sizes. For a number of species included in our study the incidence of deep hooking decreased when using non-offset circle hooks, however, these results were not consistent for all species. Our results highlight the variability in hook performance across a range of tropical demersal finfish species. The most obvious conservation benefits for both target and non-target species arise from using smaller sized hooks and non-offset circle hooks. Fishers should be encouraged to use these hook configurations to reduce the potential for post-release mortality of released fish.

Ghost fishing in the pot fishery for blue swimmer crabs Portunus pelagicus in Queensland, Australia

Blue swimmer crabs (Portunus pelagicus) are an economically important crab caught in baited traps throughout the Indo-west Pacific and Mediterranean. In Australia they are traditionally caught using rigid wire traps (approximate to pots) but there has been a recent increase in the use of collapsible pots constructed from polyethylene trawl mesh. Two experiments were conducted in Moreton Bay, Queensland, to determine the ghost fishing potential of lost crab pots on both target and bycatch species and to evaluate the differences between traditional and contemporary pot designs. A lost contemporary, collapsible trawl mesh pot will catch between 3 and 223 R pelagicus per year after the bait has been exhausted, while a traditional wire mesh pot would catch 11-74 crabs peryear. As most fishers now use the collapsible trawl mesh pots, ghost fishing mortality could be as high as 111,811-670,866 crabs per year. Bycatch retention was also higher in contemporary designs. Periods of strong winds appeared to increase the ghost fishing potential of lost pots. The use of escape gaps, larger mesh sizes and construction options that allow for the deterioration of entrance funnels to minimise ghost fishing are recommended to reduce environmental impacts.

Parameter estimation and sensitivity analysis of fat deposition models in beef steers using acslXtreme.

The Davis Growth Model (a dynamic steer growth model encompassing 4 fat deposition models) is currently being used by the phenotypic prediction program of the Cooperative Research Centre (CRC) for Beef Genetic Technologies to predict P8 fat (mm) in beef cattle to assist beef producers meet market specifications. The concepts of cellular hyperplasia and hypertrophy are integral components of the Davis Growth Model. The net synthesis of total body fat (kg) is calculated from the net energy available after accounting tor energy needs for maintenance and protein synthesis. Total body fat (kg) is then partitioned into 4 fat depots (intermuscular, intramuscular, subcutaneous, and visceral). This paper reports on the parameter estimation and sensitivity analysis of the DNA (deoxyribonucleic acid) logistic growth equations and the fat deposition first-order differential equations in the Davis Growth Model using acslXtreme (Hunstville, AL, USA, Xcellon). The DNA and fat deposition parameter coefficients were found to be important determinants of model function; the DNA parameter coefficients with days on feed >100 days and the fat deposition parameter coefficients for all days on feed. The generalized NL2SOL optimization algorithm had the fastest processing time and the minimum number of objective function evaluations when estimating the 4 fat deposition parameter coefficients with 2 observed values (initial and final fat). The subcutaneous fat parameter coefficient did indicate a metabolic difference for frame sizes. The results look promising and the prototype Davis Growth Model has the potential to assist the beef industry meet market specifications.

Hardness Methods for Testing Maize Kernels

Maize is a highly important crop to many countries around the world, through the sale of the maize crop to domestic processors and subsequent production of maize products and also provides a staple food to subsistance farms in undeveloped countries. In many countries, there have been long-term research efforts to develop a suitable hardness method that could assist the maize industry in improving efficiency in processing as well as possibly providing a quality specification for maize growers, which could attract a premium. This paper focuses specifically on hardness and reviews a number of methodologies as well as important biochemical aspects of maize that contribute to maize hardness used internationally. Numerous foods are produced from maize, and hardness has been described as having an impact on food quality. However, the basis of hardness and measurement of hardness are very general and would apply to any use of maize from any country. From the published literature, it would appear that one of the simpler methods used to measure hardness is a grinding step followed by a sieving step, using multiple sieve sizes. This would allow the range in hardness within a sample as well as average particle size and/or coarse/fine ratio to be calculated. Any of these parameters could easily be used as reference values for the development of near-infrared (NIR) spectroscopy calibrations. The development of precise NIR calibrations will provide an excellent tool for breeders, handlers, and processors to deliver specific cultivars in the case of growers and bulk loads in the case of handlers, thereby ensuring the most efficient use of maize by domestic and international processors. This paper also considers previous research describing the biochemical aspects of maize that have been related to maize hardness. Both starch and protein affect hardness, with most research focusing on the storage proteins (zeins). Both the content and composition of the zein fractions affect hardness. Genotypes and growing environment influence the final protein and starch content and. to a lesser extent, composition. However, hardness is a highly heritable trait and, hence, when a desirable level of hardness is finally agreed upon, the breeders will quickly be able to produce material with the hardness levels required by the industry.