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.

Sustained genetic control of wheat rust diseases in north-eastern Australia

Control of wheat rusts in north-eastern Australia has been based on resistance breeding since the early 1920s. It has been an enduring journey of discovery, disappointment, and achievement, which has culminated in a pool of knowledge and expertise upon which today's plant breeders can efficiently target durable resistance to the major rust diseases. This paper outlines significant advances in genetic control of rusts in the region, with particular emphasis on the invaluable role played by the University of Sydney rust control program and its influence on wheat breeding in the region and throughout Australia. This paper is part of ‘Global Landscapes in Cereal Rust Control’, see Aust. J. Agric. Res. Vol. 58, no. 6.

Genetics of heifer performance in 'wet' and 'dry' seasons and their relationships with steer performance in two tropical beef genotypes.

The genetics of heifer performance in tropical 'wet' and 'dry' seasons, and relationships with steer performance, were studied in Brahman (BRAH) and Tropical Composite (TCOMP) (50% Bos indicus, African Sanga or other tropically adapted Bos taurus; 50% non-tropically adapted Bos taurus) cattle of northern Australia. Data were from 2159 heifers (1027 BRAH, 1132 TCOMP), representing 54 BRAH and 51 TCOMP sires. Heifers were assessed after post-weaning 'wet' (ENDWET) and 'dry' (ENDDRY) seasons. Steers were assessed post-weaning, at feedlot entry, over a 70-day feed test, and after similar to 120-day finishing. Measures studied in both heifers and steers were liveweight (LWT), scanned rump fat, rib fat and M. longissimus area (SEMA), body condition score (CS), hip height (HH), serum insulin-like growth factor-I concentration (IGF-I), and average daily gains (ADG). Additional steer measures were scanned intra-muscular fat%, flight time, and daily (DFI) and residual feed intake (RFI). Uni- and bivariate analyses were conducted for combined genotypes and for individual genotypes. Genotype means were predicted for a subset of data involving 34 BRAH and 26 TCOMP sires. A meta-analysis of genetic correlation estimates examined how these were related to the difference between measurement environments for specific traits. There were genotype differences at the level of means, variances and genetic correlations. BRAH heifers were significantly (P < 0.05) faster-growing in the 'wet' season, slower-growing in the 'dry' season, lighter at ENDDRY, and taller and fatter with greater CS and IGF-I at both ENDWET and ENDDRY. Heritabilities were generally in the 20 to 60% range for both genotypes. Phenotypic and genetic variances, and genetic correlations, were commonly lower for BRAH. Differences were often explained by the long period of tropical adaptation of B. indicus. Genetic correlations were high between corresponding measures at ENDWET and ENDDRY, positive between fat and muscle measures in TCOMP but negative in BRAH (mean of 13 estimates 0.50 and -0.19, respectively), and approximately zero between steer feedlot ADG and heifer ADG in BRAH. Numerous genetic correlations between heifers and steers differed substantially from unity, especially in BRAH, suggesting there may be scope to select differently in the sexes where that would aid the differing roles of heifers and steers in production. Genetic correlations declined as measurement environments became more different, the rates of decline (environment sensitivity) sometimes differing with genotype. Similar measures (LWT, HH and ADG; IGF-I at ENDWET in TCOMP) were genetically correlated with steer DFI in heifers as in steers. Heifer SEMA was genetically correlated with steer feedlot RFI in BRAH (0.75 +/- 0.27 at ENDWET, 0.66 +/- 0.24 at ENDDRY). Selection to reduce steer RFI would reduce SEMA in BRAH heifers but otherwise have only small effects on heifers before their first joining.

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.

Genetic control of wheat quality: interactions between chromosomal regions determining protein content and composition, dough rheology, and sponge and dough baking properties

While the genetic control of wheat processing characteristics such as dough rheology is well understood, limited information is available concerning the genetic control of baking parameters, particularly sponge and dough (S&D) baking. In this study, a quantitative trait loci (QTL) analysis was performed using a population of doubled haploid lines derived from a cross between Australian cultivars Kukri x Janz grown at sites across different Australian wheat production zones (Queensland in 2001 and 2002 and Southern and Northern New South Wales in 2003) in order to examine the genetic control of protein content, protein expression, dough rheology and sponge and dough baking performance. The study highlighted the inconsistent genetic control of protein content across the test sites, with only two loci (3A and 7A) showing QTL at three of the five sites. Dough rheology QTL were highly consistent across the 5 sites, with major effects associated with the Glu-B1 and Glu-D1 loci. The Glu-D1 5 + 10 allele had consistent effects on S&D properties across sites; however, there was no evidence for a positive effect of the high dough strength Glu-B1-al allele at Glu-B1. A second locus on 5D had positive effects on S&D baking at three of five sites. This study demonstrated that dough rheology measurements were poor predictors of S&D quality. In the absence of robust predictive tests, high heritability values for S&D demonstrate that direct selection is the current best option for achieving genetic gain in this product category.

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.

Past, current and future approaches to mango genetic improvement in Australia

Improving the genetic base of cultivars that underpin commercial mango production is generally recognized as necessary for long term industry stability. Genetic improvement can take many approaches to improve cultivars, each with their own advantages and disadvantages. This paper will discuss several approaches used in the genetic improvement of mangoes in Australia, including varietal introductions, selection of monoembryonic progeny, selection within polyembryonic populations, assisted open pollination and controlled closed pollination. The current activities of the Australian National Mango Breeding Program will be outlined, and the analysis and use of hybrid phenotype data from the project for selection of next generation parents will be discussed. Some of the important traits that will enhance the competitiveness of future cultivars will be introduced and the challenges in achieving them discussed. The use of a genomics approach and its impact on future mango breeding is examined.

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.

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.

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.

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.

Negligible evidence for regional genetic population structure for two shark species Rhizoprionodon acutus (Rüppell, 1837) and Sphyrna lewini (Griffith & Smith, 1834) with contrasting biology

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

Migration of green turtles (Chelonia mydas) from Australasian feeding grounds inferred from genetic analyses

Coastal seagrass habitats in tropical and subtropical regions support aggregations of resident green turtles (Chelonia mydas) from several genetically distinct breeding populations. Migration of individuals to their respective dispersed breeding sites provides a complex pattern of migratory connectivity among nesting and feeding habitats of this species. An understanding of this pattern is important in regions where the persistence of populations is under threat from anthropogenic impacts. The present study uses mitochondrial DNA and mixed-stock analyses to assess the connectivity among seven feeding grounds across the north Australian coast and adjacent areas and 17 genetically distinct breeding populations from the Indo-Pacific region. It was hypothesised that large and geographically proximate breeding populations would dominate at nearby feeding grounds. As expected, each sampled feeding area appears to support multiple breeding populations, with two aggregations dominated by a local breeding population. Geographic distance between breeding and feeding habitat strongly influenced whether a breeding population contributed to a feeding ground (wi = 0.654); however, neither distance nor size of a breeding population was a good predictor of the extent of their contribution. The differential proportional contributions suggest the impact of anthropogenic mortality at feeding grounds should be assessed on a case-by-case basis.

Towards better management of Australia’s shark fishery: genetic analyses reveal unexpected ratios of cryptic blacktip species Carcharhinus tilstoni and C. limbatus

The common blacktip shark (Carcharhinus limbatus) and the Australian blacktip shark (C. tilstoni) are morphologically similar species that co-occur in subtropical and tropical Australia. In striking contrast to what has been previously reported, we demonstrate that the common blacktip shark is not rare in northern Australia but occurs in approximately equal frequencies with the Australian blacktip shark. Management of shark resources in northern Australia needs to take account of this new information. Species identification was performed using nucleotide sequences of the control, NADH dehydrogenase subunit 4 (ND4) and cytochrome oxidase I (COI) regions in the mitochondrial genome. The proportion of overall genetic variation (FST) between the two species was small (0.042, P < 0.01) based on allele frequencies at five microsatellite loci. We confirm that a third blacktip species (C. amblyrhynchoides, graceful shark) is closely related to C. tilstoni and C. limbatus and can be distinguished from them on the basis of mtDNA sequences from two gene regions. The Australian blacktip shark (C. tilstoni) was not encountered among 20 samples from central Indonesia that were later confirmed to be common blacktip and graceful sharks. Fisheries regulators urgently need new information on life history, population structure and morphological characters for species identification of blacktip shark species in Australia.