Bayesian model accounting for within-class biological variability in Serial Analysis of Gene Expression (SAGE)

Abstract Background An important challenge for transcript counting methods such as Serial Analysis of Gene Expression (SAGE), "Digital Northern" or Massively Parallel Signature Sequencing (MPSS), is to carry out statistical analyses that account for the within-class variability, i.e., variability due to the intrinsic biological differences among sampled individuals of the same class, and not only variability due to technical sampling error. Results We introduce a Bayesian model that accounts for the within-class variability by means of mixture distribution. We show that the previously available approaches of aggregation in pools ("pseudo-libraries") and the Beta-Binomial model, are particular cases of the mixture model. We illustrate our method with a brain tumor vs. normal comparison using SAGE data from public databases. We show examples of tags regarded as differentially expressed with high significance if the within-class variability is ignored, but clearly not so significant if one accounts for it. Conclusion Using available information about biological replicates, one can transform a list of candidate transcripts showing differential expression to a more reliable one. Our method is freely available, under GPL/GNU copyleft, through a user friendly web-based on-line tool or as R language scripts at supplemental web-site.

SAGE analysis highlights the putative role of underexpression of ribosomal proteins in GH-secreting pituitary adenomas

Background: Although the molecular pathogenesis of pituitary adenomas has been assessed by several different techniques, it still remains partially unclear. Ribosomal proteins (RPs) have been recently related to human tumorigenesis, but they have not yet been evaluated in pituitary tumorigenesis. Objective: The aim of this study was to introduce serial analysis of gene expression (SAGE), a high-throughput method, in pituitary research in order to compare differential gene expression. Methods: Two SAGE cDNA libraries were constructed, one using a pool of mRNA obtained from five GH-secreting pituitary tumors and another from three normal pituitaries. Genes differentially expressed between the libraries were further validated by real-time PCR in 22 GH-secreting pituitary tumors and in 15 normal pituitaries. Results: Computer-generated genomic analysis tools identified 13 722 and 14 993 exclusive genes in normal and adenoma libraries respectively. Both shared 6497 genes, 2188 were underexpressed and 4309 overexpressed in tumoral library. In adenoma library, 33 genes encoding RPs were underexpressed. Among these, RPSA, RPS3, RPS14, and RPS29 were validated by real-time PCR. Conclusion: We report the first SAGE library from normal pituitary tissue and GH-secreting pituitary tumor, which provide quantitative assessment of cellular transcriptome. We also validated some downregulated genes encoding RPs. Altogether, the present data suggest that the underexpression of the studied RP genes possibly collaborates directly or indirectly with other genes to modify cell cycle arrest, DNA repair, and apoptosis, leading to an environment that might have a putative role in the tumorigenesis, introducing new perspectives for further studies on molecular genesis of somatotrophinomas.

The presence of anti-citrullinated protein antibodies (ACPA) and rheumatoid factor on patients with rheumatoid arthritis (RA) does not interfere with the chance of clinical remission in a follow-up of 3 years

Autoantibodies in early rheumatoid arthritis (RA) have important diagnostic value. The association between the presence of autoantibodies against cyclic citrullinated peptide and the response to treatment is controversial. To prospectively evaluate a cohort of patients with early rheumatoid arthritis (< 12 months of symptoms) in order to determine the association between serological markers (rheumatoid factor (RF), anti-citrullinated protein antibodies) such as anti-cyclic citrullinated peptide antibodies (anti-CCP) and citrullinated anti-vimentin (anti-Sa) with the occurrence of clinical remission, forty patients diagnosed with early RA at the time of diagnosis were evaluated and followed for 3 years, in use of standardized therapeutic treatment. Demographic and clinical data were recorded, disease activity score 28 (DAS 28), as well as serology tests (ELISA) for RF (IgM, IgG, and IgA), anti-CCP (CCP2, CCP3, and CCP3.1) and anti-Sa in the initial evaluation and at 3, 6, 12, 18, 24, and 36 months of follow-up. The outcome evaluated was the percentage of patients with clinical remission, which was defined by DAS 28 lower than 2.6. Comparisons were made through the Student t test, mixed-effects regression analysis, and analysis of variance (significance level of 5%). The mean age was 45 years, and a female predominance was observed (90%). At the time of diagnosis, RF was observed in 50% of cases (RF IgA-42%, RF IgG-30%, and RF IgM-50%), anti-CCP in 50% (no difference between CCP2, CCP3, and CCP3.1) and anti-Sa in 10%. After 3 years, no change in the RF prevalence and anti-CCP was observed, but the anti-Sa increased to 17.5% (P = 0.001). The percentage of patients in remission, low, moderate, and intense disease activity, according to the DAS 28, was of 0, 0, 7.5, and 92.5% (initial evaluation) and 22.5, 7.5, 32.5, and 37.5% (after 3 years). There were no associations of the presence of autoantibodies in baseline evaluation and in serial analysis with the percentage of clinical remission during follow-up of 3 years The presence of autoantibodies in early RA has no predictive value for clinical remission in early RA.

Changes in hippocampal gene expression by 7-nitroindazole in rats submitted to forced swimming stress

Nitric oxide (NO) is an atypical neurotransmitter that has been related to the pathophysiology of major depression disorder. Increased plasma NO levels have been reported in depressed and suicidal patients. Inhibition of neuronial nitric oxide synthase (nNOS), on the other hand, induces antidepressant effects in clinical and pre-clinical trials. The mechanisms responsible for the antidepressant-like effects of nNOS inhibitors, however, are not completely understood. In this study, genomic and proteomic analyses were used to investigate the effects of the preferential nNOS inhibitor 7-nitroindazole (7-NI) on changes in global gene and protein expression in the hippocampus of rats submitted to forced swimming test (FST). Chronic treatment (14 days, i.p.) with imipramine (15 mg/kg daily) or 7-NI (60 mg/kg daily) significantly reduced immobility in the FST. Saturation curves for Serial analysis of gene expression libraries showed that the hippocampus of animals submitted to FST presented a lower number of expressed genes compared to non-FST stressed groups. Imipramine, but not 7-NI, reverted this effect. GeneGo analyses revealed that genes related to oxidative phosphorylation, apoptosis and survival controlled by HTR1A signaling and cytoskeleton remodeling controlled by Rho GTPases were significantly changed by FST. 7-NI prevented this effect. In addition, 7-NI treatment changed the expression of genes related to transcription in the cAMP response element-binding pathway. Therefore, this study suggests that changes in oxidative stress and neuroplastic processes could be involved in the antidepressant-like effects induced by nNOS inhibition.

A quantitative view of the transcriptome of Schistosoma mansoni adult-worms using SAGE

Abstract Background Five species of the genus Schistosoma, a parasitic trematode flatworm, are causative agents of Schistosomiasis, a disease that is endemic in a large number of developing countries, affecting millions of patients around the world. By using SAGE (Serial Analysis of Gene Expression) we describe here the first large-scale quantitative analysis of the Schistosoma mansoni transcriptome, one of the most epidemiologically relevant species of this genus. Results After extracting mRNA from pooled male and female adult-worms, a SAGE library was constructed and sequenced, generating 68,238 tags that covered more than 6,000 genes expressed in this developmental stage. An analysis of the ordered tag-list shows the genes of F10 eggshell protein, pol-polyprotein, HSP86, 14-3-3 and a transcript yet to be identified to be the five top most abundant genes in pooled adult worms. Whereas only 8% of the 100 most abundant tags found in adult worms of S. mansoni could not be assigned to transcripts of this parasite, 46.9% of the total ditags could not be mapped, demonstrating that the 3 sequence of most of the rarest transcripts are still to be identified. Mapping of our SAGE tags to S. mansoni genes suggested the occurrence of alternative-polyadenylation in at least 13 gene transcripts. Most of these events seem to shorten the 3 UTR of the mRNAs, which may have consequences over their stability and regulation. Conclusion SAGE revealed the frequency of expression of the majority of the S. mansoni genes. Transcriptome data suggests that alternative polyadenylation is likely to be used in the control of mRNA stability in this organism. When transcriptome was compared with the proteomic data available, we observed a correlation of about 50%, suggesting that both transcriptional and post-transcriptional regulation are important for determining protein abundance in S. mansoni. The generation of SAGE tags from other life-cycle stages should contribute to reveal the dynamics of gene expression in this important parasite.

Gene by environment QTL mapping through multiple trait analyses in blood pressure salt-sensitivity: identification of a novel QTL in rat chromosome 5

Background The genetic mechanisms underlying interindividual blood pressure variation reflect the complex interplay of both genetic and environmental variables. The current standard statistical methods for detecting genes involved in the regulation mechanisms of complex traits are based on univariate analysis. Few studies have focused on the search for and understanding of quantitative trait loci responsible for gene × environmental interactions or multiple trait analysis. Composite interval mapping has been extended to multiple traits and may be an interesting approach to such a problem. Methods We used multiple-trait analysis for quantitative trait locus mapping of loci having different effects on systolic blood pressure with NaCl exposure. Animals studied were 188 rats, the progenies of an F2 rat intercross between the hypertensive and normotensive strain, genotyped in 179 polymorphic markers across the rat genome. To accommodate the correlational structure from measurements taken in the same animals, we applied univariate and multivariate strategies for analyzing the data. Results We detected a new quantitative train locus on a region close to marker R589 in chromosome 5 of the rat genome, not previously identified through serial analysis of individual traits. In addition, we were able to justify analytically the parametric restrictions in terms of regression coefficients responsible for the gain in precision with the adopted analytical approach. Conclusion Future work should focus on fine mapping and the identification of the causative variant responsible for this quantitative trait locus signal. The multivariable strategy might be valuable in the study of genetic determinants of interindividual variation of antihypertensive drug effectiveness.

Vascular smooth muscle cells exhibit a progressive loss of rigidity with serial culture passaging

One drawback of in vitro cell culturing is the dedifferentiation process that cells experience. Smooth muscle cells (SMC) also change molecularly and morphologically with long term culture. The main objective of this study was to evaluate if culture passages interfere in vascular SMC mechanical behavior. SMC were obtained from five different porcine arterial beds. Optical magnetic twisting cytometry (OMTC) was used to characterize mechanically vascular SMC from different cultures in distinct passages and confocal microscopy/western blotting, to evaluate cytoskeleton and extracellular matrix proteins. We found that vascular SMC rigidity or viscoelastic complex modulus (G) decreases with progression of passages. A statistically significant negative correlation between G and passage was found in four of our five cultures studied. Phalloidin-stained SMC from higher passages exhibited lower mean signal intensity per cell (confocal microscopy) and quantitative western blotting analysis showed a decrease in collagen I content throughout passages. We concluded that vascular SMC progressively lose their stiffness with serial culture passaging. Thus, limiting the number of passages is essential for any experiment measuring viscoelastic properties of SMC in culture.

Longitudinal familial analysis of blood pressure involving parametric (co)variance functions

Abstract Background For analyzing longitudinal familial data we adopted a log-linear form to incorporate heterogeneity in genetic variance components over the time, and additionally a serial correlation term in the genetic effects at different levels of ages. Due to the availability of multiple measures on the same individual, we permitted environmental correlations that may change across time. Results Systolic blood pressure from family members from the first and second cohort was used in the current analysis. Measures of subjects receiving hypertension treatment were set as censored values and they were corrected. An initial check of the variance and covariance functions proposed for analyzing longitudinal familial data, using empirical semi-variogram plots, indicated that the observed trait dispersion pattern follows the assumptions adopted. Conclusion The corrections for censored phenotypes based on ordinary linear models may be an appropriate simple model to correct the data, ensuring that the original variability in the data was retained. In addition, empirical semi-variogram plots are useful for diagnosis of the (co)variance model adopted.

Performance of six diagnostic tests to screen for Chagas disease in blood banks and prevalence of Trypanosoma cruzi infection among donors with inconclusive serology screening based on the analysis of epidemiological variables

OBJECTIVE: The frequent occurrence of inconclusive serology in blood banks and the absence of a gold standard test for Chagas'disease led us to examine the efficacy of the blood culture test and five commercial tests (ELISA, IIF, HAI, c-ELISA, rec-ELISA) used in screening blood donors for Chagas disease, as well as to investigate the prevalence of Trypanosoma cruzi infection among donors with inconclusive serology screening in respect to some epidemiological variables. METHODS: To obtain estimates of interest we considered a Bayesian latent class model with inclusion of covariates from the logit link. RESULTS: A better performance was observed with some categories of epidemiological variables. In addition, all pairs of tests (excluding the blood culture test) presented as good alternatives for both screening (sensitivity > 99.96% in parallel testing) and for confirmation (specificity > 99.93% in serial testing) of Chagas disease. The prevalence of 13.30% observed in the stratum of donors with inconclusive serology, means that probably most of these are non-reactive serology. In addition, depending on the level of specific epidemiological variables, the absence of infection can be predicted with a probability of 100% in this group from the pairs of tests using parallel testing. CONCLUSION: The epidemiological variables can lead to improved test results and thus assist in the clarification of inconclusive serology screening results. Moreover, all combinations of pairs using the five commercial tests are good alternatives to confirm results.