Evaluation of reference-based two-color methods for measurement of gene expression ratios using spotted cDNA microarrays

Abstract Background Spotted cDNA microarrays generally employ co-hybridization of fluorescently-labeled RNA targets to produce gene expression ratios for subsequent analysis. Direct comparison of two RNA samples in the same microarray provides the highest level of accuracy; however, due to the number of combinatorial pair-wise comparisons, the direct method is impractical for studies including large number of individual samples (e.g., tumor classification studies). For such studies, indirect comparisons using a common reference standard have been the preferred method. Here we evaluated the precision and accuracy of reconstructed ratios from three indirect methods relative to ratios obtained from direct hybridizations, herein considered as the gold-standard. Results We performed hybridizations using a fixed amount of Cy3-labeled reference oligonucleotide (RefOligo) against distinct Cy5-labeled targets from prostate, breast and kidney tumor samples. Reconstructed ratios between all tissue pairs were derived from ratios between each tissue sample and RefOligo. Reconstructed ratios were compared to (i) ratios obtained in parallel from direct pair-wise hybridizations of tissue samples, and to (ii) reconstructed ratios derived from hybridization of each tissue against a reference RNA pool (RefPool). To evaluate the effect of the external references, reconstructed ratios were also calculated directly from intensity values of single-channel (One-Color) measurements derived from tissue sample data collected in the RefOligo experiments. We show that the average coefficient of variation of ratios between intra- and inter-slide replicates derived from RefOligo, RefPool and One-Color were similar and 2 to 4-fold higher than ratios obtained in direct hybridizations. Correlation coefficients calculated for all three tissue comparisons were also similar. In addition, the performance of all indirect methods in terms of their robustness to identify genes deemed as differentially expressed based on direct hybridizations, as well as false-positive and false-negative rates, were found to be comparable. Conclusion RefOligo produces ratios as precise and accurate as ratios reconstructed from a RNA pool, thus representing a reliable alternative in reference-based hybridization experiments. In addition, One-Color measurements alone can reconstruct expression ratios without loss in precision or accuracy. We conclude that both methods are adequate options in large-scale projects where the amount of a common reference RNA pool is usually restrictive.

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.

Sequence features responsible for intron retention in human

Abstract Background One of the least common types of alternative splicing is the complete retention of an intron in a mature transcript. Intron retention (IR) is believed to be the result of intron, rather than exon, definition associated with failure of the recognition of weak splice sites flanking short introns. Although studies on individual retained introns have been published, few systematic surveys of large amounts of data have been conducted on the mechanisms that lead to IR. Results TTo understand how sequence features are associated with or control IR, and to produce a generalized model that could reveal previously unknown signals that regulate this type of alternative splicing, we partitioned intron retention events observed in human cDNAs into two groups based on the relative abundance of both isoforms and compared relevant features. We found that a higher frequency of IR in human is associated with individual introns that have weaker splice sites, genes with shorter intron lengths, higher expression levels and lower density of both a set of exon splicing silencers (ESSs) and the intronic splicing enhancer GGG. Both groups of retained introns presented events conserved in mouse, in which the retained introns were also short and presented weaker splice sites. Conclusion Although our results confirmed that weaker splice sites are associated with IR, they showed that this feature alone cannot explain a non-negligible fraction of events. Our analysis suggests that cis-regulatory elements are likely to play a crucial role in regulating IR and also reveals previously unknown features that seem to influence its occurrence. These results highlight the importance of considering the interplay among these features in the regulation of the relative frequency of IR.

Signal transduction-related responses to phytohormones and environmental challenges in sugarcane

Abstract Background Sugarcane is an increasingly economically and environmentally important C4 grass, used for the production of sugar and bioethanol, a low-carbon emission fuel. Sugarcane originated from crosses of Saccharum species and is noted for its unique capacity to accumulate high amounts of sucrose in its stems. Environmental stresses limit enormously sugarcane productivity worldwide. To investigate transcriptome changes in response to environmental inputs that alter yield we used cDNA microarrays to profile expression of 1,545 genes in plants submitted to drought, phosphate starvation, herbivory and N2-fixing endophytic bacteria. We also investigated the response to phytohormones (abscisic acid and methyl jasmonate). The arrayed elements correspond mostly to genes involved in signal transduction, hormone biosynthesis, transcription factors, novel genes and genes corresponding to unknown proteins. Results Adopting an outliers searching method 179 genes with strikingly different expression levels were identified as differentially expressed in at least one of the treatments analysed. Self Organizing Maps were used to cluster the expression profiles of 695 genes that showed a highly correlated expression pattern among replicates. The expression data for 22 genes was evaluated for 36 experimental data points by quantitative RT-PCR indicating a validation rate of 80.5% using three biological experimental replicates. The SUCAST Database was created that provides public access to the data described in this work, linked to tissue expression profiling and the SUCAST gene category and sequence analysis. The SUCAST database also includes a categorization of the sugarcane kinome based on a phylogenetic grouping that included 182 undefined kinases. Conclusion An extensive study on the sugarcane transcriptome was performed. Sugarcane genes responsive to phytohormones and to challenges sugarcane commonly deals with in the field were identified. Additionally, the protein kinases were annotated based on a phylogenetic approach. The experimental design and statistical analysis applied proved robust to unravel genes associated with a diverse array of conditions attributing novel functions to previously unknown or undefined genes. The data consolidated in the SUCAST database resource can guide further studies and be useful for the development of improved sugarcane varieties.

Interleukin-10 production by tumor infiltrating macrophages plays a role in Human Papillomavirus 16 tumor growth

Abstract Background Human Papillomavirus, HPV, is the main etiological factor for cervical cancer. Different studies show that in women infected with HPV there is a positive correlation between lesion grade and number of infiltrating macrophages, as well as with IL-10 higher expression. Using a HPV16 associated tumor model in mice, TC-1, our laboratory has demonstrated that tumor infiltrating macrophages are M2-like, induce T cell regulatory phenotype and play an important role in tumor growth. M2 macrophages secrete several cytokines, among them IL-10, which has been shown to play a role in T cell suppression by tumor macrophages in other tumor models. In this work, we sought to establish if IL-10 is part of the mechanism by which HPV tumor associated macrophages induce T cell regulatory phenotype, inhibiting anti-tumor activity and facilitating tumor growth. Results TC-1 tumor cells do not express or respond to IL-10, but recruit leukocytes which, within the tumor environment, produce this cytokine. Using IL-10 deficient mice or blocking IL-10 signaling with neutralizing antibodies, we observed a significant reduction in tumor growth, an increase in tumor infiltration by HPV16 E7 specific CD8 lymphocytes, including a population positive for Granzyme B and Perforin expression, and a decrease in the percentage of HPV specific regulatory T cells in the lymph nodes. Conclusions Our data shows that in the HPV16 TC-1 tumor mouse model, IL-10 produced by tumor macrophages induce regulatory phenotype on T cells, an immune escape mechanism that facilitates tumor growth. Our results point to a possible mechanism behind the epidemiologic data that correlates higher IL-10 expression with risk of cervical cancer development in HPV infected women.

Immunomodulation of human monocytes following exposure to Lutzomyia intermedia saliva

Abstract Background Sand fly saliva contains potent and complex pharmacologic molecules that are able to modulate the host's hemostatic, inflammatory, and immune systems. In this study, we evaluated the effects of salivary gland sonicate (SGS) of Lutzomyia intermedia, the natural vector of Leishmania braziliensis, on monocytes obtained from the peripheral blood mononuclear cells (PBMC) of healthy volunteers. We investigated the effects of sand fly saliva on cytokine production and surface molecule expression of LPS-stimulated human monocytes uninfected or infected with L. braziliensis. Results Pre-treatment of non-infected human monocytes with L. intermedia SGS followed by LPS-stimulation led to a significant decrease in IL-10 production accompanied by a significant increase in CD86, CD80, and HLA-DR expression. Pre-treatment with SGS followed by LPS stimulation and L. braziliensis infection led to a significant increase in TNF-α, IL-6, and IL-8 production without significant alterations in co-stimulatory molecule expression. However, pre-treatment with L. intermedia SGS did not result in significant changes in the infection rate of human monocytes. Conclusion Our data indicate that L. intermedia saliva is able to modulate monocyte response, and, although this modulation is dissociated from enhanced infection with L. braziliensis, it may be associated with successful parasitism.

Reliable identification of mycobacterial species by PCR-restriction enzyme analysis (PRA)-hsp65 in a reference laboratory and elaboration of a sequence-based extended algorithm of PRA-hsp65 patterns

Abstract Background Identification of nontuberculous mycobacteria (NTM) based on phenotypic tests is time-consuming, labor-intensive, expensive and often provides erroneous or inconclusive results. In the molecular method referred to as PRA-hsp65, a fragment of the hsp65 gene is amplified by PCR and then analyzed by restriction digest; this rapid approach offers the promise of accurate, cost-effective species identification. The aim of this study was to determine whether species identification of NTM using PRA-hsp65 is sufficiently reliable to serve as the routine methodology in a reference laboratory. Results A total of 434 NTM isolates were obtained from 5019 cultures submitted to the Institute Adolpho Lutz, Sao Paulo Brazil, between January 2000 and January 2001. Species identification was performed for all isolates using conventional phenotypic methods and PRA-hsp65. For isolates for which these methods gave discordant results, definitive species identification was obtained by sequencing a 441 bp fragment of hsp65. Phenotypic evaluation and PRA-hsp65 were concordant for 321 (74%) isolates. These assignments were presumed to be correct. For the remaining 113 discordant isolates, definitive identification was based on sequencing a 441 bp fragment of hsp65. PRA-hsp65 identified 30 isolates with hsp65 alleles representing 13 previously unreported PRA-hsp65 patterns. Overall, species identification by PRA-hsp65 was significantly more accurate than by phenotype methods (392 (90.3%) vs. 338 (77.9%), respectively; p < .0001, Fisher's test). Among the 333 isolates representing the most common pathogenic species, PRA-hsp65 provided an incorrect result for only 1.2%. Conclusion PRA-hsp65 is a rapid and highly reliable method and deserves consideration by any clinical microbiology laboratory charged with performing species identification of NTM.

Expression profiling and cross-species RNA interference (RNAi) of desiccation-induced transcripts in the anhydrobiotic nematode Aphelenchus avenae

Abstract Background Some organisms can survive extreme desiccation by entering a state of suspended animation known as anhydrobiosis. The free-living mycophagous nematode Aphelenchus avenae can be induced to enter anhydrobiosis by pre-exposure to moderate reductions in relative humidity (RH) prior to extreme desiccation. This preconditioning phase is thought to allow modification of the transcriptome by activation of genes required for desiccation tolerance. Results To identify such genes, a panel of expressed sequence tags (ESTs) enriched for sequences upregulated in A. avenae during preconditioning was created. A subset of 30 genes with significant matches in databases, together with a number of apparently novel sequences, were chosen for further study. Several of the recognisable genes are associated with water stress, encoding, for example, two new hydrophilic proteins related to the late embryogenesis abundant (LEA) protein family. Expression studies confirmed EST panel members to be upregulated by evaporative water loss, and the majority of genes was also induced by osmotic stress and cold, but rather fewer by heat. We attempted to use RNA interference (RNAi) to demonstrate the importance of this gene set for anhydrobiosis, but found A. avenae to be recalcitrant with the techniques used. Instead, therefore, we developed a cross-species RNAi procedure using A. avenae sequences in another anhydrobiotic nematode, Panagrolaimus superbus, which is amenable to gene silencing. Of 20 A. avenae ESTs screened, a significant reduction in survival of desiccation in treated P. superbus populations was observed with two sequences, one of which was novel, while the other encoded a glutathione peroxidase. To confirm a role for glutathione peroxidases in anhydrobiosis, RNAi with cognate sequences from P. superbus was performed and was also shown to reduce desiccation tolerance in this species. Conclusions This study has identified and characterised the expression profiles of members of the anhydrobiotic gene set in A. avenae. It also demonstrates the potential of RNAi for the analysis of anhydrobiosis and provides the first genetic data to underline the importance of effective antioxidant systems in metazoan desiccation tolerance.

A theoretical model for estimating the margination constant of leukocytes

Abstract Background Blood leukocytes constitute two interchangeable sub-populations, the marginated and circulating pools. These two sub-compartments are found in normal conditions and are potentially affected by non-normal situations, either pathological or physiological. The dynamics between the compartments is governed by rate constants of margination (M) and return to circulation (R). Therefore, estimates of M and R may prove of great importance to a deeper understanding of many conditions. However, there has been a lack of formalism in order to approach such estimates. The few attempts to furnish an estimation of M and R neither rely on clearly stated models that precisely say which rate constant is under estimation nor recognize which factors may influence the estimation. Results The returning of the blood pools to a steady-state value after a perturbation (e.g., epinephrine injection) was modeled by a second-order differential equation. This equation has two eigenvalues, related to a fast- and to a slow-component of the dynamics. The model makes it possible to identify that these components are partitioned into three constants: R, M and SB; where SB is a time-invariant exit to tissues rate constant. Three examples of the computations are worked and a tentative estimation of R for mouse monocytes is presented. Conclusions This study establishes a firm theoretical basis for the estimation of the rate constants of the dynamics between the blood sub-compartments of white cells. It shows, for the first time, that the estimation must also take into account the exit to tissues rate constant, SB.

Divergent calcium signaling in RBCs from Tropidurus torquatus (Squamata – Tropiduridae) strengthen classification in lizard evolution

Abstract Background We have previously reported that a Teiid lizard red blood cells (RBCs) such as Ameiva ameiva and Tupinambis merianae controls intracellular calcium levels by displaying multiple mechanisms. In these cells, calcium stores could be discharged not only by: thapsigargin, but also by the Na+/H+ ionophore monensin, K+/H+ ionophore nigericin and the H+ pump inhibitor bafilomycin as well as ionomycin. Moreover, these lizards possess a P2Y-type purinoceptors that mobilize Ca2+ from intracellular stores upon ATP addition. Results Here we report, that RBCs from the tropidurid lizard Tropidurus torquatus store Ca2+ in endoplasmic reticulum (ER) pool but unlike in the referred Teiidae, these cells do not store calcium in monensin-nigericin sensitive pools. Moreover, mitochondria from T. torquatus RBCs accumulate Ca2+. Addition of ATP to a calcium-free medium does not increase the [Ca2+]c levels, however in a calcium medium we observe an increase in cytosolic calcium. This is an indication that purinergic receptors in these cells are P2X-like. Conclusion T. torquatus RBCs present different mechanisms from Teiid lizard red blood cells (RBCs), for controlling its intracellular calcium levels. At T. torquatus the ion is only stored at endoplasmic reticulum and mitochondria. Moreover activation of purinergic receptor, P2X type, was able to induce an influx of calcium from extracelullar medium. These studies contribute to the understanding of the evolution of calcium homeostasis and signaling in nucleated RBCs.

Protein secretion in Lactococcus lactis : an efficient way to increase the overall heterologous protein production

Anderson MIYOSHI, Daniela FREITAS, Luciana RIBEIRO, Jane E. GABRIEL, Sophie LECLERCQ, Maricê N. OLIVEIRA, and Valeria D. GUIMARÃES were recipients of a CAPES fellowship (project CAPESCOFECUB #319II). Luis BERMUDEZ and Sébastien NOUAILLE were recipients of a fellowship from the French Ministry of Education and Research. INRA and Région IledeFrance also financed L. BERMUDEZ and V. GUIMARAES. Cathy CHARLIER is recipient of a fellowship from INRA and Région Bretagne.

Biological evaluation of hydroxynaphthoquinones as anti-malarials

Abstract Background The hydroxynaphthoquinones have been extensively investigated over the past 50 years for their anti-malarial activity. One member of this class, atovaquone, is combined with proguanil in Malarone®, an important drug for the treatment and prevention of malaria. Methods Anti-malarial activity was assessed in vitro for a series of 3-alkyl-2-hydroxy-1,4-naphthoquinones (N1-N5) evaluating the parasitaemia after 48 hours of incubation. Potential cytotoxicity in HEK293T cells was assessed using the MTT assay. Changes in mitochondrial membrane potential of Plasmodium were measured using the fluorescent dye Mitrotracker Red CMXROS. Results Four compounds demonstrated IC50s in the mid-micromolar range, and the most active compound, N3, had an IC50 of 443 nM. N3 disrupted mitochondrial membrane potential, and after 1 hour presented an IC50ΔΨmit of 16 μM. In an in vitro cytotoxicity assay using HEK 293T cells N3 demonstrated no cytotoxicity at concentrations up to 16 μM. Conclusions N3 was a potent inhibitor of mitochondrial electron transport, had nanomolar activity against cultured Plasmodium falciparum and showed minimal cytotoxicity. N3 may serve as a starting point for the design of new hydroxynaphthoquinone anti-malarials.

Genetic diversity of Leishmania amazonensis strains isolated in northeastern Brazil as revealed by DNA sequencing, PCR-based analyses and molecular karyotyping

Abstract Background Leishmania (Leishmania) amazonensis infection in man results in a clinical spectrum of disease manifestations ranging from cutaneous to mucosal or visceral involvement. In the present study, we have investigated the genetic variability of 18 L. amazonensis strains isolated in northeastern Brazil from patients with different clinical manifestations of leishmaniasis. Parasite DNA was analyzed by sequencing of the ITS flanking the 5.8 S subunit of the ribosomal RNA genes, by RAPD and SSR-PCR and by PFGE followed by hybridization with gene-specific probes. Results ITS sequencing and PCR-based methods revealed genetic heterogeneity among the L. amazonensis isolates examined and molecular karyotyping also showed variation in the chromosome size of different isolates. Unrooted genetic trees separated strains into different groups. Conclusion These results indicate that L. amazonensis strains isolated from leishmaniasis patients from northeastern Brazil are genetically diverse, however, no correlation between genetic polymorphism and phenotype were found.

Hormone-regulated expression and distribution of versican in mouse uterine tissues

Abstract Background Remodeling of the extracellular matrix is one of the most striking features observed in the uterus during the estrous cycle and after hormone replacement. Versican (VER) is a hyaluronan-binding proteoglycan that undergoes RNA alternative splicing, generating four distinct isoforms. This study analyzed the synthesis and distribution of VER in mouse uterine tissues during the estrous cycle, in ovariectomized (OVX) animals and after 17beta-estradiol (E2) and medroxyprogesterone (MPA) treatments, either alone or in combination. Methods Uteri from mice in all phases of the estrous cycle, and animals subjected to ovariectomy and hormone replacement were collected for immunoperoxidase staining for versican, as well as PCR and quantitative Real Time PCR. Results In diestrus and proestrus, VER was exclusively expressed in the endometrial stroma. In estrus and metaestrus, VER was present in both endometrial stroma and myometrium. In OVX mice, VER immunoreaction was abolished in all uterine tissues. VER expression was restored by E2, MPA and E2+MPA treatments. Real Time PCR analysis showed that VER expression increases considerably in the MPA-treated group. Analysis of mRNA identified isoforms V0, V1 and V3 in the mouse uterus. Conclusion These results show that the expression of versican in uterine tissues is modulated by ovarian steroid hormones, in a tissue-specific manner. VER is induced in the myometrium exclusively by E2, whereas MPA induces VER deposition only in the endometrial stroma.

B cells Can Modulate the CD8 Memory T Cell after DNA Vaccination Against Experimental Tuberculosis

Abstract Background Although B cells are important as antigen presenting cells (APC) during the immune response, their role in DNA vaccination models is unknown. Methods In this study in vitro and in vivo experiments were performed to evaluate the ability of B cells to protect mice against Mycobacterium tuberculosis challenge. Results In vitro and in vivo studies showed that B cells efficiently present antigens after naked plasmid pcDNA3 encoding M. leprae 65-kDa heat shock protein (pcDNA3-Hsp65) internalization and protect B knock-out (BKO) mice against Mycobacterium tuberculosis infection. pcDNA3-Hsp65-transfected B cells adoptively transferred into BKO mice rescued the memory phenotypes and reduced the number of CFU compared to wild-type mice. Conclusions These data not only suggest that B cells play an important role in the induction of CD8 T cells but also that they improve bacterial clearance in DNA vaccine model.