STEME: efficient EM to find motifs in large data sets

MEME and many other popular motif finders use the expectation-maximization (EM) algorithm to optimize their parameters. Unfortunately, the running time of EM is linear in the length of the input sequences. This can prohibit its application to data sets of the size commonly generated by high-throughput biological techniques. A suffix tree is a data structure that can efficiently index a set of sequences. We describe an algorithm, Suffix Tree EM for Motif Elicitation (STEME), that approximates EM using suffix trees. To the best of our knowledge, this is the first application of suffix trees to EM. We provide an analysis of the expected running time of the algorithm and demonstrate that STEME runs an order of magnitude more quickly than the implementation of EM used by MEME. We give theoretical bounds for the quality of the approximation and show that, in practice, the approximation has a negligible effect on the outcome. We provide an open source implementation of the algorithm that we hope will be used to speed up existing and future motif search algorithms.     

Social Immunity in Amphibians: Evidence for Vertical Transmission of Innate Defenses

Amphibian embryos are at risk of microbial infection. Here we find evidence that innate immune defenses, both antimicrobial skin peptides and mutualistic microbiota, of adult glass frogs, Hyalinobatrachium colymbiphyllum, can be transmitted to embryos deposited on leaves above rain forest streams in Panama and can inhibit the fungal pathogen Batrachochytrium dendrobatidis. Abstract in Spanish is available at http://www.blackwell‐synergy.com/loi/btp .     

Chromosome territory arrangement and homologous pairing in nuclei of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> are predominantly random except for NOR-bearing chromosomes

Differential painting of all five chromosome pairs of Arabidopsis thaliana revealed for the first time the interphase chromosome arrangement in a euploid plant. Side-by-side arrangement of heterologous chromosome territories and homologous association of chromosomes 1, 3 and 5 (on average in 35–50% of nuclei) are in accordance with the random frequency predicted by computer simulations. Only the nucleolus organizing region (NOR)-bearing chromosome 2 and 4 homologs associate more often than randomly, since NORs mostly attach to a single nucleolus. Somatic pairing of homologous 100 kb segments occurs less frequently than homolog association, not significantly more often than expected at random and not simultaneously along the homologs. Thus, chromosome arrangement in Arabidopsis differs from that in Drosophila (characterized by somatic pairing of homologs), in spite of similar genome size, sequence organization and chromosome number. Nevertheless, in up to 31.5% of investigated Arabidopsis nuclei allelic sequences may share positions close enough for homologous recombination.Electronic Supplementary Material Supplementary material is available for this article at .     

From phenotype to genotype: a Bayesian solution

The study of biological systems commonly depends on inferring the state of a 'hidden' variable, such as an underlying genotype, from that of an 'observed' variable, such as an expressed phenotype. However, this cannot be achieved using traditional quantitative methods when more than one genetic mechanism exists for a single observable phenotype. Using a novel latent class Bayesian model, it is possible to infer the prevalence of different genetic elements in a population given a sample of phenotypes. As an exemplar, data comprising phenotypic resistance to six antimicrobials obtained from passive surveillance of Salmonella Typhimurium DT104 are analysed to infer the prevalence of individual resistance genes, as well as the prevalence of a genomic island known as SGI1 and its variants. Three competing models are fitted to the data and distinguished between using posterior predictive p-values to assess their ability to predict the observed number of unique phenotypes. The results suggest that several SGI1 variants circulate in a few fixed forms through the population from which our data were derived. The methods presented could be applied to other types of phenotypic data, and represent a useful and generic mechanism of inferring the genetic population structure of organisms.     

Captive breeding for the conservation of cichlid fishes

Cichlid fishes are by far the largest familial group of endangered vertebrates, especially the haplochromines. This paper concerns the organization and management of captive breeding of haplochromine cichlids. The setting up of a small-scale laboratory programme for the conservation of endangered species is described in terms of funding, staffing, installation and livestock husbandry. Breeding is discussed in the context of the selection of broodstock, basic reproductive biology, rearing, disease, pathological disorders and the arrangements necessary for the transfer and documentation of progeny. There are already indications in Africa and elsewhere that the dramatic decline and demise of cichlid taxa in Lake Victoria will not be an isolated phenomenon. There is no prospect that the captive breeding of cichlids can alone resolve such large-scale problems in fisheries management and ecology, or prevent the loss of taxa in nature. Nevertheless, captive breeding provides conservation options which are otherwise limited or unavailable regarding the saving of individual 'heritage' species, restocking to the wild, fundamental laboratory research and, finally, public education on the grave issue of actual and prospective mass extinctions of cichlids and other rare fishes.     

Role of Ser216 in the mechanism of action of membrane-bound lytic transglycosylase B: further evidence for substrate-assisted catalysis

Lytic transglycosylases cleave the beta-(1-->4)-glycosidic bond in the bacterial cell wall heteropolymer peptidoglycan between the N-acetylmuramic acid (MurNAc) and N-acetylglucosamine (GlcNAc) residues with the concomitant formation of a 1,6-anhydromuramoyl residue. Based on sequence alignments, Ser216 in Pseudomonas aeruginosa membrane-bound lytic transglycosylase B (MltB) was targeted for replacement with alanine to delineate its role in the enzyme's mechanism of action. The specific activity of the Ser216-->Ala MltB derivative was less than 12% of that for the wild-type enzyme, while its substrate binding affinity remained virtually unaltered. These data are in agreement with a role of Ser216 in orienting the N-acetyl group on MurNAc at the -1 subsite of MltB for its participation in a substrate-assisted mechanism of action.     

Cryptosporidium genotypes in wildlife from a new york watershed

To identify the animal sources for Cryptosporidium contamination, we genotyped Cryptosporidium spp. in wildlife from the watershed of the New York City drinking water supply, using a small-subunit rRNA gene-based PCR-restriction fragment length polymorphism analysis and DNA sequencing. A total of 541 specimens from 38 species of wildlife were analyzed. One hundred and eleven (20.5%) of the wildlife specimens were PCR positive. Altogether, 21 Cryptosporidium genotypes were found in wildlife samples, 11 of which were previously found in storm runoff in the watershed, and six of these 11 were from storm water genotypes of unknown animal origin. Four new genotypes were found, and the animal hosts for four storm water genotypes were expanded. With the exception of the cervine genotype, most genotypes were found in a limited number of animal species and have no major public health significance.     

Genetic diversity of Agrobacterium vitis as determined by DNA fingerprints of the 5'-end of the 23S rRNA gene and Random Amplified Polymorphic DNA

Sixty-nine strains of Agrobacterium vitis , the causal agent of grape crown gall, originating from different geographical regions of the USA and Europe, were characterized by fingerprint analysis of the 5'-end of the 23S rRNA gene and by Random Amplified Polymorphic DNA (RAPD) markers. For 5'-end 23S fingerprinting, amplicons were digested with Taq I, Rsa I, Ava I, Cfo I and Alu I. For RAPD analysis, three 10-mer primers were used to generate PCR products. There was a high degree of correlation between strain groupings generated by the two methods. However, more diversity was identified when groupings were based on RAPDs. For example, 28 of 29 strains having nopaline type Ti plasmids generated identical 5'-end 23S patterns but formed two distinct RAPD groups that separated strains originating from the USA and Hungary. Similarly by RAPDs, one cluster of strains carrying vitopine-type Ti plasmids could be separated into those originating in the USA and Europe. The composition of strain groups generated by 5'-end 23S and RAPDs were highly correlated with a previous fingerprint analysis of the intergenic spacer region (located between the 16S and 23S rRNA genes) and with RFLP analysis for characterizing Ti plasmids. These findings show that among Ag. vitis strains there is a high level of correlation between two regions of the rRNA operon, total genomic DNA (as determined by RAPDs) and the type of Ti plasmid they carry.