Extremely Low Microsatellite Diversity but Distinct Population Structure in a Long-Lived Threatened Species, the Australian Lungfish Neoceratodus forsteri (Dipnoi)

The Australian lungfish is a unique living representative of an ancient dipnoan lineage, listed as ‘vulnerable’ to extinction under Australia’s Environment Protection and Biodiversity Conservation Act 1999. Historical accounts indicate this species occurred naturally in two adjacent river systems in Australia, the Burnett and Mary. Current day populations in other rivers are thought to have arisen by translocation from these source populations. Early genetic work detected very little variation and so had limited power to answer questions relevant for management including how genetic variation is partitioned within and among sub-populations. In this study, we use newly developed microsatellite markers to examine samples from the Burnett and Mary Rivers, as well as from two populations thought to be of translocated origin, Brisbane and North Pine. We test whether there is significant genetic structure among and within river drainages; assign putatively translocated populations to potential source populations; and estimate effective population sizes. Eleven polymorphic microsatellite loci genotyped in 218 individuals gave an average within-population heterozygosity of 0.39 which is low relative to other threatened taxa and for freshwater fishes in general. Based on F_ST values (average over loci = 0.11) and STRUCTURE analyses, we identify three distinct populations in the natural range, one in the Burnett and two distinct populations in the Mary. These analyses also support the hypothesis that the Mary River is the likely source of translocated populations in the Brisbane and North Pine rivers, which agrees with historical published records of a translocation event giving rise to these populations. We were unable to obtain bounded estimates of effective population size, as we have too few genotype combinations, although point estimates were low, ranging from 29 - 129. We recommend that, in order to preserve any local adaptation in the three distinct populations that they be managed separately.

Characterization, development and multiplexing of microsatellite markers in three commercially exploited reef fish and their application for stock identification

Thirty-four microsatellite loci were isolated from three reef fish species; golden snapper Lutjanus johnii, blackspotted croaker Protonibea diacanthus and grass emperor Lethrinus laticaudis using a next generation sequencing approach. Both IonTorrent single reads and Illumina MiSeq paired-end reads were used, with the latter demonstrating a higher quality of reads than the IonTorrent. From the 1–1.5 million raw reads per species, we successfully obtained 10–13 polymorphic loci for each species, which satisfied stringent design criteria. We developed multiplex panels for the amplification of the golden snapper and the blackspotted croaker loci, as well as post-amplification pooling panels for the grass emperor loci. The microsatellites characterized in this work were tested across three locations of northern Australia. The microsatellites we developed can detect population differentiation across northern Australia and may be used for genetic structure studies and stock identification.

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.

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.

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.

Confirmation of QTL mapping and marker validation for partial seedling resistance to crown rot in wheat line '2-49'

We have tested the efficacy of putative microsatellite single sequence repeat (SSR) markers, previously identified in a 2-49 (Gluyas Early/Gala) × Janz doubled haploid wheat (Triticum aestivum) population, as being linked to partial seedling resistance to crown rot disease caused by Fusarium pseudograminearum. The quantitative trait loci (QTLs) delineated by these markers have been tested for linkage to resistance in an independent Gluyas Early × Janz doubled haploid population. The presence of a major QTL on chromosome 1DL (QCr.usq-1D1) and a minor QTL on chromosome 2BS (QCr.usq-2B1) was confirmed. However, a putative minor QTL on chromosome 2A was not confirmed. The QTL on 1D was inherited from Gluyas Early, a direct parent of 2-49, whereas the 2B QTL was inherited from Janz. Three other putative QTLs identified in 2-49 × Janz (on 1AL, 4BL, and 7BS) were inherited by 2-49 from Gala and were not able to be confirmed in this study. The screening of SSR markers on a small sample of elite wheat genotypes indicated that not all of the most tightly linked SSR markers flanking the major QTLs on 1D and 1A were polymorphic in all backgrounds, indicating the need for additional flanking markers when backcrossing into some elite pedigrees. Comparison of SSR haplotypes with those of other genotypes exhibiting partial crown rot resistance suggests that additional, novel sources of crown rot resistance are available.

Genetic variation within two sympatric spotted gum eucalypts exceeds between taxa variation

Population substructure and hybridization, among other factors, have the potential to cause erroneous associations in linkage disequilibrium (LD) mapping. Two closely related spotted gum eucalypts, Corymbia variegata and C. henryi (Myrtaceae) occur in sympatry in the east coast of Australia and potentially interbreed. They are morphologically similar but are distinguished as separate species based on capsule and foliage size. To determine whether they hybridize in nature and its implications for LD mapping, we investigated the level of molecular divergence between the two species at two sympatric locations separated by 300 kilometres. Very few individuals of intermediate morphology were identified, despite the two species occurring only metres apart. Analysis of genetic structure using 12 microsatellite loci showed that genetic differentiation between populations of the same species at different locations (FST = 0.07 for both species; p = 0.0001) was significantly higher than that observed between species at each location (mean FST = 0.02 and 0.04 for Cherry tree and Bunyaville respectively; p = 0.0001; all Mann-Whitney U-test p ≤ 0.01). No species-specific alleles or significant allele frequency differences were detected within a site, suggesting recurrent local gene flow between the two species. The lack of significant allele frequency differences implies no population stratification along taxonomic lines. This suggested that there is little concern for cryptic hybridization when sampling from sites of sympatry for LD mapping.

Heritability estimates for growth in the tropical abalone Haliotis asinina using microsatellites to assign parentage

The tropical abalone Haliotis asinina is a wild-caught and cultured species throughout the Indo-Pacific as well as being an emerging model species for the study of haliotids. H. asinina has the fastest recorded natural growth rate of any abalone and reaches sexual maturity within one year. As such, it is a suitable abalone species for selective breeding for commercially important traits such as rapid growth. Estimating the amount of variation in size that is attributable to heritable genetic differences can assist the development of such a selective breeding program. Here we estimated heritability for growth-related traits at 12 months of age by creating a single cohort of 84 families in a full-factorial mating design consisting of 14 sires and 6 dams. Of 500 progeny sampled, 465 were successfully assigned to their parents based on shared alleles at 5 polymorphic microsatellite loci. Using an animal model, heritability estimates were 0.48 ± 0.15 for shell length, 0.38 ± 0.13 for shell width and 0.36 ± 0.13 for weight. Genetic correlations were > 0.98 between shell parameters and weight, indicating that breeding for weight gains could be successfully achieved by selecting for shell length.

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

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

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.

Identification of quantitative trait loci for resistance to two species of root-lesion nematode (Pratylenchus thornei and P. neglectus) in wheat

Pratylenchus thornei and P. neglectus are two species of root-lesion nematode that cause substantial yield losses in wheat. No commercially available wheat variety has resistance to both species. A doubled-haploid population developed from a cross between the synthetic hexaploid wheat line CPI133872 and the bread wheat Janz was used to locate and tag quantitative trait loci (QTLs) associated with resistance to both P. thornei and P. neglectus. Wheat plants were inoculated with both species of nematode in independent replicated glasshouse trials repeated over 2 years. Known locations of wheat microsatellite markers were used to construct a framework map. After an initial single-marker analysis to detect marker-trait linkages, chromosome regions associated with putative QTLs were targetted with microsatellite markers to increase map density in the chromosome regions of interest. In total, 148 wheat microsatellite markers and 21 amplified fragment length polymorphism markers were mapped. The codominant microsatellite marker Xbarc183 on the distal end of chromosome 6DS was allelic for resistance to both P. thornei and P. neglectus. The QTL were designated QRlnt.lrc-6D.1 and QRlnn.lrc-6D.1, for the 2 traits, respectively. The allele inherited from CPI133872 explained 22.0-24.2% of the phenotypic variation for P. thornei resistance, and the allele inherited from Janz accounted for 11.3-14.0% of the phenotypic variation for P. neglectus resistance. Composite interval mapping identified markers that flank a second major QTL on chromosome 6DL (QRlnt.lrc-6D.2) that explained 8.3-13.4% of the phenotypic variation for P. thornei resistance. An additional major QTL associated with P. neglectus resistance was detected on chromosome 4DS (QRlnn.lrc-4D.1) and explained a further 10.3-15.4% of the phenotypic variation. The identification and tagging of nematode resistance genes with molecular markers will allow appropriate allele combinations to be selected, which will aid the successful breeding of wheat with dual nematode resistance.

Two sympatric spotted gum species are molecularly homogeneous.

Large fruited spotted gum eucalypt Corymbia henryi occurs sympatrically with small fruited spotted gum Corymbia citriodora subspecies variegata over a large portion of its range on the east coast of Australia. The two taxa are interfertile, have overlapping flowering times and share a common set of insect and vertebrate pollinators. Previous genetic analysis of both taxa from two geographically remote sites suggested that the two were morphotypes rather than genetically distinct species. In this study we further explore this hypothesis of genic species by expanding sampling broadly through their sympatric locations and examine local-scale spatial genetic structure in stands that differ in species and age composition. Delineation of populations at five microsatellite loci, using an individual-based approach and Bayesian modelling, as well as clustering of individuals based on allele frequencies showed the two species to be molecularly homogeneous. Genetic structure aligned largely with geographic areas of origin, and followed an isolation-by-distance model, where proximal populations were generally less differentiated than more distant ones. At the stand level, spotted gums also generally showed little structure consistent with the high levels of gene flow inferred across the species range. Disturbances in the uniformity of structuring were detected, however, and attributed to localised events giving rise to even aged stands, probably due to regeneration from a few individuals following fire.

Diagnosis of a naturally occurring dual infection of layer chickens with fowlpox virus and gallid herpesvirus 1 (infectious laryngotracheitis virus).

An outbreak of acute respiratory disease in layers was diagnosed as being of dual nature due to fowlpox and infectious laryngotracheitis using a multidisciplinary approach including virus isolation, histopathology, electron microscopy and polymerase chain reaction (PCR). The diagnosis was based on virus isolation of gallid herpesvirus 1 (GaHV-1) in chicken kidney cells and fowlpox virus (FWPV) in 9-day-old chicken embryonated eggs inoculated via the chorioallantoic membrane. The histopathology of tracheas from dead birds revealed intra-cytoplasmic and intra-nuclear inclusions suggestive of poxvirus and herpesvirus involvement. The presence of FWPV was further confirmed by electron microscopy, PCR and histology. All FWPV isolates contained the long terminal repeats of reticuloendotheliosis virus as demonstrated by PCR. GaHV-1 isolates were detected by PCR and were shown to have a different restriction fragment length polymorphism pattern when compared with the chicken embryo origin SA2 vaccine strain; however, they shared the same pattern with the Intervet chicken embryo origin vaccine strain. This is a first report of dual infection of chickens with GaHV-1 and naturally occurring FWPV with reticuloendotheliosis virus insertions. Further characterization of the viruses was carried out and the results are reported here.

Reproductive Biology of Corymbia citriodora subsp. variegata and Effective Pollination across its Native Range.

The spotted gum species complex represents a group of four eucalypt hardwood taxa that have a native range that spans the east coast of Australia, with a morphological cline from Victoria to northern Queensland. Of this group, Corymbia citriodora subsp. variegata (CCV) is widespread in south-eastern Queensland and northern New South Wales. It is currently the most commonly harvested native hardwood in Queensland. However, little basic knowledge of the reproductive biology of the species is available to inform genetic improvement and resource management programmes. Here we take an integrative approach, using both field and molecular data, to identify ecological factors important to mating patterns in native populations of CCV. Field observation of pollinator visitation and flowering phenology of 20 trees showed that foraging behaviour of pollinator guilds varies depending on flowering phenology and canopy structure. A positive effect of tree mean flowering effort was found on insect visitation, while bat visitation was predicted by tree height and by the number of trees simultaneously bearing flowers. Moreover, introduced honeybees were observed frequently, performing 73% of detected flower visits. Conversely, nectar-feeding birds and mammals were observed sporadically with lorikeets and honeyeaters each contributing to 11% of visits. Fruit bats, represented solely by the grey-headed flying fox, performed less than 2% of visits. Genotyping at six microsatellite markers in 301 seeds from 17 families sampled from four of Queensland's native forests showed that CCV displays a mixed-mating system that is mostly outcrossing (tm = 0.899 ± 0.021). Preferential effective pollination from near-neighbours was detected by means of maximum-likelihood paternity analysis with up to 16% of reproduction events resulting from selfing. Forty to 48% of fertilising pollen was also carried from longer distance (>60 m). Marked differences in foraging behaviour and visitation frequency between observed pollinator guilds suggests that the observed dichotomy of effective pollen movement in spotted gums may be due to frequent visit from introduced honeybees favouring geitonogamy and sporadic visits from honeyeaters and fruit bats resulting in potential long-distance pollinations.

Population structure, effective population size and adverse effects of stocking in the endangered Australian eastern freshwater cod Maccullochella ikei

Microsatellite markers were used to examine spatio-temporal genetic variation in the endangered eastern freshwater cod Maccullochella ikei in the Clarence River system, eastern Australia. High levels of population structure were detected. A model-based clustering analysis of multilocus genotypes identified four populations that were highly differentiated by F-statistics (FST = 0· 09 − 0· 49; P < 0· 05), suggesting fragmentation and restricted dispersal particularly among upstream sites. Hatchery breeding programmes were used to re-establish locally extirpated populations and to supplement remnant populations. Bayesian and frequency-based analyses of hatchery fingerling samples provided evidence for population admixture in the hatchery, with the majority of parental stock sourced from distinct upstream sites. Comparison between historical and contemporary wild-caught samples showed a significant loss of heterozygosity (21%) and allelic richness (24%) in the Mann and Nymboida Rivers since the commencement of stocking. Fragmentation may have been a causative factor; however, temporal shifts in allele frequencies suggest swamping with hatchery-produced M. ikei has contributed to the genetic decline in the largest wild population. This study demonstrates the importance of using information on genetic variation and population structure in the management of breeding and stocking programmes, particularly for threatened species.