Expression of target and reference genes in Daphnia magna exposed to ibuprofen

Background: Transcriptomic techniques are now being applied in ecotoxicology and toxicology to measure the impact of stressors and develop understanding of mechanisms of toxicity. Microarray technology in particular offers the potential to measure thousands of gene responses simultaneously. However, it is important that microarrays responses should be validated, at least initially, using real-time quantitative polymerase chain reaction (QPCR). The accurate measurement of target gene expression requires normalisation to an invariant internal control e. g., total RNA or reference genes. Reference genes are preferable, as they control for variation inherent in the cDNA synthesis and PCR. However, reference gene expression can vary between tissues and experimental conditions, which makes it crucial to validate them prior to application. Results: We evaluated 10 candidate reference genes for QPCR in Daphnia magna following a 24 h exposure to the non-steroidal anti-inflammatory drug (NSAID) ibuprofen (IB) at 0, 20, 40 and 80 mg IB l(-1). Six of the 10 candidates appeared suitable for use as reference genes. As a robust approach, we used a combination normalisation factor (NF), calculated using the geNorm application, based on the geometric mean of three selected reference genes: glyceraldehyde-3-phosphate dehydrogenase, ubiquitin conjugating enzyme and actin. The effects of normalisation are illustrated using as target gene leukotriene B4 12-hydroxydehydrogenase (Ltb4dh), which was upregulated following 24 h exposure to 63-81 mg IB l(-1). Conclusions: As anticipated, use of the NF clarified the response of Ltb4dh in daphnids exposed to sublethal levels of ibuprofen. Our findings emphasise the importance in toxicogenomics of finding and applying invariant internal QPCR control(s) relevant to the study conditions.

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.

Selection of endogenous genes for gene expression studies in Eucalyptus under biotic (Puccinia psidii) and abiotic (acibenzolar-S-methyl) stresses using RT-qPCR

Background: Rust caused by Puccinia psidii Winter has been limiting for the establishment of new Eucalyptus plantations, as well as for resprouting of susceptible genetic materials. Identifying host genes involved in defense responses is important to elucidate resistance mechanisms. Reverse transcription-quantitative PCR is the most common method of mRNA quantitation for gene expression analysis. This method generally employs a reference gene as an internal control to normalize results. A good endogenous control transcript shows minimal variation due to experimental conditions. Findings. We analyzed the expression of 13 genes to identify transcripts with minimal variation in leaves of 60-day-old clonal seedlings of two Eucalyptus clones (rust-resistant and susceptible) subjected to biotic (P. psidii) and abiotic (acibenzolar-S-methyl, ASM) stresses. Conclusions. For tissue samples of clones that did not receive any stimulus, a combination of the eEF2 and EglDH genes was the best control for normalization. When pathogen-inoculated and uninoculated plant samples were compared, eEF2 and UBQ together were more appropriate as normalizers. In ASM-treated and untreated leaves of both clones, transcripts of the CYP and elF4B genes combined were the ones with minimal variation. Finally, when comparing expression in both clones for ASM-treated leaves, P. psidii-inoculated leaves, ASM-treated plus P. psidii-inoculated leaves, and their respective controls, the genes with the most stable expression were EgIDH and UBQ. The chitinase gene, which is highly expressed in studies on plant resistance to phytopathogens, was used to confirm variation in gene expression due to the treatments. © 2010 Laia et al; licensee BioMed Central Ltd.

Validation of a reference control for an SYBR-Green fluorescence assay-based real-time PCR for detection of bovine herpesvirus 5 in experimentally exposed bovine embryos

The objective of this study was to optimize an internal control to improve SYBR-Green-based qPCR to amplify/detect the BoHV-5 US9 gene in bovine embryos produced invitro and experimentally exposed to the virus. We designed an SYBR-Green-based binding assay that is quick to perform, reliable, easily optimized and compares well with the published assay. Herein we demonstrated its general applicability to detect BoHV-5 US9 gene in bovine embryos produced invitro experimentally exposed to BoHV-5. In order to validate the assay, three different reference genes were tested; and the histone 2a gene was shown to be the most adequate for normalizing the qPCR reaction, by considering melting and standard curves ( p<0.05). On the other hand, no differences were found in the development of bovine embryos invitro whether they were exposed to BoHV-5 reference and field strains comparing to unexposed embryos. The developed qPCR assay may have important field applications as it provides an accurate BoHV-5 US9 gene detection using a proven reference gene and is considerably less expensive than the TaqMan qPCR currently employed in sanitary programs. © 2013 Elsevier Ltd.

Caracterização do gene ERG11 envolvido na resistência de Cryptococcus neoformans e Cryptococcus gatti as drogas azólicas e identificação de novos alvos a partir de produtos naturais

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

RNA degradation differentially affects quantitative mRNA measurements of endogenous reference genes in human placenta

It has been highlighted that RNA quality and appropriate reference gene selection is crucial for the interpretation of RT-qPCR results in human placental samples. In this context we investigated the effect of RNA degradation on the mRNA abundance of seven frequently used reference genes in 119 human placental samples. Combining RNA integrity measurements, RT-qPCR analysis and mathematical modeling we found major differences regarding the effect of RNA degradation on the measured expression levels between the different reference genes. Furthermore, we demonstrated that a modified RNA extraction method significantly improved RNA quality and consequently increased transcript levels of all reference genes.

Detection of exon deletions within an entire gene (CFTR) by relative quantification on the LightCycler

BACKGROUND: Cystic fibrosis (CF) is associated with at least 1 pathogen point sequence variant on each CFTR allele. Some symptomatic patients, however, have only 1 detectable pathogen sequence variant and carry, on the other allele, a large deletion that is not detected by conventional screening methods. METHODS: For relative quantitative real-time PCR detection of large deletions in the CFTR gene, we designed DNA-specific primers for each exon of the gene and primers for a reference gene (beta2-microglobulin). For PCR we used a LightCycler system (Roche) and calculated the gene-dosage ratio of CFTR to beta2-microglobulin. We tested the method by screening all 27 exons in 3 healthy individuals and 2 patients with only 1 pathogen sequence variant. We then performed specific deletion screenings in 10 CF patients with known large deletions and a blinded analysis in which we screened 24 individuals for large deletions by testing 8 of 27 exons. RESULTS: None of the ratios for control samples were false positive (for deletions or duplications); moreover, for all samples from patients with known large deletions, the calculated ratios for deleted exons were close to 0.5. In addition, the results from the blinded analysis demonstrated that our method can also be used for the screening of single individuals. CONCLUSIONS: The LightCycler assay allows reliable and rapid screening for large deletions in the CFTR gene and detects the copy number of all 27 exons.

Tissue-specific selection of reference genes is required for expression studies in the mouse model of maternal protein undernutrition

Suboptimal maternal nutrition during gestation results in the establishment of long-term phenotypic changes and an increased disease risk in the offspring. To elucidate how such environmental sensitivity results in physiological outcomes, the molecular characterisation of these offspring has become the focus of many studies. However, the likely modification of key cellular processes such as metabolism in response to maternal undernutrition raises the question of whether the genes typically used as reference constants in gene expression studies are suitable controls. Using a mouse model of maternal protein undernutrition, we have investigated the stability of seven commonly used reference genes (18s, Hprt1, Pgk1, Ppib, Sdha, Tbp and Tuba1) in a variety of offspring tissues including liver, kidney, heart, retro-peritoneal and inter-scapular fat, extra-embryonic placenta and yolk sac, as well as in the preimplantation blastocyst and blastocyst-derived embryonic stem cells. We find that although the selected reference genes are all highly stable within this system, they show tissue, treatment and sex-specific variation. Furthermore, software-based selection approaches rank reference genes differently and do not always identify genes which differ between conditions. Therefore, we recommend that reference gene selection for gene expression studies should be thoroughly validated for each tissue of interest. © 2011 Elsevier Inc.

Identification of superior reference genes for data normalisation of expression studies via quantitative PCR in hybrid roses (Rosa hybrida)

Background: Gene expression studies are a prerequisite for understanding the biological function of genes. Because of its high sensitivity and easy use, quantitative PCR (qPCR) has become the gold standard for gene expression quantification. To normalise qPCR measurements between samples, the most prominent technique is the use of stably expressed endogenous control genes, the so called reference genes. However, recent studies show there is no universal reference gene for all biological questions. Roses are important ornamental plants for which there has been no evaluation of useful reference genes for gene expression studies. Results: We used three different algorithms (BestKeeper, geNorm and NormFinder) to validate the expression stability of nine candidate reference genes in different rose tissues from three different genotypes of Rosa hybrida and in leaves treated with various stress factors. The candidate genes comprised the classical "housekeeping genes" (Actin, EF-1α, GAPDH, Tubulin and Ubiquitin), and genes showing stable expression in studies in Arabidopsis (PP2A, SAND, TIP and UBC). The programs identified no single gene that showed stable expression under all of the conditions tested, and the individual rankings of the genes differed between the algorithms. Nevertheless the new candidate genes, specifically, PP2A and UBC, were ranked higher as compared to the other traditional reference genes. In general, Tubulin showed the most variable expression and should be avoided as a reference gene. Conclusions: Reference genes evaluated as suitable in experiments with Arabidopsis thaliana were stably expressed in roses under various experimental conditions. In most cases, these genes outperformed conventional reference genes, such as EF1-α and Tubulin. We identified PP2A, SAND and UBC as suitable reference genes, which in different combinations may be used for normalisation in expression analyses via qPCR for different rose tissues and stress treatments. However, the vast genetic variation found within the genus Rosa, including differences in ploidy levels, might also influence expression stability of reference genes, so that future research should also consider different genotypes and ploidy levels.

Expression of apoptosis and survival genes in human memory T cell populations

RESUME Introduction: Les cellules T mémoires humaines sont classées en trois sous-populations sur la base de l'expression d'un marqueur de surface cellulaire, CD45RA, et du récepteur aux chimiokines, CCR7. Ces sous-populations, nommées cellules mémoires centrales (TcM), mémoires effectrices (TEM) et mémoires effectrices terminales (ITEM), ont des rôles fonctionnels distincts, ainsi que des capacités de prolifération et de régénération différentes. Cependant, la génération de ces différences reste encore mal comprise et on ignore les mécanismes moléculaires impliqués. Matériaux et Méthodes: Des cellules mononucléaires humaines du sang périphérique ont été séparées par cytométrie de flux selon leur expression de CD4, CD8, CD45RA et CCR7 en sous-populations de cellules CD4+ ou CD8+ naïves, TcM, TEM ou ITEM. Dans chacune de ces sous-populations, 14 gènes impliqués dans l'apoptose, la survie ou la capacité proliférative des cellules T ont été quantifiés par RT-PCR en temps réel, relativement à l'expression d'un gène de référence endogène. L'ARN provenant de 450 cellules T a été utilisé par gène et par sous-population. Les gènes analysés (cibles) comprenaient des gènes de survie (BAFF, APRIL, BAFF-R, BCMA, TACI, IL-15Rα, IL-7Rα), des gènes anti-apoptotiques (Bcl-2, BclxL, FLIP), des gènes pro-apoptotiques (Bad, Bax, Fast) et le gène anti-prolifératif, Tob. A l'aide de la méthode comparative delta-delta-CT, le taux d'expression des gènes cibles de chaque sous-population des cellules T mémoires CD4+ et CD8+, à été comparée à leur taux d'expression dans les cellules T naïves CD4+ et CD8+. Résultats: Dans les cellules CD8+, les gènes pro-apoptotiques Bax et Fast étaient surexprimés dans toutes les sous-populations mémoires, tandis que l'expression des facteurs anti-apoptotiques et de survie comme Bcl-2, APRIL et BAFF-R, étaient diminués. Ces deux tendances étaient particulièrement accentuées dans les sous-groupes des cellules mémoires TEM et TTEM. A noter que malgré le fait que leur expression était également diminuée dans les autres cellules mémoires, le facteur de survie IL-7Ra, était sélectivement surexprimé dans la sous-population de cellules TcM et l'expression d'IL-15Ra était sélectivement augmentée dans les TEM. Dans les cellules CD4+, le taux d'expression des gènes analysés était plus variable entre les sujets étudiés que dans les cellules CD8+, ne permettant pas de définir un profil d'expression spécifique. L'expression du gène de survie BAFF par contre, a été significativement augmentée dans toutes les sous-populations mémoire CD4+. Il en va de même pour l'expression d' APRIL et de BAFF-R, bien que dans moindre degré. A remarquer que l'expression du facteur anti-apoptotique Fast a été observée uniquement dans la souspopulation des TTEM. Discussion et Conclusions: Cette étude montre une nette différence entre les cellules CD8+ et CD4+, en ce qui concerne les profils d'expression des gènes impliqués dans la survie et l'apoptose des cellules T mémoires. Ceci pourrait impliquer une régulation cellulaire homéostatique distincte dans ces deux compartiments de cellules T mémoires. Dans les cellules CD8+ l'expression d'un nombre de gènes impliqués dans la survie et la protection de l'apoptose semblerait être diminuée dans les populations TEM et TTEM en comparaison à celle des sous-populations naïves et TEM, tandis que l'expression des gènes pro-apoptotiques semblerait être augmentée. Comme ceci paraît être plus accentué dans les TTEM, cela pourrait indiquer une plus grande disposition à l'apopotose dans les populations CCR7- (effectrices) et une perte de survie parallèlement à l'acquisition de capacités effectrices. Ceci parlerait en faveur d'un modèle de différentiation linéaire dans les cellules CD8+. De plus, l'augmentation sélective de l'expression d'IL-7Ra observée dans le sous-groupe de cellules mémoires TEM, et d'IL-15Ra dans celui des TEM, pourrait indiquer un moyen de sélection pour des réponses immunitaires mémoires à long terme par une réponse distincte à ces cytokines. Dans les cellules CD4+ par contre, aucun profil d'expression n'a pu être déterminé; les résultats suggèrent même une résistance relative à l'apoptose de la part des cellules mémoires. Ceci pourrait favoriser l'existence d'un modèle de différentiation plus flexible avec des possibilités d'interaction multiples. Ainsi, la surexpression sélective de BAFF, APRIL et BAFF-R dans les sous-populations individuelles des cellules mémoires pourrait être un indice de l'interaction de ces sous-groupes avec des cellules B. ABSTRACT Introduction: Based on their surface expression of the CD45 isoform and of the CCR7 chemokine receptor, memory T cells have been divided into the following three subsets: central memory (TAM), effector memory (TEM) and terminal effector memory (ITEM). Distinct functional roles and different proliferative and regenerative capacities have been attributed to each one of these subpopulations. The molecular mechanisms underlying these differences; however, remain poorly understood. Materials and Methods: According to their expression of CD4, CD8, CD45RA and CCR7, human peripheral blood mononuclear cells were sorted by flow-cytometry into CD4+ or CD8+ naïve, TAM, TEM and ITEM subsets. Using real-time PCR, the expression of 14 genes known to be involved in apoptotis, survival or proliferation of T cells was quantified separately in each individual subset, relative to an endogenous reference gene. The RNA equivalent of 450 T cells was used for each gene and subset. The target gene panel included the survival genes BAFF, APRIL, BAFF-R, BCMA, TACI, IL-15Rα and IL-7Rα, the anti-apoptotic genes Bcl2, Bcl-xL and FLIP, the pro-apoptotic genes Bad, Bax and Fast, as well as the antiproliferative gene Tob. Using the comparative CT-method, the expression of the target genes in the three memory T cell subsets of both CD4+ and CD8+ T cell populations was compared to their expression in the naïve T cells. Results: In CD8+ cells, the pro-apoptotic factors Bax and Fast were found to be upregulated in all memory T cell subsets, whereas the survival and anti-apoptotic factors Bcl-2, APRIL and BAFF-R were downregulated. These tendencies were most accentuated in TEM and TTEM subsets. Even though the survival factor IL-7Rα was also downregulated in these subsets, interestingly, it was selectively upregulated in the CD8+ TAM subset. Similarly, IL-15Rαexpression was shown to be selectively upregulated in the CD8+ TEM subset. In CD4+ cells, the expression levels of the analyzed genes showed a greater inter-individual variability than in CD8+ cells, thus suggesting the absence of any particular expression pattern for CD4+ memory T cells. However, the survival factor BAFF was found to be significantly upregulated in all CD4+ memory T cell subsets, as was also the expression of APRIL and BAFF-R, although to a lesser extent. Furthermore, it was noted that the pro-apoptotic gene Fast was only expressed in the TTEM CD4+ subset. Discussion and Conclusions: Genes involved in apoptosis and survival in human memory T cells have been shown to be expressed differently in CD8+ cells as compared to CD4+ cells, suggesting a distinct regulation of cell homeostasis in these two memory T cell compartments. The present study suggests that, in CD8+ T cells, the expression of various survival and antiapoptotic genes is downregulated in TEM and TTEM subsets, while the expression of proapoptotic genes is upregulated in comparison to the naïve and the TAM populations. These characteristics, potentially translating to a greater susceptibility to apoptosis in the CCR7- (effector) memory populations, are accentuated in the TTEM population, suggesting a loss of survival in parallel to the acquisition of effector capacities. This speaks in favour of a linear differentiation model in CD8+ T memory cells. Moreover, the observed selectively increased expression of IL-7Rα in CD8+ TAM cells - as that of IL-15Rα in CD8+ TEM cells -suggest that differential responsiveness to cytokines could confer a selection bias for distinct long-term memory cell responses. Relative to the results for CD8+ T cells, those for CD4+ T cells seem to indicate a certain resistance of the memory subsets to apoptosis, suggesting the possibility of a more flexible differentiation model with multiple checkpoints and potential interaction of CD4+ memory cells with other cells. Thus, the selective upregulation of BAFF, APRIL and BAFF-R in individual memory subsets could imply an interaction of these subsets with B cells.

Molecular Evolution of Metazoan Hypoxia-inducing Factors

Most metazoans rely on aerobic energy production, which is dependent on adequate oxygen supply. In the case of reduced oxygen supply (hypoxia), the most profound changes in gene expression are mediated by transcription factors named hypoxia-inducible factors (HIF alpha). These proteins are post-translationally regulated by prolyl-4-hydroxylase (PHD) enzymes that are direct “sensors” of cellular oxygen levels. This thesis examines the molecular evolution of metazoan HIF systems. In early metazoans the HIF system emerged from pre-existing PHD oxygen sensors and early bHLH-PAS transcription factors. In invertebrates our analysis revealed an unexpected diversity of PHD genes and HIF alpha sequence characteristics. An early branching vertebrate, the epaulette shark (Hemiscyllium ocellatum) was chosen for sequencing and hypoxia preconditioning studies of HIF alpha and PHD genes. As no quantitative PCR reference genes were available, this thesis includes the first study of reference genes in cartilaginous fish species. Applying multiple statistical analysis we also discoveredthat commonly used reference gene software may perform poorly with some data sets. Novel reference genes allowed accurate measurements of the mRNAlevels of the studied target genes. Cartilaginous fishes have three genomic duplicates of both HIF alpha and PHD genes like mammals and teleost fishes. Combining functional divergence and selection analyses it was possible to describe how sequence changes in both HIF alpha and PHD duplicates may have contributed to the differential oxygen sensitivityof HIF alphas. Additionally, novel teleost HIF-1 alpha sequences were produced and used to reveal the molecular evolution of HIF-1 alpha in this lineage rich with hypoxia tolerant species.

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

The hypoferremia that is observed during systemic inflammatory processes is mediated by hepcidin, which is a peptide that is mainly synthesized in the livers of several mammalian species. Hepcidin plays a key role in iron metabolism and in the innate immune system. It's up-regulation is particularly useful during acute inflammation, and it restricts the iron availability that is necessary for the growth of pathogenic microorganisms. In this study, the hepcidin mRNA of Equus asinus has been characterized, and the expression of donkey hepcidin in the liver has been determined. The donkey hepcidin sequence has an open reading frame (ORF) of 261 nucleotides, and the deduced corresponding protein sequence has 86 amino acids. The amino acid sequence of donkey hepcidin was most homologous to Equus caballus (98%). The mature donkey hepcidin sequence (25 amino acids) was 100% homologous to the equine mature hepcidin and has eight conserved cysteine residues that are found in all of the investigated hepcidin sequences. The expression profile of donkey hepcidin in the liver was high and was similar to the reference gene expression. The donkey hepcidin sequence was deposited in GenBankTM (HQ902884) and may be useful for additional studies on iron metabolism and the inflammatory process in this species.

Expression of TERT in precancerous gastric lesions compared to gastric cancer

The objective of this study was to determine the levels of TERT mRNA and TERT protein expression in stomach precancerous lesions such as intestinal metaplasia (IM) and gastric ulcer (GU) and compare them to gastric cancer (GC). Real-time PCR was performed to detect TERT mRNA expression levels in 35 biopsies of IM, 30 of GU, and 22 of GC and their respective normal mucosas. TERT protein was detected by immunohistochemistry in 68 samples, 34 of IM, 23 of GU, and 11 of GC. Increased TERT mRNA expression levels were observed in a significant number of cases, i.e., 46% of IM, 50% of GU, and 79% of GC. The relative mean level of TERT mRNA after normalization with the β-actin reference gene and comparison with the respective adjacent normal mucosa was slightly increased in the IM and GU groups, 2.008 ± 2.605 and 2.730 ± 4.120, respectively, but high TERT mRNA expression was observed in the GC group (17.271 ± 33.852). However, there were no statistically significant differences between the three groups. TERT protein-positive immunostaining was observed in 38% of IM, 39% of GU, and 55% of GC. No association of TERT mRNA and protein expression with Helicobacter pylori infection or other clinicopathological variables was demonstrable, except for the incomplete type vs the complete type of IM. This study confirms previous data of the high expression of both TERT mRNA and protein in gastric cancer and also demonstrates this type of changed expression in IM and GU, thus suggesting that TERT expression may be deregulated in precursor lesions that participate in the early stages of gastric carcinogenesis.

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

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

Identification of suitable internal control genes for expression studies in Coffea arabica under different experimental conditions

Background: Quantitative data from gene expression experiments are often normalized by transcription levels of reference or housekeeping genes. An inherent assumption for their use is that the expression of these genes is highly uniform in living organisms during various phases of development, in different cell types and under diverse environmental conditions. To date, the validation of reference genes in plants has received very little attention and suitable reference genes have not been defined for a great number of crop species including Coffea arabica. The aim of the research reported herein was to compare the relative expression of a set of potential reference genes across different types of tissue/organ samples of coffee. We also validated the expression profiles of the selected reference genes at various stages of development and under a specific biotic stress.Results: The expression levels of five frequently used housekeeping genes (reference genes), namely alcohol dehydrogenase (adh), 14-3-3, polyubiquitin (poly), beta-actin (actin) and glyceraldehyde-3-phosphate dehydrogenase (gapdh) was assessed by quantitative real-time RT-PCR over a set of five tissue/organ samples (root, stem, leaf, flower, and fruits) of Coffea arabica plants. In addition to these commonly used internal controls, three other genes encoding a cysteine proteinase (cys), a caffeine synthase (ccs) and the 60S ribosomal protein L7 (rpl7) were also tested. Their stability and suitability as reference genes were validated by geNorm, NormFinder and BestKeeper programs. The obtained results revealed significantly variable expression levels of all reference genes analyzed, with the exception of gapdh, which showed no significant changes in expression among the investigated experimental conditions.Conclusion: Our data suggests that the expression of housekeeping genes is not completely stable in coffee. Based on our results, gapdh, followed by 14-3-3 and rpl7 were found to be homogeneously expressed and are therefore adequate for normalization purposes, showing equivalent transcript levels in different tissue/ organ samples. Gapdh is therefore the recommended reference gene for measuring gene expression in Coffea arabica. Its use will enable more accurate and reliable normalization of tissue/organ-specific gene expression studies in this important cherry crop plant.