Approaching the taxonomic affiliation of unidentified sequences in public databases – an example from the mycorrhizal fungi

Background  

During the last few years, DNA sequence analysis has become one of the primary means of taxonomic identification of species, particularly so for species that are minute or otherwise lack distinct, readily obtainable morphological characters. Although the number of sequences available for comparison in public databases such as GenBank increases exponentially, only a minuscule fraction of all organisms have been sequenced, leaving taxon sampling a momentous problem for sequence-based taxonomic identification. When querying GenBank with a set of unidentified sequences, a considerable proportion typically lack fully identified matches, forming an ever-mounting pile of sequences that the researcher will have to monitor manually in the hope that new, clarifying sequences have been submitted by other researchers. To alleviate these concerns, a project to automatically monitor select unidentified sequences in GenBank for taxonomic progress through repeated local BLAST searches was initiated. Mycorrhizal fungi – a field where species identification often is prohibitively complex – and the much usedITSlocus were chosen as test bed.     

Distinguishing importation from diversification of quinolone-resistant <Emphasis Type="Italic">Neisseria gonorrhoeae</Emphasis> by molecular evolutionary analysis

Background  

Distinguishing the recent introduction of quinolone resistant gonococci into a population from diversification of resistant strains already in the population is important for planning effective infection control strategies. We applied molecular evolutionary analyses to DNA sequences from 9 housekeeping genes and gyrA, parC and porB of 24 quinolone resistant N. gonorrhoeae (QRNG) and 24 quinolone sensitive isolates collected in Israel during 2000–2001.     

The high-resolution NMR structure of the R21A Spc-SH3:P41 complex: Understanding the determinants of binding affinity by comparison with Abl-SH3

Background  

SH3 domains are small protein modules of 60–85 amino acids that bind to short proline-rich sequences with moderate-to-low affinity and specificity. Interactions with SH3 domains play a crucial role in regulation of many cellular processes (some are related to cancer and AIDS) and have thus been interesting targets in drug design. The decapeptide APSYSPPPPP (p41) binds with relatively high affinity to the SH3 domain of the Abl tyrosine kinase (Abl-SH3), while it has a 100 times lower affinity for the α-spectrin SH3 domain (Spc-SH3).     

What do we know about communicating risk? A brief review and suggestion for contextualising serious,but rare,risk, and the example of cox-2 selective and non-selective NSAIDs

Background  

Communicating risk is difficult. Although different methods have been proposed – using numbers, words, pictures or combinations – none has been extensively tested. We used electronic and bibliographic searches to review evidence concerning risk perception and presentation. People tend to underestimate common risk and overestimate rare risk; they respond to risks primarily on the basis of emotion rather than facts, seem to be risk averse when faced with medical interventions, and want information on even the rarest of adverse events.     

Genomic multiple sequence alignments: refinement using a genetic algorithm

Background  

Genomic sequence data cannot be fully appreciated in isolation. Comparative genomics – the practice of comparing genomic sequences from different species – plays an increasingly important role in understanding the genotypic differences between species that result in phenotypic differences as well as in revealing patterns of evolutionary relationships. One of the major challenges in comparative genomics is producing a high-quality alignment between two or more related genomic sequences. In recent years, a number of tools have been developed for aligning large genomic sequences. Most utilize heuristic strategies to identify a series of strong sequence similarities, which are then used as anchors to align the regions between the anchor points. The resulting alignment is globally correct, but in many cases is suboptimal locally. We describe a new program, GenAlignRefine, which improves the overall quality of global multiple alignments by using a genetic algorithm to improve local regions of alignment. Regions of low quality are identified, realigned using the program T-Coffee, and then refined using a genetic algorithm. Because a better COFFEE (Consistency based Objective Function For alignmEnt Evaluation) score generally reflects greater alignment quality, the algorithm searches for an alignment that yields a better COFFEE score. To improve the intrinsic slowness of the genetic algorithm, GenAlignRefine was implemented as a parallel, cluster-based program.     

Detecting overlapping coding sequences in virus genomes

Background  

Detecting new coding sequences (CDSs) in viral genomes can be difficult for several reasons. The typically compact genomes often contain a number of overlapping coding and non-coding functional elements, which can result in unusual patterns of codon usage; conservation between related sequences can be difficult to interpret – especially within overlapping genes; and viruses often employ non-canonical translational mechanisms – e.g. frameshifting, stop codon read-through, leaky-scanning and internal ribosome entry sites – which can conceal potentially coding open reading frames (ORFs).     

An EST database from saffron stigmas

Background  

Saffron (Crocus sativus L., Iridaceae) flowers have been used as a spice and medicinal plant ever since the Greek-Minoan civilization. The edible part – the stigmas – are commonly considered the most expensive spice in the world and are the site of a peculiar secondary metabolism, responsible for the characteristic color and flavor of saffron.     

An enhanced RNA alignment benchmark for sequence alignment programs

Background  

The performance of alignment programs is traditionally tested on sets of protein sequences, of which a reference alignment is known. Conclusions drawn from such protein benchmarks do not necessarily hold for the RNA alignment problem, as was demonstrated in the first RNA alignment benchmark published so far. For example, the twilight zone – the similarity range where alignment quality drops drastically – starts at 60 % for RNAs in comparison to 20 % for proteins. In this study we enhance the previous benchmark.     

The VirusBanker database uses a Java program to allow flexible searching through <Emphasis Type="Italic">Bunyaviridae</Emphasis> sequences

Background  

Viruses of the Bunyaviridae have segmented negative-stranded RNA genomes and several of them cause significant disease. Many partial sequences have been obtained from the segments so that GenBank searches give complex results. Sequence databases usually use HTML pages to mediate remote sorting, but this approach can be limiting and may discourage a user from exploring a database.     

Sleep in Kcna2 knockout mice

Background  

Shaker codes for a Drosophila voltage-dependent potassium channel. Flies carrying Shaker null or hypomorphic mutations sleep 3–4 h/day instead of 8–14 h/day as their wild-type siblings do. Shaker-like channels are conserved across species but it is unknown whether they affect sleep in mammals. To address this issue, we studied sleep in Kcna2 knockout (KO) mice. Kcna2 codes for Kv1.2, the alpha subunit of a Shaker-like voltage-dependent potassium channel with high expression in the mammalian thalamocortical system.     

RIPCAL: a tool for alignment-based analysis of repeat-induced point mutations in fungal genomic sequences

Background  

Repeat-induced point mutation (RIP) is a fungal-specific genome defence mechanism that alters the sequences of repetitive DNA, thereby inactivating coding genes. Repeated DNA sequences align between mating and meiosis and both sequences undergo C:G to T:A transitions. In most fungi these transitions preferentially affect CpA di-nucleotides thus altering the frequency of certain di-nucleotides in the affected sequences. The majority of previously published in silico analyses were limited to the comparison of ratios of pre- and post-RIP di-nucleotides in putatively RIP-affected sequences – so-called RIP indices. The analysis of RIP is significantly more informative when comparing sequence alignments of repeated sequences. There is, however, a dearth of bioinformatics tools available to the fungal research community for alignment-based RIP analysis of repeat families.     

<Emphasis Type="Italic">Geobacillus zalihae</Emphasis> sp. nov., a thermophilic lipolytic bacterium isolated from palm oil mill effluent in Malaysia

Background  

Thermophilic Bacillus strains of phylogenetic Bacillus rRNA group 5 were described as a new genus Geobacillus. Their geographical distribution included oilfields, hay compost, hydrothermal vent or soils. The members from the genus Geobacillus have a growth temperatures ranging from 35 to 78°C and contained iso-branched saturated fatty acids (iso-15:0, iso-16:0 and iso-17:0) as the major fatty acids. The members of Geobacillus have similarity in their 16S rRNA gene sequences (96.5–99.2%). Thermophiles harboring intrinsically stable enzymes are suitable for industrial applications. The quest for intrinsically thermostable lipases from thermophiles is a prominent task due to the laborious processes via genetic modification.     

Characterization of expressed class II MHC sequences in the banner-tailed kangaroo rat (<Emphasis Type="Italic">Dipodomys spectabilis</Emphasis>) reveals multiple <Emphasis Type="Italic">DRB</Emphasis> loci

Genes of the major histocompatibility complex (MHC) are exceptionally polymorphic due to the combined effects of natural and sexual selection. Most research in wild populations has focused on the second exon of a single class II locus (DRB), but complete gene sequences can provide an illuminating backdrop for studies of intragenic selection, recombination, and organization. To this end, we characterized class II loci in the banner-tailed kangaroo rat (Dipodomys spectabilis). Seven DRB-like sequences (provisionally named MhcDisp-DRB*01 through *07) were isolated from spleen cDNA and most likely comprise ≥5 loci; this multiformity is quite unlike the situation in muroid rodents such as Mus, Rattus, and Peromyscus. In silico translation revealed the presence of important structural residues for glycosylation sites, salt bonds, and CD4+ T-cell recognition. Amino-acid distances varied widely among the seven sequences (2–34%). Nuclear DNA sequences from the Disp-DRB*07 locus (∼10 kb) revealed a conventional exon/intron structure as well as a number of microsatellites and short interspersed nuclear elements (B4, Alu, and IDL-Geo subfamilies). Rates of nucleotide substitution at Disp-DRB*07 are similar in both exons and introns (π = 0.015 and 0.012, respectively), which suggests relaxed selection and may indicate that this locus is an expressed pseudogene. Finally, we performed BLASTn searches against Dipodomys ordii genomic sequences (unassembled reads) and find 90–97% nucleotide similarity between the two kangaroo rat species. Collectively, these data suggest that class II diversity in heteromyid rodents is based on polylocism and departs from the muroid architecture. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the accession numbers EU817477–EU817485.     

Role of subtyping in detecting Salmonella cross contamination in the laboratory

Background  

With the exception of M. tuberculosis, little has been published on the problems of cross-contamination in bacteriology laboratories. We performed a retrospective analysis of subtyping data from the National Salmonella Reference Laboratory (Ireland) from 2000–2007 to identify likely incidents of laboratory cross contamination.     

TRUNCATULIX – a data warehouse for the legume community

Background  

Databases for either sequence, annotation, or microarray experiments data are extremely beneficial to the research community, as they centrally gather information from experiments performed by different scientists. However, data from different sources develop their full capacities only when combined. The idea of a data warehouse directly adresses this problem and solves it by integrating all required data into one single database – hence there are already many data warehouses available to genetics. For the model legume Medicago truncatula, there is currently no such single data warehouse that integrates all freely available gene sequences, the corresponding gene expression data, and annotation information. Thus, we created the data warehouse TRUNCATULIX, an integrative database of Medicago truncatula sequence and expression data.     

Climate change and recent genetic flux in populations of <Emphasis Type="Italic">Drosophila robusta</Emphasis>

Background  

Studied since the early 1940's, chromosomal polymorphisms in the deciduous woods species Drosophila robusta have been characterized by well-defined latitudinal, longitudinal, and elevational clines, but – until at least ten years ago – stable, local population frequencies. Recent biogeographical analyses indicate that D. robusta invaded North America from southeast Asia and has persisted in eastern temperate forests for at least 20–25 my without speciating. The abundant chromosome polymorphisms found across the range of D. robusta are thus likely to be relatively ancient, having accumulated over many well known climatic cycles in North America. Sufficient long-term data are now available such that we can now gauge the rate of these evolutionary changes in natural populations due to environmental change.     

The Subviral RNA Database: a toolbox for viroids, the hepatitis delta virus and satellite RNAs research

Background  

Viroids, satellite RNAs, satellites viruses and the human hepatitis delta virus form the 'brotherhood' of the smallest known infectious RNA agents, known as the subviral RNAs. For most of these species, it is generally accepted that characteristics such as cell movement, replication, host specificity and pathogenicity are encoded in their RNA sequences and their resulting RNA structures. Although many sequences are indexed in publicly available databases, these sequence annotation databases do not provide the advanced searches and data manipulation capability for identifying and characterizing subviral RNA motifs.     

Prediction of beta-turns at over 80% accuracy based on an ensemble of predicted secondary structures and multiple alignments

Background  

β-turn is a secondary protein structure type that plays significant role in protein folding, stability, and molecular recognition. To date, several methods for prediction of β-turns from protein sequences were developed, but they are characterized by relatively poor prediction quality. The novelty of the proposed sequence-based β-turn predictor stems from the usage of a window based information extracted from four predicted three-state secondary structures, which together with a selected set of position specific scoring matrix (PSSM) values serve as an input to the support vector machine (SVM) predictor.