Laboratory information management software for genotyping workflows: Applications in high throughput crop genotyping

Background: With the advances in DNA sequencer-based technologies, it has become possible to automate several steps of the genotyping process leading to increased throughput. To efficiently handle the large amounts of genotypic data generated and help with quality control, there is a strong need for a software system that can help with the tracking of samples and capture and management of data at different steps of the process. Such systems, while serving to manage the workflow precisely, also encourage good laboratory practice by standardizing protocols, recording and annotating data from every step of the workflow Results: A laboratory information management system (LIMS) has been designed and implemented at the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) that meets the requirements of a moderately high throughput molecular genotyping facility. The application is designed as modules and is simple to learn and use. The application leads the user through each step of the process from starting an experiment to the storing of output data from the genotype detection step with auto-binning of alleles; thus ensuring that every DNA sample is handled in an identical manner and all the necessary data are captured. The application keeps track of DNA samples and generated data. Data entry into the system is through the use of forms for file uploads. The LIMS provides functions to trace back to the electrophoresis gel files or sample source for any genotypic data and for repeating experiments. The LIMS is being presently used for the capture of high throughput SSR (simple-sequence repeat) genotyping data from the legume (chickpea, groundnut and pigeonpea) and cereal (sorghum and millets) crops of importance in the semi-arid tropics. Conclusions: A laboratory information management system is available that has been found useful in the management of microsatellite genotype data in a moderately high throughput genotyping laboratory. The application with source code is freely available for academic users and can be downloaded from http://www.icrisat.org/bt-software-d-lims.htm

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.

Use of Petrifilm(R) 3M to assess coliform numbers on lamb carcasses.

Petrifilm(R) (6410) was used directly on lamb carcasses to enumerate coliforms. 10 sites on 30 carcasses were sampled at each of 4 separate meat processing establishments (works). Coliform counts obtained by this technique were statistically analysed using analysis of variance (ANOVA) to select the optimum sampling sites on the carcass and to assess contamination of the carcass by gut flora at a particular establishment. There was a large variation between sites and between works. In general, works 3 and 4 produced cleaner carcasses than works 2, which in turn was cleaner than works 1. Works 1, 2 and 4 used conventional dressing techniques and works 3 used the inverted dressing method, therefore, the coliform counts found at works 3 and 4 are achievable regardless of dressing technique. Coliform bacteria were most concentrated around the posterior pelvic rim and less prevalent at the carcass extremities. The posterior pelvic rim (sites 3 and 4) had higher (P < 0.05) coliform counts than the exterior ventral flank area (sites 5, 6, 7 and 8), which in turn had higher (P < 0.05) counts than the proximal hind and proximal fore limbs (sites 1, 2, 9 and 10) across all works. With in-line routine testing it is recommended that the majority of carcasses sampled should give coliform counts of <50 cfu/20 cm2 for sites 4 and 8. Reprinted with permission from Journal of Food Protection. Copyright held by the International Association of Food Protection, Des Moines, Iowa, USA. Authors affifiation. J.A.Guthrie & K.J.Dunlop International Food Institute of Queensland, Department of Primary Industries, Rockhampton and G.A.Saunders Veterinary Public Health Division, Livestock and Meat Authority of Queensland, Emerald.

Effect of microwave radiation on seed mortality of rubber vine (Cryptostegia grandiflora R.Br.), parthenium (Parthenium hysterophorous L.) and bellyache bush (Jatropha gossypiifolia L.)

A trial was undertaken to evaluate the effect of microwaves on seed mortality of three weed species. Seeds of rubber vine (Cryptostegia grandiflora R.Br.), parthenium (Parthenium hysterophorous L.) and bellyache bush (Jatropha gossypiifolia L.) were buried at six depths (0, 2.5, 5, 10, 20 and 40 cm) in coarse sand maintained at one of two moisture levels, oven dry or wet (field capacity), and then subjected to one of five microwave radiation durations of (0, 2, 4, 8 and 16 min). Significant interactions between soil moisture level, microwave radiation duration, seed burial depth and species were detected for mortality of seeds of all three species. Maximum seed mortality of rubber vine (88%), parthenium (67%) and bellyache bush (94%) occurred in wet soil irradiated for 16 min. Maximum seed mortality of rubber vine and bellyache bush seeds occurred in seeds buried at 2.5 cm depth whereas that of parthenium occurred in seeds buried at 10 cm depth. Maximum soil temperatures of 114.1 and 87.5°C in dry and wet soil respectively occurred at 2.5 cm depth following 16 min irradiation. Irrespective of the greater soil temperatures recorded in dry soil, irradiating seeds in wet soil generally increased seed mortality 2.9-fold compared with dry soil. Moisture content of wet soil averaged 5.7% compared with 0.1% for dry soil. Results suggest that microwave radiation has the potential to kill seeds located in the soil seed bank. However, many factors, including weed species susceptibility, determine the effectiveness of microwave radiation on buried seeds. Microwave radiation may be an alternative to conventional methods at rapidly depleting soil seed banks in the field, particularly in relatively wet soils that contain long lived weed seeds.

Predicting wood value in Queensland Caribbean pine plantations using a decision support system

A decision support system has been developed in Queensland to evaluate how changes in silvicultural regimes affect wood quality, and specifically the graded recovery of structural timber. Models of tree growth, branch architecture and wood properties were developed from data collected in routine Caribbean pine plantations and specific silvicultural experiments. These models were incorporated in software that simulates the conversion of standing trees into logs, and the logs into boards, and generates detailed data on knot location and basic density distribution. The structural grade of each board was determined by simulating the machine stress-grading process, and the predicted graded recovery provided an indicator of wood value. The decision support system improves the basis of decision-making by simulating the performance of elite genetic material under specified silvicultural regimes and by predicting links between wood quality and general stand attributes such as stocking and length of rotation.

Poisoning for production: how effective is fox baiting in south-eastern Australia?

1. The European red fox Vulpes vulpes represents a continuing threat to both livestock and native vertebrates in Australia, and is commonly managed by setting ground-level baits impregnated with 1080 (sodium fluoroacetate) poison. However, the long-term effectiveness of such control campaigns is likely to be limited due to the ability of foxes to disperse over considerable distances and to swiftly recolonize areas from where they had been removed. 2. To investigate the effectiveness of fox baiting in a production landscape, we assessed the potential for foxes to reinvade baited farm property areas within the jurisdiction of the Molong Rural Lands Protection Board (RLPB), an area of 815 000 ha on the central tablelands of New South Wales, Australia. The spatial distribution and timing of fox baiting campaigns between 1998 and 2002 was estimated from RLPB records and mapped using Geographical Information System software. The effectiveness of the control campaign was assessed on the basis of the likely immigration of foxes from non-baited farms using immigration distances calculated from published relationships between dispersal distance and home range size. 3. Few landholders undertook baiting campaigns in any given year, and the area baited was always so small that no baited property would have been sufficiently far from an unbaited property to have been immune from immigrating individuals. It is likely, therefore, that immigration onto farms negated any long-term effects of baiting operations. This study highlights some of the key deficiencies in current baiting practices in south-eastern Australia and suggests that pest management programmes should be monitored using such methods to ensure they achieve their goals.

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.

Using a supply chain approach to guide R&D

Consumers today are presented with an increasing array of products. The growing competition for consumer expenditure requires a whole of supply chain approach to maintain market share for existing cultivars and to successfully commercialise new cultivars. The supply chain needs to deliver value and satisfaction to the end customer and profitability to their members. Critical to getting the product right is developing inherent robustness into the cultivar, and developing processes and systems through the whole supply chain that maintain product quality and add value. This paper describes the approach we have used in working with supply chains in Australia and Indonesia to identify priority areas for improvement. Our experience demonstrates the need for a champion in the supply chain with significant influence and a desire to improve. The paper also describes our approach towards improving a specific supply chain to achieve successful commercialisation of a new cultivar. The cultivar was primarily selected for good production characteristics and attractive visual appeal. The performance of the fruit is being monitored from farm to retail shelf to identify points where quality is lost and practices can be improved. A targeted R&D program is investigating ways of improving production efficiency (nutrition, flowering and canopy management), maturity standards to optimise flavour, harvesting and packing practices to reduce skin damage, and ripening and handling practices to optimise shelf life. This integrated approach is based on similar approaches used to improve the performance of existing mango and avocado cultivars.

ShadowBoxer and LocusEater: programs to optimise experimental design and multiplexing strategies for genetic mark-recapture

Genetic mark–recapture requires efficient methods of uniquely identifying individuals. 'Shadows' (individuals with the same genotype at the selected loci) become more likely with increasing sample size, and bias harvest rate estimates. Finding loci is costly, but better loci reduce analysis costs and improve power. Optimal microsatellite panels minimize shadows, but panel design is a complex optimization process. locuseater and shadowboxer permit power and cost analysis of this process and automate some aspects, by simulating the entire experiment from panel design to harvest rate estimation.

How accurate are the marker orders in crop linkage maps generated from large marker datasets?

Marker ordering during linkage map construction is a critical component of QTL mapping research. In recent years, high-throughput genotyping methods have become widely used, and these methods may generate hundreds of markers for a single mapping population. This poses problems for linkage analysis software because the number of possible marker orders increases exponentially as the number of markers increases. In this paper, we tested the accuracy of linkage analyses on simulated recombinant inbred line data using the commonly used Map Manager QTX (Manly et al. 2001: Mammalian Genome 12, 930-932) software and RECORD (Van Os et al. 2005: Theoretical and Applied Genetics 112, 30-40). Accuracy was measured by calculating two scores: % correct marker positions, and a novel, weighted rank-based score derived from the sum of absolute values of true minus observed marker ranks divided by the total number of markers. The accuracy of maps generated using Map Manager QTX was considerably lower than those generated using RECORD. Differences in linkage maps were often observed when marker ordering was performed several times using the identical dataset. In order to test the effect of reducing marker numbers on the stability of marker order, we pruned marker datasets focusing on regions consisting of tightly linked clusters of markers, which included redundant markers. Marker pruning improved the accuracy and stability of linkage maps because a single unambiguous marker order was produced that was consistent across replications of analysis. Marker pruning was also applied to a real barley mapping population and QTL analysis was performed using different map versions produced by the different programs. While some QTLs were identified with both map versions, there were large differences in QTL mapping results. Differences included maximum LOD and R-2 values at QTL peaks and map positions, thus highlighting the importance of marker order for QTL mapping

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

Background: Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers. Results: The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci ( 1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions. Conclusion: The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

Applications of CLIMEX modelling leading to improved biological control

Wilmot Senaratne, Bill Palmer and Bob Sutherst recently published their paper 'Applications of CLIMEX modelling leading to improved biological control' in Proceedings of the 16th Australian Weeds Conference. They looked at three examples where modern climate matching techniques using computer software produces decisions and results than might happen using previous techniques such as climadiagrams. Assessment of climatic suitability is important at various stages of a biological control project; from initial foreign exploration, to risk assessment in preparation for the release of a particular agent, through to selection of release sites that maximise the agent´s chances of initial establishment. It is now also necessary to predict potential future distributions of both target weeds and agents under climate change.

Estimation in a multiplicative mixed model involving a genetic relationship matrix.

Genetic models partitioning additive and non-additive genetic effects for populations tested in replicated multi-environment trials (METs) in a plant breeding program have recently been presented in the literature. For these data, the variance model involves the direct product of a large numerator relationship matrix A, and a complex structure for the genotype by environment interaction effects, generally of a factor analytic (FA) form. With MET data, we expect a high correlation in genotype rankings between environments, leading to non-positive definite covariance matrices. Estimation methods for reduced rank models have been derived for the FA formulation with independent genotypes, and we employ these estimation methods for the more complex case involving the numerator relationship matrix. We examine the performance of differing genetic models for MET data with an embedded pedigree structure, and consider the magnitude of the non-additive variance. The capacity of existing software packages to fit these complex models is largely due to the use of the sparse matrix methodology and the average information algorithm. Here, we present an extension to the standard formulation necessary for estimation with a factor analytic structure across multiple environments.

Maximising the contribution of native-range studies towards the identification and prioritisation of weed biocontrol agents

Effective study in the native range to identify potential agents underpins all efforts in classical biological control of weeds. Good agents that demonstrate both a high degree of host specificity and the potential to be damaging are a very limited resource and must therefore be carefully studied and considered. The overseas component is often operationally difficult and expensive but can contribute considerably more than a list of herbivores attacking a particular target. While the principles underlying this foreign component have been understood for some time, recently developed technologies and methods can make very significant contributions to foreign studies. Molecular and genetic characterisations of both target weed and agent organism can be increasingly employed to more accurately define the identity and phylogeny of them. Climate matching and modelling software is now available and can be utilised to better select agents for particular regions of concern. Relational databases can store collection information for analysis and future enquiry while quantification of sampling effort, employment of statistical survey methods and analysis by techniques such as rarefaction curves contribute to efficient and effective searching. Obtaining good and timely identifications for discovered agent organisms is perhaps the most serious issue confronting the modern explorer. The diminishing numbers of specialist taxonomists employed at the major museums while international and national protocols demand higher standards of identity exacerbates the issue. Genetic barcoding may provide a very useful tool to overcome this problem. Native-range work also offers under-exploited opportunities for contributing towards predicting safety, abundance and efficacy of potential agents in their target environment.

Use of CLIMEX modelling to identify prospective areas for exploration to find new biological control agents for prickly acacia

Selection of biocontrol agents that are adapted to the climates in areas of intended release demands a thorough analysis of the climates of the source and release sites. We present a case study that demonstrates how use of the CLIMEX software can improve decision making in relation to the identification of prospective areas for exploration for agents to control the woody weed, prickly acacia Acacia nilotica ssp. indica in the arid areas of north Queensland.