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.

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.

Identification of genetic regions associated with black point in barley.

Black point (BP) can cause severe losses to the barley industry through downgrading and discounting of malting barley. The genetic improvement in BP resistance of barley is complex, requiring reliable screening tools, an understanding of genotype by environment interactions and an understanding of the biochemical mechanisms of melanisation involved in BP development. Thus the application of molecular markers for resistance to BP may be a useful tool for plant breeders. We have investigated the genetic regions associated with BP resistance in the barley F2 population, Valier/Binalong. Quantitative trait loci (QTLs) contributed by the resistant parent Valier, were detected on chromosomes 2HS, 2HC, 3HL, 4HL and a QTL contributed by the susceptible parent, Binalong was detected on 5HL. Three of the four QTLs were detected in two distinctly different environments. The differences observed in BP resistance between these two environments and the implications for accelerated screening are discussed. Identified SSR markers in these regions may be useful for selecting black point resistance in related breeding materials.

Understanding invasion history: genetic structure and diversity of two globally invasive plants and implications for their management

Aim: Resolving the origin of invasive plant species is important for understanding the introduction histories of successful invaders and aiding strategies aimed at their management. This study aimed to infer the number and origin(s) of introduction for the globally invasive species, Macfadyena unguis-cati and Jatropha gossypiifolia using molecular data. Location: Native range: Neotropics; Invaded range: North America, Africa, Europe, Asia, Pacific Islands and Australia. Methods: We used chloroplast microsatellites (cpSSRs) to elucidate the origin(s) of introduced populations and calculated the genetic diversity in native and introduced regions. Results: Strong genetic structure was found within the native range of M. unguis-cati, but no genetic structuring was evident in the native range of J. gossypiifolia. Overall, 27 haplotypes were found in the native range of M. unguis-cati. Only four haplotypes were found in the introduced range, with more than 96% of introduced specimens matching a haplotype from Paraguay. In contrast, 15 haplotypes were found in the introduced range of J. gossypiifolia, with all invasive populations, except New Caledonia, comprising multiple haplotypes. Main conclusions: These data show that two invasive plant species from the same native range have had vastly different introduction histories in their non-native ranges. Invasive populations of M. unguis-cati probably came from a single or few independent introductions, whereas most invasive J. gossypiifolia populations arose from multiple introductions or alternatively from a representative sample of genetic diversity from a panmictic native range. As introduced M. unguis-cati populations are dominated by a single haplotype, locally adapted natural enemies should make the best control agents. However, invasive populations of J. gossypiifolia are genetically diverse and the selection of bio-control agents will be considerably more complex.

Quantitative and molecular genetics of juvenile wood traits in radiata and slash/Caribbean pines.

The Juvenile Wood Initiative (JWI) project has been running successfully since July 2003 under a Research Agreement with FWPA and Letters of Association with the consortium partners STBA (Southern Tree Breeding Association), ArborGen and FPQ (Forestry Plantations Queensland). Over the last five and half years, JWI scientists in CSIRO, FPQ, and STBA have completed all 12 major milestones and 28 component milestones according to the project schedule. We have made benchmark progress in understanding the genetic control of wood formation and interrelationships among wood traits. The project has made 15 primary scientific findings and several results have been adopted by industry as summarized below. This progress was detailed in 10 technical reports to funding organizations and industry clients. Team scientists produced 16 scientific manuscripts (8 published, 1 in press, 2 submitted, and several others in the process of submission) and 15 conference papers or presentations. Primary Scientific Findings. The 15 major scientific findings related to wood science, inheritance and the genetic basis of juvenile wood traits are: 1. An optimal method to predict stiffness of standing trees in slash/Caribbean pine is to combine gravimetric basic density from 12 mm increment cores with a standing tree prediction of MoE using a time of flight acoustic tool. This was the most accurate and cheapest way to rank trees for breeding selection for slash/Caribbean hybrid pine. This method was also recommended for radiata pine. 2. Wood density breeding values were predicted for the first time in the STBA breeding population using a large sample of 7,078 trees (increment cores) and it was estimated that selection of the best 250 trees for deployment will produce wood density gains of 12.4%. 3. Large genetic variation for a suite of wood quality traits including density, MFA, spiral grain, shrinkage, acoustic and non-acoustic stiffness (MoE) for clear wood and standing trees were observed. Genetic gains of between 8 and 49% were predicted for these wood quality traits with selection intensity between 1 to 10% for radiata pine. 4. Site had a major effect on juvenile-mature wood transition age and the effect of selective breeding for a shorter juvenile wood formation phase was only moderate (about 10% genetic gain with 10% selection intensity, equivalent to about 2 years reduction of juvenile wood). 5. The study found no usable site by genotype interactions for the wood quality traits of density, MFA and MoE for both radiata and slash/Caribbean pines, suggesting that assessment of wood properties on one or two sites will provide reliable estimates of the genetic worth of individuals for use in future breeding. 6. There were significant and sizable genotype by environment interactions between the mainland and Tasmanian regions and within Tasmania for DBH and branch size. 7. Strong genetic correlations between rings for density, MFA and MoE for both radiata and slash/Caribbean pines were observed. This suggests that selection for improved wood properties in the innermost rings would also result in improvement of wood properties in the subsequent rings, as well as improved average performance of the entire core. 8. Strong genetic correlations between pure species and hybrid performance for each of the wood quality traits were observed in the hybrid pines. Parental performance can be used to identify the hybrid families which are most likely to have superior juvenile wood properties of the slash/Caribbean F1 hybrid in southeast Queensland. 9. Large unfavourable genetic correlations between growth and wood quality traits were a prominent feature in radiata pine, indicating that overcoming this unfavourable genetic correlation will be a major technical issue in progressing radiata pine breeding. 10. The project created the first radiata pine 18 k cDNA microarray and generated 5,952 radiata pine xylogenesis expressed sequence tags (ESTs) which assembled into 3,304 unigenes. 11. A total of 348 genes were identified as preferentially expressed genes in earlywood or latewood while a total of 168 genes were identified as preferentially expressed genes in either juvenile or mature wood. 12. Juvenile earlywood has a distinct transcriptome relative to other stages of wood development. 13. Discovered rapid decay of linkage disequilibrium (LD) in radiata pine with LD decaying to approximately 50% within 1,700 base pairs (within a typical gene). A total of 913 SNPS from sequencing 177,380 base pairs were identified for association genetic studies. 14. 149 SNPs from 44 genes and 255 SNPs from a further 51 genes (total 95 genes) were selected for association analysis with 62 wood traits, and 30 SNPs were shortlisted for their significant association with variation of wood quality traits (density, MFA and MoE) with individual significant SNPs accounting for between 1.9 and 9.7% of the total genetic variation in traits. 15. Index selection using breeding objectives was the most profitable selection method for radiata pine, but in the long term it may not be the most effective in dealing with negative genetic correlations between wood volume and quality traits. A combination of economic and biological approaches may be needed to deal with the strong adverse correlation.

Temporal genetic structure patterns in tropical maize populations under reciprocal recurrent selection.

In order to investigate the effect of long term recurrent selection on the pattern of gene diversity, thirty randomly-selected individuals from the progenitors (p) and four selection cycles (C0, C3, C6 and C11) were sampled for DNA analysis from the tropical maize (Zea mays L.) breeding populations, Atherton 1 (AT1) and Atherton 2 (AT2). Fifteen polymorphic Simple Sequence Repeat markers amplified a total of 284 and 257 alleles in AT1 and AT2 populations, respectively. Reductions in the number of alleles were observed at advanced selection cycles. About 11 and 12% of the alleles in AT1 and AT2 populations respectively, were near to fixation. However, a higher number of alleles (37% in AT1 and 33% in AT2) were close to extinction. Fisher's exact test and analysis of molecular variance (AMOVA) showed significant population differentiations. Gene diversity estimates and AMOVA revealed increased genetic differentiations at the expense of loss of heterozygosity. Population differentiations were mainly due to fixation of complementary alleles at a locus in the two breeding populations. The estimates of effective population at an advanced selection cycle were close to the population size predicted by the breeding method.

IUCN classification zones concord with, but underestimate, the population genetic structure of the zebra shark Stegostoma fasciatum in the Indo-West Pacific

The Indo-West Pacific (IWP), from South Africa in the western Indian Ocean to the western Pacific Ocean, contains some of the most biologically diverse marine habitats on earth, including the greatest biodiversity of chondrichthyan fishes. The region encompasses various densities of human habitation leading to contrasts in the levels of exploitation experienced by chondrichthyans, which are targeted for local consumption and export. The demersal chondrichthyan, the zebra shark, Stegostoma fasciatum, is endemic to the IWP and has two current regional International Union for the Conservation of Nature (IUCN) Red List classifications that reflect differing levels of exploitation: ‘Least Concern’ and ‘Vulnerable’. In this study, we employed mitochondrial ND4 sequence data and 13 microsatellite loci to investigate the population genetic structure of 180 zebra sharks from 13 locations throughout the IWP to test the concordance of IUCN zones with demographic units that have conservation value. Mitochondrial and microsatellite data sets from samples collected throughout northern Australia and Southeast Asia concord with the regional IUCN classifications. However, we found evidence of genetic subdivision within these regions, including subdivision between locations connected by habitat suitable for migration. Furthermore, parametric FST analyses and Bayesian clustering analyses indicated that the primary genetic break within the IWP is not represented by the IUCN classifications but rather is congruent with the Indonesian throughflow current. Our findings indicate that recruitment to areas of high exploitation from nearby healthy populations in zebra sharks is likely to be minimal, and that severe localized depletions are predicted to occur in zebra shark populations throughout the IWP region.

Genetic Options for nematode control

Root-lesion nematodes (Pratylenchus thornei and P. neglectus) cause severe economic loss in wheat in Australia. This project aims to develop adaptaed wheat lines with resistance and tolerance to both species. These lines will be made available to Australian wheat breeding companies for further crossing and development of resistant and tolerant wheat varieties. Sources of resistance will be synthetic hexaploid and landrace wheats from the Middle East. Suitable double haploid populations will be phenotyped for the development of molecular markers to resistance and tolerance genes. The value of resistance and tolerance will be extended to growers through collaboration in demonstration trials with NGA and ORANA and presentations at GRDC Updates.

Rapid microsatellite marker development for African mahogany ( Khaya senegalensis, Meliaceae) using next-generation sequencing and assessment of its intra-specific genetic diversity.

Khaya senegalensis (African mahogany or dry-zone mahogany) is a high-value hardwood timber species with great potential for forest plantations in northern Australia. The species is distributed across the sub-Saharan belt from Senegal to Sudan and Uganda. Because of heavy exploitation and constraints on natural regeneration and sustainable planting, it is now classified as a vulnerable species. Here, we describe the development of microsatellite markers for K. senegalensis using next-generation sequencing to assess its intra-specific diversity across its natural range, which is a key for successful breeding programs and effective conservation management of the species. Next-generation sequencing yielded 93943 sequences with an average read length of 234bp. The assembled sequences contained 1030 simple sequence repeats, with primers designed for 522 microsatellite loci. Twenty-one microsatellite loci were tested with 11 showing reliable amplification and polymorphism in K. senegalensis. The 11 novel microsatellites, together with one previously published, were used to assess 73 accessions belonging to the Australian K. senegalensis domestication program, sampled from across the natural range of the species. STRUCTURE analysis shows two major clusters, one comprising mainly accessions from west Africa (Senegal to Benin) and the second based in the far eastern limits of the range in Sudan and Uganda. Higher levels of genetic diversity were found in material from western Africa. This suggests that new seed collections from this region may yield more diverse genotypes than those originating from Sudan and Uganda in eastern Africa.

Genetic control of nodal root angle in sorghum and its implications on water extraction

Genotypic variability in root system architecture has been associated with root angle of seedlings and water extraction patterns of mature plants in a range of crops. The potential inclusion of root angle as a selection criterion in a sorghum breeding program requires (1) availability of an efficient screening method, (2) presence of genotypic variation with high heritability, and (3) an association with water extraction pattern. The aim of this study was to determine the feasibility for inclusion of nodal root angle as a selection criterion in sorghum breeding programs. A high-throughput phenotypic screen for nodal root angle in young sorghum plants has recently been developed and has been used successfully to identify significant variation in nodal root angle across a diverse range of inbred lines and a mapping population. In both cases, heritabilities for nodal root angle were high. No association between nodal root angle and plant size was detected. This implies that parental inbred lines could potentially be used to asses nodal root angle of their hybrids, although such predictability is compromised by significant interactions. To study effects of nodal root angle on water extraction patterns of mature plants, four inbred lines with contrasting nodal root angle at seedling stage were grown until at least anthesis in large rhizotrons. A consistent trend was observed that nodal root angle may affect the spatial distribution of root mass of mature plants and hence their ability to extract soil water, although genotypic differences were not significant. The potential implications of this for specific adaptation to drought stress are discussed. Results suggest that nodal root angle of young plants can be a useful selection criterion for specific drought adaptation, and could potentially be used in molecular breeding programs if QTLs for root angle can be identified. (C) 2012 Elsevier B.V. All rights reserved.

A review of the application of molecular genetics for fisheries management and conservation of sharks and rays

Since the first investigation 25 years ago, the application of genetic tools to address ecological and evolutionary questions in elasmobranch studies has greatly expanded. Major developments in genetic theory as well as in the availability, cost effectiveness and resolution of genetic markers were instrumental for particularly rapid progress over the last 10 years. Genetic studies of elasmobranchs are of direct importance and have application to fisheries management and conservation issues such as the definition of management units and identification of species from fins. In the future, increased application of the most recent and emerging technologies will enable accelerated genetic data production and the development of new markers at reduced costs, paving the way for a paradigm shift from gene to genome-scale research, and more focus on adaptive rather than just neutral variation. Current literature is reviewed in six fields of elasmobranch molecular genetics relevant to fisheries and conservation management (species identification, phylogeography, philopatry, genetic effective population size, molecular evolutionary rate and emerging methods). Where possible, examples from the Indo-Pacific region, which has been underrepresented in previous reviews, are emphasized within a global perspective. (C) 2012 The Authors Journal of Fish Biology (C) 2012 The Fisheries Society of the British Isles

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.

New SNPs for population genetic analysis reveal possible cryptic speciation of eastern Australian sea mullet (Mugil cephalus)

Sustainable management of sea mullet (Mugil cephalus) fisheries needs to account for recent observations of regional-scale differentiation. Population genetic analysis is sought to assess the situation of this ecologically and economically important fish species in eastern Australian waters. Here, we report (i) new population genetic markers [single nucleotide polymorphisms (SNPs) and potential microsatellites], (ii) first estimates of spatial genetic differentiation and (iii) prospective power tests for designing more comprehensive studies. Six DNA samples from three sampling regions (North Queensland, South Queensland and central New South Wales) on the eastern coast of Australia were used to prepare restriction site associated DNA (RAD) tag libraries from genomic DNA digested with EcoRI and MseI. A pooled sample of regional RAD tag libraries was sequenced using the Roche GS-FLX Titanium platform. A total of 172837 raw reads (17.4Mbp) were retrieved, 95500 of which were used to discover 1267 SNPs and 1417 microsatellites. A subset of 161 SNPs was validated based on 63 additional DNA samples genotyped using the Sequenom MassArray (iPLEX Gold chemistry). Altogether 92 SNPs (57%) were confirmed, with 40% of these marking fixed variants between northern and southern sampling regions. Our preliminary findings indicate a multispecies fishery stock of M. cephalus in eastern Australian waters, but suggest that strong genetic differentiation occurs north of major fishing grounds. Low potential differentiation within major fishing grounds (e.g. FST=0.0025) can be resolved with a likely power 67% by using standard sample sizes of 50 and validated subsets of available markers.

Spot form of net blotch resistance in barley is under complex genetic control

Key message: Evaluation of resistance toPyrenophora teresf.maculatain barley breeding populations via association mapping revealed a complex genetic architecture comprising a mixture of major and minor effect genes. Abstract: In the search for stable resistance to spot form of net blotch (Pyrenophora teres f. maculata, SFNB), association mapping was conducted on four independent barley (Hordeum vulgare L.) breeding populations comprising a total of 898 unique elite breeding lines from the Northern Region Barley Breeding Program in Australia for discovery of quantitative trait loci (QTL) influencing resistance at seedling and adult plant growth stages. A total of 29 significant QTL were validated across multiple breeding populations, with 22 conferring resistance at both seedling and adult plant growth stages. The remaining 7 QTL conferred resistance at either seedling (2 QTL) or adult plant (5 QTL) growth stages only. These 29 QTL represented 24 unique genomic regions, of which five were found to co-locate with previously identified QTL for SFNB. The results indicated that SFNB resistance is controlled by a large number of QTL varying in effect size with large effects QTL on chromosome 7H. A large proportion of the QTL acted in the same direction for both seedling and adult responses, suggesting that phenotypic selection for SFNB resistance performed at either growth stage could achieve adequate levels of resistance. However, the accumulation of specific resistance alleles on several chromosomes must be considered in molecular breeding selection strategies.

A molecular survey of Eimeria in chickens across Australia

Coccidiosis is a costly enteric disease of chickens caused by protozoan parasites of the genus Eimeria. Disease diagnosis and management is complicated since there are multiple Eimeria species infecting chickens and mixed species infections are common. Current control measures are only partially effective and this, combined with concerns over vaccine efficacy and increasing drug resistance, demonstrates a need for improved coccidiosis diagnosis and control. Before improvements can be made, it is important to understand the species commonly infecting poultry flocks in both backyard and commercial enterprises. The aim of this project was to conduct a survey and assessment of poultry Eimeria across Australia using genetic markers, and create a collection of isolates for each Eimeria species. A total of 260 samples (faecal or caecal) was obtained, and survey results showed that Eimeria taxa were present in 98% of commercial and 81% of backyard flocks. The distribution of each Eimeria species was widespread across Australia, with representatives of all species being found in every state and territory, and the Eimeria species predominating in commercial flocks differed from those in backyard flocks. Three operational taxonomic units also occurred frequently in commercial flocks highlighting the need to understand the impact of these uncharacterised species on poultry production. As Eimeria infections were also frequent in backyard flocks, there is a potential for backyard flocks to act as reservoirs for disease, especially as the industry moves towards free range production systems. This Eimeria collection will be an important genetic resource which is the crucial first step in the development of more sophisticated diagnostic tools and the development of new live vaccines which ultimately will provide savings to the industry in terms of more efficient coccidiosis management.