A gap-free and haplotype-resolved lemon genome provides insights into flavor synthesis and huanglongbing (HLB) tolerance

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.

Flavonoid Profiles in the Pulp of Different Lemon Cultivars and Their Antioxidant Activity Based on UPLC-Q-TOF-MS

Previous studies have indicated that there may be differences among the varieties of lemon flavonoids, but the details have not yet been made clear, which limits the comprehensive use of different cultivated lemon varieties. In this study, ultra-performance liquid chromatography-quadrupole-time-of-flight-mass spectrometry (UPLC-Q-TOF-MS) and ultraviolet-visible spectroscopy (UV-Vis) were used to investigate the types and contents of flavonoids in the flesh of the main cultivated variety (Eureka) and five common lemon varieties, as well as their in vitro antioxidant activity. A total of 21 compounds were identified, five of which were common compounds. Among them, Verna, Lisbon, and Bearss each have characteristic components that can serve as potential criteria for variety identification. Each of the six varieties of lemon has strong antioxidant activity. The antioxidant activity of different lemon varieties is related to flavonoids. Therefore, Eureka and the other five varieties of lemon are good natural antioxidants, and the cultivation and industrial production of lemons should consider the needs and selection of suitable varieties.

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.

An integrated multi-omics approach reveals polymethoxylated flavonoid biosynthesis in Citrus reticulata cv. Chachiensis

Citrus reticulata cv. Chachiensis (CRC) is an important medicinal plant, its dried mature peels named "Guangchenpi", has been used as a traditional Chinese medicine to treat cough, indigestion, and lung diseases for several hundred years. However, the biosynthesis of the crucial natural products polymethoxylated flavonoids (PMFs) in CRC remains unclear. Here, we report a chromosome-scale genome assembly of CRC with the size of 314.96 Mb and a contig N50 of 16.22 Mb. Using multi-omics resources, we discover a putative caffeic acid O-methyltransferase (CcOMT1) that can transfer a methyl group to the 3-hydroxyl of natsudaidain to form 3,5,6,7,8,3',4'-heptamethoxyflavone (HPMF). Based on transient overexpression and virus-induced gene silencing experiments, we propose that CcOMT1 is a candidate enzyme in HPMF biosynthesis. In addition, a potential gene regulatory network associated with PMF biosynthesis is identified. This study provides insights into PMF biosynthesis and may assist future research on mining genes for the biosynthesis of plant-based medicines.

Haplotype-resolved T2T reference genomes for wild and domesticated accessions shed new insights into the domestication of jujube

Chinese jujube (Ziziphus jujuba Mill.) is one of the most important deciduous tree fruits in China, with substantial economic and nutritional value. Jujube was domesticated from its wild progenitor, wild jujube (Z. jujuba var. spinosa), and both have high medicinal value. Here we report the 767.81- and 759.24-Mb haplotype-resolved assemblies of a dry-eating 'Junzao' jujube (JZ) and a wild jujube accession (SZ), using a combination of multiple sequencing strategies. Each assembly yielded two complete haplotype-resolved genomes at the telomere-to-telomere (T2T) level, and ~81.60 and 69.07 Mb of structural variations were found between the two haplotypes within JZ and SZ, respectively. Comparative genomic analysis revealed a large inversion on each of chromosomes 3 and 4 between JZ and SZ, and numerous genes were affected by structural variations, some of which were associated with starch and sucrose metabolism. A large-scale population analysis of 672 accessions revealed that wild jujube originated from the lower reaches of the Yellow River and was initially domesticated at local sites. It spread widely and was then independently domesticated at the Shanxi-Shaanxi Gorge of the middle Yellow River. In addition, we identified some new selection signals regions on genomes, which are involved in the tissue development, pollination, and other aspects of jujube tree morphology and fertilization domestication. In conclusion, our study provides high-quality reference genomes of jujube and wild jujube and new insights into the domestication history of jujube.

Transcriptome differential expression analysis of defoliation of two different lemon varieties

'Allen Eureka' is a bud variety of Eureka lemon with excellent fruiting traits. However, it suffers from severe winter defoliation that leads to a large loss of organic nutrients and seriously affects the tree's growth and development as well as the yield of the following year, and the mechanism of its response to defoliation is still unclear. In order to investigate the molecular regulatory mechanisms of different leaf abscission periods in lemon, two lemon cultivars ('Allen Eureka' and 'Yunning No. 1') with different defoliation traits were used as materials. The petiole abscission zone (AZ) was collected at three different defoliation stages, namely, the pre-defoliation stage (CQ), the mid-defoliation stage (CZ), and the post-defoliation stage (CH). Transcriptome sequencing was performed to analyze the gene expression differences between these two cultivars. A total of 898, 4,856, and 3,126 differentially expressed genes (DEGs) were obtained in CQ, CZ, and CH, respectively, and the number of DEGs in CZ was the largest. GO analysis revealed that the DEGs between the two cultivars were mainly enriched in processes related to oxidoreductase, hydrolase, DNA binding transcription factor, and transcription regulator activity in the defoliation stages. KEGG analysis showed that the DEGs were concentrated in CZ and involved plant hormone signal transduction, phenylpropanoid biosynthesis, glutathione metabolism, and alpha-linolenic acid metabolism. The expression trends of some DEGs suggested their roles in regulating defoliation in lemon. Eight gene families were obtained by combining DEG clustering analysis and weighted gene co-expression network analysis (WGCNA), including β-glucosidase, AUX/IAA, SAUR, GH3, POD, and WRKY, suggesting that these genes may be involved in the regulation of lemon leaf abscission. The above conclusions enrich the research related to lemon leaf abscission and provide reliable data for the screening of lemon defoliation candidate genes and analysis of defoliation pathways.

Analysis of huanglongbing-associated RNA-seq data reveals disturbances in biological processes within Citrus spp. triggered by Candidatus Liberibacter asiaticus infection

Introduction

Huanglongbing (HLB), a disease that's ubiquitous worldwide, wreaks havoc on the citrus industry. The primary culprit of HLB is the gram-negative bacterium Candidatus Liberibacter asiaticus (CLas) that infects the phloem, but its damaging mechanism is yet to be fully understood.

Methods and results

In this study, a multitude of tools including weighted correlation network analysis (WGCNA), protein-protein interaction (PPI) network analysis and gene expression profiling are employed to unravel the intricacies of its pathogenesis. The investigation pinpoints various central genes, such as the ethylene-responsive transcription factor 9 (ERF9) and thioredoxin reductase 1 (TrxR1), that are associated with CLas invasion and resultant disturbances in numerous biological operations. Additionally, the study uncovers a range of responses through the detection of differential expressed genes (DEGs) across different experiments. The discovery of core DEGs leads to the identification of pivotal genes such as the sieve element occlusion (SEO) and the wall-associated receptor kinase-like 15 (WAKL15). PPI network analysis highlights potential vital proteins, while GO and KEGG pathway enrichment analysis illustrate a significant impact on multiple defensive and metabolic pathways. Gene set enrichment analysis (GSEA) indicates significant alterations in biological processes such as leaf senescence and response to biotic stimuli.

Discussion

This all-encompassing approach extends valuable understanding into the pathogenesis of CLas, potentially aiding future research and therapeutic strategies for HLB.

The lemon genome and DNA methylome unveil epigenetic regulation of citric acid biosynthesis during fruit development

Citric acid gives lemons their unique flavor, which impacts their sensory traits and market value. However, the intricate process of citric acid accumulation during lemon fruit growth remains incompletely understood. Here, we achieved a chromosomal-level genome assembly for the 'Xiangshui' lemon variety, spanning 364.85 Mb across nine chromosomes. This assembly revealed 27 945 genes and 51.37% repetitive sequences, tracing the divergence from citron 2.85 million years ago. DNA methylome analysis of lemon fruits across different developmental stages revealed significant variations in DNA methylation. We observed decreased CG and CHG methylation but increased CHH methylation. Notably, the expression of RdDM pathway-related genes increased with fruit development, suggesting a connection with elevated CHH methylation, which is potentially influenced by the canonical RdDM pathway. Furthermore, we observed that elevated CHH DNA methylation within promoters significantly influenced the expression of key genes, critically contributing to vital biological processes, such as citric acid accumulation. In particular, the pivotal gene phosphoenolpyruvate carboxykinase (ClPEPCK), which regulates the tricarboxylic acid cycle, was strikingly upregulated during fruit development, concomitant with increased CHH methylation in its promoter region. Other essential genes associated with citric acid accumulation, such as the MYB transcription factor (ClPH1/4/5) and ANTHOCYANIN 1 (ClAN1), were strongly correlated with DNA methylation levels. These results strongly indicate that DNA methylation crucially orchestrates the metabolic synthesis of citric acid. In conclusion, our study revealed dynamic changes in DNA methylation during lemon fruit development, underscoring the significant role of DNA methylation in controlling the citric acid metabolic pathway.

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.

Telomere-to-telomere and gap-free genome assembly of a susceptible grapevine species (Thompson Seedless) to facilitate grape functional genomics

Grapes are globally recognized as economically significant fruit trees. Among grape varieties, Thompson Seedless holds paramount influence for fresh consumption and for extensive applications in winemaking, drying, and juicing. This variety is one of the most efficient genotypes for grape genetic modification. However, the lack of a high-quality genome has impeded effective breeding efforts. Here, we present the high-quality reference genome of Thompson Seedless with all 19 chromosomes represented as 19 contiguous sequences (N50 = 27.1 Mb) with zero gaps and prediction of all telomeres and centromeres. Compared with the previous assembly (TSv1 version), the new assembly incorporates an additional 31.5 Mb of high-quality sequenced data with annotation of a total of 30 397 protein-coding genes. We also performed a meticulous analysis to identify nucleotide-binding leucine-rich repeat genes (NLRs) in Thompson Seedless and two wild grape varieties renowned for their disease resistance. Our analysis revealed a significant reduction in the number of two types of NLRs, TIR-NB-LRR (TNL) and CC-NB-LRR (CNL), in Thompson Seedless, which may have led to its sensitivity to many fungal diseases, such as powdery mildew, and an increase in the number of a third type, RPW8 (resistance to powdery mildew 8)-NB-LRR (RNL). Subsequently, transcriptome analysis showed significant enrichment of NLRs during powdery mildew infection, emphasizing the pivotal role of these elements in grapevine's defense against powdery mildew. The successful assembly of a high-quality Thompson Seedless reference genome significantly contributes to grape genomics research, providing insight into the importance of seedlessness, disease resistance, and color traits, and these data can be used to facilitate grape molecular breeding efforts.