Estimating and evaluating the statistics of gapped local-alignment scores

We present a novel maximum-likelihood-based algorithm for estimating the distribution of alignment scores from the scores of unrelated sequences in a database search. Using a new method for measuring the accuracy of p-values, we show that our maximum-likelihood-based algorithm is more accurate than existing regression-based and lookup table methods. We explore a more sophisticated way of modeling and estimating the score distributions (using a two-component mixture model and expectation maximization), but conclude that this does not improve significantly over simply ignoring scores with small E-values during estimation. Finally, we measure the classification accuracy of p-values estimated in different ways and observe that inaccurate p-values can, somewhat paradoxically, lead to higher classification accuracy. We explain this paradox and argue that statistical accuracy, not classification accuracy, should be the primary criterion in comparisons of similarity search methods that return p-values that adjust for target sequence length.

Vibration instability in rolling mills: Modeling and experimental results

Models for the occurrence of the vibrational instability during rolling known as third octave chatter are presented and discussed. An analysis of rolling mill chatter was performed for the purpose of identifying characteristics of the vibrations and to determine any dependency on the rolling schedule. In particular, a stability criterion for the critical rolling speed is used to predict the maximum rolling speed without chatter instability on schedules from a 5 stand tandem mill rolling thin steel product. The results correlate well with measurements of critical speed occurring on the mill using a vibration monitor: This research provides significant insights into the chatter phenomena and has been used to investigate control methods for suppression of the instability.

Application of spatial analysis techniques to adjust for fertility trends and identify interplot competition in early stage sugarcane selection trials

Most sugarcane breeding programs in Australia use large unreplicated trials to evaluate clones in the early stages of selection. Commercial varieties that are replicated provide a method of local control of soil fertility. Although such methods may be useful in detecting broad trends in the field, variation often occurs on a much smaller scale. Methods such as spatial analysis adjust a plot for variability by using information from immediate neighbours. These techniques are routinely used to analyse cereal data in Australia and have resulted in increased accuracy and precision in the estimates of variety effects. In this paper, spatial analyses in which the variability is decomposed into local, natural, and extraneous components are applied to early selection trials in sugarcane. Interplot competition in cane yield and trend in sugar content were substantial in many of the trials and there were often large differences in the selections between the spatial and current method used by the Bureau of Sugar Experiment Stations. A joint modelling approach for tonnes sugar per hectare in response to fertility trends and interplot competition is recommended.

Yield response of rice (Oryza sativa L.) genotypes to different types of drought under rainfed lowlands - Part 2. Selection of drought resistant genotypes

Drought frequently reduces grain yield of rainfed lowland rice. A series of experiments were conducted in drought-prone northeast Thailand to study the magnitude and consistency of yield responses of diverse, rainfed lowland rice genotypes to drought stress environments and to examine ways to identify genotypes that confer drought resistance. One hundred and twenty-eight genotypes were grown under non-stress and four different types of drought stress conditions. The relationship of genotypic variation in yield under drought conditions to genetic yield potential, flowering time and flowering delay, and to a drought response index (DRI) that removed the effect of potential yield and flowering time on yield under stress was examined. Drought stress that developed prior to flowering generally delayed the time of flowering of genotypes, and the delay in flowering was negatively associated with grain yield, fertile panicle percentage and filled grain percentage. Genotypes with a longer delay in flowering time had extracted more water during the early drought period, and as a consequence, had higher water deficits. They were consistently associated with a larger yield reduction under drought and in one experiment with a smaller DRI. Genotypes, however, responded differently to the different drought stress conditions and there was no consistency in the DRI estimates for the different genotypes across the drought stress experiments. The results indicate that with the use of irrigated-control and drought test environments, genotypes with drought resistance can be identified by using DRI or delay in flowering. However, selections will differ depending on the type of drought condition. The inconsistency of the estimates in DRI and flowering delay across different drought conditions reflects the nature of the large genotype-by-environment interactions observed for grain yield under various types of drought in rainfed lowland conditions. (C), 2002 Elsevier Science B.V. All rights reserved.

Evaluation of various pre-treatments for the dehydration of banana and selection of suitable drying models

The thin-layer drying behaviour of bananas in a beat pump dehumidifier dryer was examined. Four pre-treatments (blanching, chilling, freezing and combined blanching and freezing) were applied to the bananas, which were dried at 50 degreesC with an air velocity of 3.1 m s(-1) and with the relative humidity of the inlet air of 10-35%. Three drying models, the simple model, the two-term exponential model and the Page model were examined. All models were evaluated using three statistical measures, correlation coefficient, root means square error, and mean absolute percent error. Moisture diffusivity was calculated based on the diffusion equation for an infinite cylindrical shape using the slope method. The rate of drying was higher for the pre-treatments involving freezing. The sample which was blanched only did not show any improvement in drying rate. In fact, a longer drying time resulted due to water absorption during blanching. There was no change in the rate for the chilled sample compared with the control. While all models closely fitted the drying data, the simple model showed greatest deviation from the experimental results. The two-term exponential model was found to be the best model for describing the drying curves of bananas because its parameters represent better the physical characteristics of the drying process. Moisture diffusivities of bananas were in the range 4.3-13.2 x 10(-10) m(2)s(-1). (C) 2002 Published by Elsevier Science Ltd.

Using crop simulation to generate genotype by environment interaction effects for sorghum in water-limited environments

Multi-environment trials (METs) used to evaluate breeding lines vary in the number of years that they sample. We used a cropping systems model to simulate the target population of environments (TPE) for 6 locations over 108 years for 54 'near-isolines' of sorghum in north-eastern Australia. For a single reference genotype, each of 547 trials was clustered into 1 of 3 'drought environment types' (DETs) based on a seasonal water stress index. Within sequential METs of 2 years duration, the frequencies of these drought patterns often differed substantially from those derived for the entire TPE. This was reflected in variation in the mean yield of the reference genotype. For the TPE and for 2-year METs, restricted maximum likelihood methods were used to estimate components of genotypic and genotype by environment variance. These also varied substantially, although not in direct correlation with frequency of occurrence of different DETs over a 2-year period. Combined analysis over different numbers of seasons demonstrated the expected improvement in the correlation between MET estimates of genotype performance and the overall genotype averages as the number of seasons in the MET was increased.

Genetic constraints on the evolution of mate recognition under natural selection

Field populations of Drosophila serrata display reproductive character displacement in cuticular hydrocarbons (CHCs) when sympatric with Drosophila birchii. We have previously shown that the naturally occurring pattern of reproductive character displacement can be experimentally replicated by exposing field allopatric populations of D. serrata to experimental sympatry with D. birchii. Here, we tested whether the repeated evolution of reproductive character displacement in natural and experimental populations was a consequence of genetic constraints on the evolution of CHCs. The genetic variance-covariance (G) matrices for CHCs were determined for populations of D. serrata that had evolved in either the presence or absence of D. birchii under field and experimental conditions. Natural selection on mate recognition under both field and experimental sympatric conditions increased the genetic variance in CHCs consistent with a response to selection based on rare alleles. A close association between G eigenstructure and the eigenstructure of the phenotypic divergence (D) matrix in natural and experimental populations suggested that G matrix eigenstructure may have determined the direction in which reproductive character displacement evolved during the reinforcement of mate recognition.

Species richness in agamid lizards: chance, body size, sexual selection or ecology?

Why does species richness vary so greatly across lineages? Traditionally, variation in species richness has been attributed to deterministic processes, although it is equally plausible that it may result from purely stochastic processes. We show that, based on the best available phylogenetic hypothesis, the pattern of cladogenesis among agamid lizards is not consistent with a random model, with some lineages having more species, and others fewer species, than expected by chance. We then use phylogenetic comparative methods to test six types of deterministic explanation for variation in species richness: body size, life history, sexual selection, ecological generalism, range size and latitude. Of eight variables we tested, only sexual size dimorphism and sexual dichromatism predicted species richness. Increases in species richness are associated with increases in sexual dichromatism but reductions in sexual size dimorphism. Consistent with recent comparative studies, we find no evidence that species richness is associated with small body size or high fecundity. Equally, we find no evidence that species richness covaries with ecological generalism, latitude or range size.

High irradiance increases organogenesis in friable callus of Caustis blakei Kuk. (Cyperaceae)

Caustis blakei is an attractive cut foliage plant harvested from the wild in Australia and marketed under the name of koala fern. Previous attempts to propagate large numbers of this plant have been unsuccessful. The effect of four light irradiances on organogenesis from compact and friable callus of C. blakei was studied for 21 wk. Both callus types produced numerous primordial shoots but many failed to develop into green plantlets. However, significantly more primordial shoots and green plantlets developed on the friable callus than on the compact callus, and significantly more green plantlets were regenerated under the higher photon irradiances of 200 and 300 mumol m(-2) s(-1) than under the lower irradiances of 100 and 150 mumol m(-2) s(-1). The compact callus produced its maximum number of green plantlets early in the experiment (after 9 wk), while the friable callus continued to produce primordial shoots and green plantlets throughout the period of the experiment, and reached its maximum production of green plantlets at 21 wk under the irradiance of 300 mumol m(-2) s(-1). Organogenesis from friable callus under high irradiance (300 mumol m(-2) s(-1)) offers an efficient propagation method for C. blakei.

Mean anisotropy of homogeneous Gaussian random fields and anisotropic norms of linear translation-invariant operators on multidimensional integer lattices

Sensitivity of output of a linear operator to its input can be quantified in various ways. In Control Theory, the input is usually interpreted as disturbance and the output is to be minimized in some sense. In stochastic worst-case design settings, the disturbance is considered random with imprecisely known probability distribution. The prior set of probability measures can be chosen so as to quantify how far the disturbance deviates from the white-noise hypothesis of Linear Quadratic Gaussian control. Such deviation can be measured by the minimal Kullback-Leibler informational divergence from the Gaussian distributions with zero mean and scalar covariance matrices. The resulting anisotropy functional is defined for finite power random vectors. Originally, anisotropy was introduced for directionally generic random vectors as the relative entropy of the normalized vector with respect to the uniform distribution on the unit sphere. The associated a-anisotropic norm of a matrix is then its maximum root mean square or average energy gain with respect to finite power or directionally generic inputs whose anisotropy is bounded above by a≥0. We give a systematic comparison of the anisotropy functionals and the associated norms. These are considered for unboundedly growing fragments of homogeneous Gaussian random fields on multidimensional integer lattice to yield mean anisotropy. Correspondingly, the anisotropic norms of finite matrices are extended to bounded linear translation invariant operators over such fields.

On continuous selection problems for multivalued mappings with the local intersection property in hyperconvex metric spaces

In the paper we present two continuous selection theorems in hyperconvex metric spaces and apply these to study xed point and coincidence point problems as well as variational inequality problems in hyperconvex metric spaces.

A structural basis for the selection of dominant alpha beta T cell receptors in antiviral immunity

We have examined the basis for immunodominant or public TCR usage in an antiviral CTL response. Residues encoded by each of the highly selected genetic elements of an immunodominant clonotype recognizing Epstein-Barr virus were critical to the antigen specificity of the receptor. Upon recognizing antigen the immunodominant TCR undergoes extensive conformational changes in the complementarity determining regions (CDRs), including the disruption of the canonical structures of the germline-encoded CDR1alpha and CDR2alpha loops to produce an enhanced fit with the HLA-peptide complex. TCR ligation induces conformational changes in the TCRalpha constant domain thought to form part of the docking site for CD3epsilon. These findings indicate that TCR immunodominance is associated with structural properties conferring receptor specificity and suggest a novel structural link between TCR ligation and intracellular signaling.

Rock fracture mean trace length estimation and confidence interval calculation using maximum likelihood methods

There has been a resurgence of interest in the mean trace length estimator of Pahl for window sampling of traces. The estimator has been dealt with by Mauldon and Zhang and Einstein in recent publications. The estimator is a very useful one in that it is non-parametric. However, despite some discussion regarding the statistical distribution of the estimator, none of the recent works or the original work by Pahl provide a rigorous basis for the determination a confidence interval for the estimator or a confidence region for the estimator and the corresponding estimator of trace spatial intensity in the sampling window. This paper shows, by consideration of a simplified version of the problem but without loss of generality, that the estimator is in fact the maximum likelihood estimator (MLE) and that it can be considered essentially unbiased. As the MLE, it possesses the least variance of all estimators and confidence intervals or regions should therefore be available through application of classical ML theory. It is shown that valid confidence intervals can in fact be determined. The results of the work and the calculations of the confidence intervals are illustrated by example. (C) 2003 Elsevier Science Ltd. All rights reserved.

Signal trait sexual dimorphism and mutual sexual selection in Drosophila serrata

The evolution of sexual dimorphism may occur when natural and sexual selection result in different optimum trait values for males and females. Perhaps the most prominent examples of sexual dimorphism occur in sexually selected traits, for which males usually display exaggerated trait levels, while females may show reduced expression of the trait. In some species, females also exhibit secondary sexual traits that may either be a consequence of a correlated response to sexual selection on males or direct sexual selection for female secondary sexual traits. In this experiment, we simultaneously measure the intersex genetic correlations and the relative strength of sexual selection on males and females for a set of cuticular hydrocarbons in Drosophila serrata. There was significant directional sexual selection on both male and female cuticular hydrocarbons: the strength of sexual selection did not differ among the sexes but males and females preferred different cuticular hydrocarbons. In contrast with many previous studies of sexual dimorphism, intersex genetic correlations were low. The evolution of sexual dimorphism in D. serrata appears to have been achieved by sex-limited expression of traits controlled by genes on the X chromosome and is likely to be in its final stages.