Murine Dendritic Cells Transcriptional Modulation upon Paracoccidioides brasiliensis Infection

Limited information is available regarding the modulation of genes involved in the innate host response to Paracoccidioides brasiliensis, the etiologic agent of paracoccidioidomycosis. Therefore, we sought to characterize, for the first time, the transcriptional profile of murine bone marrow-derived dendritic cells (DCs) at an early stage following their initial interaction with P. brasiliensis. DCs connect innate and adaptive immunity by recognizing invading pathogens and determining the type of effector T-cell that mediates an immune response. Gene expression profiles were analyzed using microarray and validated using real-time RT-PCR and protein secretion studies. A total of 299 genes were differentially expressed, many of which are involved in immunity, signal transduction, transcription and apoptosis. Genes encoding the cytokines IL-12 and TNF-alpha, along with the chemokines CCL22, CCL27 and CXCL10, were up-regulated, suggesting that P. brasiliensis induces a potent proinflammatory response in DCs. In contrast, pattern recognition receptor (PRR)-encoding genes, particularly those related to Toll-like receptors, were down-regulated or unchanged. This result prompted us to evaluate the expression profiles of dectin-1 and mannose receptor, two other important fungal PRRs that were not included in the microarray target cDNA sequences. Unlike the mannose receptor, the dectin-1 receptor gene was significantly induced, suggesting that this beta-glucan receptor participates in the recognition of P. brasiliensis. We also used a receptor inhibition assay to evaluate the roles of these receptors in coordinating the expression of several immune-related genes in DCs upon fungal exposure. Altogether, our results provide an initial characterization of early host responses to P. brasiliensis and a basis for better understanding the infectious process of this important neglected pathogen.

BPIFB1 (LPLUNC1) is upregulated in cystic fibrosis lung disease

Although the biology the PLUNC (recently renamed BPI fold, BPIF) family of secreted proteins is poorly understood, multiple array based studies have suggested that some are differentially expressed in lung diseases. We have examined the expression of BPIFB1 (LPLUNC1), the prototypic two-domain containing family member, in lungs from CF patients and in mouse models of CF lung disease. BPIFB1 was localized in CF lung samples along with BPIFA1, MUC5AC, CD68 and NE and directly compared to histologically normal lung tissues and that of bacterial pneumonia. We generated novel antibodies to mouse BPIF proteins to conduct similar studies on ENaC transgenic (ENaC-Tg) mice, a model for CF-like lung disease. Small airways in CF demonstrated marked epithelial staining of BPIFB1 in goblet cells but staining was absent from alveolar regions. BPIFA1 and BPIFB1 were not co-localised in the diseased lungs. In ENaC-Tg mice there was strong staining of both proteins in the airways and luminal contents. This was most marked for BPIFB1 and was noted within 2 weeks of birth. The two proteins were present in distinct cells within epithelium. BPIFB1 was readily detected in BAL from ENaC-Tg mice but was absent from wild-type mice. Alterations in the expression of BPIF proteins is associated with CF lung disease in humans and mice. It is unclear if this elevation of protein production, which results from phenotypic alteration of the cells within the diseased epithelium, plays a role in the pathogenesis of the disease.

A Bayesian Approach for Decision Making on the Identification of Genes with Different Expression Levels: An Application to Escherichia coli Bacterium Data

A common interest in gene expression data analysis is to identify from a large pool of candidate genes the genes that present significant changes in expression levels between a treatment and a control biological condition. Usually, it is done using a statistic value and a cutoff value that are used to separate the genes differentially and nondifferentially expressed. In this paper, we propose a Bayesian approach to identify genes differentially expressed calculating sequentially credibility intervals from predictive densities which are constructed using the sampled mean treatment effect from all genes in study excluding the treatment effect of genes previously identified with statistical evidence for difference. We compare our Bayesian approach with the standard ones based on the use of the t-test and modified t-tests via a simulation study, using small sample sizes which are common in gene expression data analysis. Results obtained report evidence that the proposed approach performs better than standard ones, especially for cases with mean differences and increases in treatment variance in relation to control variance. We also apply the methodologies to a well-known publicly available data set on Escherichia coli bacterium.

Metabolic engineering of beta-oxidation in Penicillium chrysogenum for improved semi-synthetic cephalosporin biosynthesis

Industrial production of semi-synthetic cephalosporins by Penicillium chrysogenum requires supplementation of the growth media with the side-chain precursor adipic acid. In glucose-limited chemostat cultures of P. chrysogenum, up to 88% of the consumed adipic acid was not recovered in cephalosporinrelated products, but used as an additional carbon and energy source for growth. This low efficiency of side-chain precursor incorporation provides an economic incentive for studying and engineering the metabolism of adipic acid in P. cluysogenum. Chemostat-based transcriptome analysis in the presence and absence of adipic acid confirmed that adipic acid metabolism in this fungus occurs via beta-oxidation. A set of 52 adipate-responsive genes included six putative genes for acyl-CoA oxidases and dehydrogenases, enzymes responsible for the first step of beta-oxidation. Subcellular localization of the differentially expressed acyl-CoA oxidases and dehydrogenases revealed that the oxidases were exclusively targeted to peroxisomes, while the dehydrogenases were found either in peroxisomes or in mitochondria. Deletion of the genes encoding the peroxisomal acyl-CoA oxidase Pc20g01800 and the mitochondrial acyl-CoA dehydrogenase Pc20g07920 resulted in a 1.6- and 3.7-fold increase in the production of the semi-synthetic cephalosporin intermediate adipoyl-6-APA, respectively. The deletion strains also showed reduced adipate consumption compared to the reference strain, indicating that engineering of the first step of beta-oxidation successfully redirected a larger fraction of adipic acid towards cephalosporin biosynthesis. (C) 2012 Elsevier Inc. All rights reserved.

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.

Molecular determinants of caste differentiation in the highly eusocial honeybee Apis mellifera

Abstract Background In honeybees, differential feeding of female larvae promotes the occurrence of two different phenotypes, a queen and a worker, from identical genotypes, through incremental alterations, which affect general growth, and character state alterations that result in the presence or absence of specific structures. Although previous studies revealed a link between incremental alterations and differential expression of physiometabolic genes, the molecular changes accompanying character state alterations remain unknown. Results By using cDNA microarray analyses of >6,000 Apis mellifera ESTs, we found 240 differentially expressed genes (DEGs) between developing queens and workers. Many genes recorded as up-regulated in prospective workers appear to be unique to A. mellifera, suggesting that the workers' developmental pathway involves the participation of novel genes. Workers up-regulate more developmental genes than queens, whereas queens up-regulate a greater proportion of physiometabolic genes, including genes coding for metabolic enzymes and genes whose products are known to regulate the rate of mass-transforming processes and the general growth of the organism (e.g., tor). Many DEGs are likely to be involved in processes favoring the development of caste-biased structures, like brain, legs and ovaries, as well as genes that code for cytoskeleton constituents. Treatment of developing worker larvae with juvenile hormone (JH) revealed 52 JH responsive genes, specifically during the critical period of caste development. Using Gibbs sampling and Expectation Maximization algorithms, we discovered eight overrepresented cis-elements from four gene groups. Graph theory and complex networks concepts were adopted to attain powerful graphical representations of the interrelation between cis-elements and genes and objectively quantify the degree of relationship between these entities. Conclusion We suggest that clusters of functionally related DEGs are co-regulated during caste development in honeybees. This network of interactions is activated by nutrition-driven stimuli in early larval stages. Our data are consistent with the hypothesis that JH is a key component of the developmental determination of queen-like characters. Finally, we propose a conceptual model of caste differentiation in A. mellifera based on gene-regulatory networks.

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.

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.

Global gene expression under nitrogen starvation in Xylella fastidiosa: contribution of the σ54 regulon

Abstract Background Xylella fastidiosa, a Gram-negative fastidious bacterium, grows in the xylem of several plants causing diseases such as citrus variegated chlorosis. As the xylem sap contains low concentrations of amino acids and other compounds, X. fastidiosa needs to cope with nitrogen limitation in its natural habitat. Results In this work, we performed a whole-genome microarray analysis of the X. fastidiosa nitrogen starvation response. A time course experiment (2, 8 and 12 hours) of cultures grown in defined medium under nitrogen starvation revealed many differentially expressed genes, such as those related to transport, nitrogen assimilation, amino acid biosynthesis, transcriptional regulation, and many genes encoding hypothetical proteins. In addition, a decrease in the expression levels of many genes involved in carbon metabolism and energy generation pathways was also observed. Comparison of gene expression profiles between the wild type strain and the rpoN null mutant allowed the identification of genes directly or indirectly induced by nitrogen starvation in a σ54-dependent manner. A more complete picture of the σ54 regulon was achieved by combining the transcriptome data with an in silico search for potential σ54-dependent promoters, using a position weight matrix approach. One of these σ54-predicted binding sites, located upstream of the glnA gene (encoding glutamine synthetase), was validated by primer extension assays, confirming that this gene has a σ54-dependent promoter. Conclusions Together, these results show that nitrogen starvation causes intense changes in the X. fastidiosa transcriptome and some of these differentially expressed genes belong to the σ54 regulon.

Identification of protein-coding and non-coding RNA expression profiles in CD34+ and in stromal cells in refractory anemia with ringed sideroblasts

Abstract Background Myelodysplastic syndromes (MDS) are a group of clonal hematological disorders characterized by ineffective hematopoiesis with morphological evidence of marrow cell dysplasia resulting in peripheral blood cytopenia. Microarray technology has permitted a refined high-throughput mapping of the transcriptional activity in the human genome. Non-coding RNAs (ncRNAs) transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression, and in the regulation of exon-skipping and intron retention. Characterization of ncRNAs in progenitor cells and stromal cells of MDS patients could be strategic for understanding gene expression regulation in this disease. Methods In this study, gene expression profiles of CD34+ cells of 4 patients with MDS of refractory anemia with ringed sideroblasts (RARS) subgroup and stromal cells of 3 patients with MDS-RARS were compared with healthy individuals using 44 k combined intron-exon oligoarrays, which included probes for exons of protein-coding genes, and for non-coding RNAs transcribed from intronic regions in either the sense or antisense strands. Real-time RT-PCR was performed to confirm the expression levels of selected transcripts. Results In CD34+ cells of MDS-RARS patients, 216 genes were significantly differentially expressed (q-value ≤ 0.01) in comparison to healthy individuals, of which 65 (30%) were non-coding transcripts. In stromal cells of MDS-RARS, 12 genes were significantly differentially expressed (q-value ≤ 0.05) in comparison to healthy individuals, of which 3 (25%) were non-coding transcripts. Conclusions These results demonstrated, for the first time, the differential ncRNA expression profile between MDS-RARS and healthy individuals, in CD34+ cells and stromal cells, suggesting that ncRNAs may play an important role during the development of myelodysplastic syndromes.

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.

Renal proteome in mice with different susceptibilities to fluorosis

A/J and 129P3/J mouse strains have different susceptibilities to dental fluorosis due to their genetic backgrounds. They also differ with respect to several features of fluoride (F) metabolism and metabolic handling of water. This study was done to determine whether differences in F metabolism could be explained by diversities in the profile of protein expression in kidneys. Weanling, male A/J mice (susceptible to dental fluorosis, n = 18) and 129P3/J mice (resistant, n = 18) were housed in pairs and assigned to three groups given low-F food and drinking water containing 0, 10 or 50 ppm [F] for 7 weeks. Renal proteome profiles were examined using 2D-PAGE and LC-MS/MS. Quantitative intensity analysis detected between A/J and 129P3/J strains 122, 126 and 134 spots differentially expressed in the groups receiving 0, 10 and 50 ppmF, respectively. From these, 25, 30 and 32, respectively, were successfully identified. Most of the proteins were related to metabolic and cellular processes, followed by response to stimuli, development and regulation of cellular processes. In F-treated groups, PDZK-1, a protein involved in the regulation of renal tubular reabsorption capacity was down-modulated in the kidney of 129P3/J mice. A/J and 129P3/J mice exhibited 11 and 3 exclusive proteins, respectively, regardless of F exposure. In conclusion, proteomic analysis was able to identify proteins potentially involved in metabolic handling of F and water that are differentially expressed or even not expressed in the strains evaluated. This can contribute to understanding the molecular mechanisms underlying genetic susceptibility to dental fluorosis, by indicating key-proteins that should be better addressed in future studies

Genes involved in thoracic exoskeleton formation during the pupal-to-adult molt in a social insect model, Apis mellifera

Background: The insect exoskeleton provides shape, waterproofing, and locomotion via attached somatic muscles. The exoskeleton is renewed during molting, a process regulated by ecdysteroid hormones. The holometabolous pupa transforms into an adult during the imaginal molt, when the epidermis synthe3sizes the definitive exoskeleton that then differentiates progressively. An important issue in insect development concerns how the exoskeletal regions are constructed to provide their morphological, physiological and mechanical functions. We used whole-genome oligonucleotide microarrays to screen for genes involved in exoskeletal formation in the honeybee thoracic dorsum. Our analysis included three sampling times during the pupal-to-adult molt, i.e., before, during and after the ecdysteroid-induced apolysis that triggers synthesis of the adult exoskeleton. Results: Gene ontology annotation based on orthologous relationships with Drosophila melanogaster genes placed the honeybee differentially expressed genes (DEGs) into distinct categories of Biological Process and Molecular Function, depending on developmental time, revealing the functional elements required for adult exoskeleton formation. Of the 1,253 unique DEGs, 547 were upregulated in the thoracic dorsum after apolysis, suggesting induction by the ecdysteroid pulse. The upregulated gene set included 20 of the 47 cuticular protein (CP) genes that were previously identified in the honeybee genome, and three novel putative CP genes that do not belong to a known CP family. In situ hybridization showed that two of the novel genes were abundantly expressed in the epidermis during adult exoskeleton formation, strongly implicating them as genuine CP genes. Conserved sequence motifs identified the CP genes as members of the CPR, Tweedle, Apidermin, CPF, CPLCP1 and Analogous-to-Peritrophins families. Furthermore, 28 of the 36 muscle-related DEGs were upregulated during the de novo formation of striated fibers attached to the exoskeleton. A search for cis-regulatory motifs in the 5′-untranslated region of the DEGs revealed potential binding sites for known transcription factors. Construction of a regulatory network showed that various upregulated CP- and muscle-related genes (15 and 21 genes, respectively) share common elements, suggesting co-regulation during thoracic exoskeleton formation. Conclusions: These findings help reveal molecular aspects of rigid thoracic exoskeleton formation during the ecdysteroid-coordinated pupal-to-adult molt in the honeybee.

A proteomic analysis of the functional effects of fatty acids in NIH 3T3 fibroblasts

Previous studies have demonstrated that long chain fatty acids influence fibroblast function at sub-lethal concentrations. This study is the first to assess the effects of oleic, linoleic or palmitic acids on protein expression of fibroblasts, as determined by standard proteomic techniques. The fatty acids were not cytotoxic at the concentration used in this work as assessed by membrane integrity, DNA fragmentation and the MTT assay but significantly increased cell proliferation. Subsequently, a proteomic analysis was performed using two dimensional difference gel electrophoresis (2D-DIGE) and MS based identification. Cells treated with 50 μM oleic, linoleic or palmitic acid for 24 h were associated with 24, 22, 16 spots differentially expressed, respectively. Among the identified proteins, α-enolase and far upstream element binding protein 1 (FBP-1) are of importance due to their function in fibroblast-associated diseases. However, modulation of α-enolase and FBP-1 expression by fatty acids was not validated by the Western blot technique.

High-throughput sequencing of small RNA transcriptomes reveals critical biological features targeted by microRNAs in cell models used for squamous cell cancer research

Abstract Background The implication of post-transcriptional regulation by microRNAs in molecular mechanisms underlying cancer disease is well documented. However, their interference at the cellular level is not fully explored. Functional in vitro studies are fundamental for the comprehension of their role; nevertheless results are highly dependable on the adopted cellular model. Next generation small RNA transcriptomic sequencing data of a tumor cell line and keratinocytes derived from primary culture was generated in order to characterize the microRNA content of these systems, thus helping in their understanding. Both constitute cell models for functional studies of microRNAs in head and neck squamous cell carcinoma (HNSCC), a smoking-related cancer. Known microRNAs were quantified and analyzed in the context of gene regulation. New microRNAs were investigated using similarity and structural search, ab initio classification, and prediction of the location of mature microRNAs within would-be precursor sequences. Results were compared with small RNA transcriptomic sequences from HNSCC samples in order to access the applicability of these cell models for cancer phenotype comprehension and for novel molecule discovery. Results Ten miRNAs represented over 70% of the mature molecules present in each of the cell types. The most expressed molecules were miR-21, miR-24 and miR-205, Accordingly; miR-21 and miR-205 have been previously shown to play a role in epithelial cell biology. Although miR-21 has been implicated in cancer development, and evaluated as a biomarker in HNSCC progression, no significant expression differences were seen between cell types. We demonstrate that differentially expressed mature miRNAs target cell differentiation and apoptosis related biological processes, indicating that they might represent, with acceptable accuracy, the genetic context from which they derive. Most miRNAs identified in the cancer cell line and in keratinocytes were present in tumor samples and cancer-free samples, respectively, with miR-21, miR-24 and miR-205 still among the most prevalent molecules at all instances. Thirteen miRNA-like structures, containing reads identified by the deep sequencing, were predicted from putative miRNA precursor sequences. Strong evidences suggest that one of them could be a new miRNA. This molecule was mostly expressed in the tumor cell line and HNSCC samples indicating a possible biological function in cancer. Conclusions Critical biological features of cells must be fully understood before they can be chosen as models for functional studies. Expression levels of miRNAs relate to cell type and tissue context. This study provides insights on miRNA content of two cell models used for cancer research. Pathways commonly deregulated in HNSCC might be targeted by most expressed and also by differentially expressed miRNAs. Results indicate that the use of cell models for cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. Additionally, one new miRNA-like molecule with a potential role in cancer was identified in the cell lines and clinical samples.