# The Distribution of Coumarins and Furanocoumarins in Citrus Species Closely Matches Citrus Phylogeny and Reflects the Organization of Biosynthetic Pathways

Audray Dugrand-Judek, Alexandre Olry, Alain Hehn, Gilles Costantino, Patrick Ollitrault, Yann Froelicher, Frédéric Bourgaud

<figure><img src="image_1.png" /></figure>

## To cite this version:

Audray Dugrand-Judek, Alexandre Olry, Alain Hehn, Gilles Costantino, Patrick Ollitrault, et al.. The Distribution of Coumarins and Furanocoumarins in Citrus Species Closely Matches Citrus Phylogeny and Reflects the Organization of Biosynthetic Pathways. PLoS ONE, 2015, 10 (11), e01475936. doi:10.1371/journal.pone.0142757. This version is licensed under a Creative Commons Attribution 4.0 International License.

HAL Id: hal-01475936

https://hal.univ-lorraine.fr/hal-01475936v1

Submitted on 28 Feb 2017

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

RESEARCH ARTICLE

# The Distribution of Coumarins and Furanocoumarins in *Citrus* Species Closely Matches *Citrus* Phylogeny and Reflects the Organization of Biosynthetic Pathways

Audray Dugrand-Judek¹,², Alexandre Olry¹,², Alain Hehn¹,², Gilles Costantino³, Patrick Ollitrault⁴,⁵, Yann Froelicher⁶*, Frédéric Bourgaud¹,²*

1 Université de Lorraine, UMR 1121 Laboratoire Agronomie et Environnement Nancy-Colmar, 2 avenue de la forêt de Haye, TSA 40602, 54518, Vandœuvre-lès-Nancy, France, 2 INRA, UMR 1121 Université de Lorraine, UMR 1121 Laboratoire Agronomie et Environnement Nancy-Colmar, 2 avenue de la forêt de Haye, TSA 40602, 54518, Vandœuvre-lès-Nancy, France, 3 INRA, UMR AGAP, Station INRA, F-20230, San Giuliano, France, 4 Centro de Protección Vegetal y Biotecnología, Instituto Valenciano de Investigaciones Agrarias (Ivia), 46113 Moncada, Valencia, Spain, 5 CIRAD, UMR AGAP, Station de Roujol, 97170 Petit-Bourg, Guadeloupe, France, 6 CIRAD, UMR AGAP, Station INRA, F-20230, San Giuliano, France

* frederic.bourgaud@univ-lorraine.fr (FB); froelicher@cirad.fr (YF)

## Abstract

*Citrus* plants are able to produce defense compounds such as coumarins and furanocoumarins to cope with herbivorous insects and pathogens. In humans, these chemical compounds are strong photosensitizers and can interact with medications, leading to the "grapefruit juice effect". Removing coumarins and furanocoumarins from food and cosmetics imply additional costs and might alter product quality. Thus, the selection of *Citrus* cultivars displaying low coumarin and furanocoumarin contents constitutes a valuable alternative. In this study, we performed ultra-performance liquid chromatography coupled with mass spectrometry analyses to determine the contents of these compounds within the peel and the pulp of 61 *Citrus* species representative of the genetic diversity all *Citrus*. Generally, *Citrus* peel contains larger diversity and higher concentrations of coumarin/furanocoumarin than the pulp of the same fruits. According to the chemotypes found in the peel, *Citrus* species can be separated into 4 groups that correspond to the 4 ancestral taxa (pummelos, mandarins, citrons and papedas) and extended with their respective secondary species descendants. Three of the 4 ancestral taxa (pummelos, citrons and papedas) synthesize high amounts of these compounds, whereas mandarins appear practically devoid of them. Additionally, all ancestral taxa and their hybrids are logically organized according to the coumarin and furanocoumarin pathways described in the literature. This organization allows hypotheses to be drawn regarding the biosynthetic origin of compounds for which the biogenesis remains unresolved. Determining coumarin and furanocoumarin contents is also helpful for hypothesizing the origin of *Citrus* species for which the phylogeny is presently not firmly established. Finally, this work also notes favorable hybridization schemes that will lead to low coumarin and furanocoumarin contents, and we propose to

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Dugrand-Judek A, Olry A, Hehn A, Costantino G, Ollitrault P, Froelicher Y, et al. (2015) The Distribution of Coumarins and Furanocoumarins in *Citrus* Species Closely Matches *Citrus* Phylogeny and Reflects the Organization of Biosynthetic Pathways. PLoS ONE 10(11): e0142757. doi:10.1371/journal.pone.0142757

**Editor:** Chunxian Chen, USDA/ARS, UNITED STATES

**Received:** July 21, 2015

**Accepted:** October 23, 2015

**Published:** November 11, 2015

**Copyright:** © 2015 Dugrand-Judek et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** The authors have no support or funding to report.

**Competing Interests:** The authors have declared that no competing interests exist.
select mandarins and Ichang papeda as *Citrus* varieties for use in creating species devoid of these toxic compounds in future breeding programs.

**Introduction**

Citrus originate from the tropical and subtropical regions of Southeast Asia and their culture was initiated in India and China during the first millenary BC [1]. Currently, citrus are cultivated between the latitudes 40°N and 40°S in tropical and subtropical regions [2]. In 2013, citrus crops represented the most important fruit produced in the world, with more than 135 million tons [3]. Citrus belong to the Rutaceae family, which is composed of approximately 160 genera and 1900 species [4]; true citrus belong to the Aurantioideae subfamily, the Citreae tribe and the Citrineae subtribe. True citrus include 6 genera: *Citrus*, *Fortunella*, *Poncirus*, *Microcitrus*, *Eremocitrus* and *Clymenia* [5]. Citrus taxonomy remains controversial due to a long cultivation history, complex reproductive biology and somatic bud mutation. Swingle and Reece (1967) [5] and Tanaka (1977) [6], which recognize 16 and 162 species, respectively, remain the two major classification systems currently used. However, this last decade, molecular analyses have provided decisive information regarding *Citrus* domestication and the relations between various cultivated species of *Citrus* [7–14]. Four ancestral taxa, *Citrus medica* L. (citron), *Citrus reticulata* Blanco (mandarin), *Citrus maxima* (Burm.) Merr. (pummelo), and *Citrus micrantha* Wester (papeda), have been identified as the ancestors of all cultivated Citrus [8,10,12,14]. Among the four basic horticultural groups, no evidence of interspecific introgression was found in pummelos, citrons, and *C. micrantha*. In contrast, evidence of introgression by *C. maxima* was found in sweet mandarin by sequencing and resequencing [13,14]. However, ‘Cleopatra’, ‘Sunki’ and ‘Shekwasha’ mandarin, which share the acidic mandarin cytoplasm defined by Froelicher et al. (2011) [10] from mitochondrial data and confirmed by Carbonell-Caballero et al. (2015) [15] with chloroplast data, present no evidence of introgression [14]. The secondary species *Citrus sinensis* (L.) Osb. (sweet orange), *Citrus aurantium* L. (sour orange), *Citrus paradisi* Macf. (grapefruit), *Citrus limon* (L.) Burm. (lemon), and *Citrus aurantifolia* (Christm.) Swing. (lime) arose from the hybridization of the 4 ancestral taxa [8,10,12,14]. In fact, *C. sinensis*, *C. aurantium*, *C. paradisi*, and tangors had a 2 ancestral taxa admixture structure (*C. reticulata* and *C. maxima*) with variable contributions. The origins of *C. limon* and *C. aurantifolia* (Swingle and Reece taxonomy) are more complex and come from 2 or 3 taxa hybridization. Similarly, *C. aurantifolia* originates from different taxa hybridization (2 to 4 taxa) composed systematically of citron and papeda taxa. *Citrus amblycarpa* secondary species displayed approximately 50% contributions from *C. reticulata* and *C. micrantha* [14].

Coumarins (benzo-α-pyrones) constitute a class of secondary metabolites commonly found in higher plants [16]. In a limited number of plant families, such as Rutaceae, Moraceae, Apiaceae and Fabaceae, umbelliferone (7-hydroxycoumarin), a ubiquitous coumarin of higher plants, can undergo subsequent biochemical modifications corresponding to a prenylation step at C6 or C8 followed by the closure of a furan ring. These modifications generate furanocoumarins, a new subclass of compounds involved in plant defense against pathogens [16] and, more generally, in plant environmental adaptation [16,17]. Citrus, as members of the Rutaceae, can synthesize both coumarins and furanocoumarins and are rich in dimethylallylated and/or geranylated compounds (see Fig 1) such as bergamottin (5-geranyloxypsoralen), aurapten (7-geranyloxycoumarin) or imperatorin (8-dimethylallyloxypsoralen).
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Schematic representation of the coumarin and furanocoumarin pathways. Solid lines indicate direct enzymatic steps; dashed lines represent multiple enzymatic steps. Coumarins are framed in orange, furanocoumarins of the bergapten cluster are framed in blue, furanocoumarins of the xanthotoxin cluster are framed in red, and furanocoumarins of the isopimpinellin cluster are framed in green.</figcaption></figure>

doi:10.1371/journal.pone.0142757.g001

In parallel to their ecological functions in plants, coumarins and furanocoumarins can be deleterious for humans. These compounds are potential photosensitizers that can cause severe phytophotodermatitis after either skin contact [18] or ingestion [19] followed by sun UV exposure. This photosensitization property is a notable problem with citrus essential oils such as bergamot oil because they are extensively used in perfumes and cosmetics. In addition, furanocoumarin monomers [20] and dimers such as paradisins [21] are known to inhibit cytochrome P450 (CYP) from different families [22,23]. Notably, interactions between furanocoumarins and human CYP3A4, a CYP found in intestinal enterocytes and in liver hepatocytes, have generated many concerns because this P450 (i) represents half of all human P450 enzyme pools
and (ii) plays a major role in the metabolism of drugs altered by oxidation [24]. CYP3A4 furanocoumarin-mediated inhibition increases the bioavailability of possibly more than 85 medicines, among which 43 can lead to severe adverse effects, such as heart rhythm disorders or respiratory depression [25]. This interaction was first discovered with grapefruit juice consumption in the beginning of the 1990s [26] and is often called the "grapefruit juice effect". Since then, grapefruit consumption steadily declined from approximately 2.4 million tons in 1991 to approximately 1.75 million tons in 2009 in the developed countries aware of this effect [27]. Yet, citrus fruits also display beneficial effects on health, such as antioxidant properties [28], neuroprotective effects [29] or the reduction of cardiovascular risk [30].

In relation to the undesirable effects of furanocoumarins in humans, different studies have aimed at removing these toxicants from grapefruit juices or citrus essential oils. For example, furanocoumarins can be removed by different extraction and absorption processes [31], degraded by heat treatment [32] or ultraviolet irradiation [33], or removed by the use of autoclaved fungi [34]. However, these treatments imply additional costs and might alter the qualities of grapefruit juice or essential oils. Because the furanocoumarin content in plants is determined by the occurrence of specialized genes [35–38] and is therefore highly heritable, an alternative solution for lowering the furanocoumarin concentration in citrus would consist of developing new breeding programs that would allow access to new citrus varieties with low or zero furanocoumarin content.

Thus far, several studies have attempted to identify and quantitate coumarins and furanocoumarins in citrus tissue [39–42]. However, none of these studies was sufficiently comprehensive regarding the number of chemical compounds and/or citrus varieties that were studied. Thus, a comprehensive study regarding coumarin and furanocoumarin distribution among citrus varieties in terms of both the diversity of chemicals and quantities of compounds synthesized is necessary for building up future breeding projects.

Because our objective was to perform quantitative analyses of already known coumarins and furanocoumarins in citrus, we limited our study to the identification and quantitation of compounds for which authentic standards were available. Under these premises, the present study reports the chemical diversity represented by 27 coumarins and furanocoumarins in 61 *Citrus* varieties from 21 different species (according to Swingle's classification) that are representative of the genetic diversity of all *Citrus*. Our study focused on the sole *Citrus* genus, as it is the largest of the 6 true citrus genera in terms of species number, the most cultivated and the most used by the food and perfume industries. The 4 ancestral taxa were represented by 7 pumelos, 11 mandarins (8 sweet mandarins and 3 acidic mandarins), 3 citrons and 3 papedas. Secondary species were also investigated, with 4 sweet oranges, 4 sour oranges, 3 grapefruits, 4 small mandarin hybrids, 7 lemons, 10 limes, 1 bergamot and 4 other hybrids: Halimii (Mountain), Yuzu, Latipes (Khasi) and Amblycarpa (Nasnaran mandarin). Both citrus peel and pulp were analyzed because the pulp is used for human consumption and because the peel can be used as a flavoring or fragrance ingredient in food and perfumes.

The coumarin and furanocoumarin contents found in citrus peel and pulp samples are compared, and the chemical diversity is discussed with respect to the genetic diversity of the studied varieties and to the organization of the biosynthetic pathways. With reference to the *Citrus* diversity obtained from molecular markers, we highlight the genetic crossings that may have resulted in low coumarin and furanocoumarin content varieties. Finally, the *Citrus* varieties with low coumarin and furanocoumarin content identified in this study will constitute invaluable genetic resources for future breeding programs, promoting new citrus species devoid of these toxic compounds.
# Materials and Methods

## Plant material

Citrus fruits (except *C. micrantha*) were harvested at the *Citrus* Biological Resource Center (CRB CITRUS, INRA-CIRAD, France) at San Giuliano (Corsica, France). *C. micrantha* fruits were harvested at Bachès nursery located in Eus (France). This study included 61 *Citrus* varieties representative of the *Citrus* diversity described in Table 1 and in S1 Table. For each *Citrus* variety, 5 fruits were harvested from the tree.

## Sample preparation

This study was performed on citrus peel and pulp samples. All fruits were harvested and sampled as described previously [39]. Briefly, citrus peel samples consisted of the flavedo and albedo tissues, without exceeding 4 mm for the albedo part, as its thickness greatly varies from one species to the other. Citrus pulp samples only consisted of the juice vesicles, without the walls separating the segments, which were excluded. The only difference with the citrus peel samples is that the citrus pulp was ground for only 2 minutes in the ball mill. Fresh samples were weighed, frozen in liquid nitrogen and stored in a freezer at -80°C until extraction was performed.

## Extraction of compounds

The extraction of the compounds from the peel and pulp samples was identical to the method described by Dugrand *et al.* (2013) [39]. In brief, frozen samples were lyophilized and ground in a ball mill. Two extractions of the compounds were conducted in 80:20 methanol HPLC-grade/water solutions. The citrus extracts obtained were evaporated and then resuspended in a 75:25 methanol HPLC-grade/ultrapure water solution.

## UPLC-MS analyses and data processing

UPLC-MS analyses and data processing were performed with the same equipment and methods described by Dugrand *et al.* (2013) [39]. In short, the UPLC separation of the coumarins and fur-anocoumarins was performed on a C18 reversed-phase column thermostated at 40°C. The mobile phase consisted of 0.1% formic acid in ultrapure water (A) and 0.1% formic acid in meth-anol (B). The elution gradient was (A:B; v/v): 90:10 at 0 min, 80:20 at 0.74 min, 40:60 at 5.88 min, 10:90 at 10 min, 0:100 between 12 and 16 min, and 90:10 from 16.01 to 20 min. The total analysis lasted 20 min. The flow rate was set at 0.2 mL/min and the injection volume was 3 µL. Concern-ing the MS system, we used a dual ion source (DUIS). The inlet, desolvation line and heating block temperatures were respectively set at 350, 250 and 400°C. The voltage of the capillary was set at 4.5 kV. The detection of the compounds was performed in single ion monitoring mode (SIM). Data were acquired and processed on LabSolution software version 5.52 sp2 (Shimadzu).

## Method validation

Method validation has been performed for citrus peel samples in our previous work [39] but has not been performed for citrus pulp samples. Full validation of a method includes checking the linearity, limits of detection (LOD) and quantitation (LOQ), specificity, accuracy, precision and robustness. However, several of these parameters were unchanged between the 2 types of samples studied here (peel vs. pulp) and did not need to be repeated. The measurements that did not need to be repeated include the following: (i) the linearity, which concerns the calibra-tion solutions and not the samples themselves; (ii) the accuracy, which was performed as required in a complex matrix devoid of the studied compounds (here, on tomato leaf samples
<table><caption>Table 1. List of the accessions investigated in the study and their phylogenetic constitution.</caption><thead><tr><th></th><th>Swingle's classification system</th><th>Phylogenetic constitution</th><th>Number of accessions</th></tr></thead><tbody><tr><td colspan="4"><strong>Horticultural group of ancestral species</strong></td></tr><tr><td>Papeda</td><td>C. micrantha Wester</td><td>C. micrantha</td><td>2</td></tr><tr><td>Citron</td><td>C. medica L.</td><td>C. medica</td><td>3</td></tr><tr><td>Pummelo</td><td>C. maxima (Burm.) Merr.</td><td>C. maxima</td><td>7</td></tr><tr><td>Mandarin</td><td>C. reticulata Blanco</td><td>C. reticulata</td><td>11</td></tr><tr><td colspan="4"><strong>Horticultural group of secondary species</strong></td></tr><tr><td>Small mandarin hybrid</td><td>C. reticulata x C. sinensis<br/>C. reticulata x C. sinensis?<br/>C. reticulata Blanco<br/>(C. reticulata x C. sinensis) x C. tangerina</td><td>C. reticulata x ((C. maxima x C. reticulata) x C. reticulata)<br/>C. reticulata x ((C. maxima x C. reticulata) x C. reticulata)?<br/>?</td><td>1<br/>1<br/>1<br/>1</td></tr><tr><td>Sweet orange</td><td>C. sinensis (L.) Osbeck</td><td>(C. maxima x C. reticulata) x C. reticulata</td><td>4</td></tr><tr><td>Grapefruit</td><td>C. paradisi Macf.</td><td>C. maxima x ((C. maxima x C. reticulata) x C. reticulata)</td><td>3</td></tr><tr><td>Sour orange</td><td>C. aurantium L.</td><td>C. maxima x C. reticulata</td><td>4</td></tr><tr><td>Lemon</td><td>C. limon (L.) Burm. f.</td><td>(C. maxima x C. reticulata) x C. medica<br/>(C. maxima x C. reticulata) x C. medica or ((C. maxima x C. reticulata)x C. reticulata) x C. medica<br/>C. reticulata x C. medica</td><td>2<br/>1<br/>4</td></tr><tr><td>Lime</td><td>C. aurantifolia (Christm.) Swing.</td><td>C. micrantha x C. medica<br/>2 (C. micrantha x C. medica) x C. medica<br/>((C. maxima x C. reticulata) x C. medica) x 2 (C. micrantha x C. medica)<br/>(C. maxima x C. reticulata) x C. medica</td><td>4<br/>1<br/>2<br/>2<br/>2</td></tr><tr><td>Bergamot</td><td>C. aurantifolia var Bergamia</td><td>((C. maxima x C. reticulata) x C. medica) x (C. maxima x C. reticulata)</td><td>1</td></tr><tr><td>Amblycarpa</td><td>C. reticulata Blanco hybrid</td><td>C. micrantha x C. reticulata</td><td>1</td></tr><tr><td>Mandarin x citron hybrid</td><td>C. reticulata Blanco var. austera Swing. x C. medica L.</td><td>C. reticulata x C. medica</td><td>1</td></tr><tr><td>Ichang</td><td>C. ichangensis Swing.</td><td>?</td><td>1</td></tr><tr><td>Halimii</td><td>C. halimii B.C. Stone (Tan. Class.)</td><td>?</td><td>1</td></tr><tr><td>Yuzu</td><td>C. ichang austera ex. Tan. hybrid</td><td>?</td><td>1</td></tr><tr><td>Latipes</td><td>C. latipes (Swing.) Tan.</td><td>?</td><td>1</td></tr></tbody></table>

doi:10.1371/journal.pone.0142757.t001

because they do not contain coumarins and furanocoumarins); (iii) the precision; and (iv) the robustness. These latter elements only depended on the equipment. In contrast, the LOD, LOQ, and specificity depended on the matrix containing the chemical compounds and had to be checked again on citrus pulp samples. The LOD, LOQ and specificity were determined using a lime (var. Tahiti, *C. aurantifolia* (Christm.) Swing.) because it contains both high diversity and quantities of coumarins and furanocoumarins in the pulp.

## Limits of detection and quantitation in citrus pulp

The LOD and LOQ of a lime pulp sample were determined using the same method as Dugrand and colleagues (2013) [39]. Because lime pulp does not contain all 27 of the investigated
compounds, this sample was enriched in coumarins and furanocoumarins to confirm the ability to measure the LOD and LOQ for all the studied compounds. The data are presented in S2 Table. The LOD and LOQ values were generally low. Only 6',7'-dihydroxybergamottin, 8-geranyloxypsoralen, byakangelicin and cnidicin had higher limit values. Notably, the LOD and LOQ values in the citrus peel were presented in our previous work [39].

## Specificity

A lime pulp sample was prepared and injected into the UPLC-MS system to check its coumarin and furanocoumarin contents. Based on this analysis, we decided to determine the specificity of 6 quantifiable compounds: aurapten, bergapten, oxypeucedanin, oxypeucedanin hydrate, 8-geranyloxypsoralen and isopimpinellin. For this test, 4 additional lime pulp samples (from the same pulp powder) were spiked with these 6 coumarins and furanocoumarins at 25, 50, 75 and 100% of their initial concentrations. S2 Table shows the equations linking the compounds concentrations, the percentages of the initial concentrations added and the linear correlation coefficients. No matrix effect was observed for the coumarins and furanocoumarins in pulp because the linear correlation coefficients of the equations were above 0.99.

## Statistical analyses

The data were processed using R statistical software (http://www.R-project.org). The gplots R package and the heatmap function were used to build the heatmaps. The qualitative data (0 for the absence of a given compound and 1 for the presence of this compound) were used to create a profile matrix of the 27 coumarins and furanocoumarins detected by UPLC-MS within the 61 *Citrus* varieties. The quantitative data (means of 5 biological replicates expressed as mg.kg⁻¹ fresh weight) were used to build (1) Neighbor Joining (NJ) analyses based on the Range Manhattan dissimilarity index (R pvclust package and pvclust function), (2) Principal components analyses (PCAs) where data were transformed as log₁₀(1 + x) (R Ade4 package and dudi.pca function) and where accessions of the ancestral taxa were used as active individuals to define the PCA axes and accessions of the secondary species were additional individuals.

## Results

### General coumarin and furanocoumarin diversity in *Citrus*

In total, 5 coumarins and 21 furanocoumarins were quantitatively determined in the peel and pulp samples of 61 *Citrus* varieties. In addition, 5-geranyloxy-7-methoxycoumarin was unequivocally determined when present in a sample but could not be quantitated due to the insufficient amount of standard available to establish calibration curves. The coumarin and furanocoumarin contents in the citrus peel and pulp samples are presented in Figs 2 and 3, respectively, and are highly variable between the 61 *Citrus* varieties analyzed (see Supporting information, S1 File for complete quantitative results). Generally, citrus fruit peel has broader diversity of compounds than the corresponding pulp, as shown by the total numbers of coumarins and furanocoumarins found in the peel (649 occurrences) and pulp (561 occurrences) samples of the 61 studied *Citrus* varieties. Similar to the diversity of compounds, higher compound concentrations were generally found in peel than in pulp (Figs 2 and 3). Exceptions to this latter rule, when the pulp appears more concentrated than the peel, were observed in the cases of Deep Red and Pink pummelo and Bouquettier de Nice and Chinotto sour oranges for which the pulp is exceptionally rich in furanocoumarins from the bergapten group (Tables A and G in S1 File). Another notable exception was bergamottin, a bergapten-geranylated derivative, which is regularly found more concentrated in the pulp of pummelos (Deep Red,
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Coumarin and furanocoumarin quantities (mmol kg⁻¹ fresh weight) in the peel extracts of the 61 *Citrus* species investigated. Coumarins are represented in orange; furanocoumarins of the bergapten cluster, in blue; furanocoumarins of the xanthotoxin cluster, in red; and furanocoumarins of the isopimpinellin cluster, in green. Ancestral taxa and secondary species are highlighted in the bottom of the graph. Sweet mandarins are presented in bright red, and acidic mandarins are illustrated in deep red.</figcaption></figure>

doi:10.1371/journal.pone.0142757.g002

Chandler) and their hybrid grapefruits (Duncan, Marsh) and sour oranges (Maroc, Bouquettier de Nice, Granito) (Tables A, E and G in S1 File). Based on the results presented in Figs 2 and 3, the highest producers of coumarins and furanocoumarins are 2 papedas: Micrantha (10518 and 1118 mmol kg⁻¹ fresh weight in peel and pulp, respectively) and Combava (6136 and 1189 mmol kg⁻¹ fresh weight in peel and pulp, respectively). Citrons and some of their descendants (particularly Tahiti, Bears, Alemow, Mexican, Excelsa, Giant and Key limes), pummelos and
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Coumarin and furanocoumarin quantities (mmol kg⁻¹ fresh weight) in the pulp extracts of the 61 *Citrus* species investigated. Coumarins are represented in orange; furanocoumarins of the bergapten cluster, in blue; furanocoumarins of the xanthotoxin cluster, in red; and furanocoumarins of the isopimpinellin cluster, in green. Ancestral taxa and secondary species are highlighted in the bottom of the graph. Sweet mandarins are presented in bright red, and acidic mandarins are illustrated in deep red.</figcaption></figure>

doi:10.1371/journal.pone.0142757.g003

some of their descendants (grapefruit, sour oranges and bergamot) and Nasnaran also synthesize important amounts of coumarins and furanocoumarins with respect to their peel and pulp contents. Other citron descendants (particularly those for which Micrantha is not related to, e. g., Eureka and Meyer lemons and Marrakech, Palestine, Rangpur and Brazil sweet limes) contain intermediate levels in the peel, compared with extremely low levels in the pulp. Other varieties such as mandarins, small mandarin hybrids, sweet oranges, some limes, some lemons,
Ichang papeda and Mountain display low amounts of coumarins and furanocoumarins in both their peel and pulp. As examples, the peel/pulp of Bendiguangju satsuma, Murcott tangor and Shamouti sweet orange present extremely low coumarin and furanocoumarin contents, with 1.85/0.97, 2.10/1.24, and 2.02/1.15 mmol kg⁻¹ fresh weight, respectively.

Heatmap representations associated with hierarchical ascendant classifications (HACs) were chosen to establish classifications and correlations between the chemotypes and *Citrus* varieties. Heatmaps were constructed based on the sole presence or absence (0 for the absence of the compound and 1 for the detection of the compound) of each of the coumarins/furanocoumarins (5-geranyloxy-7-methoxycoumarin included). This analysis allowed the chemotypes observed for the peel (Fig 4) and for the pulp (Fig 5) in all *Citrus* varieties to be summarized.

Based on the largest compound diversity found in peel ( heatmap, Fig 4A ), coumarins and furanocoumarins were classified by HACs into 4 different subgroups ( dendrogram, Fig 4B ).

Group I is composed of 3 coumarins (aurapten, limettin, and 5-geranyloxy-7-methoxycou-
marin), bergapten and all its dimethylallylated derivatives (isoimperatorin, oxypeucedanin,
and oxypeucedanin-hydrate), and the geranylated compounds bergamottin and isopimpinell-
lin. This group includes the most widespread coumarins/furanocoumarins in the *Citrus* species
studied, as exemplified by oxypeucedanin, which was ubiquitously found in the peel of all 61
investigated varieties.

Group II is composed of xanthotoxin dimethylallylated compounds (imperatorin, heracle-
nin, and heraclenol) associated with isopimpinellin dimethylallylated derivatives, where
DMAPP is attached to the oxygen bound to C8 (phellopterin, byakangelicol, byakangelicin,
and cnidicin).

Group III includes parent compounds such as psoralen, bergaptol, xanthotoxol and umbel-
liferone, which subsequently generate additional diversified compounds in the coumarin/fura-
nocoumarin pathways. Group III also contains cnidilin, a double DMAPP-derivative of
isopimpinellin, and 8-geranyloxypsoralen. Group III includes some of the most scarcely found
coumarins/furanocoumarins among the studied *Citrus* varieties.

Group IV is composed of two coumarins, osthol and epoxyaurapten, which are associated with two bergapten geranylated compounds (epoxybergamottin and 6-7-dihydroxybergamot- tin). Compounds from group IV are less common in *Citrus* varieties and are restricted to some horticultural groups.

The distribution of coumarins and furanocoumarins in citrus pulp ( heatmap Fig 5A ) signifi-
cantly differs from that previously described in the peel ( heatmap, Fig 4A ). In pulp, coumarin/
furanocoumarin classifications are only limited to 3 groups ( dendrogram, Fig 5B ), denoted as
I', II', and III' to avoid confusion with the groups described previously in Fig 5.

Group I' is composed of the coumarins osthol and epoxyaurapten and the furanocoumarins psoralen, the 5 other compounds from the bergapten cluster (prenylated and geranylated forms) and 8-geranyloxypsoralen. Most of the compounds from group I' are geranylated forms.

Group II' contains the coumarin umbelliferone, the furanocoumarin oxypeucedanin
hydrate, the compounds from the xanthotoxin cluster (except 8-geranyloxypsoralen) and the
prenylated compounds of the isopimpinellin cluster. These coumarins and furanocoumarins
from group II' are primarily dimethylallylated compounds.

Group III' includes compounds that can be found in many *Citrus* fruits: the coumarins
5-geranyloxy-7-methoxycoumarin, aurapten and limettin and the furanocoumarins oxypeuce-
danin, bergapten and isopimpinellin. Compounds from group III' are the most widespread in
citrus pulp samples.
<figure><img src="image_8.png" /><figcaption>**Fig 4.** Heatmap displaying the coumarin and furanocoumarin profiles in the peel of the 61 *Citrus* species investigated. A: heatmap; B: dendrogram of compounds; C: dendrogram of *Citrus* varieties. A red square highlights the presence of a given compound, while a yellow square shows its absence. Coumarins are represented in orange; furanocoumarins of the bergapten cluster, in blue; furanocoumarins of the xanthotoxin cluster, in red; and furanocoumarins of the isopimpinellin cluster, in green.</figcaption></figure>

doi:10.1371/journal.pone.0142757.g004

## Organization of *Citrus* diversity based on coumarin and furanocoumarin contents in peel

Based on the largest compound diversity found in peel (heatmap, Fig 4A), *Citrus* varieties can be classified by HACs into 4 clades (dendrogram, Fig 4C).
<figure><img src="image_9.png" /><figcaption>**Fig 5.** Heatmap displaying the coumarin and furanocoumarin profiles in the pulp of the 61 *Citrus* species investigated. A: heatmap; B: dendrogram of compounds; C: dendrogram of *Citrus* varieties. A red square highlights the presence of a given compound, while a yellow square shows its absence. Coumarins are represented in orange; furanocoumarins of the bergapten cluster, in blue; furanocoumarins of the xanthotoxin cluster, in red; and furanocoumarins of the isopimpinellin cluster, in green.</figcaption></figure>

doi:10.1371/journal.pone.0142757.g005

Clade 1 gathers *Citrus* from the pummelo group, i.e., 7 pummelos and pummelo hybrids (3 grapefruits, 4 sour oranges, and 1 bergamot), which are associated with 2 papedas (Combava, Micrantha) and 2 other hybrids: Yuzu and Nasnaran. These varieties primarily synthesize compounds from the groups I and IV, including all the coumarins biosynthetically originating
from the parent compound umbelliferone (osthol, epoxyaurapten, aurapten, limettin, 5-geranyloxy, and 7-methoxycoumarin; Fig 1) and furanocoumarins from the bergapten cluster: bergapten, bergaptol and its geranylated derivatives (bergamottin, epoxybergamottin and 6',7'-dihydroxybergamottin), and its dimethylallylated derivatives (isoimperatorin, oxypeucedanin and oxypeucedanin hydrate; Fig 4A).

Clade 2 is characterized by the presence of limes (except Palestine), the Khasi hybrid and the Eureka lemon. Compared to the other *Citrus* groups, clade 2 exhibits the greatest coumarin/furanocoumarin diversity as shown in the heatmap representation (Fig 4A), where a large assortment of prenylated coumarins and furanocoumarins derived from bergapten, xanthotoxin and isopimpinellin parent compounds (all 4 groups represented) is displayed.

Clade 3 is composed of the citron group, i.e., the 3 citrons, the lemons (except Eureka) and the Palestine sweet lime. This clade also includes the mandarin × citron hybrid. Fruit peels from clade 3 are characterized by the occurrence of prenylated coumarins associated with prenylated furanocoumarins from the bergapten, xanthotoxin and isopimpinellin clusters (Fig 4A). This clade 3 chemotype is quite comparable to that found in clade 2, although coumarins and furanocoumarins from the previously described groups III and IV are less represented in clade 3.

Clade 4 is composed of all the mandarins (sweet and acidic), all the small mandarin hybrids, all the sweet oranges, Ichang papeda and the Mountain hybrid. Clade 4 includes species that are characterized by limited chemotype diversity, primarily restricted to compounds associated previously with group I (prenylated coumarins and furanocoumarins from the bergapten cluster, Fig 4A).

In addition to qualitative analyses, we used the quantitative results of coumarins and furanocoumarins contents in citrus peel to analyze the implication of phylogenetic origins on coumarins and furanocoumarins contents. NJ analyses and PCAs were performed for this purpose. The four ancestral taxa (mandarins, pummelos, papedas and citrons) could be perfectly separated in the NJ analysis based on the quantitative data of coumarins and furanocoumarins in citrus peel (Fig 6). On this NJ analysis, bootstrap values are quite high (91 for mandarins, 71 for pummelos and papedas and 79 for citrons; Fig 6), which demonstrates that it is a strong representation of the structure of citrus diversity. Based on this strong structure of coumarins and furanocoumarins contents between the four ancestral taxa, we analyzed the coumarins and furanocoumarins contents of the secondary species according to their phylogenomic contents by PCA analysis, using ancestral taxa accessions as active individuals and secondary species accessions as additional individuals.

When using the accessions of the four ancestral taxa, the very high coumarins and furanocoumarins contents of papedas (*C. micrantha* and *C. hystrix*) and some of their hybrids resulted in a very strong impact of these accessions on the PCA axes definition (S1 Fig), even though the quantitative data were transformed as log₁₀(1 + x). As a result, the resolution of the differentiation between the three other basic taxa was low and thus the relative position of secondary species difficult to infer (S1 Fig). Therefore, to analyze the organization of coumarins and furanocoumarins contents of the three other ancestral taxa (pummelos, mandarins and citrons) and their respective hybrids, we performed a PCA analysis with only the corresponding accessions (Fig 7). The structure of the diversity of citrus varieties is well represented on the first two axes of this PCA with a value above 67% (Fig 7). Pummelos, mandarins and citrons are well separated on the PCA and the secondary species are well distributed among these ancestral taxa (Fig 7A). For example, varieties with a *C. maxima* × *C. reticulata* admixture structure are placed on an axis between these 2 taxa. In the same way, the citrus varieties with a *C. medica* × *C. reticulata* admixture structure genome are distributed between these 2 ancestral taxa. Pummelos, grapefruits and sour oranges are characterized by their peel content in
<figure><img src="image_10.png" /><figcaption>Fig 6. NJ analysis of the 23 citrus varieties belonging to the 4 ancestral taxa based on the coumarins and furanocoumarins contents in peel expressed as mg.kg⁻¹ fresh weight. Numbers in black represent the bootstrap probability values. The colors correspond to the phylogenetic constitution of the varieties and are indicated in Figs 2 and 3.</figcaption></figure>

doi:10.1371/journal.pone.0142757.g006

dimethylallylated coumarins (aurapten, epoxyaurapten and osthol) and dimethylallylated furanocoumarins from the bergapten cluster essentially (bergamottin, epoxybergamottin and 6', 7'-dihydroxybergamottin) (Fig 7B). In contrast, citrons, lemons and limes are characterized by their peel content in limettin and in dimethylallylated furanocoumarins from the xanthotoxin cluster (imperatorin, heraclenin and heraclenol) and from the isopimpinellin cluster (phellopterin, byakangelicol and byakangelicin) (Fig 7B). Mandarins, sweet oranges and mandarin hybrids are poor producers of coumarins and furanocoumarins, even though they produce isopimpinellin and bergapten (Fig 7B).

The PCA with the four basic taxa (S1 Fig) reveals the position of papeda hybrids that are placed between papedas and citrons. The coumarins and furanocoumarins contents of Alemow
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**Fig 7.** PCA representation of the citrus varieties belonging to the 3 ancestral taxa pummelos, mandarins and citrons and their respective hybrids based on the coumarins and furanocoumarins contents in peel expressed as mg.kg⁻¹ fresh weight. A: Distribution of the observations (citrus varieties); B: Distribution of the variables (coumarins and furanocoumarins). Ancestral taxa were used as active individuals (circles) to build the PCA, while secondary species are supplementary individuals (squares). Quantitative data were transformed as log₁₀(1 + x). The colors correspond to the phylogenetic constitution of the varieties and are indicated in Figs 2 and 3.</figcaption></figure>

doi:10.1371/journal.pone.0142757.g007
and Mexican limes place them closer to citrons than papedas, while Nasnaran is closer to papedas than to mandarins (S1 Fig).

## Organization of *Citrus* diversity based on coumarin and furanocoumarin contents in pulp

The identification of coumarins and furanocoumarins in the pulp allowed 3 clades (Fig 5C) to be distinguished instead of the 4 described previously based on the peel content (Fig 4C). For clarity between Figs 5C and 4C, these 3 new clades are nominated 1', 2', and 3'.

Clade 1' is composed of the 7 pummelos, the 4 sour oranges, the 3 grapefruits, the Combava, the Bergamot, the Nasnaran mandarin and the Eureka lemon. Clade 1' is composed of the same varieties found previously in clade 1 (Fig 4C), with the exception of Yuzu and Micrantha, which are found only in clade 1, and the Eureka lemon, which is restricted to clade 1'. The pulp of these varieties primarily contains prenylated coumarins and prenylated derivatives from the bergapten cluster classified in groups I' and III'.

Clade 2' includes all the sweet and acidic mandarins, all the small mandarin hybrids, all the sweet oranges, the lemons (except Eureka), 2 limes (Brazil and Palestine), the 2 citrons for which the pulp compartment could be investigated (NB: due to the absence of pulp in Buddha's Hand citron fruits, no investigation could be conducted in the pulp of this third citron variety), Ichang, Mountain, Yuzu and the mandarin × citron hybrid. Clade 2' is composed of varieties classified previously in clades 1 (Yuzu), 2 (Brazil sweet lime), 3 (citrons, lemons and limes) and 4 (mandarins, small mandarin hybrids, sweet oranges, Mountain and Ichang) based on the peel chemotype. Varieties from clade 2' show the narrowest diversity in coumarins and furanocoumarins in the fruit pulp, primarily restricted to the compounds classified in group III'.

Clade 3' contains the limes (except Brazil and Palestine), the Micrantha papeda and Khasi. These varieties are associated with clades 1 (Micrantha) and 2 (limes and Khasi) established previously from the peel chemotype. The pulp observed from clade 3' varieties shows the most diversified chemotypes, with coumarins and furanocoumarins belonging to groups I', II', and III'. In particular, dimethylallylated compounds from the xanthotoxin cluster (i.e., imperatorin, heraclenin, and heraclenol), as well as the parent compound xanthotoxol and its methoxylated derivative xanthotoxin, are exclusively restricted to some of the clade 3' fruit pulp.

As we did before, the quantitative data of coumarins and furanocoumarins contents in citrus pulp were used to build a NJ analysis between the four ancestral taxa. This NJ analysis displays a weak representation of the structure of citrus diversity, with low bootstrap values for the main branching and incomplete differentiation between the four basic taxa (S2 Fig). Therefore, citrus pulp coumarins and furanocoumarins contents are not strongly related with phylogenetic differentiation. Thus we do not perform an analysis to study the relationship of coumarins and furanocoumarins contents on secondary species and their phylogenetic origin.

## Discussion

### Peel presents the richest coumarin/furanocoumarin content

The results obtained from the peel and pulp of the 61 studied varieties clearly demonstrate that the envelope of the fruit is the richest part in terms of both molecular diversity (Fig 4 vs. Fig 5) and compound concentration (Figs 2 and 3, Tables A to L in S1 File). Peel prominence regarding coumarin and furanocoumarin contents must be related to the presence of lysogenic cavities in the flavedo (i.e., pigmented pericarp of the peel), where large quantities of essential oils and polyphenolics are stored. The occurrence of specialized cells involved in the biosynthesis of coumarins and furanocoumarins has been hypothesized in young grapefruits (*C. paradisi*
'Duncan') [43]: epithelial cells surrounding the secretory cavities show a high level of tran-
scripts potentially committed to the biosynthetic routes, which most likely reflects their special-
ization for coumarin and furanocoumarin biosynthesis. Future studies aimed at classifying
*Citrus* varieties based on coumarin/furanocoumarin chemotypes should be limited to the peel
compartment only.

Coumarin and furanocoumarin diversity found in *Citrus* reflects the
organization of the biosynthetic pathways

Qualitative studies are useful to understand the organization of the biosynthetic pathways, as they reflect the ability of varieties to produce different metabolites. HACs based on the presence/absence of coumarins/furanocoumarins in peel allowed 4 categories of compounds to be distinguished (Fig 4B, groups I to IV). Each group is logically structured with coumarins and furanocoumarins engaged in the same cascade of enzymatic reactions (i.e., groups of substrates and products either leading to the same reaction end product or originating from the same precursor compound).

Group I is composed of one dimethylallylated coumarin (aurapten), one geranylated coumarin (5-geranyloxy, 7-methoxycoumarin) and one dimethoxy coumarin (limettin). All these products originate from the biosynthetic transformation of the parent compound umbelliferone. In accordance with the rule that substrates and products from the same biosynthetic route are regrouped, group I is composed of the 3 prenylated derivatives that form after the dimethylallylation of bergaptol (i.e., isoimperatorin, oxypeucedanin, and oxypeucedanin hydrate) and one prenylated derivative (bergamottin) that forms after the geranylation of the same bergaptol.

Notably, group I has the C5, C8 dimethoxylated furanocoumarin isopimpinellin and the C5 mono-methoxylated bergapten, both of which are ubiquitously found in *Citrus*, whereas xanthotoxin, the C8-isomer of bergapten that occurs in group III, is primarily restricted to *Citrus* from clade 2 (limes and papedas). At present, the route leading from psoralen to isopimpinellin is not established, and whether isopimpinellin is formed after C5 methoxylation of xanthotoxin, C8 methoxylation of bergapten, or after both is unknown [16,35]. Based on our data regarding the strong co-occurrence of isopimpinellin and bergapten, isopimpinellin is likely derived from bergapten in *Citrus*.

Interestingly, the 3 compounds that form after O-geranylation of bergaptol (bergamottin → epoxybergamottin → 6', 7', dihydroxybergamottin) are not clustered in the same group: bergamottin (i.e., 5 geranyloxy-psoralen), the first compound, is located in group I, and the latter 2 epoxy derivatives are in group IV. Symmetrically, group IV contains epoxyaurapten, another epoxy derivative for which the parent geranyloxy compound (aurapten) is located in group I. Because compounds from group I are represented in almost all *Citrus* accessions and because compounds from group IV are principally limited to pummelos and their hybrids, it is tempting to hypothesize that the enzymes responsible for the synthesis of O-geranylated compounds (i.e., aromatic O-prenyltransferases) might be ubiquitously spread in *Citrus*, whereas the enzymes in charge of the synthesis of the two epoxydes (epoxybergamottin and epoxyaurapten) and of 6',7',dihydroxybergamottin might be more scarcely distributed and restricted to the pummelo group (clade 1). Group IV also includes osthol, the unique C-prenylated compound investigated in this study. Notably, the distribution of this C-prenylated structure is much more limited than other O-prenylated coumarins (aurapten, 5-geranyloxy, and 7-methoxycoumarin), although it is difficult to draw general conclusions regarding the distribution C vs. O prenyltransferases due to the lack of C-prenylated compounds analyzed in this study.

Group II is exclusively composed of *O*-prenylated compounds characterized by a dimethy-
lallyloxy chain bound at C8 of the furanocoumarin core structure. These compounds include
imperatorin, heraclenin, and heraclenol from the xanthotoxin cluster and cnidicin, phellopterin, byakangelicol, and byakangelicin from the isopimpinellin cluster. They are principally found in *Citrus* clade 2 (limes and papedas) and clade 3 (citrons, lemons, and sour limes). Whether a single aromatic O-prenyltransferase might operate on multiple substrates such as xanthotoxol and 8-hydroxybergapten or whether different enzymes are engaged in the synthesis of imperatorin and phellopterin remain unknown. This point is presently difficult to address because the number of aromatic prenyltransferases described from plants thus far is limited to only seven. Notably, five of these seven aromatic prenyltransferases appear to be extremely specific regarding the prenyl donors that they accept (DMAPP, GPP) [44-46,38], whereas 2 prenyltransferases described recently from hop plant show affinity for both DMAPP and GPP [47]. Four of these enzymes also demonstrate a high specificity for prenyl acceptors, as they only tolerate a single aromatic entity as a substrate [44-46,38], including a C-prenyltransferase described recently from lemon that is involved in the synthesis of 8-geranylumbelliferone [46]. A 5th enzyme isolated from *Sophora flavescens* accepts more promiscuous aromatic substrates, as it is able to prenylate 2 structurally close isoflavones: genistein and its methylated derivative biochanin A [45]. Consistent with this *Sophora flavescens* enzyme, the two hop prenyltransferases described recently show affinity for different substrates such as acylphloroglucinols and naringenin chalcone [47].

Group III primarily contains intermediate compounds that are subsequently transformed into prenylated/methoxylated/epoxydated derivatives such as psoralen, xanthotoxol, bergaptol, xanthotoxin, and umbelliferone. These intermediate compounds are the most scarcely found among citrus fruit peel (Fig 4A) because they are rapidly converted into more chemically diverse end products. Group III also includes cnidilin and 8-geranyloxypsoralen, two prenylated compounds rarely found in *Citrus*.

## Coumarin and furanocoumarin chemical diversity in relation to the recent advances in *Citrus* phylogeny

Based on molecular markers and sequencing, 4 ancestral taxa, *C. maxima*, *C. medica*, *C. reticulata* and *C. micrantha*, have been identified in *Citrus* [8,10,12,14]. These basic taxa also could be well distinguished based on the coumarins and furanocoumarins contents in their peel (Fig 6). Ancestral taxa present no evidence of interspecific introgression, except in sweet mandarin, which presents approximately 10% introgression of *C. maxima* [13,14]. However, the acidic mandarin 'Cleopatra', 'Sunki' and 'Shekwasha' present no evidence of introgression [14]. The coumarin and furanocoumarin chemotypes observed in fruit peel and pulp were compared within the 4 ancestral taxa and were relatively homogeneous inside the taxa, except in papedas. *C. maxima* presents relatively high coumarin and furanocoumarin contents characterized by the presence of all coumarins and bergapten-related compounds (bergapten, bergaptol and its geranylated derivatives, bergamottin, epoxybergamottin, and 6',7'-dihydroxybergamottin) (Figs 2 and 3, Table A in S1 File). *C. reticulata* presents the lowest coumarin and furanocoumarin contents (Figs 2 and 3, Table D in S1 File). This taxon is characterized by the total absence of xanthotoxin in peel and pulp. No difference was observed between sweet and acidic mandarin. *C. medica* is characterized by the occurrence of prenylated coumarins associated with prenylated furanocoumarins from bergapten, xanthotoxin and isopimpinellin clusters (Figs 2 and 3, Table B in S1 File). The Papeda taxon is not cultivated, except Combava, and remains the least known of all ancestral taxa. Only few genetic resources for this group are available in the different germplasm banks in the world, and extreme profiles can be observed within this group. Micrantha and Combava present the highest furanocoumarin content (from bergapten, xanthotoxin and isopimpinellin) (Figs 2 and 3, Table C in S1 File). However, we
must keep in mind that Micrantha was harvested outside Corsica, where all other samples were
collected, and that different pedo-climatic conditions might influence coumarin and furano-
coumarin synthesis [48]. In contrast, Ichang papeda, also considered a true papeda, has very
low coumarin and furanocoumarin contents (Figs 2 and 3, Table C in S1 File). In different
studies [8,15], C. *ichangensis* (Ichang papeda) is grouped with the wild mandarins C. *tachibana*
and C. *mangshanensis*. Unfortunately, these wild mandarins were unavailable to our study
because these genotypes are only present in Asia. However, these species are considered geneti-
cally distinct from other mandarins although they clustered with them. Our results are concor-
dant with these studies, as C. *Ichangensis* is localized with all the mandarin accessions.

All other species result from hybridization between these four ancestral taxa. The coumarins and furanocoumarins found in these secondary species display contrasting molecular diversity associated with variable quantitative traits. Many of these species, such as C. sinensis, C. aurantium, C. paradisi, and tangors, originate from C. reticulata x C. maxima hybridization, with variable contributions of each taxon. For the tangors, the contribution of C. reticulata is higher (from 63% to 89%) [49]. C. sinensis presents a contribution of 60% of C. reticulata and 40% of C. maxima, C. aurantium 50% and 50%; and C. paradisi approximately 40% and 60% respectively [49]. Mandarin hybrids (tangors) and sweet orange, closely related to mandarins, are logically present in the coumarin/furanocoumarin mandarin clade (clade IV, Fig 4C). In the same way, these secondary species are closer to the mandarin cluster on the PCA based on the citrus peel content in coumarins and furanocoumarins (Fig 7A). This taxonomical group is of special interest because it contains the lowest coumarin and furanocoumarin producers amongst all the studied Citrus species. In contrast, C. aurantium and C. paradisi have similar coumarin and furanocoumarin profiles compared with C. maxima (Fig 4C, clade I and Fig 7A). All these results show a strong concordance between phytochemical and molecular data.

*C. limon* (Swingle and Reece taxonomy) is more complex and originates from 2 or 3 taxa hybridization. *C. limon* can result from a direct hybridization between *C. medica* and *C. reticulata* ('Volkamer lemon', 'Rough lemon' and 'Rangpur' lime) or more complex hybridization between 3 ancestral taxa, *C. medica*, *C. reticulata*, and *C. maxima* ('Meyer', 'Eureka' lemons and Marrakech lime) [49]. Most of these species are present in the same clade 3 as citron (*C. medica*, Fig 4C) or are placed on an axis between the citron and mandarin clusters on the PCA based on the quantitative results in citrus peel (Fig 7A).

Similarly, *C. aurantifolia* originates from different taxa hybridization (2 to 4 taxa) composed systematically of citron and papeda taxons. 'Alemow', 'Excelsa', Mexican' and Giant key limes display extremely similar phylogenomic origins, with close to 50% contributions from *C. medica* and *C. micrantha* (14). *C. aurantifolia* 'Palestine' and 'Brazil' sweet lime have phylogenetic structures similar to *C. limon* ('Meyer', 'Eureka' lemons, and 'Marrakech' lime). *C. aurantifolia* 'Palestine' and 'Brazil' sweet lime contains close to 50% *C. medica* and displayed 2 other ancestral taxa, *C. reticulata* and *C. maxima*, with a greater contribution from *C. reticulata* (32-34%) than *C. maxima* (15-20% in 'Palestinian' sweet lime) [14,49]. Coppenrhad lime is a triploid originating from a (*C. micrantha* x *C. medica*) x *C. medica* hybridization with a diploid gamete of the *C. micrantha* x *C. medica* parent [49]. Tahiti and Bearss lime, which are also triploid, have contributions from the 4 ancestral taxa (1/6 *C. maxima*, 1/3 *C. medica*, 1/3 *C. micrantha* and 1/6 *C. reticulata*) [49]. All of these secondary species containing Papeda, except Palestine sweet lime, are grouped in the same clade 2 Fig 4C), indicating that our phytochemical data are consistent with the recent molecular data cited previously. The mandarin × citron hybrid is a simple hybrid of *C. medica* x *C. reticulata* [49] as Rough lemon, Rangpur lime and Volkamer lemon. This hybrid has similar profiles as the coumarin and furanocoumarin citron profiles and is grouped in clade 3 (Fig 4C).
Bergamot is considered to have emerged from a lemon X sour orange crossing and contains 40% of pummelo [49], which is consistent with its position in the pummelo clade. In the same way, bergamot is placed between pummelos and mandarins, close to the sour oranges, on the PCA based on the quantitative results of coumarins and furanocoumarins in peel (Fig 7A). With reference to our survey, this finding could be due to the persistence of some chemical compounds that might have been inherited from pummelo such as osthol, epoxyaurapten and bergamottin (Figs 4 and 7).

The origin of the Mountain hybrid remains undetermined; however, this hybrid has the same maternal origin as papeda [10]. Yuzu might be derived from the hybridization between Ichang papeda and a mandarin [5]. These 2 hybrids are present in the mandarin clade, which also contains the Ichang papeda displaying extremely low coumarin and furanocoumarin contents. Thus, based on our phytochemical data, these hybrids are most likely derived from Ichang papeda and/or mandarin.

The Khasi hybrid displays coumarin and furanocoumarin profiles similar to Micrantha papeda, allowing us to draw the hypothesis of affiliation with this cultivar.

Nasnaran mandarin is a direct mandarin-papeda hybrid [14] with a papeda mitotype [10]. This hybrid presents a strong papeda fragrance, although it looks similar to a mandarin. Here, its profile is more closely related to papedas, particularly Combava, which might suggest a parental link between them.

In brief, all secondary species resulting from a cross directly involving the ancestral taxa *C. maxima* (sour oranges and grapefruits), *C. reticulata* (sweet oranges and some tangors), *C. medica* (some lemons) or *C. micrantha* (some limes and Nasnaran mandarin) present coumarin and furanocoumarin profiles similar to the parent ancestral taxa, with a prevalence of the chemotypes found in pummelos, citrons and papedas over mandarins and, identically, a prevalence of the chemotypes found in papedas over citrons. These findings constitute valuable information for future breeding programs that will allow the coumarin and furanocoumarin contents in *Citrus* species to be reduced.

## Breeding programs for selection against coumarins and furanocoumarins in *Citrus*: Global and targeted approaches

We observed previously that the lowest coumarin and furanocoumarin contents are found in sweet and acidic mandarins and most of their hybrids (including sweet oranges and small mandarin hybrids), Ichang papeda and the 2 hybrids Yuzu and Mountain. In contrast, pummelos, citrons and papedas, which contain the highest quantities of coumarins and furanocoumarins, generate hybrids that are themselves rich in these compounds. This finding led us to hypothesize that breeding programs aimed at decreasing the toxic coumarin and furanocoumarin contents in *Citrus* should preferentially involve mandarins or Ichang papeda. *C. sinensis* is a good example of the possibility of obtaining a low coumarin and furanocoumarin hybrid from *C. reticulata*, although it has a relatively high genetic background from pummelo (40%) [14]. This hypothesis is in agreement with a study evaluating the inhibition of the CYP3A4 by immature *Citrus* fruit extracts [50]. In this work, *C. paradisi* was shown to be the most effective inhibitor of CYP3A4, while *C. unshiu* displayed relatively weak inhibition of the cytochrome most likely due to low furanocoumarin content. Other *Citrus* fruits were tested; their capacities to inhibit CYP3A4 are ordered as follows: *C. paradisi* > *C. limon* > *C. aurantium* > *C. grandis* (also called *C. maxima*) > *C. reticulata* > *C. sinensis* > *C. unshiu* [50]. This order is globally consistent with our quantitative results, considering that the coumarin and furanocoumarin contents vary greatly from one species to another and depend on environmental conditions.
Pummelos are high producers of geranylated compounds from the bergapten cluster: bergamottin, epoxybergamottin and 6',7'-dihydroxybergamottin. These chemical compounds are phototoxic and are strong inhibitors of CYP3A4 [20,51,52]. Sour oranges and grapefruits are direct descendants of pummelos that also produce these 3 compounds, and their other parents (sweet mandarin and sweet orange, respectively) generally do not produce them. Therefore, we can hypothesize that crosses involving pummelos and mandarins or sweet oranges usually include the toxic furanocoumarins bergamottin, epoxybergamottin and 6',7'-dihydroxybergamottin (Fig 8). This is in accordance with a recent statement issued by Xu and coworkers, who confirmed that pummelos as parents have a definite contribution to the furanocoumarins content (bergapten, bergamottin, epoxybergamottin and 6',7'-dihydroxybergamottin) of grapefruits [42]. However, sweet oranges and small mandarin hybrids that have a greater contribution of *C. reticulata* in their genotype have lost the capacity to produce these 3 compounds (Fig 8). This finding suggests that the genes encoding the enzyme machinery involved in the biosynthesis of bergamottin, epoxybergamottin and 6',7'-dihydroxybergamottin in pummelos might be dominant and transmitted to their direct hybrid progenies. This hypothesis of dominant genes is consistent with molecular data regarding coumarin/furanocoumarin biosynthesis available to date, where the capacity to synthesize such compounds is given by very specific enzymes encoded by specialized genes such as dioxygenases, prenyltransferases, and cytochrome P450s [35–38,46].

<figure><img src="image_13.png" /><figcaption>**Fig 8.** Schematic representation of coumarin/furanocoumarin inheritance from the 4 ancestral taxa (pummelo, mandarin, citron, and papeda) in the cultivated *Citrus* species. Parental relations between species are illustrated by arrows. Thick arrows represent high chemotype similarities between hybrids and their ascendants. Same colors were used as in Figs 4 and 5 for representing *Citrus* species. Coumarins and furanocoumarins are indicated in the same color as the ancestral taxon that transmitted them to their descendants. The symbol "+" indicates the probable transmission of compounds in the secondary species, while the symbol "-" indicates their disappearance.</figcaption></figure>

doi:10.1371/journal.pone.0142757.g008
Crosses involving citrons are interesting, with the objective to primarily remove the toxic epoxybergamottin and 6',7'-dihydroxybergamottin introduced by pummelos (see lemon progenies: Eureka, Meyer and Marrakech) or by *C. micrantha* (see lime progenies) (Fig 8). However, bergamottin and the phototoxic xanthotoxin [53] from Micrantha are maintained in some limes, and citrons bring other phototoxic compounds to their descendants such as heraclenin and heraclenol [51,54], which are absent in pummelos and Micrantha (Fig 8). Acidic mandarins might be a solution to prevent the transmission of heraclenin and heraclenol from citrons because among the 4 *C. reticulata* × *C. medica* lemon hybrids investigated (Rough, Rangpur, Yellow Rangpur and Volkamer), only the Volkamer lemon produces these 2 compounds (Fig 8). Alternatively, citrons can also bring the phototoxic limettin and oxypeucedanin hydrate [51] to lemons (see Eureka, Meyer and Marrakech lemons progenies), 2 chemical compounds absent in sour oranges (Fig 8). This problem cannot be solved by the use of mandarins as the parental line (see Rough, Rangpur and Volkamer lemons progenies) because they all contain limettin and sometimes oxypeucedanin hydrate (Fig 8).

In conclusion, this study highlights the different chemotypes found in *Citrus* species regarding coumarin and furanocoumarin diversity and content. Three of 4 ancestral taxa (pummelos, citrons and papedas) synthesize high amounts of these compounds, whereas mandarins, from the last ancestral taxon, appear practically devoid of them. As for hybrid species, their corresponding chemotypes appear inherited from the 4 ancestral taxa, with strong hypotheses of a prevalence of pummelos, citrons and papedas over mandarins and of papedas over citrons. Compared to other species, mandarins are of invaluable interest for *Citrus* breeding programs aimed at removing toxic furanocoumarins in *Citrus*. Particularly, sweet oranges and small mandarins hybrids, which have a higher contribution of mandarin in their genotype, do not produce the toxic compounds bergamottin, epoxybergamottin and 6',7'-dihydroxybergamottin. An alternative strategy for the future would lie in the identification of an upstream gene involved in the first steps leading to the furanocoumarin pathway. Such gene characterization should allow gene editing strategies to be performed and would constitute an efficient solution for removing these toxic compounds from *Citrus* species.

**Supporting Information**

S1 Fig. PCA representation of the 61 citrus varieties based on the coumarins and furanocoumarins contents in peel expressed as mg.kg⁻¹ fresh weight. Ancestral taxa were used as active individuals (circles) to build the PCA, while secondary species are supplementary individuals (squares). Quantitative data were transformed as log₁₀(1 + x). The colors correspond to the phylogenetic constitution of the varieties and are indicated in Figs 2 and 3.
(TIF)

S2 Fig. NJ analysis of the 23 citrus varieties belonging to the 4 ancestral taxa based on the coumarins and furanocoumarins contents in pulp expressed as mg.kg⁻¹ fresh weight. Numbers in black represent the bootstrap probability values. The colors correspond to the phylogenetic constitution of the varieties and are indicated in Figs 2 and 3.
(TIF)

S1 File. Concentrations (in mg kg⁻¹ fresh weight ± standard deviation) of the coumarins and the furanocoumarins in the peel and pulp extracts of the 61 examined citrus species.
(PDF)

S1 Table. List of the 61 citrus species investigated in this study and their phylogenetic con-
stitutions.
(PDF)
S2 Table. Limits of detection and quantitation of the coumarins and furanocoumarins in the citrus pulp and specificity (equation and coefficient of determination, r²) of the UPLC-MS method.

(PDF)

## Author Contributions

Conceived and designed the experiments: YF FB. Performed the experiments: AD AO GC. Analyzed the data: AD AO AH PO YF FB. Contributed reagents/materials/analysis tools: AD AO GC PO YF. Wrote the paper: AD PO YF FB.

## References

1. Praloran J-C. Citrus. Paris: G.P. Maisonneuve et Larose; 1971.

2. Gmitter FG, Chen C, Rao MN, Soneji JR. Citrus Fruits. In: Kole C, editor. Fruits and Nuts. Springer Berlin Heidelberg; 2007. pp. 265–279. Available: http://link.springer.com.bases-doc.univ-lorraine.fr/chapter/10.1007/978-3-540-34533-6_14.

3. Nelson DR. Cytochrome P450 homepage [Internet]. 2015 [cited 27 Apr 2015]. Available: http://dnelson.uthsc.edu/CytochromeP450.html.

4. Groppo M, Pirani JR, Salatino MLF, Blanco SR, Kallunki JA. Phylogeny of Rutaceae based on twonon-coding regions from cpDNA. Am J Bot. 2008; 95: 985–1005. doi: 10.3732/ajb.2007313 PMID: 21632420

5. Swingle W, Reece P. The botany of Citrus and its wild relatives. In: Reuther W, Webber H, Batchelor L, editors. The citrus industry. Berkeley, CA: University of California; 1967. pp. 190–430.

6. Tanaka T. Fundamental discussion of Citrus classification. Stud Citrol. 1977; 14: 1–6.

7. Federici CT, Fang DQ, Scora RW, Roosevelt ML. Phylogenetic relationships within the genus Citrus (Rutaceae) and related genera as revealed by RFLP and RAPD analysis. Theor Appl Genet. 1998; 96: 812–822. doi: 10.1007/s001220050807

8. Nicolosi E, Deng ZN, Gentile A, Malfa SL, Continella G, Tribulato E. Citrus phylogeny and genetic origin of important species as investigated by molecular markers. Theor Appl Genet. 2000; 100: 1155–1166. doi: 10.1007/s001220051419

9. Barkley NA, Roosevelt ML, Krueger RR, Federici CT. Assessing genetic diversity and population structure in a citrus germplasm collection utilizing simple sequence repeat markers (SSRs). Theor Appl Genet. 2006; 112: 1519–1531. PMID: 16699791

10. Froelicher Y, Mouhaya W, Bassene J-B, Costantino G, Kamiri M, Luro F, et al. New universal mitochon-drial PCR markers reveal new information on maternal citrus phylogeny. Tree Genet Genomes. 2011; 7: 49–61. doi: 10.1007/s11295-010-0314-x

11. Ollitrault P, Terol J, Garcia-Lor A, Bérard A, Chauveau A, Froelicher Y, et al. SNP mining in C. clemen-tina BAC end sequences; transferability in the Citrus genus (Rutaceae), phylogenetic inferences and perspectives for genetic mapping. BMC Genomics. 2012; 13: 13. doi: 10.1186/1471-2164-13-13 PMID: 22233093

12. Garcia-Lor A, Curk F, Snoussi-Trifa H, Morillon R, Ancillo G, Luro F, et al. A nuclear phylogenetic analy-sis: SNPs, indels and SSRs deliver new insights into the relationships in the “true citrus fruit trees” group (Citrinae, Rutaceae) and the origin of cultivated species. Ann Bot. 2013; 111: 1–19. doi: 10.1093/aob/mcs227 PMID: 23104641

13. Wu GA, Prochnik S, Jenkins J, Salse J, Hellsten U, Murat F, et al. Sequencing of diverse mandarin, pummelo and orange genomes reveals complex history of admixture during citrus domestication. Nat Biotechnol. 2014; 32: 656–662. doi: 10.1038/nbt.2906 PMID: 24908277

14. Curk F, Ancillo G, Ollitrault F, Perrier X, Jacquemoud-Collet J-P, Garcia-Lor A, et al. Nuclear species-diagnostic SNP markers mined from 454 amplicon sequencing reveal admixture genomic structure of modern citrus varieties. PloS One. 2015; 10: e0125628. doi: 10.1371/journal.pone.0125628 PMID: 25973611

15. Carbonell-Caballero J, Alonso R, Ibañez V, Terol J, Talon M, Dopazo J. A Phylogenetic Analysis of 34 Chloroplast Genomes Elucidates the Relationships between Wild and Domestic Species within the Genus Citrus. Mol Biol Evol. 2015; doi: 10.1093/molbev/msv082

16. Bourgaud F, Hehn A, Larbat R, Doerper S, Gontier E, Kellner S, et al. Biosynthesis of coumarins in plants: a major pathway still to be unravelled for cytochrome P450 enzymes. Phytochem Rev. 2006; 5: 293–308.
17. Berenbaum MR, Nitao JK, Zangerl AR. Adaptive significance of furanocoumarin diversity in Pastinaca sativa (Apiaceae). J Chem Ecol. 1991; 17: 207–15. doi: 10.1007/BF00994434 PMID: 24258446

18. Robl M, Robl R, Marinoni LP, Abagge KT, Carvalho VO. Assemble the puzzle: bizarre-looking lesions. Arch Child. 2013; 98: 915.

19. Calka O, Akdeniz N, Metin A, Behcet L. Phototoxic dermatitis due to Chenopodium album in a mother and son. Contact Dermatitis. 2005; 53: 58–60. PMID: 15982237

20. Messer A, Raquet N, Lohr C, Schrenk D. Major furocoumarins in grapefruit juice II: phototoxicity, photo-genotoxicity, and inhibitory potency vs. cytochrome P450 3A4 activity. Food Chem Toxicol. 2012; 50: 756–760. doi: 10.1016/j.fct.2011.11.023 PMID: 22155270

21. Ohta T, Maruyama T, Nagahashi M, Miyamoto Y, Hosoi S, Kiuchi F, et al. Paradisin C: a new CYP3A4 inhibitor from grapefruit juice. Tetrahedron. 2002; 58: 6631–6635.

22. Gravot A, Larbat R, Hehn A, Lièvre K, Gontier E, Goergen J-L, et al. Cinnamic acid 4-hydroxylase mechanism-based inactivation by psoralen derivatives: cloning and characterization of a C4H from a psoralen producing plant _Ruta graveolens_ exhibiting low sensitivity to psoralen inactivation. Arch Biochem Biophys. 2004; 422: 71–80. PMID: 14725859

23. Lin H, Kenaan C, Hollenberg PF. Identification of the residue in Human CYP3A4 that is covalently modified by bergamottin and the reactive intermediate that contributes to the grapefruit juice effect. Drug Metab Dispos. 2012; 40: 998–1006. doi: 10.1124/dmd.112.044560 PMID: 22344702

24. Wilkinson GR. Drug metabolism and variability among patients in drug response. N Engl J Med. 2005; 352: 2211–2221. doi: 10.1056/NEJMra032424 PMID: 15917386

25. Bailey DG, Dresser G, Arnold JM. Grapefruit-medication interactions: forbidden fruit or avoidable consequences? Can Med Assoc J. 2013; 185: 309–16.

26. Bailey DG, Spence JD, Munoz C, Arnold JM. Interaction of citrus juices with felodipine and nifedipine. Lancet. 1991; 337: 268–9. PMID: 1671113

27. UNCTAD. Citrus consumption [Internet]. 2011. Available: http://www.unctad.info/en/Infocomm/.

28. Liu Y, Heying E, Tanumihardjo SA. History, global distribution, and nutritional importance of Citrus fruits. Compr Rev Food Sci Food Saf. 2012; 11: 530–545.

29. Hwang S-L, Shih P-H, Yen G-C. Neuroprotective effects of citrus flavonoids. J Agric Food Chem. 2012; 60: 877–885. doi: 10.1021/jf204452y PMID: 22224368

30. Buscemi S, Rosafio G, Arcoleo G, Mattina A, Canino B, Montana M, et al. Effects of red orange juice intake on endothelial function and inflammatory markers in adult subjects with increased cardiovascular risk. Am J Clin Nutr. 2012; 95: 1089–95. doi: 10.3945/ajcn.111.031088 PMID: 22492368

31. Paine MF, Criss AB, Watkins PB. Two major grapefruit juice components differ in time to onset of intestinal CYP3A4 inhibition. J Pharmacol Exp Ther. 2005; 312: 1151–1160. PMID: 15485894

32. Uesawa Y, Mohri K. The use of heat treatment to eliminate drug interactions due to grapefruit juice. Biol Pharm Bull. 2006; 29: 2274–2278. PMID: 17077528

33. Uesawa Y, Mohri K. UV-irradiated grapefruit juice loses pharmacokinetic interaction with nifedipine in rats. Biol Pharm Bull. 2006; 29: 1286–1289. PMID: 16755036

34. Myung K, Narciso JA, Manthey JA. Removal of furanocoumarins in grapefruit juice by edible fungi. J Agric Food Chem. 2008; 56: 12064–12068. doi: 10.1021/jf802713g PMID: 19012403

35. Larbat R, Kellner S, Specker S, Hehn A, Gontier E, Hans J, et al. Molecular cloning and functional characterization of psoralen synthase, the first committed monooxygenase of furanocoumarin biosynthesis. J Biol Chem. 2007; 282: 542–554. PMID: 17068340

36. Larbat R, Hehn A, Hans J, Schneider S, Jugde H, Schneider B, et al. Isolation and functional characterization of CYP71AJ4 encoding for the first P450 monooxygenase of angular furanocoumarin biosynthesis. J Biol Chem. 2009; 284: 4776–4785. doi: 10.1074/jbc.M807351200 PMID: 19098286

37. Vialart G, Hehn A, Olry A, Ito K, Krieger C, Larbat R, et al. A 2-oxoglutarate-dependent dioxygenase from Ruta graveolens L. exhibits p-coumaroyl CoA 2'-hydroxylase activity (C2'H): a missing step in the synthesis of umbelliferone in plants. Plant J. 2012; 70: 460–470. doi: 10.1111/j.1365-313X.2011.04879.x PMID: 22168819

38. Karamat F, Olry A, Munakata R, Koeduka T, Sugiyama A, Paris C, et al. A coumarin-specific prenyl-transferase catalyzes the crucial biosynthetic reaction for furanocoumarin formation in parsley. Plant J Cell Mol Biol. 2014; 77: 627–638. doi: 10.1111/tpj.12409

39. Dugrand A, Olry A, Duval T, Hehn A, Froelicher Y, Bourgaud F. Coumarin and furanocoumarin quantitation in citrus peel via ultraperformance liquid chromatography coupled with mass spectrometry (UPLC-MS). J Agric Food Chem. 2013; 61: 10677–10684. doi: 10.1021/jf402763t PMID: 24117278
40. Vandermolen KM, Cech NB, Paine MF, Oberlies NH. Rapid quantitation of furanocoumarins and flavonoids in grapefruit juice using ultra-performance liquid chromatography. Phytochem Anal. 2013; 24: 654-660. doi: 10.1002/pca.2449 PMID: 23780830

41. Durand-Hulak M, Dugrand A, Duval T, Bidel LPR, Jay-Allemand C, Froelicher Y, et al. Mapping the genetic and tissular diversity of 64 phenolic compounds in Citrus species using a UPLC-MS approach. Ann Bot. 2015; 115: 861-877. doi: 10.1093/aob/mcv012 PMID: 25757470

42. Xu J, Ma L, Jiang D, Zhu S, Yan F, Xie Y, et al. Content evaluation of 4 furanocoumarin monomers in various citrus germplasms. Food Chem. 2015; 187: 75-81. doi: 10.1016/j.foodchem.2015.04.007 PMID: 25977000

43. Voo SS, Grimes HD, Lange BM. Assessing the biosynthetic capabilities of secretory glands in Citrus peel. Plant Physiol. 2012; 159: 81-94. doi: 10.1104/pp.112.194233 PMID: 22452856

44. Sasaki K, Mito K, Ohara K, Yamamoto H, Yazaki K. Cloning and characterization of naringenin 8-phenyltransferase, a flavonoid-specific prenyltransferase of Sophora flavescens. Plant Physiol. 2008; 146: 1075-1084. doi: 10.1104/pp.107.110544 PMID: 18218974

45. Sasaki K, Tsurumaru Y, Yamamoto H, Yazaki K. Molecular characterization of a membrane-bound prenyltransferase specific for isoflavone from Sophora flavescens. J Biol Chem. 2011; 286: 24125-24134. doi: 10.1074/jbc.M111.244426 PMID: 21576242

46. Munakata R, Inoue T, Koeduka T, Karamat F, Olry A, Sugiyama A, et al. Molecular cloning and characterization of a geranyl diphosphate-specific aromatic prenyltransferase from lemon. Plant Physiol. 2014; 166: 80-90. doi: 10.1104/pp.114.246892 PMID: 25077796

47. Li H, Ban Z, Qin H, Ma L, King AJ, Wang G. A heteromeric membrane-bound prenyltransferase complex from hop catalyzes three sequential aromatic prenylations in the bitter acid pathway. Plant Physiol. 2015; 167: 650-659. doi: 10.1104/pp.114.253682 PMID: 25564559

48. Bourgaud F, Brunel M, Guckert A, Forlot P. Effect of nitrogen nutrition and environmental conditions on the production of pharmaceutically useful metabolites by Psoralea cinerea. Eur J Agron. 1992; 1: 37-43.

49. Curk F. Citrus gene pool organization and nuclear genomic interspecific admixture of cultivated citrus. 2014.

50. Fujita T, Kawase A, Niwa T, Tomohiro N, Masuda M, Matsuda H, et al. Comparative evaluation of 12 immature citrus fruit extracts for the inhibition of cytochrome P450 isoform activities. Biol Pharm Bull. 2008; 31: 925-930. PMID: 18451520

51. Raquet N, Schrenk D. Application of the equivalency factor concept to the phototoxicity and -genotoxicity of furocoumarin mixtures. Food Chem Toxicol Int J Publ Br Ind Biol Res Assoc. 2014; 68: 257-266.

52. Wangensteen H, Molden E, Christensen H, Malterud KE. Identification of epoxybergamottin as a CYP3A4 inhibitor in grapefruit peel. Eur J Clin Pharmacol. 2003; 58: 663-668. PMID: 12610742

53. Pathak MA, Fellman JH, Kaufman KD. The effect of structural alterations on the erythemal activity of furocoumarins: psoralens. J Invest Dermatol. 1960; 35: 165-183. PMID: 13733305

54. Colombain M, Goll V, Muyard F, Girard C, Bévalot F, Richert L. A bioassay using the human hepato-blastoma cell line HepG2 for detecting phototoxicity of furocoumarins. Planta Med. 2001; 67: 644-646. doi: 10.1055/s-2001-17364 PMID: 11582543