A simulated annealing algorithm for finding consensus sequences

Motivation: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. Results: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.

Cleavage of hemoglobin by hookworm cathespin D aspartic proteases and its potential contribution to host specificity

Hookworms routinely reach the gut of nonpermissive hosts but fail to successfully feed, develop, and reproduce. To investigate the effects of host-parasite coevolution on the ability of hookworms to feed in nonpermissive hosts, we cloned and expressed aspartic proteases from canine and human hookworms. We show here that a cathepsin D-like protease from the canine hookworm Ancylosotoma caninum (Ac-APR-1) and the orthologous protease from the human hookworm Necator americanus (Na-APR-1) are expressed in the gut and probably exert their proteolytic activity extracellularly. Both proteases were detected immunologically and enzymatically in somatic extracts of adult worms. The two proteases were expressed in baculovirus, and both cleaved human and dog hemoglobin (Hb) in vitro. Each protease digested Hb from its permissive host between twofold (whole molecule) and sixfold (synthetic peptides) more efficiently than Hb from the nonpermissive host, despite the two proteases' having identical residues lining their active site clefts. Furthermore, both proteases cleaved Hb at numerous distinct sites and showed different substrate preferences. The findings suggest that the paradigm of matching the molecular structure of the food source within a host to the molecular structure of the catabolic proteases of the parasite is an important contributing factor for host-parasite compatibility and host species range.

Evidence for an independent radiation of endosymbiotic litostome ciliates within Australian marsupial herbivores

The recent discovery of isotrichid-like ciliates occurring as endosymbionts in macropodid marsupials posed interesting questions in regard to both their phyletic origin (all previous records confined to eutherian mammals) and their morphological evolution (Australian forms possibly representing missing links between previously described genera). The SSU rRNA gene was sequenced for three species (Dasytricha dehorityi, D. dogieli, and Batricha tasmaniensis) and aligned against representatives of all major ciliate classes. The Australian species did not group with the other isotrichid species but instead formed an independent radiation. Discrepancies between recent global phylogenies of the phylum Ciliophora were examined by manipulation of the aligned sequence data set. Sources of conflict between these studies did not stem from differences in outgroup choice or phylogenetic reconstruction methods. Differences in the application of confidence limits and primary sequence alignment have probably resulted in the reporting of spurious associations which are not supported by more conservative confidence or alignment methodology. At present, the ciliate subphylum Intramacro-nucleata is an unresolved polytomy which may be due to deficiencies in the SSU rRNA gene sequence dataset or indicate that the ciliates radiated into their extant classes by rapid burst-like evolution. (C) 2001 academic Press.

Total synthesis and absolute stereochemistry of plakortone D

The first total synthesis of plakortone D is described and thereby establishes the structure and absolute stereochemistry of the most biologically active member of the marine-derived plakortone family. The sterically congested bicyclic lactone core results from a Pd(II)-induced hydroxycyclization−carbonylation−lactonization sequence on an enediol whose chirality was installed by AD-technology. Attachment of the side chain, also constructed using AD-methodology, was achieved by using a modified Julia coupling. The described approach enables acquisition of other plakortones and analogues, in the correct (natural) stereochemical series.

Cloning of a differentially expressed tropomyosin isoform from cultured rabbit aortic smooth muscle cells

The four known tropomyosin genes have highly conserved DNA and amino acid sequences, and at least 18 isoforms are generated by alternative RNA splicing in muscle and non-muscle cells. No rabbit tropomyosin nucleotide sequences are known, although protein sequences for alpha- and beta-tropomyosin expressed by rabbit skeletal muscle have been described. Subtractive hybridisation was used to select for genes differentially expressed in rabbit aortic smooth muscle cells (SMC), during the change in cell phenotype in primary culture that is characterised by a loss of cytoskeletal filaments and contractile proteins. This led to the cloning of a tropomyosin gene predominantly expressed in rabbit SMC during this change. The full-length cDNA clone, designated rabbit TM-beta, contains an open reading frame of 284 amino acids, 5' untranslated region (UTR) of I 17 base pairs and 3' UTR of 79 base pairs. It is closely related to the beta-gene isoforms in other species, with the highest homology in DNA and protein sequences to the human fibroblast isoform TM-1 (91.7% identity in 1035 bp and 93.3% identity in the entire 284 amino acid sequence of the protein), It differs from rabbit skeletal muscle P-tropomyosin (81.7% homology at the protein level) mainly in two regions at amino acids 189-213 and 258-283 suggesting alternative splicing of exons 6a for 6b and 9d for 9a. Since this TM-P gene was the only gene strongly enough expressed in SMC changing phenotype to be observed by the subtractive hybridisation screen, it likely plays a significant role in this process. (C) 2002 Published by Elsevier Science Ltd.

The nexus of value chain integration and e-business applications on public sector agriculture R&D management

We examine the potential impact of interconnectivity of value chain partnerships through electronic means (e-business practices) on the management of Public Sector Agriculture R&D in Australia. We review the changing forms of managing research and development, the forces driving these changes, and R&D processes that are theoretically consistent with the move towards value chain involvement and the increase in active constituents in Public Sector Agriculture R&D. We then explore the potential of emerging e-business models to change the patterns of inter-connectivity, speed and omnipresence of partners in the value chain. Three e-business R&D management practices are identified that provide the prerequisite flexibility necessary to take advantage of opportunistic markets. These R&D business practices are: compressing R&D to reduce time to market, fostering co-development to enter a market at the last moment and building flexible products that allow adjustment at the last possible moment. Some fundamental reallocation of existing resources will be required to meet these markets. Implications of these e-business practices for R&D management are discussed.

Characterization of mouse profilaggrin: Evidence for nuclear engulfment and translocation of the profilaggrin B-domain during epidermal differentiation

Filaggrin is a keratin filament associated protein that is expressed in granular layer keratinocytes and derived by sequential proteolysis from a polyprotein precursor termed profilaggrin. Depending on the species, each profilaggrin molecule contains between 10 and 20 filaggrin subunits organized as tandem repeats with a calcium-binding domain at the N-terminal end. We now report the characterization of the complete mouse gene. The structural organization of the mouse gene is identical to the human profilaggrin gene and consists of three exons with a 4 kb intron within the 5' noncoding region and a 1.7 kb intron separating the sequences encoding the calcium-binding EF-hand motifs. A processed pseudogene was found embedded within the second intron. The third and largest exon encodes the second EF-hand, a basic domain (designated the B-domain) followed by 12 filaggrin repeats and a unique C-terminal tail domain. A polyclonal anti-body raised against the conceptually translated sequence of the B-domain specifically stained keratohyalin granules and colocalized with a filaggrin antibody in granular layer cells. In upper granular layer cells, B-domain containing keratohyalin granules were in close apposition to the nucleus and, in some cells, appeared to be completely engulfed by the nucleus. In transition layer cells, B-domain staining was evident in the nucleus whereas filaggrin staining remained cytoplasmic. Nuclear staining of the B-domain was also observed in primary mouse keratinocytes induced to differentiate. This study has also revealed significant sequence homology between the mouse and human promoter sequences and in the calcium-binding domain but the remainder of the protein-coding region shows substantial divergence.

Improved detection of lacZ reporter gene expression in transgenic epithelia by immunofluorescence microscopy

The bacterial lacZ gene is commonly used as a reporter for the in vivo analysis of gene regulation in transgenic mice. However, several laboratories have reported poor detection of beta-galactosidase (the lacZ gene product) using histochemical techniques, particularly in skin. Here we report the difficulties we encountered in assessing lacZ expression in transgenic keratinocytes using classic X-gal histochemical protocols in tissues shown to express the transgene by mRNA in situ hybridization. We found that lacZ reporter gene expression could be reliably detected in frozen tissue sections by immunofluorescence analysis using a beta-galactosidase-specific antibody. Moreover, we were able to localize both transgene and endogenous gene products simultaneously using double-label immunofluorescence. Our results suggest that antibody detection of beta-galactosidase should be used to verify other assays of lacZ expression, particularly where low expression levels are suspected or patchy expression is observed.

Urine D-arabinitol/L-arabinitol ratio in diagnosing Candida infection in patients with haematological malignancy and HIV infection

Adult patients with hematologic malignancies along with HIV infected patients were prospectively studied to determine the performance of urine D-arabinitol/L-arabinitol (DA/LA) ratio in diagnosing invasive candidiasis. Ten evaluable febrile neutropenic patients had proven invasive candidiasis and elevated DA/LA ratios were found in 5. Invasive candidiasis with normal DA/LA ratios was most frequently due to Candida krusei infection. This Candida species is a non-producer of arabinitol. Only 4 of 81 febrile neutropenic patients given either antifungal prophylaxis or empiric antifungal treatment had elevated DA/LA ratios. Only 1 of 15 HIV positive patients with either oropharyngeal or esophageal candidiasis had elevated DA/LA ratios. Widespread use of fluconazole prophylaxis in bone marrow transplantation patients at the study hospital has led to an increased prevalence of C. krusei infection. This is the likely reason for the low sensitivity of the test in proven and suspected invasive Candida infections reported here. (C) 2002 Elsevier Science Inc. All rights reserved.

Alterations in Plasmodium falciparum genotypes during sequential infections suggest the presence of strain specific immunity

Many of the asexual stage Plasmodium falciparum proteins that are the targets of host protective responses are markedly polymorphic. The full repertoire of diversity is not defined for any antigen. Most studies have focused on the genes encoding merozoite surface proteins 1 and 2 (MSP1, MSP2). We explored the extent of diversity of some of the less studied merozoite surface antigens and analyzed the degree of complexity of malaria field isolates by deriving nucleotide sequences of several antigens. We have determined the genotype of apical membrane antigen 1 (AMA1) in a group of 30 field samples, collected over 29 months, from individuals living in an area of intense malaria transmission in Irian Jaya, identifying 14 different alleles. AMA1 genotyping was combined with previously determined MSP2 typing. AMA1 had the greatest power in distinguishing between isolates but methodological problems, especially when mixed infections are present, suggest it is not an ideal typing target. MSP1, MSP3, and glutamate-rich protein genotypes were also determined from a smaller group of samples, and all results were combined to derive an extended antigenic haplotype. Within this subset of 10 patients, nine different genotypes could be discerned; however, five patients were all infected with the same strain. This strain was present in individuals from two separate villages and was still present 12 months later. This strain was predominant at the first time point but had disappeared at the fourth time point. This significant change in malaria genotypes could be due to strain-specific immunity developing in this population.

Cytochrome P450(cin) (CYP176A), isolation, expression, and characterization

Cytochromes P450 are members of a superfamily of hemoproteins involved in the oxidative metabolism of various physiologic and xenobiotic compounds in eukaryotes and prokaryotes. Studies on bacterial P450s, particularly those involved in monoterpene oxidation, have provided an integral contribution to our understanding of these proteins, away from the problems encountered with eukaryotic forms. We report here a novel cytochrome P450 (P450(cin), CYP176A1) purified from a strain of Citrobacter braakii that is capable of using cineole 1 as its sole source of carbon and energy. This enzyme has been purified to homogeneity and the amino acid sequences of three tryptic peptides determined. By using this information, a PCR-based cloning strategy was developed that allowed the isolation of a 4-kb DNA fragment containing the cytochrome P450(cin) gene (cinA). Sequencing revealed three open reading frames that were identified on the basis of sequence homology as a cytochrome P450, an NADPH-dependent flavodoxin/ferrodoxin reductase, and a flavodoxin. This arrangement suggests that P450(cin) may be the first isolated P450 to use a flavodoxin as its natural redox partner. Sequencing also identified the unprecedented substitution of a highly conserved, catalytically, important active site threonine with an asparagine residue. The P450 gene was subcloned and heterologously expressed in Escherichia coli at similar to2000 nmol/liter of original culture, and purification was achieved by standard protocols. Postulating the native E. coli flavodoxin/flavodoxin reductase system might mimic the natural redox partners of P450,in, it was expressed in E. coli in the presence of cineole 1. A product was formed in vivo that was tentatively identified by gas chromatography-mass spectrometry as 2-hydroxycineole 2. Examination of P450(cin) by UV-visible spectroscopy revealed typical spectra characteristic of P450s, a high affinity for cineole 1 (K-D = 0.7 mum), and a large spin state change of the heme iron associated with binding of cineole 1. These facts support the hypothesis that cineole 1 is the natural substrate for this enzyme and that P450(cin) catalyzes the initial monooxygenation of cineole 1 biodegradation. This constitutes the first characterization of an enzyme involved in this pathway.

The ultrastructure of Amylovorax dehorityi comb. nov and erection of the Amylovoracidae fam. nov (Ciliophora : Trichostomatia)

The ultrastructural features of the holotrichous ciliates inhabiting macropodid maruspials were investigated to resolve their morphological similarity to other trichostome ciliates with observed differences in their small subunit rRNA gene sequences. The ultrastructure of Amylovorax dehorityi nov. comb. (formerly Dasytricha dehorityi) was determined by transmission electron microscopy. The somatic kineties are composed of monokinetids whose microtubules show a typical litostome pattern. The somatic cortex is composed of ridges which separate kinety rows, granular ectoplasm and a basal layer of hydrogenosomes lining the tela corticalis. The vestibulum is an invagination of the pellicle lined down one side with kineties (invaginated extensions of the somatic kineties); transverse tubules line the surface of the vestibulum and small nematodesmata surround it forming a cone-like network of struts. Cytoplasmic organelles include hydrogenosomes, irregularly shaped contractile vacuoles surrounded by a sparse spongioplasm, food vacuoles containing bacteria and large numbers of starch granules. This set of characteristics differs sufficiently from those of isotrichids and members of the genus Dasytricha to justify the erection of a new genus (Amylovorax) and a new family (Amylovoracidae). Dasytricha dehorityi, D. dogieli and D. mundayi are reassigned to the new genus Amylovorax and a new species A. quokka is erected. While the gross morphological similarities between Amylovorax and Dasytricha may be explained by convergent evolution, ultrastructural features indicate that these two genera have probably diverged independently from haptorian ancestors by successive reduction of the cortical and vestibular support structures.