Exploring complex fitness surfaces: Multiple ornamentation and polymorphism in male guppies

The sexual ornamentation used by male guppies to attract females comprises many components, each of which varies considerably among males. Although natural and sexual selection have been shown to contribute to divergence among populations in male sexual ornaments, the role of sexual selection in maintaining polymorphism within populations is less clear. We used both parametric quadratic regression and nonparametric projection pursuit regression techniques to reveal the major axes of non-linear sexual selection on male ornaments. We visualized the fitness surfaces defined by these axes using thin-plate splines to allow a direct comparison of the two methodologies. Identification of the major axes of selection and their visualization was critical in determining the form and strength of nonlinear selection. Both types of analysis revealed fitness surfaces comprising three peaks, suggesting that there is more than one way to make an attractive guppy. Disruptive selection may be an important process underlying the presence of multiple sexual ornaments and may contribute to the maintenance of the high levels of polymorphism in male sexual ornaments found in guppy populations.

Further evidence for the reliance of catalysis by rabbit muscle pyruvate kinase upon isomerization of the ternary complex between enzyme and products

Isothermal calorimetry has been used to examine the effect of thermodynamic non-ideality on the kinetics of catalysis by rabbit muscle pyruvate kinase as the result of molecular crowding by inert cosolutes. The investigation, designed to detect substrate-mediated isomerization of pyruvate kinase, has revealed a 15% enhancement of maximal velocity by supplementation of reaction mixtures with 0.1 M proline, glycine or sorbitol. This effect of thermodynamic non-ideality implicates the existence of a substrate-induced conformational change that is governed by a minor volume decrease and a very small isomerization constant; and hence, substantiates earlier inferences that the rate-determining step in pyruvate kinase kinetics is isomerization of the ternary enzyme product complex rather than the release of products. (C) 2003 Elsevier Science B.V. All rights reserved.

A partial velocity approach to subcycling structural dynamics

Subcycling, or the use of different timesteps at different nodes, can be an effective way of improving the computational efficiency of explicit transient dynamic structural solutions. The method that has been most widely adopted uses a nodal partition. extending the central difference method, in which small timestep updates are performed interpolating on the displacement at neighbouring large timestep nodes. This approach leads to narrow bands of unstable timesteps or statistical stability. It also can be in error due to lack of momentum conservation on the timestep interface. The author has previously proposed energy conserving algorithms that avoid the first problem of statistical stability. However, these sacrifice accuracy to achieve stability. An approach to conserve momentum on an element interface by adding partial velocities is considered here. Applied to extend the central difference method. this approach is simple. and has accuracy advantages. The method can be programmed by summing impulses of internal forces, evaluated using local element timesteps, in order to predict a velocity change at a node. However, it is still only statistically stable, so an adaptive timestep size is needed to monitor accuracy and to be adjusted if necessary. By replacing the central difference method with the explicit generalized alpha method. it is possible to gain stability by dissipating the high frequency response that leads to stability problems. However. coding the algorithm is less elegant, as the response depends on previous partial accelerations. Extension to implicit integration, is shown to be impractical due to the neglect of remote effects of internal forces acting across a timestep interface. (C) 2002 Elsevier Science B.V. All rights reserved.

Dynamics of Gray Matter Loss in Alzheimer's Disease

We detected and mapped a dynamically spreading wave of gray matter loss in the brains of patients with Alzheimer's disease (AD). The loss pattern was visualized in four dimensions as it spread over time from temporal and limbic cortices into frontal and occipital brain regions, sparing sensorimotor cortices. The shifting deficits were asymmetric (left hemisphere >right hemisphere) and correlated with progressively declining cognitive status ( p 15% loss). The maps distinguished different phases of AD and differentiated AD from normal aging. Local gray matter loss rates (5.3 +/- 2.3% per year in AD v 0.9 +/- 0.9% per year in controls) were faster in the left hemisphere ( p < 0.029) than the right. Transient barriers to disease progression appeared at limbic/frontal boundaries. This degenerative sequence, observed in vivo as it developed, provides the first quantitative, dynamic visualization of cortical atrophic rates in normal elderly populations and in those with dementia.

Speciation and phylogeography in Caridina indistincta, a complex of freshwater shrimps from Australian heathland streams

The Caridina indistincta complex is a group of closely related atyid shrimps that inhabit coastal freshwater streams throughout north-eastern Australia. Using mitochondrial DNA sequence data (cytochrome oxidase 1, CO1), we (1) inferred the timing of speciation in the C. indistincta group and (2) examined the intraspecific phylogeographic patterns within the group. Assuming a shrimp-specific rate of CO1 evolution, the level of sequence divergence among species suggests that speciation took place during the Miocene epoch. Within one widespread mainland species, phylogeographic patterns suggest strong geographic 'regionalisation' of mtDNA lineages that are most likely of Pleistocene origin. By contrast, another species comprises two highly divergent mtDNA lineages that occur in sympatry. We suggest that although Pleistocene sea-level regressions appear important in generating population-level phylogeographic patterns, these events were largely unimportant in the formation of species in this group.

Epitope map for a growth hormone receptor agonist monoclonal antibody, MAb 263

Monoclonal antibody (MAb) 263 is a widely used monoclonal antibody that recognizes the extracellular domain (ECD) of the GH receptor. It has been shown to act as a GH agonist both in vitro and in vivo, and we report here that it must be divalent to exert its effect on the full-length receptor. To understand the mechanism of its agonist action, we have determined the precise epitope for this antibody using a novel random PCR mutagenesis approach together with expression screening in yeast. A library of 5200 clones of rabbit GH receptor ECD mutants were screened both with MAb 263 and with an anticarboxy-tag antibody to verify complete ECD expression. Sequencing for clones that expressed complete ECD but were not MAb 263 positive identified 20 epitope residues distributed in a discontinuous manner throughout the ECD. The major part of the epitope, as revealed after mapping onto the crystal structure model of the ECD molecule, was located on the side and upper portion of domain 1, particularly within the D - E strand disulfide loop 79 - 96. Molecular dynamics docking of an antibody of the same isotype as MAb 263 was used to dock the bivalent antibody to the 1528-Angstrom(2) epitope and to visualize the likely consequences of MAb binding. The minimized model enables the antibody to grasp two receptors in a pincer-like movement from opposite sides, facilitating alignment of the receptor dimerization domains in a manner similar to, but not identical with, GH.

The Lake Tekapo experiment (LTEX): An investigation of atmospheric boundary layer processes in complex terrain

A research program on atmospheric boundary layer processes and local wind regimes in complex terrain was conducted in the vicinity of Lake Tekapo in the southern Alps of New Zealand, during two 1-month field campaigns in 1997 and 1999. The effects of the interaction of thermal and dynamic forcing were of specific interest, with a particular focus on the interaction of thermal forcing of differing scales. The rationale and objectives of the field and modeling program are described, along with the methodology used to achieve them. Specific research aims include improved knowledge of the role of surface forcing associated with varying energy balances across heterogeneous terrain, thermal influences on boundary layer and local wind development, and dynamic influences of the terrain through channeling effects. Data were collected using a network of surface meteorological and energy balance stations, radiosonde and pilot balloon soundings, tethered balloon and kite-based systems, sodar, and an instrumented light aircraft. These data are being used to investigate the energetics of surface heat fluxes, the effects of localized heating/cooling and advective processes on atmospheric boundary layer development, and dynamic channeling. A complementary program of numerical modeling includes application of the Regional Atmospheric Modeling System (RAMS) to case studies characterizing typical boundary layer structures and airflow patterns observed around Lake Tekapo. Some initial results derived from the special observation periods are used to illustrate progress made to date. In spite of the difficulties involved in obtaining good data and undertaking modeling experiments in such complex terrain, initial results show that surface thermal heterogeneity has a significant influence on local atmospheric structure and wind fields in the vicinity of the lake. This influence occurs particularly in the morning. However, dynamic channeling effects and the larger-scale thermal effect of the mountain region frequently override these more local features later in the day.

Influence of hysteresis on tidal capillary fringe dynamics in a well-sorted sand

Nielsen and Perrochet [Adv. Water Resour. 23 (2000) 503] presented experimental data for cyclic water movement in the vadose zone above an oscillating watertable. The response of the watertable to cyclic forcing was characterised by the ratios of the forcing head to watertable amplitudes and their associated phase lag. They found that their non-hysteretic Richards' equation model failed to represent the observed behaviour of these parameters. This paper explores the effect on the simulated capillary fringe dynamics (in terms of these parameters) of including varying degrees of hysteresis in the moisture retention curve used in a numerical model of their experiment. It is clear that hysteresis can indeed account for observed discrepancies between simulation and experiment and that the effect of hysteresis varies with the frequency of oscillation. The use of a single-valued mean retention curve, as advocated by some authors, fails to provide a match between the simulated and observed behaviour of the Nielsen and Perrochet parameters, but is shown to be adequate for predicting time-averaged soil moisture profiles. (C) 2003 Elsevier Ltd. All rights reserved.

Genomovar diversity amongst Burkholderia cepacia complex isolates from an Australian adult cystic fibrosis unit

In this study, a combination of recA-based PCR assays and 16S rDNA restriction fragment length polymorphism (RFLP) analysis was used to determine the genomovar diversity of clinical Burkholderia cepacia complex isolates. Twenty-eight isolates were prospectively collected from patients attending a large Australian adult cystic fibrosis (CF) unit, 22 isolates were referred from other Australian CF units and a further eight isolates originated from patients without CF. The 28 prospectively collected isolates were distributed amongst the following genomovars: Burkholderia cepacia genomovar I (28.6%), Burkholderia multivorans (21.4%), Burkholderia cepacia genomovar III (39.3%), Burkholderia vietnamiensis (3.6%) and Burkholderia ambifaria (7.1%). The results of this study highlight the usefulness of 16S rDNA RFLP typing for the identification of other Burkholderia spp. and non-fermenting gram-negative bacteria.

Genotype and environment effects on dynamics of harvest index during grain filling in sorghum

An approach based on a linear rate of increase in harvest index (141) with time after anthesis has been used as a simple means-to predict grain growth and yield in many crop simulation models. When applied to diverse situations, however, this approach has been found to introduce significant error in grain yield predictions. Accordingly, this study was undertaken to examine the stability of the HI approach for yield prediction in sorghum [Sorghum bicolor (L.) Moench]. Four field experiments were conducted under nonlimiting water. and N conditions. The experiments were sown at times that ensured a broad range in temperature and radiation conditions. Treatments consisted of two population densities and three genotypes varying in maturity. Frequent sequential harvests were used to monitor crop growth, yield, and the dynamics of 111. Experiments varied greatly in yield and final HI. There was also a tendency for lower HI with later maturity. Harvest index dynamics also varied among experiments and, to a lesser extent, among treatments within experiments. The variation was associated mostly with the linear rate of increase in HI and timing of cessation of that increase. The average rate of HI increase was 0.0198 d(-1), but this was reduced considerably (0.0147) in one experiment that matured in cool conditions. The variations found in IN dynamics could be largely explained by differences in assimilation during grain filling and remobilization of preanthesis assimilate. We concluded that this level of variation in HI dynamics limited the general applicability of the HI approach in yield prediction and suggested a potential alternative for testing.

Phylogeny of the Sporobolus indicus complex, based on internal transcribed spacer (ITS) sequences

The entire internal transcribed spacer ( ITS) region, including the 5.8S subunit of the nuclear ribosomal DNA ( rDNA), was sequenced by direct double-stranded sequencing of polymerase chain reaction (PCR) amplified fragments. The study included 40 Sporobolus ( Family Poaceae, subfamily Chloridoideae) seed collections from 14 putative species ( all 11 species from the S. indicus complex and three Australian native species). These sequences, along with those from two out-group species [ Pennisetum alopecuroides ( L.) Spreng. and Heteropogon contortus ( L.) P. Beauv. ex Roemer & Schultes, Poaceae, subfamily Panicoideae], were analysed by the parsimony method (PAUP; version 4.0b4a) to infer phylogenetic relationships among these species. The length of the ITS1, 5.8S subunit and ITS2 region were 222, 164 and 218 base pairs ( bp), respectively, in all species of the S. indicus complex, except for the ITS2 region of S. diandrus P. Beauv. individuals, which was 217 bp long. Of the 624 characters included in the analysis, 245 ( 39.3%) of the 330 variable sites contained potential phylogenetic information. Differences in sequences among the members of the S. pyramidalis P. Beauv., S. natalensis (Steud.) Dur & Schinz and S. jacquemontii Kunth. collections were 0%, while differences ranged from 0 to 2% between these and other species of the complex. Similarly, differences in sequences among collections of S. laxus B. K. Simon, S. sessilis B. K. Simon, S. elongatus R. Br. and S. creber De Nardi were 0%, compared with differences of 1-2% between these four species and the rest of the complex. When comparing S. fertilis ( Steud.) Clayton and S. africanus (Poir.) Robyns & Tourney, differences between collections ranged from 0 to 1%. Parsimony analysis grouped all 11 species of the S. indicus complex together, indicating a monophyletic origin. For the entire data set, pair-wise distances among members of the S. indicus complex varied from 0.00 to 1.58%, compared with a range of 20.08-21.44% among species in the complex and the Australian native species studied. A strict consensus phylogenetic tree separated 11 species of the S. indicus complex into five major clades. The phylogeny, based on ITS sequences, was found to be congruent with an earlier study on the taxonomic relationship of the weedy Sporobolus grasses revealed from random amplified polymorphic DNA ( RAPD). However, this cladistic analysis of the complex was not in agreement with that created on past morphological analyses and therefore gives a new insight into the phylogeny of the S. indicus complex.

SPAnDX: a process-based population dynamics model to explore management and climate change impacts on an invasive alien plant, Acacia nilotica

This paper describes a process-based metapopulation dynamics and phenology model of prickly acacia, Acacia nilotica, an invasive alien species in Australia. The model, SPAnDX, describes the interactions between riparian and upland sub-populations of A. nilotica within livestock paddocks, including the effects of extrinsic factors such as temperature, soil moisture availability and atmospheric concentrations of carbon dioxide. The model includes the effects of management events such as changing the livestock species or stocking rate, applying fire, and herbicide application. The predicted population behaviour of A. nilotica was sensitive to climate. Using 35 years daily weather datasets for five representative sites spanning the range of conditions that A. nilotica is found in Australia, the model predicted biomass levels that closely accord with expected values at each site. SPAnDX can be used as a decision-support tool in integrated weed management, and to explore the sensitivity of cultural management practices to climate change throughout the range of A. nilotica. The cohort-based DYMEX modelling package used to build and run SPAnDX provided several advantages over more traditional population modelling approaches (e.g. an appropriate specific formalism (discrete time, cohort-based, process-oriented), user-friendly graphical environment, extensible library of reusable components, and useful and flexible input/output support framework). (C) 2003 Published by Elsevier Science B.V.

Synthesis of optically complex core-shell colloidal suspensions: Pathways to multiplexed biological screening

The ability to generate enormous random libraries of DNA probes via split-and-mix synthesis on solid supports is an important biotechnological application of colloids that has not been fully utilized to date. To discriminate between colloid-based DNA probes each colloidal particle must be 'encoded' so it is distinguishable from all other particles. To this end, we have used novel particle synthesis strategies to produce large numbers of optically encoded particle suitable for DNA library synthesis. Multifluorescent particles with unique and reproducible optical signatures (i.e., fluorescence and light-scattering attributes) suitable for high-throughput flow cytometry have been produced. In the spectroscopic study presented here, we investigated the optical characteristics of multi-fluorescent particles that were synthesized by coating silica 'core' particles with up to six different fluorescent dye shells alternated with non-fluorescent silica 'spacer' shells. It was observed that the diameter of the particles increased by up to 20% as a result of the addition of twelve concentric shells and that there was a significant reduction in fluorescence emission intensities from inner shells as an increasing number of shells were deposited.

The preparation of zinc(II) and cadmium(II) complexes of the pentadentate N3S2 ligand formed from 2,6-diacetylpyridine and S-benzyldithiocarbazate (H2SNNNS) and the X-ray crystal structure of the novel dimeric [Zn2(SNNNS)2] complex

The pentadentate chelating agent, 2,6-diacetylpyridinebis(S-benzyldithiocarbazate) (H2SNNNS) reacts with zinc(II) and cadmium(II) ions forming stable complexes of empirical formula, [M(SNNNS)] (M=Zn2+, Cd2+; SNNNS2 =doubly deprotonated anionic form of the Schiff base). These complexes have been characterized by a variety of physico-chemical techniques. IR and H-1 NMR spectral evidence indicate that the Schiff base coordinates to the zinc(II) and cadmium(II) ions via the pyridine nitrogen atoms, the azomethine nitrogen atoms and the mercaptide sulfur atoms. The crystal and molecular structure of the zinc(II) complex has been determined by X-ray diffraction. The complex is a dimer in which the pyridine nitrogen atom,the azomethine nitrogen atom and the thiolate sulfur atom from one ligand coordinate to one of the zinc(II) ions whereas the azomethine and thiolate sulfur atoms from another ligand complete pentacoordination around the zinc(II) ion, the ligands being coordinated in their deprotonated forms. The coordination geometry about each zinc(II) can be considered as intermediate between a square-pyramid and trigonal-bipyramid. The cadmium(II) complex is also assigned with a dimeric structure. (C) 2003 Elsevier Ltd. All rights reserved.