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.

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.

Global gene expression profiling of oral cavity cancers suggests molecular heterogeneity within anatomic subsites

Abstract Background Oral squamous cell carcinoma (OSCC) is a frequent neoplasm, which is usually aggressive and has unpredictable biological behavior and unfavorable prognosis. The comprehension of the molecular basis of this variability should lead to the development of targeted therapies as well as to improvements in specificity and sensitivity of diagnosis. Results Samples of primary OSCCs and their corresponding surgical margins were obtained from male patients during surgery and their gene expression profiles were screened using whole-genome microarray technology. Hierarchical clustering and Principal Components Analysis were used for data visualization and One-way Analysis of Variance was used to identify differentially expressed genes. Samples clustered mostly according to disease subsite, suggesting molecular heterogeneity within tumor stages. In order to corroborate our results, two publicly available datasets of microarray experiments were assessed. We found significant molecular differences between OSCC anatomic subsites concerning groups of genes presently or potentially important for drug development, including mRNA processing, cytoskeleton organization and biogenesis, metabolic process, cell cycle and apoptosis. Conclusion Our results corroborate literature data on molecular heterogeneity of OSCCs. Differences between disease subsites and among samples belonging to the same TNM class highlight the importance of gene expression-based classification and challenge the development of targeted therapies.

Adjunct effect of the antimicrobial photodynamic therapy to an association of non-surgical and surgical periodontal treatment in modulation of gene expression

Background: This study has evaluated the effect of antimicrobial photodynamic therapy (aPDT) used in conjunction with non-surgical and surgical periodontal treatment (PT) in modulating gene expression during periodontal wound healing. Methods: Fifteen patients with chronic periodontitis, presenting bilaterally lower molars with class III furcation lesions and scheduled for extraction, were selected. In initial therapy, scaling and root planing (SRP) was performed in the Control Group (CG), while SRP + aPDT were performed in the Test Group (TG). 45 days later, flap surgery plus SRP, and flap surgery plus SRP + aPDT were performed in the CG and TG, respectively. At 21 days post-surgery, the newly formed granulation tissue was collected, and Real-time PCR evaluated the expression of the genes: tumor necrosis factor-?, interleukin-1?, interleukin-4, interleukin-10, matrix metalloproteinase-2 (MMP-2), tissue inhibitor of metalloproteinase-2 (TIMP-2), osteoprotegerin (OPG), receptor activator of nuclear factor- ?B ligand (RANKL), type I collagen, alkaline phosphatase, osteopontin, osteocalcin, and bone sialoprotein. Results: There were statistically significant differences between the groups in relation to mRNA levels for MMP-2 (TG = 3.26 ± 0.89; CG = 4.23 ± 0.97; p = 0.01), TIMP-2/MMP-2 ratio (TG = 0.91 ± 0.34; CG = 0.73 ± 0.32; p = 0.04), OPG (TG = 0.84 ± 0.45; CG = 0.30 ± 0.26; p = 0.001), and OPG/RANKL ratio (TG = 0.60 ± 0.86; CG = 0.23 ± 0.16; p = 0.04), favoring the TG. Conclusion: The present data suggest that the aPDT associated to nonsurgical and surgical periodontal therapy may modulate the extracellular matrix and bone remodeling by up regulating the TIMP- 2/MMP-2 and OPG/RANKL mRNA ratio, but the clinical relevance needs to be evaluated in further studies.

Independent of ErbB1 gene copy number, EGF stimulates migration but is not associated with cell proliferation in non-small cell lung cancer

Abstract Background Lung cancer often exhibits molecular changes, such as the overexpression of the ErbB1 gene. ErbB1 encodes epidermal growth factor receptor (EGFR), a tyrosine kinase receptor, involved mainly in cell proliferation and survival. EGFR overexpression has been associated with more aggressive disease, poor prognosis, low survival rate and low response to therapy. ErbB1 amplification and mutation are associated with tumor development and are implicated in ineffective treatment. The aim of the present study was to investigate whether the ErbB1 copy number affects EGFR expression, cell proliferation or cell migration by comparing two different cell lines. Methods The copies of ErbB1 gene was evaluated by FISH. Immunofluorescence and Western blotting were performed to determine location and expression of proteins mentioned in the present study. Proliferation was studied by flow cytometry and cell migration by wound healing assay and time lapse. Results We investigated the activation and function of EGFR in the A549 and HK2 lung cancer cell lines, which contain 3 and 6 copies of ErbB1, respectively. The expression of EGFR was lower in the HK2 cell line. EGFR was activated after stimulation with EGF in both cell lines, but this activation did not promote differences in cellular proliferation when compared to control cells. Inhibiting EGFR with AG1478 did not modify cellular proliferation, confirming previous data. However, we observed morphological alterations, changes in microfilament organization and increased cell migration upon EGF stimulation. However, these effects did not seem to be consequence of an epithelial-mesenchymal transition. Conclusion EGFR expression did not appear to be associated to the ErbB1 gene copy number, and neither of these aspects appeared to affect cell proliferation. However, EGFR activation by EGF resulted in cell migration stimulation in both cell lines.