Improved statistical analysis of low abundance phenomena in bimodal bacterial populations.

Accurate detection of subpopulation size determinations in bimodal populations remains problematic yet it represents a powerful way by which cellular heterogeneity under different environmental conditions can be compared. So far, most studies have relied on qualitative descriptions of population distribution patterns, on population-independent descriptors, or on arbitrary placement of thresholds distinguishing biological ON from OFF states. We found that all these methods fall short of accurately describing small population sizes in bimodal populations. Here we propose a simple, statistics-based method for the analysis of small subpopulation sizes for use in the free software environment R and test this method on real as well as simulated data. Four so-called population splitting methods were designed with different algorithms that can estimate subpopulation sizes from bimodal populations. All four methods proved more precise than previously used methods when analyzing subpopulation sizes of transfer competent cells arising in populations of the bacterium Pseudomonas knackmussii B13. The methods' resolving powers were further explored by bootstrapping and simulations. Two of the methods were not severely limited by the proportions of subpopulations they could estimate correctly, but the two others only allowed accurate subpopulation quantification when this amounted to less than 25% of the total population. In contrast, only one method was still sufficiently accurate with subpopulations smaller than 1% of the total population. This study proposes a number of rational approximations to quantifying small subpopulations and offers an easy-to-use protocol for their implementation in the open source statistical software environment R.

Low abundance of long-tongued pollinators leads to pollen limitation in four specialized hawkmoth-pollinated plants in the Atlantic Rain forest, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury

Abstract Background The sequencing of the D.melanogaster genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of Drosophila protein-coding genes contain one or more alternative exons. A recent transcription map of the Drosophila embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of Drosophila transcriptome. Results Bioinformatic analysis of 1,303 Drosophila ORESTES clusters identified 68 sequences derived from unannotated regions in the current Drosophila genome version (4.3). Of these, a set of 38 was analysed by polyA+ northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The SP212 gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this locus is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury. Conclusion Using the ORESTES methodology we identified 17 novel exons from low abundance Drosophila transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.

Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources

Microarrays have established as instrumental for bacterial detection, identification, and genotyping as well as for transcriptomic studies. For gene expression analyses using limited numbers of bacteria (derived from in vivo or ex vivo origin, for example), RNA amplification is often required prior to labeling and hybridization onto microarrays. Evaluation of the fidelity of the amplification methods is crucial for the robustness and reproducibility of microarray results. We report here the first utilization of random primers and the highly processive Phi29 phage polymerase to amplify material for transcription profiling analyses. We compared two commercial amplification methods (GenomiPhi and MessageAmp kits) with direct reverse-transcription as the reference method, focusing on the robustness of mRNA quantification using either microarrays or quantitative RT-PCR. Both amplification methods using either poly-A tailing followed by in vitro transcription, or direct strand displacement polymerase, showed appreciable linearity. Strand displacement technique was particularly affordable compared to in vitro transcription-based (IVT) amplification methods and consisted in a single tube reaction leading to high amplification yields. Real-time measurements using low-, medium-, and highly expressed genes revealed that this simple method provided linear amplification with equivalent results in terms of relative messenger abundance as those obtained by conventional direct reverse-transcription.

Contribution of APOBEC3G/F Activity to the Development of Low-Abundance Drug-Resistant HIV-1 variants.

Plasma drug-resistant minority HIV-1 variants (DRMV) increase the risk of virological failure to first-line NNRTI antiretroviral therapy (ART). The origin of DRMVs in ART-naive patients, however, remains unclear. In a large pan-European case-control study investigating the clinical relevance of pre-existing DRMVs using 454 pyrosequencing, the six most prevalent plasma DRMVs detected corresponded to G-to-A nucleotide mutations (V90I, V106I, V108I, E138K, M184I and M230I). Here, we evaluated if such DRMVs could have emerged from APOBEC3G/F activity. Out of 236 ART-naïve evaluated subjects, APOBEC3G/F hypermutation signatures were detected in plasma viruses of 14 (5.9%) individuals. Samples with minority E138K, M184I, and M230I mutations, but not those with V90I, V106I, or V108I were significantly associated with APOBEC3G/F activity (Fisher's p<0.005), defined as presence of >0.5% of sample sequences with an APOBEC3G/F signature. Mutations E138K, M184I and M230I co-occurred in the same sequence as APOBEC3G/F signatures in 3/9 (33%), 5/11 (45%) and 4/8 (50%) of samples, respectively; such linkage was not found for V90I, V106I or V108I. In-frame STOP codons were observed in 1.5% of all clonal sequences; 14.8% of them co-occurred with APOBEC3G/F signatures. APOBEC3G/F-associated E138K, M184I and M230I appeared within clonal sequences containing in-frame STOP codons in 2/3 (66%), 5/5 (100%) and 4/4 (100%) of the samples. In a reanalysis of the parent case-control study, presence of APOBEC3G/F signatures was not associated with virological failure. In conclusion, the contribution of APOBEC3G/F editing to the development of DRMVs is very limited and does not affect the efficacy of NNRTI ART.

The low abundance of U7 snRNA is partly determined by its Sm binding site.

In transient expression studies after DNA transfection of HeLa cells, the mouse U7 gene produces only approximately 30% of the RNA produced by a mouse U1b gene. This difference persists even when the transfected genes have all their 5' and 3' flanking sequences exchanged suggesting a post-transcriptional effect. When the special U7 Sm binding site is mutated to a consensus derived from the major snRNAs (Sm-opt), the U7 RNA level increases 4- to 5-fold, whereas no RNA is detected from a U7 gene with a non-functional Sm binding site (Sm-mut). Moreover, U1b genes with the U7 Sm binding site yield reduced RNA levels. The Sm-opt site also alters the cellular behaviour of the corresponding U7 snRNA. It accumulates to a higher level in the nucleus than wild type U7 RNA, and is better immunoprecipitable with anti-Sm antibodies. Injection experiments in Xenopus oocytes indicate that the U7 genes with either Sm-opt or Sm-mut sites produce similar amounts of RNA as wild type U7, but that they differ in opposing ways in the processing of precursors to mature size U7 snRNA and in nuclear accumulation. However, in reconstitution experiments using Xenopus oocytes, we show that U7 Sm-opt RNA, despite its efficient nuclear accumulation, is not active in 3' processing of histone pre-mRNA, whereas wild type U7 RNA is assembled into functional snRNPs, which correctly process histone pre-mRNA substrate. This suggests a functional importance of the special U7 Sm sequence.

Rapid affinity purification of erythropoietin from biological samples using disposable monoliths.

Identification of post-translational modifications of proteins in biological samples often requires access to preanalytical purification and concentration methods. In the purification step high or low molecular weight substances can be removed by size exclusion filters, and high abundant proteins can be removed, or low abundant proteins can be enriched, by specific capturing tools. In this paper is described the experience and results obtained with a recently emerged and easy-to-use affinity purification kit for enrichment of the low amounts of EPO found in urine and plasma specimens. The kit can be used as a pre-step in the EPO doping control procedure, as an alternative to the commonly used ultrafiltration, for detecting aberrantly glycosylated isoforms. The commercially available affinity purification kit contains small disposable anti-EPO monolith columns (6 ?L volume, Ø7 mm, length 0.15 mm) together with all required buffers. A 24-channel vacuum manifold was used for simultaneous processing of samples. The column concentrated EPO from 20 mL urine down to 55 ?L eluate with a concentration factor of 240 times, while roughly 99.7% of non-relevant urine proteins were removed. The recoveries of Neorecormon (epoetin beta), and the EPO analogues Aranesp and Mircera applied to buffer were high, 76%, 67% and 57%, respectively. The recovery of endogenous EPO from human urine was 65%. High recoveries were also obtained when purifying human, mouse and equine EPO from serum, and human EPO from cerebrospinal fluid. Evaluation with the accredited EPO doping control method based on isoelectric focusing (IEF) showed that the affinity purification procedure did not change the isoform distribution for rhEPO, Aranesp, Mircera or endogenous EPO. The kit should be particularly useful for applications in which it is essential to avoid carry-over effects, a problem commonly encountered with conventional particle-based affinity columns. The encouraging results with EPO propose that similar affinity monoliths, with the appropriate antibodies, should constitute useful tools for general applications in sample preparation, not only for doping control of EPO and other hormones such as growth hormone and insulin but also for the study of post-translational modifications of other low abundance proteins in biological and clinical research, and for sample preparation prior to in vitro diagnostics.

Combinatorial peptide ligand libraries and plant proteomics: a winning strategy at a price

The mechanism of action and properties of a solid-phase ligand library made of hexapeptides (combinatorial peptide ligand libraries or CPLL), for capturing the "hidden proteome", i.e. the low- and very low-abundance proteins constituting the vast majority of species in any proteome, as applied to plant tissues, are reviewed here. Plant tissues are notoriously recalcitrant to protein extraction and to proteome analysis. Firstly, rigid plant cell walls need to be mechanically disrupted to release the cell content and, in addition to their poor protein yield, plant tissues are rich in proteases and oxidative enzymes, contain phenolic compounds, starches, oils, pigments and secondary metabolites that massively contaminate protein extracts. In addition, complex matrices of polysaccharides, including large amount of anionic pectins, are present. All these species compete with the binding of proteins to the CPLL beads, impeding proper capture and identification / detection of low-abundance species. When properly pre-treated, plant tissue extracts are amenable to capture by the CPLL beads revealing thus many new species among them low-abundance proteins. Examples are given on the treatment of leaf proteins, of corn seed extracts and of exudate proteins (latex from Hevea brasiliensis). In all cases, the detection of unique gene products via CPLL capture is at least twice that of control, untreated sample.

Combinatorial peptide ligand libraries and plant proteomics: a winning strategy at a price

The mechanism of action and properties of a solid-phase ligand library made of hexapeptides (combinatorial peptide ligand libraries or CPLL, for capturing the "hidden proteome", i.e. the low- and very low-abundance proteins Constituting the vast majority of species in any proteome. as applied to plant tissues, are reviewed here. Plant tissues are notoriously recalcitrant to protein extraction and to proteome analysis, Firstly, rigid plant cell walls need to be mechanically disrupted to release the cell content and, in addition to their poor protein yield, plant tissues are rich in proteases and oxidative enzymes, contain phenolic Compounds, starches, oils, pigments and secondary metabolites that massively contaminate protein extracts. In addition, complex matrices of polysaccharides, including large amount of anionic pectins, are present. All these species compete with the binding of proteins to the CPLL beads, impeding proper capture and identification I detection of low-abundance species. When properly pre-treated, plant tissue extracts are amenable to capture by the CPLL beads revealing thus many new species among them low-abundance proteins. Examples are given on the treatment of leaf proteins, of corn seed extracts and of exudate proteins (latex from Hevea brasiliensis). In all cases, the detection of unique gene products via CPLL Capture is at least twice that of control, untreated sample. (c) 2008 Elsevier B.V. All rights reserved.

A preliminary comparative analysis of the H. sapiens saliva proteome between cysteine fractionation,performic acid oxidation LC/MALDI/MS/MS and LC/ESI/MS/MS

We present a comparative study between LC/MALDI/MS/MS and LC/ESI/MS/MS. Diagnostic biomarkers in saliva have been identified for monitoring caries, periodontitis, oral cancer, salivary gland diseases, and systemic disorders e.g. hepatitis and HIV[1]. Saliva is similar to serum in that there are a small number of highly abundant proteins and many low abundance proteins. There are 35 previously identified salivary proteins [1-4]. We prepared a representative sample of cysteine containing peptides and oxidised them to improve their fragmentation under MALDI conditions. In total 20 proteins were identified with 6 been identified by both methods. Surprisingly there was little overlap in the peptides used to identify the proteins between the two methods

Separating the wheat from the chaff: unbiased filtering of background tandem mass spectra improves protein identification

Only a small fraction of spectra acquired in LC-MS/MS runs matches peptides from target proteins upon database searches. The remaining, operationally termed background, spectra originate from a variety of poorly controlled sources and affect the throughput and confidence of database searches. Here, we report an algorithm and its software implementation that rapidly removes background spectra, regardless of their precise origin. The method estimates the dissimilarity distance between screened MS/MS spectra and unannotated spectra from a partially redundant background library compiled from several control and blank runs. Filtering MS/MS queries enhanced the protein identification capacity when searches lacked spectrum to sequence matching specificity. In sequence-similarity searches it reduced by, on average, 30-fold the number of orphan hits, which were not explicitly related to background protein contaminants and required manual validation. Removing high quality background MS/MS spectra, while preserving in the data set the genuine spectra from target proteins, decreased the false positive rate of stringent database searches and improved the identification of low-abundance proteins.

Identifizierung neuer potentieller Biomarker für das Kolonkarzinom sowie dessen Vorstufen

Zur Detektion von neuen potentiellen Biomarkerkandidaten für das KRK sowie dessen Vorstufen wurde Blutplasma von gesunden Kontrollen und von Patienten mit kolorektalen Adenomen und Karzinomen untersucht. Mit Hilfe immunoaffinitäts-basierter Chromatographiesäulen wurde das Blutplasma mit einer automatisierten Methode von hoch und moderat abundanten Proteinen entfernt und daraufhin massenspektrometrisch analysiert. Hierfür wurde eine labelfreie quantitative und datenunabhängige Methode angewendet. Nach Auswertung der Daten kristallisierten sich mehrere potentielle Biomarkerkandidaten heraus, wobei das Augenmerk auf solchen Proteinen lag, die sowohl bei Adenompatienten als auch bei Karzinompatienten vermehrt im Plasma vorlagen. Nach Western Blot Experimenten und einer weiteren Validierung mittels eines ELISAs an einem größeren Patientenset zeigte sich, dass die beiden Proteine AMBP und CRP als potentielle Biomarker zur Detektion von Adenomen und Karzinomen herangezogen werden konnten. Receiver Operating Charakteristik (ROC)-Analysen ergaben einen AUC-Wert von 0,79 für AMBP und 0,77 für CRP für das Erkennen von Adenomen und Karzinomen. Dies resultierte in einer Sensitivität von 72,5% für AMBP bzw. 60% für CRP bei einer Spezifität von 80%. Durch die Kombination beider Markerproteine erhielt man sogar eine noch höhere Sensitivität von 80% bei 80% Spezifität für das Aufdecken kolorektaler Neoplasien. Ebenso zeigte sich in der histochemischen Analyse, dass AMBP direkt von intestinalen Epithelzellen exprimiert wird.

Stimulation of mouse mammary tumor virus superantigen expression by an intragenic enhancer.

The mechanisms regulating expression of mouse mammary tumor virus (MMTV)-encoded superantigens from the viral sag gene are largely unknown, due to problems with detection and quantification of these low-abundance proteins. To study the expression and regulation of the MMTV sag gene, we have developed a sensitive and quantitative reporter gene assay based on a recombinant superantigen-human placental alkaline phosphatase fusion protein. High sag-reporter expression in Ba/F3, an early B-lymphoid cell line, depends on enhancers in either of the viral long terminal repeats (LTRs) and is largely independent of promoters in the 5' LTR. The same enhancer region is also required for general expression of MMTV genes from the 5' LTR. The enhancer was mapped to a 548-bp fragment of the MMTV LTR lying within sag and shown to be sufficient to stimulate expression from a heterologous simian virus 40 promoter. No enhancer activity of the MMTV LTR was observed in XC sarcoma cells, which are permissive for MMTV. Our results demonstrate a major role for this enhancer in MMTV gene expression in early B-lymphoid cells.

Proteomic analysis of k-casein micro-heterogeneity

The proteome of bovine milk is dominated by just six gene products that constitute approximately 95% of milk protein. Nonetheless, over 150 protein spots can be readily detected following two-dimensional electrophoresis of whole milk. Many of these represent isoforms of the major gene products produced through extensive posttranslational modification. Peptide mass fingerprinting of in-gel tryptic digests (using matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) in reflectron mode with alpha-cyano-4-hydroxycinnamic acid as the matrix) identified 10 forms of K-casein with isoelectric point (pl) values from 4.47 to 5.81, but could not distinguish between them. MALDI-TOF MS in linear mode, using sinapinic acid as the matrix, revealed a large tryptic peptide (mass > 5990 Da) derived from the C-terminus that contained all the known sites of genetic variance, phosphorylation and glycosylation. Two genetic variants present as singly or doubly phosphorylated forms could be distinguished using mass data alone. Glycoforms containing a single acidic tetrasaccharide were also identified. The differences in electrophoretic mobility of these isoforms were consistent with the addition of the acidic groups. While more extensively glycosylated forms were also observed, substantial loss of N-acetylneuraminic acid from the glycosyl group was evident in the MALDI spectra such that ions corresponding to the intact glycopeptide were not observed and assignment of the glycoforms was not possible. However, by analysing the pl shifts observed on the two-dimensional gels in conjunction with the MS data, the number of N-acetylneuraminic acid residues, and hence the glycoforms present, could be determined.