A comparison of methods for differential expression analysis of RNA-seq data.

BACKGROUND: Finding genes that are differentially expressed between conditions is an integral part of understanding the molecular basis of phenotypic variation. In the past decades, DNA microarrays have been used extensively to quantify the abundance of mRNA corresponding to different genes, and more recently high-throughput sequencing of cDNA (RNA-seq) has emerged as a powerful competitor. As the cost of sequencing decreases, it is conceivable that the use of RNA-seq for differential expression analysis will increase rapidly. To exploit the possibilities and address the challenges posed by this relatively new type of data, a number of software packages have been developed especially for differential expression analysis of RNA-seq data. RESULTS: We conducted an extensive comparison of eleven methods for differential expression analysis of RNA-seq data. All methods are freely available within the R framework and take as input a matrix of counts, i.e. the number of reads mapping to each genomic feature of interest in each of a number of samples. We evaluate the methods based on both simulated data and real RNA-seq data. CONCLUSIONS: Very small sample sizes, which are still common in RNA-seq experiments, impose problems for all evaluated methods and any results obtained under such conditions should be interpreted with caution. For larger sample sizes, the methods combining a variance-stabilizing transformation with the 'limma' method for differential expression analysis perform well under many different conditions, as does the nonparametric SAMseq method.

Temporal analysis of prostaglandin F2 alpha receptor, caspase 3, and cyclooxygenase 2 messenger RNA expression and prostaglandin F2 alpha receptor and cyclooxygenase 2 protein expression in endometrial tissue from multiparous Nelore (Bos taurus indicus) cows treated with cloprostenol sodium during puerperium

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

Selection of suitable housekeeping genes for expression analysis in glioblastoma using quantitative RT-PCR

Background: Considering the broad variation in the expression of housekeeping genes among tissues and experimental situations, studies using quantitative RT-PCR require strict definition of adequate endogenous controls. For glioblastoma, the most common type of tumor in the central nervous system, there was no previous report regarding this issue. Results: Here we show that amongst seven frequently used housekeeping genes TBP and HPRT1 are adequate references for glioblastoma gene expression analysis. Evaluation of the expression levels of 12 target genes utilizing different endogenous controls revealed that the normalization method applied might introduce errors in the estimation of relative quantities. Genes presenting expression levels which do not significantly differ between tumor and normal tissues can be considered either increased or decreased if unsuitable reference genes are applied. Most importantly, genes showing significant differences in expression levels between tumor and normal tissues can be missed. We also demonstrated that the Holliday Junction Recognizing Protein, a novel DNA repair protein over expressed in lung cancer, is extremely over-expressed in glioblastoma, with a median change of about 134 fold. Conclusion: Altogether, our data show the relevance of previous validation of candidate control genes for each experimental model and indicate TBP plus HPRT1 as suitable references for studies on glioblastoma gene expression.

Expression analysis of putative vitellogenin and lipophorin receptors in honey bee (Apis mellifera L.) queens and workers

Two members of the low density lipoprotein receptor (LDLR) family were identified as putative orthologs for a vitellogenin receptor (Amvgr) and a lipophorin receptor (Amlpr) in the Apis mellifera genome. Both receptor sequences have the structural motifs characteristic of LDLR family members and show a high degree of similarity with sequences of other insects. RT-PCR analysis of Amvgr and Amlpr expression detected the presence of both transcripts in different tissues of adult female (ovary, fat body, midgut, head and specifically hypopharyngeal gland), as well as in embryos. In the head RNA samples we found two variant forms of AmLpR: a full length one and a shorter one lacking 29 amino acids in the O-linked sugar domain. In ovaries the expression levels of the two honey bee LDLR members showed opposing trends: whereas Amvgr expression was upregulated as the ovaries became activated, Amlpr transcript levels gradually declined. In situ hybridization analysis performed on ovaries detected Amvgr mRNA exclusively in germ line cells and corroborated the qPCR results showing an increase in Amvgr gene expression concomitant with follicle growth. (C) 2008 Elsevier Ltd. All rights reserved.

RNA-seq analysis of the Quercus suber root response to drought

Dissertação de mestrado em Plant Molecular Biology, Biotechnology and Bioentrepreneurship

Small nucleolar RNA expression profiling identifies potential prognostic markers in peripheral T-cell lymphoma.

Peripheral T-cell lymphoma (PTCL) is a rare, heterogeneous type of non-Hodgkin lymphoma (NHL) that, in general, is associated with a poor clinical outcome. Therefore, a current major challenge is the discovery of new prognostic tools for this disease. In the present study, a cohort of 122 patients with PTCL was collected from a multicentric T-cell lymphoma consortium (TENOMIC). We analyzed the expression of 80 small nucleolar RNAs (snoRNAs) using high-throughput quantitative PCR. We demonstrate that snoRNA expression analysis may be useful in both the diagnosis of some subtypes of PTCL and the prognostication of both PTCL-not otherwise specified (PTCL-NOS; n = 26) and angio-immunoblastic T-cell lymphoma (AITL; n = 46) patients treated with chemotherapy. Like miRNAs, snoRNAs are globally down-regulated in tumor cells compared with their normal counterparts. In the present study, the snoRNA signature was robust enough to differentiate anaplastic large cell lymphoma (n = 32) from other PTCLs. For PTCL-NOS and AITL, we obtained 2 distinct prognostic signatures with a reduced set of 3 genes. Of particular interest was the prognostic value of HBII-239 snoRNA, which was significantly over-expressed in cases of AITL and PTCL-NOS that had favorable outcomes. Our results suggest that snoRNA expression profiles may have a diagnostic and prognostic significance for PTCL, offering new tools for patient care and follow-up.

Identification of reference genes for expression analysis by real-time quantitative PCR in drought-stressed soybean

The objective of this work was to validate, by quantitative PCR in real time (RT-qPCR), genes to be used as reference in studies of gene expression in soybean in drought-stressed trials. Four genes commonly used in soybean were evaluated: Gmβ-actin, GmGAPDH, GmLectin and GmRNAr18S. Total RNA was extracted from six samples: three from roots in a hydroponic system with different drought intensities (0, 25, 50, 75 and 100 minutes of water stress), and three from leaves of plants grown in sand with different soil moistures (15, 5 and 2.5% gravimetric humidity). The raw cycle threshold (Ct) data were analyzed, and the efficiency of each primer was calculated for an overall analysis of the Ct range among the different samples. The GeNorm application was used to evaluate the best reference gene, according to its stability. The GmGAPDH was the least stable gene, with the highest mean values of expression stability (M), and the most stable genes, with the lowest M values, were the Gmβ-actin and GmRNAr18S, when both root and leaves samples were tested. These genes can be used in RT-qPCR as reference gene for expression analysis.

Differential gene expression analysis of iodide-treated rat thyroid follicular cell line PCCl3

The inhibitory effect of supraphysiological iodide concentrations on thyroid hormone synthesis (Wolff - Chaikoff effect) and on thyrocyte proliferation is largely known as iodine autoregulation. However, the molecular mechanisms by which iodide modulates thyroid function remain unclear. In this paper, we analyze the transcriptome profile of the rat follicular cell lineage PCCl3 under untreated and treated conditions with 10 (- 3) M sodium iodide (NaI). Serial analysis of gene expression (SAGE) revealed 84 transcripts differentially expressed in response to iodide (p <= 0.001). We also showed that iodide excess inhibits the expression of essential genes for thyroid differentiation: Tshr, Nis, Tg, and Tpo. Relative expression of 14 of 20 transcripts selected by SAGE was confirmed by real-time PCR. Considering the key role of iodide organification in thyroid physiology, we also observed that both the oxidized form of iodide and iodide per se are responsible for gene expression modulation in response to iodide excess. (c) 2008 Elsevier Inc. All rights reserved.

Global Gene Expression Analysis during Sporulation of the Aquatic Fungus Blastocladiella emersonii

The Blastocladiella emersonii life cycle presents a number of drastic biochemical and morphological changes, mainly during two cell differentiation stages: germination and sporulation. To investigate the transcriptional changes taking place during the sporulation phase, which culminates with the production of the zoospores, motile cells responsible for the dispersal of the fungus, microarray experiments were performed. Among the 3,773 distinct genes investigated, a total of 1,207 were classified as differentially expressed, relative to time zero of sporulation, at at least one of the time points analyzed. These results indicate that accurate transcriptional control takes place during sporulation, as well as indicating the necessity for distinct molecular functions throughout this differentiation process. The main functional categories overrepresented among upregulated genes were those involving the microtubule, the cytoskeleton, signal transduction involving Ca(2+), and chromosome organization. On the other hand, protein biosynthesis, central carbon metabolism, and protein degradation were the most represented functional categories among downregulated genes. Gene expression changes were also analyzed in cells sporulating in the presence of subinhibitory concentrations of glucose or tryptophan. Data obtained revealed overexpression of microtubule and cytoskeleton transcripts in the presence of glucose, probably causing the shape and motility problems observed in the zoospores produced under this condition. In contrast, the presence of tryptophan during sporulation led to upregulation of genes involved in oxidative stress, proteolysis, and protein folding. These results indicate that distinct physiological pathways are involved in the inhibition of sporulation due to these two classes of nutrient sources.

Sequencing and expression analysis of hepcidin mRNA in donkey (Equus asinus) liver

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

Gene expression analysis reveals novel altered genes in nasal polyps

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

Identification of reference genes for expression analysis by real-time quantitative PCR in drought-stressed soybean

O objetivo deste trabalho foi validar, pela técnica de PCR quantitativo em tempo real (RT-qPCR) genes para serem utilizados como referência em estudos de expressão gênica em soja, em ensaios de estresse hídrico. Foram avaliados quatro genes comumente utilizados em soja: Gmβ-actin, GmGAPDH, GmLectin e GmRNAr18S. O RNA total foi extraído de seis amostras: três amostras de raízes em sistema de hidroponia com diferentes intensidades de déficit hídrico (0, 25, 50, 75 e 100 minutos de estresse hídrico), e três amostras de folhas de plantas cultivadas em areia com diferentes umidades do solo (15, 5 e 2,5% de umidade gravimétrica). Os dados brutos do intervalo cycle threshold (Ct) foram analisados, e a eficiência de cada iniciador foi calculada para uma analise da Ct entre as diferentes amostras. A aplicação do programa GeNorm foi utilizada para a avaliação dos melhores genes de referência, de acordo com a estabilidade. O GmGAPDH foi o gene menos estável, com o maior valor médio de estabilidade de expressão (M), e os genes mais estáveis, com menor valor de M, foram o Gmβ-actin e GmRNAr18S, tanto nas amostras de raízes como nas de folhas. Estes genes podem ser usados em RT-qPCR como gens de referência para análises de expressão gênica.

Epigenetic and expression analysis of TRAIL-R2 and BCL2: on the TRAIL to knowledge of apoptosis in ovarian tumors

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Evaluation of Toll-Like Receptor 2 and 4 RNA Expression and the Cytokine Profile in Postmenopausal Women with Metabolic Syndrome

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

rpb2 is a reliable reference gene for quantitative gene expression analysis in the dermatophyte Trichophyton rubrum

The selection of reference genes used for data normalization to quantify gene expression by real-time PCR amplifications (qRT-PCR) is crucial for the accuracy of this technique. In spite of this, little information regarding such genes for qRT-PCR is available for gene expression analyses in pathogenic fungi. Thus, we investigated the suitability of eight candidate reference genes in isolates of the human dermatophyte Trichophyton rubrum subjected to several environmental challenges, such as drug exposure, interaction with human nail and skin, and heat stress. The stability of these genes was determined by geNorm, NormFinder and Best-Keeper programs. The gene with the most stable expression in the majority of the conditions tested was rpb2 (DNA-dependent RNA polymerase II), which was validated in three T. rubrum strains. Moreover, the combination of rpb2 and chs1 (chitin synthase) genes provided for the most reliable qRT-PCR data normalization in T. rubrum under a broad range of biological conditions. To the best of our knowledge this is the first report on the selection of reference genes for qRT-PCR data normalization in dermatophytes and the results of these studies should permit further analysis of gene expression under several experimental conditions, with improved accuracy and reliability.