Raw data for project "Development of a DNA based aging technique for use in fisheries assessments" Fisheries and Development Corporation project 2007/033

The primary aim of this study was to determine the relationship between telomere length and age in a range of marine invertebrates including abalone (Haliotis spp) oysters (Saccostrea glomerata), spiny lobsters (Sagmariasus verreauxi formerly Jasus verreauxi and Jasus edwardsii) and school prawns (Metapenaeus macleayi). Additionally, this relationship was studied in a vertebrate organism using the freshwater fish Silver perch (Bidyanus bidyanus). Telomere length differences between tissues were also examined in some species such as Saccostrea glomerata, Sagmariasus verreauxi and Bidyanus bidyanus. In some cases cultured specimens of known age were used and this is quoted in the spreadsheets. For other wild-caught specimens where age was not known, size was used as a proxy for age. This may be a broad size class, or be determined by shell size or carapace length depending on the organism. Each spreadsheet contains raw data of telomere length estimates from Terminal Restriction Fragment Assays (TRF) for various individuals of each species including appropriate details such as age or size and tissue. Telomere length estimates are given in base pairs (bp). In most cases replicate experiments were conducted on groups of samples three times but on a small number of occasions only two replicate experiments were conducted. Further description of the samples can be found in final report of FRDC 2007/033. The arithmetic average for each individual (sample ID) across the two or three replicate experiments is also given. Bidyanus bidyanus (SilverPerch) Two sheets are contained within. a) Comparison of telomere length between different tissues (heart, liver and muscle) within the three year old age class - two replicate experiments were conducted. b) Comparison of telomere length between fish of different but known ages (0.25, 1, 2, and 3 years old) in each of three tissues, heart, liver and muscle – three replicate experiments were conducted per tissue. Haliotis spp (Abalone species) Three species were tested. H. asinina Telomere length was compared in two age classes-11 month and 18 month old abalone using muscle tissue from the foot. Within gel-variation was also estimated using a single sample run three times on one gel (replicate experiment). H. laevigata x H. rubra hybrids Telomere length was compared in three known age classes – two, three and four years old using muscle tissue from the foot. H. rubra Telomere length was compared in a range of different sized abalone using muscle tissue from the foot. Shell size is also given for each abalone Saccostrea glomerata Three sheets are contained within the file. a) Samples came from Moreton Bay Queensland in 2007. Telomere length was compared in two tissues (gill and mantle) of oysters in three age groups (1, 3 and 4 years) b) Samples came from Moreton Bay Queensland in 2009. Telomere length was compared in three age classes using DNA from gill tissue only c) Samples came from Wallis Lake, New South Wales. Telomere length was estimated from whole body minus the shell from 1 year old oysters, gill tissue of 3 age classes (1.5 years, 3 and 4 years), mantle tissue of two age classes (3 and 4 years). Sagmariasus verreauxi (formerly Jasus verreauxi) Telomere length was estimated from abdomen tissue of puerulus, gill and muscle tissue of 3 year old, large and very large size classes of lobsters. Jasus edwardsii Telomere length was measured in two size classes of lobsters- adults of varying sizes using muscle tissue and puerulus using tissues from the abdomen minus the exoskeleton. Metapenaeus macleayi Telomere length was measured in three size classes of school prawns adults. Muscle tissue was used, minus the exoskeleton.

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.

Assessment and rationalization of genetic diversity of Papua New Guinea taro (Colocasia esculenta) using SSR DNA fingerprinting

Taro (Colocasia esculenta) accessions were collected from 15 provinces of Papua New Guinea (PNG). The collection, totalling 859 accessions was collated for characterization and a core collection of 81 accessions (10%) was established on the basis of characterization data generated on 30 agro-morphological descriptors, and DNA fingerprinting using seven SSR primers. The selection of accessions was based on cluster analysis of the morphological data enabling initial selection of 20% accessions. The 20% sample was then reduced and rationalized to 10% based on molecular data generated by SSR primers. This represents the first national core collection of any species established in PNG based on molecular markers. The core has been integrated with core from other Pacific Island countries, contributing to a Pacific regional core collection, which is conserved in vitro in the South Pacific Regional Germplasm Centre at Fiji. The core collection is a valuable resource for food security of the South Pacific region and is currently being utilized by the breeding programmes of small Pacific Island countries to broaden the genetic base of the crop.

Concordant molecular and phenotypic data delineate new taxonomy and conservation priorities for the endangered mountain yellow-legged frog

The mountain yellow-legged frog Rana muscosa sensu lato, once abundant in the Sierra Nevada of California and Nevada, and the disjunct Transverse Ranges of southern California, has declined precipitously throughout its range, even though most of its habitat is protected. The species is now extinct in Nevada and reduced to tiny remnants in southern California, where as a distinct population segment, it is classified as Endangered. Introduced predators (trout), air pollution and an infectious disease (chytridiomycosis) threaten remaining populations. A Bayesian analysis of 1901 base pairs of mitochondrial DNA confirms the presence of two deeply divergent clades that come into near contact in the Sierra Nevada. Morphological studies of museum specimens and analysis of acoustic data show that the two major mtDNA clades are readily differentiated phenotypically. Accordingly, we recognize two species, Rana sierrae, in the northern and central Sierra Nevada, and R. muscosa, in the southern Sierra Nevada and southern California. Existing data indicate no range overlap. These results have important implications for the conservation of these two species as they illuminate a profound mismatch between the current delineation of the distinct population segments (southern California vs. Sierra Nevada) and actual species boundaries. For example, our study finds that remnant populations of R. muscosa exist in both the southern Sierra Nevada and the mountains of southern California, which may broaden options for management. In addition, despite the fact that only the southern California populations are listed as Endangered, surveys conducted since 1995 at 225 historic (1899-1994) localities from museum collections show that 93.3% (n=146) of R. sierrae populations and 95.2% (n=79) of R. muscosa populations are extinct. Evidence presented here underscores the need for revision of protected population status to include both species throughout their ranges.

An identification tool for the Australian weedy Sporobolus species based on random amplified polymorphic DNA (RAPD) profiles

Based on morphological features alone, there is considerable difficulty in identifying the 5 most economically damaging weed species of Sporobolus [viz. S. pyramidalis P. Beauv., S. natalensis (Steud.) Dur and Schinz, S. fertilis (Steud.) Clayton, S. africanus (Poir.) Robyns and Tourney, and S. jacquemontii Kunth.] found in Australia. A polymerase chain reaction (PCR)-based random amplified polymorphic DNA (RAPD) technique was used to create a series of genetic markers that could positively identify the 5 major weeds from the other less damaging weedy and native Sporobolus species. In the initial RAPD profiling experiment, using arbitrarily selected primers and involving 12 species of Sporobolus, 12 genetic markers were found that, when used in combination, could consistently identify the 5 weedy species from all others. Of these 12 markers, the most diagnostic were UBC51490 for S. pyramidalis and S. natalensis; UBC43310.2000.2100 for S. fertilis and S. africanus; and ORA20850 and UBC43470 for S. jacquemontii. Species-specific markers could be found only for S. jacquemontii. In an effort to understand why there was difficulty in obtaining species-specific markers for some of the weedy species, a RAPD data matrix was created using 40 RAPD products. These 40 products amplified by 6 random primers from 45 individuals belonging to 12 species, were then subjected to numerical taxonomy and multivariate system (NTSYS pc version 1.70) analysis. The RAPD similarity matrix generated from the analysis indicated that S. pyramidalis was genetically more similar to S. natalensis than to other species of the 'S. indicus complex'. Similarly, S. jacquemontii was more similar to S. pyramidalis, and S. fertilis was more similar to S. africanus than to other species of the complex. Sporobolus pyramidalis, S. jacquemontii, S. africanus, and S. creber exhibited a low within-species genetic diversity, whereas high genetic diversity was observed within S. natalensis, S. fertilis, S. sessilis, S. elongates, and S. laxus. Cluster analysis placed all of the introduced species (major and minor weedy species) into one major cluster, with S. pyramidalis and S. natalensis in one distinct subcluster and S. fertilis and S. africanus in another. The native species formed separate clusters in the phenograms. The close genetic similarity of S. pyramidalis to S. natalensis, and S. fertilis to S. africanus may explain the difficulty in obtaining RAPD species-specific markers. The importance of these results will be within the Australian dairy and beef industries and will aid in the development of integrated management strategy for these weeds.

Mitochondrial DNA supports the identification of two endangered river sharks (Glyphis glyphis and Glyphis garricki) across northern Australia

The river sharks (genus Glyphis) are a small group of poorly known sharks occurring in tropical rivers and estuarine waters across northern Australia, south-east Asia and the subcontinent. The taxonomy of the genus has long been unclear due to very few individuals having been caught and examined, resulting in a paucity of data regarding their distribution, biology and ecology. Only recently has attention focussed on the two Australian species, G. glyphis and G. garricki. This study is a result of a rare opportunity to collate the few samples that have been collected from these species and the bull shark Carcharhinus leucas, which shares an overlapping range. These samples were analysed using the DNA barcoding approach (cox1 mitochondrial gene), compared with six other species of carcharhinids and evaluated in light of the current taxonomic classification. Nine species-specific nucleotide differences were found between G. glyphis and G. garricki and no intra-specific variation provides strong support for the separation into distinct species. Significant differences were also observed at the inter-generic level, with Glyphis forming a distinct clade from Carcharhinus. This study provides the basis for future molecular studies required to better address conservation issues confronting G. glyphis and G. garricki in Australia.

DNA bar-coding reveals source and patterns of Thaumastocoris peregrinus invasions in South Africa and South America.

Thaumastocoris peregrinus is a recently introduced invertebrate pest of non-native Eucalyptus plantations in the Southern Hemisphere. It was first reported from South Africa in 2003 and in Argentina in 2005. Since then, populations have grown explosively and it has attained an almost ubiquitous distribution over several regions in South Africa on 26 Eucalyptus species. Here we address three key questions regarding this invasion, namely whether only one species has been introduced, whether there were single or multiple introductions into South Africa and South America and what the source of the introduction might have been. To answer these questions, bar-coding using mitochondrial DNA (COI) sequence diversity was used to characterise the populations of this insect from Australia, Argentina, Brazil, South Africa and Uruguay. Analyses revealed three cryptic species in Australia, of which only T. peregrinus is represented in South Africa and South America. Thaumastocoris peregrinus populations contained eight haplotypes, with a pairwise nucleotide distance of 0.2-0.9% from seventeen locations in Australia. Three of these haplotypes are shared with populations in South America and South Africa, but the latter regions do not share haplotypes. These data, together with the current distribution of the haplotypes and the known direction of original spread in these regions, suggest that at least three distinct introductions of the insect occurred in South Africa and South America before 2005. The two most common haplotypes in Sydney, one of which was also found in Brisbane, are shared with the non-native regions. Sydney populations of T. peregrinus, which have regularly reached outbreak levels in recent years, might thus have served as source of these three distinct introductions into other regions of the Southern Hemisphere.

DNA-based identification of Quambalaria pitereka causing severe leaf blight of Corymbia citriodora in China

Quambalaria spp. include serious plant pathogens, causing leaf and shoot blight of Corymbia and Eucalyptus spp. In this study, a disease resembling Quambalaria leaf blight was observed on young Corymbia citriodora trees in a plantation in the Guangdong Province of China. Comparisons of rDNA sequence data showed that the causal agent of the disease is Q. pitereka. This study provides the first report of Quambalaria leaf blight from China, and it is also the first time that this pathogen has been found on trees outside the native range of Eucalypts.

Development of a DNA based aging technique for use in fisheries assessments. Final report, Australian Fisheries Research and Development Corporation Project 2007/033.

The productivity of a fisheries resource can be quantified from estimates of recruitment, individual growth and natural and fisheries-related mortality, assuming the spatial extent of the resource has been quantified and there is minimal immigration or emigration. The sustainability of a fisheries resource is facilitated by management controls such as minimum and maximum size limits and total allowable catch. Minimum size limits are often set to allow individuals the opportunity to reproduce at least once before the chance of capture. Total allowable catches are a proportion of the population biomass, which is estimated based on known reproduction, recruitment, mortality and growth rates. In some fisheries, however, management actions are put in place without quantification of the resource through the stock assessment process. This occurs because species-specific information, for example individual growth, may not be available. In these circumstances, management actions need to be precautionary to protect against future resource collapse, but this often means that the resource is lightly exploited. Consequently, the productivity of the resource is not fully realised. Australia’s most valuable fisheries are invertebrate fisheries (Australian Department of Agriculture Fisheries and Forestry, 2008). For example, Australian fisheries (i.e. excluding aquaculture) production of crustaceans (largely prawns, rock lobster and crab) was 41,000 tonnes in 2006/7, worth $778 million. Production from mollusc (largely abalone, scallops, oysters and squid) fisheries was 39,000 tonnes, worth $502 million. Together, in 2006/7 crustacean and mollusc fisheries represented 58% of the total value of Australian wild fisheries production. Sustainable management of Australia’s invertebrate fisheries is frustrated by the lack of data on species-specific growth rates. This project investigated a new method to estimate age, and hence individual growth rates, in invertebrate fisheries species. The principle behind the new aging method was that telomeres (i.e. DNA end-caps of chromosomes) get shorter as an individual gets older. We studied commercial crustacean and molluscan species. A vertebrate fish species (silver perch, Bidyanus bidyanus) was used as a control to standardise our work against the literature. We found a clear relationship between telomere length and shell size for temperate abalone (Haliotis rubra). Further research is recommended before the method can be implemented to assist management of wildharvested abalone populations. Age needs to be substituted for shell size in the relationship and it needs to be studied for abalone from several regions. This project showed that telomere length declined with increasing age in Sydney rock oysters (Saccostrea glomerata) and was affected by regional variation. A relationship was not apparent between telomere length and age (or size as a surrogate for age) for crustacean species (school prawns, Metapenaeus macleayi; eastern rock lobster, Sagmariasus verreauxi; southern rock lobster, Jasus edwardsii; and spanner crabs, Ranina ranina). For school prawns, there was no difference between telomere length in males and females. Further research is recommended, however, as telomeric DNA from crustaceans was difficult to analyse using the terminal restriction fragment (TRF) assay. Telomere lengths of spanner crabs and lobsters were at the upper limit of resolution of the assay used and results were affected by degradation and possible contamination of telomeric DNA. It is possible that telomere length is an indicator of remaining lifespan in molluscan and crustacean individuals, as suggested for some vertebrate species (e.g. Monaghan, 2010). Among abalone of similar shell size and among lobster pueruli, there was evidence of individuals having significantly longer or shorter telomeres than the group average. At a population level, this may be a surrogate for estimates of future natural mortality, which may have usefulness in the management of those populations. The method used to assay telomere length (terminal restriction fragment assay) performed adequately for most species, but it was too expensive and time-consuming to be considered a useful tool for gathering information for fisheries management. Research on alternative methods is strongly recommended.

Comparative analysis of Tritrichomonas foetus (Riedmuller, 1928) cat genotype, T. foetus (Riedmuller, 1928) cattle genotype and Tritrichomonas suis (Davaine, 1875) at 10 DNA loci

The parasitic protists in the genus Tritrichomonas cause significant disease in domestic cattle and cats. To assess the genetic diversity of feline and bovine isolates of Tritrichomonas foetus (Riedmuller, 1928) Wenrich and Emmerson, 1933, we used 10 different genetic regions, namely the protein coding genes of cysteine proteases 1,2 and 4-9 (CP1, 2, 4-9) involved in the pathogenesis of the disease caused by the parasite. The cytosolic malate dehydrogenase 1 (MDH1) and internal transcribed spacer region 2 of the rDNA unit (ITS2) were included as additional markers. The gene sequences were compared with those of Tritrichomonas suis (Davaine. 1875) Morgan and Hawkins, 1948 and Tritrichomonas mobilensis Culberson et al., 1986. The study revealed 100% identity for all 10 genes among all feline isolates (=T. foetus cat genotype), 100% identity among all bovine isolates (=T. foetus cattle genotype) and a genetic distinctness of 1% between the cat and cattle genotypes of T. foetus. The cattle genotype of T. foetus was 100% identical to T. suis at nine loci (CP1, 2,4-8, ITS2, MDH1). At CP9, three out of four T. suis isolates were identical to the T. foetus cattle genotype, while the T. suis isolate SUI-H3B sequence contained a single unique nucleotide substitution. Tritrichomonas mobilensis was 0.4% and 0.7% distinct from the cat and cattle genotypes of T. foetus, respectively. The genetic differences resulted in amino acid changes in the CP genes, most pronouncedly in CP2, potentially providing a platform for elucidation of genotype-specific host-pathogen interactions of T. foetus. On the basis of this data we judge T. suis and T. foetus to be subjective synonyms. For the first time, on objective nomenclatural grounds, the authority of T. suis is given to Davaine, 1875, rather than the commonly cited Gruby and Delafond, 1843. To maintain prevailing usage of T. foetus, we are suppressing the senior synomym T. suis Davaine, 1875 according to Article 23.9, because it has never been used as a valid name after 1899 and T. foetus is widely discussed as the cause of bovine trichomonosis. Thus bovine, feline and porcine isolates should all be given the name T. foetus. This promotes the stability of T. foetus for the veterinary and economically significant venereal parasite causing bovine trichomonosis. (C) 2012 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.