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.

Selecting for increased barley grain size

In this study, 120–144 commercial varieties and breeding lines were assessed for grain size attributes including plump grain (>2.8 mm) and retention (>2.5 mm+>2.8 mm). Grain samples were produced from replicated trials at 25 sites across four years. Climatic conditions varied between years as well as between sites. Several of the trial sites were irrigated while the remaining were produced under dryland conditions. A number of the dryland sites suffered severe drought stress. The grain size data was analysed for genetic (G), environmental (E) and genotype by environment (G×E) interactions. All analyses included maturity as a covariate. The genetic effect on grain size was greater than environmental or maturity effects despite some sites suffering terminal moisture stress. The model was used to calculate heritability values for each site used in the study. These values ranged from 89 to 98% for plump grain and 88 to 96% for retention. The results demonstrated that removing the sources of non-heritable variation, such as maturity and field effects, can improve genetic estimates of the retention and plump grain fractions. By partitioning all variance components, and thereby having more robust estimates of genetic differences, plant breeders can have greater confidence in selecting barley genotypes which maintain large, stable grain size across a range of environments.

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.

Optimising between and within family selection for harvest weight in prawns with restricted harvest size.

Manual grading of prawns restricts the number that can be harvested. A restricted harvest size places a limit on the opposing within family and between family sources of selection pressure. A simulation study with inbreeding constrained at 0.5% per generation, a harvest size of 2000, heritability of 0.3, common family environmental effect of 0.1, indicates that maximum response to selection could be achieved with as few as 40 families. Increasing the number of families above 80 may reduce total selection response. It is important to be aware that increasing the number of families may not always yield a greater genetic response.

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)

Bayesian mixture model estimation of aerosol particle size distributions

In this paper, we examine approaches to estimate a Bayesian mixture model at both single and multiple time points for a sample of actual and simulated aerosol particle size distribution (PSD) data. For estimation of a mixture model at a single time point, we use Reversible Jump Markov Chain Monte Carlo (RJMCMC) to estimate mixture model parameters including the number of components which is assumed to be unknown. We compare the results of this approach to a commonly used estimation method in the aerosol physics literature. As PSD data is often measured over time, often at small time intervals, we also examine the use of an informative prior for estimation of the mixture parameters which takes into account the correlated nature of the parameters. The Bayesian mixture model offers a promising approach, providing advantages both in estimation and inference.

The effect of 3-5-6 TPA on fruit drop and fruit size in the lychee (Litchi Chinensis) cultivars 'Fay Zee Siu' (feizixiao) , 'Kaimana', 'Kwai Mai Pink', 'Souey Tung' and 'Tai So' (Mauritus)

Fruit drop can cause major yield losses in Australian lychee orchards, the severity varying with cultivar and season. Research in China, South Africa and Israel has demonstrated the potential for synthetic auxins used as foliar sprays to reduce fruit drop in lychee. Trials tested the efficacy of the synthetic auxin 3-5-6 trichloro-2-phridyl-oxyacetic acid (TPA) applied as a foliar spray at 50 ppm on fruit drop and fruit size on the cultivars ‘Fay Zee Siu’, ‘Kaimana’, ‘Kwai Mai Pink’, ‘Souey Tung’ and ‘Tai So’. TPA reduced fruit drop when applied to fruit greater than 12 mm in length but increased fruit drop when fruit were smaller. Fruit size at the time of application had less effect on the response than the level of natural fruit drop. When natural fruit drop was high, TPA significantly reduced it; by up to 18.7 in ‘Fay Zee Siu’, 37.1 in ‘Kaimana’, 39.8 in ‘Kwai Mai Pink’, 15.1 in ‘Souey Tung’ and 7.7 in ‘Tai So’. TPA was less effective when natural fruit drop was low. TPA increased the number of large fruit and frequently increased the number of small fruit at harvest. The small fruit were associated with an increase in the retention of fruit with poorly developed (chicken tongue) seed. Average fruit size was generally larger (up to 12.7 in ‘Souey Tung’ and 22 in ‘Tai So’) with TPA applications.

What are the driving factors influencing the size distribution of airborne synthetic clay particles emitted from a jet milling process?

In the field of workplace air quality, measuring and analyzing the size distribution of airborne particles to identify their sources and apportion their contribution has become widely accepted, however, the driving factors that influence this parameter, particularly for nanoparticles (< 100 nm), have not been thoroughly determined. Identification of driving factors, and in turn, general trends in size distribution of emitted particles would facilitate the prediction of nanoparticles’ emission behavior and significantly contribute to their exposure assessment. In this study, a comprehensive analysis of the particle number size distribution data, with a particular focus on the ultrafine size range of synthetic clay particles emitted from a jet milling machine was conducted using the multi-lognormal fitting method. The results showed relatively high contribution of nanoparticles to the emissions in many of the tested cases, and also, that both surface treatment and feed rate of the machine are significant factors influencing the size distribution of the emitted particles of this size. In particular, applying surface treatments and increasing the machine feed rate have the similar effect of reducing the size of the particles, however, no general trend was found in variations of size distribution across different surface treatments and feed rates. The findings of our study demonstrate that for this process and other activities, where no general trend is found in the size distribution of the emitted airborne particles due to dissimilar effects of the driving factors, each case must be treated separately in terms of workplace exposure assessment and regulations.

Factors influencing seasonal and monthly changes in the group size of chital or axis deer in southern India

Chital or axis deer (Axis axis) form fluid groups that change in size temporally and in relation to habitat. Predictions of hypotheses relating animal density, rainfall, habitat structure, and breeding seasonality, to changes in chital group size were assessed simultaneously using multiple regression models of monthly data collected over a 2 yr period in Guindy National Park, in southern India. Over 2,700 detections of chital groups were made during four seasons in three habitats (forest, scrubland and grassland). In scrubland and grassland, chital group size was positively related to animal density, which increased with rainfall. This suggests that in these habitats, chital density increases in relation to food availability, and group sizes increase due to higher encounter rate and fusion of groups. The density of chital in forest was inversely related to rainfall, but positively to the number of fruiting tree species and availability of fallen litter, their forage in this habitat. There was little change in mean group size in the forest, although chital density more than doubled during the dry season and summer. Dispersion of food items or the closed nature of the forest may preclude formation of larger groups. At low densities, group sizes in all three habitats were similar. Group sizes increased with chital density in scrubland and grassland, but more rapidly in the latter—leading to a positive relationship between openness and mean group size at higher densities. It is not clear, however, that this relationship is solely because of the influence of habitat structure. The rutting index (monthly percentage of adult males in hard antler) was positively related to mean group size in forest and scrubland, probably reflecting the increase in group size due to solitary males joining with females during the rut. The fission-fusion system of group formation in chital is thus interactively influenced by several factors. Aspects that need further study, such as interannual variability, are highlighted.

Evolution of polyandry by reduction in progeny number variance in structured populations

When there is a variation in the quality of males in a population, multiple mating can lead to an increase in the genetic fitness of a female by reducing the variance of the progeny number. The extent of selective advantage obtainable by this process is investigated for a population subdivided into structured demes. It is seen that for a wide range of model parameters (deme size, distribution of male quality, local resource level), multiple mating leads to a considerable increase in the fitness. Frequency-dependent selection or a stable coexistence between polyandry and monandry can also result when the possible costs involved in multiple mating are taken into account.

On Finding the Natural Number of Topics with Latent Dirichlet Allocation: Some Observations

It is important to identify the ``correct'' number of topics in mechanisms like Latent Dirichlet Allocation(LDA) as they determine the quality of features that are presented as features for classifiers like SVM. In this work we propose a measure to identify the correct number of topics and offer empirical evidence in its favor in terms of classification accuracy and the number of topics that are naturally present in the corpus. We show the merit of the measure by applying it on real-world as well as synthetic data sets(both text and images). In proposing this measure, we view LDA as a matrix factorization mechanism, wherein a given corpus C is split into two matrix factors M-1 and M-2 as given by C-d*w = M1(d*t) x Q(t*w).Where d is the number of documents present in the corpus anti w is the size of the vocabulary. The quality of the split depends on ``t'', the right number of topics chosen. The measure is computed in terms of symmetric KL-Divergence of salient distributions that are derived from these matrix factors. We observe that the divergence values are higher for non-optimal number of topics - this is shown by a `dip' at the right value for `t'.

Design of speaker identification schemes for large number of speakers (A)

Design of speaker identification schemes for a small number of speakers (around 10) with a high degree of accuracy in controlled environment is a practical proposition today. When the number of speakers is large (say 50–100), many of these schemes cannot be directly extended, as both recognition error and computation time increase monotonically with population size. The feature selection problem is also complex for such schemes. Though there were earlier attempts to rank order features based on statistical distance measures, it has been observed only recently that the best two independent measurements are not the same as the combination in two's for pattern classification. We propose here a systematic approach to the problem using the decision tree or hierarchical classifier with the following objectives: (1) Design of optimal policy at each node of the tree given the tree structure i.e., the tree skeleton and the features to be used at each node. (2) Determination of the optimal feature measurement and decision policy given only the tree skeleton. Applicability of optimization procedures such as dynamic programming in the design of such trees is studied. The experimental results deal with the design of a 50 speaker identification scheme based on this approach.

H I deficiency in cluster spiral galaxies - Dependence on galaxy size

From the available H I data on spiral galaxies in three rich Abell clusters and the Virgo Cluster, it is shown that galaxies with medium to large optical sizes tend to be more severely deficient in atomic hydrogen than the small galaxies. This is so both in terms of the fractional number of galaxies that are deficient and the amount of gas lost by a galaxy. The fraction of H I-deficient galaxies increases with size over most of the size range, saturating or dropping only for the largest galaxies. A comparative study is made of various currently accepted gas removal mechanisms, namely those which are a result of galaxy-intracluster medium interactions, e.g., ram pressure stripping, as well as those due to galaxy-galaxy interactions, i.e., collisions and tidal interactions. It is shown that, with the exception of tidal interactions, all of these mechanisms would produce a size dependence in H I deficiency that is the opposite of that observed. That is, the gas in the largest galaxies would be the least affected by these mechanisms. However, if there is significant mass segregation, these processes may give the trends observed in the size dependence of H I deficiency.

Fatigue crack propagation model and size effect in concrete using dimensional analysis

It is well known that fatigue in concrete causes excessive deformations and cracking leading to structural failures. Due to quasi-brittle nature of concrete and formation of a fracture process zone, the rate of fatigue crack growth depends on a number of parameters, such as, the tensile strength, fracture toughness, loading ratio and most importantly the structural size. In this work, an analytical model is proposed for estimating the fatigue crack growth in concrete by using the concepts of dimensional analysis and including the above parameters. Knowing the governed and the governing parameters of the physical problem and by using the concepts of self-similarity, a relationship is obtained between different parameters involved. It is shown that the proposed fatigue law is able to capture the size effect in plain concrete and agrees well with different experimental results. Through a sensitivity analysis, it is shown that the structural size plays a dominant role followed by loading ratio and the initial crack length in fatigue crack propagation. (C) 2010 Elsevier Ltd. All rights reserved.

Rainbow Connection Number and Radius

Rainbow connection number, rc(G), of a connected graph G is the minimum number of colours needed to colour its edges, so that every pair of vertices is connected by at least one path in which no two edges are coloured the same. In this note we show that for every bridgeless graph G with radius r, rc(G) <= r(r+2). We demonstrate that this bound is the best possible for rc(G) as a function of r, not just for bridgeless graphs, but also for graphs of any stronger connectivity. It may be noted that, for a general 1-connected graph G, rc(G) can be arbitrarily larger than its radius (K_{1,n} for instance). We further show that for every bridgeless graph G with radius r and chordality (size of a largest induced cycle) k, rc(G) <= rk. Hitherto, the only reported upper bound on the rainbow connection number of bridgeless graphs is 4n/5 - 1, where n is order of the graph [Caro et al., 2008]