MapScript: A map algebra programming language incorporating neighborhood analysis

Map algebra is a data model and simple functional notation to study the distribution and patterns of spatial phenomena. It uses a uniform representation of space as discrete grids, which are organized into layers. This paper discusses extensions to map algebra to handle neighborhood operations with a new data type called a template. Templates provide general windowing operations on grids to enable spatial models for cellular automata, mathematical morphology, and local spatial statistics. A programming language for map algebra that incorporates templates and special processing constructs is described. The programming language is called MapScript. Example program scripts are presented to perform diverse and interesting neighborhood analysis for descriptive, model-based and processed-based analysis.

A technique for modular logic program refinement

A refinement calculus provides a method for transforming specifications to executable code, maintaining the correctness of the code with respect to its specification. In this paper we introduce modules into a logic programming refinement calculus. Modules allow data types to be grouped together with sets of procedures that manipulate the data types. By placing restrictions on the way a program uses a module, we develop a technique for refining the module so that it uses a more efficient representation of the data type.

Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

Instantaneous outbursts in underground coal mines: An overview and association with coal type

Instantaneous outbursts in underground coal mines have occurred in at least 16 countries, involving both methane (CH4) and carbon dioxide (CO2). The precise mechanisms of an instantaneous outburst are still unresolved but must consider the effects of stress, gas content and physico-mechanical properties of the coal. Other factors such as mining methods (e.g., development heading into the coal seam) and geological features (e.g., coal seam disruptions from faulting) can combine to exacerbate the problem. Prediction techniques continue to be unreliable and unexpected outburst incidents resulting in fatalities are a major concern for underground coal operations. Gas content thresholds of 9 m(3)/t for CH4 and 6 m(3)/t for CO2 are used in the Sydney Basin, to indicate outburst-prone conditions, but are reviewed on an individual mine basis and in mixed as situations. Data on the sorption behaviour of Bowen Basin coals from Australia have provided an explanation for the conflicting results obtained by coal face desorption indices used for outburst-proneness assessment. A key factor appears to be different desorption rates displayed by banded coals, which is supported by both laboratory and mine-site investigations. Dull coal bands with high fusinite and semifusinite contents tend to display rapid desorption from solid coal, for a given pressure drop. The opposite is true for bright coal bands with high vitrinite contents and dull coal bands with high inertodetrinite contents. Consequently, when face samples of dull, fusinite-or semifusinite-rich coal of small particle size are taken for desorption testing, much gas has already escaped and low readings result. The converse applies for samples taken from coal bands with high vitrinite and/or inertodetrinite contents. In terms of outburst potential, it is the bright, vitrinite-rich and the dull, inertodetrinite-rich sections of a coal seam that appear to be more outburst-prone. This is due to the ability of the solid coal to retain gas, even after pressure reduction, creating a gas content gradient across the coal face sufficient to initiate an outburst. Once the particle size of the coal is reduced, rapid gas desorption can then take place. (C) 1998 Elsevier Science.

Structure-activity relationships of omega-conotoxins MVIIA, MVIIC and 14 loop splice hybrids at N and P/Q-type calcium channels

The omega-conotoxins are a set of structurally related, four-loop, six cysteine containing peptides, that have a range of selectivities for different subtypes of the voltage-sensitive calcium channel (VSCC). To investigate the basis of the selectivity displayed by these peptides, we have studied the binding affinities of two naturally occurring omega-conotoxins, MVIIA and MVIIC and a series of 14 MVIIA/MVIIC loop hybrids using radioligand binding assays for N and P/Q-type Ca2+ channels in rat brain tissue. A selectivity profile was developed from the ratio of relative potencies at N-type VSCCs (using [I-125]GVIA radioligand binding assays) and P/Q-type VSCCs (using [I-125]MVIIC radioligand binding assays). in these peptides, loops 2 and 4 make the greatest contribution to VSCC subtype selectivity, while the effects of loops 1 and 3 are negligible. Peptides with homogenous combinations of loop 2 and 4 display clear selectivity preferences, while those with heterogeneous combinations of loops 2 and 4 are less discriminatory. H-1 NMR spectroscopy revealed that the global folds of MVIIA, MVIIC and the 14 loop hybrid peptides were similar; however, several differences in local structure were identified. Based on the binding data and the 3D structures of MVIIA, GVIA and MVIIC, we have developed a preliminary pharmacophore based on the omega-conotoxin residues most Likely to interact with the N-type VSCC. (C) 1999 Academic Press.

Functional implications of the human T-lymphotropic virus type 1 transmembrane glycoprotein helical hairpin structure

Retrovirus entry into cells follows receptor binding by the surface exposed envelope glycoprotein (Env) subunit (SU), which triggers the membrane fusion activity of the transmembrane (TM) protein. TM protein fragments expressed in the absence of SU adopt helical hairpin structures comprising a central coiled coil, a region of chain reversal containing a disulfide-bonded loop, and a C-terminal segment that packs onto the exterior of the coiled coil in an antiparallel manner. Here we used in vitro mutagenesis to test the functional role of structural elements observed in a model helical hairpin, gp21 of human T-lymphotropic virus type 1. Membrane fusion activity requires the stabilization of the N and C termini of the central coiled coil by a hydrophobic N cap and a small hydrophobic core, respectively. A conserved Gly-Gly hinge motif preceding the disulfide-bonded loop, a salt bridge that stabilizes the chain reversal region, and interactions between the C-terminal segment and the coiled coil are also critical for fusion activity. Our data support a model whereby the chain reversal region transmits a conformational signal from receptor-bound SU to induce the fusion-activated helical hairpin conformation of the TM protein.

Multiple conformational epitopes are recognized by natural and induced immunity to the E7 protein of human papilloma virus type 16 in man

The reactivity of sera from patients with cervical cancer with the E7 protein of human papilloma virus type 16 (HPV16) was estimated using a novel non-radioactive immunoprecipitation assay and four established protein-and peptide-based immunoassays. Six of 14 sera from patients with cervical cancer and 1 of 10 sera from healthy laboratory staff showed repeated reactivity with E7 in at least one assay. Four of the 7 reactive sera were consistently reactive in more than one assay, but only one was reactive in all four assays. Following immunization with E7, 2 of 5 patients with cervical cancer had increased E7-specific reactivity, measurable in one or more assays. No single assay was particularly sensitive for E7 reactivity, or predictive of cervical cancer. Mapping of E7 reactivity to specific E7 peptides was unsuccessful, suggesting that natural or induced E7 reactivity in human serum is commonly directed to conformational epitopes of E7, These results suggest that each assay employed with is study measures a different aspect of E7 reactivity, and that various reactivities to E7 may manifest following HPV infection or immunization. This finding is of significance for monitoring of E7 immunotherapy and for serological screening for cervical cancer. Copyright (C) 2000 S.Karger, AG. Basel.

Conformations of platypus venom C-type natriuretic peptide in aqueous solution and sodium dodecyl sulfate micelles

Nuclear magnetic resonance spectroscopy was used to investigate the conformations of the platypus venom C-type natriuretic peptide A (OvCNPa) in aqueous solutions and in solutions containing sodium dodecyl sulfate (SDS) micelles. The chemically synthesized OvCNPa showed a substantial decrease in flexibility in aqueous solution at 10 degreesC, allowing the observation of medium- and long-range nuclear Overhauser enhancement (NOE) connectivities. Three-dimensional structures calculated using these data showed flexible and reasonably well-defined regions, the locations of which were similar in the two solvents. In aqueous solution, the linear part that spans residues 3-14 was basically an extended conformation while the cyclic portion, defined by residues 23-39, contained a series of beta-turns. The overall shape of the cyclic portion was similar to that observed for an atrial natriuretic peptide (ANP) variant in aqueous solution. OvCNPa adopted a different conformation in SDS micelles wherein the N-terminal region, defined by residues 2-10, was more compact, characterised by turns and a helix, while the cyclic region had turns and an overall shape that was fundamentally different from those structures observed in aqueous solution. The hydrophobic cluster, situated at the centre of the ring of the structure in aqueous solution, was absent in the structure in the presence of SDS micelles. Thus, OvCNPa interacts with SDS micelles and can possibly form ion-channels in cell membranes. (C) 2002 Elsevier Science Ltd. All rights reserved.

Spatial data mining for enhanced soil map modelling

The principle of using induction rules based on spatial environmental data to model a soil map has previously been demonstrated Whilst the general pattern of classes of large spatial extent and those with close association with geology were delineated small classes and the detailed spatial pattern of the map were less well rendered Here we examine several strategies to improve the quality of the soil map models generated by rule induction Terrain attributes that are better suited to landscape description at a resolution of 250 m are introduced as predictors of soil type A map sampling strategy is developed Classification error is reduced by using boosting rather than cross validation to improve the model Further the benefit of incorporating the local spatial context for each environmental variable into the rule induction is examined The best model was achieved by sampling in proportion to the spatial extent of the mapped classes boosting the decision trees and using spatial contextual information extracted from the environmental variables.

Image calibration to like-values in mapping shallow water quality from multitemporal data

The applicability of image calibration to like-values in mapping water quality parameters from multitemporal images is explored, Six sets of water samples were collected at satellite overpasses over Moreton Bay, Brisbane, Australia. Analysis of these samples reveals that waters in this shallow bay are mostly TSS-dominated, even though they are occasionally dominated by chlorophyll as well. Three of the images were calibrated to a reference image based on invariant targets. Predictive models constructed from the reference image were applied to estimating total suspended sediment (TSS) and Secchi depth from another image at a discrepancy of around 35 percent. Application of the predictive model for TSS concentration to another image acquired at a time of different water types resulted in a discrepancy of 152 percent. Therefore, image calibration to like-values could be used to reliably map certain water quality parameters from multitemporal TM images so long as the water type under study remains unchanged. This method is limited in that the mapped results could be rather inaccurate if the water type under study has changed considerably. Thus, the approach needs to be refined in shallow water from multitemporal satellite imagery.

How well does B-type natriuretic peptide predict death and cardiac events in patients with heart failure: systematic review

Objective To assess how well B-type natriuretic peptide (BNP) predicts prognosis in patients with heart failure. Design Systematic review of studies assessing BNP for prognosis m patients with heart failure or asymptomatic patients. Data sources Electronic searches of Medline and Embase from January 1994 to March 2004 and reference lists of included studies. Study selection and data extraction We included all studies that estimated the relation between BNP measurement and the risk of death, cardiac death, sudden death, or cardiovascular event in patients with heart failure or asymptomatic patients, including initial values and changes in values in response to treatment. Multivariable models that included both BNP and left ventricular ejection fraction as predictors were used to compare the prognostic value of each variable. Two reviewers independently selected studies and extracted data. Data synthesis 19 studies used BNP to estimate the relative risk of death or cardiovascular events in heart failure patients and five studies in asymptomatic patients. In heart failure patients, each 100 pg/ml increase was associated with a 35% increase in the relative risk of death. BNP was used in 35 multivariable models of prognosis. In nine of the models, it was the only variable to reach significance-that is, other variables contained no prognostic information beyond that of BNP. Even allowing for the scale of the variables, it seems to be a strong indicator of risk. Conclusion Although systematic reviews of prognostic studies have inherent difficulties, including die possibility of publication bias, the results of the studies in this review show that BNP is a strong prognostic indicator for both asymptomatic patients mid for patients with heart failure at all stages of disease.

E7 oncoprotein of human papillomavirus type 16 expressed constitutively in the epidermis has no effect on E7-specific B- or Th-repertoires or on the immune response induced or sustained after immunization with E7 protein

A line of FVB (H-2(q)) mice transgenic for the E6/E7 open reading frames of Human Papillomavirus type 16 driven from the alpha-A crystallin promoter expresses E7 mRNA in lens and skin epithelium. E7 protein is detectable in adult skin, coinciding with the development or inflammatory skin disease, which progresses to papillomata and squamous carcinomata in some mice. By examining the outcome of parenteral immunization with E7 protein, we sought to determine whether endogenous expression of E7 in skin had induced a preexisting immune outcome, i.e., specific immunity or tolerance, or whether the mice remain naive (''ignorant'') to E7. Our data show that the antibody response to defined E7 B-epitopes, the proliferative response to Th epitopes, and the delayed-type hypersensitivity (DTH) response to whole E7 did not differ between groups or young and old E6/E7 transgenic mice (likely having different degrees of lifetime exposure to E7 protein) or between E6/E7-transgenic and nontransgenic parental strain control mice. Although an E7-specific CTL response could not be induced in the H-2(q) background of these mice, incorporation of a D-b allele into the genome allowed comparison of D-b-restricted CTL responses in E6/E7 transgenic and nontransgenic mice. Experiments indicated that the E7-immunization-induced CTL response did not differ significantly between E6/E7 transgenic and nontransgenic mice. We interpret these results to indicate that in spite of expression of E7 protein in adult skin, E6/E7 transgenic mice remain immunologically naive (ignorant) of E7 epitopes presented by immunization. (C) 1997 Academic Press.

Autopathogenic T helper cell type 1 (Th1) and protective Th2 clones differ in their recognition of the autoantigenic peptide of myelin proteolipid protein

We previously generated a panel of T helper cell 1 (Th1) clones specific for an encephalitogenic peptide of myelin proteolipid protein (PLP) peptide 139-151 (HSLGKWLGHPDKF) that induces experimental autoimmune encephalomyelitis (EAE) upon adoptive transfer. In spite of the differences in their T cell receptor (TCR) gene usage, all these Th1 clones required W144 as the primary and most critical TCR contact residue for the activation. In this study, we determined the TCR contact residues of a panel of Th2/Th0 clones specific for the PLP peptide 139-151 generated either by immunization with the PLP 139-151 peptide with anti-B7-1 antibody or by immunization with an altered peptide Q144. Using alanine-substituted peptide analogues of the native PLP peptide, we show that the Th2 clones have shifted their primary contact residue to the NH2-terminal end of the peptide. These Th2 cells do not show any dependence on the W144, but show a critical requirement for L141/G142 as their major TCR contact residue. Thus, in contrast with the Th1 clones that did not proliferate to A144-substituted peptide, the Th2 clones tolerated a substitution at position 144 and proliferated to A144 peptide. This alternative A144 reactive repertoire appears to have a critical role in the regulation of autoimmune response to PLP 139-151 because preimmunization with A144 to expand the L141/G142-reactive repertoire protects mice from developing EAE induced with the native PLP 139-151 peptide. These data suggest that a balance between two different T cell repertoires specific for same autoantigenic epitope can determine disease phenotype, i.e., resistance or susceptibility to an autoimmune disease.

The secondary passive film for type 304 stainless steel in 0.5 M H2SO4

An examination has been carried out of the secondary passive film on Type 304 stainless steel in 0.5 M H2SO4. The characterization techniques used were electrochemical (potentiodynamic; potentiostatic, and film reduction experiments) and surface analytical. A bilayer model for the secondary passive film is proposed. It appears that next to the metal, there is a modified passive film which controls the electrochemical response; i.e., governs the current for any applied potential. On top of this modified passive film, the experimental data are consistent with a ''porous'' corrosion-product film which adds to the total film thickness but has little influence on the electrochemical response. The composition of the secondary passive film corresponds most probably to a mixed Fe/Cr oxide/hydroxide enriched in Cr3+, With a composition similar to a primary passive film.