Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Origin and Rapid Evolution of a Novel Murine Erythroleukemia Virus of the Spleen Focus-Forming Virus Family

The Friend spleen focus-forming virus (SFFV) env gene encodes a glycoprotein with apparent M_r of 55,000 that binds to erythropoietin receptors (EpoR) to stimulate erythroblastosis. A retroviral vector that does not encode any Env glycoprotein was packaged into retroviral particles and was coinjected into mice in the presence of a nonpathogenic helper virus. Although most mice remained healthy, one mouse developed splenomegaly and polycythemia at 67 days; the virus from this mouse reproducibly caused the same symptoms in secondary recipients by 2 to 3 weeks postinfection. This disease, which was characterized by extramedullary erythropoietin-independent erythropoiesis in the spleens and livers, was also reproduced in long-term bone marrow cultures. Viruses from the diseased primary mouse and from secondary recipients converted an erythropoietin-dependent cell line (BaF3/EpoR) into factor-independent derivatives but had no effect on the interleukin-3-dependent parental BaF3 cells. Most of these factor-independent cell clones contained a major Env-related glycoprotein with an M_r of 60,000. During further in vivo passaging, a virus that encodes an M_r-55,000 glycoprotein became predominant. Sequence analysis indicated that the ultimate virus is a new SFFV that encodes a glycoprotein of 410 amino acids with the hallmark features of classical gp55s. Our results suggest that SFFV-related viruses can form in mice by recombination of retroviruses with genomic and helper virus sequences and that these novel viruses then evolve to become increasingly pathogenic.The independently isolated Friend and Rauscher erythroleukemia viruses (18, 48) are complexes of a replication competent murine leukemia virus (MuLV) and a replication-defective spleen focus-forming virus (SFFV) (39, 42, 47). The SFFVs encode Env glycoproteins (gp55) that are inefficiently processed to form larger cell surface derivatives (gp55^p) (19). The gp55 binds to erythropoietin receptors (EpoR) to cause erythroblast proliferation and splenomegaly in susceptible mice. Evidence has suggested that the critical mitogenic interaction occurs exclusively on cell surfaces (7, 16).SFFVs are structurally closely related to mink cell focus-inducing viruses (MCFs) (2, 5, 10, 50), a class of replication-competent murine retroviruses that has a broad host range termed polytropic (15, 21). Although MCFs are not inherited as replication-competent intact proviruses, the mouse genome contains numerous dispersed polytropic env gene sequences (8, 17, 27). MCFs apparently readily form de novo by recombination when ecotropic host range MuLVs replicate in mice (14, 15, 26, 43). MCF env genes typically are hybrid recombinants that contain a 5′ polytropic-specific region and a 3′ ecotropic-specific portion (26). They encode a gPr90 Env glycoprotein that is cleaved by partial proteolysis to form the products gp70 surface (SU) glycoprotein plus p15E transmembrane (TM) protein (32, 39, 47). In addition, MCFs often differ from ecotropic MuLVs in their long terminal repeat (LTR) sequences (8, 20, 26, 28, 29, 45).Based on their sequences, SFFVs could have derived from MCFs by a 585-base deletion and by a single-base addition in the ecotropic-specific portion of the env gene (10). Evidence suggests that both the 585-bp deletion and the frameshift mutation probably contribute to SFFV pathogenesis (3, 49). Several pathogenic differences among SFFV strains have also been ascribed to amino acid sequence differences in the ecotropic-specific portion of the Env glycoproteins (9, 41).This report describes the origin and rapid stepwise evolution of a new SFFV. This new pathogenic virus initially formed in a mouse that had been injected with an ecotropic strain of MuLV in the presence of a retroviral vector that does not encode any Env glycoprotein. The mouse developed erythroleukemia, splenomegaly, and polycythemia after a long lag phase. At that time the spleen contained viruses with env genes that were able to activate EpoR. Serial passage of this initial virus isolate resulted in selection of a novel SFFV that encodes a gp55 glycoprotein of 410 amino acids. This experimental system provides a method for isolating new SFFVs and for mapping the stages in their genesis.     

Increased Lipid Accumulation in the Chlamydomonas reinhardtii sta7-10 Starchless Isoamylase Mutant and Increased Carbohydrate Synthesis in Complemented Strains

The accumulation of bioenergy carriers was assessed in two starchless mutants of Chlamydomonas reinhardtii (the sta6 [ADP-glucose pyrophosphorylase] and sta7-10 [isoamylase] mutants), a control strain (CC124), and two complemented strains of the sta7-10 mutant. The results indicate that the genetic blockage of starch synthesis in the sta6 and sta7-10 mutants increases the accumulation of lipids on a cellular basis during nitrogen deprivation relative to that in the CC124 control as determined by conversion to fatty acid methyl esters. However, this increased level of lipid accumulation is energetically insufficient to completely offset the loss of cellular starch that is synthesized by CC124 during nitrogen deprivation. We therefore investigated acetate utilization and O₂ evolution to obtain further insights into the physiological adjustments utilized by the two starchless mutants in the absence of starch synthesis. The results demonstrate that both starchless mutants metabolize less acetate and have more severely attenuated levels of photosynthetic O₂ evolution than CC124, indicating that a decrease in overall anabolic processes is a significant physiological response in the starchless mutants during nitrogen deprivation. Interestingly, two independent sta7-10:STA7 complemented strains exhibited significantly greater quantities of cellular starch and lipid than CC124 during acclimation to nitrogen deprivation. Moreover, the complemented strains synthesized significant quantities of starch even when cultured in nutrient-replete medium.Microalgae are able to efficiently convert sunlight, water, and CO₂ into a variety of products suitable for renewable energy applications, including H₂, carbohydrates, and lipids (11, 12, 16, 21, 38, 41, 44). The unicellular green alga Chlamydomonas reinhardtii has emerged as a model organism for studying algal physiology, photosynthesis, metabolism, nutrient stress, and the synthesis of bioenergy carriers (12, 15, 19, 24, 32). During acclimation to nitrogen deprivation, C. reinhardtii cells accumulate significant quantities of starch and form lipid bodies (4, 5, 8, 26, 28, 30, 34, 43, 46, 48). Despite the significance of these products in algal physiology and in biofuels applications, the metabolic, enzymatic, and regulatory mechanisms controlling the partitioning of metabolites into these distinct carbon stores in algae are poorly understood. Several C. reinhardtii starch mutants with various phenotypic changes in starch content and structure have been isolated (2,–4). Two of these, the sta6 and sta7 mutants, contain single-gene disruptions that result in “starchless” phenotypes with severely attenuated levels of starch granule accumulation (2, 4, 34, 39, 40, 48).The disrupted loci in the two isolated starchless mutants are distinct and each mutant has a unique phenotype (7, 40). In the sta6 mutant, the small, catalytic subunit of ADP-glucose pyrophosphorylase (AGPase-SS) is disrupted (2, 4, 48), and this mutant accumulates less than 1% of the starch observed in wild-type (WT) cells under conditions of nitrogen deprivation. The sta7 mutant contains a disrupted isoamylase gene (7, 8, 10, 39, 40) and also has severely attenuated levels of starch, but it accumulates a soluble glycogen-like product (4, 9). In this study, we conducted an examination of the unique physiological acclimations that are utilized by these mutants to adapt to the loss of starch synthesis. As the genetic lesions in these two mutants are distinct and block starch synthesis via two very different mechanisms, we investigated the physiological consequences of starch inhibition in both of these mutants from a holistic bioenergy perspective, which included photosynthetic parameters and the overall yields of lipids and carbohydrates, the two primary bioenergy carriers in C. reinhardtii. Specifically, we examined whether the inability to synthesize starch would result in the accumulation of additional lipid, alter cellular growth or cell size, affect acetate utilization, and/or influence photosynthetic O₂ evolution. Our data indicate that both the sta6 (BAFJ5) and sta7 (sta7-10) mutants accumulate more lipid than the CC124 control during nitrogen deprivation. However, the additional lipid does not completely offset the loss of starch synthesis from a complete energetic perspective. Increased lipid accumulation during nitrogen stress has also been reported for a variety of starch mutants in recent papers (26, 27, 46). A significant feature in both of the starchless mutants studied here is that O₂ evolution and acetate utilization are diminished during nitrogen stress, which is undesirable from an overall bioenergy perspective. Remarkably, complementation of sta7-10 with genomic DNA encoding the wild-type isoamylase gene resulted in cells that were larger than those of the sta6, sta7-10, and CC124 strains, exhibited the highest total lipid levels during nitrogen deprivation, and overaccumulated starch even in nutrient-replete medium.     

Differentiation of Lactobacillus Species by Molecular Typing

A total of 64 type, reference, clinical, health food, and stock isolates of microaerophilic Lactobacillus species were examined by restriction fragment length polymorphisms. Of particular interest were members of six of the eight species most commonly recovered from the vaginas of healthy premenopausal women, namely, Lactobacillus jensenii, L. casei, L. rhamnosus, L. acidophilus, L. plantarum, and L. fermentum. Six main groupings were identified on the basis of ribotyping. This technique was able to classify fresh isolates to the species level. In the case of the ribotype A grouping for L. rhamnosus, differences between strains were evident by chromosome typing (chromotyping). Many isolates did not possess plasmids. Six L. rhamnosus strains isolated from four different health food products appeared to be identical to L. rhamnosus ATCC 21052. The molecular typing system is useful for identifying and differentiating Lactobacillus isolates. Studies of strains of potential importance to the urogenital flora should include molecular characterization as a means of comparing genetic traits with those of strains whose characteristics associated with colonization and antagonism against pathogens have been defined.Lactobacilli colonizing human tissues have long been considered important for the maintenance of a healthy gastrointestinal tract (4) and urogenital tract (20, 21, 22, 25). The disruption of the lactobacillus flora has been associated with many urogenital infections, and as these afflict over 150 million women worldwide each year, this area of study is an important one. Indeed, many patients resort to taking health food products containing lactobacilli as a means of trying to maintain a healthy intestinal (and in some cases vaginal) flora.The typing of lactobacilli has generally been conducted by cell and colony morphology and biochemical tests. These techniques type bacteria based on their ability to ferment sugars and produce acids such as lactic acid and acetic acid (9). Unfortunately, these typing methods are not completely accurate, and strains which show intermediate characteristics are frequently encountered (9, 33).Many studies emphasize that the classification of lactobacilli is unsatisfactory and does not reflect the real phylogenetic relatedness of different strains and species (6, 19, 30). Several new genetic and chemotaxonomic approaches have been used during the last 14 years with an aim of improving the classification and identification of lactobacilli: for example, analysis of plasmid content (18), sodium dodecyl sulfate-polyacrylamide gel electrophoresis patterns of whole-cell protein (19) and of total soluble cell protein (8), sequencing of rRNA (5, 6, 19), restriction endonuclease fingerprinting (14, 30), and DNA-DNA hybridization (6, 17, 19, 27). All of these approaches have improved the taxonomic knowledge of the generic and suprageneric relationships of lactobacilli. However, no analysis of urogenital or health food isolates of lactobacilli by molecular typing has been reported.Plasmid typing of Lactobacillus strains has been suggested as a taxonomic tool in a number of reports (27, 31), but there is evidence to suggest that it is not very effective (2, 12, 15, 28), since plasmids can be absent or unstable.Chromosome typing (restriction endonuclease fingerprinting of chromosomal DNA) (chromotyping) has been applied to the discrimination of strains of lactobacilli and has been found more specific and reproducible than plasmid content analysis (14). However, one of the disadvantages of chromotyping is that comparing electrophoretic patterns consisting of up to 100 bands is difficult. Hence, chromotyping alone may not be widely used for typing large numbers of strains.A sensitive but rapid method for species differentiation involves the use of ribotyping (1, 27), which combines Southern hybridization of chromosomal DNA fingerprints with the use of Escherichia coli rRNA probes, thereby discriminating between various species and individual strains of lactobacilli. The method involves the separation and identification of the rRNA genes present within the bacterial genome, which exist in various copy numbers (10, 13), making it possible to delineate species based on differences in the restriction fragment length polymorphisms of the rRNA genes. Of particular importance is the fact that specific regions of the rRNA genes have remained well conserved because of their functional importance, thus allowing the detection of a broad range of bacteria with 16S and 23S rRNA of E. coli as probes.The aim of this study was to develop a method to test the efficacy of ribotyping of Lactobacillus type and reference strains and to use this method to characterize a number of clinical, health food, and laboratory isolates.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Genomic Plasticity of the Human Fungal Pathogen Candida albicans

The genomic plasticity of Candida albicans, a commensal and common opportunistic fungal pathogen, continues to reveal unexpected surprises. Once thought to be asexual, we now know that the organism can generate genetic diversity through several mechanisms, including mating between cells of the opposite or of the same mating type and by a parasexual reduction in chromosome number that can be accompanied by recombination events (2, 12, 14, 53, 77, 115). In addition, dramatic genome changes can appear quite rapidly in mitotic cells propagated in vitro as well as in vivo. The detection of aneuploidy in other fungal pathogens isolated directly from patients (145) and from environmental samples (71) suggests that variations in chromosome organization and copy number are a common mechanism used by pathogenic fungi to rapidly generate diversity in response to stressful growth conditions, including, but not limited to, antifungal drug exposure. Since cancer cells often become polyploid and/or aneuploid, some of the lessons learned from studies of genome plasticity in C. albicans may provide important insights into how these processes occur in higher-eukaryotic cells exposed to stresses such as anticancer drugs.The purpose of this review is to describe the tools used to detect genome changes, to highlight recent advances in our understanding of large-scale chromosome changes that arise in Candida albicans, and to discuss the role of specific stresses in eliciting these genome changes. The types of genomic diversity that have been characterized suggest that C. albicans can undergo extreme genomic changes in order to survive stresses in the human host. We propose that C. albicans and other pathogens may have evolved mechanisms not only to tolerate but also to generate large-scale genetic variation as a means of adaptation.C. albicans is a polymorphic yeast with a 16-Mb (haploid) genome organized in 8 diploid chromosomes (140, 154, 203). The C. albicans genome displays a very high degree of plasticity. This plasticity includes the types of genomic changes frequently observed with cancer cells, including gross chromosomal rearrangements, aneuploidy, and loss of heterozygosity (reviewed in references 100, 117, and 157). Similar to somatic cancer cells, C. albicans reproduces primarily through asexual clonal division (65, 84). Nonetheless, it has retained much of the machinery needed for mating and meiosis (189), yet meiosis has never been observed (13, 120).C. albicans has two mating-type-like (MTL) alleles, MTLa and MTLα (76). The MTL locus is on the left arm of chromosome 5 (Chr5), approximately 80 kbp from the centromere. Most C. albicans isolates are heterozygous for the MTL locus, but approximately 3 to 10% of clinical isolates are naturally homozygous at MTL (104, 108). Mating can occur between strains carrying the opposite MTL locus, and most strains that were found to be naturally MTL homozygous are mating competent (104, 108). MTL-homozygous strains were also constructed from MTL-heterozygous strains by deletion of either the MTLa or MTLα locus (77) or by selection for Chr5 loss on sorbose (87, 115).Mating between these diploid strains of opposite mating type can occur both in vitro (115) and in vivo (77, 97). The products are tetraploid and do not undergo a conventional meiotic reduction in ploidy (12, 120). Rather, they undergo random loss of multiple chromosomes, a process termed “concerted chromosome loss,” until they reach a near-diploid genome content (2, 12, 53, 85). A subset of these cells also undergoes multiple gene conversion events reminiscent of meiotic recombination, and most remain trisomic for one to several chromosomes (53). While mating and concerted chromosome loss have been induced in the laboratory, the role of the parasexual cycle during the host-pathogen interaction and in the response to stresses, such as exposure to antifungal drugs, remains unclear. The prevailing model is that adaptive mutations (such as those that occur with the acquisition of drug resistance) evolve through somatic events, including point mutations, recombination, gene conversion, loss of heterozygosity, and/or aneuploidy (13).