MOE and MOR assessment technologies for improving graded recovery of exotic pines in Australia.

This project was designed to provide the structural softwood processing industry with the basis for improved green and dry grading to allow maximise MGP grade yields, consistent product performance and reduced processing costs. To achieve this, advanced statistical techniques were used in conjunction with state-of-the-art property measurement systems. Specifically, the project aimed to make two significant steps forward for the Australian structural softwood industry: • assessment of technologies, both existing and novel, that may lead to selection of a consistent, reliable and accurate device for the log yard and green mill. The purpose is to more accurately identify and reject material that will not make a minimum grade of MGP10 downstream; • improved correlation of grading MOE and MOR parameters in the dry mill using new analytical methods and a combination of devices. The three populations tested were stiffness-limited radiata pine, strength-limited radiata pine and Caribbean pine. Resonance tests were conducted on logs prior to sawmilling, and on boards. Raw data from existing in-line systems were captured for the green and dry boards. The dataset was analysed using classical and advanced statistical tools to provide correlations between data sets and to develop efficient strength and stiffness prediction equations. Stiffness and strength prediction algorithms were developed from raw and combined parameters. Parameters were analysed for comparison of prediction capabilities using in-line parameters, off-line parameters and a combination of in-line and off-line parameters. The results show that acoustic resonance techniques have potential for log assessment, to sort for low stiffness and/or low strength, depending on the resource. From the log measurements, a strong correlation was found between the average static MOE of the dried boards within a log and the predicted value. These results have application in segregating logs into structural and non-structural uses. Some commercial technologies are already available for this application such as Hitman LG640. For green boards it was found that in-line and laboratory acoustic devices can provide a good prediction of dry static MOE and moderate prediction for MOR.There is high potential for segregating boards at this stage of processing. Grading after the log breakdown can improve significantly the effectiveness of the mill. Subsequently, reductions in non-structural volumes can be achieved. Depending on the resource it can be expected that a 5 to 8 % reduction in non structural boards won’t be dried with an associated saving of $70 to 85/m3. For dry boards, vibration and a standard Metriguard CLT/HCLT provided a similar level of prediction on stiffness limited resource. However, Metriguard provides a better strength prediction in strength limited resources (due to this equipment’s ability to measure local characteristics). The combination of grading equipment specifically for stiffness related predictors (Metriguard or vibration) with defect detection systems (optical or X-ray scanner) provides a higher level of prediction, especially for MOR. Several commercial technologies are already available for acoustic grading on board such those from Microtec, Luxscan, Falcon engineering or Dynalyse AB for example. Differing combinations of equipment, and their strategic location within the processing chain, can dramatically improve the efficiency of the mill, the level of which will vary depending of the resource. For example, an initial acoustic sorting on green boards combined with an optical scanner associated with an acoustic system for grading dry board can result in a large reduction of the proportion of low value low non-structural produced. The application of classical MLR on several predictors proved to be effective, in particular for MOR predictions. However, the usage of a modern statistics approach(chemometrics tools) such as PLS proved to be more efficient for improving the level of prediction. Compared to existing technologies, the results of the project indicate a good improvement potential for grading in the green mill, ahead of kiln drying and subsequent cost-adding processes. The next stage is the development and refinement of systems for this purpose.

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.

Grey Mackerel (Scomberomorus semifasciatus) microsatellite genotypes (nine loci) representing twelve regions along Australia's northern coastline from the eastern coast of Queensland to the western coast of Western Australia.

Scomberomorus semifasciatus is an Australian endemic found in tropical, coastal waters from eastern to western Australia. Commercial and recreational exploitation is common and regulated by state-based authorities. This study used mitochondrial DNA sequence and microsatellite markers to elucidate the population structure of Scomberomorus semifasciatus collected from twelve, equidistant sampling locations. Samples (n=544) were genotyped with nine microsatellite loci, and 353 were sequenced for d-loop (384 bp) and ATP (800bp) mitochondrial DNA gene regions. Combined interpretation of microsatellite and mtDNA data identified four genetic stocks of S. semifasciatus: Western Australia, northwest coast of the Northern Territory, Gulf of Carpentaria and the east coast of Queensland. Connectivity among stocks across northern Australia from the Northern Territory to the east coast of Queensland was high, but in contrast, there was a clear genetic break between populations in Western Australia compared to the rest of northern Australia. This indicates a restriction to gene flow possibly associated with suboptimal habitat along the Kimberley coast (northwestern Australia). The appropriate scale of management for this species corresponds to the jurisdictions of the three Australian states, except that the Gulf of Carpentaria stock should be co-managed by authorities in Queensland and Northern Territory.

Fisheries Queensland Long Term Monitoring Program Underwater Visual Census of twenty-four reefs on the Great Barrier Reef between 1999 and 2004.

In 1999, the Department of Employment, Economic Development and Innovation (DEEDI), Fisheries Queensland undertook a new initiative to collect long term monitoring data of various important stocks including reef fish. This data and monitoring manual for the reef fish component of that program which was based on Underwater Visual Census methodology of 24 reefs on the Great Barrier Reef between 1999 and 2004. Data was collected using six 50m x 5m transects at 4 sites on 24 reefs. Benthic cover type was also recorded for 10m of each transect. The attached Access Database contains 5 tables being: SITE DETAILS TABLE Survey year Data entry complete REF survey site ID Site # (1-4) Location (reef name) Site Date (date surveyed) Observer 1 (3 initials to identify who estimated fish lengths and recorded benthic cover) TRANSECT DETAILS Survey ID Transect Number (1-6) Time (the transect was surveyed) Visibility (in metres) Minimum Depth surveyed (m) Maximum Depth surveyed (m) Percent of survey completed (%) Comments SUBSTRATE Survey ID Transect Number (1-6) then % cover of each of eth following categories of benthic cover types Dead Coral Live Coral Soft Coral Rubble Sand Sponge Algae Sea Grass Other COORDINATES (over survey sites) from -14 38.792 to -19 44.233 and from 145 21.507 to 149 55.515 SIGHTINGS ID Survey ID Transect Number (1-6) CAAB Code Scientific Name Reef Fish Length (estimated Fork Length of fish; -1 = unknown or not recorded) Outside Transect (if a fish was observed outside a transect -1 was recorded) Morph Code (F = footballer morph for Plectropomus laevis, S = Spawning colour morph displayed)

Uncertainty in length measurements of live coral trout: implications for compliance and enforcement of minimum legal size limits.

Report on evidence of shrinkage of live coral trout during professional fishing operations on the Great Barrier Reef in 2000. Excel data includes the following fields: Column A. Fish (fish number from 1 -24) Column B. Bin (1-8, container the fish was held in during the experiment) Column C. Measure (1-7, number of the measurement of each fish) Column D. Observer (1 or 2, making the measurement) Column E. Time 2 Column F. Time (time of the day the measurement was made) Column G. FL (Fork Length) Column H. TL (Total Length) Column I. Difference (difference in length between measures) Column J. Order Column K. Temperature (surface water temp under the boat)

Hierarchical Bayesian modelling of plant pest invasions with human-mediated dispersal

Hierarchical Bayesian models can assimilate surveillance and ecological information to estimate both invasion extent and model parameters for invading plant pests spread by people. A reliability analysis framework that can accommodate multiple dispersal modes is developed to estimate human-mediated dispersal parameters for an invasive species. Uncertainty in the observation process is modelled by accounting for local natural spread and population growth within spatial units. Broad scale incursion dynamics are based on a mechanistic gravity model with a Weibull distribution modification to incorporate a local pest build-up phase. The model uses Markov chain Monte Carlo simulations to infer the probability of colonisation times for discrete spatial units and to estimate connectivity parameters between these units. The hierarchical Bayesian model with observational and ecological components is applied to a surveillance dataset for a spiralling whitefly (Aleurodicus dispersus) invasion in Queensland, Australia. The model structure provides a useful application that draws on surveillance data and ecological knowledge that can be used to manage the risk of pest movement.

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.

Southern scallop (Pecten fumatus) microsatellite genotypes (fifteen loci) from eighteen locations along the coastline of south-eastern Australia

Microsatellite genotypes of individual scallops

Water and thermal regimes for field pea in Australia and their implications for breeding.

There is a large gap between the refined approaches to characterise genotypes and the common use of location and season as a coarse surrogate for environmental characterisation of breeding trials. As a framework for breeding, the aim of this paper is quantifying the spatial and temporal patterns of thermal and water stress for field pea in Australia. We compiled a dataset for yield of the cv. Kaspa measured in 185 environments, and investigated the associations between yield and seasonal patterns of actual temperature and modelled water stress. Correlations between yield and temperature indicated two distinct stages. In the first stage, during crop establishment and canopy expansion before flowering, yield was positively associated with minimum temperature. Mean minimum temperature below similar to 7 degrees C suggests that crops were under suboptimal temperature for both canopy expansion and radiation-use efficiency during a significant part of this early growth period. In the second stage, during critical reproductive phases, grain yield was negatively associated with maximum temperature over 25 degrees C. Correlations between yield and modelled water supply/demand ratio showed a consistent pattern with three phases: no correlation at early stages of the growth cycle, a progressive increase in the association that peaked as the crop approached the flowering window, and a progressive decline at later reproductive stages. Using long-term weather records (1957-2010) and modelled water stress for 104 locations, we identified three major patterns of water deficit nation wide. Environment type 1 (ET1) represents the most favourable condition, with no stress during most of the pre-flowering phase and gradual development of mild stress after flowering. Type 2 is characterised by increasing water deficit between 400 degree-days before flowering and 200 degree-days after flowering and rainfall that relieves stress late in the season. Type 3 represents the more stressful condition with increasing water deficit between 400 degree-days before flowering and maturity. Across Australia, the frequency of occurrence was 24% for ET1, 32% for ET2 and 43% for ET3, highlighting the dominance of the most stressful condition. Actual yield averaged 2.2 t/ha for ET1, 1.9 t/ha for ET2 and 1.4 t/ha for ET3, and the frequency of each pattern varied substantially among locations. Shifting from a nominal (i.e. location and season) to a quantitative (i.e. stress type) characterisation of environments could help improving breeding efficiency of field pea in Australia.

Queensland Coastal Wetland Resources: Round Hill Head to Tin Can Inlet

The wetland resources of the Queensland coastline have been mapped by the Resource Condition and Trend Unit, Fisheries Group, Department of Primary Industries Queensland. This process is being undertaken in order to provide a baseline dataset for Fish Habitat Area (FHA) declaration, Ramsar site nomination and continued monitoring of these important fish habitats. This report summarises the results of the mapping undertaken from Round Hill Head to Tin Can Inlet. The study was undertaken in order to: 1. document and map the coastal wetland communities from Round Hill Head (24°S) to Tin Can Inlet (26°S); 2. document levels of existing disturbance to and protection of the wetlands; 3. examine existing recreational and commercial fisheries in the region; and 4. evaluate the conservation values of the areas investigated from the viewpoint of fisheries productivity and as habitat for important and/or threatened species.

Raw data for "Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson)."

Abstract of Macbeth, G. M., Broderick, D., Buckworth, R. & Ovenden, J. R. (In press, Feb 2013). Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson). G3: Genes, Genomes and Genetics. Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. But, pitfalls are predicted when the effective size is large, as estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes around 5000. This was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest standard deviation of Ne estimates occurred when low frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or non-breeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were above 24,000.

Long-term repeated burning reduces Lantana camara regeneration in a dry eucalypt forest

Previous short-term studies predict that the use of fire to manage lantana (Lantana camara) may promote its abundance. We tested this prediction by examining long-term recruitment patterns of lantana in a dry eucalypt forest in Australia from 1959 to 2007 in three fire frequency treatments: repeated annual burning, repeated triennial burning and long unburnt. The dataset was divided into two periods (1959–1972, 1974–2007) due to logging that occurred at the study site between 1972 and 1974 and the establishment of the triennial burn treatment in 1973. Our results showed that repeated burning decreased lantana regeneration under an annual burn regime in the pre- and post-logging periods and maintained low levels of regeneration in the triennial burn compartment during the post-logging period. In the absence of fire, lantana recruitment exhibited a dome-shaped response over time, with the total population peaking in 1982 before declining to 2007. In addition to fire regime, soil pH and carbon to nitrogen ratio, the density of taller conspecifics and the interaction between rainfall and fire regime were found to influence lantana regeneration change over time. The results suggest that the reported positive association between fire disturbance and abundance of lantana does not hold for all forest types and that fire should be considered as part of an integrated weed management strategy for lantana in more fire-tolerant ecosystems.

A Bayesian hurdle model for analysis of an insect resistance monitoring database

Motivated by the analysis of the Australian Grain Insect Resistance Database (AGIRD), we develop a Bayesian hurdle modelling approach to assess trends in strong resistance of stored grain insects to phosphine over time. The binary response variable from AGIRD indicating presence or absence of strong resistance is characterized by a majority of absence observations and the hurdle model is a two step approach that is useful when analyzing such a binary response dataset. The proposed hurdle model utilizes Bayesian classification trees to firstly identify covariates and covariate levels pertaining to possible presence or absence of strong resistance. Secondly, generalized additive models (GAMs) with spike and slab priors for variable selection are fitted to the subset of the dataset identified from the Bayesian classification tree indicating possibility of presence of strong resistance. From the GAM we assess trends, biosecurity issues and site specific variables influencing the presence of strong resistance using a variable selection approach. The proposed Bayesian hurdle model is compared to its frequentist counterpart, and also to a naive Bayesian approach which fits a GAM to the entire dataset. The Bayesian hurdle model has the benefit of providing a set of good trees for use in the first step and appears to provide enough flexibility to represent the influence of variables on strong resistance compared to the frequentist model, but also captures the subtle changes in the trend that are missed by the frequentist and naive Bayesian models. © 2014 Springer Science+Business Media New York.

Rapid genome wide mapping of phosphine resistance loci by a simple regional averaging analysis in the red flour beetle, Tribolium castaneum

Background Next-generation sequencing technology is an important tool for the rapid, genome-wide identification of genetic variations. However, it is difficult to resolve the ‘signal’ of variations of interest and the ‘noise’ of stochastic sequencing and bioinformatic errors in the large datasets that are generated. We report a simple approach to identify regional linkage to a trait that requires only two pools of DNA to be sequenced from progeny of a defined genetic cross (i.e. bulk segregant analysis) at low coverage (<10×) and without parentage assignment of individual SNPs. The analysis relies on regional averaging of pooled SNP frequencies to rapidly scan polymorphisms across the genome for differential regional homozygosity, which is then displayed graphically. Results Progeny from defined genetic crosses of Tribolium castaneum (F4 and F19) segregating for the phosphine resistance trait were exposed to phosphine to select for the resistance trait while the remainders were left unexposed. Next generation sequencing was then carried out on the genomic DNA from each pool of selected and unselected insects from each generation. The reads were mapped against the annotated T. castaneum genome from NCBI (v3.0) and analysed for SNP variations. Since it is difficult to accurately call individual SNP frequencies when the depth of sequence coverage is low, variant frequencies were averaged across larger regions. Results from regional SNP frequency averaging identified two loci, tc_rph1 on chromosome 8 and tc_rph2 on chromosome 9, which together are responsible for high level resistance. Identification of the two loci was possible with only 5-7× average coverage of the genome per dataset. These loci were subsequently confirmed by direct SNP marker analysis and fine-scale mapping. Individually, homozygosity of tc_rph1 or tc_rph2 results in only weak resistance to phosphine (estimated at up to 1.5-2.5× and 3-5× respectively), whereas in combination they interact synergistically to provide a high-level resistance >200×. The tc_rph2 resistance allele resulted in a significant fitness cost relative to the wild type allele in unselected beetles over eighteen generations. Conclusion We have validated the technique of linkage mapping by low-coverage sequencing of progeny from a simple genetic cross. The approach relied on regional averaging of SNP frequencies and was used to successfully identify candidate gene loci for phosphine resistance in T. castaneum. This is a relatively simple and rapid approach to identifying genomic regions associated with traits in defined genetic crosses that does not require any specialised statistical analysis.