The database search problem: a question of rational decision making

This paper applies probability and decision theory in the graphical interface of an influence diagram to study the formal requirements of rationality which justify the individualization of a person found through a database search. The decision-theoretic part of the analysis studies the parameters that a rational decision maker would use to individualize the selected person. The modeling part (in the form of an influence diagram) clarifies the relationships between this decision and the ingredients that make up the database search problem, i.e., the results of the database search and the different pairs of propositions describing whether an individual is at the source of the crime stain. These analyses evaluate the desirability associated with the decision of 'individualizing' (and 'not individualizing'). They point out that this decision is a function of (i) the probability that the individual in question is, in fact, at the source of the crime stain (i.e., the state of nature), and (ii) the decision maker's preferences among the possible consequences of the decision (i.e., the decision maker's loss function). We discuss the relevance and argumentative implications of these insights with respect to recent comments in specialized literature, which suggest points of view that are opposed to the results of our study.

Fourmidable: a database for ant genomics.

BACKGROUND: Fourmidable is an infrastructure to curate and share the emerging genetic, molecular, and functional genomic data and protocols for ants. DESCRIPTION: The Fourmidable assembly pipeline groups nucleotide sequences into clusters before independently assembling each cluster. Subsequently, assembled sequences are annotated via Interproscan and BLAST against general and insect-specific databases. Gene-specific information can be retrieved using gene identifiers, searching for similar sequences or browsing through inferred Gene Ontology annotations. The database will readily scale as ultra-high throughput sequence data and sequences from additional species become available. CONCLUSION: Fourmidable currently houses EST data from two ant species and microarray gene expression data for one of these. Fourmidable is publicly available at http://fourmidable.unil.ch.

Visualization of RNA polymerase II ternary transcription complexes formed in vitro on a Xenopus laevis vitellogenin gene.

Stable ternary transcription complexes assembled in vitro, using a HeLa whole-cell extract, have been isolated and visualized by electron microscopy. The formation of these stable complexes on the DNA fragment used as template, the 5' end region of the Xenopus laevis vitellogenin gene B2, depends on factors present in the whole-cell extract, RNA polymerase II and at least two nucleotides. Interestingly, bending in the DNA fragment was frequently observed at the binding site of RNA polymerase II. Dinucleotides that can prime initiation within a short sequence of approximately 10 contiguous nucleotides centered around the initiation site used in vivo, also favour the formation of stable complexes. In addition, pre-initiation complexes were isolated and it was shown that factors in the extract involved in their formation are more abundant than the RNA polymerase II molecules available for binding. The possible implication of this observation relative to the in vivo situation is discussed.

CATMA: a complete Arabidopsis GST database.

The Complete Arabidopsis Transcriptome Micro Array (CATMA) database contains gene sequence tag (GST) and gene model sequences for over 70% of the predicted genes in the Arabidopsis thaliana genome as well as primer sequences for GST amplification and a wide range of supplementary information. All CATMA GST sequences are specific to the gene for which they were designed, and all gene models were predicted from a complete reannotation of the genome using uniform parameters. The database is searchable by sequence name, sequence homology or direct SQL query, and is available through the CATMA website at http://www.catma.org/.

RNA/DNA co-analysis from blood stains--results of a second collaborative EDNAP exercise.

A second collaborative exercise on RNA/DNA co-analysis for body fluid identification and STR profiling was organized by the European DNA Profiling Group (EDNAP). Six human blood stains, two blood dilution series (5-0.001 μl blood) and, optionally, bona fide or mock casework samples of human or non-human origin were analyzed by the participating laboratories using a RNA/DNA co-extraction or solely RNA extraction method. Two novel mRNA multiplexes were used for the identification of blood: a highly sensitive duplex (HBA, HBB) and a moderately sensitive pentaplex (ALAS2, CD3G, ANK1, SPTB and PBGD). The laboratories used different chemistries and instrumentation. All of the 18 participating laboratories were able to successfully isolate and detect mRNA in dried blood stains. Thirteen laboratories simultaneously extracted RNA and DNA from individual stains and were able to utilize mRNA profiling to confirm the presence of blood and to obtain autosomal STR profiles from the blood stain donors. The positive identification of blood and good quality DNA profiles were also obtained from old and compromised casework samples. The method proved to be reproducible and sensitive using different analysis strategies. The results of this collaborative exercise involving a RNA/DNA co-extraction strategy support the potential use of an mRNA based system for the identification of blood in forensic casework that is compatible with current DNA analysis methodology.

Detection of live and antibiotic-killed bacteria by quantitative real-time PCR of specific fragments of rRNA.

Assessing bacterial viability by molecular markers might help accelerate the measurement of antibiotic-induced killing. This study investigated whether rRNA could be suitable for this purpose. Cultures of penicillin-susceptible and penicillin-tolerant (Tol1 mutant) Streptococcus gordonii were exposed to mechanistically different penicillin and levofloxacin. Bacterial survival was assessed by viable counts and compared to quantitative real-time PCR amplification of either the 16S rRNA genes or the 16S rRNA, following reverse transcription. Penicillin-susceptible S. gordonii lost > or =4 log(10) CFU/ml of viability over 48 h of penicillin treatment. In comparison, the Tol1 mutant lost < or =1 log(10) CFU/ml. Amplification of a 427-bp fragment of 16S rRNA genes yielded amplicons that increased proportionally to viable counts during bacterial growth but did not decrease during drug-induced killing. In contrast, the same 427-bp fragment amplified from 16S rRNA paralleled both bacterial growth and drug-induced killing. It also differentiated between penicillin-induced killing of the parent and the Tol1 mutant (> or =4 log(10) CFU/ml and < or =1 log(10) CFU/ml, respectively) and detected killing by mechanistically unrelated levofloxacin. Since large fragments of polynucleotides might be degraded faster than smaller fragments, the experiments were repeated by amplifying a 119-bp region internal to the original 427-bp fragment. The amount of 119-bp amplicons increased proportionally to viability during growth but remained stable during drug treatment. Thus, 16S rRNA was a marker of antibiotic-induced killing, but the size of the amplified fragment was critical for differentiation between live and dead bacteria.

Prevalence of Neural Tube Defects (NTD) in Younger Mothers in Europe 2000-2008: Analysis of the EUROCAT Database - Report to Bayer Schiering Pharma

M.C. Addor is included in the Eurocat Working Group

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.

Permanent.Plot.ch - a database for Swiss permanent vegetation plots

The aim of the Permanent.Plot.ch project is the conservation of historical data about permanent plots in Switzerland and the monitoring of vegetation in a context of environmental changes (mainly climate and land use). Permanent plots are currently being recognized as valuable tools to monitor long-term effects of environmental changes on vegetation. Often used in short studies (3 to 5 years), they are generally abandoned at the end of projects. However, their full potential might only be revealed after 10 or more years, once the location is lost. For instance, some of the oldest permanent plots in Switzerland (first half of the 20th century) were nearly lost, although they are now very valuable data. The Permanent.Plot.ch national database (GIVD ID EU-CH-001), by storing historical and recent data, will allow to ensuring future access to data from permanent vegetation plots. As the database contains some private data, it is not directly available on internet but an overview of the data can be downloaded from internet (http://www.unil.ch/ppch) and precise data are available on request.

Genome-wide search reveals a novel GacA-regulated small RNA in Pseudomonas species.

BACKGROUND: Small RNAs (sRNAs) are widespread among bacteria and have diverse regulatory roles. Most of these sRNAs have been discovered by a combination of computational and experimental methods. In Pseudomonas aeruginosa, a ubiquitous Gram-negative bacterium and opportunistic human pathogen, the GacS/GacA two-component system positively controls the transcription of two sRNAs (RsmY, RsmZ), which are crucial for the expression of genes involved in virulence. In the biocontrol bacterium Pseudomonas fluorescens CHA0, three GacA-controlled sRNAs (RsmX, RsmY, RsmZ) regulate the response to oxidative stress and the expression of extracellular products including biocontrol factors. RsmX, RsmY and RsmZ contain multiple unpaired GGA motifs and control the expression of target mRNAs at the translational level, by sequestration of translational repressor proteins of the RsmA family. RESULTS: A combined computational and experimental approach enabled us to identify 14 intergenic regions encoding sRNAs in P. aeruginosa. Eight of these regions encode newly identified sRNAs. The intergenic region 1698 was found to specify a novel GacA-controlled sRNA termed RgsA. GacA regulation appeared to be indirect. In P. fluorescens CHA0, an RgsA homolog was also expressed under positive GacA control. This 120-nt sRNA contained a single GGA motif and, unlike RsmX, RsmY and RsmZ, was unable to derepress translation of the hcnA gene (involved in the biosynthesis of the biocontrol factor hydrogen cyanide), but contributed to the bacterium's resistance to hydrogen peroxide. In both P. aeruginosa and P. fluorescens the stress sigma factor RpoS was essential for RgsA expression. CONCLUSION: The discovery of an additional sRNA expressed under GacA control in two Pseudomonas species highlights the complexity of this global regulatory system and suggests that the mode of action of GacA control may be more elaborate than previously suspected. Our results also confirm that several GGA motifs are required in an sRNA for sequestration of the RsmA protein.

MyD88 and TLR9 dependent immune responses mediate resistance to Leishmania guyanensis infections, irrespective of Leishmania RNA virus burden.

Infections with Leishmania parasites of the Leishmania Viannia subgenus give rise to both localized cutaneous (CL), and metastatic leishmaniasis. Metastasizing disease forms including disseminated (DCL) and mutocutaneous (MCL) leishmaniasis result from parasitic dissemination and lesion formation at sites distal to infection and have increased inflammatory responses. The presence of Leishmania RNA virus (LRV) in L. guyanensis parasites contributes to the exacerbation of disease and impacts inflammatory responses via activation of TLR3 by the viral dsRNA. In this study we investigated other innate immune response adaptor protein modulators and demonstrated that both MyD88 and TLR9 played a crucial role in the development of Th1-dependent healing responses against L. guyanensis parasites regardless of their LRV status. The absence of MyD88- or TLR9-dependent signaling pathways resulted in increased Th2 associated cytokines (IL-4 and IL-13), which was correlated with low transcript levels of IL-12p40. The reliance of IL-12 was further confirmed in IL12AB-/- mice, which were completely susceptible to infection. Protection to L. guyanensis infection driven by MyD88- and TLR9-dependent immune responses arises independently to those induced due to high LRV burden within the parasites.

Recent misconceptions about the "database search problem": a probabilistic analysis using Bayesian networks

This paper analyses and discusses arguments that emerge from a recent discussion about the proper assessment of the evidential value of correspondences observed between the characteristics of a crime stain and those of a sample from a suspect when (i) this latter individual is found as a result of a database search and (ii) remaining database members are excluded as potential sources (because of different analytical characteristics). Using a graphical probability approach (i.e., Bayesian networks), the paper here intends to clarify that there is no need to (i) introduce a correction factor equal to the size of the searched database (i.e., to reduce a likelihood ratio), nor to (ii) adopt a propositional level not directly related to the suspect matching the crime stain (i.e., a proposition of the kind 'some person in (outside) the database is the source of the crime stain' rather than 'the suspect (some other person) is the source of the crime stain'). The present research thus confirms existing literature on the topic that has repeatedly demonstrated that the latter two requirements (i) and (ii) should not be a cause of concern.

Genomic Study of RNA Polymerase II and III SNAP(c)-Bound Promoters Reveals a Gene Transcribed by Both Enzymes and a Broad Use of Common Activators.

SNAP(c) is one of a few basal transcription factors used by both RNA polymerase (pol) II and pol III. To define the set of active SNAP(c)-dependent promoters in human cells, we have localized genome-wide four SNAP(c) subunits, GTF2B (TFIIB), BRF2, pol II, and pol III. Among some seventy loci occupied by SNAP(c) and other factors, including pol II snRNA genes, pol III genes with type 3 promoters, and a few un-annotated loci, most are primarily occupied by either pol II and GTF2B, or pol III and BRF2. A notable exception is the RPPH1 gene, which is occupied by significant amounts of both polymerases. We show that the large majority of SNAP(c)-dependent promoters recruit POU2F1 and/or ZNF143 on their enhancer region, and a subset also recruits GABP, a factor newly implicated in SNAP(c)-dependent transcription. These activators associate with pol II and III promoters in G1 slightly before the polymerase, and ZNF143 is required for efficient transcription initiation complex assembly. The results characterize a set of genes with unique properties and establish that polymerase specificity is not absolute in vivo.

Mynodbcsv: Lightweight Zero-Config Database Solution for Handling Very Large CSV Files.

Volumes of data used in science and industry are growing rapidly. When researchers face the challenge of analyzing them, their format is often the first obstacle. Lack of standardized ways of exploring different data layouts requires an effort each time to solve the problem from scratch. Possibility to access data in a rich, uniform manner, e.g. using Structured Query Language (SQL) would offer expressiveness and user-friendliness. Comma-separated values (CSV) are one of the most common data storage formats. Despite its simplicity, with growing file size handling it becomes non-trivial. Importing CSVs into existing databases is time-consuming and troublesome, or even impossible if its horizontal dimension reaches thousands of columns. Most databases are optimized for handling large number of rows rather than columns, therefore, performance for datasets with non-typical layouts is often unacceptable. Other challenges include schema creation, updates and repeated data imports. To address the above-mentioned problems, I present a system for accessing very large CSV-based datasets by means of SQL. It's characterized by: "no copy" approach - data stay mostly in the CSV files; "zero configuration" - no need to specify database schema; written in C++, with boost [1], SQLite [2] and Qt [3], doesn't require installation and has very small size; query rewriting, dynamic creation of indices for appropriate columns and static data retrieval directly from CSV files ensure efficient plan execution; effortless support for millions of columns; due to per-value typing, using mixed text/numbers data is easy; very simple network protocol provides efficient interface for MATLAB and reduces implementation time for other languages. The software is available as freeware along with educational videos on its website [4]. It doesn't need any prerequisites to run, as all of the libraries are included in the distribution package. I test it against existing database solutions using a battery of benchmarks and discuss the results.

Molecular basis of messenger RNA recognition by the specific bacterial repressing clamp RsmA/CsrA.

Proteins of the RsmA/CsrA family are global translational regulators in many bacterial species. We have determined the solution structure of a complex formed between the RsmE protein, a member of this family from Pseudomonas fluorescens, and a target RNA encompassing the ribosome-binding site of the hcnA gene. The RsmE homodimer with its two RNA-binding sites makes optimal contact with an 5'-A/UCANGGANGU/A-3' sequence in the mRNA. When tightly gripped by RsmE, the ANGGAN core folds into a loop, favoring the formation of a 3-base-pair stem by flanking nucleotides. We validated these findings by in vivo and in vitro mutational analyses. The structure of the complex explains well how, by sequestering the Shine-Dalgarno sequence, the RsmA/CsrA proteins repress translation.