A nuclear phylogenetic analysis: 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

Pangenome analysis provides insight into the evolution of the orange subfamily and a key gene for citric acid accumulation in citrus fruits

The orange subfamily (Aurantioideae) contains several Citrus species cultivated worldwide, such as sweet orange and lemon. The origin of Citrus species has long been debated and less is known about the Aurantioideae. Here, we compiled the genome sequences of 314 accessions, de novo assembled the genomes of 12 species and constructed a graph-based pangenome for Aurantioideae. Our analysis indicates that the ancient Indian Plate is the ancestral area for Citrus-related genera and that South Central China is the primary center of origin of the Citrus genus. We found substantial variations in the sequence and expression of the PH4 gene in Citrus relative to Citrus-related genera. Gene editing and biochemical experiments demonstrate a central role for PH4 in the accumulation of citric acid in citrus fruits. This study provides insights into the origin and evolution of the orange subfamily and a regulatory mechanism underpinning the evolution of fruit taste.

Whole-genome sequencing and variant discovery of Citrus reticulata "Kinnow" from Pakistan

Citrus is a source of nutritional and medicinal advantages, cultivated worldwide with major groups of sweet oranges, mandarins, grapefruits, kumquats, lemons and limes. Pakistan produces all major citrus groups with mandarin (Citrus reticulata) being the prominent group that includes local commercial cultivars Feutral's Early, Dancy, Honey, and Kinnow. The present study designed to understand the genetic architecture of this unique variety of Citrus reticulata 'Kinnow.' The whole-genome resequencing and variant calling was performed to map the genomic variability that might be responsible for its particular characteristics like taste, seedlessness, juice content, thickness of peel, and shelf-life. A total of 139,436,350 raw sequence reads were generated with 20.9 Gb data in Fastq format having 98% effectiveness and 0.2% base call error rate. Overall, 3,503,033 SNPs, 176,949 MNPs, 323,287 INS, and 333,083 DEL were identified using the GATK4 variant calling pipeline against Citrus clementina. Furthermore, g:Profiler was applied for annotating the newly found variants, harbor genes/transcripts and their involved pathways. A total of 73,864 transcripts harbors 4,336,352 variants, most of the observed variants were predicted in non-coding regions and 1009 transcripts were found well annotated by different databases. Out of total aforementioned transcripts, 588 involved in biological processes, 234 in molecular functions and 167 transcripts in cellular components. In a nutshell, 18,153 high impact variants and 216 genic variants found in the current study, which may be used after its functional validation for marker-assisted breeding programs of "Kinnow" to propagate its valued traits for the improvement of contemporary citrus varieties in the region.

Comprehensive genome-wide identification and transferability of chromosome-specific highly variable microsatellite markers from citrus species

Citrus species among the most important and widely consumed fruit in the world due to Vitamin C, essential oil glands, and flavonoids. Highly variable simple sequence repeats (SSR) markers are one of the most informative and versatile molecular markers used in perennial tree genetic research. SSR survey of Citrus sinensis and Citrus maxima were identified perfect SSRs spanning nine chromosomes. Furthermore, we categorized all SSR motifs into three major classes based on their tract lengths. We designed and validated a class I SSRs in the C. sinensis and C. maxima genome through electronic polymerase chain reaction (ePCR) and found 83.89% in C. sinensis and 78.52% in C. maxima SSRs producing a single amplicon. Then, we selected extremely variable SSRs (> 40 nt) from the ePCR-verified class I SSRs and in silico validated across seven draft genomes of citrus, which provided us a subset of 84.74% in C. sinensis and 77.53% in C. maxima highly polymorphic SSRs. Out of these, 129 primers were validated on 24 citrus genotypes through wet-lab experiment. We found 127 (98.45%) polymorphic HvSSRs on 24 genotypes. The utility of the developed HvSSRs was demonstrated by analysing genetic diversity of 181 citrus genotypes using 17 HvSSRs spanning nine citrus chromosomes and were divided into 11 main groups through 17 HvSSRs. These chromosome-specific SSRs will serve as a powerful genomic tool used for future QTL mapping, molecular breeding, investigation of population genetic diversity, comparative mapping, and evolutionary studies among citrus and other relative genera/species.

Reuse of Food Waste: The Chemical Composition and Health Properties of Pomelo (Citrus maxima) Cultivar Essential Oils

The aim of the present study is to investigate the chemical profile, antioxidant activity, carbohydrate-hydrolysing enzyme inhibition, and hypolipidemic effect of essential oils (EOs) extracted from Sicilian Citrus maxima (pomelo) flavedo. Using gas-chromatography-mass spectrometry analysis (GC-MS) we analysed the Eos of five cultivars of C. maxima, namely, ‘Chadock’, ‘Maxima’, ‘Pyriformis’, ‘Terracciani’, and ‘Todarii’, and their blends. The antioxidant activity was performed by using a multi-target approach using 2,2′-Azino-Bis-3-Ethylbenzothiazoline-6-Sulfonic acid (ABTS), 2,2-Diphenyl-1-picrylhydrazyl (DPPH), ferric reducing ability power (FRAP), and β-carotene bleaching tests. The α-amylase, α-glucosidase, and lipase-inhibitory activities were also assessed. GC-MS analyses revealed D-limonene as the main monoterpene hydrocarbon in all cultivars, albeit with different percentages in the range of 21.72–71.13%. A good content of oxygenated monoterpenes was detected for all cultivars, especially for ‘Todarii’. The analysis of the principal components (PCA), and related clusters (HCA), was performed to find chemo-diversity among the analysed samples. EOs from ‘Chadock’ and ‘Maxima’ were statistically similar to each other, and they differed from P3 in the smaller amount of sesquiterpene hydrocarbons, while the oils from ‘Terracciani’ and ‘Todarii’ were found to be chemically and statistically different. ‘Chadock’ EO was the most active to scavenge radicals (IC₅₀ values of 22.24 and 27.23 µg/mL in ABTS and DPPH tests, respectively). ‘Terracciani’ EO was the most active against both lipase and α-amylase, whereas the blends obtained by the combination (1:1 v/v) of C. maxima ‘Maxima’ + ‘Todarii’ were the most active against α-glucosidase. Generally, the blends did not exert a unique behaviour in potentiating or reducing the bioactivity of the pomelo EOs.

Potential application of simple easy-to-use insertion-deletion (InDel) markers in citrus cultivar identification

We previously developed insertion-deletion (InDel) markers that distinguish three genotypes (two homozygous and one heterozygous) of diverse citrus cultivars. These InDel markers were codominant and could be clearly detected by using simple agarose gel electrophoresis. We sought to establish a method for cultivar identification using these 28 InDel markers to genotype 31 citrus cultivars. The results revealed that a minimum of 6 markers were required to identify individuals using the three-genotype classification method. Furthermore, we found that a simple method for distinguishing between two genotypes (homozygous and heterozygous) could be used to identify individuals using a minimum of 7 markers. Our findings provide a basis for the development of simple and rapid citrus cultivar identification methods.

Chemical Variability of Peel and Leaf Essential Oils in the Citrus Subgenus Papeda (Swingle) and Few Relatives

The Papeda Citrus subgenus includes several species belonging to two genetically distinct groups, containing mostly little-exploited wild forms of citrus. However, little is known about the potentially large and novel aromatic diversity contained in these wild citruses. In this study, we characterized and compared the essential oils obtained from peels and leaves from representatives of both Papeda groups, and three related hybrids. Using a combination of GC, GC-MS, and ¹³C-NMR spectrometry, we identified a total of 60 compounds in peel oils (PO), and 76 compounds in leaf oils (LO). Limonene was the major component in almost all citrus PO, except for C. micrantha and C. hystrix, where β-pinene dominated (around 35%). LO composition was more variable, with different major compounds among almost all samples, except for two citrus pairs: C. micrantha/C. hystrix and two accessions of C. ichangensis. In hybrid relatives, the profiles were largely consistent with their Citrus/Papeda parental lineage. This high chemical diversity, not only among the sections of the subgenus Papeda, but also between species and even at the intraspecific level, suggests that Papeda may be an important source of aroma diversity for future experimental crosses with field crop species.

Strategies to Produce Grapefruit-Like Citrus Varieties With a Low Furanocoumarin Content and Distinctive Flavonoid Profiles

Pummelos and hybrids, such as grapefruits, have high furanocoumarin and low flavonoid contents. Furanocoumarins interact negatively with certain drugs, while flavonoids are antioxidant compounds with health benefits. To obtain new grapefruit-like varieties with low furanocoumarin and high flavonoid contents, diploid and triploid hybrid populations from crosses between diploid and tetraploid "Clemenules" clementine and diploid "Pink" pummelo were recovered and analyzed. With regard to furanocoumarins, triploids produce less bergapten, bergamottin and 6,7-DHB than diploids. Regarding flavonoids, triploids yielded more eriocitrin, narirutin, hesperidin and neohesperidin than diploids, whereas no differences were observed in neoeriocitrin and naringin. These results indicate that, the strategy to recover triploid hybrids by 4x × 2x crosses is more appropriate than the recovery of diploid hybrids by 2x × 2x crosses for obtaining grapefruit-like varieties of citrus with lower furanocoumarin and higher flavonoid contents.

Genomic insights into citrus domestication and its important agronomic traits

Citrus originated in Southeast Asia, and it has become one of the most important fruit crops worldwide. Citrus has a long and obscure domestication history due to its clonal propagation, long life cycle, wide sexual compatibility, and complex genetic background. As the genomic information of both wild and cultivated citrus becomes available, their domestication history and underlying traits or genes are becoming clear. This review outlines the genomic features of wild and cultivated species. We propose that the reduction of citric acid is a critical trait for citrus domestication. The genetic model representing the change during domestication may be associated with a regulatory complex known as WD-repeat-MYB-bHLH-WRKY (WMBW), which is involved in acidification and anthocyanin accumulation. The reduction in or loss of anthocyanins may be due to a hitchhiking effect of fruit acidity selection, in which mutation occurs in the common regulator of these two pathways in some domesticated types. Moreover, we have summarized the domestication traits and candidate genes for breeding purposes. This review represents a comprehensive summary of the genes controlling key traits of interest, such as acidity, metabolism, and disease resistance. It also sheds light on recent advances in early flowering from transgenic studies and provides a new perspective for fast breeding of citrus. Our review lays a foundation for future research on fruit acidity, flavor, and disease resistance in citrus.

Integration of mandarin (Citrus reticulata) cytogenetic map with its genome sequence

Citrus is an extremely important genus in terms of world fruit production. Despite its economic importance and the small genome sizes of its species (2n = 18, 1C = 430 ± 68 Mbp), entire genomic assemblies have only recently become available for some of its representatives. Together with the previous CMA/DAPI banding and fluorescence in situ hybridization (FISH) in the group, these data are important for understanding the complex relationships between its species and for assisting breeding programs. To anchor genomic data with the cytogenetic map of mandarin (Citrus reticulata), the parental species of several economically important hybrids such as sweet orange and clementine, 18 BAC (bacterial artificial chromosome) clones were used. Eleven clementine BACs were positioned by BAC-FISH, doubling the number of chromosome markers so far available for BAC-FISH in citrus. Additionally, six previously mapped BACs were end-sequenced, allowing, together with one BAC previously sequenced, their assignment to scaffolds and the subsequent integration of chromosomes and the genome assembly. This study therefore established correlations between mandarin scaffolds and chromosomes, allowing further structural genomic and comparative study with the sweet orange genome, as well as insights into the chromosomal evolution of the group.

Male and female inheritance patterns in tetraploid 'Moncada' mandarin

KEY MESSAGE

Tetraploid `Moncada´ mandarin, used as male and female in interploidy hybridizations, displays mainly tetrasomic inheritance for most LGs, with slight variations according to the direction of the crossing. Triploid-breeding programs in citrus are key tool to develop seedless cultivars. Obtaining triploid citrus hybrids may be achieved through different strategies, such as the exploitation of female unreduced gamete in crosses between diploid parents and diploid by tetraploid sexual hybridizations, in which tetraploid genotypes can be used as male or female parents. Genetic configuration of triploid populations from interploid crosses greatly depends on the chromosomic segregation mode of the tetraploid parent used. Here, we have analyzed the inheritance of the tetraploid 'Moncada' mandarin and compared the genetic structures of the resulting gametes when used as male and as female parent. The preferential chromosome pairing rate is calculated from the parental heterozygosity restitution (PHR) of codominant molecular markers, indicating the proportion between disomic and tetrasomic segregation. Tetraploid 'Moncada' both as female and male parent largely exhibited tetrasomic segregation. However, as female parent, one linkage group (LG8) showed intermediate segregation with tendency towards tetrasomic inheritance, while another linkage group (LG4) evidenced a clear intermediate segregation. On the other hand, when used as male parent two linkage groups (LG5 and LG6) showed values that fit an intermediate inheritance model with tetrasomic tendency. Significant doubled reduction (DR) rates were observed in five linkage groups as female parent, and in six linkage groups as male parent. The new knowledge generated here will serve to define crossing strategies in citrus improvement programs to efficiently obtain new varieties of interest in the global fresh consumption market.

Resistance to 'Candidatus Liberibacter asiaticus,' the Huanglongbing Associated Bacterium, in Sexually and/or Graft-Compatible Citrus Relatives

Huanglongbing (HLB) is the most destructive, yet incurable disease of citrus. Finding sources of genetic resistance to HLB-associated 'Candidatus Liberibacter asiaticus' (Las) becomes strategic to warrant crop sustainability, but no resistant Citrus genotypes exist. Some Citrus relatives of the family Rutaceae, subfamily Aurantioideae, were described as full-resistant to Las, but they are phylogenetically far, thus incompatible with Citrus. Partial resistance was indicated for certain cross-compatible types. Moreover, other genotypes from subtribe Citrinae, sexually incompatible but graft-compatible with Citrus, may provide new rootstocks able to restrict bacterial titer in the canopy. Use of seedlings from monoembryonic species and inconsistencies in previous reports likely due to Las recalcitrance encouraged us to evaluate more accurately these Citrus relatives. We tested for Las resistance a diverse collection of graft-compatible Citrinae species using an aggressive and consistent challenge-inoculation and evaluation procedure. Most Citrinae species examined were either susceptible or partially resistant to Las. However, Eremocitrus glauca and Papua/New Guinea Microcitrus species as well as their hybrids and those with Citrus arose here for the first time as full-resistant, opening the way for using these underutilized genotypes as Las resistance sources in breeding programs or attempting using them directly as possible new Las-resistant Citrus rootstocks or interstocks.

Analysis of Crossover Events and Allele Segregation Distortion in Interspecific Citrus Hybrids by Single Pollen Genotyping

In citrus, a classical method of studying crossovers and segregation distortion (SD) is the genetic analysis of progenies. A new strategy combining fluorescence-activated cell sorting and whole genome amplification of haploid pollen nuclei with a large set of molecular markers, offers the opportunity to efficiently determine the frequency of crossovers and the identification of SD without the need to generate segregating populations. Here we have analyzed meiotic crossover events in a pollen nuclei population from "Eureka" lemon and the allelic SD was evaluated in a pollen nuclei population from a clementine × sweet orange hybrid ("CSO"). Data obtained from the "CSO" pollen nuclei population were compared to those obtained from genotyping of a segregating population ("RTSO") arising from a hand-made sexual hybridization between diploid non apomictic selected tangor (mandarin × sweet orange; "RTO" tangor) as female parent pollinated with "CSO" tangor as male parent. The analysis of crossovers rates on chromosome 1 revealed the presence of up to five crossovers events on one arm and four on the corresponding other arm, with an average of 1.97 crossovers per chromosome while no crossover events were observed in five "Eureka" lemon pollen nuclei. The rate of SD observed in "CSO" pollen nuclei (13.8%) was slightly lower than that recovered in the "RTSO" population (20.7%). In the pollen nuclei population, SD was found on linkage group (LG) 2, while the "RTSO" population showed SD on LGs 2 and 7. Potential male gametic selection mechanisms were distinguished in pollen grains, while in the population, mechanisms of gametophytic selection and/or zygotic selection were observed. This methodology is a very useful tool to facilitate research focused on the reproductive biology of citrus and study the mechanisms that affect crossovers and SD.

Genetic, morphological and chemical investigations reveal the genetic origin of Pompia (C. medica tuberosa Risso & Poiteau) - An old endemic Sardinian citrus fruit

Citrus fruits have been introduced to the Mediterranean area from Asia for centuries and spontaneous crosses have generated several hybrid forms, some of which have had agricultural or industrial success while others have remained niche food or ornamental products, or have disappeared. Pompia (C. medica tuberosa Risso & Poiteau) is an old endemic citrus fruit from Sardinia of unknown genetic origin. Initial phenotypic and molecular characterizations revealed a high degree of similarity with lemon (C. limon (L.) Burm.) and citron (C. medica L.). To identify the ancestors of Pompia, 70 citrus species of the Citrus genus were genotyped with 36 codominant molecular markers (SSR and InDel) of nuclear and cytoplasmic genomes. Diversity analysis and allelic comparisons between each citrus species at each locus indicated that Pompia resembles lemon and limonette of Marrakech, i.e. the result of a cross between sour orange (C. aurantium L.) and citron, where citron was the pollinator. Two Italian citron varieties were identified as potential male parents, i.e. Diamante and Common Poncire. However, we were unable to differentiate varieties of sour oranges because varietal diversification in this horticultural group resulted from DNA sequence variations that SSR or InDel markers could not reveal. Rhob el Arsa and Poncire de Collioure were found to be two synonyms of Pompia. Pompia appeared to be equally distinct from citron, lemon and sour orange based on the overall analysis of the fruit, leaf and seed phenotype, and juice chemical composition. At the leaf level, the Pompia essential oil (EO) composition is close to that of citron whereas the zest is much closer to that of sour orange.

Genome sequencing and CRISPR/Cas9 gene editing of an early flowering Mini-Citrus (Fortunella hindsii)

Hongkong kumquat (Fortunella hindsii) is a wild citrus species characterized by dwarf plant height and early flowering. Here, we identified the monoembryonic F. hindsii (designated as 'Mini-Citrus') for the first time and constructed its selfing lines. This germplasm constitutes an ideal model for the genetic and functional genomics studies of citrus, which have been severely hindered by the long juvenility and inherent apomixes of citrus. F. hindsii showed a very short juvenile period (~8 months) and stable monoembryonic phenotype under cultivation. We report the first de novo assembled 373.6 Mb genome sequences (Contig-N50 2.2 Mb and Scaffold-N50 5.2 Mb) for F. hindsii. In total, 32 257 protein-coding genes were annotated, 96.9% of which had homologues in other eight Citrinae species. The phylogenomic analysis revealed a close relationship of F. hindsii with cultivated citrus varieties, especially with mandarin. Furthermore, the CRISPR/Cas9 system was demonstrated to be an efficient strategy to generate target mutagenesis on F. hindsii. The modifications of target genes in the CRISPR-modified F. hindsii were predominantly 1-bp insertions or small deletions. This genetic transformation system based on F. hindsii could shorten the whole process from explant to T1 mutant to about 15 months. Overall, due to its short juvenility, monoembryony, close genetic background to cultivated citrus and applicability of CRISPR, F. hindsii shows unprecedented potentials to be used as a model species for citrus research.

Phylogenetic Origin of Primary and Secondary Metabolic Pathway Genes Revealed by C. maxima and C. reticulata Diagnostic SNPs

Modern cultivated Citrus species and varieties result from interspecific hybridization between four ancestral taxa. Among them, Citrus maxima and Citrus reticulata, closely associated with the pummelo and mandarin horticultural groups, respectively, were particularly important as the progenitors of sour and sweet oranges (Citrus aurantium and Citrus sinensis), grapefruits (Citrus paradisi), and hybrid types resulting from modern breeding programs (tangors, tangelos, and orangelos). The differentiation between the four ancestral taxa and the phylogenomic structure of modern varieties widely drive the phenotypic diversity's organization. In particular, strong phenotypic differences exist in the coloration and sweetness and represent important criteria for breeders. In this context, focusing on the genes of the sugar, carotenoid, and chlorophyll biosynthesis pathways, the aim of this work was to develop a set of diagnostic single-nucleotide polymorphism (SNP) markers to distinguish the ancestral haplotypes of C. maxima and C. reticulata and to provide information at the intraspecific diversity level (within C. reticulata or C. maxima). In silico analysis allowed the identification of 3,347 SNPs from selected genes. Among them, 1,024 were detected as potential differentiation markers between C. reticulata and C. maxima. A total of 115 SNPs were successfully developed using a competitive PCR technology. Their transferability among all Citrus species and the true citrus genera was very good, with only 0.87% of missing data. The ancestral alleles of the SNPs were identified, and we validated the usefulness of the developed markers for tracing the ancestral haplotype in large germplasm collections and sexually recombined progeny issued from the C. reticulata/C. maxima admixture gene pool. These markers will pave the way for targeted association studies based on ancestral haplotypes.

Genotyping by sequencing can reveal the complex mosaic genomes in gene pools resulting from reticulate evolution: a case study in diploid and polyploid citrus

BACKGROUND AND AIMS

Reticulate evolution, coupled with reproductive features limiting further interspecific recombinations, results in admixed mosaics of large genomic fragments from the ancestral taxa. Whole-genome sequencing (WGS) data are powerful tools to decipher such complex genomes but still too costly to be used for large populations. The aim of this work was to develop an approach to infer phylogenomic structures in diploid, triploid and tetraploid individuals from sequencing data in reduced genome complexity libraries. The approach was applied to the cultivated Citrus gene pool resulting from reticulate evolution involving four ancestral taxa, C. maxima, C. medica, C. micrantha and C. reticulata.

METHODS

A genotyping by sequencing library was established with the restriction enzyme ApeKI applying one base (A) selection. Diagnostic single nucleotide polymorphisms (DSNPs) for the four ancestral taxa were mined in 29 representative varieties. A generic pipeline based on a maximum likelihood analysis of the number of read data was established to infer ancestral contributions along the genome of diploid, triploid and tetraploid individuals. The pipeline was applied to 48 diploid, four triploid and one tetraploid citrus accessions.

KEY RESULTS

Among 43 598 mined SNPs, we identified a set of 15 946 DSNPs covering the whole genome with a distribution similar to that of gene sequences. The set efficiently inferred the phylogenomic karyotype of the 53 analysed accessions, providing patterns for common accessions very close to that previously established using WGS data. The complex phylogenomic karyotypes of 21 cultivated citrus, including bergamot, triploid and tetraploid limes, were revealed for the first time.

CONCLUSIONS

The pipeline, available online, efficiently inferred the phylogenomic structures of diploid, triploid and tetraploid citrus. It will be useful for any species whose reproductive behaviour resulted in an interspecific mosaic of large genomic fragments. It can also be used for the first generations of interspecific breeding schemes.

Facing Climate Change: Biotechnology of Iconic Mediterranean Woody Crops

The Mediterranean basin is especially sensitive to the adverse outcomes of climate change and especially to variations in rainfall patterns and the incidence of extremely high temperatures. These two concurring adverse environmental conditions will surely have a detrimental effect on crop performance and productivity that will be particularly severe on woody crops such as citrus, olive and grapevine that define the backbone of traditional Mediterranean agriculture. These woody species have been traditionally selected for traits such as improved fruit yield and quality or alteration in harvesting periods, leaving out traits related to plant field performance. This is currently a crucial aspect due to the progressive and imminent effects of global climate change. Although complete genome sequence exists for sweet orange (Citrus sinensis) and clementine (Citrus clementina), olive tree (Olea europaea) and grapevine (Vitis vinifera), the development of biotechnological tools to improve stress tolerance still relies on the study of the available genetic resources including interspecific hybrids, naturally occurring (or induced) polyploids and wild relatives under field conditions. To this respect, post-genomic era studies including transcriptomics, metabolomics and proteomics provide a wide and unbiased view of plant physiology and biochemistry under adverse environmental conditions that, along with high-throughput phenotyping, could contribute to the characterization of plant genotypes exhibiting physiological and/or genetic traits that are correlated to abiotic stress tolerance. The ultimate goal of precision agriculture is to improve crop productivity, in terms of yield and quality, making a sustainable use of land and water resources under adverse environmental conditions using all available biotechnological tools and high-throughput phenotyping. This review focuses on the current state-of-the-art of biotechnological tools such as high throughput -omics and phenotyping on grapevine, citrus and olive and their contribution to plant breeding programs.

Assessing Ploidy Level Analysis and Single Pollen Genotyping of Diploid and Euploid Citrus Genotypes by Fluorescence-Activated Cell Sorting and Whole-Genome Amplification

Flow cytometry is widely used to determine genome size and ploidy level in plants. This technique, when coupled with fluorescence-activated cell sorting (FACS), whole genome amplification and genotyping (WGA), opens up new opportunities for genetic studies of individualized nuclei. This strategy was used to analyze the genetic composition of single pollen nuclei of different citrus species. The flow cytometry and microscope observations allowed us to differentiate the populations of pollen nuclei present in the diploid and euploid genotypes analyzed, showing that citrus has binuclear pollen. We have identified in the "CSO" tangor an additional nuclei population composed by the vegetative plus generative nuclei. Genotyping of this nuclei population revealed that vegetative and generative nuclei show the same genetic configuration. In addition, we have demonstrated the presence of unreduced gametes in the diploid genotype "Mexican lime." Genomic amplification is a robust method for haploid nuclei genotyping with several molecular markers, whereas in diploid nuclei using heterozygous markers showed a bias towards one of the two alleles, limiting the use of this tool in this type of nuclei. We further discuss the importance and applications of single pollen genotyping in citrus genetic studies.

Molecular Characterization and Stress Tolerance Evaluation of New Allotetraploid Somatic Hybrids Between Carrizo Citrange and Citrus macrophylla W. rootstocks

Polyploidy is one of the main forces that drives the evolution of plants and provides great advantages for breeding. Somatic hybridization by protoplast fusion is used in citrus breeding programs. This method allows combining the whole parental genomes in a single genotype, adding complementary dominant characters, regardless of parental heterozygosity. It also contributes to surpass limitations imposed by reproductive biology and quickly generates progenies that combine the required traits. Two allotetraploid somatic hybrids recovered from the citrus rootstocks-Citrus macrophylla (CM) and Carrizo citrange (CC)-were characterized for morphology, genome composition using molecular markers (SNP, SSR, and InDel), and their tolerance to iron chlorosis, salinity, and Citrus tristeza virus (CTV). Both hybrids combine the whole parental genomes even though the loss of parental alleles was detected in most linkage groups. Mitochondrial genome was inherited from CM in both the hybrids, whereas recombination was observed for chloroplastic genome. Thus, somatic hybrids differ from each other in their genome composition, indicating that losses and rearrangements occurred during the fusion process. Both inherited the tolerance to stem pitting caused by CTV from CC, are tolerant to iron chlorosis such as CM, and have a higher tolerance to salinity than the sensitive CC. These hybrids have potential as improved rootstocks to grow citrus in areas with calcareous and saline soils where CTV is present, such as the Mediterranean region. The provided knowledge on the effects of somatic hybridization on the genome composition, anatomy, and physiology of citrus rootstocks will be key for breeding programs that aim to address current and future needs of the citrus industry.

Endogenous pararetrovirus sequences are widely present in Citrinae genomes

Endogenous pararetroviruses (EPRVs) are characterized in several plant genomes and their biological effects have been reported. In this study, hundreds of EPRV segments were identified in six Citrinae genomes. A total of 1034 EPRV segments were identified in the genomes of sweet orange, 2036 in pummelo, 598 in clementine mandarin, 752 in Ichang papeda, 2060 in citron and 245 in atalantia. Genomic analysis indicated that EPRV segments tend to cluster as hot spots in the genomes, particularly on chromosome 2 and 5. Large numbers of simple repeats and transposable elements were identified in the 2-kb flanking regions of the EPRV segments. Comparative genomic analysis and PCR experiments showed that there are highly conserved EPRV segments and species-specific EPRV segments between the Citrinae genomes. Phylogenetic analysis suggested that the integration events of EPRVs could initiate in a common progenitor of Citrinae species and repeatedly occur during the Citrinae divergence.

Origin of C. latifolia and C. aurantiifolia triploid limes: the preferential disomic inheritance of doubled-diploid 'Mexican' lime is consistent with an interploid hybridization hypothesis

Background and Aims

Two main types of triploid limes are produced worldwide. The 'Tahiti' lime type (Citrus latifolia) is predominant, while the 'Tanepao' type (C. aurantiifolia) is produced to a lesser extent. Both types result from natural interspecific hybridization involving a diploid gamete of C. aurantiifolia 'Mexican' lime type (itself a direct interspecific C. micrantha × C. medica hybrid). The meiotic behaviour of a doubled-diploid 'Mexican' lime, the interspecific micrantha/medica recombination and the resulting diploid gamete structures were analysed to investigate the possibility that 'Tahiti' and 'Tanepao' varieties are derived from natural interploid hybridization.

Methods

A population of 85 tetraploid hybrids was established between a doubled-diploid clementine and a doubled-diploid 'Mexican' lime and used to infer the genotypes of 'Mexican' lime diploid gametes. Meiotic behaviour was studied through combined segregation analysis of 35 simple sequenbce repeat (SSR) and single nucleotide polymorphismn (SNP) markers covering the nine citrus chromosomes and cytogenetic studies. It was supplemented by pollen viability assessment.

Key Results

Pollen viability of the doubled-diploid Mexican lime (64 %) was much higher than that of the diploid. On average, 65 % of the chromosomes paired as bivalents and 31.4 % as tetravalents. Parental heterozygosity restitution ranged from 83 to 99 %. Disomic inheritance with high preferential pairing values was deduced for three chromosomes. Intermediate inheritances, with disomic trend, were found for five chromosomes, and an intermediate inheritance was observed for one chromosome. The average effective interspecific recombination rate was low (1.2 cM Mb-1).

Conclusion

The doubled-diploid 'Mexican' lime had predominantly disomic segregation, producing interspecific diploid gamete structures with high C. medica/C. micrantha heterozygosity, compatible with the phylogenomic structures of triploid C. latifolia and C. aurantiifolia varieties. This disomic trend limits effective interspecific recombination and diversity of the diploid gamete population. Interploid reconstruction breeding using doubled-diploid lime as one parent is a promising approach for triploid lime diversification.

Genotyping by sequencing reveals the interspecific C. maxima / C. reticulata admixture along the genomes of modern citrus varieties of mandarins, tangors, tangelos, orangelos and grapefruits

The mandarin horticultural group is an important component of world citrus production for the fresh fruit market. This group formerly classified as C. reticulata is highly polymorphic and recent molecular studies have suggested that numerous cultivated mandarins were introgressed by C. maxima (the pummelos). C. maxima and C. reticulata are also the ancestors of sweet and sour oranges, grapefruit, and therefore of all the "small citrus" modern varieties (mandarins, tangors, tangelos) derived from sexual hybridization between these horticultural groups. Recently, NGS technologies have greatly modified how plant evolution and genomic structure are analyzed, moving from phylogenetics to phylogenomics. The objective of this work was to develop a workflow for phylogenomic inference from Genotyping By Sequencing (GBS) data and to analyze the interspecific admixture along the nine citrus chromosomes for horticultural groups and recent varieties resulting from the combination of the C. reticulata and C. maxima gene pools. A GBS library was established from 55 citrus varieties, using the ApekI restriction enzyme and selective PCR to improve the read depth. Diagnostic polymorphisms (DPs) of C. reticulata/C. maxima differentiation were identified and used to decipher the phylogenomic structure of the 55 varieties. The GBS approach was powerful and revealed 30,289 SNPs and 8,794 Indels with 12.6% of missing data. 11,133 DPs were selected covering the nine chromosomes with a higher density in genic regions. GBS combined with the detection of DPs was powerful for deciphering the "phylogenomic karyotypes" of cultivars derived from admixture of the two ancestral species after a limited number of interspecific recombinations. All the mandarins, mandarin hybrids, tangelos and tangors analyzed displayed introgression of C. maxima in different parts of the genome. C. reticulata/C. maxima admixture should be a major component of the high phenotypic variability of this germplasm opening up the way for association studies based on phylogenomics.

Unreduced Megagametophyte Production in Lemon Occurs via Three Meiotic Mechanisms, Predominantly Second-Division Restitution

Unreduced (2n) gametes have played a pivotal role in polyploid plant evolution and are useful for sexual polyploid breeding in various species, particularly for developing new seedless citrus varieties. The underlying mechanisms of 2n gamete formation were recently revealed for Citrus reticulata but remain poorly understood for other citrus species, including lemon (C. limon [L.] Burm. f.). Here, we investigated the frequency and causal meiotic mechanisms of 2n megagametophyte production in lemon. We genotyped 48progeny plants of two lemon genotypes, "Eureka Frost" and "Fino", using 16 Simple Sequence Repeat (SSR) and 18 Single Nucleotide Polymorphism (SNP) markers to determine the genetic origin of the progenies and the underlying mechanisms for 2n gamete formation. We utilized a maximum-likelihood method based on parental heterozygosity restitution (PHR) of centromeric markers and analysis of PHR patterns along the chromosome. The frequency of 2n gamete production was 4.9% for "Eureka Frost" and 8.3% for "Fino", with three meiotic mechanisms leading to 2n gamete formation. We performed the maximum-likelihood method at the individual level via centromeric marker analysis, finding that 88% of the hybrids arose from second-division restitution (SDR), 7% from first-division restitution (FDR) or pre-meiotic doubling (PRD), and 5% from post-meiotic genome doubling (PMD). The pattern of PHR along LG1 confirmed that SDR is the main mechanism for 2n gamete production. Recombination analysis between markers in this LG revealed partial chiasma interference on both arms. We discuss the implications of these restitution mechanisms for citrus breeding and lemon genetics.

Bacterial phytopathogen infection disrupts belowground plant indirect defense mediated by tritrophic cascade

Plants can defend themselves against herbivores through activation of defensive pathways and attraction of third-trophic-level predators and parasites. Trophic cascades that mediate interactions in the phytobiome are part of a larger dynamic including the pathogens of the plant itself, which are known to greatly influence plant defenses. As such, we investigated the impact of a phloem-limited bacterial pathogen, Candidatus Liberibacter asiaticus (CLas), in cultivated citrus rootstock on a well-studied belowground tritrophic interaction involving the attraction of an entomopathogenic nematode (EPN), Steinernema diaprepesi, to their root-feeding insect hosts, Diaprepes abbreviatus larvae. Using belowground olfactometers, we show how CLas infection interferes with this belowground interaction by similarly inducing the release of a C12 terpene, pregeijerene, and disconnecting the association of the terpene with insect presence. D. abbreviatus larvae that were not feeding but in the presence of a CLas-infected plant were more likely to be infected by EPN than those near uninfected plants. Furthermore, nonfeeding larvae associated with CLas-infected plants were just as likely to be infected by EPN as those near noninfected plants with D. abbreviatus larval damage. Larvae of two weevil species, D. abbreviatus and Pachnaeus litus, were also more attracted to plants with infection than to uninfected plants. D. abbreviatus larvae were most active when exposed to pregeijerene at a concentration of 0.1 μg/μl. We attribute this attraction to CLas-infected plants to the same signal previously thought to be a herbivore-induced plant volatile specifically induced by root-feeding insects, pregeijerene, by assessing volatiles collected from the roots of infected plants and uninfected plants with and without feeding D. abbreviatus. Synthesis. Phytopathogens can influence the structuring of soil communities extending to the third trophic level. Field populations of EPN may be less effective at host-finding using pregeijerene as a cue in citrus grove agroecosystems with high presence of CLas infection.

Anticancer Potential of Citrus Juices and Their Extracts: A Systematic Review of Both Preclinical and Clinical Studies

Background: During the last decades, a huge body of evidence has been accumulated suggesting that Citrus fruits and their juices might have a role in preventing many diseases including cancer.

Objective: To summarize the numerous evidences on the potential of Citrus juices and their extracts as anticancer agents.

Data sources: A systematic review of articles written in English using MEDLINE (1946-present), EMBASE (1974-present) and Web of Sciences (1970-present) was performed independently by two reviewers. Search terms included Citrus, Citrus aurantifolia, Citrus sinensis, Citrus paradisi, Citrus fruits, Citrus fruits extract, cancer, neoplasm, neoplasia, tumor, metastasis, carcinogenesis, proliferation. The last search was performed on March 16th, 2017.

Study selection: Study selection and systematic review were carried out in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement. Prior to the beginning of the review, Authors defined a checklist for inclusion criteria, thus including articles which meet the following: (i) published on peer-reviewed scientific journals; (ii) Citrus juice used alone; (iii) extracts derived from Citrus juice; (iii) for preclinical studies, an exposure time to Citrus juices and their extracts more than 24 h. Reviews, meta-analyses, conference abstracts and book chapters were excluded.

Data extraction: Three reviewers independently performed the extraction of articles.

Data synthesis: 22 papers met our inclusion criteria and were eligible for inclusion in the final review. According to the kind of study, the selected ones were further divided in preclinical (n = 20) and observational (n = 2) studies.

Conclusion: The studies discussed in this review strongly corroborate the role of Citrus juices and their derivatives as potential resource against cancer.

Changes in Anthocyanin Production during Domestication of Citrus

Mandarin (Citrus reticulata), citron (Citrus medica), and pummelo (Citrus maxima) are important species of the genus Citrus and parents of the interspecific hybrids that constitute the most familiar commercial varieties of Citrus: sweet orange, sour orange, clementine, lemon, lime, and grapefruit. Citron produces anthocyanins in its young leaves and flowers, as do species in genera closely related to Citrus, but mandarins do not, and pummelo varieties that produce anthocyanins have not been reported. We investigated the activity of the Ruby gene, which encodes a MYB transcription factor controlling anthocyanin biosynthesis, in different accessions of a range of Citrus species and in domesticated cultivars. A white mutant of lemon lacks functional alleles of Ruby, demonstrating that Ruby plays an essential role in anthocyanin production in Citrus Almost all the natural variation in pigmentation by anthocyanins in Citrus species can be explained by differences in activity of the Ruby gene, caused by point mutations and deletions and insertions of transposable elements. Comparison of the allelic constitution of Ruby in different species and cultivars also helps to clarify many of the taxonomic relationships in different species of Citrus, confirms the derivation of commercial varieties during domestication, elucidates the relationships within the subgenus Papeda, and allows a new genetic classification of mandarins.

Genomic Analysis of Terpene Synthase Family and Functional Characterization of Seven Sesquiterpene Synthases from Citrus sinensis

Citrus aroma and flavor, chief traits of fruit quality, are derived from their high content in essential oils of most plant tissues, including leaves, stems, flowers, and fruits. Accumulated in secretory cavities, most components of these oils are volatile terpenes. They contribute to defense against herbivores and pathogens, and perhaps also protect tissues against abiotic stress. In spite of their importance, our understanding of the physiological, biochemical, and genetic regulation of citrus terpene volatiles is still limited. The availability of the sweet orange (Citrus sinensis L. Osbeck) genome sequence allowed us to characterize for the first time the terpene synthase (TPS) family in a citrus type. CsTPS is one of the largest angiosperm TPS families characterized so far, formed by 95 loci from which just 55 encode for putative functional TPSs. All TPS angiosperm families, TPS-a, TPS-b, TPS-c, TPS-e/f, and TPS-g were represented in the sweet orange genome, with 28, 18, 2, 2, and 5 putative full length genes each. Additionally, sweet orange β-farnesene synthase, (Z)-β-cubebene/α-copaene synthase, two β-caryophyllene synthases, and three multiproduct enzymes yielding β-cadinene/α-copaene, β-elemene, and β-cadinene/ledene/allo-aromandendrene as major products were identified, and functionally characterized via in vivo recombinant Escherichia coli assays.

RAD-Seq analysis of typical and minor Citrus accessions, including Bhutanese varieties

We analyzed the reduced-representation genome sequences of Citrus species by double-digest restriction site-associated DNA sequencing (ddRAD-Seq) using 44 accessions, including typical and minor accessions, such as Bhutanese varieties. The results of this analysis using typical accessions were consistent with previous reports that citron, papeda, pummelo, and mandarin are ancestral species, and that most Citrus species are derivatives or hybrids of these four species. Citrus varieties often reproduce asexually and heterozygosity is highly conserved within each variety. Because this approach could readily detect conservation of heterozygosity, it was able to discriminate citrus varieties such as satsuma mandarin from closely related species. Thus, this method provides an inexpensive way to protect citrus varieties from unintended introduction and to prevent the provision of incorrect nursery stocks to customers. One Citrus variety in Bhutan was morphologically similar to Mexican lime and was designated as Himalayan lime. The current analysis confirmed the previous proposition that Mexican lime is a hybrid between papeda and citron, and also suggested that Himalayan lime is a probable hybrid between mandarin and citron. In addition to Himalayan lime, current analysis suggested that several accessions were formed by previously undescribed combinations.

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

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

Towards a molecular taxonomic key of the Aurantioideae subfamily using chloroplastic SNP diagnostic markers of the main clades genotyped by competitive allele-specific PCR

Background

Chloroplast DNA is a primary source of molecular variations for phylogenetic analysis of photosynthetic eukaryotes. However, the sequencing and analysis of multiple chloroplastic regions is difficult to apply to large collections or large samples of natural populations. The objective of our work was to demonstrate that a molecular taxonomic key based on easy, scalable and low-cost genotyping method should be developed from a set of Single Nucleotide Polymorphisms (SNPs) diagnostic of well-established clades. It was applied to the Aurantioideae subfamily, the largest group of the Rutaceae family that includes the cultivated citrus species.

Results

The publicly available nucleotide sequences of eight plastid genomic regions were compared for 79 accessions of the Aurantioideae subfamily to search for SNPs revealing taxonomic differentiation at the inter-tribe, inter-subtribe, inter-genus and interspecific levels. Diagnostic SNPs (DSNPs) were found for 46 of the 54 clade levels analysed. Forty DSNPs were selected to develop KASPar markers and their taxonomic value was tested by genotyping 108 accessions of the Aurantioideae subfamily. Twenty-seven markers diagnostic of 24 clades were validated and they displayed a very high rate of transferability in the Aurantioideae subfamily (only 1.2 % of missing data on average). The UPGMA from the validated markers produced a cladistic organisation that was highly coherent with the previous phylogenetic analysis based on the sequence data of the eight plasmid regions. In particular, the monophyletic origin of the “true citrus” genera plus Oxanthera was validated. However, some clarification remains necessary regarding the organisation of the other wild species of the Citreae tribe.

Conclusions

We validated the concept that with well-established clades, DSNPs can be selected and efficiently transformed into competitive allele-specific PCR markers (KASPar method) allowing cost-effective highly efficient cladistic analysis in large collections at subfamily level. The robustness of this genotyping method is an additional decisive advantage for network collaborative research. The availability of WGS data for the main “true citrus” species should soon make it possible to develop a set of DSNP markers allowing very fine resolution of this very important horticultural group.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0426-x) contains supplementary material, which is available to authorized users.

Patterns of genetic structure and evidence of gene flow among Tunisian Citrus species based on informative nSSR markers

This study investigates the extent of genetic diversity, phylogenetic relationships and the amount of gene flow among Tunisian Citrus species based on a set of 15 informative nuclear SSR molecular markers. Genotyping data highlighted an allelic richness among Tunisian Citrus species and has allowed the detection of 168 alleles among them 104.19 were effective. The partition of the total genetic diversity (HT=0.832) showed that the highest amount of variation within the Citrus species is HS=0.550, while the relative amount of the between-species genetic diversity GST does not exceed 0.338. This pattern of genetic structure was supported by low-to-moderate FST pairwise values and the presence of a gene flow (Nm) among the eight Citrus species. The lowest genetic differentiation was revealed between the species C. sinensis and C. insitorum (FST=0.111, Nm=1.99), while the highest genetic differentiation was recorded between the species C. aurantifolia and C. paradisi (FST=0.367, Nm=0.43). The established Neighbor Joining analysis showed that all genotypes were widely discriminated and clearly pooled according to their species of origin, with minor exceptions.

Phylogenetic origin of limes and lemons revealed by cytoplasmic and nuclear markers

BACKGROUND AND AIMS

The origin of limes and lemons has been a source of conflicting taxonomic opinions. Biochemical studies, numerical taxonomy and recent molecular studies suggested that cultivated Citrus species result from interspecific hybridization between four basic taxa (C. reticulata,C. maxima,C. medica and C. micrantha). However, the origin of most lemons and limes remains controversial or unknown. The aim of this study was to perform extended analyses of the diversity, genetic structure and origin of limes and lemons.

METHODS

The study was based on 133 Citrus accessions. It combined maternal phylogeny studies based on mitochondrial and chloroplastic markers, and nuclear structure analysis based on the evaluation of ploidy level and the use of 123 markers, including 73 basic taxa diagnostic single nucleotide polymorphism (SNP) and indel markers.

KEY RESULTS

The lime and lemon horticultural group appears to be highly polymorphic, with diploid, triploid and tetraploid varieties, and to result from many independent reticulation events which defined the sub-groups. Maternal phylogeny involves four cytoplasmic types out of the six encountered in the Citrus genus. All lime and lemon accessions were highly heterozygous, with interspecific admixture of two, three and even the four ancestral taxa genomes. Molecular polymorphism between varieties of the same sub-group was very low.

CONCLUSIONS

Citrus medica contributed to all limes and lemons and was the direct male parent for the main sub-groups in combination with C. micrantha or close papeda species (for C. aurata, C. excelsa, C. macrophylla and C. aurantifolia--'Mexican' lime types of Tanaka's taxa), C. reticulata(for C. limonia, C. karna and C. jambhiri varieties of Tanaka's taxa, including popular citrus rootstocks such as 'Rangpur' lime, 'Volkamer' and 'Rough' lemons), C. aurantium (for C. limetta and C. limon--yellow lemon types--varieties of Tanaka's taxa) or the C. maxima × C. reticulate hybrid (for C. limettioides--'Palestine sweet' lime types--and C. meyeri). Among triploid limes, C. latifolia accessions ('Tahiti' and 'Persian' lime types) result from the fertilization of a haploid ovule of C. limon by a diploid gamete of C. aurantifolia, while C. aurantifolia triploid accessions ('Tanepao' lime types and 'Madagascar' lemon) probably result from an interspecific backcross (a diploid ovule of C. aurantifolia fertilized by C. medica). As limes and lemons were vegetatively propagated (apomixis, horticultural practices) the intra-sub-group phenotypic diversity results from asexual variations.

Isolation and Functional Characterization of a Lycopene β-cyclase Gene Promoter from Citrus

Lycopene β-cyclases are key enzymes located at the branch point of the carotenoid biosynthesis pathway. However, the transcriptional regulatory mechanisms of LCYb1 in citrus with abundant carotenoid accumulation are still unclear. To understand the molecular basis of CsLCYb1 expression, we isolated and functionally characterized the 5' upstream sequences of CsLCYb1 from citrus. The full-length CsLCYb1 promoter and a series of its 5' deletions were fused to the β-glucuronidase (GUS) reporter gene and transferred into different plants (tomato, Arabidopsis and citrus callus) to test the promoter activities. The results of all transgenic species showed that the 1584 bp upstream region from the translational start site displayed maximal promoter activity, and the minimal promoter containing 746 bp upstream sequences was sufficient for strong basal promoter activity. Furthermore, the CsLCYb1 promoter activity was developmentally and tissue-specially regulated in transgenic Arabidopsis, and it was affected by multiple hormones and environmental cues in transgenic citrus callus under various treatments. Finer deletion analysis identified an enhancer element existing as a tandem repeat in the promoter region between -574 to -513 bp and conferring strong promoter activity. The copy numbers of the enhancer element differed among various citrus species, leading to the development of a derived simple sequence repeat marker to distinguish different species. In conclusion, this study elucidates the expression characteristics of the LCYb1 promoter from citrus and further identifies a novel enhancer element required for the promoter activity. The characterized promoter fragment would be an ideal candidate for genetic engineering and seeking of upstream trans-acting elements.

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

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 select mandarins and Ichang papeda as Citrus varieties for use in creating species devoid of these toxic compounds in future breeding programs.

Exploring the diversity in Citrus fruit colouration to decipher the relationship between plastid ultrastructure and carotenoid composition

Differentiation of new and characteristic plastid ultrastructures during ripening of citrus fruits in both peel and pulp appears to be strongly correlated with the content and complement of carotenoids. Most of the species of the Citrus genus display a wide range in fruit colouration due to differences in carotenoids; however, how this diversity is related and may contribute to plastid differentiation and ultrastructure is currently unknown. To that end, carotenoid profile and plastid ultrastructure were compared in peel and pulp of three sweet oranges: the ordinary orange-coloured Navel, rich in β,β-xanthophylls, the yellow Pinalate mutant with an elevated content of colourless carotenes and reduced β,β-xanthophylls, and the red-fleshed Cara Cara with high concentration of colourless carotenes and lycopene in the pulp; and two grapefruits: the white Marsh, with low carotenoid content, and the red Star Ruby, accumulating upstream carotenes and lycopene. The most remarkable differences in plastid ultrastructure among varieties were detected in the pulp at full colour, coinciding with major differences in carotenoid composition. Accumulation of lycopene in Cara Cara and Star Ruby pulp was associated with the presence of needle-like crystals in the plastids, while high content of upstream carotenes in Pinalate pulp was related to the development of a novel plastid type with numerous even and round vesicles. The presence of plastoglobuli was linked to phytoene and xanthophyll accumulation, suggesting these structures as the main sites for the accumulation of these pigments. Peel chromoplasts were richer in membranes compared to pulp chromoplasts, reflecting their different biogenesis. In summary, differences in carotenoid composition and accumulation of unusual carotenoids are mirrored by the development of diverse and novel chromoplast types, revealing the plasticity of these organelles to rearrange carotenoids inside different structures to allow massive accumulation and thus contributing to the chemical stability of the carotenoids.

Parentage Reconstruction in Eucalyptus nitens Using SNPs and Microsatellite Markers: A Comparative Analysis of Marker Data Power and Robustness

Pedigree reconstruction using molecular markers enables efficient management of inbreeding in open-pollinated breeding strategies, replacing expensive and time-consuming controlled pollination. This is particularly useful in preferentially outcrossed, insect pollinated Eucalypts known to suffer considerable inbreeding depression from related matings. A single nucleotide polymorphism (SNP) marker panel consisting of 106 markers was selected for pedigree reconstruction from the recently developed high-density Eucalyptus Infinium SNP chip (EuCHIP60K). The performance of this SNP panel for pedigree reconstruction in open-pollinated progenies of two Eucalyptus nitens seed orchards was compared with that of two microsatellite panels with 13 and 16 markers respectively. The SNP marker panel out-performed one of the microsatellite panels in the resolution power to reconstruct pedigrees and out-performed both panels with respect to data quality. Parentage of all but one offspring in each clonal seed orchard was correctly matched to the expected seed parent using the SNP marker panel, whereas parentage assignment to less than a third of the expected seed parents were supported using the 13-microsatellite panel. The 16-microsatellite panel supported all but one of the recorded seed parents, one better than the SNP panel, although there was still a considerable level of missing and inconsistent data. SNP marker data was considerably superior to microsatellite data in accuracy, reproducibility and robustness. Although microsatellites and SNPs data provide equivalent resolution for pedigree reconstruction, microsatellite analysis requires more time and experience to deal with the uncertainties of allele calling and faces challenges for data transferability across labs and over time. While microsatellite analysis will continue to be useful for some breeding tasks due to the high information content, existing infrastructure and low operating costs, the multi-species SNP resource available with the EuCHIP60k, opens a whole new array of opportunities for high-throughput, genome-wide or targeted genotyping in species of Eucalyptus.

Nuclear species-diagnostic SNP markers mined from 454 amplicon sequencing reveal admixture genomic structure of modern citrus varieties

Most cultivated Citrus species originated from interspecific hybridisation between four ancestral taxa (C. reticulata, C. maxima, C. medica, and C. micrantha) with limited further interspecific recombination due to vegetative propagation. This evolution resulted in admixture genomes with frequent interspecific heterozygosity. Moreover, a major part of the phenotypic diversity of edible citrus results from the initial differentiation between these taxa. Deciphering the phylogenomic structure of citrus germplasm is therefore essential for an efficient utilization of citrus biodiversity in breeding schemes. The objective of this work was to develop a set of species-diagnostic single nucleotide polymorphism (SNP) markers for the four Citrus ancestral taxa covering the nine chromosomes, and to use these markers to infer the phylogenomic structure of secondary species and modern cultivars. Species-diagnostic SNPs were mined from 454 amplicon sequencing of 57 gene fragments from 26 genotypes of the four basic taxa. Of the 1,053 SNPs mined from 28,507 kb sequence, 273 were found to be highly diagnostic for a single basic taxon. Species-diagnostic SNP markers (105) were used to analyse the admixture structure of varieties and rootstocks. This revealed C. maxima introgressions in most of the old and in all recent selections of mandarins, and suggested that C. reticulata × C. maxima reticulation and introgression processes were important in edible mandarin domestication. The large range of phylogenomic constitutions between C. reticulata and C. maxima revealed in mandarins, tangelos, tangors, sweet oranges, sour oranges, grapefruits, and orangelos is favourable for genetic association studies based on phylogenomic structures of the germplasm. Inferred admixture structures were in agreement with previous hypotheses regarding the origin of several secondary species and also revealed the probable origin of several acid citrus varieties. The developed species-diagnostic SNP marker set will be useful for systematic estimation of admixture structure of citrus germplasm and for diverse genetic studies.

Network analysis of postharvest senescence process in citrus fruits revealed by transcriptomic and metabolomic profiling

Citrus (Citrus spp.), a nonclimacteric fruit, is one of the most important fruit crops in global fruit industry. However, the biological behavior of citrus fruit ripening and postharvest senescence remains unclear. To better understand the senescence process of citrus fruit, we analyzed data sets from commercial microarrays, gas chromatography-mass spectrometry, and liquid chromatography-mass spectrometry and validated physiological quality detection of four main varieties in the genus Citrus. Network-based approaches of data mining and modeling were used to investigate complex molecular processes in citrus. The Citrus Metabolic Pathway Network and correlation networks were constructed to explore the modules and relationships of the functional genes/metabolites. We found that the different flesh-rind transport of nutrients and water due to the anatomic structural differences among citrus varieties might be an important factor that influences fruit senescence behavior. We then modeled and verified the citrus senescence process. As fruit rind is exposed directly to the environment, which results in energy expenditure in response to biotic and abiotic stresses, nutrients are exported from flesh to rind to maintain the activity of the whole fruit. The depletion of internal substances causes abiotic stresses, which further induces phytohormone reactions, transcription factor regulation, and a series of physiological and biochemical reactions.

Maximum-likelihood method identifies meiotic restitution mechanism from heterozygosity transmission of centromeric loci: application in citrus

Polyploidisation is a key source of diversification and speciation in plants. Most researchers consider sexual polyploidisation leading to unreduced gamete as its main origin. Unreduced gametes are useful in several crop breeding schemes. Their formation mechanism, i.e., First-Division Restitution (FDR) or Second-Division Restitution (SDR), greatly impacts the gametic and population structures and, therefore, the breeding efficiency. Previous methods to identify the underlying mechanism required the analysis of a large set of markers over large progeny. This work develops a new maximum-likelihood method to identify the unreduced gamete formation mechanism both at the population and individual levels using independent centromeric markers. Knowledge of marker-centromere distances greatly improves the statistical power of the comparison between the SDR and FDR hypotheses. Simulating data demonstrated the importance of selecting markers very close to the centromere to obtain significant conclusions at individual level. This new method was used to identify the meiotic restitution mechanism in nineteen mandarin genotypes used as female parents in triploid citrus breeding. SDR was identified for 85.3% of 543 triploid hybrids and FDR for 0.6%. No significant conclusions were obtained for 14.1% of the hybrids. At population level SDR was the predominant mechanisms for the 19 parental mandarins.

A Phylogenetic Analysis of 34 Chloroplast Genomes Elucidates the Relationships between Wild and Domestic Species within the Genus Citrus

Citrus genus includes some of the most important cultivated fruit trees worldwide. Despite being extensively studied because of its commercial relevance, the origin of cultivated citrus species and the history of its domestication still remain an open question. Here, we present a phylogenetic analysis of the chloroplast genomes of 34 citrus genotypes which constitutes the most comprehensive and detailed study to date on the evolution and variability of the genus Citrus. A statistical model was used to estimate divergence times between the major citrus groups. Additionally, a complete map of the variability across the genome of different citrus species was produced, including single nucleotide variants, heteroplasmic positions, indels (insertions and deletions), and large structural variants. The distribution of all these variants provided further independent support to the phylogeny obtained. An unexpected finding was the high level of heteroplasmy found in several of the analyzed genomes. The use of the complete chloroplast DNA not only paves the way for a better understanding of the phylogenetic relationships within the Citrus genus but also provides original insights into other elusive evolutionary processes, such as chloroplast inheritance, heteroplasmy, and gene selection.

Mapping the genetic and tissular diversity of 64 phenolic compounds in Citrus species using a UPLC-MS approach

BACKGROUND AND AIMS

Phenolic compounds contribute to food quality and have potential health benefits. Consequently, they are an important target of selection for Citrus species. Numerous studies on this subject have revealed new molecules, potential biosynthetic pathways and linkage between species. Although polyphenol profiles are correlated with gene expression, which is responsive to developmental and environmental cues, these factors are not monitored in most studies. A better understanding of the biosynthetic pathway and its regulation requires more information about environmental conditions, tissue specificity and connections between competing sub-pathways. This study proposes a rapid method, from sampling to analysis, that allows the quantitation of multiclass phenolic compounds across contrasting tissues and cultivars.

METHODS

Leaves and fruits of 11 cultivated citrus of commercial interest were collected from adult trees grown in an experimental orchard. Sixty-four phenolic compounds were simultaneously quantified by ultra-high-performance liquid chromatography coupled with mass spectrometry.

KEY RESULTS

Combining data from vegetative tissues with data from fruit tissues improved cultivar classification based on polyphenols. The analysis of metabolite distribution highlighted the massive accumulation of specific phenolic compounds in leaves and the external part of the fruit pericarp, which reflects their involvement in plant defence. The overview of the biosynthetic pathway obtained confirmed some regulatory steps, for example those catalysed by rhamnosyltransferases. The results suggest that three other steps are responsible for the different metabolite profiles in 'Clementine' and 'Star Ruby' grapefruit.

CONCLUSIONS

The method described provides a high-throughput method to study the distribution of phenolic compounds across contrasting tissues and cultivars in Citrus, and offers the opportunity to investigate their regulation and physiological roles. The method was validated in four different tissues and allowed the identification and quantitation of 64 phenolic compounds in 20 min, which represents an improvement over existing methods of analysing multiclass polyphenols.

De novo transcriptome assembly of pummelo and molecular marker development

Pummelo (Citrus grandis) is an important fruit crop worldwide because of its nutritional value. To accelerate the pummelo breeding program, it is essential to obtain extensive genetic information and develop relative molecular markers. Here, we obtained a 12-Gb transcriptome dataset of pummelo through a mixture of RNA from seven tissues using Illumina pair-end sequencing, assembled into 57,212 unigenes with an average length of 1010 bp. The annotation and classification results showed that a total of 39,584 unigenes had similar hits to the known proteins of four public databases, and 31,501 were classified into 55 Gene Ontology (GO) functional sub-categories. The search for putative molecular markers among 57,212 unigenes identified 10,276 simple sequence repeats (SSRs) and 64,720 single nucleotide polymorphisms (SNPs). High-quality primers of 1174 SSR loci were designed, of which 88.16% were localized to nine chromosomes of sweet orange. Of 100 SSR primers that were randomly selected for testing, 87 successfully amplified clear banding patterns. Of these primers, 29 with a mean PIC (polymorphic information content) value of 0.52 were effectively applied for phylogenetic analysis. Of the 20 SNP primers, 14 primers, including 54 potential SNPs, yielded target amplifications, and 46 loci were verified via Sanger sequencing. This new dataset will be a valuable resource for molecular biology studies of pummelo and provides reliable information regarding SNP and SSR marker development, thus expediting the breeding program of pummelo.

Next generation haplotyping to decipher nuclear genomic interspecific admixture in Citrus species: analysis of chromosome 2

BACKGROUND

The most economically important Citrus species originated by natural interspecific hybridization between four ancestral taxa (Citrus reticulata, Citrus maxima, Citrus medica, and Citrus micrantha) and from limited subsequent interspecific recombination as a result of apomixis and vegetative propagation. Such reticulate evolution coupled with vegetative propagation results in mosaic genomes with large chromosome fragments from the basic taxa in frequent interspecific heterozygosity. Modern breeding of these species is hampered by their complex heterozygous genomic structures that determine species phenotype and are broken by sexual hybridisation. Nevertheless, a large amount of diversity is present in the citrus gene pool, and breeding to allow inclusion of desirable traits is of paramount importance. However, the efficient mobilization of citrus biodiversity in innovative breeding schemes requires previous understanding of Citrus origins and genomic structures. Haplotyping of multiple gene fragments along the whole genome is a powerful approach to reveal the admixture genomic structure of current species and to resolve the evolutionary history of the gene pools. In this study, the efficiency of parallel sequencing with 454 methodology to decipher the hybrid structure of modern citrus species was assessed by analysis of 16 gene fragments on chromosome 2.

RESULTS

454 amplicon libraries were established using the Fluidigm array system for 48 genotypes and 16 gene fragments from chromosome 2. Haplotypes were established from the reads of each accession and phylogenetic analyses were performed using the haplotypic data for each gene fragment. The length of 454 reads and the level of differentiation between the ancestral taxa of modern citrus allowed efficient haplotype phylogenetic assignations for 12 of the 16 gene fragments. The analysis of the mixed genomic structure of modern species and cultivars (i) revealed C. maxima introgressions in modern mandarins, (ii) was consistent with previous hypotheses regarding the origin of secondary species, and (iii) provided a new picture of the evolution of chromosome 2.

CONCLUSIONS

454 sequencing was an efficient strategy to establish haplotypes with significant phylogenetic assignations in Citrus, providing a new picture of the mixed structure on chromosome 2 in 48 citrus genotypes.

Utility of indels for species-level identification of a biologically complex plant group: a study with intergenic spacer in Citrus

The Consortium of Barcode of Life plant working group proposed to use the defined portion of plastid genes rbcL and matK either singly or in combination as the standard DNA barcode for plants. But DNA barcode based identification of biologically complex plant groups are always a challenging task due to the occurrence of natural hybridization. Here, we examined the use of indels polymorphism in trnH-psbA and trnL-trnF sequences for rapid species identification of citrus. DNA from young leaves of selected citrus species were isolated and matK gene (~800 bp) and trnH-psbA spacer (~450 bp) of Chloroplast DNA was amplified for species level identification. The sequences within the group taxa of Citrus were aligned using the ClustalX program. With few obvious misalignments were corrected manually using the similarity criterion. We identified a 54 bp inverted repeat or palindrome sequence (27-80 regions) and 6 multi residues indel coding regions. Large inverted repeats in cpDNA provided authentication at the higher taxonomic levels. These diagnostics indel marker from trnH-psbA were successful in identifying different species (5 out of 7) within the studied Citrus except Citrus limon and Citrus medica. These two closely related species are distinguished through the 6 bp deletion in trnL-trnF. This study demonstrated that the indel polymorphism based approach easily characterizes the Citrus species and the same may be applied in other complex groups. Likewise other indels occurring intergenic spacer of chloroplast regions may be tested for rapid identification of other secondary citrus species.

Sequencing of diverse mandarin, pummelo and orange genomes reveals complex history of admixture during citrus domestication

Cultivated citrus are selections from, or hybrids of, wild progenitor species whose identities and contributions to citrus domestication remain controversial. Here we sequence and compare citrus genomes--a high-quality reference haploid clementine genome and mandarin, pummelo, sweet-orange and sour-orange genomes--and show that cultivated types derive from two progenitor species. Although cultivated pummelos represent selections from one progenitor species, Citrus maxima, cultivated mandarins are introgressions of C. maxima into the ancestral mandarin species Citrus reticulata. The most widely cultivated citrus, sweet orange, is the offspring of previously admixed individuals, but sour orange is an F1 hybrid of pure C. maxima and C. reticulata parents, thus implying that wild mandarins were part of the early breeding germplasm. A Chinese wild 'mandarin' diverges substantially from C. reticulata, thus suggesting the possibility of other unrecognized wild citrus species. Understanding citrus phylogeny through genome analysis clarifies taxonomic relationships and facilitates sequence-directed genetic improvement.

Genetic diversity and population structure of the Tibetan poplar (Populus szechuanica var. tibetica) along an altitude gradient

Background

The Tibetan poplar (Populus szechuanica var. tibetica Schneid), which is distributed at altitudes of 2,000-4,500 m above sea level, is an ecologically important species of the Qinghai-Tibet Plateau and adjacent areas. However, the genetic adaptations responsible for its ability to cope with the harsh environment remain unknown.

Results

In this study, a total of 24 expressed sequence tag microsatellite (EST-SSR) markers were used to evaluate the genetic diversity and population structure of Tibetan poplars along an altitude gradient. The 172 individuals were of genotypes from low-, medium- and high-altitude populations, and 126 alleles were identified. The expected heterozygosity (H_E) value ranged from 0.475 to 0.488 with the highest value found in low-altitude populations and the lowest in high-altitude populations. Genetic variation was low among populations, indicating a limited influence of altitude on microsatellite variation. Low genetic differentiation and high levels of gene flow were detected both between and within the populations along the altitude gradient. An analysis of molecular variance (AMOVA) showed that 6.38% of the total molecular variance was attributed to diversity between populations, while 93.62% variance was associated with differences within populations. There was no clear correlation between genetic variation and altitude, and a Mantel test between genetic distance and altitude resulted in a coefficient of association of r = 0.001, indicating virtually no correlation.

Conclusion

Microsatellite genotyping results showing genetic diversity and low differentiation suggest that extensive gene flow may have counteracted local adaptations imposed by differences in altitude. The genetic analyses carried out in this study provide new insight for conservation and optimization of future arboriculture.