A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

Genetic diversity, aggressiveness and metalaxyl sensitivity of Pythium aphanidermatum populations infecting cucumber in Oman

Seventy three isolates of Pythium aphanidermatum obtained from cucumber from four different regions of Oman and 16 isolates of muskmelon from the Batinah region in Oman were characterized for aggressiveness, sensitivity to metalaxyl and genetic diversity using AFLP fingerprinting. Twenty isolates of P. aphanidermatum from diverse hosts from different countries were also included in the study. Most isolates from Oman were found to be aggressive on cucumber seedlings and all were highly sensitive to metalaxyl (EC50 < 0•80 µg mL−1). Isolates from cucumber and muskmelon were as aggressive as each other on both hosts (P > 0.05), which implies a lack of host specialization in P. aphanidermatum on these two hosts in Oman. AFLP analysis of all isolates using four primer-pair combinations resolved 152 bands, of which 61 (~40%) were polymorphic. Isolates of P. aphanidermatum from Oman and other countries exhibited high genetic similarity (mean = 94.1%) and produced 59 different AFLP profiles. Analysis of molecular variance indicated that most AFLP variation among populations of P. aphanidermatum in Oman was associated with geographical regions (FST = 0.118; P < 0.0001), not hosts (FST = -0.004; P = 0.4323). These data were supported by the high rate of recovery (24%) of identical phenotypes between cucumber and muskmelon fields in the same region as compared to the low recovery (10%) across regions in Oman, which suggests more frequent movement of Pythium inoculum among muskmelon and cucumber fields in the same region compared to movement across geographically separated regions. However, recovering clones among regions and different countries may imply circulation of Pythium inoculum via common sources in Oman and also intercontinental spread of isolates.

Genetic diversity, aggressiveness and metalaxyl sensitivity of Pythium spinosum infecting cucumber in Oman

A total of 24 isolates of Pythium spinosum from cucumber obtained from five regions in Oman were characterized for genetic diversity using amplified fragment length polymorphism (AFLP) fingerprinting and three isolates from the Netherlands, South Africa and Japan were included for comparison. Isolates from Oman were also characterized for aggressiveness on cucumber seedlings and sensitivity to metalaxyl. Identity of all isolates was confirmed using sequences of the internal transcribed spacer (ITS) region of the ribosomal DNA (rDNA), which showed more than 99% nucleotide similarity among all isolates. Using six primer-pair combinations, AFLP fingerprinting resolved 295 AFLP markers of which 193 were polymorphic among isolates from other countries and only six were polymorphic among isolates of P. spinosum from Oman. Seven different AFLP phenotypes of P. spinosum were recovered in Oman; two of them were found to contain over 79% of isolates and one was recovered from all regions in Oman. Phenotypes from Oman showed very high (?99%) levels of genetic similarity to each other compared to moderate (mean =53%) levels of genetic similarity with phenotypes from other countries. In addition, all isolates from Oman were found to be highly sensitive to metalaxyl and all were aggressive on cucumber seedlings at 25°C. The high genetic similarity among phenotypes of P. spinosum in Oman as well as recovering two major clones across regions may suggest that P. spinosum has been recently introduced in Oman via a common source.

Transcriptome analysis of Neotyphodium and Epichloë grass endophytes

Large-scale gene discovery has been performed for the grass fungal endophytes Neotyphodium coenophialum, Neotyphodium lolii, and Epichloë festucae. The resulting sequences have been annotated by comparison with public DNA and protein sequence databases and using intermediate gene ontology annotation tools. Endophyte sequences have also been analysed for the presence of simple sequence repeat and single nucleotide polymorphism molecular genetic markers. Sequences and annotation are maintained within a MySQL database that may be queried using a custom web interface. Two cDNA-based microarrays have been generated from this genome resource. They permit the interrogation of 3806 Neotyphodium genes (NchipTM microarray), and 4195 Neotyphodium and 920 Epichloë genes (EndoChipTM microarray), respectively. These microarrays provide tools for high-throughput transcriptome analysis, including genome-specific gene expression studies, profiling of novel endophyte genes, and investigation of the host grass–symbiont interaction. Comparative transcriptome analysis in Neotyphodium and Epichloë was performed

Assessing for genetic and environmental effects on ruminant feed quality in barley ( Hordeum vulgare )

Grain samples from a combined intermediate and advanced stage barley breeding trial series, grown at two sites in two consecutive years were assessed for detailed grain quality and ruminant feed quality. The results indicated that there were significant genetic and environmental effects for “feed” traits as measured using grain hardness, acid detergent fibre (ADF), starch and in-sacco dry matter digestibility (ISDMD) assays. In addition, there was strong genotypic discrimination for the regressed feed performance traits, namely Net Energy (NE) and Average Daily Gain (ADG). There was considerable variation in genetic correlations for all traits based on variance from the cultivars used, sites or laboratory processing effects. There was a high level of heritability ranging from 89% to 88% for retention, 60% to 80% for protein and 56% to 68% for ADF. However, there were only low to moderate levels of heritability for the feed traits, with starch 30–39%, ISDMD 55–63%, ADF 56–68%, particle size 47–73%, 31–48% NE and ADG 44–51%. These results suggest that there were real differences in the feed performance of barleys and that selection for cattle feed quality is potentially a viable option for breeding programs.

Applications of pedigree-based genome mapping in wheat and barley breeding programs.

The aim of the pedigree-based genome mapping project is to investigate and develop systems for implementing marker assisted selection to improve the efficiency of selection and increase the rate of genetic gain in breeding programs. Pedigree-based whole genome marker application provides a vehicle for incorporating marker technologies into applied breeding programs by bridging the gap between marker-trait association and marker implementation. We report on the development of protocols for implementation of pedigree-based whole genome marker analysis in breeding programs within the Australian northern winter cereals region. Examples of applications from the Queensland DPI&F wheat and barley breeding programs are provided, commenting on the use of microsatellites and other types of molecular markers for routine genomic analysis, the integration of genotypic, phenotypic and pedigree information for targeted wheat and barley lines, the genomic impacts of strong selection pressure in case study pedigrees, and directions for future pedigree-based marker development and analysis.

Using immunology and resistant sheep to beat the fly

New methods for controlling blowfly strike will be needed when mulesing is phased out and the availability or efficacy of insecticides for control of fly strike decreases. The Australian Sheep Industry CRC has pursued two approaches for the development of new methods to help control blowfly strike. In the first, genetic resistance of sheep to survival and growth of blowfly larvae was examined. Resistance to growth of larvae was heritable (0.29 ± 0.22). The trait was not associated with resistance to internal parasites, nor was it influenced by wool characteristics such as fibre diameter or coefficient of variation of fibre diameter. This new trait differs from resistance to fly strike associated with resistance to fleece rot. Because measurement of the trait is labour intensive, gene markers or correlated measures are needed before it will be suitable for industry adoption. The second approach examined the impact of larval products on the immmune system of the sheep. Larvae suppress the sheep immune system and thereby limit the ability of the sheep to reject the larvae. The immunosuppresive agent is being purified and strategies to abolish its activity are being explored. Abolition of immunosuppression would create opportunities for the development of new vaccines againts blowfly strike.

The genetic diversity of ampeloviruses in Australian pineapples and their association with mealybug wilt disease

Pineapple mealybug wilt-associated virus 1 (PMWaV-1), 2 (PMWaV-2) and -3 (PMWaV-3) have been detected in Australian commercial pineapple crops, along with a previously undescribed ampelovirus, for which the name Pineapple mealybug wilt-associated virus 5 (PMWaV-5) is proposed. Partial sequences extending from open reading frame 1b through to the heat shock protein homologue were obtained for PMWaV-1, -3 and -5. Phylogenetic analyses of selected regions of these sequences indicated that PMWaV-5 is a distinct species and most closely related to PMWaV-1. The amino acid sequence variation observed in the RNA-dependent RNA polymerase region of PMWaV-1 isolates was 95.8–98.4% and of PMWaV-3 isolates was 92.2–99.5%. In surveys of mealybug wilt disease (MWD) affected crops, none of the four viruses was clearly associated with the disease at all survey sites. A statistically significant association (P < 0.001) between the presence of PMWaV-2 and symptoms was observed at one survey site (site 3), but the virus was at a low incidence at the remaining three survey sites. By contrast, although PMWaV-1 and -3 were equally distributed between symptomless and MWD-affected plants at site 3, there was a statistically significant (P < 0.001) association between each of these two viruses and MWD at sites 1 and 4. At site 2, there was a statistically significant (P < 0.001) association only between PMWaV-3 and MWD. PMWaV-1 was the most commonly found of the four viruses and conversely PMWaV-5 was only occasionally found. Australian isolates of PMWaV-1, -2 and -3 were transmitted by the mealybug species Dysmicoccus brevipes.

The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus)

This study compares estimates of the census size of the spawning population with genetic estimates of effective current and long-term population size for an abundant and commercially important marine invertebrate, the brown tiger prawn (Penaeus esculentus). Our aim was to focus on the relationship between genetic effective and census size that may provide a source of information for viability analyses of naturally occurring populations. Samples were taken in 2001, 2002 and 2003 from a population on the east coast of Australia and temporal allelic variation was measured at eight polymorphic microsatellite loci. Moments-based and maximum-likelihood estimates of current genetic effective population size ranged from 797 to 1304. The mean long-term genetic effective population size was 9968. Although small for a large population, the effective population size estimates were above the threshold where genetic diversity is lost at neutral alleles through drift or inbreeding. Simulation studies correctly predicted that under these experimental conditions the genetic estimates would have non-infinite upper confidence limits and revealed they might be overestimates of the true size. We also show that estimates of mortality and variance in family size may be derived from data on average fecundity, current genetic effective and census spawning population size, assuming effective population size is equivalent to the number of breeders. This work confirms that it is feasible to obtain accurate estimates of current genetic effective population size for abundant Type III species using existing genetic marker technology.

The genetic structure of Australasian green turtles (Chelonia mydas): Exploring the geographical scale of genetic exchange

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

Achievements in forest tree improvement in Australia and New Zealand 6: Genetic improvement and conservation of Araucaria cunninghamii in Queensland

Araucaria cunninghamii (hoop pine) typically occurs as an emergent tree over subtropical and tropical rainforests, in a discontinuous distribution that extends from West Irian Jaya at about 0°30'S, through the highlands of Indonesian New Guinea and Papua New Guinea, along the east coast of Australia from 11°39'S in Queensland to 30°35'S in northern New South Wales. Plantations established in Queensland since the 1920s now total about 44000 ha, and constitute the primary source for the continuing supply of hoop pine quality timber and pulpwood, with a sustainable harvest exceeding 440 000 m3 y-1. Establishment of these managed plantations allowed logging of all native forests of Araucaria species (hoop pine and bunya pine, A. bidwillii) on state-owned lands to cease in the late 1980s, and the preservation of large areas of araucarian forest types within a system of state-owned and managed reserves. The successful plantation program with this species has been strongly supported by genetic improvement activities since the late 1940s - through knowledge of provenance variation and reproductive biology, the provision of reliable sources of improved seed, and the capture of substantial genetic gains in traits of economic importance (for example growth, stem straightness, internode length and spiral grain). As such, hoop pine is one of the few tropical tree species that, for more than half a century, has been the subject of continuous genetic improvement. The history of commercialisation and genetic improvement of hoop pine provides an excellent example of the dual economic and conservation benefits that may be obtained in tropical tree species through the integration of gene conservation and genetic improvement with commercial plantation development. This paper outlines the natural distribution and reproductive biology of hoop pine, describes the major achievements of the genetic improvement program in Queensland over the past 50+ y, summarises current understanding of the genetic variation and control of key selection traits, and outlines the means by which genetic diversity in the species is being conserved.

Domestication to crop improvement: Genetic resources for Sorghum and Saccharum (Andropogoneae).

Background: Both sorghum (Sorghum bicolor) and sugarcane (Saccharum officinarum) are members of the Andropogoneae tribe in the Poaceae and are each other's closest relatives amongst cultivated plants. Both are relatively recent domesticates and comparatively little of the genetic potential of these taxa and their wild relatives has been captured by breeding programmes to date. This review assesses the genetic gains made by plant breeders since domestication and the progress in the characterization of genetic resources and their utilization in crop improvement for these two related species. Genetic Resources: The genome of sorghum has recently been sequenced providing a great boost to our knowledge of the evolution of grass genomes and the wealth of diversity within S. bicolor taxa. Molecular analysis of the Sorghum genus has identified close relatives of S. bicolor with novel traits, endosperm structure and composition that may be used to expand the cultivated gene pool. Mutant populations (including TILLING populations) provide a useful addition to genetic resources for this species. Sugarcane is a complex polyploid with a large and variable number of copies of each gene. The wild relatives of sugarcane represent a reservoir of genetic diversity for use in sugarcane improvement. Techniques for quantitative molecular analysis of gene or allele copy number in this genetically complex crop have been developed. SNP discovery and mapping in sugarcane has been advanced by the development of high-throughput techniques for ecoTILLING in sugarcane. Genetic linkage maps of the sugarcane genome are being improved for use in breeding selection. The improvement of both sorghum and sugarcane will be accelerated by the incorporation of more diverse germplasm into the domesticated gene pools using molecular tools and the improved knowledge of these genomes.

QTL analysis of ergot resistance in sorghum

Sorghum ergot, caused predominantly by Claviceps africana Frederickson, Mantle, de Milliano, is a significant threat to the sorghum industry worldwide. The objectives of this study were firstly, to identify molecular markers linked to ergot resistance and to two pollen traits, pollen quantity (PQ) and pollen viability (PV), and secondly, to assess the relationship between the two pollen traits and ergot resistance in sorghum. A genetic linkage map of sorghum RIL population R931945-2-2 x IS 8525 (resistance source) was constructed using 303 markers including 36 SSR, 117 AFLP™, 148 DArT™ and two morphological trait loci. Composite interval mapping identified nine, five, and four QTL linked to molecular markers for percentage ergot infection (PCERGOT), PQ and PV, respectively, at a LOD >2.0. Co-location/linkage of QTL were identified on four chromosomes while other QTL for the three traits mapped independently, indicating that both pollen and non pollen-based mechanisms of ergot resistance were operating in this sorghum population. Of the nine QTL identified for PCERGOT, five were identified using the overall data set while four were specific to the group data sets defined by temperature and humidity. QTL identified on SBI-02 and SBI-06 were further validated in additional populations. This is the first report of QTL associated with ergot resistance in sorghum. The markers reported herein could be used for marker-assisted selection for this important disease of sorghum.

Identification of QTL for sugar-related traits in a sweet x grain sorghum (Sorghum bicolor L. Moench) recombinant inbred population

QTL for stem sugar-related and other agronomic traits were identified in a converted sweet (R9188) × grain (R9403463-2-1) sorghum population. QTL analyses were conducted using phenotypic data for 11 traits measured in two field experiments and a genetic map comprising 228 SSR and AFLP markers grouped into 16 linkage groups, of which 11 could be assigned to the 10 sorghum chromosomes (SBI-01 to SBI-10). QTL were identified for all traits and were generally co-located to five locations (SBI-01, SBI-03, SBI-05, SBI-06 and SBI-10). QTL alleles from R9188 were detected for increased sucrose content and sugar content on SBI-01, SBI-05 and SBI-06. R9188 also contributed QTL alleles for increased Brix on SBI-05 and SBI-06, and increased sugar content on SBI-03. QTL alleles from R9403463-2-1 were found for increased sucrose content and sucrose yield on SBI-10, and increased glucose content on SBI-07. QTL alleles for increased height, later flowering and greater total dry matter yield were located on SBI-01 of R9403463-2-1, and SBI-06 of R9188. QTL alleles for increased grain yield from both R9403463-2-1 and R9188 were found on SBI-03. As an increase in stem sugars is an important objective in sweet sorghum breeding, the QTL identified in this study could be further investigated for use in marker-assisted selection of sweet sorghum.

An assessment of the genetic relationship between sweet and grain sorghums, within Sorghum bicolor ssp. bicolor (L.) Moench, using AFLP markers

Compared to grain sorghums, sweet sorghums typically have lower grain yield and thick, tall stalks which accumulate high levels of sugar (sucrose, fructose and glucose). Unlike commercial grain sorghum (S. bicolor ssp. bicolor) cultivars, which are usually F1 hybrids, commercial sweet sorghums were selected as wild accessions or have undergone limited plant breeding. Although all sweet sorghums are classified within S. bicolor ssp. bicolor, their genetic relationship with grain sorghums is yet to be investigated. Ninety-five genotypes, including 31 sweet sorghums and 64 grain sorghums, representing all five races within the subspecies bicolor, were screened with 277 polymorphic amplified fragment length polymorphism (AFLP) markers. Cluster analysis separated older sweet sorghum accessions (collected in mid 1800s) from those developed and released during the early to mid 1900s. These groups were emphasised in a principle component analysis of the results such that sweet sorghum lines were largely distinguished from the others, particularly by a group of markers located on sorghum chromosomes SBI-08 and SBI-10. Other studies have shown that QTL and ESTs for sugar-related traits, as well as for height and anthesis, map to SBI-10. Although the clusters obtained did not group clearly on the basis of racial classification, the sweet sorghum lines often cluster with grain sorghums of similar racial origin thus suggesting that sweet sorghum is of polyphyletic origin within S. bicolor ssp. bicolor.