Population genetics of the frog-killing fungus Batrachochytrium dendrobatidis

Global amphibian decline by chytridiomycosis is a major environmental disaster that has been attributed to either recent fungal spread or environmental change that promotes disease. Here, we present a population genetic comparison of Batrachochytrium dendrobatidis isolates from an intensively studied region of frog decline, the Sierra Nevada of California. In support of a novel pathogen, we find low diversity, no amphibian-host specificity, little correlation between fungal genotype and geography, local frog extirpation by a single fungal genotype, and evidence of human-assisted fungus migration. In support of endemism, at a local scale, we find some diverse, recombining populations. Therefore neither epidemic spread nor endemism alone explains this particular amphibian decline. Recombination raises the possibility of resistant sporangia and a mechanism for rapid spread as well as persistence that could greatly complicate global control of the pathogen.

Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes

Despite international protection of white sharks Carcharodon carcharias, important conservation parameters such as abundance, population structure and genetic diversity are largely unknown. The tissue of 97 predominately juvenile white sharks sampled from spatially distant eastern and southwestern Australian coastlines was sequenced for the mitochondrial DNA (mtDNA) control region and genotyped with 6 nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (F-ST = 0.142, p < 0.0001), implying female reproductive philopatry. This concurs with recent satellite and acoustic tracking findings which suggest the sustained presence of discrete east coast nursery areas. Furthermore, population subdivision was found between the same regions with biparentally inherited micro satellite markers (F-ST = 0.009, p < 0.05), suggesting that males may also exhibit some degree of reproductive philopatry; 5 sharks captured along the east coast had mtDNA haplotypes that resembled western Indian Ocean sharks more closely than Australian/New Zealand sharks, suggesting that transoceanic dispersal, or migration resulting in breeding, may occur sporadically. Our most robust estimate of contemporary genetic effective population size was low and close to thresholds at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least 1, possibly 2, orders of magnitude below our historical effective size estimates. Population decline could expose these genetically isolated populations to detrimental genetic effects. Regional Australian white shark conservation management units should be implemented until genetic population structure, size and diversity can be investigated in more detail.

Diversity in the sunflower: Puccinia helianthi pathosystem in Australia

Sunflower rust caused by Puccinia helianthi is the most important disease of sunflower in Australia with the potential to cause significant yield losses in susceptible hybrids. Rapid and frequent virulence changes in the rust fungus population limit the effective lifespan of commercial cultivars and impose constant pressure on breeding programs to identify and deploy new sources of resistance. This paper contains a synopsis of virulence data accumulated over 25 years, and more recent studies of genotypic diversity and sexual recombination. We have used this synopsis, generated from both published and unpublished data, to propose the origin, evolution and distribution of new pathotypes of P. helianthi. Virulence surveys revealed that diverse pathotypes of P. helianthi evolve in wild sunflower populations, most likely because sexual recombination and subsequent selection of recombinant pathotypes occurs there. Wild sunflower populations provide a continuum of genetically heterogeneous hosts on which P. helianthi can potentially complete its sexual cycle under suitable environmental conditions. Population genetics analysis of a worldwide collection of P. helianthi indicated that Australian isolates of the pathogen are more diverse than non-Australian isolates. Additionally, the presence of the same pathotype in different genotypic backgrounds supported evidence from virulence data that sexual recombination has occurred in the Australian population of P. helianthi at some time. A primary aim of the work described was to apply our knowledge of pathotype evolution to improve resistance in sunflower to sunflower rust. Molecular markers were identified for a number of previously uncharacterised sunflower rust R-genes. These markers have been used to detect resistance genes in breeding lines and wild sunflower germplasm. A number of virulence loci that do not recombine were identified in P. helianthi. The resistance gene combinations corresponding to these virulence loci are currently being introgressed with breeding lines to generate hybrids with durable resistance to sunflower rust.

Rapid Global Expansion of the Fungal Disease Chytridiomycosis into Declining and Healthy Amphibian Populations.

The fungal disease chytridiomycosis, caused by Batrachochytrium dendrobatidis, is enigmatic because it occurs globally in both declining and apparently healthy (non-declining) amphibian populations. This distribution has fueled debate concerning whether, in sites where it has recently been found, the pathogen was introduced or is endemic. In this study, we addressed the molecular population genetics of a global collection of fungal strains from both declining and healthy amphibian populations using DNA sequence variation from 17 nuclear loci and a large fragment from the mitochondrial genome. We found a low rate of DNA polymorphism, with only two sequence alleles detected at each locus, but a high diversity of diploid genotypes. Half of the loci displayed an excess of heterozygous genotypes, consistent with a primarily clonal mode of reproduction. Despite the absence of obvious sex, genotypic diversity was high (44 unique genotypes out of 59 strains). We provide evidence that the observed genotypic variation can be generated by loss of heterozygosity through mitotic recombination. One strain isolated from a bullfrog possessed as much allelic diversity as the entire global sample, suggesting the current epidemic can be traced back to the outbreak of a single clonal lineage. These data are consistent with the current chytridiomycosis epidemic resulting from a novel pathogen undergoing a rapid and recent range expansion. The widespread occurrence of the same lineage in both healthy and declining populations suggests that the outcome of the disease is contingent on environmental factors and host resistance.

Simulating the evolution of glyphosate resistance in grains farming in northern Australia.

Background and Aims: The evolution of resistance to herbicides is a substantial problem in contemporary agriculture. Solutions to this problem generally consist of the use of practices to control the resistant population once it evolves, and/or to institute preventative measures before populations become resistant. Herbicide resistance evolves in populations over years or decades, so predicting the effectiveness of preventative strategies in particular relies on computational modelling approaches. While models of herbicide resistance already exist, none deals with the complex regional variability in the northern Australian sub-tropical grains farming region. For this reason, a new computer model was developed. Methods: The model consists of an age- and stage-structured population model of weeds, with an existing crop model used to simulate plant growth and competition, and extensions to the crop model added to simulate seed bank ecology and population genetics factors. Using awnless barnyard grass (Echinochloa colona) as a test case, the model was used to investigate the likely rate of evolution under conditions expected to produce high selection pressure. Key Results: Simulating continuous summer fallows with glyphosate used as the only means of weed control resulted in predicted resistant weed populations after approx. 15 years. Validation of the model against the paddock history for the first real-world glyphosate-resistant awnless barnyard grass population shows that the model predicted resistance evolution to within a few years of the real situation. Conclusions: This validation work shows that empirical validation of herbicide resistance models is problematic. However, the model simulates the complexities of sub-tropical grains farming in Australia well, and can be used to investigate, generate and improve glyphosate resistance prevention strategies.

White shark (Carcharodon carcharias) microsatellite genotypes (six loci) from regions along Australia's southern and eastern coastline.

Despite international protection of white sharks (Carcharodon carcharias), important conservation parameters such as abundance, population structure and genetic diversity are largely unknown. The tissue of 97 predominately juvenile white sharks sampled from spatially distant eastern and southwestern Australian coastlines was sequenced for the mitochondrial DNA (mtDNA) control region and genotyped with six nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (FST = 0.142, p < 0.001), implying female natal philopatry. This concords with recent satellite and acoustic tracking findings which suggest the sustained presence of discrete east coast nursery areas. Furthermore, population subdivision was found between the same regions with biparentally inherited microsatellite markers (FST = 0.009, p <0.05), suggesting that males may also exhibit some degree of reproductive philopatry. Five sharks captured along the east coast had mtDNA haplotypes that resembled western Indian Ocean sharks more closely than Australian/New Zealand sharks, suggesting that transoceanic dispersal or migration resulting in breeding may occur sporadically. Our most robust estimate of contemporary genetic effective population size was low and below the threshold at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least one, possibly two orders of magnitude below our historical effective size estimates. Further population decline could expose these genetically isolated populations to detrimental genetic effects. Regional Australian white shark conservation management units should be implemented until genetic population structure, size and diversity can be investigated in more detail. Reference: Blower, D. C., Pandolfi, J. M., Gomez-Cabrera, M. del C., Bruce, B. D. & Ovenden, J. R. (In press - April 2012). Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes. Marine Ecology Progress Series.

Biology, ecology and conservation of the Mobulidae

The Mobulidae are zooplanktivorous elasmobranchs comprising two recognized species of manta rays (Manta spp.) and nine recognized species of devil rays (Mobula spp.). They are found circumglobally in tropical, subtropical and temperate coastal waters. Although mobulids have been recorded for over 400 years, critical knowledge gaps still compromise the ability to assess the status of these species. On the basis of a review of 263 publications, a comparative synthesis of the biology and ecology of mobulids was conducted to examine their evolution, taxonomy, distribution, population trends, movements and aggregation, reproduction, growth and longevity, feeding, natural mortality and direct and indirect anthropogenic threats. There has been a marked increase in the number of published studies on mobulids since c. 1990, particularly for the genus Manta, although the genus Mobula remains poorly understood. Mobulid species have many common biological characteristics although their ecologies appear to be species-specific, and sometimes region-specific. Movement studies suggest that mobulids are highly mobile and have the potential to rapidly travel large distances. Fishing pressure is the major threat to many mobulid populations, with current levels of exploitation in target fisheries unlikely to be sustainable. Advances in the fields of population genetics, acoustic and satellite tracking, and stable-isotope and fatty-acid analyses will provide new insights into the biology and ecology of these species. Future research should focus on the uncertain taxonomy of mobulid species, the degree of overlap between their large-scale movement and human activities such as fisheries and pollution, and the need for management of inter-jurisdictional fisheries in developing nations to ensure their long-term sustainability. Closer collaboration among researchers worldwide is necessary to ensure standardized sampling and modelling methodologies to underpin global population estimates and status.

Random Mating Between Two Widely Divergent Mitochondrial Lineages of Cryptolestes ferrugineus (Coleoptera: Laemophloeidae): A Test of Species Limits in a Phosphine-Resistant Stored Product Pest

Effective pest management relies on accurate delimitation of species and, beyond this, on accurate species identification. Mitochondrial COI sequences are useful for providing initial indications in delimiting species but, despite acknowledged limitations in the method, many studies involving COI sequences and species problems remain unresolved. Here we illustrate how such impasses can be resolved with microsatellite and nuclear sequence data, to assess more directly the amount of gene flow between divergent lineages. We use a population genetics approach to test for random mating between two 8 ± 2% divergent COI lineages of the rusty grain beetle, Cryptolestes ferrugineus (Stephens). This species has become strongly resistant to phosphine, a fumigant used worldwide for disinfesting grain. The possibility of cryptic species would have significant consequences for resistance management, especially if resistance was confined to one mitochondrial lineage. We find no evidence of restricted gene flow or nonrandom mating across the two COI lineages of these beetles, rather we hypothesize that historic population structure associated with early Pleistocene climate changes likely contributed to divergent lineages within this species.

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.

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.

Genetic population structure of red snappers (Lutjanus malabaricus Bloch & Schneider, 1801 and Lutjanus erythropterus Bloch, 1790) in central and eastern Indonesia and northern Australia

The genetic population structure of red snapper Lutjanus malabaricus and Lutjanus erythropterus in eastern Indonesia and northern Australia was investigated by allozyme electrophoresis and sequence variation in the control region of mtDNA. Samples were collected from eight sites in Indonesia and four sites in northern Australia for both species. A total of 13 allozyme loci were scored. More variable loci were observed in L. malabaricus than in L. erythropterus. Sequence variation in the control region (left domain) of the mitochondrial genome was assessed by RFLP and direct sequencing. MtDNA haplotype diversity was high (L. erythropterus, 0.95 and L. malabaricus, 0.97), as was intraspecific sequence divergence, (L. erythropterus, 0.0-12.5% and L. malabaricus, 0.0-9.5%). The pattern of mtDNA haplotype frequencies grouped both species into two broad fisheries stocks with a genetic boundary either between Kupang and Sape (L. malabaricus) or between Kupang and Australian Timor Sea (L. erythropertus). The allozyme analyses revealed similar boundaries for L. erythropterus. Seven allozymes stocks compared to two mtDNA stocks of L. malabaricus including Ambon, which was not sampled with mtDNA, however, were reported. Possible reasons for differences in discrimination between the methods include: i) increased power of multiple allozyme loci over the single mtDNA locus, ii) insufficient gene sampling in the mtDNA control region and iii) relative evolutionary dynamics of nuclear (allozyme loci) and mitochondrial DNA in these taxa. Allozyme and haplotype data did not distinguish separate stocks among the four Australian locations nor the central Indonesian (Bali and Sape locations) for both L. malabaricus and L. erythropterus.

Population genetic structure of the brown tiger prawn, Penaeus esculentus, in tropical northern Australia

Eight polymorphic microsatellite loci were analysed in six population samples from four locations of the Australian endemic brown tiger prawn, Penaeus esculentus. Tests of Hardy-Weinberg equilibrium were generally in accord with expectations, with only one locus, in two samples, showing significant deviations. Three samples were taken in different years from the Exmouth Gulf. These showed no significant heterogeneity, and it was concluded that they were from a single panmictic population. A sample from Shark Bay, also on the west coast of Australia, showed barely detectable differentiation from Exmouth Gulf (F (ST) = 0 to 0.0014). A northeast sample from the Gulf of Carpentaria showed low (F (ST) = 0.008) but significant differentiation from Moreton Bay, on the east coast. However, Exmouth Gulf/Shark Bay samples were well differentiated from the Gulf of Carpentaria/Moreton Bay (F (ST) = 0.047-0.063). The data do not fit a simple isolation by distance model. It is postulated that the east-west differentiation largely reflects the isolation of east and west coast populations that occurred at the last glacial maximum when there was a land bridge between north-eastern Australia and New Guinea.

Mapping quantitative trait loci for resistance to Pratylenchus thornei from synthetic hexaploid wheat in the International Triticeae Mapping Initiative (ITMI) population

Root-lesion nematode (Pratylenchus thornei) is a serious pathogen of wheat in many countries. The International Triticeae Mapping Initiative (ITMI) population of recombinant inbred lines (RILs) was assessed for resistance to P. thornei to determine the chromosome locations of the resistance genes. The ITMI population is derived from a cross between the resistant synthetic hexaploid wheat W-7984 and a susceptible bread wheat cultivar Opata 85. Two years of phenotypic data for resistance to P. thornei were obtained in replicated glasshouse trials. Quantitative trait locus (QTL) analysis was performed using available segregation and map data for 114 RILs. A QTL on chromosome 6DS showed consistent effects for reduced nematode numbers (partial resistance) across years and accounted for 11% and 23% of the phenotypic variation. A second QTL for P. thornei resistance on chromosome 2BS accounted for an additional 19% and 5%. Restriction fragment length polymorphism (RFLP) and simple sequence repeat (SSR) markers associated with the QTLs are physically located in regions rich in major genes at the distal ends of the short chromosome arms of 6D and 2B. SSR markers with potential for marker-assisted selection of P. thornei resistance effective in different genetic backgrounds have been identified.

Population structure of South African and Australian Pyrenophora teres isolates.

There are two recognized forms of the disease net blotch of barley: the net form caused by Pyrenophora teres f. teres (PTT) and the spot form caused by P. teres f. maculata (PTM). In this study, amplified fragment length polymorphism analysis was used to investigate the genetic diversity and population structure of 60 PTT and 64 PTM isolates collected across Australia (66 isolates) and in the south-western Cape of South Africa (58 isolates). For comparison, P. tritici-repentis, Exserohilum rostratum and Bipolaris sorokiniana samples were also included in the analyses. Both distance-and model-based cluster analyses separated the PTT and PTM isolates into two strongly divergent genetic groups. Significant variation was observed both among the South African and Australian populations of PTT and PTM and among sampling locations for the PTT samples. Results suggest that sexual reproduction between the two forms is unlikely and that reproduction within the PTT and PTM groups occurs mainly asexually.

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.