Conditional simulation of random fields by successive residuals

This paper presents a new approach to the LU decomposition method for the simulation of stationary and ergodic random fields. The approach overcomes the size limitations of LU and is suitable for any size simulation. The proposed approach can facilitate fast updating of generated realizations with new data, when appropriate, without repeating the full simulation process. Based on a novel column partitioning of the L matrix, expressed in terms of successive conditional covariance matrices, the approach presented here demonstrates that LU simulation is equivalent to the successive solution of kriging residual estimates plus random terms. Consequently, it can be used for the LU decomposition of matrices of any size. The simulation approach is termed conditional simulation by successive residuals as at each step, a small set (group) of random variables is simulated with a LU decomposition of a matrix of updated conditional covariance of residuals. The simulated group is then used to estimate residuals without the need to solve large systems of equations.

Development and characterization of polymorphic microsatellite markers in taro (Colocasia esculenta)

Microsatellite-containing sequences were isolated from enriched genomic libraries of taro (Colocasia esculenta (L.) Schott). The sequencing of 269 clones yielded 77 inserts containing repeat motifs. The majority of these (81.7%) were dinucleotide or trinucleotide repeats. The GT/CA repeat motif was the most common, accounting for 42% of all repeat types. From a total of 43 primer pairs designed, 41 produced markers within the expected size range. Sixteen (39%) were polymorphic when screened against a restricted set of taro genotypes from Southeast Asia and Oceania, with an average of 3.2 alleles detected on each locus. These markers represent a useful resource for taro germplasm management, genome mapping, and marker-assisted selection.

Solution Structures of the cis- and trans-Pro30 Isomers of a Novel 38-Residue Toxin from the Venom of Hadronyche Infensa sp. that Contains a Cystine-Knot Motif within Its Four Disulfide Bonds

The primary sequence and three-dimensional structure of a novel peptide toxin isolated from the Australian funnel-web spider Hadronyche infensa sp. is reported. ACTX-HI:OB4219 contains 38 amino acids, including eight-cysteine residues that form four disulfide bonds. The connectivities of these disulfide bonds were previously unknown but have been unambiguously determined in this study. Three of these disulfide bonds are arranged in an inhibitor cystine-knot (ICK) motif, which is observed in a range of other disulfide-rich peptide toxins. The motif incorporates an embedded ring in the structure formed by two of the disulfides and their connecting backbone segments penetrated by a third disulfide bond. Using NMR spectroscopy, we determined that despite the isolation of a single native homologous product by RP-HPLC, ACTX-HI:OB4219 possesses two equally populated conformers in solution. These two conformers were determined to arise from cis/trans isomerization of the bond preceding Pro30. Full assignment of the NMR spectra for both conformers allowed for the calculation of their structures, revealing, the presence of a triple-stranded antiparallel sheet consistent with the inhibitor cystine-knot (ICK) motif.

The Three-dimensional Solution Structure of NaD1, a New Floral Defensin from Nicotiana alata and its Application to a Homology Model of the Crop Defense Protein alfAFP

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of NaD1, a novel antifungal and insecticidal protein isolated from the flowers of Nicotiana alata. NaD1 is a basic, cysteine-rich protein of 47 residues and is the first example of a plant defensin from flowers to be characterized structurally. Its three-dimensional structure consists of an a-helix and a triple-stranded anti-parallel beta-sheet that are stabilized by four intramolecular disulfide bonds. NaD1 features all the characteristics of the cysteine-stabilized up motif that has been described for a variety of proteins of differing functions ranging from antibacterial insect defensins and ion channel-perturbing scorpion toxins to an elicitor of the sweet taste response. The protein is biologically active against insect pests, which makes it a potential candidate for use in crop protection. NaD1 shares 31% sequence identity with alfAFP, an antifungal protein from alfalfa that confers resistance to a fungal pathogen in transgenic potatoes. The structure of NaD1 was used to obtain a homology model of alfAFP, since NaD1 has the highest level of sequence identity with alfAFP of any structurally characterized antifungal defensin. The structures of NaD1 and alfAFP were used in conjunction with structure - activity data for the radish defensin Rs-AFP2 to provide an insight into structure-function relationships. In particular, a putative effector site was identified in the structure of NaD1 and in the corresponding homology model of alfAFP. (C) 2002 Elsevier Science Ltd. All rights reserved.

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.

Outcome variables for osteoarthritis clinical trials: the OMERACT-OARSI set of responder criteria

Improvement in analysis and reporting results of osteoarthritis (OA) clinical trials has been recently obtained because of harmonization and standardization of the selection of outcome variables (OMERACT 3 and OARSI). Moreover, OARSI has recently proposed the OARSI responder criteria. This composite index permits presentation of results of symptom modifying clinical trials in OA based on individual patient responses (responder yes/no). The 2 organizations (OMERACT and OARSI) established. a task force aimed at evaluating: (1) the variability of observed placebo and active treatment effects using the OARSI responder criteria; and (2) the possibility of proposing a simplified set of criteria. The conclusions of the task force were presented and discussed during the OMERACT 6 conference, where a simplified set of responder criteria (OMERACT-OARSI set of criteria) was proposed.

Skolem-type difference sets for cycle systems

Cyclic m-cycle systems of order v are constructed for all m greater than or equal to 3, and all v = 1(mod 2m). This result has been settled previously by several authors. In this paper, we provide a different solution, as a consequence of a more general result, which handles all cases using similar methods and which also allows us to prove necessary and sufficient conditions for the existence of a cyclic m-cycle system of K-v - F for all m greater than or equal to 3, and all v = 2(mod 2m).

Quasiequilibrium problem in H-spaces

In the present paper, we study the quasiequilibrium problem and generalized quasiequilibrium problem of generalized quasi-variational inequality in H-spaces by a new method. Some new equilibrium existence theorems are given. Our results are different from corresponding given results or contain some recent results as their special cases. (C) 2003 Elsevier Science Ltd. All rights reserved.

The modeling of turbulent reactive flows based on multiple mapping conditioning

A new modeling approach-multiple mapping conditioning (MMC)-is introduced to treat mixing and reaction in turbulent flows. The model combines the advantages of the probability density function and the conditional moment closure methods and is based on a certain generalization of the mapping closure concept. An equivalent stochastic formulation of the MMC model is given. The validity of the closuring hypothesis of the model is demonstrated by a comparison with direct numerical simulation results for the three-stream mixing problem. (C) 2003 American Institute of Physics.

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.

Prospects for whole genome linkage disequilibrium mapping in domestic dog breeds

Linkage disequilibrium (LD) mapping is commonly used as a fine mapping tool in human genome mapping and has been used with some success for initial disease gene isolation in certain isolated inbred human populations. An understanding of the population history of domestic dog breeds suggests that LID mapping could be routinely utilized in this species for initial genome-wide scans. Such an approach offers significant advantages over traditional linkage analysis. Here, we demonstrate, using canine copper toxicosis in the Bedlington terrier as the model, that LID mapping could be reasonably expected to be a useful strategy in low-resolution, genome-wide scans in pure-bred dogs. Significant LID was demonstrated over distances up to 33.3 cM. It is very unlikely, for a number of reasons discussed, that this result could be extrapolated to the rest of the genome. It is, however, consistent with the expectation given the population structure of canine breeds and, in this breed at least, with the hypothesis that it may be possible to utilize LID in a genome-wide scan. In this study, LD mapping confirmed the location of the copper toxicosis in Bedlington terrier gene (CT-BT) and was able to do so in a population that was refractory to traditional linkage analysis.

Solution structure of CnErg1 (Ergtoxin), a HERG specific scorpion toxin

The three-dimensional structure of chemically synthesized CnErg1 (Ergtoxin), which specifically blocks HERG (human ether-a-go-go-related gene) K+ channels, was determined by nuclear magnetic resonance spectroscopy. CnErg1 consists of a triple-stranded beta-sheet and an a-helix, as is typical of K+ channel scorpion toxins. The peptide structure differs from the canonical structures in that the first beta-strand is shorter and is nearer to the second beta-strand rather than to the third beta-strand on the C-terminus. There is also a large hydrophobic patch on the surface of the toxin, surrounding a central lysine residue, Lys13. We postulate that this hydrophobic patch is likely to form part of the binding surface of the toxin. (C) 2003 Published by Elsevier Science B.V. on behalf of the Federation of European Biochemical Societies.