A spatially explicit multi-competitor coexistence model of penaeid (shrimp) distribution on the Australian Great Barrier Reef

A spatially explicit multi-competitor coexistence model was developed for meta-populations of prawns (shrimp) occupying habitat patches across the Great Barrier Reef, where dispersal was localised and dispersal rates varied between species. Prawns were modelled as individuals moving to and from patches or cells according to pre-set decision rules. The landscape was simulated as a matrix of cells with each cell having a spatially explicit survival index for each species. Mixed species prawn assemblages moved over this simplified spatially explicit landscape. A low level of chronic random environmental disturbance was assumed (cyclone and tropical storm damage) with additional acute spatially confined disturbance due to commercial trawling, modelled as an increase in mortality affecting inter-specific competition. The general form of the results was for increased disturbance to favour good-colonising "generalist" species at the expense of good-competitor "specialists". Increasing fishing mortality (local patch extinctions) combined with poor colonising ability resulted in low equilibrium abundance for even the best competitor, while in the same circumstances the poorest competitor but best coloniser could have the highest equilibrium abundance. This mimics the switch from high-value prawn species to lower-value prawn species as trawl effort increases, reflected in historic catch and effort logbook data and reported anecdotaly from the north Queensland trawl fleet. To match the observed distribution and behaviour of prawn assemblages, a combination inter-species competition, a spatially explicit landscape, and a defined pattern of disturbance (trawling) was required. Modelling this combination could simulate not only general trends in spatial distribution of each of prawn species but also localised concentrations observed in the survey data

Using simulation techniques to investigate methods to determine resistance of helminths to anthelmintic treatment

The widespread and increasing resistance of internal parasites to anthelmintic control is a serious problem for the Australian sheep and wool industry. As part of control programmes, laboratories use the Faecal Egg Count Reduction Test (FECRT) to determine resistance to anthelmintics. It is important to have confidence in the measure of resistance, not only for the producer planning a drenching programme but also for companies investigating the efficacy of their products. The determination of resistance and corresponding confidence limits as given in anthelmintic efficacy guidelines of the Standing Committee on Agriculture (SCA) is based on a number of assumptions. This study evaluated the appropriateness of these assumptions for typical data and compared the effectiveness of the standard FECRT procedure with the effectiveness of alternative procedures. Several sets of historical experimental data from sheep and goats were analysed to determine that a negative binomial distribution was a more appropriate distribution to describe pre-treatment helminth egg counts in faeces than a normal distribution. Simulated egg counts for control animals were generated stochastically from negative binomial distributions and those for treated animals from negative binomial and binomial distributions. Three methods for determining resistance when percent reduction is based on arithmetic means were applied. The first was that advocated in the SCA guidelines, the second similar to the first but basing the variance estimates on negative binomial distributions, and the third using Wadley’s method with the distribution of the response variate assumed negative binomial and a logit link transformation. These were also compared with a fourth method recommended by the International Co-operation on Harmonisation of Technical Requirements for Registration of Veterinary Medicinal Products (VICH) programme, in which percent reduction is based on the geometric means. A wide selection of parameters was investigated and for each set 1000 simulations run. Percent reduction and confidence limits were then calculated for the methods, together with the number of times in each set of 1000 simulations the theoretical percent reduction fell within the estimated confidence limits and the number of times resistance would have been said to occur. These simulations provide the basis for setting conditions under which the methods could be recommended. The authors show that given the distribution of helminth egg counts found in Queensland flocks, the method based on arithmetic not geometric means should be used and suggest that resistance be redefined as occurring when the upper level of percent reduction is less than 95%. At least ten animals per group are required in most circumstances, though even 20 may be insufficient where effectiveness of the product is close to the cut off point for defining resistance.

The risk of pollen-mediated gene flow from exotic Corymbia plantations into native Corymbia populations in Australia

Sectors of the forest plantation industry in Australia are set to expand in the near future using species or hybrids of the spotted gums (Corymbia, Section Politaria). Plantations of these taxa have already been introduced across temperate and subtropical Australia, representing locally exotic introductions from native stands in Queensland and New South Wales. A literature review was undertaken to provide insights into the potential for pollen-mediated gene flow from these plantations into native populations. Three factors suggest that such gene flow is likely; (1) interspecific hybridisation within the genus has frequently been recorded, including between distantly related species from different sections, (2) apparent high levels of vertebrate pollinator activity may result in plantation pollen being moved over hundreds of kilometres, (3) much of the plantation estate is being established among closely related taxa and therefore few barriers to gene flow are expected. Across Australia, 20 of the 100 native Corymbia taxa were found to have regional level co-occurrence with plantations. These were located most notably within regions of north-east New South Wales and south-east Queensland, however, co-occurrence was also found in south-west Western Australia and eastern Victoria. The native species found to have co-occurrence were then assessed for the presence of reproductive barriers at each step in the process of gene flow that may reduce the number of species at risk even further. The available data suggest three risk categories exist for Corymbia. The highest risk was for gene flow from plantations of spotted gums to native populations of spotted gums. This was based on the expected limited existence of pre- and post-zygotic barriers, substantial long-distance pollen dispersal and an apparent broad period of flowering in Corymbia citriodora subsp. variegata plantations. The following risk category focussed on gene flow from Corymbia torelliana × C. c. variegata hybrid plantations into native C. c. variegata, as the barriers associated with the production and establishment of F1 hybrids have been circumvented. For the lowest risk category, Corymbia plantations may present a risk to other non-spotted gum species, however, further investigation of the particular cross-combinations is required. A list of research directions is provided to better quantify these risks. Empirical data will need to be combined within a risk assessment framework that will not only estimate the likelihood of exotic gene flow, but also consider the conservation status/value of the native populations. In addition, the potential impacts of pollen flow from plantations will need to be weighed up against their various economic and environmental benefits.

Relationship Between Hardness Genes and Quality in Barley (Hordeum vulgare)

Barley (Hordeum vulgare) genotypes were sequenced for polymorphism in the hardness genes, these being the three hordoindoline (hin a, hin b1 and hin b2) genes. The variation in haplotype was determined by sequencing for single nucleotide polymorphisms (SNPs). Polymorphism between each gene was then compared to grain hardness (three methods), malt quality characteristics (hot water extract and friability) and cattle feed quality. Two haplotypes were found in a set of forty barley genotypes. For hin a, two alleles were present, namely hin a1 and hin a2. However, there was no specific hin a allele that was associated with grain hardness, malt and feed quality. Barley has two hin b genes, namely hin b1 and hin b2, and the genotypes tested here had one of two alleles for each gene. However, there were no obvious effects on hardness or quality from either of these hin b alleles. Unlike wheat, where a clear relationship has been demonstrated between a number of SNPs in the wheat hardness genes and quality (soft or hard wheat), there was no such relationship for barley. Despite the wide range in hardness, malt and feed quality, there were only two haplotypes for each of the hin a, hin b1 and hin b2 genes and there was no clear relationship between grain hardness, malt or feed quality. The genotypes used in this study demonstrated that there was a low level of polymorphism in hardness genes in current commercial varieties as well as breeding lines and these polymorphisms had no impact on quality.

Short-term responses of leaf growth rate to water deficit scale up to whole-plant and crop levels: An integrated modelling approach in maize.

Physiological and genetic studies of leaf growth often focus on short-term responses, leaving a gap to whole-plant models that predict biomass accumulation, transpiration and yield at crop scale. To bridge this gap, we developed a model that combines an existing model of leaf 6 expansion in response to short-term environmental variations with a model coordinating the development of all leaves of a plant. The latter was based on: (1) rates of leaf initiation, appearance and end of elongation measured in field experiments; and (2) the hypothesis of an independence of the growth between leaves. The resulting whole-plant leaf model was integrated into the generic crop model APSIM which provided dynamic feedback of environmental conditions to the leaf model and allowed simulation of crop growth at canopy level. The model was tested in 12 field situations with contrasting temperature, evaporative demand and soil water status. In observed and simulated data, high evaporative demand reduced leaf area at the whole-plant level, and short water deficits affected only leaves developing during the stress, either visible or still hidden in the whorl. The model adequately simulated whole-plant profiles of leaf area with a single set of parameters that applied to the same hybrid in all experiments. It was also suitable to predict biomass accumulation and yield of a similar hybrid grown in different conditions. This model extends to field conditions existing knowledge of the environmental controls of leaf elongation, and can be used to simulate how their genetic controls flow through to yield.

Isolation and characterization of novel microsatellite markers and their application for diversity assessment in cultivated groundnut (Arachis hypogaea)

Background: Cultivated peanut or groundnut (Arachis hypogaea L.) is the fourth most important oilseed crop in the world, grown mainly in tropical, subtropical and warm temperate climates. Due to its origin through a single and recent polyploidization event, followed by successive selection during breeding efforts, cultivated groundnut has a limited genetic background. In such species, microsatellite or simple sequence repeat (SSR) markers are very informative and useful for breeding applications. The low level of polymorphism in cultivated germplasm, however, warrants a need of larger number of polymorphic microsatellite markers for cultivated groundnut. Results: A microsatellite- enriched library was constructed from the genotype TMV2. Sequencing of 720 putative SSR-positive clones from a total of 3,072 provided 490 SSRs. 71.2% of these SSRs were perfect type, 13.1% were imperfect and 15.7% were compound. Among these SSRs, the GT/CA repeat motifs were the most common (37.6%) followed by GA/CT repeat motifs (25.9%). The primer pairs could be designed for a total of 170 SSRs and were optimized initially on two genotypes. 104 (61.2%) primer pairs yielded scorable amplicon and 46 (44.2%) primers showed polymorphism among 32 cultivated groundnut genotypes. The polymorphic SSR markers detected 2 to 5 alleles with an average of 2.44 per locus. The polymorphic information content (PIC) value for these markers varied from 0.12 to 0.75 with an average of 0.46. Based on 112 alleles obtained by 46 markers, a phenogram was constructed to understand the relationships among the 32 genotypes. Majority of the genotypes representing subspecies hypogaea were grouped together in one cluster, while the genotypes belonging to subspecies fastigiata were grouped mainly under two clusters. Conclusion. Newly developed set of 104 markers extends the repertoire of SSR markers for cultivated groundnut. These markers showed a good level of PIC value in cultivated germplasm and therefore would be very useful for germplasm analysis, linkage mapping, diversity studies and phylogenetic relationships in cultivated groundnut as well as related Arachis species.

DArT markers: Diversity analyses and mapping in Sorghum bicolor

The sequential nature of gel-based marker systems entails low throughput and high costs per assay. Commonly used marker systems such as SSR and SNP are also dependent on sequence information. These limitations result in high cost per data point and significantly limit the capacity of breeding programs to obtain sufficient return on investment to justify the routine use of marker-assisted breeding for many traits and particularly quantitative traits. Diversity Arrays Technology (DArT™) is a cost effective hybridisation-based marker technology that offers a high multiplexing level while being independent of sequence information. This technology offers sorghum breeding programs an alternative approach to whole-genome profiling. We report on the development, application, mapping and utility of DArT™ markers for sorghum germplasm. Results: A genotyping array was developed representing approximately 12,000 genomic clones using PstI+BanII complexity with a subset of clones obtained through the suppression subtractive hybridisation (SSH) method. The genotyping array was used to analyse a diverse set of sorghum genotypes and screening a Recombinant Inbred Lines (RIL) mapping population. Over 500 markers detected variation among 90 accessions used in a diversity analysis. Cluster analysis discriminated well between all 90 genotypes. To confirm that the sorghum DArT markers behave in a Mendelian manner, we constructed a genetic linkage map for a cross between R931945-2-2 and IS 8525 integrating DArT and other marker types. In total, 596 markers could be placed on the integrated linkage map, which spanned 1431.6 cM. The genetic linkage map had an average marker density of 1/2.39 cM, with an average DArT marker density of 1/3.9 cM. Conclusion: We have successfully developed DArT markers for Sorghum bicolor and have demonstrated that DArT provides high quality markers that can be used for diversity analyses and to construct medium-density genetic linkage maps. The high number of DArT markers generated in a single assay not only provides a precise estimate of genetic relationships among genotypes, but also their even distribution over the genome offers real advantages for a range of molecular breeding and genomics applications.

Effects of fishing on tropical reef associated shark populations on the Great Barrier Reef

Three data sets were examined to define the level of interaction of reef associated sharks with the commercial Coral Reef Fin Fish Fishery within the Great Barrier Reef (GBR). Data were examined from fishery logbooks, an observer program within the fishery and a fishery-independent survey conducted as part of the Effects of Line Fishing (ELF) Experiment. The majority of the identified catch was comprised of grey reef (62-72%), whitetip reef (16-29%) and blacktip reef (6-13%) sharks. Logbook data revealed spatially and temporally variable landings of shark from the GBR. Catch per unit effort (CPUE) through time was stable for the period from 1989 to 2006 with no evidence of increase or decline. Data from observer and ELF data sets indicated no differences in CPUE among regions. The ELF data set demonstrated that CPUE was higher in Marine National Park zones (no fishing) when compared to General Use zones (open to fishing). The ongoing and consistent catches of reef sharks in the fishery and effectiveness of no-fishing zones suggest that management zones within the GBR Marine Park are effective at protecting a portion of the reef shark population from exploitation.

Annual forecasting of the Australian macadamia crop - integrating tree census data with statistical climate-adjustment models

To facilitate marketing and export, the Australian macadamia industry requires accurate crop forecasts. Each year, two levels of crop predictions are produced for this industry. The first is an overall longer-term forecast based on tree census data of growers in the Australian Macadamia Society (AMS). This data set currently accounts for around 70% of total production, and is supplemented by our best estimates of non-AMS orchards. Given these total tree numbers, average yields per tree are needed to complete the long-term forecasts. Yields from regional variety trials were initially used, but were found to be consistently higher than the average yields that growers were obtaining. Hence, a statistical model was developed using growers' historical yields, also taken from the AMS database. This model accounted for the effects of tree age, variety, year, region and tree spacing, and explained 65% of the total variation in the yield per tree data. The second level of crop prediction is an annual climate adjustment of these overall long-term estimates, taking into account the expected effects on production of the previous year's climate. This adjustment is based on relative historical yields, measured as the percentage deviance between expected and actual production. The dominant climatic variables are observed temperature, evaporation, solar radiation and modelled water stress. Initially, a number of alternate statistical models showed good agreement within the historical data, with jack-knife cross-validation R2 values of 96% or better. However, forecasts varied quite widely between these alternate models. Exploratory multivariate analyses and nearest-neighbour methods were used to investigate these differences. For 2001-2003, the overall forecasts were in the right direction (when compared with the long-term expected values), but were over-estimates. In 2004 the forecast was well under the observed production, and in 2005 the revised models produced a forecast within 5.1% of the actual production. Over the first five years of forecasting, the absolute deviance for the climate-adjustment models averaged 10.1%, just outside the targeted objective of 10%.

Responses in carcass lean pH at 24 hours post-mortem in lines of Large White pigs selected for growth rate on a fixed ration over a set time

A restricted maximum likelihood analysis applied to an animal model showed no significant differences (P > 0.05) in pH value of the longissimus dorsi measured at 24 h post-mortem (pH24) between high and low lines of Large White pigs selected over 4 years for post-weaning growth rate on restricted feeding. Genetic and phenotypic correlations between pH24 and production and carcass traits were estimated using all performance testing records combined with the pH24 measurements (5.05-7.02) on slaughtered animals. The estimate of heritability for pH24 was moderate (0.29 ± 0.18). Genetic correlations between pH24 and production or carcass composition traits, except for ultrasonic backfat (UBF), were not significantly different from zero. UBF had a moderate, positive genetic correlation with pH24 (0.24 ± 0.33). These estimates of genetic correlations affirmed that selection for increased growth rate on restricted feeding is likely to result in limited changes in pH24 and pork quality since the selection does not put a high emphasis on reduced fatness.

The effects of high stocking rates on milk production from dryland and irrigated Mediterranean pastures

An experiment using herds of similar to 20 cows (farmlets) assessed the effects of high stocking rates on production and profitability of feeding systems based on dryland and irrigated perennial ryegrass-based pastures in a Mediterranean environment in South Australia over 4 years. A target level of milk production of 7000 L/cow.year was set, based on predicted intakes of 2.7 t DM/cow.year as concentrates, pasture intakes from 1.5 to 2.7 t/cow.year and purchased fodder. In years 1 and 2, up to 1.5 t DM/cow.year of purchased fodder was used and in years 3 and 4 the amounts were increased if necessary to enable levels of milk production per cow to be maintained at target levels. Cows in dryland farmlets calved in March to May inclusive and were stocked at 2.5, 2.9, 3.3, 3.6 and 4.1 cows/ha, while those in irrigated farmlets calved in August to October inclusive and were stocked at 4.1, 5.2, 6.3 and 7.4 cows/ha. In the first 2 years, when inputs of purchased fodder were limited, milk production per cow was reduced with higher stocking rates (P < 0.01), but in years 3 and 4 there were no differences. Mean production was 7149 kg/cow.year in years 1 and 2, and 8162 kg/cow.year in years 3 and 4. Production per hectare was very closely related to stocking rate in all years (P < 0.01), increasing from 18 to 34 t milk/ha.year for dryland farmlets (1300 to 2200 kg milk solids/ha) and from 30 to 60 t milk/ha.year for irrigated farmlets (2200 to 4100 kg milk solids/ha). Almost all of these increases were attributed to the increases in grain and purchased fodder inputs associated with the increases in stocking rate. Net pasture accumulation rates and pasture harvest were generally not altered with stocking rate, though as stocking rate increased there was a change to more of the pasture being grazed and less conserved in both dryland and irrigated farmlets. Total pasture harvest averaged similar to 8 and 14 t DM/ha.year for dryland and irrigated pastures, respectively. An exception was at the highest stocking rate under irrigation, where pugging during winter was associated with a 14% reduction in annual pasture growth. There were several indications that these high stocking rates may not be sustainable without substantial changes in management practice. There were large and positive nutrient balances and associated increases in soil mineral content (P < 0.01), especially for phosphorus and nitrate nitrogen, with both stocking rate and succeeding years. Levels under irrigation were considerably higher (up to 90 and 240 mg/kg of soil for nitrate nitrogen and phosphorus, respectively) than under dryland pastures (60 and 140 mg/kg, respectively). Soil organic carbon levels did not change with stocking rate, indicating a high level of utilisation of forage grown. Weed ingress was also high (to 22% DM) in all treatments and especially in heavily stocked irrigated pastures during winter. It was concluded the higher stocking rates used exceeded those that are feasible for Mediterranean pastures in this environment and upper levels of stocking are suggested to be 2.5 cows/ha for dryland pastures and 5.2 cows/ha for irrigated pastures. To sustain these suggested stocking rates will require further development of management practices to avoid large increases in soil minerals and weed invasion of pastures.

Action learning in partnership with Landcare and catchment management groups to support increased pasture sowings in southern inland Queensland

The incorporation of sown pastures as short-term rotations into the cropping systems of northern Australia has been slow. The inherent chemical fertility and physical stability of the predominant vertisol soils across the region enabled farmers to grow crops for decades without nitrogen fertiliser, and precluded the evolution of a crop–pasture rotation culture. However, as less fertile and less physically stable soils were cropped for extended periods, farmers began to use contemporary farming and sown pasture technologies to rebuild and maintain their soils. This has typically involved sowing long-term grass and grass–legume pastures on the more marginal cropping soils of the region. In partnership with the catchment management authority, the Queensland Murray–Darling Committee (QMDC) and Landcare, a pasture extension process using the LeyGrain™ package was implemented in 2006 within two Grain & Graze projects in the Maranoa-Balonne and Border Rivers catchments in southern inland Queensland. The specific objectives were to increase the area sown to high quality pasture and to gain production and environmental benefits (particularly groundcover) through improving the skills of producers in pasture species selection, their understanding and management of risk during pasture establishment, and in managing pastures and the feed base better. The catalyst for increasing pasture sowings was a QMDC subsidy scheme for increasing groundcover on old cropping land. In recognising a need to enhance pasture knowledge and skills to implement this scheme, the QMDC and Landcare producer groups sought the involvement of, and set specific targets for, the LeyGrain workshop process. This is a highly interactive action learning process that built on the existing knowledge and skills of the producers. Thirty-four workshops were held with more than 200 producers in 26 existing groups and with private agronomists. An evaluation process assessed the impact of the workshops on the learning and skill development by participants, their commitment to practice change, and their future intent to sow pastures. The results across both project catchments were highly correlated. There was strong agreement by producers (>90%) that the workshops had improved knowledge and skills regarding the adaptation of pasture species to soils and climates, enabling a better selection at the paddock level. Additional strong impacts were in changing the attitudes of producers to all aspects of pasture establishment, and the relative species composition of mixtures. Producers made a strong commitment to practice change, particularly in managing pasture as a specialist crop at establishment to minimise risk, and in the better selection and management of improved pasture species (particularly legumes and the use of fertiliser). Producers have made a commitment to increase pasture sowings by 80% in the next 5 years, with fourteen producers in one group alone having committed to sow an additional 4893 ha of pasture in 2007–08 under the QMDC subsidy scheme. The success of the project was attributed to the partnership between QMDC and Landcare groups who set individual workshop targets with LeyGrain presenters, the interactive engagement processes within the workshops themselves, and the follow-up provided by the LeyGrain team for on-farm activities.

Two sympatric spotted gum species are molecularly homogeneous.

Large fruited spotted gum eucalypt Corymbia henryi occurs sympatrically with small fruited spotted gum Corymbia citriodora subspecies variegata over a large portion of its range on the east coast of Australia. The two taxa are interfertile, have overlapping flowering times and share a common set of insect and vertebrate pollinators. Previous genetic analysis of both taxa from two geographically remote sites suggested that the two were morphotypes rather than genetically distinct species. In this study we further explore this hypothesis of genic species by expanding sampling broadly through their sympatric locations and examine local-scale spatial genetic structure in stands that differ in species and age composition. Delineation of populations at five microsatellite loci, using an individual-based approach and Bayesian modelling, as well as clustering of individuals based on allele frequencies showed the two species to be molecularly homogeneous. Genetic structure aligned largely with geographic areas of origin, and followed an isolation-by-distance model, where proximal populations were generally less differentiated than more distant ones. At the stand level, spotted gums also generally showed little structure consistent with the high levels of gene flow inferred across the species range. Disturbances in the uniformity of structuring were detected, however, and attributed to localised events giving rise to even aged stands, probably due to regeneration from a few individuals following fire.

Refined Global Analysis of Bemisia tabaci (Hemiptera: Sternorrhyncha: Aleyrodoidea: Aleyrodidae) Mitochondrial Cytochrome Oxidase 1 to Identify Species Level Genetic Boundaries.

Identifying species boundaries within morphologically indistinguishable cryptic species complexes is often contentious. For the whitefly Bemisia tabaci (Gennadius) (Hemiptera: Sternorrhyncha: Aleyrodoidea: Aleyrodidae), the lack of a clear understanding about the genetic limits of the numerous genetic groups and biotypes so far identified has resulted in a lack of consistency in the application of the terms, the approaches use to apply them and in our understanding of what genetic structure within B. tabaci means. Our response has been to use mitochondrial gene cytochrome oxidase one to consider how to clearly and consistently define genetic separation. Using Bayesian phylogenetic analysis and analysis of sequence pairwise divergence we found a considerably higher to number of genetic groups than had been previously determined with two breaks in the distribution, one at 11% and another at 3.5%. At >11% divergence, 11 distinct groups were resolved, whereas at >3.5% divergence 24 groups were identified. Consensus sequences for each of these groups were determined and were shown to be useful in the correct assignment of sequences of unknown origin. The 3.5% divergence bound is consistent with species level separations in other insect taxa and Suggests that B. tabaci is it cryptic species composed of at least 24 distinct species. We further show that the placement of Bemesia atriplex (Froggatt) within the B. tabaci in, group adds further weight to the argument for species level separation within B. tabaci. This new analysis, which constructs consensus sequences and uses these its a standard against which unknown sequences call be compared, provides for the first time it consistent means of identifying the genetic hounds of each species with it high degree of certainty.

Development of a multi-locus sequence typing scheme for avian isolates of Pasteurella multocida.

A total of 63 isolates of Pasteurella multocida from Australian poultry, all associated with fowl cholera outbreaks, and three international reference strains, representing the three subspecies within P. multocida were used to develop a multi-locus sequence typing scheme. Primers were designed for conserved regions of seven house-keeping enzymes - adk, est, gdh, mdh, pgi, pmi and zwf - and internal fragments of 570-784 bp were sequenced for all isolates and strains. The number of alleles at the different loci ranged from 11 to 20 and a total of 29 allelic profiles or sequence types were recognised amongst the 66 strains. There was a strong concordance between the MLST data and the existing multi-locus enzyme electrophoresis and ribotyping data. When used to study a sub-set of isolates with a known detailed epidemiological history, the MLST data matched the results given by restriction endonuclease analysis, pulsed-field gel electrophoresis, ribotyping and REP-PCR. The MLST scheme provides a high level of resolution and is an excellent tool for studying the population structure and epidemiology of P. multocida.