Comparative phylogeographic summary statistics for testing simultaneous vicariance

Testing for simultaneous vicariance across comparative phylogeographic data sets is a notoriously difficult problem hindered by mutational variance, the coalescent variance, and variability across pairs of sister taxa in parameters that affect genetic divergence. We simulate vicariance to characterize the behaviour of several commonly used summary statistics across a range of divergence times, and to characterize this behaviour in comparative phylogeographic datasets having multiple taxon-pairs. We found Tajima's D to be relatively uncorrelated with other summary statistics across divergence times, and using simple hypothesis testing of simultaneous vicariance given variable population sizes, we counter-intuitively found that the variance across taxon pairs in Nei and Li's net nucleotide divergence (pi(net)), a common measure of population divergence, is often inferior to using the variance in Tajima's D across taxon pairs as a test statistic to distinguish ancient simultaneous vicariance from variable vicariance histories. The opposite and more intuitive pattern is found for testing more recent simultaneous vicariance, and overall we found that depending on the timing of vicariance, one of these two test statistics can achieve high statistical power for rejecting simultaneous vicariance, given a reasonable number of intron loci (> 5 loci, 400 bp) and a range of conditions. These results suggest that components of these two composite summary statistics should be used in future simulation-based methods which can simultaneously use a pool of summary statistics to test comparative the phylogeographic hypotheses we consider here.

DNA markers linked to yield, yield components, and morphological traits in autotetraploid lucerne (Medicago sativa L.)

We have mapped and identifed DNA markers linked to morphology, yield, and yield components of lucerne, using a backcross population derived from winter-active parents. The high-yielding and recurrent parent, D, produced individual markers that accounted for up to 18% of total yield over 6 harvests, at Gatton, south-eastern Queensland. The same marker, AC/TT8, was consistently identified at each individual harvest, and in individual harvests accounted for up to 26% of the phenotypic variation for yield. This marker was located in linkage group 2 of the D map, and several other markers positively associated with yield were consistently identified in this linkage group. Similarly, markers negatively associated with yield were consistently identified in the W116 map, W116 being the low-yielding parent. Highly significant positive correlations were observed between total yield and yield for harvests 1-6, and between total yield and stem length, tiller number, leaf yield/plant, leaf yield/5 stems, stem yield/plant, and stem yield/5 stems. Highly significant QTL were located for all these characters as well as for leaf shape and pubescence.

Structure and ultrastructure of the silk glands and spinneret of Helicoverpa armigera (Hubner) (Lepidoptera : Noctuidae)

This study provides comprehensive documentation of silk production in the pest moth Helicoverpa armigera from gland secretion to extrusion of silk thread. The structure of the silk glands, accessory structures and extrusion apparatus are reported. The general schema of the paired silk glands follows that found for Lepidoptera. Morphology of the duct, silk press, muscle attachments and spigot are presented as a three-dimensional reconstruction and the cuticular crescent-shaped profile of the silk press is demonstrated in both open and closed forms with attendant muscle blocks, allowing advances in our knowledge of how the silk press functions to regulate the extrusion of silk. Growth of the spigot across instars is documented showing a distinctive developmental pattern for this extrusion device. Its shape and structure are related to use and load-bearing activity. (c) 2005 Elsevier Ltd. All rights reserved.

The effect of box shape on the dynamic properties of proteins simulated under periodic boundary conditions

The effect of the box shape on the dynamic behavior of proteins simulated under periodic boundary conditions is evaluated. In particular, the influence of simulation boxes defined by the near-densest lattice packing (NDLP) in conjunction with rotational constraints is compared to that of standard box types without these constraints. Three different proteins of varying size, shape, and secondary structure content were examined in the study. The statistical significance of differences in RMSD, radius of gyration, solvent-accessible surface, number of hydrogen bonds, and secondary structure content between proteins, box types, and the application or not of rotational constraints has been assessed. Furthermore, the differences in the collective modes for each protein between different boxes and the application or not of rotational constraints have been examined. In total 105 simulations were performed, and the results compared using a three-way multivariate analysis of variance (MANOVA) for properties derived from the trajectories and a three-way univariate analysis of variance (ANOVA) for collective modes. It is shown that application of roto-translational constraints does not have a statistically significant effect on the results obtained from the different simulations. However, the choice of simulation box was found to have a small (5-10%), but statistically significant effect on the behavior of two of the three proteins included in the study. (c) 2005 Wiley Periodicals, Inc.

Vehicle routing and crew scheduling for metropolitan mail distribution at Australia Post

This paper presents a new multi-depot combined vehicle and crew scheduling algorithm, and uses it, in conjunction with a heuristic vehicle routing algorithm, to solve the intra-city mail distribution problem faced by Australia Post. First we describe the Australia Post mail distribution problem and outline the heuristic vehicle routing algorithm used to find vehicle routes. We present a new multi-depot combined vehicle and crew scheduling algorithm based on set covering with column generation. The paper concludes with a computational investigation examining the affect of different types of vehicle routing solutions on the vehicle and crew scheduling solution, comparing the different levels of integration possible with the new vehicle and crew scheduling algorithm and comparing the results of sequential versus simultaneous vehicle and crew scheduling, using real life data for Australia Post distribution networks.

Characterisation of the signalment, clinical and survival characteristics of 41 cats with mast cell neoplasia

Mast cell tumours (MCTs) are relatively common tumours of cats, and are the second most common cutaneous tumours in cats in the USA. While the primary splenic form of the disease is far less common, it is usually associated with more severe clinical signs. Signalment, clinical and survival characteristics of mast cell neoplasia were characterised in 41 cats. The most common tumour location was cutaneous/ subcutaneous head and trunk. Stage la was the most common tumour stage at first diagnosis (n = 20), followed by stage 4 (both stage 4a and stage 4b; n = 10). Of 22 cats that underwent excisional biopsy, mast cell neoplasia recurred in four cats during the study period. Three of the 41 cats presented with simultaneous cutaneous and either splenic or lymph node tumours. A comparison between cats with only cutaneous tumours (n = 30) and those with tumours involving the spleen or lymph nodes (n = 11) showed longer survival times for the cutaneous-only group (P = 0.031). Twelve of the 41 cats died of mast cell neoplasia during the study period. When a subgroup of cats with only cutaneous tumours (no lymph node or visceral involvement) were divided according to whether there were multiple (five or more) tumours (n = 6) or a single tumour (n = 19), cats with single tumours survived longer than those with multiple tumours (P = 0.001). Solitary cutaneous feline MCTs without spread to the lymph nodes usually manifest as benign disease with a relatively protracted course. However, multiple cutaneous tumours, recurrent tumours and primary splenic disease should receive a guarded prognosis due to the relatively short median survival times associated with these forms of the disease. (C) 2006 ESFM and AAFR Published by Elsevier Ltd. All rights reserved.

Modelling broccoli development, yield and quality

Broccoli is a vegetable crop of increasing importance in Australia, particularly in south-east Queensland and farmers need to maintain a regular supply of good quality broccoli to meet the expanding market. A predictive model of ontogeny, incorporating climatic data including frost risk, would enable farmers to predict harvest maturity date and select appropriate cultivar – sowing date combinations. To develop procedures for predicting ontogeny, yield and quality, field studies using three cultivars, ‘Fiesta’, ‘Greenbelt’ and ‘Marathon’, were sown on eight dates from 11 March to 22 May 1997, and grown under natural and extended (16 h) photoperiods at the University of Queensland, Gatton Campus. Cultivar, rather than the environment, mainly determined head quality attributes of head shape and branching angle. Yield and quality were not influenced by photoperiod. A better understanding of genotype and environmental interactions will help farmers optimise yield and quality, by matching cultivars with time of sowing. The estimated base and optimum temperature for broccoli development were 0°C and 20 °C, respectively, and were consistent across cultivars, but thermal time requirements for phenological intervals were cultivar specific. Differences in thermal time requirement from floral initiation to harvest maturity between cultivars were small and of little importance, but differences in thermal time requirement from emergence to floral initiation were large. Sensitivity to photoperiod and solar radiation was low in the three cultivars used. This research has produced models to assist broccoli farmers in crop scheduling and cultivar selection in south-east Queensland.

Micromanipulation of chloroplasts using optical tweezers

This paper describes experiments using optical tweezers to probe chloroplast arrangement, shape and consistency in cells of living leaf tissue and in suspension. Dual optical tweezers provided two-point contact on a single chloroplast or two-point contact on two adhered chloroplasts for manipulation in suspension. Alternatively, a microstirrer consisting of a birefringent particle trapped in an elliptically polarized laser trap was used to induce motion and tumbling of a selected chloroplast suspended in a solution. We demonstrate that displacement of chloroplasts inside the cell is extremely difficult, presumably due to chloroplast adhesion to the cytoskeleton and connections between organelles. The study also confirms that the chloroplasts are very thin and extremely cup-shaped with a concave inner surface and a convex outer surface.

Effect of multiple transverse modes in self-mixing sensors based on vertical-cavity surface-emitting lasers

We investigate the effect of coexisting transverse modes on the operation of self-mixing sensors based on vertical-cavity surface-emitting lasers (VCSELs). The effect of multiple transverse modes on the measurement of displacement and distance were examined by simulation and in laboratory experiment. The simulation model shows that the periodic change in the shape and magnitude of the self-mixing signal with modulation current can be properly explained by the different frequency-modulation coefficients of the respective transverse modes in VCSELs. The simulation results are in excellent agreement with measurements performed on single-mode and multimode VCSELs and on self-mixing sensors based on these VCSELs.

XSophe-Sophe-XeprView (R). A computer simulation software suite (v. 1.1.3) for the analysis of continuous wave EPR spectra

The XSophe-Sophe-XeprView((R)) computer simulation software suite enables scientists to easily determine spin Hamiltonian parameters from isotropic, randomly oriented and single crystal continuous wave electron paramagnetic resonance (CW EPR) spectra from radicals and isolated paramagnetic metal ion centers or clusters found in metalloproteins, chemical systems and materials science. XSophe provides an X-windows graphical user interface to the Sophe programme and allows: creation of multiple input files, local and remote execution of Sophe, the display of sophelog (output from Sophe) and input parameters/files. Sophe is a sophisticated computer simulation software programme employing a number of innovative technologies including; the Sydney OPera HousE (SOPHE) partition and interpolation schemes, a field segmentation algorithm, the mosaic misorientation linewidth model, parallelization and spectral optimisation. In conjunction with the SOPHE partition scheme and the field segmentation algorithm, the SOPHE interpolation scheme and the mosaic misorientation linewidth model greatly increase the speed of simulations for most spin systems. Employing brute force matrix diagonalization in the simulation of an EPR spectrum from a high spin Cr(III) complex with the spin Hamiltonian parameters g(e) = 2.00, D = 0.10 cm(-1), E/D = 0.25, A(x) = 120.0, A(y) = 120.0, A(z) = 240.0 x 10(-4) cm(-1) requires a SOPHE grid size of N = 400 (to produce a good signal to noise ratio) and takes 229.47 s. In contrast the use of either the SOPHE interpolation scheme or the mosaic misorientation linewidth model requires a SOPHE grid size of only N = 18 and takes 44.08 and 0.79 s, respectively. Results from Sophe are transferred via the Common Object Request Broker Architecture (CORBA) to XSophe and subsequently to XeprView((R)) where the simulated CW EPR spectra (1D and 2D) can be compared to the experimental spectra. Energy level diagrams, transition roadmaps and transition surfaces aid the interpretation of complicated randomly oriented CW EPR spectra and can be viewed with a web browser and an OpenInventor scene graph viewer.

Determining forest structural attributes using an inverted geometric-optical model in mixed eucalypt forests, Southeast Queensland, Australia

The Montreal Process indicators are intended to provide a common framework for assessing and reviewing progress toward sustainable forest management. The potential of a combined geometrical-optical/spectral mixture analysis model was assessed for mapping the Montreal Process age class and successional age indicators at a regional scale using Landsat Thematic data. The project location is an area of eucalyptus forest in Emu Creek State Forest, Southeast Queensland, Australia. A quantitative model relating the spectral reflectance of a forest to the illumination geometry, slope, and aspect of the terrain surface and the size, shape, and density, and canopy size. Inversion of this model necessitated the use of spectral mixture analysis to recover subpixel information on the fractional extent of ground scene elements (such as sunlit canopy, shaded canopy, sunlit background, and shaded background). Results obtained fron a sensitivity analysis allowed improved allocation of resources to maximize the predictive accuracy of the model. It was found that modeled estimates of crown cover projection, canopy size, and tree densities had significant agreement with field and air photo-interpreted estimates. However, the accuracy of the successional stage classification was limited. The results obtained highlight the potential for future integration of high and moderate spatial resolution-imaging sensors for monitoring forest structure and condition. (C) Elsevier Science Inc., 2000.

Quasistationarity of continuous-time Markov chains with positive drift

We shall study continuous-time Markov chains on the nonnegative integers which are both irreducible and transient, and which exhibit discernible stationarity before drift to infinity sets in. We will show how this 'quasi' stationary behaviour can be modelled using a limiting conditional distribution: specifically, the limiting state probabilities conditional on not having left 0 for the last time. By way of a dual chain, obtained by killing the original process on last exit from 0, we invoke the theory of quasistationarity for absorbing Markov chains. We prove that the conditioned state probabilities of the original chain are equal to the state probabilities of its dual conditioned on non-absorption, thus allowing us to establish the simultaneous existence and then equivalence, of their limiting conditional distributions. Although a limiting conditional distribution for the dual chain is always a quasistationary distribution in the usual sense, a similar statement is not possible for the original chain.

Central projections of Drosophila sensory neurons in the transition from embryo to larva

Sensory axons of different sensory modalities project into typical domains within insect ganglia. Tactile and gustatory axons project into a ventral layer of neuropil and proprioceptive afferents, including chordotonal axone, into an intermediate or dorsal layer. Here, we describe the central projections of sensory neurons in the first instar Drosophila larva, relating them to the projection of the same sensory afferents in the embryo and to sensory afferents of similar type in other insects. Several neurons show marked morphologic changes in their axon terminals in the transition between the embryo and larva. During a short morphogenetic period late in embryogenesis, the axon terminals of the dorsal bipolar dendrite stretch receptor change their shape and their distribution within the neuromere. In the larva, external sense organ neurons (es) project their axons into a ventral layer of neuropil. Chordotonal sensory neurons (ch) project into a slightly more dorsal region that is comparable to their projection in adults. The multiple dendrite (md) neurons show two distinctive classes of projection. One group of md neurons projects into the ventral-most neuropil region, the same region into which es neurons project. Members of this group are related by lineage to es neurons or share a requirement for expression of the same proneural gene during development. Other md neurons project into a more dorsal region. Sensory receptors projecting into dorsal neuropil possibly provide proprioceptive feedback from the periphery to central motorneurons and are candidates for future genetic and cellular analysis of simple neural circuitry. J. Comp. Neurol. 425:34-44, 2000. (C) 2000 Wiley-Liss, Inc.

Improved VDC casting of aluminium alloys through an understanding of the surface properties of the molten metal

Vertical direct chill (VDC) casting of aluminium alloys is a mature process that has evolved over many decades through gradual change to both equipment design and casting practice. Today, air-pressurised, continuous lubrication, hot top mould systems with advanced station automation are selected as the process of choice for producing extrusion billet. Specific sets of operating parameters are employed on these stations for each alloy and size combination to produce optimal billet quality. The designs and parameters are largely derived from past experience and accumulated know-how. Recent experimental work at the University of Queensland has concentrated on understanding the way in which the surface properties of liquid aluminium alloys, e.g., surface tension, wetting angle and oxide skin strength, influence the size and shape of the naturally-stab le meniscus for a given alloy, temperature and atmosphere. The wide range of alloy-and condition-dependent values measured has led to the consideration of how these properties impact the stability of the enforced molten metal meniscus within the hot top mould cavity. The actual shape and position of the enforced meniscus is controlled by parameters such as the upstream conduction distance (UCD) from sub-mould cooling and the molten metal head. The degree of deviation of this actual meniscus from the predicted stable meniscus is considered to be a key driver in surface defect formation. This paper reports on liquid alloy property results and proposes how this knowledge might be used to better design VDC mould systems and casting practices.

A new wavelet-based method for the solution of the population balance equation

A new wavelet-based method for solving population balance equations with simultaneous nucleation, growth and agglomeration is proposed, which uses wavelets to express the functions. The technique is very general, powerful and overcomes the crucial problems of numerical diffusion and stability that often characterize previous techniques in this area. It is also applicable to an arbitrary grid to control resolution and computational efficiency. The proposed technique has been tested for pure agglomeration, simultaneous nucleation and growth, and simultaneous growth and agglomeration. In all cases, the predicted and analytical particle size distributions are in excellent agreement. The presence of moving sharp fronts can be addressed without the prior investigation of the characteristics of the processes. (C) 2001 Published by Elsevier Science Ltd.