Geographic differentiation of polymorphism in the Plasmodium falciparum malaria vaccine candidate gene SERA5

SERA5 is regarded as a promising malaria vaccine candidate of the most virulent human malaria parasite Plasmodium falciparum. SERA5 is a 120 kDa abundantly expressed blood-stage protein containing a papain-like protease. Since substantial polymorphism in blood-stage vaccine candidates may potentially limit their efficacy, it is imperative to fully investigate polymorphism of the SERA5 gene (sera5). In this study, we performed evolutionary and population genetic analysis of sera5. The level of inter-species divergence (kS = 0.076) between P. falciparum and Plasmodium reichenowi, a closely related chimpanzee malaria parasite is comparable to that of housekeeping protein genes. A signature of purifying selection was detected in the proenzyme and enzyme domains. Analysis of 445 near full-length P. falciparum sera5 sequences from nine countries in Africa, Southeast Asia, Oceania and South America revealed extensive variations in the number of octamer repeat (OR) and serine repeat (SR) regions as well as substantial level of single nucleotide polymorphism (SNP) in non-repeat regions (2562 bp). Remarkably, a 14 amino acid sequence of SERA5 (amino acids 59-72) that is known to be the in vitro target of parasite growth inhibitory antibodies was found to be perfectly conserved in all 445 worldwide isolates of P. falciparum evaluated. Unlike other major vaccine target antigen genes such as merozoite surface protein-1, apical membrane antigen-1 or circumsporozoite protein, no strong evidence for positive selection was detected for SNPs in the non-repeat regions of sera5. A biased geographical distribution was observed in SNPs as well as in the haplotypes of the sera5 OR and SR regions. In Africa, OR- and SR-haplotypes with low frequency (<5%) and SNPs with minor allele frequency (<5%) were abundant and were mostly continent-specific. Consistently, significant genetic differentiation, assessed by the Wright's fixation index (FST) of inter-population variance in allele frequencies, was detected for SNPs and both OR- and SR-haplotypes among almost all parasite populations. The exception was parasite populations between Tanzania and Ghana, suggesting frequent gene flow in Africa. The present study points to the importance of investigating whether biased geographical distribution for SNPs and repeat variants in the OR and SR regions affect the reactivity of human serum antibodies to variants. (C) 2011 Elsevier Ltd. All rights reserved.

Cellular and humoral immune responses against the Plasmodium vivax MSP-119 malaria vaccine candidate in individuals living in an endemic area in north-eastern Amazon region of Brazil

Abstract Background Plasmodium vivax merozoite surface protein-1 (MSP-1) is an antigen considered to be one of the leading malaria vaccine candidates. PvMSP-1 is highly immunogenic and evidences suggest that it is target for protective immunity against asexual blood stages of malaria parasites. Thus, this study aims to evaluate the acquired cellular and antibody immune responses against PvMSP-1 in individuals naturally exposed to malaria infections in a malaria-endemic area in the north-eastern Amazon region of Brazil. Methods The study was carried out in Paragominas, Pará State, in the Brazilian Amazon. Blood samples were collected from 35 individuals with uncomplicated malaria. Peripheral blood mononuclear cells were isolated and the cellular proliferation and activation was analysed in presence of 19 kDa fragment of MSP-1 (PvMSP-119) and Plasmodium falciparum PSS1 crude antigen. Antibodies IgE, IgM, IgG and IgG subclass and the levels of TNF, IFN-γ and IL-10 were measured by enzyme-linked immunosorbent assay. Results The prevalence of activated CD4+ was greater than CD8+ T cells, in both ex-vivo and in 96 h culture in presence of PvMSP-119 and PSS1 antigen. A low proliferative response against PvMSP-119 and PSS1 crude antigen after 96 h culture was observed. High plasmatic levels of IFN-γ and IL-10 as well as lower TNF levels were also detected in malaria patients. However, in the 96 h supernatant culture, the dynamics of cytokine responses differed from those depicted on plasma assays; in presence of PvMSP-119 stimulus, higher levels of TNF were noted in supernatant 96 h culture of malaria patient’s cells while low levels of IFN-γ and IL-10 were verified. High frequency of malaria patients presenting antibodies against PvMSP-119 was evidenced, regardless class or IgG subclass.PvMSP-119-induced antibodies were predominantly on non-cytophilic subclasses. Conclusions The results presented here shows that PvMSP-119 was able to induce a high cellular activation, leading to production of TNF and emphasizes the high immunogenicity of PvMSP-119 in naturally exposed individuals and, therefore, its potential as a malaria vaccine candidate.

Identification of Leishmania infantum chagasi proteins in urine of patients with visceral leishmaniasis: a promising antigen discovery approach of vaccine candidates

Visceral leishmaniasis (VL) is a serious lethal parasitic disease caused by Leishmania donovani in Asia and by Leishmania infantum chagasi in southern Europe and South America. VL is endemic in 47 countries with an annual incidence estimated to be 500 000 cases. This high incidence is due in part to the lack of an efficacious vaccine. Here, we introduce an innovative approach to directly identify parasite vaccine candidate antigens that are abundantly produced in vivo in humans with VL. We combined RP-HPLC and mass spectrometry and categorized three L. infantum chagasi proteins, presumably produced in spleen, liver and bone marrow lesions and excreted in the patients urine. Specifically, these proteins were the following: Li-isd1 (XP_001467866.1), Li-txn1 (XP_001466642.1) and Li-ntf2 (XP_001463738.1). Initial vaccine validation studies were performed with the rLi-ntf2 protein produced in Escherichia coli mixed with the adjuvant BpMPLA-SE. This formulation stimulated potent Th1 response in BALB/c mice. Compared to control animals, mice immunized with Li-ntf2+ BpMPLA-SE had a marked parasite burden reduction in spleens at 40 days post-challenge with virulent L. infantum chagasi. These results strongly support the proposed antigen discovery strategy of vaccine candidates to VL and opens novel possibilities for vaccine development to other serious infectious diseases.

Plasmodium vivax Duffy binding protein: baseline antibody responses and parasite polymorphisms in a well-consolidated settlement of the Amazon Region

Objective To investigate risk factors associated with the acquisition of antibodies against Plasmodium vivax Duffy binding protein (PvDBP) a leading malaria vaccine candidate in a well-consolidated agricultural settlement of the Brazilian Amazon Region and to determine the sequence diversity of the PvDBP ligand domain (DBPII) within the local malaria parasite population. Methods Demographic, epidemiological and clinical data were collected from 541 volunteers using a structured questionnaire. Malaria parasites were detected by conventional microscopy and PCR, and blood collection was used for antibody assays and molecular characterisation of DBPII. Results The frequency of malaria infection was 7% (6% for P. vivax and 1% for P. falciparum), with malaria cases clustered near mosquito breeding sites. Nearly 50% of settlers had anti-PvDBP IgG antibodies, as detected by enzyme-linked immunosorbent assay (ELISA) with subjects age being the only strong predictor of seropositivity to PvDBP. Unexpectedly, low levels of DBPII diversity were found within the local malaria parasites, suggesting the existence of low gene flow between P. vivax populations, probably due to the relative isolation of the studied settlement. Conclusion The recognition of PvDBP by a significant proportion of the community, associated with low levels of DBPII diversity among local P. vivax, reinforces the variety of malaria transmission patterns in communities from frontier settlements. Such studies should provide baseline information for antimalarial vaccines now in development.

Head-to-Head Comparison of Three Vaccination Strategies Based on DNA and Raw Insect-Derived Recombinant Proteins against Leishmania

Parasitic diseases plague billions of people among the poorest, killing millions annually, and causing additional millions of disability-adjusted life years lost. Leishmaniases affect more than 12 million people, with over 350 million people at risk. There is an urgent need for efficacious and cheap vaccines and treatments against visceral leishmaniasis (VL), its most severe form. Several vaccination strategies have been proposed but to date no head-to-head comparison was undertaken to assess which is the best in a clinical model of the disease. We simultaneously assayed three vaccination strategies against VL in the hamster model, using KMPII, TRYP, LACK, and PAPLE22 vaccine candidate antigens. Four groups of hamsters were immunized using the following approaches: 1) raw extracts of baculovirus-infected Trichoplusia ni larvae expressing individually one of the four recombinant proteins (PROT); 2) naked pVAX1 plasmids carrying the four genes individually (DNA); 3) a heterologous prime-boost (HPB) strategy involving DNA followed by PROT (DNA-PROT); and 4) a Control including empty pVAX1 plasmid followed by raw extract of wild-type baculovirus-infected T. ni larvae. Hamsters were challenged with L. infantum promastigotes and maintained for 20 weeks. While PROT vaccine was not protective, DNA vaccination achieved protection in spleen. Only DNA-PROT vaccination induced significant NO production by macrophages, accompanied by a significant parasitological protection in spleen and blood. Thus, the DNA-PROT strategy elicits strong immune responses and high parasitological protection in the clinical model of VL, better than its corresponding naked DNA or protein versions. Furthermore, we show that naked DNA coupled with raw recombinant proteins produced in insect larvae biofactories -the cheapest way of producing DNA-PROT vaccines-is a practical and cost-effective way for potential "off the shelf" supplying vaccines at very low prices for the protection against leishmaniases, and possibly against other parasitic diseases affecting the poorest of the poor.

Genetic relatedness of Plasmodium falciparum isolates and the origin of allelic diversity at the merozoite surface protein-1 (MSP-1) locus in Brazil and Vietnam

Abstract Background Despite the extensive polymorphism at the merozoite surface protein-1 (MSP-1) locus of Plasmodium falciparum, that encodes a major repetitive malaria vaccine candidate antigen, identical and nearly identical alleles frequently occur in sympatric parasites. Here we used microsatellite haplotyping to estimate the genetic distance between isolates carrying identical and nearly identical MSP-1 alleles. Methods We analyzed 28 isolates from hypoendemic areas in north-western Brazil, collected between 1985 and 1998, and 23 isolates obtained in mesoendemic southern Vietnam in 1996. MSP-1 alleles were characterized by combining PCR typing with allele-specific primers and partial DNA sequencing. The following single-copy microsatellite markers were typed : Polyα, TA42 (only for Brazilian samples), TA81, TA1, TA87, TA109 (only for Brazilian samples), 2490, ARAII, PfG377, PfPK2, and TA60. Results The low pair-wise average genetic distance between microsatellite haplotypes of isolates sharing identical MSP-1 alleles indicates that epidemic propagation of discrete parasite clones originated most identical MSP-1 alleles in parasite populations from Brazil and Vietnam. At least one epidemic clone propagating in Brazil remained relatively unchanged over more than one decade. Moreover, we found no evidence that rearrangements of MSP-1 repeats, putatively created by mitotic recombination events, generated new alleles within clonal lineages of parasites in either country. Conclusion Identical MSP-1 alleles originated from co-ancestry in both populations, whereas nearly identical MSP-1 alleles have probably appeared independently in unrelated parasite lineages.

Therapeutic DNA Vaccine Encoding Peptide P10 against Experimental Paracoccidioidomycosis

Paracoccidioidomycosis (PCM), caused by Paracoccidioides brasiliensis, is the most prevalent invasive fungal disease in South America. Systemic mycoses are the 10th most common cause of death among infectious diseases in Brazil and PCM is responsible for more than 50% of deaths due to fungal infections. PCM is typically treated with sulfonamides, amphotericin B or azoles, although complete eradication of the fungus may not occur and relapsing disease is frequently reported. A 15-mer peptide from the major diagnostic antigen gp43, named P10, can induce a strong T-CD4+ helper-1 immune response in mice. The TEPITOPE algorithm and experimental data have confirmed that most HLA-DR molecules can present P10, which suggests that P10 is a candidate antigen for a PCM vaccine. In the current work, the therapeutic efficacy of plasmid immunization with P10 and/or IL-12 inserts was tested in murine models of PCM. When given prior to or after infection with P. brasiliensis virulent Pb 18 isolate, plasmid-vaccination with P10 and/or IL-12 inserts successfully reduced the fungal burden in lungs of infected mice. In fact, intramuscular administration of a combination of plasmids expressing P10 and IL-12 given weekly for one month, followed by single injections every month for 3 months restored normal lung architecture and eradicated the fungus in mice that were infected one month prior to treatment. The data indicate that immunization with these plasmids is a powerful procedure for prevention and treatment of experimental PCM, with the perspective of being also effective in human patients.

Exome analysis of HIV patients submitted to dendritic cells therapeutic vaccine reveals an association of CNOT1 gene with response to the treatment

INTRODUCTION: With the aim of searching genetic factors associated with the response to an immune treatment based on autologous monocyte-derived dendritic cells pulsed with autologous inactivated HIV, we performed exome analysis by screening more than 240,000 putative functional exonic variants in 18 HIV-positive Brazilian patients that underwent the immune treatment. METHODS: Exome analysis has been performed using the ILLUMINA Infinium HumanExome BeadChip. zCall algorithm allowed us to recall rare variants. Quality control and SNP-centred analysis were done with GenABEL R package. An in-house implementation of the Wang method permitted gene-centred analysis. RESULTS: CCR4-NOT transcription complex, subunit 1 (CNOT1) gene (16q21), showed the strongest association with the modification of the response to the therapeutic vaccine (p=0.00075). CNOT1 SNP rs7188697 A/G was significantly associated with DC treatment response. The presence of a G allele indicated poor response to the therapeutic vaccine (p=0.0031; OR=33.00; CI=1.74-624.66), and the SNP behaved in a dominant model (A/A vs. A/G+G/G p=0.0009; OR=107.66; 95% CI=3.85-3013.31), being the A/G genotype present only in weak/transient responders, conferring susceptibility to poor response to the immune treatment. DISCUSSION: CNOT1 is known to be involved in the control of mRNA deadenylation and mRNA decay. Moreover, CNOT1 has been recently described as being involved in the regulation of inflammatory processes mediated by tristetraprolin (TTP). The TTP-CCR4-NOT complex (CNOT1 in the CCR4-NOT complex is the binding site for TTP) has been reported as interfering with HIV replication, through post-transcriptional control. Therefore, we can hypothesize that genetic variation occurring in the CNOT1 gene could impair the TTP-CCR4-NOT complex, thus interfering with HIV replication and/or host immune response. CONCLUSIONS: Being aware that our findings are exclusive to the 18 patients studied with a need for replication, and that the genetic variant of CNOT1 gene, localized at intron 3, has no known functional effect, we propose a novel potential candidate locus for the modulation of the response to the immune treatment, and open a discussion on the necessity to consider the host genome as another potential variant to be evaluated when designing an immune therapy study

Pandemic unadjuvanted influenza A (H1N1) vaccine in dermatomyositis and polymyositis: Immunogenicity independent of therapy and no harmful effect in disease

The goal of the present study was to evaluate the influence of the influenza A H1N1/2009 vaccine on dermatomyositis/polymyositis (DM/PM) disease parameters and the potential deleterious effect of therapy on immune response. Thirty-seven DM and 21 PM patients (Bohan and Peter's criteria) were gender- and age-matched to 116 healthy controls. Seroprotection, seroconversion, the geometric mean titers (GMTs) and the factor increase (FI) in the GMTs were calculated. Disease safety was determined from a muscle enzyme analysis and the DM/PM scores [patient's visual analog scale (VAS), physician's VAS, manual muscle strength (MMT-8)] evaluated pre- and post-vaccination. The mean age (43.1 +/- 9.9 vs. 43.8 +/- 8.4 years, p = 0.607) and gender distribution (p = 1.00) were comparable between the patients and controls. After 21 days, seroconversion (p = 0.394), seroprotection (p = 0.08), GMT (p = 0.573) and the FI in the GMT (p = 0.496) were similar in both groups. The disease and muscle parameters remained stable throughout the study, including the creatine kinase (p = 0.20) and aldolase levels (p = 0.98), the physicians' VAS (p = 1.00), the patients' VAS (p = 1.00) and the MMT-8 (p = 1.00). Regarding the influence of treatment, the seroconversion rates were comparable between the controls and patients undergoing treatment with glucocorticoid (GC) (p = 0.969), GC >0.5 mg/kg/day (p = 0.395) and GC + immunosuppressors (p = 0.285). Vaccine-related adverse events were mild and similar in the DM/PM and control groups (p > 0.05). Our data support the administration of the pandemic influenza A H1N1/2009 vaccination in DM/PM, as we found no short-term harmful effects related to the disease itself and adequate immunogenicity in spite of therapy. Further studies are necessary to identify any long-term adverse effects in patients with these diseases.(c) 2012 Elsevier Ltd. All rights reserved.

TNF blockers show distinct patterns of immune response to the pandemic influenza A H1N1 vaccine in inflammatory arthritis patients

Methods. One hundred and twenty patients (RA, n = 41; AS, n = 57; PsA, n = 22) on anti-TNF agents (monoclonal, n = 94; soluble receptor, n = 26) were compared with 116 inflammatory arthritis patients under DMARDs and 117 healthy controls. Seroprotection, seroconversion (SC), geometric mean titre, factor increase in geometric mean titre and adverse events were evaluated 21 days after vaccination. Results. After immunization, SC rates (58.2% vs 74.3%, P = 0.017) were significantly lower in SpA patients receiving anti-TNF therapy, whereas no difference was observed in RA patients receiving this therapy compared with healthy controls (P = 0.067). SpA patients receiving mAbs (infliximab/adalimumab) had a significantly lower SC rate compared with healthy controls (51.6% vs 74.3%, P = 0.002) or those on DMARDs (51.6% vs 74.7%, P = 0.005), whereas no difference was observed for patients on etanercept (86.7% vs 74.3%, P = 0.091). Further analysis of non-seroconverting and seroconverting SpA patients revealed that the former group had a higher mean age (P = 0.003), a higher frequency of anti-TNF (P = 0.031) and mAbs (P = 0.001) and a lower frequency of MTX (P = 0.028). In multivariate logistic regression, only older age (P = 0.015) and mAb treatment (P = 0.023) remained significant factors for non-SC in SpA patients. Conclusion. This study revealed a distinct disease pattern of immune response to the pandemic influenza vaccine in inflammatory arthritis patients receiving anti-TNF agents, illustrated by a reduced immunogenicity solely in SpA patients using mAbs. Trial Registration: ClinicalTrials.gov, ext-link-type="uri" xlink:href="www.clinicaltrials.gov" xmlns:xlink="http://www.w3.org/1999/xlink">www.clinicaltrials.gov, NCT01151644.

Naturally acquired antibodies to Plasmodium vivax blood-stage vaccine candidates (PvMSP-1(19) and PvMSP-3 alpha(359-798)) and their relationship with hematological features in malaria patients from the Brazilian Amazon

An important step when designing a vaccine is identifying the antigens that function as targets of naturally acquired antibodies. We investigated specific antibody responses against two Plasmodium vivax vaccine candidates, PvMSP-1(19) and PvMSP-3 alpha(359-798). Moreover, we assessed the relationship between these antibodies and morbidity parameters. PvMSP-1(19) was the most immunogenic antigen and the frequency of responders to this protein tended to increase in P. vivax patients with higher parasitemia. For both antigens, IgG antibody responses tended to be lower in patients who had experienced their first bout of malaria. Furthermore, anemic patients presented higher IgG antibody responses to PvMSP-3 alpha(359-798). Since the humoral response involves a number of antibodies acting simultaneously on different targets, we performed a Principal Component Analysis (PCA). Anemic patients had, on average, higher first principal component scores (IgG1/IgG2/IgG3/IgG4 anti-MSP3 alpha), which were negatively correlated with hemoglobin levels. Since antibodies against PfMSP-3 have been strongly associated with clinical protection, we cannot exclude the possibility of a dual role of PvMSP-3 specific antibodies in both immunity and pathogenesis of vivax malaria. Our results confirm the high immunogenicity of the conserved C terminus of PvMSP-1 and points to the considerable immunogenicity of polymorphic PvMSP-3 alpha(359-798) during natural infection. (C) 2012 Institut Pasteur. Published by Elsevier Masson SAS. All rights reserved.

Broad and Cross-Clade CD4(+) T-Cell Responses Elicited by a DNA Vaccine Encoding Highly Conserved and Promiscuous HIV-1 M-Group Consensus Peptides

T-cell based vaccine approaches have emerged to counteract HIV-1/AIDS. Broad, polyfunctional and cytotoxic CD4(+) T-cell responses have been associated with control of HIV-1 replication, which supports the inclusion of CD4(+) T-cell epitopes in vaccines. A successful HIV-1 vaccine should also be designed to overcome viral genetic diversity and be able to confer immunity in a high proportion of immunized individuals from a diverse HLA-bearing population. In this study, we rationally designed a multiepitopic DNA vaccine in order to elicit broad and cross-clade CD4(+) T-cell responses against highly conserved and promiscuous peptides from the HIV-1 M-group consensus sequence. We identified 27 conserved, multiple HLA-DR-binding peptides in the HIV-1 M-group consensus sequences of Gag, Pol, Nef, Vif, Vpr, Rev and Vpu using the TEPITOPE algorithm. The peptides bound in vitro to an average of 12 out of the 17 tested HLA-DR molecules and also to several molecules such as HLA-DP, -DQ and murine IA(b) and IA(d). Sixteen out of the 27 peptides were recognized by PBMC from patients infected with different HIV-1 variants and 72% of such patients recognized at least 1 peptide. Immunization with a DNA vaccine (HIVBr27) encoding the identified peptides elicited IFN-gamma secretion against 11 out of the 27 peptides in BALB/c mice; CD4(+) and CD8(+) T-cell proliferation was observed against 8 and 6 peptides, respectively. HIVBr27 immunization elicited cross-clade T-cell responses against several HIV-1 peptide variants. Polyfunctional CD4(+) and CD8(+) T cells, able to simultaneously proliferate and produce IFN-gamma and TNF-alpha, were also observed. This vaccine concept may cope with HIV-1 genetic diversity as well as provide increased population coverage, which are desirable features for an efficacious strategy against HIV-1/AIDS.

Effects of extrusion, lipid concentration and purity on physico-chemical and biological properties of cationic liposomes for gene vaccine applications

We developed cationic liposomes containing DNA through a conventional process involving steps of (i) preformation of liposomes, (ii) extrusion, (iii) drying and rehydration and (iv) DNA complexation. Owing to its high prophylactic potentiality against tuberculosis, which had already been demonstrated in preclinical assays, we introduced modifications into the conventional process towards getting a simpler and more economical process for further scale-up. Elimination of the extrusion step, increasing the lipid concentration (from 16 to 64 mM) of the preformed liposomes and using good manufacturing practice bulk lipids (96-98% purity) instead of analytical grade purity lipids (99.9-100%) were the modifications studied. The differences in the physico-chemical properties, such as average diameter, zeta potential, melting point and morphology of the liposomes prepared through the modified process, were not as significant for the biological properties, such as DNA loading on the cationic liposomes, and effective immune response in mice after immunisation as the control liposomes prepared through the conventional process. Beneficially, the modified process increased productivity by 22% and reduced the cost of raw material by 75%.

Homeobox gene expression profile indicates HOXA5 as a candidate prognostic marker in oral squamous cell carcinoma

The search for molecular markers to improve diagnosis, individualize treatment and predict behavior of tumors has been the focus of several studies. This study aimed to analyze homeobox gene expression profile in oral squamous cell carcinoma (OSCC) as well as to investigate whether some of these genes are relevant molecular markers of prognosis and/or tumor aggressiveness. Homeobox gene expression levels were assessed by microarrays and qRT-PCR in OSCC tissues and adjacent non-cancerous matched tissues (margin), as well as in OSCC cell lines. Analysis of microarray data revealed the expression of 147 homeobox genes, including one set of six at least 2-fold up-regulated, and another set of 34 at least 2-fold down-regulated homeobox genes in OSCC. After qRT-PCR assays, the three most up-regulated homeobox genes (HOXA5, HOXD10 and HOXD11) revealed higher and statistically significant expression levels in OSCC samples when compared to margins. Patients presenting lower expression of HOXA5 had poorer prognosis compared to those with higher expression (P=0.03). Additionally, the status of HOXA5, HOXD10 and HOXD11 expression levels in OSCC cell lines also showed a significant up-regulation when compared to normal oral keratinocytes. Results confirm the presence of three significantly upregulated (>4-fold) homeobox genes (HOXA5, HOXD10 and HOXD11) in OSCC that may play a significant role in the pathogenesis of these tumors. Moreover, since lower levels of HOXA5 predict poor prognosis, this gene may be a novel candidate for development of therapeutic strategies in OSCC.

Gene expression profile during coffee fruit development and identification of candidate markers for phenological stages

The objective of this work was to identify genes that could be used as suitable markers for molecular recognition of phenological stages during coffee (Coffea arabica) fruit development. Four cultivars were evaluated as to their differential expression of genes associated to fruit development and maturation processes. Gene expression was characterized by both semi-quantitative and quantitative RT-PCR, in fruit harvested at seven different developmental stages, during three different seasons. No size polymorphisms or differential expression were observed among the cultivars for the evaluated genes; however, distinct expression profiles along fruit development were determined for each gene. Four out of the 28 evaluated genes exhibited a regular expression profile in all cultivars and harvest seasons, and, therefore, they were validated as candidate phenological markers of coffee fruit. The gene a-galactosidase can be used as a marker of green stage, caffeine synthase as a marker of transition to green and yellowish-green stages, and isocitrate lyase and ethylene receptor 3 as markers of late maturation.