Genomics of the origin and evolution of Citrus.

The genus Citrus, comprising some of the most widely cultivated fruit crops worldwide, includes an uncertain number of species. Here we describe ten natural citrus species, using genomic, phylogenetic and biogeographic analyses of 60 accessions representing diverse citrus germ plasms, and propose that citrus diversified during the late Miocene epoch through a rapid southeast Asian radiation that correlates with a marked weakening of the monsoons. A second radiation enabled by migration across the Wallace line gave rise to the Australian limes in the early Pliocene epoch. Further identification and analyses of hybrids and admixed genomes provides insights into the genealogy of major commercial cultivars of citrus. Among mandarins and sweet orange, we find an extensive network of relatedness that illuminates the domestication of these groups. Widespread pummelo admixture among these mandarins and its correlation with fruit size and acidity suggests a plausible role of pummelo introgression in the selection of palatable mandarins. This work provides a new evolutionary framework for the genus Citrus.

Comparative transcriptomics of wild and commercial Citrus during early ripening reveals how domestication shaped fruit gene expression.

BACKGROUND: Interspecific hybridizations and admixtures were key in Citrus domestication, but very little is known about their impact at the transcriptomic level. To determine the effects of genome introgressions on gene expression, the transcriptomes of the pulp and flavedo of three pure species (citron, pure mandarin and pummelo) and four derived domesticated genetic admixtures (sour orange, sweet orange, lemon and domesticated mandarin) have been analyzed at color break. RESULTS: Many genes involved in relevant physiological processes for domestication, such sugar/acid metabolism and carotenoid/flavonoid synthesis, were differentially expressed among samples. In the low-sugar, highly acidic species lemon and citron, many genes involved in sugar metabolism, the TCA cycle and GABA shunt displayed a reduced expression, while the P-type ATPase CitPH5 and most subunits of the vacuolar ATPase were overexpressed. The red-colored species and admixtures were generally characterized by the overexpression in the flavedo of specific pivotal genes involved in the carotenoid biosynthesis, including phytoene synthase, ζ-carotene desaturase, ß-lycopene cyclase and CCD4b, a carotenoid cleavage dioxygenase. The expression patterns of many genes involved in flavonoid modifications, especially the flavonoid and phenylpropanoid O-methyltransferases showed extreme diversity. However, the most noticeable differential expression was shown by a chalcone synthase gene, which catalyzes a key step in the biosynthesis of flavonoids. This chalcone synthase was exclusively expressed in mandarins and their admixed species, which only expressed the mandarin allele. In addition, comparisons between wild and domesticated mandarins revealed that the major differences between their transcriptomes concentrate in the admixed regions. CONCLUSION: In this work we present a first study providing broad evidence that the genome introgressions that took place during citrus domestication largely shaped gene expression in their fruits.

Regulation of carotenoid and chlorophyll pools in hesperidia, anatomically unique fruits found only in Citrus.

Domesticated citrus varieties are woody perennials and interspecific hybrid crops of global economic and nutritional importance. The citrus fruit "hesperidium" is a unique morphological innovation not found in any other plant lineage. Efforts to improve the nutritional quality of the fruit are predicated on understanding the underlying regulatory mechanisms responsible for fruit development, including temporal control of chlorophyll degradation and carotenoid biosynthesis. Here, we investigated the molecular basis of the navel orange (Citrus sinensis) brown flavedo mutation, which conditions flavedo that is brown instead of orange. To overcome the limitations of using traditional genetic approaches in citrus and other woody perennials, we developed a strategy to elucidate the underlying genetic lesion. We used a multi-omics approach to collect data from several genetic sources and plant chimeras to successfully decipher this mutation. The multi-omics strategy applied here will be valuable in driving future gene discovery efforts in citrus as well as in other woody perennial plants. The comparison of transcriptomic and genomic data from multiple genotypes and plant sectors revealed an underlying lesion in the gene encoding STAY-GREEN (SGR) protein, which simultaneously regulates carotenoid biosynthesis and chlorophyll degradation. However, unlike SGR of other plant species, we found that the carotenoid and chlorophyll regulatory activities could be uncoupled in the case of certain SGR alleles in citrus and thus we propose a model for the molecular mechanism underlying the brown flavedo phenotype. The economic and nutritional value of citrus makes these findings of wide interest. The strategy implemented, and the results obtained, constitute an advance for agro-industry by driving opportunities for citrus crop improvement.

Single-nucleotide mosaicism in citrus: Estimations of somatic mutation rates and total number of variants.

Most of the hundreds of citrus varieties are derived from spontaneous mutations. We characterized the dynamics of single-nucleotide mosaicism in a 36-yr-old clementine (Citrus ×clementina hort. ex Tanaka) tree, a commercial citrus whose vegetative behavior is known in detail. Whole-genome sequencing identified 73 reliable somatic mutations, 48% of which were transitions from G/C to A/T, suggesting ultraviolet (UV) exposure as mutagen. The mutations accumulated in sectorized areas of the tree in a nested hierarchy determined by the branching pattern, although some variants detected in the basal parts were also found in the new growth and were fixed in some branches and leaves of much younger age. The estimate of mutation rates in our tree was 4.4 × 10-10 bp-1 yr-1 , a rate in the range reported in other perennials. Assuming a perfect configuration and taking advantage of previous counts on the number of total leaves of typical clementine trees, these mutation determinations allowed to estimate for the first time the total number of variants present in a standard adult tree (1,500-5,000) and the somatic mutations generated in a typical leaf flush (0.92-1.19). From an evolutionary standpoint, the sectoral distribution of somatic mutations and the habit of periodic foliar renewal of long-lived plants appear to increase genetic heterogeneity and, therefore, the adaptive role of somatic mutations reducing the mutational load and providing fitness benefits.

Transcriptome analysis of the pulp of citrus fruitlets suggests that domestication enhanced growth processes and reduced chemical defenses increasing palatability.

To identify key traits brought about by citrus domestication, we have analyzed the transcriptomes of the pulp of developing fruitlets of inedible wild Ichang papeda (Citrus ichangensis), acidic Sun Chu Sha Kat mandarin (C. reticulata) and three palatable segregants of a cross between commercial Clementine (C. x clementina) and W. Murcott (C. x reticulata) mandarins, two pummelo/mandarin admixtures of worldwide distribution. RNA-seq comparison between the wild citrus and the ancestral sour mandarin identified 7267 differentially expressed genes, out of which 2342 were mapped to 117 KEGG pathways. From the remaining genes, a set of 2832 genes was functionally annotated and grouped into 45 user-defined categories. The data suggest that domestication promoted fundamental growth processes to the detriment of the production of chemical defenses, namely, alkaloids, terpenoids, phenylpropanoids, flavonoids, glucosinolates and cyanogenic glucosides. In the papeda, the generation of energy to support a more active secondary metabolism appears to be dependent upon upregulation of glycolysis, fatty acid degradation, Calvin cycle, oxidative phosphorylation, and ATP-citrate lyase and GABA pathways. In the acidic mandarin, downregulation of cytosolic citrate degradation was concomitant with vacuolar citrate accumulation. These changes affected nitrogen and carbon allocation in both species leading to major differences in organoleptic properties since the reduction of unpleasant secondary metabolites increases palatability while acidity reduces acceptability. The comparison between the segregants and the acidic mandarin identified 357 transcripts characterized by the occurrence in the three segregants of additional downregulation of secondary metabolites and basic structural cell wall components. The segregants also showed upregulation of genes involved in the synthesis of methyl anthranilate and furaneol, key substances of pleasant fruity aroma and flavor, and of sugar transporters relevant for sugar accumulation. Transcriptome and qPCR analysis in developing and ripe fruit of a set of genes previously associated with citric acid accumulation, demonstrated that lower acidity is linked to downregulation of these regulatory genes in the segregants. The results suggest that the transition of inedible papeda to sour mandarin implicated drastic gene expression reprograming of pivotal pathways of the primary and secondary metabolism, while palatable mandarins evolved through progressive refining of palatability properties, especially acidity.

Shaping the biology of citrus: I. Genomic determinants of evolution.

We performed genomic analyses on wild species of the genus Citrus to identify major determinants of evolution. The most notable effect occurred on the pathogen-defense genes, as observed in many other plant genera. The gene space was also characterized by changes in gene families intimately related to relevant biochemical properties of citrus fruit, such as pectin modifying enzymes, HDR (4-hydroxy-3-methylbut-2-enyl diphosphate reductase) genes, and O-methyltransferases. Citrus fruits are highly abundant on pectins and secondary metabolites such as terpenoids and flavonoids, the targets of these families. Other gene types under positive selection, expanded through tandem duplications and retained as triplets from whole genome duplications, codified for purple acid phosphatases and MATE-efflux proteins. Although speciation has not been especially rapid in the genus, analyses of selective pressure at the codon level revealed that the extant species evolved from the ancestral citrus radiation show signatures of pervasive adaptive evolution and is therefore potentially responsible for the vast phenotypic differences observed among current species.

Shaping the biology of citrus: II. Genomic determinants of domestication.

We performed genomic analyses on species and varieties of the genus Citrus to identify several determinants of domestication, based on the pattern of pummelo [Citrus maxima (Burr. f) Merr] and mandarin (Citrus reticulata Blanco) admixture into the ancestral genome, as well as population genetic tests at smaller scales. Domestication impacted gene families regulating pivotal components of citrus flavor (such as acidity) because in edible mandarin varieties, chromosome areas with negative Tajimas values were enriched with genes associated with the regulation of citric acid. Detection of sweeps in edible mandarins that diverged from wild relatives indicated that domestication reduced chemical defenses involving cyanogenesis and alkaloid synthesis, thus increasing palatability. Also, a cluster of SAUR genes in domesticated mandarins derived from the pummelo genome appears to contain candidate genes controlling fruit size. Similarly, conserved stretches of pure mandarin areas were likely important as well for domestication, as, for example, a fragment in chromosome 1 that is involved in the apomictic reproduction of most edible mandarins. Interestingly, our results also support the hypothesis that various genes subject to selective pressure during evolution or derived from whole genome duplication events later became potential targets of domestication.

Reprogramming of Retrotransposon Activity during Speciation of the Genus Citrus.

Speciation of the genus Citrus from a common ancestor has recently been established to begin ∼8 Ma during the late Miocene, a period of major climatic alterations. Here, we report the changes in activity of Citrus LTR retrotransposons during the process of diversification that gave rise to the current Citrus species. To reach this goal, we analyzed four pure species that diverged early during Citrus speciation, three recent admixtures derived from those species and an outgroup of the Citrus clade. More than 30,000 retrotransposons were grouped in ten linages. Estimations of LTR insertion times revealed that retrotransposon activity followed a species-specific pattern of change that could be ascribed to one of three different models. In some genomes, the expected pattern of gradual transposon accumulation was suddenly arrested during the radiation of the ancestor that gave birth to the current Citrus species. The individualized analyses of retrotransposon lineages showed that in each and every species studied, not all lineages follow the general pattern of the species itself. For instance, in most of the genomes, the retrotransposon activity of elements from the SIRE lineage reached its highest level just before Citrus speciation, while for Retrofit elements, it has been steadily growing. Based on these observations, we propose that Citrus retrotransposons may respond to stressful conditions driving speciation as a part of the genetic response involved in adaptation. This proposal implies that the evolving conditions of each species interact with the internal regulatory mechanisms of the genome controlling the proliferation of mobile elements.