High-breakdown linear discriminant analysis

The classification rules of linear discriminant analysis are defined by the true mean vectors and the common covariance matrix of the populations from which the data come. Because these true parameters are generally unknown, they are commonly estimated by the sample mean vector and covariance matrix of the data in a training sample randomly drawn from each population. However, these sample statistics are notoriously susceptible to contamination by outliers, a problem compounded by the fact that the outliers may be invisible to conventional diagnostics. High-breakdown estimation is a procedure designed to remove this cause for concern by producing estimates that are immune to serious distortion by a minority of outliers, regardless of their severity. In this article we motivate and develop a high-breakdown criterion for linear discriminant analysis and give an algorithm for its implementation. The procedure is intended to supplement rather than replace the usual sample-moment methodology of discriminant analysis either by providing indications that the dataset is not seriously affected by outliers (supporting the usual analysis) or by identifying apparently aberrant points and giving resistant estimators that are not affected by them.

Mean anisotropy of homogeneous Gaussian random fields and anisotropic norms of linear translation-invariant operators on multidimensional integer lattices

Sensitivity of output of a linear operator to its input can be quantified in various ways. In Control Theory, the input is usually interpreted as disturbance and the output is to be minimized in some sense. In stochastic worst-case design settings, the disturbance is considered random with imprecisely known probability distribution. The prior set of probability measures can be chosen so as to quantify how far the disturbance deviates from the white-noise hypothesis of Linear Quadratic Gaussian control. Such deviation can be measured by the minimal Kullback-Leibler informational divergence from the Gaussian distributions with zero mean and scalar covariance matrices. The resulting anisotropy functional is defined for finite power random vectors. Originally, anisotropy was introduced for directionally generic random vectors as the relative entropy of the normalized vector with respect to the uniform distribution on the unit sphere. The associated a-anisotropic norm of a matrix is then its maximum root mean square or average energy gain with respect to finite power or directionally generic inputs whose anisotropy is bounded above by a≥0. We give a systematic comparison of the anisotropy functionals and the associated norms. These are considered for unboundedly growing fragments of homogeneous Gaussian random fields on multidimensional integer lattice to yield mean anisotropy. Correspondingly, the anisotropic norms of finite matrices are extended to bounded linear translation invariant operators over such fields.

Predicting residue-wise contact orders in proteins by support vector regression

Background: The residue-wise contact order (RWCO) describes the sequence separations between the residues of interest and its contacting residues in a protein sequence. It is a new kind of one-dimensional protein structure that represents the extent of long-range contacts and is considered as a generalization of contact order. Together with secondary structure, accessible surface area, the B factor, and contact number, RWCO provides comprehensive and indispensable important information to reconstructing the protein three-dimensional structure from a set of one-dimensional structural properties. Accurately predicting RWCO values could have many important applications in protein three-dimensional structure prediction and protein folding rate prediction, and give deep insights into protein sequence-structure relationships. Results: We developed a novel approach to predict residue-wise contact order values in proteins based on support vector regression (SVR), starting from primary amino acid sequences. We explored seven different sequence encoding schemes to examine their effects on the prediction performance, including local sequence in the form of PSI-BLAST profiles, local sequence plus amino acid composition, local sequence plus molecular weight, local sequence plus secondary structure predicted by PSIPRED, local sequence plus molecular weight and amino acid composition, local sequence plus molecular weight and predicted secondary structure, and local sequence plus molecular weight, amino acid composition and predicted secondary structure. When using local sequences with multiple sequence alignments in the form of PSI-BLAST profiles, we could predict the RWCO distribution with a Pearson correlation coefficient (CC) between the predicted and observed RWCO values of 0.55, and root mean square error (RMSE) of 0.82, based on a well-defined dataset with 680 protein sequences. Moreover, by incorporating global features such as molecular weight and amino acid composition we could further improve the prediction performance with the CC to 0.57 and an RMSE of 0.79. In addition, combining the predicted secondary structure by PSIPRED was found to significantly improve the prediction performance and could yield the best prediction accuracy with a CC of 0.60 and RMSE of 0.78, which provided at least comparable performance compared with the other existing methods. Conclusion: The SVR method shows a prediction performance competitive with or at least comparable to the previously developed linear regression-based methods for predicting RWCO values. In contrast to support vector classification (SVC), SVR is very good at estimating the raw value profiles of the samples. The successful application of the SVR approach in this study reinforces the fact that support vector regression is a powerful tool in extracting the protein sequence-structure relationship and in estimating the protein structural profiles from amino acid sequences.

Modifying insect population age structure to control vector-borne disease

Age is a critical determinant of the ability of most arthropod vectors to transmit a range of human pathogens. This is due to the fact that most pathogens require a period of extrinsic incubation in the arthropod host before pathogen transmission can occur. This developmental period for the pathogen often comprises a significant proportion of the expected lifespan of the vector. As such, only a small proportion of the population that is oldest contributes to pathogen transmission. Given this, strategies that target vector age would be expected to obtain the most significant reductions in the capacity of a vector population to transmit disease. The recent identification of biological agents that shorten vector lifespan, such as Wolbachia, entomopathogenic fungi and densoviruses, offer new tools for the control of vector-borne diseases. Evaluation of the efficacy of these strategies under field conditions will be possible due to recent advances in insect age-grading techniques. Implementation of all of these strategies will require extensive field evaluation and consideration of the selective pressures that reductions in vector longevity may induce on both vector and pathogen.

Evidence of a spotted fever-like rickettsia and a potential new vector from northeastern Australia

A spotted fever-like rickettsia was identified in a Hemaphysalis tick by polymerase chain reaction (PCR) amplification and sequencing of the 16S rDNA, ompA, and ompB genes. A comparison of these nucleotide sequences with those of other spotted fever group (SFG) rickettsiae revealed that the Hemaphysalis tick rickettsia was distinct from other previously reported strains. Phylogenetic analysis based on both ompA and ompB also indicates that the strain’s closest relatives are the agents of Thai tick typhus (Rickettsia honei strain TT-118) and Flinders Island spotted fever (R. honei). This study represents the first report of an R. honei-like agent from a Hemaphysalis tick in Australia and of a spotted fever group rickettsia from Cape York Peninsula, Queensland.

Field prevalence of Wolbachia in the mosquito vector Aedes albopictus

The endosymbiotic bacteria in the genus Wolbachia have been proposed as a potential candidate to deliver pathogen-blocking genes into natural populations of medically important insects. The successful application of Wolbachia in insect vector control depends on the ability of the agent to successfully invade and maintain itself at high frequency under field conditions. Here, we evaluated the prevalence of Wolbachia infections in a field population of the Wolbachia-superinfected mosquito Aedes albopictus. A field prevalence of 100% (n = 1,016) was found in a single population in eastern Thailand via polymerase chain reaction (PCR) testing of Wolbachia both from individual parent females and their corresponding F1 offspring. This is the first report of accurate Wolbachia prevalence in a field population of an insect disease vector. The prevalence of superinfection was estimated to be 99.41%. All single-infected individual mosquitoes (n = 6) were found to harbor group A Wolbachia. For this particular population, none was found to be single-infected with group B Wolbachia. Our results also show that PCR testing of field materials alone without checking F1 offspring overestimated the natural prevalence of single infection. Thus, the confirmation of infection status by means of F1 offspring was critical to the accurate estimates of Wolbachia prevalence under field conditions.

Prospects for control of African trypanosomiasis by tsetse vector manipulation

The extensive antigenic variation phenomena African trypanosomes display in their mammalian host have hampered efforts to develop effective vaccines against trypanosomiasis. Human disease management aims largely to treat infected hosts by chemotherapy, whereas control of animal diseases relies on reducing tsetse populations as well as on drug therapy. The control strategies for animal diseases are carried out and financed by livestock owners, who have an obvious economic incentive. Sustaining largely insecticide-based control at a local level and relying on drugs for treatment of infected hosts for a disease for which there is no evidence of acquired immunity could prove extremely costly in the long run. It is more likely that a combination of several methods in an integrated, phased and area-wide approach would be more effective in controlling these diseases and subsequently improving agricultural output. New approaches that are environmentally acceptable, efficacious and affordable are clearly desirable for control of various medically and agriculturally important insects including tsetse. Here, Serap Aksoy and colleagues discuss molecular genetic approaches to modulate tsetse vector competence.

Semliki Forest virus as an expression vector in insect cell lines

Studies were undertaken to determine if replication-deficient Semliki Forest virus expression vectors could be successfully used to express foreign gene constructs in insect cell lines. Using green fluorescent protein (GFP) as a marker we recorded infection levels of nearly 100% in the Aedes albopictus cell lines C6/36 and Aa23T, as well as in the Ae. aegypti cell line MOS20. The virus was capable of infecting an Anopheles gambiae cell line MOS55. The amount of GFP protein produced in each cell line was quantified. Northern analysis of viral transcription revealed the presence of novel transcripts in Aa23T, C6/36, and MOS55 cell lines, but not in the BHK or MOS20. The initial characterization of these transcripts is described.

Control of vector-borne disease by genetic manipulation of insect populations: technological requirements and research priorities

The possibility of controlling vector-borne disease through the development and release of transgenic insect vectors has recently gained popular support and is being actively pursued by a number of research laboratories around the world. Several technical problems must be solved before such a strategy could be implemented: genes encoding refractory traits (traits that render the insect unable to transmit the pathogen) must be identified, a transformation system for important vector species has to be developed, and a strategy to spread the refractory trait into natural vector populations must be designed. Recent advances in this field of research make it seem likely that this technology will be available in the near future. In this paper we review recent progress in this area as well as argue that care should be taken in selecting the most appropriate disease system with which to first attempt this form of intervention. Much attention is currently being given to the application of this technology to the control of malaria, transmitted by Anopheles gambiae in Africa. While malaria is undoubtedly the most important vector-borne disease in the world and its control should remain an important goal, we maintain that the complex epidemiology of malaria together with the intense transmission rates in Africa may make it unsuitable for the first application of this technology. Diseases such as African trypanosomiasis, transmitted by the tsetse fly, or unstable malaria in India may provide more appropriate initial targets to evaluate the potential of this form of intervention.

Modification of arthropod vector competence via symbiotic bacteria

Some of the world's most devastating diseases are transmitted by arthropod vectors. Attempts to control these arthropods are currently being challenged by the widespread appearance of insecticide resistance. It is therefore desirable to develop alternative strategies to complement existing methods of vector control. In this review, Charles Beard, Scott O'Neill, Robert Tesh, Frank Richards and Serap Aksoy present an approach for introducing foreign genes into insects in order to confer refractoriness to vector populations, ie. the inability to transmit disease-causing agents. This approach aims to express foreign anti-parasitic or anti-viral gene products in symbiotic bacteria harbored by insects. The potential use of naturally occurring symbiont-based mechanisms in the spread of such refractory phenotypes is also discussed.

Metabolite mean transit times in the liver as predicted by various models of hepatic elimination

Predicted area under curve (AUC), mean transit time (MTT) and normalized variance (CV2) data have been compared for parent compound and generated metabolite following an impulse input into the liver, Models studied were the well-stirred (tank) model, tube model, a distributed tube model, dispersion model (Danckwerts and mixed boundary conditions) and tanks-in-series model. It is well known that discrimination between models for a parent solute is greatest when the parent solute is highly extracted by the liver. With the metabolite, greatest model differences for MTT and CV2 occur when parent solute is poorly extracted. In all cases the predictions of the distributed tube, dispersion, and tasks-in-series models are between the predictions of the rank and tube models. The dispersion model with mixed boundary conditions yields identical predictions to those for the distributed tube model (assuming an inverse gaussian distribution of tube transit times). The dispersion model with Danckwerts boundary conditions and the tanks-in series models give similar predictions to the dispersion (mixed boundary conditions) and the distributed tube. The normalized variance for parent compound is dependent upon hepatocyte permeability only within a distinct range of permeability values. This range is similar for each model but the order of magnitude predicted for normalized variance is model dependent. Only for a one-compartment system is the MIT for generated metabolite equal to the sum of MTTs for the parent compound and preformed metabolite administered as parent.

Optimising fed-batch production of recombinant proteins using the baculovirus expression vector system

Fed-batch culture can offer significant improvement in recombinant protein production compared to batch culture in the baculovirus expression vector system (BEVS), as shown by Nguyen et al. (1993) and Bedard et al. (1994) among others. However, a thorough analysis of fed-batch culture to determine its limits in improving recombinant protein production over batch culture has yet to be performed. In this work, this issue is addressed by the optimisation of single-addition fed-batch culture. This type of fed-batch culture involves the manual addition of a multi-component nutrient feed to batch culture before infection with the baculovirus. The nutrient feed consists of yeastolate ultrafiltrate, lipids, amino acids, vitamins, trace elements, and glucose, which were added to batch cultures of Spodoptera frugiperda (Sf9) cells before infection with a recombinant Autographa californica nuclear polyhedrosis virus (Ac-NPV) expressing beta-galactosidase (beta-Gal). The fed-batch production of beta-Gal was optimised using response surface methods (RSM). The optimisation was performed in two stages, starting with a screening procedure to determine the most important variables and ending with a central-composite experiment to obtain a response surface model of volumetric beta-Gal production. The predicted optimum volumetric yield of beta-Gal in fed-batch culture was 2.4-fold that of the best yields in batch culture. This result was confirmed by a statistical analysis of the best fed-batch and batch data (with average beta-Gal yields of 1.2 and 0.5 g/L, respectively) obtained from this laboratory. The response surface model generated can be used to design a more economical fed-batch operation, in which nutrient feed volumes are minimised while maintaining acceptable improvements in beta-Gal yield. (C) 1998 John Wiley & Sons, Inc.

The velocity of probability transport in quantum mechanics

In his study of the 'time of arrival' problem in the nonrelativistic quantum mechanics of a single particle, Allcock [1] noted that the direction of the probability flux vector is not necessarily the same as that of the mean momentum of a wave packet, even when the packet is composed entirely of plane waves with a common direction of momentum. Packets can be constructed, for example for a particle moving under a constant force, in which probability flows for a finite time in the opposite direction to the momentum. A similar phenomenon occurs for the Dirac electron. The maximum amount of probabilitiy backflow which can occur over a given time interval can be calculated in each case.

Gene transfer and expression of a non-viral polycation-based vector in CD4(+) cells

CD4-selective targeting of an antibody-polycation-DNA complex was investigated The complex was synthesized with the anti-CD4 monoclonal antibody B-F5, polylysine(268) (pLL) and either the pGL3 control vector containing the luciferase reporter gene or the pGeneGrip vector containing the green fluorescent protein (GFP) gene. B-F5-pLL-DNA complexes inhibited the binding of I-125-B-F5 to CD4(+) Jurkat cells, while complexes synthesised either without B-F5 or using a non-specific mouse IgG1 antibody had little or no effect Expression of the luciferase reporter gene was achieved in Jurkat cells using the B-F5-pLL-pGL3 complex and was enhanced in the presence of PMA. Negligible luciferase activity was defected with the non-specific antibody complex in Jurkat cells or with the B-F5-pLL-pGL3 complex in the CD4(-) K-562 cells. Using complexes synthesised with the pGeneGrip vector, the transfection efficiency in Jurkat and K-562 cells was examined using confocal microscopy. More than 95% of Jurkat cells expressed GFP and the level of this expression was markedly enhanced by PMA. Negligible GFP expression was seen in K-562 cells or when B-F5 was replaced by a nonspecific antibody. Using flow cytometry, fluorescein-labelled complex showed specific targeting to CD4(+) cells in a mixed cell population from human peripheral blood. These studies demonstrate the selective transfection of CD4(+) T-lymphoid cells using a polycation-based gene delivery system. The complex may provide a means of delivering anti-HIV gene therapies to CD4(+) cells in vivo.

Transport properties of strongly correlated metals: A dynamical mean-field approach

The temperature dependence of the transport properties of the metallic phase of a frustrated Hubbard model on the hypercubic lattice at half-filling is calculated. Dynamical mean-held theory, which maps the Hubbard model onto a single impurity,Anderson model that is solved self-consistently, and becomes exact in the limit of large dimensionality, is used. As the temperature increases there is a smooth crossover from coherent Fermi liquid excitations at low temperatures to incoherent excitations at high temperatures. This crossover leads to a nonmonotonic temperature dependence for the resistance, thermopower, and Hall coefficient, unlike in conventional metals. The resistance smoothly increases from a quadratic temperature dependence at low temperatures to large values which can exceed the Mott-Ioffe-Regel value ha/e(2) (where a is a lattice constant) associated with mean free paths less than a lattice constant. Further signatures of the thermal destruction of quasiparticle excitations are a peak in the thermopower and the absence of a Drude peak in the optical conductivity. The results presented here are relevant to a wide range of strongly correlated metals, including transition metal oxides, strontium ruthenates, and organic metals.