Decreasing the number of false positives in sequence classification

Abstract Background A large number of probabilistic models used in sequence analysis assign non-zero probability values to most input sequences. To decide when a given probability is sufficient the most common way is bayesian binary classification, where the probability of the model characterizing the sequence family of interest is compared to that of an alternative probability model. We can use as alternative model a null model. This is the scoring technique used by sequence analysis tools such as HMMER, SAM and INFERNAL. The most prevalent null models are position-independent residue distributions that include: the uniform distribution, genomic distribution, family-specific distribution and the target sequence distribution. This paper presents a study to evaluate the impact of the choice of a null model in the final result of classifications. In particular, we are interested in minimizing the number of false predictions in a classification. This is a crucial issue to reduce costs of biological validation. Results For all the tests, the target null model presented the lowest number of false positives, when using random sequences as a test. The study was performed in DNA sequences using GC content as the measure of content bias, but the results should be valid also for protein sequences. To broaden the application of the results, the study was performed using randomly generated sequences. Previous studies were performed on aminoacid sequences, using only one probabilistic model (HMM) and on a specific benchmark, and lack more general conclusions about the performance of null models. Finally, a benchmark test with P. falciparum confirmed these results. Conclusions Of the evaluated models the best suited for classification are the uniform model and the target model. However, the use of the uniform model presents a GC bias that can cause more false positives for candidate sequences with extreme compositional bias, a characteristic not described in previous studies. In these cases the target model is more dependable for biological validation due to its higher specificity.

Kinematic distances of pre-main-sequence stars in the Lupus star-forming region

To reliably determine the main physical parameters (masses and ages) of young stars, we must know their distances. While the average distance to nearby star-forming regions (<300 pc) is often known, the distances to individual stars are usually unknown. Individual distances to members of young moving groups can be derived from their radial velocities and proper motions using the convergent-point strategy. We investigate the kinematic properties of the Lupus moving group with the primary objective of deriving individual distances to all group members.

Genome characterization through a mathematical model of the genetic code: an analysis of the whole chromosome 1 of A. thaliana

The objective of this work is to characterize the genome of the chromosome 1 of A.thaliana, a small flowering plants used as a model organism in studies of biology and genetics, on the basis of a recent mathematical model of the genetic code. I analyze and compare different portions of the genome: genes, exons, coding sequences (CDS), introns, long introns, intergenes, untranslated regions (UTR) and regulatory sequences. In order to accomplish the task, I transformed nucleotide sequences into binary sequences based on the definition of the three different dichotomic classes. The descriptive analysis of binary strings indicate the presence of regularities in each portion of the genome considered. In particular, there are remarkable differences between coding sequences (CDS and exons) and non-coding sequences, suggesting that the frame is important only for coding sequences and that dichotomic classes can be useful to recognize them. Then, I assessed the existence of short-range dependence between binary sequences computed on the basis of the different dichotomic classes. I used three different measures of dependence: the well-known chi-squared test and two indices derived from the concept of entropy i.e. Mutual Information (MI) and Sρ, a normalized version of the “Bhattacharya Hellinger Matusita distance”. The results show that there is a significant short-range dependence structure only for the coding sequences whose existence is a clue of an underlying error detection and correction mechanism. No doubt, further studies are needed in order to assess how the information carried by dichotomic classes could discriminate between coding and noncoding sequence and, therefore, contribute to unveil the role of the mathematical structure in error detection and correction mechanisms. Still, I have shown the potential of the approach presented for understanding the management of genetic information.

Real-time polymerase chain-reaction detection of pathogens is feasible to supplement the diagnostic sequence for urinary tract infections

To evaluate, in a prospective pilot study, the feasibility of identifying pathogens in urine using real-time polymerase chain reaction (PCR), and to compare the results with the conventional urine culture-based procedures.

Multimodal approach in the use of clinical scoring, morphological MRI and biochemical T2-mapping and diffusion-weighted imaging in their ability to assess differences between cartilage repair tissue after microfracture therapy and matrix-associated autologous chondrocyte transplantation: a pilot study

OBJECTIVE: The aim of the present pilot study is to show initial results of a multimodal approach using clinical scoring, morphological magnetic resonance imaging (MRI) and biochemical T2-relaxation and diffusion-weighted imaging (DWI) in their ability to assess differences between cartilage repair tissue after microfracture therapy (MFX) and matrix-associated autologous chondrocyte transplantation (MACT). METHOD: Twenty patients were cross-sectionally evaluated at different post-operative intervals from 12 to 63 months after MFX and 12-59 months after MACT. The two groups were matched by age (MFX: 36.0+/-10.4 years; MACT: 35.1+/-7.7 years) and post-operative interval (MFX: 32.6+/-16.7 months; MACT: 31.7+/-18.3 months). After clinical evaluation using the Lysholm score, 3T-MRI was performed obtaining the MR observation of cartilage repair tissue (MOCART) score as well as T2-mapping and DWI for multi-parametric MRI. Quantitative T2-relaxation was achieved using a multi-echo spin-echo sequence; semi-quantitative diffusion-quotient (signal intensity without diffusion-weighting divided by signal intensity with diffusion weighting) was prepared by a partially balanced, steady-state gradient-echo pulse sequence. RESULTS: No differences in Lysholm (P=0.420) or MOCART (P=0.209) score were observed between MFX and MACT. T2-mapping showed lower T2 values after MFX compared to MACT (P=0.039). DWI distinguished between healthy cartilage and cartilage repair tissue in both procedures (MFX: P=0.001; MACT: P=0.007). Correlations were found between the Lysholm and the MOCART score (Pearson: 0.484; P=0.031), between the Lysholm score and DWI (Pearson:-0.557; P=0.011) and a trend between the Lysholm score and T2 (Person: 0.304; P=0.193). CONCLUSION: Using T2-mapping and DWI, additional information could be gained compared to clinical scoring or morphological MRI. In combination clinical, MR-morphological and MR-biochemical parameters can be seen as a promising multimodal tool in the follow-up of cartilage repair.

Improvement of cardiac function with device-based diaphragmatic stimulation in chronic heart failure patients: the randomized, open-label, crossover Epiphrenic II Pilot Trial.

AIMS Device-based pacing-induced diaphragmatic stimulation (PIDS) may have therapeutic potential for chronic heart failure (HF) patients. We studied the effects of PIDS on cardiac function and functional outcomes. METHODS AND RESULTS In 24 chronic HF patients with CRT, an additional electrode was attached to the left diaphragm. Randomized into two groups, patients received the following PIDS modes for 3 weeks in a different sequence: (i) PIDS off (control group); (ii) PIDS 0 ms mode (PIDS simultaneously with ventricular CRT pulse); or (iii) PIDS optimized mode (PIDS with optimized delay to ventricular CRT pulse). For PIDS optimization, acoustic cardiography was used. Effects of each PIDS mode on dyspnoea, power during exercise testing, and LVEF were assessed. Dyspnoea improved with the PIDS 0 ms mode (P = 0.057) and the PIDS optimized mode (P = 0.034) as compared with the control group. Maximal power increased from median 100.5 W in the control group to 104.0 W in the PIDS 0 ms mode (P = 0.092) and 109.5 W in the PIDS optimized mode (P = 0.022). Median LVEF was 33.5% in the control group, 33.0% in the PIDS 0 ms mode, and 37.0% in the PIDS optimized mode (P = 0.763 and P = 0.009 as compared with the control group, respectively). PIDS was asymptomatic in all patients. CONCLUSION PIDS improves dyspnoea, working capacity, and LVEF in chronic HF patients over a 3 week period in addition to CRT. This pilot study demonstrates proof of principle of an innovative technology which should be confirmed in a larger sample. TRIAL REGISTRATION NCT00769678.

The sequence comparison index: A novel method for comparing proteins and proteomes

Historically morphological features were used as the primary means to classify organisms. However, the age of molecular genetics has allowed us to approach this field from the perspective of the organism's genetic code. Early work used highly conserved sequences, such as ribosomal RNA. The increasing number of complete genomes in the public data repositories provides the opportunity to look not only at a single gene, but at organisms' entire parts list. ^ Here the Sequence Comparison Index (SCI) and the Organism Comparison Index (OCI), algorithms and methods to compare proteins and proteomes, are presented. The complete proteomes of 104 sequenced organisms were compared. Over 280 million full Smith-Waterman alignments were performed on sequence pairs which had a reasonable expectation of being related. From these alignments a whole proteome phylogenetic tree was constructed. This method was also used to compare the small subunit (SSU) rRNA from each organism and a tree constructed from these results. The SSU rRNA tree by the SCI/OCI method looks very much like accepted SSU rRNA trees from sources such as the Ribosomal Database Project, thus validating the method. The SCI/OCI proteome tree showed a number of small but significant differences when compared to the SSU rRNA tree and proteome trees constructed by other methods. Horizontal gene transfer does not appear to affect the SCI/OCI trees until the transferred genes make up a large portion of the proteome. ^ As part of this work, the Database of Related Local Alignments (DaRLA) was created and contains over 81 million rows of sequence alignment information. DaRLA, while primarily used to build the whole proteome trees, can also be applied shared gene content analysis, gene order analysis, and creating individual protein trees. ^ Finally, the standard BLAST method for analyzing shared gene content was compared to the SCI method using 4 spirochetes. The SCI system performed flawlessly, finding all proteins from one organism against itself and finding all the ribosomal proteins between organisms. The BLAST system missed some proteins from its respective organism and failed to detect small ribosomal proteins between organisms. ^

Searching for the ideal DNA recognition code of covalent ligands

The combitiatorial approach restriction endonuclease protection selection and amplification REPSA was successfully used to determine ideal DNA interactions sites of covalent ligands. Unlike most other combinatorial methods, REPSA is based on inhibition of enzymatic cleavage by specific ligand-DNA complexes, which enables identification of binding sites of various ligands. However, the inherent nature of this technique posses a problem during selection of binding sites of covalent ligands. By modifying the technique according to the nature of the ligand, we demonstrate the flexibility of REPSA in identifying the preferred binding sites for monocovalent ligands, topoisomerase I and tallimustine, and the bicovalent ligand topoisomerase II. From among the preferred binding sites, we identified the consensus binding sequence of camptothecin induced topoisomerase I cleavage as ‘aGWT/Gc’, and tallimustine consensus sequences as ‘GTTCTA’ and ‘TTTTTTC’. We have shown for the first time that preferential binding of tallimustine occurs at sequences not previously reported. Furthermore, our data indicate that tallimustine is a novel DNA minor groove, guanine-specific alkylating agent. ^ Additionally, we have demonstrated in vivo that sequence-specific covalent DNA-binding small molecules have the ability to regulate transcription by inhibiting RNA polymerase II. Tallimustine, binding to its preferred sequences located in the 5′ untranslated region were an effective impediment for transcribing polymerase II. The ability of covalent binding small molecules to target predetermined DNA sequences located downstream of the promoter suggests a general approach for regulation of gene expression. ^

Secondary mineral assemblages in a volcanic sequence of ODP Hole 104-642E

Mineralogical and geochemical analyses of alteration products from upper and lower volcanic series recovered during ODP Leg 104 reveal variations both in composition and order of crystallization of clay minerals vesicles and voids filling and replacing glass. These results provide information about successive alteration stages of rocks and interlayered volcaniclastic sediments. The first stage, related to initial basalt-seawater interaction, is characterized by development of Fe-smectites, especially Fe-rich saponite. A second stage of intermittently superimposed subaerial weathering is marked by iron-oxides-halloysite-kaolinite formation. The third episode, interpreted as hydrothermal on the basis of O-isotopic data, is defined by postburial coprecipitation of Fe-poor, Mg-rich saponite and celadonite. A distinct final and pervasive hydrothermal stage, occurring mainly in the lower series and dominated by Al-smectites-zeolites assemblage, indicates changes toward a more reducing alteration environment.

Stable isotope, seismic sequence boundaries and bioevents in ODP Hole 166-1006A

During ODP Leg 166, the recovery of cores from a transect of drill sites across the Bahamas margin from marginal to deep basin environments was an essential requirement for the study of the response of the sedimentary systems to sea-level changes. A detailed biostratigraphy based on planktonic foraminifera was performed on ODP Hole 1006A for an accurate stratigraphic control. The investigated late middle Miocene-early Pliocene sequence spans the interval from about 12.5 Ma (Biozone N12) to approximately 4.5 Ma (Biozone N19). Several bioevents calibrated with the time scale of Berggren et al. (1995a,b) were identified. The ODP Site 1006 benthic oxygen isotope stratigraphy can be correlated to the corresponding deep-water benthic oxygen isotope curve from ODP Site 846 in the Eastern Equatorial Pacific (Shackleton et al., 1995. Proc. ODP Sci. Res. 138, 337-356), which was orbitally tuned for the entire Pliocene into the latest Miocene at 6.0 Ma. The approximate stratigraphic match of the isotopic signals from both records between 4.5 and 6.0 Ma implies that the paleoceanographic signal from the Bahamas is not simply a record of regional variations but, indeed, represents glacio-eustatic fluctuations. The ODP Site 1006 oxygen and carbon isotope record, based on benthic and planktonic foraminifera, was used to define paleoceanographic changes on the margin, which could be tied to lithostratigraphic events on the Bahamas carbonate platform using seismic sequence stratigraphy. The oxygen isotope values show a general cooling trend from the middle to late Miocene, which was interrupted by a significant trend towards warmer sea-surface temperatures (SST) and associated sea-level rise with decreased ice volume during the latest Miocene. This trend reached a maximum coincident with the Miocene/Pliocene boundary. An abrupt cooling in the early Pliocene then followed the warming which continued into the earliest Pliocene. The late Miocene paleoceanographic evolution along the Bahamas margin can be observed in the ODP Site 1006 delta13C values, which support other evidence for the beginning of the closure of the Panama gateway at 8 Ma followed by a reduced intermediate water supply of water from the Pacific into the Caribbean at about 5 Ma. A general correlation of lower sedimentation rates with the major seismic sequence boundaries (SSBs) was observed. Additionally, the SSBs are associated with transitions towards more positive oxygen isotope excursions. This observed correspondence implies that the presence of a SSB, representing a density impedance contrast in the sedimentary sequence, may reflect changes in the character of the deposited sediment during highstands versus those during lowstands. However, not all of the recorded oxygen isotope excursions correspond to SSBs. The absence of a SSB in association with an oxygen isotope excursion indicates that not all oxygen isotope sea-level events impact the carbonate margin to the same extent, or maybe even represent equivalent sea-level fluctuations. Thus, it can be tentatively concluded that SSBs produced on carbonate margins do record sea-level fluctuations but not every sea-level fluctuation is represented by a SSB in the sequence stratigraphic record.

Late Pliocene dinoflagellate cyst and diatom analysis from a high resolution sequence in DSDP Site 594, Chatham Rise, south west Pacific

A high resolution mixed carbonate and siliciclastic sequence from DSDP Site 594 contains a detailed record of climate change in the late Pliocene. The sequence can be accurately dated by the LAD of Nitzschia weaveri, the LAD of Thalassiosira insigna, the LAD of T. vulnifica and the LAD of T. kolbei diatom datums. Carbonate content and delta18O signatures provide added resolution and place the sequence between isotope stage 100 and 92. The sequence contains well-preserved and diverse dinoflagellate cyst floras. Use of principal component (PCA) and canonical correspondence analyses (CCA) identifies changes in the assemblages that principally reflect warming and cooling trends. Species association with warmer climates included Impagidinium patulum, I. paradoxum and I. sp. cf. paradoxum while those from cooler climates include Invertecysta tabulata and I. velorum. CCA is shown to be a valuable method of determining the past environmental preferences of extinct species such as I. tabulata.