Spatial structure and habitat variation in a grasshopper hybrid zone

A hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) in northern Spain has been analyzed for variation in morphology and ecology. These species are readily distinguished by the number of stridulatory pegs on the hind femur. Both sexes are fully winged and inhabit disturbed habitats throughout the study area. We develop a maximum-likelihood approach to fitting a two-dimensional cline to geographical variation in quantitative traits and for estimating associations of population mean with local habitat. This method reveals a cline in peg number approximately 30 km south of the Picos de Europa Mountains that shows substantial deviations in population mean compared with the expectations of simple tension zone models. The inclusion of variation in local vegetation in the model explains a significant proportion of the residual variation in peg number, indicating that habitat-genotype associations contribute to the observed spatial pattern. However, this association is weak, and a number of populations continue to show strong deviations in mean even after habitat is included in the final model. These outliers may be the result of long-distance colonization of sites distant from the cline center or may be due to a patchy pattern of initial contact during postglacial expansion. As well as contrasting with the smooth hybrid zones described for Chorthippus parallelus, this situation also contrasts with the mosaic hybrid zones observed in Gryllus crickets and in parts of the hybrid zone between Bombina toad species, where habitat-genotype associations account for substantial amounts of among-site variation.

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.

EphA receptors and ephrin-A ligands exhibit highly regulated spatial and temporal expression patterns in the developing olfactory system

The spatiotemporal expression patterns of the chemorepulsive EphA receptors, EphA4 and EphA7, and three ephrins-A2, A4 and A5, were examined in the developing rat primary olfactory system. Unlike the visual system that has simple and stable gradients of Ephs and ephrins, the olfactory system demonstrates complex spatiotemporal expression patterns of these molecules. Using immunohistochemistry, we demonstrate that expression of these molecules is dynamic and tightly regulated both within and between different cell types. We reveal restricted targeting of these proteins within subcellular compartments of some neurons. EphA4, ephrin-A2 and ephrin-A5 were expressed by primary olfactory axons during the embryonic formation of the olfactory nerve. There were no gradients in expression along the rostrocaudal or ventrodorsal axes in the nasal cavity and olfactory bulb. However, during the early neonatal period, axons expressing different levels of ephrin-A5 sorted out and terminated in a subpopulation of glomeruli that were mosaically dispersed throughout the bulb. The expression of EphA4 and ephrin-A2 was dramatically down-regulated on all axons during the early neonatal period of glomerular formation. The uniform co-expression of receptors and ligands before glomerular formation suggests they play a generic role in axon-axon interactions in the olfactory nerve and nerve fibre layer. In contrast, loss of EphA4 from axons during glomerular formation may facilitate the interaction of ephrin-A5 with Eph receptors on target cells in the bulb. While EphA4, EphA5 and EphA7 are not mosaically expressed by bulbar neurons, other Eph receptors may have expression patterns complementary to the ephrin-A5-positive subpopulation of glomeruli. (C) 2002 Elsevier Science B.V. All rights reserved.

Spatial and temporal patterns in Australian wheat yield and their relationship with ENSO

Spatial and temporal variability in wheat production in Australia is dominated by rainfall occurrence. The length of historical production records is inadequate, however, to analyse spatial and temporal patterns conclusively. In this study we used modelling and simulation to identify key spatial patterns in Australian wheat yield, identify groups of years in the historical record in which spatial patterns were similar, and examine association of those wheat yield year groups with indicators of the El Nino Southern Oscillation (ENSO). A simple stress index model was trained on 19 years of Australian Bureau of Statistics shire yield data (1975-93). The model was then used to simulate shire yield from 1901 to 1999 for all wheat-producing shires. Principal components analysis was used to determine the dominating spatial relationships in wheat yield among shires. Six major components of spatial variability were found. Five of these represented near spatially independent zones across the Australian wheatbelt that demonstrated coherent temporal (annual) variability in wheat yield. A second orthogonal component was required to explain the temporal variation in New South Wales. The principal component scores were used to identify high- and low-yielding years in each zone. Year type groupings identified in this way were tested for association with indicators of ENSO. Significant associations were found for all zones in the Australian wheatbelt. Associations were as strong or stronger when ENSO indicators preceding the wheat season (April-May phases of the Southern Oscillation Index) were used rather than indicators based on classification during the wheat season. Although this association suggests an obvious role for seasonal climate forecasting in national wheat crop forecasting, the discriminatory power of the ENSO indicators, although significant, was not strong. By examining the historical years forming the wheat yield analog sets within each zone, it may be possible to identify novel climate system or ocean-atmosphere features that may be causal and, hence, most useful in improving seasonal forecasting schemes.

Spatial distribution of vectors of Ross River virus and Barmah Forest virus on Russell Island, Moreton Bay, Queensland

We used a network of 20 carbon dioxide- and octenol-supplemented light traps to sample adult mosquitoes throughout Russell Island in southern Moreton Bay, south-east Queensland. Between February and April 2001, an estimated 1365 564 adult female mosquitoes were collected. In contrast to an average catch of 9754 female mosquitoes per trap night on Russell Island, reference traps set on Macleay Island and on the mainland returned average catches of 3172 and 222, respectively. On Russell Island, Ochlerotatus vigilax (Skuse), Coquillettidia linealis (Skuse), Culex annulirostris Skuse and Verrallina funerea (Theobald), known or suspected vectors of Ross River (RR) and/or Barmah Forest (BF) viruses, comprised 89.6% of the 25 taxa collected. When the spatial distributions of the above species were mapped and analysed using local spatial statistics, all were found to be present in highest numbers towards the southern end of the island during most of the 7 weeks. This indicated the presence of more suitable adult harbourage sites and/or suboptimal larval control efficacy. As immature stages and the breeding habitat of Cq. linealis are as yet undescribed, this species in particular presents a considerable impediment to proposed development scenarios. The method presented here of mapping the numbers of mosquitoes throughout a local government area allows specific areas that have high vector numbers to be defined.

Temporal and spatial expression of fibroblast growth factor receptor 4 isoforms in murine tissues

Fibroblast growth factor receptors (FGFRs) undergo highly regulated spatial and temporal changes of expression during development. This study describes the use of quantitative reverse transcriptase-polymerase chain reaction and immunochemistry to assess the changes in expression of FGFR4 as compared to its FGFR4-17a and -17b isoforms in mouse tissues, from early embryogenesis through to adulthood. Compared to FGFR4, the expression of the isoforms is more restricted at all developmental stages tested. The reverse transcriptase-polymerase chain reaction demonstrated that FGFR4 is expressed in more tissue types than either of its isoforms: it was found predominantly in lung, liver, brain, skeletal muscle and kidney, whereas the FGFR4-17a form was detected in lung and skeletal muscle, and the FGFR4-17b form only in lung, liver, skeletal muscle and kidney. Immunohistochemistry confirmed strong FGFR4-17b expression in the postnatal lung. When combined, the results suggest that FGFR4 variants play important roles particularly in lung and skeletal muscle development.

Spatial and temporal changes in multiple hormone groups during lateral bud release shortly following apex decapitation of chickpea (Cicer arietinum) seedlings

Although the co-ordination of promotive root-sourced cytokinin (CK) and inhibitory shoot apex-sourced auxin (IAA) is central to all current models on lateral bud dormancy release, control by those hormones alone has appeared inadequate in many studies. Thus it was hypothesized that the IAA : CK model is the central control but that it must be considered within the relevant timeframe leading to lateral bud release and against a backdrop of interactions with other hormone groups. Therefore, IAA and a wide survey of cytokinins (CKs), were examined along with abscisic acid (ABA) and polyamines (PAs) in released buds, tissue surrounding buds and xylem sap at 1 and 4 h after apex removal, when lateral buds of chickpea are known to break dormancy. Three potential lateral bud growth inhibitors, IAA, ABA and cis-zeatin 9-riboside (ZR), declined sharply in the released buds and xylem following decapitation. This is in contrast to potential dormancy breaking CKs like trans-ZR and trans-zeantin 9-riboside 5'phosphate (ZRMP), which represented the strongest correlative changes by increasing 3.5-fold in xylem sap and 22-fold in buds. PAs had not changed significantly in buds or other tissues after 4 h, so they were not directly involved in the breaking of bud dormancy. Results from the xylem and surrounding tissues indicated that bud CK increases resulted from a combination synthesis in the bud and selective loading of CK nucleotides into the xylem from the root.

Direct visualization of Ras proteins in spatially distinct cell surface microdomains

Localization of signaling complexes to specific micro-domains coordinates signal transduction at the plasma membrane. Using immunogold electron microscopy of plasma membrane sheets coupled with spatial point pattern analysis, we have visualized morphologically featureless microdomains including lipid rafts, in situ and at high resolution. We find that an inner-plasma membrane lipid raft marker displays cholesterol-dependent clustering in microdomains with a mean diameter of 44 nm that occupy 35% of the cell surface. Cross-linking an outer-leaflet raft protein results in the redistribution of inner leaflet rafts, but they retain their modular structure. Analysis of Ras microlocalization shows that inactive H-ras is distributed between lipid rafts and a cholesterol-independent micro-domain. Conversely, activated H-ras and K-ras reside predominantly in nonoverlapping, cholesterol-independent microdomains. Galectin-1 stabilizes the association of activated H-ras with these nonraft microdomains, whereas K-ras clustering is supported by farnesylation, but not geranylgeranylation. These results illustrate that the inner plasma membrane comprises a complex mosaic of discrete microdomains. Differential spatial localization within this framework can likely account for the distinct signal outputs from the highly homologous Ras proteins.

Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Linking physiological and architectural models of cotton

Despite the strong influence of plant architecture on crop yield, most crop models either ignore it or deal with it in a very rudimentary way. This paper demonstrates the feasibility of linking a model that simulates the morphogenesis and resultant architecture of individual cotton plants with a crop model that simulates the effects of environmental factors on critical physiological processes and resulting yield in cotton. First the varietal parameters of the models were made concordant. Then routines were developed to allocate the flower buds produced each day by the crop model amongst the potential positions generated by the architectural model. This allocation is done according to a set of heuristic rules. The final weight of individual bolls and the shedding of buds and fruit caused by water, N, and C stresses are processed in a similar manner. Observations of the positions of harvestable fruits, both within and between plants, made under a variety of agronomic conditions that had resulted in a broad range of plant architectures were compared to those predicted by the model with the same environmental inputs. As illustrated by comparisons of plant maps, the linked models performed reasonably well, though performance of the fruiting point allocation and shedding algorithms could probably be improved by further analysis of the spatial relationships of retained fruit. (C) 2002 Elsevier Science Ltd. All rights reserved.

Knowledge-based system on optimum design of liquid retaining structures with genetic algorithms

This paper delineates the development of a prototype hybrid knowledge-based system for the optimum design of liquid retaining structures by coupling the blackboard architecture, an expert system shell VISUAL RULE STUDIO and genetic algorithm (GA). Through custom-built interactive graphical user interfaces under a user-friendly environment, the user is directed throughout the design process, which includes preliminary design, load specification, model generation, finite element analysis, code compliance checking, and member sizing optimization. For structural optimization, GA is applied to the minimum cost design of structural systems with discrete reinforced concrete sections. The design of a typical example of the liquid retaining structure is illustrated. The results demonstrate extraordinarily converging speed as near-optimal solutions are acquired after merely exploration of a small portion of the search space. This system can act as a consultant to assist novice designers in the design of liquid retaining structures.