S-locus F-Box Protein as Pollen S determinant Targets Non-self S-RNase underlying Self-Incompatibility in Citrus.

Self-incompatibility (SI) is a crucial mechanism that prevents self-fertilization and inbreeding in flowering plants. Citrus exhibits SI regulated by a polymorphic S-locus containing an S-RNase gene and multiple S-locus F-box (SLF) genes. It has been documented that S-RNase functions as the pistil S determinant, but there is no direct evidence that the SLFs closely linked with S-RNase function as pollen S determinants in Citrus. This study assembled the genomes of two pummelo (Citrus grandis) plants and obtained three novel complete and well-annotated S-haplotypes and isolated 36 SLF or SLF-like alleles on the S-loci. Phylogenetic analysis of 138 SLFs revealed that the SLFs were classified into 12 types, including six types with divergent or missing alleles. Furthermore, transformation experiments verified that the conserved S6-SLF7a protein can lead the transition of SI to self-compatibility (SC) by recognizing non-self S8-RNase in 'Mini-Citrus' plants (S7S8 and S8S29, Fortunella hindsii), a model plant for citrus gene function studies. In vitro assays demonstrated interactions between SLFs of different S haplotypes and the Skp1-Cullin1-F-box subunit CgSSK1 protein. This study provides direct evidence that SLF controls the pollen function in Citrus, demonstrating its role in the 'non-self-recognition' SI system.

Acetylome reprograming participates in the establishment of fruit metabolism during polyploidization in citrus.

Polyploidization leads to novel phenotypes and is a major force in evolution. However, the relationship between the evolution of new traits and variations in the post-translational modifications (PTM) of proteins during polyploidization has not been studied. Acetylation of lysine residues is a common protein PTM that plays a critical regulatory role in central metabolism. To test whether changes in metabolism in citrus fruit is associated with the reprogramming of lysine acetylation (Kac) in non-histone proteins during allotetraploidization, we performed a global acetylome analysis of fruits from a synthetic allotetraploid citrus and its diploid parents. A total of 4,175 Kac sites were identified on 1,640 proteins involved in a wide range of fruit traits. In the allotetraploid, parental dominance (i.e. resemblance to one of the two parents) in specific fruit traits, such as fruit acidity and flavonol metabolism, was highly associated with parental Kac level dominance in pertinent enzymes. This association is due to Kac-mediated regulation of enzyme activity. Moreover, protein Kac probably contributes to the discordance between the transcriptomic and proteomic variations during allotetraploidization. The acetylome reprogramming can be partially explained by the expression pattern of several lysine deacetylases (KDACs). Overexpression of silent information regulator 2 (CgSRT2) and histone deacetylase 8 (CgHDA8) diverted metabolic flux from primary metabolism to secondary metabolism and partially restored a metabolic status to the allotetraploid, which expressed attenuated levels of CgSRT2 and CgHDA8. Additionally, KDAC inhibitor treatment greatly altered metabolism in citrus fruit. Collectively, these findings reveal the important role of acetylome reprogramming in trait evolution during polyploidization.

Transcriptomes and DNA methylomes in apomictic cells delineate nucellar embryogenesis initiation in citrus.

Citrus nucellar poly-embryony (NPE) is a mode of sporophytic apomixis that asexual embryos formed in the seed through adventitious embryogenesis from the somatic nucellar cells. NPE allows clonal propagation of rootstocks, but it impedes citrus cross breeding. To understand the cellular processes involved in NPE initiation, we profiled the transcriptomes and DNA methylomes in laser microdissection captured citrus apomictic cells. In apomictic cells, ribosome biogenesis and protein degradation were activated, whereas auxin polar transport was repressed. Reactive oxygen species (ROS) accumulated in the poly-embryonic ovules, and response to oxidative stress was provoked. The global DNA methylation level, especially that of CHH context, was decreased, whereas the methylation level of the NPE-controlling key gene CitRWP was increased. A C2H2 domain-containing transcription factor gene and CitRWP co-expressed specifically in apomictic cells may coordinate to initiate NPE. The activated embryogenic development and callose deposition processes indicated embryogenic fate of nucellar embryo initial (NEI) cells. In our working model for citrus NPE initiation, DNA hyper-methylation may activate transcription of CitRWP, which increases C2H2 expression and ROS accumulation, triggers epigenetic regulation and regulates cell fate transition and NEI cell identity in the apomictic cells.

Identification of a delayed leaf greening gene from a mutation of pummelo.

Delayed greening of young leaves is an unusual phenomenon of plants in nature. Citrus are mostly evergreen tree species. Here, a natural mutant of "Guanxi" pummelo (Citrus maxima), which shows yellow leaves at the young stage, was characterized to identify the genes underlying the trait of delayed leaf greening in plants. A segregating population with this mutant as the seed parent and a normal genotype as the pollen parent was generated. Two DNA pools respectively from the leaves of segregating seedlings with extreme phenotypes of normal leaf greening and delayed leaf greening were collected for sequencing. Bulked segregant analysis (BSA) and InDel marker analysis demonstrated that the delayed leaf greening trait is governed by a 0.3 Mb candidate region on chromosome 6. Gene expression analysis further identified a key candidate gene (Citrus Delayed Greening gene 1, CDG1) in the 0.3 Mb region, which showed significantly differential expression between the genotypes with delayed and normal leaf greening phenotypes. There was a 67 bp InDel region difference in the CDG1 promoter and the InDel region contains a TATA-box element. Confocal laser-scanning microscopy revealed that the CDG1-GFP fusion protein signals were co-localized with the chloroplast signals in the protoplasts. Overexpression of CDG1 in tobacco and Arabidopsis led to the phenotype of delayed leaf greening. These results suggest that the CDG1 gene is involved in controlling the delayed leaf greening phenotype with important functions in chloroplast development.

The miR399-CsUBC24 Module Regulates Reproductive Development and Male Fertility in Citrus.

MicroRNA399 (miR399) regulates phosphate homeostasis in plants by down-regulating the expression of PHOSPHATE2 (PHO2, or UBC24 encoding the ubiquitin-conjugating E2 enzyme). We previously identified CsmiR399a.1 in a small RNA sequencing screen of a male-sterile somatic cytoplasmic hybrid (or cybrid) of pummelo (Citrus grandis). Here, we report that miR399 affects reproductive development and male fertility in citrus. Down-regulation of CsmiR399a.1 using a short tandem target mimic (STTM) led to abnormal floral development, inhibition of anther dehiscence, and decreased pollen fertility. When grown in inorganic phosphate (Pi)-sufficient conditions, CsmiR399a.1-STTM plants had lower total phosphorus content in their leaves than the wild type and showed typical symptoms of Pi deficiency. In CsmiR399a.1-STTM plants, the expression of genes involved in starch metabolism and Pi homeostasis was significantly different than in the wild type. Thus, we conclude that miR399-STTM mimicked Pi deficiency, disturbed starch metabolism, and was responsible for pollen grain collapse in the transgenic lines. We identified CsUBC24, a citrus homolog of Arabidopsis (Arabidopsis thaliana) AtUBC24 (PHO2), as a target of CsmiR399a.1 that physically interacts with the floral development regulators SEPALLATA family (CsSEP1.1, CsSEP1.2, and CsSEP3) and the anther dehiscence regulator INDUCER OF CBF EXPRESSION1 (CsICE1). We hypothesize that CsUBC24 downregulates the CsSEPs, which disrupts the floral meristem identity regulatory network and leads to developmental abnormalities in flowers. By interacting with CsICE1, CsUBC24 disturbs stomate function on the anther surface, which inhibits anther dehiscence. These findings indicate that a miR399-based mechanism influences both reproductive development and male fertility in citrus.

Glycemic index, glycemic load, and glycemic response to pomelo in patients with type 2 diabetes.

Food intake has a great influence on blood glucose in patients with diabetes. This study was to determine the glycemic index (GI) and glycemic load (GL) of a particular pomelo named Majia pomelo and its effects on postprandial glucose (PPG) in patients with type 2 diabetes (T2D). Twenty healthy subjects and 20 T2D patients (controlled on lifestyle measures and/or metformin) were tested on 2 separate days with 50 g of glucose and 50 g equivalent of carbohydrates from Majia pomelo for GI measurement. To test effects of Majia pomelo on PPG, 19 hospitalized T2D patients (controlled on insulin therapy) were selected for a 9-day study. The dose of insulin for each patient was adjusted on the first 3 days. A total of 100 mg Majia pomelo was consumed per meal in the last 3 tested days. Blood glucose was measured to evaluate the glycemic excursions. The GIs for Majia pomelo in healthy individuals and T2D patients were 78.34±1.88 and 72.15±1.95 respectively. The value of GL was as low as 4.23 in diabetic patients with serving size of 100 g pomelo, indicting Majia pomelo as a high GI but low GL fruit. Consumption of Majia pomelo in hospitalized T2D patients did not cause significant glucose fluctuation. It was concluded that high GI pomelo can serve as a low GL fruit if it is consumed with a limited daily amount and thus can be supplied to diabetic patients. These results may mean more varieties of food choices for T2D patients.

Genome-scale mRNA and small RNA transcriptomic insights into initiation of citrus apomixis.

Nucellar embryony (NE) is an adventitious form of apomixis common in citrus, wherein asexual embryos initiate directly from nucellar cells surrounding the embryo sac. NE enables the fixation of desirable agronomic traits and the production of clonal offspring of virus-free rootstock, but impedes progress in hybrid breeding. In spite of the great importance of NE in citrus breeding and commercial production, little is understood about the underlying molecular mechanisms. In this study, the stages of nucellar embryo initiation (NEI) were determined for two polyembryonic citrus cultivars via histological observation. To explore the genes and regulatory pathways involved in NEI, we performed mRNA-seq and sRNA-seq analyses of ovules immediately prior to and at stages during NEI in the two pairs of cultivars. A total of 305 differentially expressed genes (DEGs) were identified between the poly- and monoembryonic ovules. Gene ontology (GO) analysis revealed that several processes are significantly enriched based on DEGs. In particular, response to stress, and especially response to oxidative stress, was over-represented in polyembryonic ovules. Nearly 150 miRNAs, comprising ~90 conserved and ~60 novel miRNAs, were identified in the ovules of either cultivar pair. Only two differentially expressed miRNAs (DEMs) were identified, of which the novel miRN23-5p was repressed whereas the targets accumulated in the polyembryonic ovules. This integrated study on the transcriptional and post-transcriptional regulatory profiles between poly- and monoembryonic citrus ovules provides new insights into the mechanism of NE, which should contribute to revealing the regulatory mechanisms of plant apomixis.

Fragile Sites of 'Valencia' Sweet Orange (Citrus sinensis) Chromosomes Are Related with Active 45s rDNA.

Citrus sinensis chromosomes present a morphological differentiation of bands after staining by the fluorochromes CMA and DAPI, but there is still little information on its chromosomal characteristics. In this study, the chromosomes in 'Valencia' C. sinensis were analyzed by fluorescence in situ hybridization (FISH) using telomere DNA and the 45S rDNA gene as probes combining CMA/DAPI staining, which showed that there were two fragile sites in sweet orange chromosomes co-localizing at distended 45S rDNA regions, one proximally locating on B-type chromosome and the other subterminally locating on D-type chromosome. While the chromosomal CMA banding and 45S rDNA FISH mapping in the doubled haploid line of 'Valencia' C. sinensis indicated six 45S rDNA regions, four were identified as fragile sites as doubled comparing its parental line, which confirmed the cytological heterozygosity and chromosomal heteromorphisms in sweet orange. Furthermore, Ag-NOR identified two distended 45S rDNA regions to be active nucleolar organizing regions (NORs) in diploid 'Valencia' C. sinensis. The occurrence of quadrivalent in meiosis of pollen mother cells (PMCs) in 'Valencia' sweet orange further confirmed it was a chromosomal reciprocal translocation line. We speculated this chromosome translocation was probably related to fragile sites. Our data provide insights into the chromosomal characteristics of the fragile sites in 'Valencia' sweet orange and are expected to facilitate the further investigation of the possible functions of fragile sites.

Prediction and functional analysis of the sweet orange protein-protein interaction network.

BACKGROUND: Sweet orange (Citrus sinensis) is one of the most important fruits world-wide. Because it is a woody plant with a long growth cycle, genetic studies of sweet orange are lagging behind those of other species. RESULTS: In this analysis, we employed ortholog identification and domain combination methods to predict the protein-protein interaction (PPI) network for sweet orange. The K-nearest neighbors (KNN) classification method was used to verify and filter the network. The final predicted PPI network, CitrusNet, contained 8,195 proteins with 124,491 interactions. The quality of CitrusNet was evaluated using gene ontology (GO) and Mapman annotations, which confirmed the reliability of the network. In addition, we calculated the expression difference of interacting genes (EDI) in CitrusNet using RNA-seq data from four sweet orange tissues, and also analyzed the EDI distribution and variation in different sub-networks. CONCLUSIONS: Gene expression in CitrusNet has significant modular features. Target of rapamycin (TOR) protein served as the central node of the hormone-signaling sub-network. All evidence supported the idea that TOR can integrate various hormone signals and affect plant growth. CitrusNet provides valuable resources for the study of biological functions in sweet orange.

L-Ascorbic acid metabolism during fruit development in an ascorbate-rich fruit crop chestnut rose (Rosa roxburghii Tratt).

Chestnut rose (Rosa roxburghii Tratt) is a fruit crop that contains unusually high levels of l-ascorbic acid (AsA; ∼1300 mg 100g(-1) FW). To explore the mechanisms underlying AsA metabolism, we investigated the distribution and abundance of AsA during fruit development. We also analyzed gene expression patterns, enzyme activities, and content of metabolites related to AsA biosynthesis and recycling. AsA first accumulated during late fruit development and continued to accumulate during ripening, with the highest accumulation rate near fruit maturity. The redox state of AsA in fruit was also enhanced during late fruit development, while leaf and other tissues had much lower levels of AsA and the redox state of AsA was lower. In mature fruit, AsA was mainly distributed in the cytoplasm of the mesocarp. Correlation analysis suggested that the gene expression patterns, enzyme activities, and related metabolite concentrations involved in the l-galactose pathway showed relatively high correlations with the accumulation rate of AsA. The gene expression pattern and activity of dehydroascorbate reductase (DHAR, EC 1.8.5.1) correlated strongly with AsA concentration, possibly indicating the crucial role of DHAR in the accumulation of high levels of AsA in chestnut rose fruit. Over expression of DHAR in Arabidopsis significantly increased the reduced AsA content and redox state. This was more effective than over expression of the l-galactose pathway gene GDP-d-mannose-3,5-epimerase (EC 5.1.3.18). These findings will enhance understanding of the molecular mechanisms regulating accumulation of AsA in chestnut rose.

iTRAQ-based quantitative proteomics analysis revealed alterations of carbohydrate metabolism pathways and mitochondrial proteins in a male sterile cybrid pummelo.

Comprehensive and quantitative proteomic information on citrus floral bud is significant for understanding male sterility of the cybrid pummelo (G1+HBP) with nuclear genome of HBP and foreign mitochondrial genome of G1. Scanning electron microscopy and transmission electron microscopy analyses of the anthers showed that the development of pollen wall in G1+HBP was severely defective with a lack of exine and sporopollenin formation. Proteomic analysis was used to identify the differentially expressed proteins between male sterile G1+HBP and fertile type (HBP) with the aim to clarify their potential roles in anther development and male sterility. On the basis of iTRAQ quantitative proteomics, we identified 2235 high-confidence protein groups, 666 of which showed differentially expressed profiles in one or more stages. Proteins up- or down-regulated in G1+HBP were mainly involved in carbohydrate and energy metabolism (e.g., pyruvate dehydrogenase, isocitrate dehydrogenase, ATP synthase, and malate dehydrogenase), nucleotide binding (RNA-binding proteins), protein synthesis and degradation (e.g., ribosome proteins and proteasome subunits). Additionally, the proteins located in mitochondria also showed changed expression patterns. These findings provide a valuable inventory of proteins involved in floral bud development and contribute to elucidate the mechanism of cytoplasmic male sterility in the cybrid pummelo.

Citrus sinensis annotation project (CAP): a comprehensive database for sweet orange genome.

Citrus is one of the most important and widely grown fruit crop with global production ranking firstly among all the fruit crops in the world. Sweet orange accounts for more than half of the Citrus production both in fresh fruit and processed juice. We have sequenced the draft genome of a double-haploid sweet orange (C. sinensis cv. Valencia), and constructed the Citrus sinensis annotation project (CAP) to store and visualize the sequenced genomic and transcriptome data. CAP provides GBrowse-based organization of sweet orange genomic data, which integrates ab initio gene prediction, EST, RNA-seq and RNA-paired end tag (RNA-PET) evidence-based gene annotation. Furthermore, we provide a user-friendly web interface to show the predicted protein-protein interactions (PPIs) and metabolic pathways in sweet orange. CAP provides comprehensive information beneficial to the researchers of sweet orange and other woody plants, which is freely available at http://citrus.hzau.edu.cn/.

Genome-wide characterization and expression analysis of genetic variants in sweet orange.

Heterozyosity is an important feature of many plant genomes, and is related to heterosis. Sweet orange, a highly heterozygous species, is thought to have originated from an inter-species hybrid between pummelo and mandarin. To investigate the heterozygosity of the sweet orange genome and examine how this heterozygosity affects gene expression, we characterized the genome of Valencia orange for single nucleotide variations (SNVs), small insertions and deletions (InDels) and structural variations (SVs), and determined their functional effects on protein-coding genes and non-coding sequences. Almost half of the genes containing large-effect SNVs and InDels were expressed in a tissue-specific manner. We identified 3542 large SVs (>50 bp), including deletions, insertions and inversions. Most of the 296 genes located in large-deletion regions showed low expression levels. RNA-Seq reads and DNA sequencing reads revealed that the alleles of 1062 genes were differentially expressed. In addition, we detected approximately 42 Mb of contigs that were not found in the reference genome of a haploid sweet orange by de novo assembly of unmapped reads, and annotated 134 protein-coding genes within these contigs. We discuss how this heterozygosity affects the quality of genome assembly. This study advances our understanding of the genome architecture of sweet orange, and provides a global view of gene expression at heterozygous loci.

Evaluation of codon biology in citrus and Poncirus trifoliata based on genomic features and frame corrected expressed sequence tags.

Citrus, as one of the globally important fruit trees, has been an object of interest for understanding genetics and evolutionary process in fruit crops. Meta-analyses of 19 Citrus species, including 4 globally and economically important Citrus sinensis, Citrus clementina, Citrus reticulata, and 1 Citrus relative Poncirus trifoliata, were performed. We observed that codons ending with A- or T- at the wobble position were preferred in contrast to C- or G- ending codons, indicating a close association with AT richness of Citrus species and P. trifoliata. The present study postulates a large repertoire of a set of optimal codons for the Citrus genus and P. trifoliata and demonstrates that GCT and GGT are evolutionary conserved optimal codons. Our observation suggested that mutational bias is the dominating force in shaping the codon usage bias (CUB) in Citrus and P. trifoliata. Correspondence analysis (COA) revealed that the principal axis [axis 1; COA/relative synonymous codon usage (RSCU)] contributes only a minor portion (∼10.96%) of the recorded variance. In all analysed species, except P. trifoliata, Gravy and aromaticity played minor roles in resolving CUB. Compositional constraints were found to be strongly associated with the amino acid signatures in Citrus species and P. trifoliata. Our present analysis postulates compositional constraints in Citrus species and P. trifoliata and plausible role of the stress with GC3 and coevolution pattern of amino acid.

The draft genome of sweet orange (Citrus sinensis).

Oranges are an important nutritional source for human health and have immense economic value. Here we present a comprehensive analysis of the draft genome of sweet orange (Citrus sinensis). The assembled sequence covers 87.3% of the estimated orange genome, which is relatively compact, as 20% is composed of repetitive elements. We predicted 29,445 protein-coding genes, half of which are in the heterozygous state. With additional sequencing of two more citrus species and comparative analyses of seven citrus genomes, we present evidence to suggest that sweet orange originated from a backcross hybrid between pummelo and mandarin. Focused analysis on genes involved in vitamin C metabolism showed that GalUR, encoding the rate-limiting enzyme of the galacturonate pathway, is significantly upregulated in orange fruit, and the recent expansion of this gene family may provide a genomic basis. This draft genome represents a valuable resource for understanding and improving many important citrus traits in the future.

Comparative transcript profiling of a male sterile cybrid pummelo and its fertile type revealed altered gene expression related to flower development.

Male sterile and seedless characters are highly desired for citrus cultivar improvement. In our breeding program, a male sterile cybrid pummelo, which could be considered as a variant of male fertile pummelo, was produced by protoplast fusion. Herein, ecotopic stamen primordia initiation and development were detected in this male sterile cybrid pummelo. Histological studies revealed that the cybrid showed reduced petal development in size and width, and retarded stamen primordia development. Additionally, disorganized cell proliferation was also detected in stamen-like structures (fused to petals and/or carpel). To gain new insight into the underlying mechanism, we compared, by RNA-Seq analysis, the nuclear gene expression profiles of floral buds of the cybrid with that of fertile pummelo. Gene expression profiles which identified a large number of differentially expressed genes (DEGs) between the two lines were captured at both petal primordia and stamen primordia distinguishable stages. For example, nuclear genes involved in nucleic acid binding and response to hormone synthesis and metabolism, genes required for floral bud identification and expressed in particular floral whorls. Furthermore, in accordance with flower morphology of the cybrid, expression of PISTILLATA (PI) was reduced in stamen-like structures, even though it was restricted to correct floral whorls. Down-regulated expression of APETALA3 (AP3) coincided with that of PI. These finding indicated that, due to their whorl specific effects in flower development, citrus class-B MADS-box genes likely constituted 'perfect targets' for CMS retrograde signaling, and that dysfunctional mitochondria seemed to cause male sterile phenotype in the cybrid pummelo.

Identification and comparative profiling of miRNAs in an early flowering mutant of trifoliate orange and its wild type by genome-wide deep sequencing.

MicroRNAs (miRNAs) are a new class of small, endogenous RNAs that play a regulatory role in various biological and metabolic processes by negatively affecting gene expression at the post-transcriptional level. While the number of known Arabidopsis and rice miRNAs is continuously increasing, information regarding miRNAs from woody plants such as citrus remains limited. Solexa sequencing was performed at different developmental stages on both an early flowering mutant of trifoliate orange (precocious trifoliate orange, Poncirus trifoliata L. Raf.) and its wild-type in this study, resulting in the obtainment of 141 known miRNAs belonging to 99 families and 75 novel miRNAs in four libraries. A total of 317 potential target genes were predicted based on the 51 novel miRNAs families, GO and KEGG annotation revealed that high ranked miRNA-target genes are those implicated in diverse cellular processes in plants, including development, transcription, protein degradation and cross adaptation. To characterize those miRNAs expressed at the juvenile and adult development stages of the mutant and its wild-type, further analysis on the expression profiles of several miRNAs through real-time PCR was performed. The results revealed that most miRNAs were down-regulated at adult stage compared with juvenile stage for both the mutant and its wild-type. These results indicate that both conserved and novel miRNAs may play important roles in citrus growth and development, stress responses and other physiological processes.

A global view of gene activity at the flowering transition phase in precocious trifoliate orange and its wild-type [Poncirus trifoliata (L.) Raf.] by transcriptome and proteome analysis.

Most of what we know about the molecular genetics of flowering time regulation comes from studies in the model plants. However, little is known about the regulation of flowering transition in perennial species or in species with particular growth habits compared with model plants. Here comparative transcriptome and proteome profiling of spring shoots was performed on an early flowering trifoliate orange mutant (precocious trifoliate orange, Poncirus trifoliata) and its wild-type. A total of 19,215 read-mapped genes were observed in two genotypes by RNA-Seq. Of these, 1450 and 1159 genes were specifically observed in the mutant and wild-type libraries, respectively. There were 355 genes that were expressed differently in the two genotypes. A total of 1664 proteins were identified by the iTRAQ technique, and transcript and protein profiles were parallel across the time course for 50% of the comparisons made, but divergent patterns were also observed, indicative of post-transcriptional events. In addition, a global survey of messenger RNA splicing events identified 16,343 splice junctions among 12,688 genes and showed that alternative 3' splice is the most prevalent form of alternative splicing. We further identify 5698 novel transcribed regions that are not overlapping with annotated citrus transcriptome in two genotypes. Understanding of the regulation of flowering transition in citrus can help in the development of new genetic or management strategies to improve fruit production.

Comparative transcript profiling of gene expression between seedless Ponkan mandarin and its seedy wild type during floral organ development by suppression subtractive hybridization and cDNA microarray.

BACKGROUND: Seedlessness is an important agronomic trait for citrus, and male sterility (MS) is one main cause of seedless citrus fruit. However, the molecular mechanism of citrus seedlessness remained not well explored. RESULTS: An integrative strategy combining suppression subtractive hybridization (SSH) library with cDNA microarray was employed to study the underlying mechanism of seedlessness of a Ponkan mandarin seedless mutant (Citrus reticulata Blanco). Screening with custom microarray, a total of 279 differentially expressed clones were identified, and 133 unigenes (43 contigs and 90 singletons) were obtained after sequencing. Gene Ontology (GO) distribution based on biological process suggested that the majority of differential genes are involved in metabolic process and respond to stimulus and regulation of biology process; based on molecular function they function as DNA/RNA binding or have catalytic activity and oxidoreductase activity. A gene encoding male sterility-like protein was highly up-regulated in the seedless mutant compared with the wild type, while several transcription factors (TFs) such as AP2/EREBP, MYB, WRKY, NAC and C2C2-GATA zinc-finger domain TFs were down-regulated. CONCLUSION: Our research highlighted some candidate pathways that participated in the citrus male gametophyte development and could be beneficial for seedless citrus breeding in the future.

Transcriptional profiling of genes involved in embryogenic, non-embryogenic calluses and somatic embryogenesis of Valencia sweet orange by SSH-based microarray.

Somatic embryogenesis (SE) is a most promising technology that is used for in vitro germplasm conservation and genetic improvement via biotechnological approaches in citrus. Herein, three suppression subtractive hybridization (SSH) libraries were constructed using calluses of Citrus sinensis cv. 'Valencia' to explore the molecular mechanisms that underlie the SE in citrus. A total of 880 unisequences were identified by microarray screening based on these three SSH libraries. Gene ontology analysis of the differentially expressed genes indicated that nucleolus associated regulation and biogenesis processes, hormone signal transduction, and stress factors might be involved in SE. Transcription factors might also play an important role. LEC1/B3 domain regulatory network genes (LEC1, L1L, FUS3, ABI3, and ABI5) were isolated in citrus SE. Some new transcription factors associated with citrus SE, like a B3 domain containing gene and HB4, were identified. To understand the influence of these isolated genes on SE competence, their expression profiles were compared among callus lines of seven citrus cultivars with different SE competence. The expression dynamics suggested that these genes could be necessary for the SE initiation and might play a role in embryogenic competence maintenance in different cultivars. On the basis of gene expression profiles, an overview of major physiological and biosynthesis processes at different developmental stages during citrus SE is presented. For the first time, these data provide a global resource for transcriptional events important for SE in citrus, and the specific genes offer new information for further investigation on citrus SE maintenance and development.