Cloning and expression analysis of S-RNase homologous gene in Citrus reticulata Blanco cv. Wuzishatangju

Downregulated expression of S2-RNase attenuates self-incompatibility in "Guiyou No. 1" pummelo

Self-incompatibility (SI) substantially restricts the yield and quality of citrus. Therefore, breeding and analyzing self-compatible germplasm is of great theoretical and practical significance for citrus. Here, we focus on the mechanism of a self-compatibility mutation in 'Guiyou No. 1' pummelo (Citrus maxima), which is a spontaneous mutant of 'Shatian' pummelo (Citrus maxima, self-incompatibility). The rate of fruit set and the growth of pollen tubes in the pistil confirmed that a spontaneous mutation in the pistil is responsible for the self-compatibility of 'Guiyou No. 1'. Segregation ratios of the S genotype in F₁ progeny, expression analysis, and western blotting validated that the reduced levels of S₂-RNase mRNA contribute to the loss of SI in 'Guiyou No. 1'. Furthermore, we report a phased assembly of the 'Guiyou No. 1' pummelo genome and obtained two complete and well-annotated S haplotypes. Coupled with an analysis of SV variations, methylation levels, and gene expression, we identified a candidate gene (CgHB40), that may influence the regulation of the S₂-RNase promoter. Our data provide evidence that a mutation that affects the pistil led to the loss of SI in 'Guiyou No. 1' by influencing a poorly understood mechanism that affects transcriptional regulation. This work significantly advances our understanding of the genetic basis of the SI system in citrus and provides information on the regulation of S-RNase genes.

SKP1-like protein, CrSKP1-e, interacts with pollen-specific F-box proteins and assembles into SCF-type E3 complex in 'Wuzishatangju' (Citrus reticulata Blanco) pollen

S-ribonuclease (S-RNase)-based self-incompatibility (SI) mechanisms have been extensively studied in Solanaceae, Rosaceae and Plantaginaceae. S-RNase-based SI is controlled by two closely related genes, S-RNase and S-locus F-box (SLF), located at a polymorphic S-locus. In the SI system, the SCF-type (SKP1-CUL1-F-box-RBX1) complex functions as an E3 ubiquitin ligase complex for ubiquitination of non-self S-RNase. Pummelo (Citrus grandis) and several mandarin cultivars are suggested to utilize an S-RNase-based SI system. However, the molecular mechanism of the non-S-factors involved in the SI reaction is not straightforward in Citrus. To investigate the SCF-type E3 complex responsible for the SI reaction in mandarin, SLF, SKP1-like and CUL1 candidates potentially involved in the SI reaction of ‘Wuzishatangju’ (Citrus reticulata Blanco) were identified based on the genome-wide identification and expression analyses. Sixteen pollen-specific F-box genes (CrFBX1-CrFBX16), one pollen-specific SKP1-like gene (CrSKP1-e) and two CUL1 genes (CrCUL1A and CrCUL1B) were identified and cloned from ‘Wuzishatangju’. Yeast two-hybrid (Y2H) and in vitro binding assays showed that five CrFBX proteins could bind to CrSKP1-e, which is an ortholog of SSK1 (SLF-interacting-SKP1-like), a non-S-factor responsible for the SI reaction. Luciferase complementation imaging (LCI) and in vitro binding assays also showed that CrSKP1-e interacts with the N-terminal region of both CrCUL1A and CrCUL1B. These results indicate that CrSKP1-e may serve as a functional member of the SCF-type E3 ubiquitin ligase complex in ‘Wuzishatangju’.

Self-(In)compatibility Systems: Target Traits for Crop-Production, Plant Breeding, and Biotechnology

Self-incompatibility (SI) mechanisms prevent self-fertilization in flowering plants based on specific discrimination between self- and non-self pollen. Since this trait promotes outcrossing and avoids inbreeding it is a widespread mechanism of controlling sexual plant reproduction. Growers and breeders have effectively exploited SI as a tool for manipulating domesticated crops for thousands of years. However, only within the past thirty years have studies begun to elucidate the underlying molecular features of SI. The specific S-determinants and some modifier factors controlling SI have been identified in the sporophytic system exhibited by Brassica species and in the two very distinct gametophytic systems present in Papaveraceae on one side and in Solanaceae, Rosaceae, and Plantaginaceae on the other. Molecular level studies have enabled SI to SC transitions (and vice versa) to be intentionally manipulated using marker assisted breeding and targeted approaches based on transgene integration, silencing, and more recently CRISPR knock-out of SI-related factors. These scientific advances have, in turn, provided a solid basis to implement new crop production and plant breeding practices. Applications of self-(in)compatibility include widely differing objectives such as crop yield and quality improvement, marker-assisted breeding through SI genotyping, and development of hybrids for overcoming intra- and interspecific reproductive barriers. Here, we review scientific progress as well as patented applications of SI, and also highlight future prospects including further elucidation of SI systems, deepening our understanding of SI-environment relationships, and new perspectives on plant self/non-self recognition.

Overexpression of a Novel ROP Gene from the Banana (MaROP5g) Confers Increased Salt Stress Tolerance

Rho-like GTPases from plants (ROPs) are plant-specific molecular switches that are crucial for plant survival when subjected to abiotic stress. We identified and characterized 17 novel ROP proteins from Musa acuminata (MaROPs) using genomic techniques. The identified MaROPs fell into three of the four previously described ROP groups (Groups II⁻IV), with MaROPs in each group having similar genetic structures and conserved motifs. Our transcriptomic analysis showed that the two banana genotypes tested, Fen Jiao and BaXi Jiao, had similar responses to abiotic stress: Six genes (MaROP-3b, -5a, -5c, -5f, -5g, and -6) were highly expressed in response to cold, salt, and drought stress conditions in both genotypes. Of these, MaROP5g was most highly expressed in response to salt stress. Co-localization experiments showed that the MaROP5g protein was localized at the plasma membrane. When subjected to salt stress, transgenic Arabidopsis thaliana overexpressing MaROP5g had longer primary roots and increased survival rates compared to wild-type A. thaliana. The increased salt tolerance conferred by MaROP5g might be related to reduced membrane injury and the increased cytosolic K⁺/Na⁺ ratio and Ca²⁺ concentration in the transgenic plants as compared to wild-type. The increased expression of salt overly sensitive (SOS)-pathway genes and calcium-signaling pathway genes in MaROP5g-overexpressing A. thaliana reflected the enhanced tolerance to salt stress by the transgenic lines in comparison to wild-type. Collectively, our results suggested that abiotic stress tolerance in banana plants might be regulated by multiple MaROPs, and that MaROP5g might enhance salt tolerance by increasing root length, improving membrane injury and ion distribution.

Reproduction in woody perennial Citrus: an update on nucellar embryony and self-incompatibility

Review on citrus reproduction. Citrus is one of the most important and widely grown fruit crops. It possesses several special reproductive characteristics, such as nucellar embryony and self-incompatibility. The special phenomenon of nucellar embryony in citrus, also known as the polyembryony, is a kind of sporophytic apomixis. During the past decade, the emergence of novel technologies and the construction of multiple citrus reference genomes have facilitated rapid advances to our understanding of nucellar embryony. Indeed, several research teams have preliminarily determined the genetic basis of citrus apomixis. On the other hand, the phenomenon of self-incompatibility that promotes genetic diversity by rejecting self-pollen and accepting non-self-pollen is difficult to study in citrus because the long juvenile period of citrus presents challenges to identifying candidate genes that control this phenomenon. In this review, we focus on advances to our understanding of reproduction in citrus from the last decade and discuss priorities for the coming decade.

Temperatures during flower bud development affect pollen germination, self-incompatibility reaction and early fruit development of clementine (Citrus clementina Hort. ex Tan.)

One of the key environmental factors affecting plant reproductive systems is temperature. Characterising such effects is especially relevant for some commercially important genera such as Citrus. In this genus, failure of fertilisation results in parthenocarpic fruit development and seedlessness, which is a much-prized character. Here, we characterise the effects of temperature on flower and ovary development, and on pollen-pistil interactions in 'Comune' clementine (Citrus clementina Hort. ex Tan.). We examine flower bud development, in vitro pollen germination and pollen-pistil interaction at different temperatures (15, 20, 25 or 30 °C). These temperatures span the range from 'cold' to 'hot' weather during the flowering season in many citrus-growing regions. Temperature had a strong effect on flower and ovary development, pollen germination, and pollen tube growth kinetics. In particular, parthenocarpic fruit development (indicated by juice vesicle growth) was initiated early if flowers were exposed to warmer temperatures during anthesis. Exposure to different temperatures during flower bud development also alters expression of the self-incompatibility reaction. This affects the point in the pistil at which pollen tube growth is arrested and confirms the role of sub- and supra-optimal temperatures in determining the numbers of pollen tubes reaching the ovary.

Comparative transcriptional survey between self-incompatibility and self-compatibility in Citrus reticulata Blanco

Seedlessness is an excellent economical trait, and self-incompatibility (SI) is one of important factors resulting in seedless fruit in Citrus. However, SI molecular mechanism in Citrus is still unclear. In this study, RNA-Seq technology was used to identify differentially expressed genes related to SI reaction of 'Wuzishatangju' (Citrus reticulata Blanco). A total of 35.67GB raw RNA-Seq data was generated and was de novo assembled into 50,364 unigenes with an average length of 897bp and N50 value of 1549. Twenty-three candidate unigenes related to SI were analyzed using qPCR at different tissues and stages after self- and cross-pollination. Seven pollen S genes (Unigene0050323, Unigene0001060, Unigene0004230, Unigene0004222, Unigene0012037, Unigene0048889 and Unigene0004272), three pistil S genes (Unigene0019191, Unigene0040115, Unigene0036542) and three genes (Unigene0038751, Unigene0031435 and Unigene0029897) associated with the pathway of ubiquitin-mediated proteolysis were identified. Unigene0031435, Unigene0038751 and Unigene0029897 are probably involved in SI reaction of 'Wuzishatangju' based on expression analyses. The present study provides a new insight into the molecular mechanism of SI in Citrus at the transcriptional level.

Genome-wide identification and functional analysis of S-RNase involved in the self-incompatibility of citrus

S-RNase-based self-incompatibility is found in Solanaceae, Rosaceae, and Scrophulariaceae, and is the most widespread mechanism that prevents self-fertilization in plants. Although 'Shatian' pummelo (Citrus grandis), a traditional cultivated variety, possesses the self-incompatible trait, the role of S-RNases in the self-incompatibility of 'Shatian' pummelo is poorly understood. To identify genes associated with self-incompatibility in citrus, we identified 16 genes encoding homologs of ribonucleases in the genomes of sweet orange (Citrus sinensis) and clementine mandarin (Citrus clementine). We preliminarily distinguished S-RNases from S-like RNases with a phylogenetic analysis that classified these homologs into three groups, which is consistent with the previous reports. Expression analysis provided evidence that CsRNS1 and CsRNS6 are S-like RNase genes. The expression level of CsRNS1 was increased during fruit development. The expression of CsRNS6 was increased during the formation of embryogenic callus. In contrast, we found that CsRNS3 possessed several common characteristics of the pistil determinant of self-incompatibility: it has an alkaline isoelectric point (pI), harbors only one intron, and is specifically expressed in style. We obtained a cDNA encoding CgRNS3 from 'Shatian' pummelo and found that it is homolog to CsRNS3 and that CgRNS3 exhibited the same expression pattern as CsRNS3. In an in vitro culture system, the CgRNS3 protein significantly inhibited the growth of self-pollen tubes from 'Shatian' pummelo, but after a heat treatment, this protein did not significantly inhibit the elongation of self- or non-self-pollen tubes. In conclusion, an S-RNase gene, CgRNS3, was obtained by searching the genomes of sweet orange and clementine for genes exhibiting sequence similarity to ribonucleases followed by expression analyses. Using this approach, we identified a protein that significantly inhibited the growth of self-pollen tubes, which is the defining property of an S-RNase.

CrWSKP1, an SKP1-like Gene, Is Involved in the Self-Incompatibility Reaction of "Wuzishatangju" (Citrus reticulata Blanco)

Plant S-phase kinase-associated protein 1 (SKP1) genes play crucial roles in plant development and differentiation. However, the role of SKP1 in citrus is unclear. Herein, we described a novel SKP1-like gene, designated as CrWSKP1, from "Wuzishatangju" (Citrus reticulata Blanco). The cDNA sequence of CrWSKP1 is 779 base pairs (bp) and contains an open reading frame (ORF) of 477 bp. The genomic sequence of the CrWSKP1 gene is 1296 bp with two exons and one intron. CrWSKP1 has high identity with SKP1-like genes from other plant species within two conserved regions. Approximately 85% of pollen tubes of self-pollinated CrWSKP1 transgenic tobaccos became twisted at four days after self-pollination. Pollen tube numbers of self-pollinated CrWSKP1 transformants entering into ovules were significantly fewer than that of the control. Seed number of self-pollinated CrWSKP1 transformants was significantly reduced. These results suggested that the CrWSKP1 is involved in the self-incompatibility (SI) reaction of "Wuzishatangju".

No evidence for Fabaceae Gametophytic self-incompatibility being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes

BACKGROUND

Fabaceae species are important in agronomy and livestock nourishment. They have a long breeding history, and most cultivars have lost self-incompatibility (SI), a genetic barrier to self-fertilization. Nevertheless, to improve legume crop breeding, crosses with wild SI relatives of the cultivated varieties are often performed. Therefore, it is fundamental to characterize Fabaceae SI system(s). We address the hypothesis of Fabaceae gametophytic (G)SI being RNase based, by recruiting the same S-RNase lineage gene of Rosaceae, Solanaceae or Plantaginaceae SI species.

RESULTS

We first identify SSK1 like genes (described only in species having RNase based GSI), in the Trifolium pratense, Medicago truncatula, Cicer arietinum, Glycine max, and Lupinus angustifolius genomes. Then, we characterize the S-lineage T2-RNase genes in these genomes. In T. pratense, M. truncatula, and C. arietinum we identify S-RNase lineage genes that in phylogenetic analyses cluster with Pyrinae S-RNases. In M. truncatula and C. arietinum genomes, where large scaffolds are available, these sequences are surrounded by F-box genes that in phylogenetic analyses also cluster with S-pollen genes. In T. pratense the S-RNase lineage genes show, however, expression in tissues not involved in GSI. Moreover, levels of diversity are lower than those observed for other S-RNase genes. The M. truncatula and C. arietinum S-RNase and S-pollen like genes phylogenetically related to Pyrinae S-genes, are also expressed in tissues other than those involved in GSI. To address if other T2-RNases could be determining Fabaceae GSI, here we obtained a style with stigma transcriptome of Cytisus striatus, a species that shows significant difference on the percentage of pollen growth in self and cross-pollinations. Expression and polymorphism analyses of the C. striatus S-RNase like genes revealed that none of these genes, is the S-pistil gene.

CONCLUSION

We find no evidence for Fabaceae GSI being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes. There is no evidence that T2-RNase lineage genes could be determining GSI in C. striatus. Therefore, to characterize the Fabaceae S-pistil gene(s), expression analyses, levels of diversity, and segregation analyses in controlled crosses are needed for those genes showing high expression levels in the tissues where GSI occurs.

Molecular characterization and expression analysis of S1 self-incompatibility locus-linked pollen 3.15 gene in Citrus reticulata

Gametophytic self-incompatibility (GSI) is controlled by a highly polymorphic locus called the S-locus, which is an important factor that can result in seedless fruit in Citrus. The S1 self-incompatibility locus-linked pollen 3.15 gene (S1-3.15 ) belongs to a type of S locus gene. The role of S1-3.15 in the SI reaction of Citrus has not yet been reported. In this study, full-length sequences of cDNA and DNA encoding the S1-3.15 gene, referred to as CrS1-3.15 , were isolated from 'Wuzishatangju' (Self-incompatibility, SI) and 'Shatangju' (Self-compatibility, SC). The predicted amino acid sequences of CrS1-3.15 between 'Wuzishatangju' and 'Shatangju' differ by only three amino acids. Compared to 'Wuzishatangju', three bases were substituted in the genomic DNA of CrS1-3.15 from 'Shatangju'. Southern blot results showed that one copy of CrS1-3.15 existed in the genomic DNA of both 'Wuzishatangju' and 'Shatangju'. The expression level of the CrS1-3.15 gene in the ovaries of 'Shatangju' was approximately 60-fold higher than that in the ovaries of 'Wuzishatangju'. When 'Wuzishatangju' was cross-pollinated, the expression of CrS1-3.15 was upregulated in the ovaries at 3 d, and the highest expression levels were detected in the ovaries at 6 d after cross-pollination of 'Wuzishatangju' × 'Shatangju'. To obtain the CrS1-3.15 protein, the full-length cDNA of CrS1-3.15 genes from 'Wuzishatangju' and 'Shatangju' was successfully expressed in Pichia pastoris. Pollen germination frequency of 'Wuzishatangju' was inhibited significantly with increasing CrS1-3.15 protein concentrations from SI 'Wuzishatangju'.

Identifying differentially expressed genes in pollen from self-incompatible "Wuzishatangju" and self-compatible "Shatangju" mandarins

Self-incompatibility (SI) is one of the important factors that can result in seedless fruit in Citrus. However, the molecular mechanism of SI in Citrus is not yet clear. In this study, two suppression subtractive hybridization (SSH) libraries (forward, F and reverse, R) were constructed to isolate differentially expressed genes in pollen from "Wuzishatangju" (SI) and "Shatangju" (self-compatibility, SC) mandarins. Four hundred and sixty-eight differentially expressed cDNA clones from 2077 positive clones were sequenced and identified. Differentially expressed ESTs are possibly involved in the SI reaction of "Wuzishatangju" by regulating pollen development, kinase activity, ubiquitin pathway, pollen-pistil interaction, and calcium ion binding. Twenty five SI candidate genes were obtained, six of which displayed specific expression patterns in various organs and stages after self- and cross-pollination. The expression level of the F-box gene (H304) and S1 (F78) in the pollen of "Wuzishatangju" was 5-fold higher than that in "Shatangju" pollen. The F-box gene, S1, UBE2, UBE3, RNaseHII, and PCP were obviously up-regulated in pistils at 3 d after self-pollination of "Wuzishatangju", approximately 3-, 2-, 10-, 5-, 5-, and 2-fold higher, respectively than that at the same stage after cross-pollination of "Wuzishatangju" × "Shatangju" pistils. The potential involvement of these genes in the pollen SI reaction of "Wuzishatangju" is discussed.

Molecular characterization and expression analysis of ubiquitin-activating enzyme E1 gene in Citrus reticulata

Ubiquitin-activating enzyme E1 (UBE1) catalyzes the first step in the ubiquitination reaction, which targets a protein for degradation via a proteasome pathway. UBE1 plays an important role in metabolic processes. In this study, full-length cDNA and DNA sequences of UBE1 gene, designated CrUBE1, were obtained from 'Wuzishatangju' (self-incompatible, SI) and 'Shatangju' (self-compatible, SC) mandarins. 5 amino acids and 8 bases were different in cDNA and DNA sequences of CrUBE1 between 'Wuzishatangju' and 'Shatangju', respectively. Southern blot analysis showed that there existed only one copy of the CrUBE1 gene in genome of 'Wuzishatangju' and 'Shatangju'. The temporal and spatial expression characteristics of the CrUBE1 gene were investigated using semi-quantitative RT-PCR (SqPCR) and quantitative real-time PCR (qPCR). The expression level of the CrUBE1 gene in anthers of 'Shatangju' was approximately 10-fold higher than in anthers of 'Wuzishatangju'. The highest expression level of CrUBE1 was detected in pistils at 7days after self-pollination of 'Wuzishatangju', which was approximately 5-fold higher than at 0 h. To obtain CrUBE1 protein, the full-length cDNA of CrUBE1 genes from 'Wuzishatangju' and 'Shatangju' were successfully expressed in Pichia pastoris. Pollen germination frequency of 'Wuzishatangju' was significantly inhibited with increasing of CrUBE1 protein concentrations from 'Wuzishatangju'.

Identification of differentially expressed genes in pistils from self-incompatible Citrus reticulata by suppression subtractive hybridization

Self-incompatibility (SI) is one important factor that can result in Citrus seedlessness. However, the molecular mechanism of SI in Citrus is not clear yet. To isolate the pistil's SI-related genes, a suppression subtractive hybridization library was constructed using mature pistils of 'Wuzishatangju' mandarin (SI) as the tester and mature pistils of 'Shatangju' mandarin (self-compatibility, SC) as the driver. 229 differentially expressed cDNA clones from 967 positive clones were sequenced and identified. Differentially expressed ESTs are possibly involved in the SI reaction of 'Wuzishatangju' through a regulating signaling pathway, serine/threonine phosphatase activity, receptor kinase, embryonic development, gibberellin stimulus, or transcription. 11 out of 36 SI candidate genes displayed different expression patterns in various tissues and stages after self- and cross-pollination of 'Wuzishatangju'. The expression of CaBP (WY65), a senescence-protease (WY372), an unknown gene (WY283), and a WRKY (WY17) were up-regulated in the styles of 'Wuzishatangju' while higher expression of WY190 was observed in styles of 'Shatangju'. Highest expression levels of WY65, WY372, an annexin (WY598), the zinc-finger protein (WY376), a C2-protein (WY291), and an unknown gene (WY318) were detected in styles at 3 days after self-pollination of 'Wuzishatangju' while lowest levels were observed in styles at 3 days after cross-pollination of 'Wuzishatangju' × 'Shatangju'. The potential involvement of these genes in the SI reaction is discussed.

Comparative transcriptome analysis of stylar canal cells identifies novel candidate genes implicated in the self-incompatibility response of Citrus clementina

BACKGROUND

Reproductive biology in citrus is still poorly understood. Although in recent years several efforts have been made to study pollen-pistil interaction and self-incompatibility, little information is available about the molecular mechanisms regulating these processes. Here we report the identification of candidate genes involved in pollen-pistil interaction and self-incompatibility in clementine (Citrus clementina Hort. ex Tan.). These genes have been identified comparing the transcriptomes of laser-microdissected stylar canal cells (SCC) isolated from two genotypes differing for self-incompatibility response ('Comune', a self-incompatible cultivar and 'Monreal', a self- compatible mutation of 'Comune').

RESULTS

The transcriptome profiling of SCC indicated that the differential regulation of few specific, mostly uncharacterized transcripts is associated with the breakdown of self-incompatibility in 'Monreal'. Among them, a novel F-box gene showed a drastic up-regulation both in laser microdissected stylar canal cells and in self-pollinated whole styles with stigmas of 'Comune' in concomitance with the arrest of pollen tube growth. Moreover, we identify a non-characterized gene family as closely associated to the self-incompatibility genetic program activated in 'Comune'. Three different aspartic-acid rich (Asp-rich) protein genes, located in tandem in the clementine genome, were over-represented in the transcriptome of 'Comune'. These genes are tightly linked to a DELLA gene, previously found to be up-regulated in the self-incompatible genotype during pollen-pistil interaction.

CONCLUSION

The highly specific transcriptome survey of the stylar canal cells identified novel genes which have not been previously associated with self-pollen rejection in citrus and in other plant species. Bioinformatic and transcriptional analyses suggested that the mutation leading to self-compatibility in 'Monreal' affected the expression of non-homologous genes located in a restricted genome region. Also, we hypothesize that the Asp-rich protein genes may act as Ca2+ "entrapping" proteins, potentially regulating Ca2+ homeostasis during self-pollen recognition.

Reference genes for accurate transcript normalization in citrus genotypes under different experimental conditions

Real-time reverse transcription PCR (RT-qPCR) has emerged as an accurate and widely used technique for expression profiling of selected genes. However, obtaining reliable measurements depends on the selection of appropriate reference genes for gene expression normalization. The aim of this work was to assess the expression stability of 15 candidate genes to determine which set of reference genes is best suited for transcript normalization in citrus in different tissues and organs and leaves challenged with five pathogens (Alternaria alternata, Phytophthora parasitica, Xylella fastidiosa and Candidatus Liberibacter asiaticus). We tested traditional genes used for transcript normalization in citrus and orthologs of Arabidopsis thaliana genes described as superior reference genes based on transcriptome data. geNorm and NormFinder algorithms were used to find the best reference genes to normalize all samples and conditions tested. Additionally, each biotic stress was individually analyzed by geNorm. In general, FBOX (encoding a member of the F-box family) and GAPC2 (GAPDH) was the most stable candidate gene set assessed under the different conditions and subsets tested, while CYP (cyclophilin), TUB (tubulin) and CtP (cathepsin) were the least stably expressed genes found. Validation of the best suitable reference genes for normalizing the expression level of the WRKY70 transcription factor in leaves infected with Candidatus Liberibacter asiaticus showed that arbitrary use of reference genes without previous testing could lead to misinterpretation of data. Our results revealed FBOX, SAND (a SAND family protein), GAPC2 and UPL7 (ubiquitin protein ligase 7) to be superior reference genes, and we recommend their use in studies of gene expression in citrus species and relatives. This work constitutes the first systematic analysis for the selection of superior reference genes for transcript normalization in different citrus organs and under biotic stress.