Somatic variations led to the selection of acidic and acidless orange cultivars

Global characterization of somatic mutations and DNA methylation changes during vegetative propagation in strawberries

Somatic mutations arise and accumulate during tissue culture and vegetative propagation, potentially affecting various traits in horticultural crops, but their characteristics are still unclear. Here, somatic mutations in regenerated woodland strawberry derived from tissue culture of shoot tips under different conditions and 12 cultivated strawberry individuals are analyzed by whole genome sequencing. The mutation frequency of single nucleotide variants is significantly increased with increased hormone levels or prolonged culture time in the range of 3.3 × 10-8-3.0 × 10-6 mutations per site. CG methylation shows a stable reduction (0.71%-8.03%) in regenerated plants, and hypoCG-DMRs are more heritable after sexual reproduction. A high-quality haplotype-resolved genome is assembled for the strawberry cultivar "Beni hoppe." The 12 "Beni hoppe" individuals randomly selected from different locations show 4731-6005 mutations relative to the reference genome, and the mutation frequency varies among the subgenomes. Our study has systematically characterized the genetic and epigenetic variants in regenerated woodland strawberry plants and different individuals of the same strawberry cultivar, providing an accurate assessment of somatic mutations at the genomic scale and nucleotide resolution in plants.

Diversity and Evolution of NLR Genes in Citrus Species

NLR genes are crucial components of the effector-triggered immunity (ETI) system, responsible for recognizing pathogens and initiating immune responses. Although NLR genes in many plant species have been extensively studied, the diversity of NLR genes in citrus remains largely unknown. Our analysis revealed significant variations in the copy numbers of NLR genes among these species. Gene duplication and recombination were identified as the major driving forces behind this diversity. Additionally, horizontal gene transfer (HGT) emerged as the principal mechanism responsible for the increase in NLR gene copy number in A. buxifolia. The citrus NLR genes were classified into four categories: TIR-NBS-LRR (TNL), CC-NBS-LRR (CNL), RPW8-NBS-LRR (RNL), and NL. Our findings indicate that TNL, RNL, and CNL genes originated from NL genes through the acquisition of TIR and RPW8 domains, along with CC motifs, followed by the random loss of corresponding domains. Phylogenetic analysis suggested that citrus NLR genes originated alongside the species and underwent adaptive evolution, potentially playing crucial roles in the global colonization of citrus. This study provides important insights into the diversity of citrus NLR genes and serves as a foundational dataset for future research aimed at breeding disease-resistant citrus varieties.

Dynamic Analysis of the Fruit Sugar-Acid Profile in a Fresh-Sweet Mutant and Wild Type in 'Shatangju' (Citrus reticulata cv.)

Citrate is a major determinant of fruit flavor quality. Currently, citrus species and/or varieties with significant alterations in citrate level have greatly advanced the molecular basis of citrate accumulation in fruit. However, in-depth dissections of the molecular mechanism specific to citrate accumulation are still limited due to the lack of mutants, especially within one single variety. In this study, a fresh-sweet 'Shatangju' mutant (Citrus reticulata cv.) was obtained during a survey of citrus resources in Guangdong, China, and the phenotype, fruit morphology, and primary flavor profiles were comparatively analyzed. Unlike the wild-type 'Shatangju' (WT), the mutant (MT) material exhibited a dwarfed and multi-branched tree shape, delayed flowering and fruit ripening at maturity, a prolonged fruit tree-retention time, and a decreased single fruit weight at maturity. Dynamic measurement of the metabolite levels further suggested that the contents and fluctuation patterns of vitamin C, malate, quinate, and oxalate showed no obvious difference between MT and MT fruits, while the citrate level in MT fruits significantly decreased over various developmental stages, ranging from 0.356 to 1.91 mg g-1 FW. In addition, the accumulation patterns of the major soluble sugars (sucrose, fructose, and glucose), as well as the sugar/acid ratio, were also altered in MT fruits during development. Taken together, this study provides a novel acid-free 'Shatangju' mutant, which can serve as a powerful tool for the research of fruit flavor quality, especially for the comprehensive understanding of the molecular mechanism of citrate accumulation in fruits.

The genomes of Australian wild limes

Australian wild limes occur in highly diverse range of environments and are a unique genetic resource within the genus Citrus. Here we compare the haplotype-resolved genome assemblies of six Australian native limes, including four new assemblies generated using PacBio HiFi and Hi-C sequencing data. The size of the genomes was between 315 and 391 Mb with contig N50s from 29.5 to 35 Mb. Gene completeness of the assemblies was estimated to be from 98.4 to 99.3% and the annotations from 97.7 to 98.9% based upon BUSCO, confirming the high contiguity and completeness of the assembled genomes. High collinearity was observed among the genomes and the two haplotype assemblies for each species. Gene duplication and evolutionary analysis demonstrated that the Australian citrus have undergone only one ancient whole-genome triplication event during evolution. The highest number of species-specific and expanded gene families were found in C. glauca and they were primarily enriched in purine, thiamine metabolism, amino acids and aromatic amino acids metabolism which might help C. glauca to mitigate drought, salinity, and pathogen attacks in the drier environments in which this species is found. Unique genes related to terpene biosynthesis, glutathione metabolism, and toll-like receptors in C. australasica, and starch and sucrose metabolism genes in both C. australis and C. australasica might be important candidate genes for HLB tolerance in these species. Expanded gene families were not lineage specific, however, a greater number of genes related to plant-pathogen interactions, predominantly disease resistant protein, was found in C. australasica and C. australis.

Comparative Analysis of Transposable Elements in the Genomes of Citrus and Citrus-Related Genera

Transposable elements (TEs) significantly contribute to the evolution and diversity of plant genomes. In this study, we explored the roles of TEs in the genomes of Citrus and Citrus-related genera by constructing a pan-genome TE library from 20 published genomes of Citrus and Citrus-related accessions. Our results revealed an increase in TE content and the number of TE types compared to the original annotations, as well as a decrease in the content of unclassified TEs. The average length of TEs per assembly was approximately 194.23 Mb, representing 41.76% (Murraya paniculata) to 64.76% (Citrus gilletiana) of the genomes, with a mean value of 56.95%. A significant positive correlation was found between genome size and both the number of TE types and TE content. Consistent with the difference in mean whole-genome size (39.83 Mb) between Citrus and Citrus-related genera, Citrus genomes contained an average of 34.36 Mb more TE sequences than Citrus-related genomes. Analysis of the estimated insertion time and half-life of long terminal repeat retrotransposons (LTR-RTs) suggested that TE removal was not the primary factor contributing to the differences among genomes. These findings collectively indicate that TEs are the primary determinants of genome size and play a major role in shaping genome structures. Principal coordinate analysis (PCoA) of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) identifiers revealed that the fragmented TEs were predominantly derived from ancestral genomes, while intact TEs were crucial in the recent evolutionary diversification of Citrus. Moreover, the presence or absence of intact TEs near the AdhE superfamily was closely associated with the bitterness trait in the Citrus species. Overall, this study enhances TE annotation in Citrus and Citrus-related genomes and provides valuable data for future genetic breeding and agronomic trait research in Citrus.

Cas12a RNP-mediated co-transformation enables transgene-free multiplex genome editing, long deletions, and inversions in citrus chromosome

Introduction

Citrus canker, caused by Xanthomonas citri subsp. citri (Xcc), is a devastating disease worldwide. Previously, we successfully generated canker-resistant Citrus sinensis cv. Hamlin lines in the T0 generation. This was achieved through the transformation of embryogenic protoplasts using the ribonucleoprotein (RNP) containing Cas12a and one crRNA to edit the canker susceptibility gene, CsLOB1, which led to small indels.

Methods

Here, we transformed embryogenic protoplasts of Hamlin with RNP containing Cas12a and three crRNAs.

Results

Among the 10 transgene-free genome-edited lines, long deletions were obtained in five lines. Additionally, inversions were observed in three of the five edited lines with long deletions, but not in any edited lines with short indel mutations, suggesting long deletions maybe required for inversions. Biallelic mutations were observed for each of the three target sites in four of the 10 edited lines when three crRNAs were used, demonstrating that transformation of embryogenic citrus protoplasts with Cas12a and three crRNAs RNP can be very efficient for multiplex editing. Our analysis revealed the absence of off-target mutations in the edited lines. These cslob1 mutant lines were canker- resistant and no canker symptoms were observed after inoculation with Xcc and Xcc growth was significantly reduced in the cslob1 mutant lines compared to the wild type plants.

Discussion

Taken together, RNP (Cas12a and three crRNAs) transformation of embryogenic protoplasts of citrus provides a promising solution for transgene-free multiplex genome editing with high efficiency and for deletion of long fragments.

Volatile metabolomics and transcriptomics analyses provide insights into the mechanism of volatile changes during fruit development of 'Ehime 38' (Citrus reticulata) and its bud mutant

Volatile compounds are important determinants affecting fruit flavor. Previous study has identified a bud mutant of 'Ehime 38' (Citrus reticulata) with different volatile profile. However, the volatile changes between WT and MT during fruit development and underlying mechanism remain elusive. In this study, a total of 35 volatile compounds were identified in the pulps of WT and MT at five developmental stages. Both varieties accumulated similar and the highest levels of volatiles at stage S1, and showed a downward trend as the fruit develops. However, the total volatile contents in the pulps of MT were 1.4-2.5 folds higher than those in WT at stages S2-S5, which was mainly due to the increase in the content of d-limonene. Transcriptomic and RT-qPCR analysis revealed that most genes in MEP pathway were positively correlated with the volatile contents, of which DXS1 might mainly contribute to the elevated volatiles accumulation in MT by increasing the flux into the MEP pathway. Moreover, temporal expression analysis indicated that these MEP pathway genes functioned at different developmental stages. This study provided comprehensive volatile metabolomics and transcriptomics characterizations of a citrus mutant during fruit development, which is valuable for fruit flavor improvement in citrus.

Cytochrome P450 CitCYP97B modulates carotenoid accumulation diversity by hydroxylating β-cryptoxanthin in Citrus

Carotenoids in plant foods provide health benefits by functioning as provitamin A. One of the vital provitamin A carotenoids, β-cryptoxanthin, is typically plentiful in citrus fruit. However, little is known about the genetic basis of β-cryptoxanthin accumulation in citrus. Here, we performed a widely targeted metabolomic analysis of 65 major carotenoids and carotenoid derivatives to characterize carotenoid accumulation in Citrus and determine the taxonomic profile of β-cryptoxanthin. We used data from 81 newly sequenced representative accessions and 69 previously sequenced Citrus cultivars to reveal the genetic basis of β-cryptoxanthin accumulation through a genome-wide association study. We identified a causal gene, CitCYP97B, which encodes a cytochrome P450 protein whose substrate and metabolic pathways in land plants were undetermined. We subsequently demonstrated that CitCYP97B functions as a novel monooxygenase that specifically hydroxylates the β-ring of β-cryptoxanthin in a heterologous expression system. In planta experiments provided further evidence that CitCYP97B negatively regulates β-cryptoxanthin content. Using the sequenced Citrus accessions, we found that two critical structural cis-element variations contribute to increased expression of CitCYP97B, thereby altering β-cryptoxanthin accumulation in fruit. Hybridization/introgression appear to have contributed to the prevalence of two cis-element variations in different Citrus types during citrus evolution. Overall, these findings extend our understanding of the regulation and diversity of carotenoid metabolism in fruit crops and provide a genetic target for production of β-cryptoxanthin-biofortified products.

Chromosome-Scale, De Novo, Phased Genome Assemblies of Three Australian Limes: Citrus australasica, C. inodora, and C. glauca

Huanglongbing (HLB) is a severe citrus disease worldwide. Wild Australian limes like Citrus australasica, C. inodora, and C. glauca possess beneficial HLB resistance traits. Individual trees of the three taxa were extensively used in a breeding program for over a decade to introgress resistance traits into commercial-quality citrus germplasm. We generated high-quality, phased, de novo genome assemblies of the three Australian limes using PacBio long-read sequencing. The genome assembly sizes of the primary and alternate haplotypes were determined for C. australasica (337 Mb/335 Mb), C. inodora (304 Mb/299 Mb), and C. glauca (376 Mb/379 Mb). The nine chromosome-scale scaffolds included 86-91% of the genome sequences generated. The integrity and completeness of the assembled genomes were estimated to be at 97.2-98.8%. Gene annotation studies identified 25,461 genes in C. australasica, 27,665 in C. inodora, and 30,067 in C. glauca. Genes belonging to 118 orthogroups were specific to Australian lime genomes compared to other citrus genomes analyzed. Significantly fewer canonical resistance (R) genes were found in C. inodora and C. glauca (319 and 449, respectively) compared to C. australasica (576), C. clementina (579), and C. sinensis (651). Similar patterns were observed for other gene families associated with potential HLB resistance, including Phloem protein 2 (PP2) and Callose synthase (CalS) genes predicted in the Australian lime genomes. The genomic information on Australian limes developed in the present study will help understand the genetic basis of HLB resistance.

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.

Seasonal drought promotes citrate accumulation in citrus fruit through the CsABF3-activated CsAN1-CsPH8 pathway

Plenty of rainfall but unevenly seasonal distribution happens regularly in southern China. Seasonal drought from summer to early autumn leads to citrus fruit acidification, but how seasonal drought regulates citrate accumulation remains unknown. Herein, we employed a set of physiological, biochemical, and molecular approaches to reveal that CsABF3 responds to seasonal drought stress and modulates citrate accumulation in citrus fruits by directly regulating CsAN1 and CsPH8. Here, we demonstrated that irreversible acidification of citrus fruits is caused by drought lasting for > 30 d during the fruit enlargement stage. We investigated the transcriptome characteristics of fruits affected by drought and corroborated the pivotal roles of a bHLH transcription factor (CsAN1) and a P3A-ATPase gene (CsPH8) in regulating citrate accumulation in response to drought. Abscisic acid (ABA)-responsive element binding factor 3 (CsABF3) was upregulated by drought in an ABA-dependent manner. CsABF3 activated CsAN1 and CsPH8 expression by directly and specifically binding to the ABA-responsive elements (ABREs) in the promoters and positively regulated citrate accumulation. Taken together, this study sheds new light on the regulatory module ABA-CsABF3-CsAN1-CsPH8 responsible for citrate accumulation under drought stress, which advances our understanding of quality formation of citrus fruit.

Multi-omics analysis provides new insights into the changes of important nutrients and fructose metabolism in loquat bud sport mutant

Bud sport is a common and stable somatic variation in perennial fruit trees, and often leads to significant modification of fruit traits and affects the breeding value. To investigate the impact of bud sport on the main metabolites in the fruit of white-fleshed loquat, we conducted a multi-omics analysis of loquat fruits at different developmental stages of a white-fleshed bud sport mutant of Dongting loquat (TBW) and its wild type (TBY). The findings from the detection of main fruit quality indices and metabolites suggested that bud sport resulted in a reduction in the accumulation of carotenoids, fructose, titratable acid and terpenoids at the mature stage of TBW, while leading to the accumulation of flavonoids, phenolic acids, amino acids and lipids. The comparably low content of titratable acid further enhances the balanced and pleasent taste profile of TBW. Expression patterns of differentially expressed genes involved in fructose metabolism exhibited a significant increase in the expression level of S6PDH (EVM0006243, EVM0044405) prior to fruit maturation. The comparison of protein sequences and promoter region of S6PDH between TBY and TBW revealed no structural variations that would impact gene function or expression, indicating that transcription factors may be responsible for the rapid up-regulation of S6PDH before maturation. Furthermore, correlation analysis helped to construct a comprehensive regulatory network of fructose metabolism in loquat, including 23 transcription factors, six structural genes, and nine saccharides. Based on the regulatory network and existing studies, it could be inferred that transcription factors such as ERF, NAC, MYB, GRAS, and bZIP may promote fructose accumulation in loquat flesh by positively regulating S6PDH. These findings improve our understanding of the nutritional value and breeding potential of white-fleshed loquat bud sport mutant, as well as serve as a foundation for exploring the genes and transcription factors that regulate fructose metabolism in loquat.

Generation of transgene-free canker-resistant Citrus sinensis cv. Hamlin in the T0 generation through Cas12a/CBE co-editing

Citrus canker disease affects citrus production. This disease is caused by Xanthomonas citri subsp. citri (Xcc). Previous studies confirmed that during Xcc infection, PthA4, a transcriptional activator like effector (TALE), is translocated from the pathogen to host plant cells. PthA4 binds to the effector binding elements (EBEs) in the promoter region of canker susceptibility gene LOB1 (EBEPthA4-LOBP) to activate its expression and subsequently cause canker symptoms. Previously, the Cas12a/CBE co-editing method was employed to disrupt EBEPthA4-LOBP of pummelo, which is highly homozygous. However, most commercial citrus cultivars are heterozygous hybrids and more difficult to generate homozygous/biallelic mutants. Here, we employed Cas12a/CBE co-editing method to edit EBEPthA4-LOBP of Hamlin (Citrus sinensis), a commercial heterozygous hybrid citrus cultivar grown worldwide. Binary vector GFP-p1380N-ttLbCas12a:LOBP1-mPBE:ALS2:ALS1 was constructed and shown to be functional via Xcc-facilitated agroinfiltration in Hamlin leaves. This construct allows the selection of transgene-free regenerants via GFP, edits ALS to generate chlorsulfuron-resistant regenerants as a selection marker for genome editing resulting from transient expression of the T-DNA via nCas9-mPBE:ALS2:ALS1, and edits gene(s) of interest (i.e., EBEPthA4-LOBP in this study) through ttLbCas12a, thus creating transgene-free citrus. Totally, 77 plantlets were produced. Among them, 8 plantlets were transgenic plants (#HamGFP1 - #HamGFP8), 4 plantlets were transgene-free (#HamNoGFP1 - #HamNoGFP4), and the rest were wild type. Among 4 transgene-free plantlets, three lines (#HamNoGFP1, #HamNoGFP2 and #HamNoGFP3) contained biallelic mutations in EBEpthA4, and one line (#HamNoGFP4) had homozygous mutations in EBEpthA4. We achieved 5.2% transgene-free homozygous/biallelic mutation efficiency for EBEPthA4-LOBP in C. sinensis cv. Hamlin, compared to 1.9% mutation efficiency for pummelo in a previous study. Importantly, the four transgene-free plantlets and 3 transgenic plantlets that survived were resistant against citrus canker. Taken together, Cas12a/CBE co-editing method has been successfully used to generate transgene-free canker-resistant C. sinensis cv. Hamlin in the T0 generation via biallelic/homozygous editing of EBEpthA4 of the canker susceptibility gene LOB1.

A tetraploid-dominated cytochimera developed from a natural bud mutant of the nonapomictic mandarin variety 'Orah'

Nonapomictic citrus tetraploids are desirable in citrus breeding for the production of triploid, seedless varieties, and polyploid rootstocks. However, only a few lines have been reported, and they were all generated using chemical methods. A 2x + 4 × cytochimera of the nonapomictic citrus variety 'Orah' mandarin, which developed from a bud mutant, was found due to its morphology differing from that of diploid plants and characterised via ploidy analysis combining flow cytometry and chromosome observation. The chimaera was stable, and there were 1.86-1.90 times as tetraploid cells as diploid cells. Anatomical structure observation revealed that the 'Orah' chimaera may be a periclinal chimaera with diploid cells in the L1 layer and tetraploid cells in the L2 and L3 layers. The chimaera showed some typical traits of polyploid plants, including thicker shoots, wider and thicker leaves, larger flowers and fruits, and fewer but larger seeds in fruits than in diploid plants. Almost all the seeds of the chimaera were monoembryonic. Most of the self-pollinated progenies of the chimaera were identified as tetraploids, and some triploid, pentaploid, and hexaploid plants were found. As a female, the chimaera produced allotriploids when crossed with Australian finger lime. In addition, 6 plants developed from polyembryonic seeds of the chimaera were identified as sexual tetraploid progenies with low-level recombinant genomes. Therefore, the 'Orah' 2x + 4 × chimaera can be used as a female parent to produce hybrid triploid and tetraploid citrus plants with high efficiency. Identification of the chimaera demonstrated that tetraploid citrus plants, especially nonapomictic varieties, can be generated from shoot bud mutants.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01456-x.

Chemical and genetic basis of orange flavor

Sweet orange (Citrus sinensis) exhibits limited genetic diversity and high susceptibility to Huanglongbing (HLB). Breeding HLB-tolerant orange-like hybrids is in dire need. However, our understanding of the key compounds responsible for orange flavor and their genetic regulation remains elusive. Evaluating 179 juice samples, including oranges, mandarins, Poncirus trifoliata, and hybrids, distinct volatile compositions were found. A random forest model predicted untrained samples with 78% accuracy and identified 26 compounds crucial for orange flavor. Notably, seven esters differentiated orange from mandarin flavor. Cluster analysis showed six esters with shared genetic control. Differential gene expression analysis identified C. sinensis alcohol acyltransferase 1 (CsAAT1) responsible for ester production in orange. Its activity was validated through overexpression assays. Phylogeny revealed the functional allele was inherited from pummelo. A SNP-based DNA marker in the coding region accurately predicted phenotypes. This study enhances our understanding of orange flavor compounds and their biosynthetic pathways and expands breeding options for orange-like cultivars.

Development and Assessment of SNP Genotyping Arrays for Citrus and Its Close Relatives

Rapid advancements in technologies provide various tools to analyze fruit crop genomes to better understand genetic diversity and relationships and aid in breeding. Genome-wide single nucleotide polymorphism (SNP) genotyping arrays offer highly multiplexed assays at a relatively low cost per data point. We report the development and validation of 1.4M SNP Axiom® Citrus HD Genotyping Array (Citrus 15AX 1 and Citrus 15AX 2) and 58K SNP Axiom® Citrus Genotyping Arrays for Citrus and close relatives. SNPs represented were chosen from a citrus variant discovery panel consisting of 41 diverse whole-genome re-sequenced accessions of Citrus and close relatives, including eight progenitor citrus species. SNPs chosen mainly target putative genic regions of the genome and are accurately called in both Citrus and its closely related genera while providing good coverage of the nuclear and chloroplast genomes. Reproducibility of the arrays was nearly 100%, with a large majority of the SNPs classified as the most stringent class of markers, "PolyHighResolution" (PHR) polymorphisms. Concordance between SNP calls in sequence data and array data average 98%. Phylogenies generated with array data were similar to those with comparable sequence data and little affected by 3 to 5% genotyping error. Both arrays are publicly available.

Phased genomics reveals hidden somatic mutations and provides insight into fruit development in sweet orange

Although revisiting the discoveries and implications of genetic variations using phased genomics is critical, such efforts are still lacking. Somatic mutations represent a crucial source of genetic diversity for breeding and are especially remarkable in heterozygous perennial and asexual crops. In this study, we focused on a diploid sweet orange (Citrus sinensis) and constructed a haplotype-resolved genome using high fidelity (HiFi) reads, which revealed 10.6% new sequences. Based on the phased genome, we elucidate significant genetic admixtures and haplotype differences. We developed a somatic detection strategy that reveals hidden somatic mutations overlooked in a single reference genome. We generated a phased somatic variation map by combining high-depth whole-genome sequencing (WGS) data from 87 sweet orange somatic varieties. Notably, we found twice as many somatic mutations relative to a single reference genome. Using these hidden somatic mutations, we separated sweet oranges into seven major clades and provide insight into unprecedented genetic mosaicism and strong positive selection. Furthermore, these phased genomics data indicate that genomic heterozygous variations contribute to allele-specific expression during fruit development. By integrating allelic expression differences and somatic mutations, we identified a somatic mutation that induces increases in fruit size. Applications of phased genomics will lead to powerful approaches for discovering genetic variations and uncovering their effects in highly heterozygous plants. Our data provide insight into the hidden somatic mutation landscape in the sweet orange genome, which will facilitate citrus breeding.

Vacuolar proteomic analysis reveals tonoplast transporters for accumulation of citric acid and sugar in citrus fruit

Vacuole largely dictates the fruit taste and flavor, as most of the sugars and organic acids are stored in the vacuoles of the fruit. However, difficulties associated with vacuole separation severely hinder identification and characterization of vacuolar proteins in fruit species. In this study, we established an effective approach for separating vacuoles and successfully purified vacuolar protein from six types of citrus fruit with varying patterns of sugar and organic acid contents. By using label-free LC-MS/MS proteomic analysis, 1443 core proteins were found to be associated with the essential functions of vacuole in citrus fruit. Correlation analysis of metabolite concentration with proteomic data revealed a transporter system for the accumulation of organic acid and soluble sugars in citrus. Furthermore, we characterized the physiological roles of selected key tonoplast transporters, ABCG15, Dict2.1, TMT2, and STP7 in the accumulation of citric acid and sugars. These findings provide a novel perspective and practical solution for investigating the transporters underlying the formation of citrus taste and flavor.

Soluble sugar and organic acid composition and flavor evaluation of Chinese cherry fruits

Chinese cherry is an economically important fruit crop native to China. Flavor quality is greatly influenced by compositions of soluble sugars and organic acids. To better understand the flavor quality of Chinese cherry, we determined sugar and acid components in thirty-eight landrace and cultivar collections, and two wild resources using the HPLC method. Glucose and fructose were the main components, accounting for 85.91% of soluble sugars. Malic acid was the predominant organic acid, with an average proportion of 65.73% of total acids. Correlation and PCA analysis revealed seven key indicators for evaluating fruit flavor. Compared with wild Chinese cherry, the cultivated collections exhibited higher levels of soluble sugars, especially fructose, and lower levels of organic acid, particularly malic acid in fruits. Finally, we have established grading criteria for seven flavor indicators in Chinese cherry. Our study provides valuable references for identifying flavor compounds and improving flavor quality of Chinese cherry.

GIFTdb: a useful gene database for plant fruit traits improving

Fruit traits are critical determinants of plant fitness, resource diversity, productive and quality. Gene regulatory networks in plants play an essential role in determining fruit traits, such as fruit size, yield, firmness, aroma and other important features. Many research studies have focused on elucidating the associated signaling pathways and gene interaction mechanism to better utilize gene resources for regulating fruit traits. However, the availability of specific database of genes related to fruit traits for use by the plant research community remains limited. To address this limitation, we developed the Gene Improvements for Fruit Trait Database (GIFTdb, http://giftdb.agroda.cn). GIFTdb contains 35 365 genes, including 896 derived from the FR database 1.0, 305 derived from 30 882 articles from 2014 to 2021, 236 derived from the Universal Protein Resource (UniProt) database, and 33 928 identified through homology analysis. The database supports several aided analysis tools, including signal transduction pathways, gene ontology terms, protein-protein interactions, DNAWorks, Basic Local Alignment Search Tool (BLAST), and Protein Subcellular Localization Prediction (WoLF PSORT). To provide information about genes currently unsupported in GIFTdb, potential fruit trait-related genes can be searched based on homology with the supported genes. GIFTdb can provide valuable assistance in determining the function of fruit trait-related genes, such as MYB306-like, by conducting a straightforward search. We believe that GIFTdb will be a valuable resource for researchers working on gene function annotation and molecular breeding to improve fruit traits.

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.

Unlikely heroes on the long and winding road to potato inbreeding

Conversion of potato from a tetraploid, heterozygous, vegetatively propagated crop to a diploid F1 hybrid, propagated via botanical seed, would constitute a considerable advance for global agriculture, but faces multiple challenges. One such challenge is the difficulty in inbreeding potato, which involves purging deleterious alleles from its genome. This commentary discusses possible reasons for this difficulty and highlights a recent sequence-based effort to classify SNP variation, in potato germplasm, according to its deleterious potential. Tools and strategies connected to this database may facilitate development of F1 hybrids.

Chromosome-level genome assembly and population genomics of Robinia pseudoacacia reveal the genetic basis for its wide cultivation

Urban greening provides important ecosystem services and ideal places for urban recreation and is a serious consideration for municipal decision-makers. Among the tree species cultivated in urban green spaces, Robinia pseudoacacia stands out due to its attractive flowers, fragrances, high trunks, wide adaptability, and essential ecosystem services. However, the genomic basis and consequences of its wide-planting in urban green spaces remains unknown. Here, we report the chromosome-level genome assembly of R. pseudoacacia, revealing a genome size of 682.4 Mb and 33,187 protein-coding genes. More than 99.3% of the assembly is anchored to 11 chromosomes with an N50 of 59.9 Mb. Comparative genomic analyses among 17 species reveal that gene families related to traits favoured by urbanites, such as wood formation, biosynthesis, and drought tolerance, are notably expanded in R. pseudoacacia. Our population genomic analyses further recover 11 genes that are under recent selection. Ultimately, these genes play important roles in the biological processes related to flower development, water retention, and immunization. Altogether, our results reveal the evolutionary forces that shape R. pseudoacacia cultivated for urban greening. These findings also present a valuable foundation for the future development of agronomic traits and molecular breeding strategies for R. pseudoacacia.

Genome-wide characterization of long terminal repeat retrotransposons provides insights into trait evolution of four cucurbit species

Cucurbits are a diverse plant family that includes economically important crops, such as cucumber, watermelon, melon, and pumpkin. Knowledge of the roles that long terminal repeat retrotransposons (LTR-RTs) have played in diversification of cucurbit species is limited; to add to understanding of the roles of LTR-RTs, we assessed their distributions in four cucurbit species. We identified 381, 578, 1086, and 623 intact LTR-RTs in cucumber (Cucumis sativus L. var. sativus cv. Chinese Long), watermelon (Citrullus lanatus subsp. vulgaris cv. 97103), melon (Cucumis melo cv. DHL92), and Cucurbita (Cucurbita moschata var. Rifu), respectively. Among these LTR-RTs, the Ale clade of the Copia superfamily was the most abundant in all the four cucurbit species. Insertion time and copy number analysis revealed that an LTR-RT burst occurred approximately 2 million years ago in cucumber, watermelon, melon, and Cucurbita, and may have contributed to their genome size variation. Phylogenetic and nucleotide polymorphism analyses suggested that most LTR-RTs were formed after species diversification. Analysis of gene insertions by LTR-RTs revealed that the most frequent insertions were of Ale and Tekay and that genes related to dietary fiber synthesis were the most commonly affected by LTR-RTs in Cucurbita. These results increase our understanding of LTR-RTs and their roles in genome evolution and trait characterization in cucurbits.

Chromosome-level assemblies of cultivated water chestnut Trapa bicornis and its wild relative Trapa incisa

Water chestnut (Trapa L.) is a floating-leaved aquatic plant with high edible and medicinal value. In this study, we presented chromosome-level genome assemblies of cultivated large-seed species Trapa bicornis and its wild small-seed relative Trapa incisa by using PacBio HiFi long reads and Hi-C technology. The T. bicornis and T. incisa assemblies consisted of 479.90 Mb and 463.97 Mb contigs with N50 values of 13.52 Mb and 13.77 Mb, respectively, and repeat contents of 62.88% and 62.49%, respectively. A total of 33,306 and 33,315 protein-coding genes were predicted in T. bicornis and T. incisa assemblies, respectively. There were 159,232 structural variants affecting more than 11 thousand genes detected between the two genomes. The phylogenetic analysis indicated that the lineage leading to Trapa was diverged from the lineage to Sonneratia approximately 23 million years ago. These two assemblies provide valuable resources for future evolutionary and functional genomic research and molecular breeding of water chestnut.

Anthocyanins distribution, transcriptional regulation, epigenetic and post-translational modification in fruits

Anthocyanins have indispensable functions in plant resistance, human health, and fruit coloring, which arouse people's favorite. It has been reported that anthocyanins are widely found in fruits, and can be affected by numerous factors. In this review, we systematically summarize anthocyanin functions, classifications, distributions, biosynthesis, decoration, transportation, transcriptional regulation, DNA methylation, and post-translational regulation in fruits.

Comparative transcriptome profiling of susceptible and tolerant citrus species at early and late stage of infection by "Candidatus Liberibacter asiaticus"

Citrus Huanglongbing (HLB), caused by "Candidatus Liberibacter asiaticus" (CLas), is the most destructive disease threatening global citrus industry. Most commercial cultivars were susceptible to HLB, although some showed tolerant to HLB phenotypically. Identifying tolerant citrus genotypes and understanding the mechanism correlated with tolerance to HLB is essential for breeding citrus variety tolerance/resistance to HLB. In this study, the graft assay with CLas-infected bud were performed in four citrus genotypes, including Citrus reticulata Blanco, C. sinensis, C. limon, and C. maxima. HLB tolerance was observed in C. limon and C. maxima, while C. Blanco and C. sinensis were susceptible to HLB. The time-course transcriptomic analysis revealed a significant variation in differentially expressed genes (DEGs) related to HLB between susceptible and tolerant cultivar group at early and late infection stage. Functional analysis of DEGs indicated that the activation of genes involved in SA-mediated defense response, PTI, cell wall associated immunity, endochitinase, phenylpropanoid and alpha-linolenic/linoleic lipid metabolism played an important in the tolerance of C. limon and C. maxima to HLB at early infection stage. In addition, the overactive plant defense combined with the stronger antibacterial activity (antibacterial secondary and lipid metabolism) and the suppression of pectinesterase were contributed to the long-term tolerance to HLB in C. limon and C. maxima at late infection stage. Particularly, the activation of ROS scavenging genes (catalases and ascorbate peroxidases) could help to reduce HLB symptoms in tolerant cultivars. In contrast, the overexpression of genes involved in oxidative burst and ethylene metabolism, as well as the late inducing of defense related genes could lead to the early HLB symptom development in susceptible cultivars at early infection stage. The weak defense response and antibacterial secondary metabolism, and the induce of pectinesterase were responsible for sensitivity to HLB in C. reticulata Blanco and C. sinensis at late infection stage. This study provided new insights into the tolerance/sensitivity mechanism against HLB and valuable guidance for breeding of HLB-tolerant/resistant cultivars.

Weighted gene coexpression correlation network analysis reveals the potential molecular regulatory mechanism of citrate and anthocyanin accumulation between postharvest 'Bingtangcheng' and 'Tarocco' blood orange fruit

BACKGROUND

Organic acids and anthocyanins are the most important compounds for the flavor and nutritional quality of citrus fruit. However, there are few reports on the involvement of co-regulation of citrate and anthocyanin metabolism. Here, we performed a comparative transcriptome analysis to elucidate the genes and pathways involved in both citrate and anthocyanin accumulation in postharvest citrus fruit with 'Tarocco' blood orange (TBO; high accumulation) and 'Bingtangcheng' sweet orange (BTSO; low accumulation).

RESULTS

A robust core set of 825 DEGs were found to be temporally associated with citrate and anthocyanin accumulation throughout the storage period through transcriptome analysis. Further according to the results of weighted gene coexpression correlation network analysis (WGCNA), the turquoise and brown module was highly positively correlated with both of the content of citrate and anthocyanin, and p-type ATPase (PH8), phosphoenolpyruvate carboxylase kinase (PEPCK), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3'-hydroxylase (F3'H) and glutathione S transferase (GST) were considered key structural genes. Moreover, MYB family transcription factor (PH4), Zinc finger PHD-type transcription factor (CHR4, HAC12), Zinc finger SWIM-type transcription factor (FAR1) and Zinc finger C3H1-type transcription factor (ATC3H64) were considered hub genes related to these structural genes. Further qRT-PCR analysis verified that these transcription factors were highly expressed in TBO fruit and their expression profiles were significantly positively correlated with the structural genes of citrate and anthocyanin metabolism as well as the content of citrate and anthocyanin content.

CONCLUSIONS

The findings suggest that the CHR4, FAR1, ATC3H64 and HAC12 may be the new transcription regulators participate in controlling the level of citrate and anthocyanin in postharvest TBO fruit in addition to PH4. These results may providing new insight into the regulation mechanism of citrate and anthocyanin accumulation in citrus fruit.

Somatic Mutations in Fruit Trees: Causes, Detection Methods, and Molecular Mechanisms

Somatic mutations are genetic changes that occur in non-reproductive cells. In fruit trees, such as apple, grape, orange, and peach, somatic mutations are typically observed as "bud sports" that remain stable during vegetative propagation. Bud sports exhibit various horticulturally important traits that differ from those of their parent plants. Somatic mutations are caused by internal factors, such as DNA replication error, DNA repair error, transposable elements, and deletion, and external factors, such as strong ultraviolet radiation, high temperature, and water availability. There are several methods for detecting somatic mutations, including cytogenetic analysis, and molecular techniques, such as PCR-based methods, DNA sequencing, and epigenomic profiling. Each method has its advantages and limitations, and the choice of method depends on the research question and the available resources. The purpose of this review is to provide a comprehensive understanding of the factors that cause somatic mutations, techniques used to identify them, and underlying molecular mechanisms. Furthermore, we present several case studies that demonstrate how somatic mutation research can be leveraged to discover novel genetic variations. Overall, considering the diverse academic and practical value of somatic mutations in fruit crops, especially those that require lengthy breeding efforts, related research is expected to become more active.

Genome-wide identification, characterization, and expression profile ofNBS-LRRgene family in sweet orange (Citrussinensis)

BACKGROUND

The NBS-LRR (nucleotide-binding site-leucine-rich repeat gene) gene family, known as the plant R (resistance) gene family with the most members, plays a significant role in plant resistance to various external adversity stresses. The NBS-LRR gene family has been researched in many plant species. Citrus is one of the most vital global cash crops, the number one fruit group, and the third most traded agricultural product world wild. However, as one of the largest citrus species, a comprehensive study of the NBS-LRR gene family has not been reported on sweet oranges.

METHODS

In this study, NBS-LRR genes were identified from the Citrus sinensis genome (v3.0), with a comprehensive analysis of this gene family performed, including phylogenetic analysis, gene structure, cis-acting element of a promoter, and chromosomal localization, among others. The expression pattern of NBS-LRR genes was analyzed when sweet orange fruits were infected by Penicillium digitatum, employing experimental data from our research group. It first reported the expression patterns of NBS-LRR genes under abiotic stresses, using three transcript data from NCBI (National Center for Biotechnology Information).

RESULTS

In this study, 111 NBS-LRR genes were identified in the C. sinensis genome (v3.0) and classified into seven subfamilies according to their N-terminal and C-terminal domains. The phylogenetic tree results indicate that genes containing only the NBS structural domain are more ancient in the sweet orange NBS-LRR gene family. The chromosome localization results showed that 111 NBS-LRR genes were distributed unevenly on nine chromosomes, with the most genes distributed on chromosome 1. In addition, we identified a total of 18 tandem duplication gene pairs in the sweet orange NBS-LRR gene family, and based on the Ka/Ks ratio, all of the tandem duplication genes underwent purifying selection. Transcriptome data analysis showed a significant number of NBS-LRR genes expressed under biotic and abiotic stresses, and some reached significantly different levels of expression. It indicates that the NBS-LRR gene family is vital in resistance to biotic and abiotic stresses in sweet oranges.

CONCLUSION

Our study provides the first comprehensive framework on the NBS-LRR family of genes, which provides a basis for further in-depth studies on the biological functions of NBS-LRR in growth, development, and response to abiotic stresses in sweet orange.

Citrus genomic resources unravel putative genetic determinants of Huanglongbing pathogenicity

Citrus HLB caused by Candidatus Liberibacter asiaticus is a pathogen-triggered immune disease. Here, we identified putative genetic determinants of HLB pathogenicity by integrating citrus genomic resources to characterize the pan-genome of accessions that differ in their response to HLB. Genome-wide association mapping and analysis of allele-specific expression between susceptible, tolerant, and resistant accessions further refined candidates underlying the response to HLB. We first developed a phased diploid assembly of Citrus sinensis 'Newhall' genome and produced resequencing data for 91 citrus accessions that differ in their response to HLB. These data were combined with previous resequencing data from 356 accessions for genome-wide association mapping of the HLB response. Genes determinants for HLB pathogenicity were associated with host immune response, ROS production, and antioxidants. Overall, this study has provided a significant resource of citrus genomic data and identified candidate genes to be further explored to understand the genetic determinants of HLB pathogenicity.

Molecular cytogenetic map visualizes the heterozygotic genome and identifies translocation chromosomes in Citrus sinensis

Citrus sinensis is the most cultivated and economically valuable Citrus species in the world, whose genome has been assembled by three generation sequencings. However, chromosome recognition remains a problem due to the small size of chromosomes, and difficulty in differentiating between pseudo and real chromosomes because of a highly heterozygous genome. Here, we employ fluorescence in situ hybridization (FISH) with 9 chromosome painting probes, 30 oligo pools, and 8 repetitive sequences to visualize 18 chromosomes. Then, we develop an approach to identify each chromosome in one cell through single experiment of oligo-FISH and Chromoycin A3 (CMA) staining. By this approach, we construct a high-resolution molecular cytogenetic map containing the physical positions of CMA banding and 38 sequences of FISH including centromere regions, which enable us to visualize significant differences between homologous chromosomes. Based on the map, we locate several highly repetitive sequences on chromosomes and estimate sizes and copy numbers of each site. In particular, we discover the translocation regions of chromosomes 4 and 9 in C. sinensis "Valencia." The high-resolution molecular cytogenetic map will help improve understanding of sweet orange genome assembly and also provide a fundamental reference for investigating chromosome evolution and chromosome engineering for genetic improvement in Citrus.

High level of somatic mutations detected in a diploid banana wild relative Musa basjoo

Plants are thought to lack an early segregating germline and often retain both asexual and sexual reproduction, both of which may allow somatic mutations to enter the gametes or clonal progeny, and thereby impact plant evolution. It is yet unclear how often these somatic mutations occur during plant development and what proportion is transmitted to their sexual or cloned offspring. Asexual "seedless" propagation has contributed greatly to the breeding in many fruit crops, such as citrus, grapes and bananas. Whether plants in these lineages experience substantial somatic mutation accumulation is unknown. To estimate the somatic mutation accumulation and inheritance among a clonal population of plant, here we assess somatic mutation accumulation in Musa basjoo, a diploid banana wild relative, using 30 whole-genome resequenced samples collected from five structures, including leaves, sheaths, panicle, roots and underground rhizome connecting three clonal individuals. We observed 18.5 high proportion de novo somatic mutations on average between each two adjacent clonal suckers, equivalent to ~ 2.48 × 10^-8 per site per asexual generation, higher than the per site per sexual generation rates (< 1 × 10^-8) reported in Arabidopsis and peach. Interestingly, most of these inter-ramet somatic mutations were shared simultaneously in different tissues of the same individual with a high level of variant allele fractions, suggesting that these somatic mutations arise early in ramet development and that each individual may develop only from a few apical stem cells. These results thus suggest substantial mutation accumulation in a wild relative of banana. Our work reveals the significance of somatic mutation in Musa basjoo genetics variations and contribute to the trait improvement breeding of bananas and other asexual clonal crops.

A chromosome-level phased genome enabling allele-level studies in sweet orange: a case study on citrus Huanglongbing tolerance

Sweet orange originated from the introgressive hybridizations of pummelo and mandarin resulting in a highly heterozygous genome. How alleles from the two species cooperate in shaping sweet orange phenotypes under distinct circumstances is unknown. Here, we assembled a chromosome-level phased diploid Valencia sweet orange (DVS) genome with over 99.999% base accuracy and 99.2% gene annotation BUSCO completeness. DVS enables allele-level studies for sweet orange and other hybrids between pummelo and mandarin. We first configured an allele-aware transcriptomic profiling pipeline and applied it to 740 sweet orange transcriptomes. On average, 32.5% of genes have a significantly biased allelic expression in the transcriptomes. Different cultivars, transgenic lineages, tissues, development stages, and disease status all impacted allelic expressions and resulted in diversified allelic expression patterns in sweet orange, but particularly citrus Huanglongbing (HLB) shifted the allelic expression of hundreds of genes in leaves and calyx abscission zones. In addition, we detected allelic structural mutations in an HLB-tolerant mutant (T19) and a more sensitive mutant (T78) through long-read sequencing. The irradiation-induced structural mutations mostly involved double-strand breaks, while most spontaneous structural mutations were transposon insertions. In the mutants, most genes with significant allelic expression ratio alterations (≥1.5-fold) were directly affected by those structural mutations. In T19, alleles located at a translocated segment terminal were upregulated, including CsDnaJ, CsHSP17.4B, and CsCEBPZ. Their upregulation is inferred to keep phloem protein homeostasis under the stress from HLB and enable subsequent stress responses observed in T19. DVS will advance allelic level studies in citrus.

Pan-mitogenomics reveals the genetic basis of cytonuclear conflicts in citrus hybridization, domestication, and diversification

Although interactions between the cytoplasmic and nuclear genomes occurred during diversification of many plants, the evolutionary conflicts due to cytonuclear interactions are poorly understood in crop breeding. Here, we constructed a pan-mitogenome and identified chimeric open reading frames (ORFs) generated by extensive structural variations (SVs). Meanwhile, short reads from 184 accessions of citrus species were combined to construct three variation maps for the nuclear, mitochondrial, and chloroplast genomes. The population genomic data showed discordant topologies between the cytoplasmic and nuclear genomes because of differences in mutation rates and levels of heteroplasmy from paternal leakage. An analysis of species-specific SVs indicated that mitochondrial heteroplasmy was common and that chloroplast heteroplasmy was undetectable. Interestingly, we found a prominent divergence in the mitogenomes and the highest genetic load in the, which may provide the basis for cytoplasmic male sterility (CMS) and thus influence the reshuffling of the cytoplasmic and nuclear genomes during hybridization. Using cytoplasmic replacement experiments, we identified a type of species-specific CMS in mandarin related to two chimeric mitochondrial genes. Our analyses indicate that cytoplasmic genomes from mandarin have rarely been maintained in hybrids and that paternal leakage produced very low levels of mitochondrial heteroplasmy in mandarin. A genome-wide association study (GWAS) provided evidence for three nuclear genes that encode pentatricopeptide repeat (PPR) proteins contributing to the cytonuclear interactions in the Citrus genus. Our study demonstrates the occurrence of evolutionary conflicts between cytoplasmic and nuclear genomes in citrus and has important implications for genetics and breeding.

miR156 regulates somatic embryogenesis by modulating starch accumulation in citrus

Somatic embryogenesis (SE) is a major regeneration approach for in vitro cultured tissues of plants, including citrus. However, SE capability is difficult to maintain, and recalcitrance to SE has become a major obstacle to plant biotechnology. We previously reported that miR156-SPL modules regulate SE in citrus callus. However, the downstream regulatory pathway of the miR156-SPL module in SE remains unclear. In this study, we found that transcription factors CsAGL15 and CsFUS3 bind to the CsMIR156A promoter and activate its expression. Suppression of csi-miR156a function leads to up-regulation of four target genes, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (CsSPL) genes, and reduction of SE efficiency. In the short tandem target mimic (STTM)-miR156a overexpression callus (MIM156), the number of amyloplasts and starch content were significantly reduced, and genes involved in starch synthesis and transport were down-regulated. csi-miR172d was down-regulated, whereas the target genes, CsTOE1.1 and CsTOE1.2, which inhibit the expression of starch biosynthesis genes, were up-regulated. In our working model, CsAGL15 and CsFUS3 activate csi-miR156a, which represses CsSPLs and further regulates csi-miR172d and CsTOEs, thus altering starch accumulation in callus cells and regulating SE in citrus. This study elucidates the pathway of miR156-SPLs and miR172-TOEs-mediated regulation of SE, and provides new insights into enhancing SE capability in citrus.

Structural variation and parallel evolution of apomixis in citrus during domestication and diversification

Apomixis, or asexual seed formation, is prevalent in Citrinae via a mechanism termed nucellar or adventitious embryony. Here, multiple embryos of a maternal genotype form directly from nucellar cells in the ovule and can outcompete the developing zygotic embryo as they utilize the sexually derived endosperm for growth. Whilst nucellar embryony enables the propagation of clonal plants of maternal genetic constitution, it is also a barrier to effective breeding through hybridization. To address the genetics and evolution of apomixis in Citrinae, a chromosome-level genome of the Hongkong kumquat (Fortunella hindsii) was assembled following a genome-wide variation map including structural variants (SVs) based on 234 Citrinae accessions. This map revealed that hybrid citrus cultivars shelter genome-wide deleterious mutations and SVs into heterozygous states free from recessive selection, which may explain the capability of nucellar embryony in most cultivars during Citrinae diversification. Analyses revealed that parallel evolution may explain the repeated origin of apomixis in different genera of Citrinae. Within Fortunella, we found that apomixis of some varieties originated via introgression. In apomictic Fortunella, the locus associated with apomixis contains the FhRWP gene, encoding an RWP-RK domain-containing protein previously shown to be required for nucellar embryogenesis in Citrus. We found the heterozygous SV in the FhRWP and CitRWP promoters from apomictic Citrus and Fortunella, due to either two or three miniature inverted transposon element (MITE) insertions. A transcription factor, FhARID, encoding an AT-rich interaction domain-containing protein binds to the MITEs in the promoter of apomictic varieties, which facilitates induction of nucellar embryogenesis. This study provides evolutionary genomic and molecular insights into apomixis in Citrinae and has potential ramifications for citrus breeding.

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

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

An apple somatic mutation of delayed fruit maturation date is primarily caused by a retrotransposon insertion-associated large deletion

Somatic mutations may alter important traits in tree fruits, such as fruit color, size and maturation date. Autumn Gala (AGala), a somatic mutation from apple cultivar Gala, matures 4 weeks later than Gala. To understand the mechanisms underlying the delayed maturation, RNA-seq analyses were conducted with fruit sampled at 13 (Gala) and 16 (AGala) time-points during their growth and development. Weighted gene co-expression network analysis (WGCNA) of 23 372 differentially expressed genes resulted in 25 WGCNA modules. Of these, modules 1 (r = -0.98, P = 2E-21) and 2 (r = -0.52, P = 0.004), which were suppressed in AGala, were correlated with fruit maturation date. Surprisingly, 77 of the 152 member genes in module 1 were harbored in a 2.8-Mb genomic region on chromosome 6 that was deleted and replaced by a 10.7-kb gypsy-like retrotransposon (Gy-36) from chromosome 7 in AGala. Among the 77 member genes, MdACT7 was the most suppressed (by 10.5-fold) in AGala due to a disruptive 2.5-kb insertion in coding sequence. Moreover, MdACT7 is the exclusive apple counterpart of Arabidopsis ACT7 known of essential roles in plant development, and the functional allele MdACT7, which was lost to the deletion in AGala, was associated with early fruit maturation in 268 apple accessions. Overexpressing alleles MdACT7 and Mdact7 in an Arabidopsis act7 line showed that MdACT7 largely rescued its stunted growth and delayed initial flowering while Mdact7 did not. Therefore, the 2.8-Mb hemizygous deletion is largely genetically causal for fruit maturation delay in AGala, and the total loss of MdACT7 might have contributed to the phenotype.

Multi-omics analyses reveal MdMYB10 hypermethylation being responsible for a bud sport of apple fruit color

Apple bud sports offer a rich resource for clonal selection of numerous elite cultivars. The accumulation of somatic mutations as plants develop may potentially impact the emergence of bud sports. Previous studies focused on somatic mutation in the essential genes associated with bud sports. However, the rate and function of genome-wide somatic mutations that accumulate when a bud sport arises remain unclear. In this study, we identified a branch from a 10-year-old tree of the apple cultivar 'Oregon Spur II' as a bud sport. The mutant branch showed reduced red coloration on fruit skin. Using this plant material, we assembled a high-quality haplotype reference genome consisting of 649.61 Mb sequences with a contig N50 value of 2.04 Mb. We then estimated the somatic mutation rate of the apple tree to be 4.56 × 10 -8 per base per year, and further identified 253 somatic single-nucleotide polymorphisms (SNPs), including five non-synonymous SNPs, between the original type and mutant samples. Transcriptome analyses showed that 69 differentially expressed genes between the original type and mutant fruit skin were highly correlated with anthocyanin content. DNA methylation in the promoter of five anthocyanin-associated genes was increased in the mutant compared with the original type as determined using DNA methylation profiling. Among the genetic and epigenetic factors that directly and indirectly influence anthocyanin content in the mutant apple fruit skin, the hypermethylated promoter of MdMYB10 is important. This study indicated that numerous somatic mutations accumulated at the emergence of a bud sport from a genome-wide perspective, some of which contribute to the low coloration of the bud sport.