The draft genome of sweet orange (Citrus sinensis)

Genome-wide identification and expression analysis of the GAD family reveal their involved in embryogenesis and hormones responses in Dimocarpus longan Lour

Glutamate decarboxylase (GAD) is involved in GABA metabolism and plays an essential regulatory role in plant growth, abiotic stresses, and hormone response. This study investigated the expression mechanism of the GAD family during longan early somatic embryogenesis (SE) and identified 6 GAD genes based on the longan genome. Homology analysis indicated that DlGAD genes had a closer relationship with dicotyledonous plants. The analysis of cis-acting elements in the promoter region suggests that the GAD genes were associated with various stress responses and hormones. RNA sequencing (RNA-Seq) and the qRT-PCR data indicated that most DlGAD genes were highly expressed in the incomplete compact pro-embryogenic cultures (ICpEC) and upregulated in longan embryogenic callus (EC) after treatments with 2,4-D, high temperature (35 °C), IAA, and ABA. Moreover, the RNA-Seq analysis also revealed that DlGADs exhibit different expression patterns in various tissues and organs. The subcellular localization results showed that DlGAD5 was localized in the cytoplasm, suggesting that it played a role in the cytoplasm. Transient overexpression of DlGAD5 enhanced the expression levels of DlGADs and increased the activity of glutamate decarboxylase in longan embryogenic callus (EC), while the content of glutamic acid decreased. Thus, the DlGAD gene can play an important role in the early somatic embryogenesis of longan by responding to hormones such as IAA and ABA. DlGAD5 can affect the growth and development of longan by stimulating the expression of the DlGAD gene family, thereby increasing the GAD activity in the early SE of longan, participating in hormone synthesis and signaling pathways.

Identification of DEMETER-like DNA demethylase gene family in citrus and their role in drought stress-adaptive responses

DEMETER-Like DNA demethylases (DMLs) are epigenetic regulators of many developmental and biological processes in plants. No comprehensive information about the DML gene family in citrus is available to date. Here, a total of three DML genes in the genomes of Citrus sinensis (named CsDML1-3) and C. clementina (named CcDML1-3) were identified and analyzed. They encode hydrophilic and relatively large proteins, with prediction of nuclear localization, containing the conserved domains and motifs typical of plant DMLs. Protein interaction network analysis suggested that they interact primarily with proteins related to the maintenance of DNA methylation and remodeling of chromatin. Analysis of their promoter regions led to the identification of several cis-acting regulatory elements involved in stress response, including drought, heat and cold stresses. The presence of several miRNA targets and potential phosphorylation sites suggest that their expression is also regulated at post-transcriptional and post-translational levels. RNA-Seq data and quantitative real-time PCR analysis showed a low and drought-regulated gene expression of the citrus DMLs in different plant tissues. CsDML1 and CsDML3 were also differentially regulated by deficit irrigation in fruits at different developmental stages, with a positive and significant correlation found between CsDML1 and PHYTOENE SYNTHASE (PSY) and between CsDML3 and ATP CITRATE LYASEs (ACLs) and ZETA-CAROTENE DESATURASE (ZDS) gene expression. These results indicate that the citrus DMLs are potentially functional enzymes involved in developmental processes and drought stress-adaptive responses, providing a useful reference for further investigation of their functions and applications on the citrus improvement.

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.

Systematic analysis and functional verification of citrus ascorbate peroxidases reveal that CsAPX01 and CsAPX02 negatively regulate citrus bacterial canker through the hydrogen peroxide regulation

Ascorbate peroxidases (APXs) are antioxidant enzymes that play vital roles in redox homeostasis in plants. Citrus is susceptible to infection by Xanthomonas citri subsp. citri (Xcc), resulting in citrus bacterial canker (CBC). The present study used bioinformatic and expression analyses to investigate the APX family in Citrus sinensis. Bioinformatic research revealed the chromosomal locations, phylogeny, gene structure, promoter elements, functional domains, conserved motifs, and most likely physicochemical properties of the sequences. Six APXs clustered in three groups were identified, with each protein containing a single peroxidase domain. The promoter regions contained a variety of transcription factor-binding and hormone-response components. Xcc infection induced different CsAPX01 and CsAPX02 expressions in the CBC-susceptible Wanjincheng and CBC-resistant Kumquat varieties. Subcellular localization and transient expression showed that CsAPX01 and CsAPX02 were expressed in the cytoplasm and nucleus and had hydrogen peroxide (H2O2)-scavenging activity. Virus-induced gene silencing (VIGS) of CsAPX01 and CsAPX02 resulted in strong resistance to CBC and H2O2 bursts without effects on the plant phenotype. The current study focused on investigating and characterizing the citrus APX family. It was found that CsAPX01 and CsAPX02 exacerbated CBC by altering the balance of H2O2. These findings emphasize the importance of APXs in enhancing plant resistance to pathogens.

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.

Genome-wide characterization of RR gene family members in Zanthoxylum armatum and the subsequent functional characterization of the C-type RR

Response Regulators (RRs) are crucial regulators in plant development and stress responses, comprising A-type, B-type, C-type, and pseudo-RR subfamilies. However, previous studies have often focused on specific subfamilies, which restricts our understanding of the complete RR gene family. In this study, we conducted a comprehensive analysis of 63 RR members from Zanthoxylum armatum, using phylogenetic relationships, motif composition, cis-acting elements, gene duplication and collinearity analyses. Segmental repeats among ZaRR genes enhanced the various environmental adaptabilities of Z. armatum, and the B-type ZaRR exhibited significant collinearity with the RRs in P. trichocarpa and C. sinensis. Cis-element analysis indicated ZaRRs play a significant role in abiotic stress and phytohormone pathways, particularly in light, drought, cold, abscisic acid (ABA) and salicylic acid (SA) responses. Abundant Ethylene Response Factor (ERF) and reproduction-associated binding sites in ZaRR promoters suggested their roles in stress and reproductive processes. A-type ZaRRs were implicated in plant vegetative and reproductive growth, whereas B-type ZaRRs contributed to both growth and stress responses. C-type ZaRRs were associated with plant reproductive growth, whereas pseudo-RRs may function in plant stress responses, such as water logging, cold, and response to ethylene (ETH), SA, and jasmonic acid (JA). Ectopic expression of ZaRR24, a C-type RR, inhibits growth, induces early flowering, and shortens fruit length in Arabidopsis. ZaRR24 overexpression also affected the expression of A- and B-type RRs, as well as floral meristem and organ identity genes. These findings establish a solid and comprehensive foundation for RR gene research in Z. armatum, and provide a platform for investigating signal transduction in other woody plants.

Elucidation of mechanisms underlying active oxygen burst in Citrus sinensis after Diaporthe citri infection using transcriptome analysis

Introduction

Reactive oxygen species (ROS) generation is a common disease defense mechanism in plants. However, it is unclear whether Citrus host activates defense response against Diaporthe citri causing citrus melanose disease by producing ROS, and the underlying molecular mechanisms are unknown.

Methods

DAB staining and RNA-Seq technology were used to compare the active oxygen burst and differential gene expression, respectively, in uninfected and infected Citrus sinensis leaves at different time points during D. citri infection in vivo. The functions of CsRBOH (a significant DEG) were confirmed in N. benthamiana through the Agrobacterium-mediated transient expression system.

Results

DAB staining indicated that C. sinensis initiated defense against D. citri infection within 24 h by generating ROS. Illumina sequencing revealed 25,557 expressed genes of C. sinensis. The most upregulated DEGs (n = 1,570) were identified 72 h after fungal inoculation (sample denoted as CD72). In the CD72 vs. Cs (samples at 0 h after fungal inoculation) comparison, the KEGG pathway category with the highest number of genes (n = 62) and most significant enrichment was Protein processing in endoplasmic reticulum, followed by Glutathione metabolism and MAPK signaling pathway-plant. GO analysis revealed that the DEGs of CD72 vs. Cs related to active oxygen burst and chitin recognition were significantly grouped into the regulation of biological processes and molecular functions, with GO terms including response to ROS, response to fungus, and oxidoreductase activity. Remarkably, CsRBOH was significantly enriched in the GO and KEGG analyses, and its expression pattern in qRT-PCR and DAB staining results were consistent. Among the 63 ROS-related DEGs, HSP genes and genes associated with the peroxidase family were highly significant as revealed by protein-protein interaction networks. Furthermore, ROS accumulation, cell death, and upregulation of defense-related genes were observed in N. benthamiana leaves with CsRBOH expressed through the Agrobacterium-mediated transient expression system.

Conclusion

Our findings suggested that C. sinensis activates CsRBOH and ROS-related genes, leading to ROS accumulation to resist the invasion by D. citri. This study laid the foundation for future research on molecular mechanisms and breeding of C. sinensis cultivars resistant to citrus melanose.

Triumphs of genomic-assisted breeding in crop improvement

Conventional breeding approaches have played a significant role in meeting the food demand remarkably well until now. However, the increasing population, yield plateaus in certain crops, and limited recombination necessitate using genomic resources for genomics-assisted crop improvement programs. As a result of advancements in the next-generation sequence technology, GABs have developed dramatically to characterize allelic variants and facilitate their rapid and efficient incorporation in crop improvement programs. Genomics-assisted breeding (GAB) has played an important role in harnessing the potential of modern genomic tools, exploiting allelic variation from genetic resources and developing cultivars over the past decade. The availability of pangenomes for major crops has been a significant development, albeit with varying degrees of completeness. Even though adopting these technologies is essentially determined on economic grounds and cost-effective assays, which create a wealth of information that can be successfully used to exploit the latent potential of crops. GAB has been instrumental in harnessing the potential of modern genomic resources and exploiting allelic variation for genetic enhancement and cultivar development. GAB strategies will be indispensable for designing future crops and are expected to play a crucial role in breeding climate-smart crop cultivars with higher nutritional value.

Longitudinal Transcriptomic, Proteomic, and Metabolomic Response of Citrus sinensis to Diaphorina citri Inoculation of Candidatus Liberibacter asiaticus

Huanglongbing (HLB) is a fatal citrus disease that is currently threatening citrus varieties worldwide. One putative causative agent, Candidatus Liberibacter asiaticus (CLas), is vectored by Diaphorina citri, known as the Asian citrus psyllid (ACP). Understanding the details of CLas infection in HLB disease has been hindered by its Candidatus nature and the inability to confidently detect it in diseased trees during the asymptomatic stage. To identify early changes in citrus metabolism in response to inoculation of CLas using its natural psyllid vector, leaves from Madam Vinous sweet orange (Citrus sinensis (L.) Osbeck) trees were exposed to CLas-positive ACP or CLas-negative ACP and longitudinally analyzed using transcriptomics (RNA sequencing), proteomics (liquid chromatography-tandem mass spectrometry; data available in Dryad: 10.25338/B83H1Z), and metabolomics (proton nuclear magnetic resonance). At 4 weeks postexposure (wpe) to psyllids, the initial HLB plant response was primarily to the ACP and, to a lesser extent, the presence or absence of CLas. Additionally, analysis of 4, 8, 12, and 16 wpe identified 17 genes and one protein as consistently differentially expressed between leaves exposed to CLas-positive ACP versus CLas-negative ACP. This study informs identification of early detection molecular targets and contributes to a broader understanding of vector-transmitted plant pathogen interactions.

Citrus psorosis virus 24K protein inhibits the processing of miRNA precursors by interacting with components of the biogenesis machinery

Sweet orange (Citrus sinensis) is one of the most important fruit crops worldwide. Virus infections in this crop can interfere with cellular processes, causing dramatic economic losses. By performing RT-qPCR analyses, we demonstrated that citrus psorosis virus (CPsV)-infected orange plants exhibited higher levels of unprocessed microRNA (miRNA) precursors than healthy plants. This result correlated with the reported reduction of mature miRNAs species. The protein 24K, the CPsV suppressor of RNA silencing (VSR), interacts with miRNA precursors in vivo. Thus, this protein becomes a candidate responsible for the increased accumulation of unprocessed miRNAs. We analyzed 24K RNA-binding and protein-protein interaction domains and described patterns of its subcellular localization. We also showed that 24K colocalizes within nuclear D-bodies with the miRNA biogenesis proteins DICER-LIKE 1 (DCL1), HYPONASTIC LEAVES 1 (HYL1), and SERRATE (SE). According to the results of bimolecular fluorescence complementation and co-immunoprecipitation assays, the 24K protein interacts with HYL1 and SE. Thus, 24K may inhibit miRNA processing in CPsV-infected citrus plants by direct interaction with the miRNA processing complex. This work contributes to the understanding of how a virus can alter the regulatory mechanisms of the host, particularly miRNA biogenesis and function.IMPORTANCESweet oranges can suffer from disease symptoms induced by virus infections, thus resulting in drastic economic losses. In sweet orange plants, CPsV alters the accumulation of some precursors from the regulatory molecules called miRNAs. This alteration leads to a decreased level of mature miRNA species. This misregulation may be due to a direct association of one of the viral proteins (24K) with miRNA precursors. On the other hand, 24K may act with components of the cell miRNA processing machinery through a series of predicted RNA-binding and protein-protein interaction domains.

Haplotype-resolved genome assembly for tetraploid Chinese cherry (Prunus pseudocerasus) offers insights into fruit firmness

Chinese cherry (Prunus pseudocerasus) holds considerable importance as one of the primary stone fruit crops in China. However, artificially improving its traits and genetic analysis are challenging due to lack of high-quality genomic resources, which mainly result from difficulties associated with resolving its tetraploid and highly heterozygous genome. Herein, we assembled a chromosome-level, haplotype-resolved genome of the cultivar 'Zhuji Duanbing', comprising 993.69 Mb assembled into 32 pseudochromosomes using PacBio HiFi, Oxford Nanopore, and Hi-C. Intra-haplotype comparative analyses revealed extensive intra-genomic sequence and expression consistency. Phylogenetic and comparative genomic analyses demonstrated that P. pseudocerasus was a stable autotetraploid species, closely related to wild P. pusilliflora, with the two diverging ~18.34 million years ago. Similar to other Prunus species, P. pseudocerasus underwent a common whole-genome duplication event that occurred ~139.96 million years ago. Because of its low fruit firmness, P. pseudocerasus is unsuitable for long-distance transportation, thereby restricting its rapid development throughout China. At the ripe fruit stage, P. pseudocerasus cv. 'Zhuji Duanbing' was significantly less firm than P. avium cv. 'Heizhenzhu'. The difference in firmness is attributed to the degree of alteration in pectin, cellulose, and hemicellulose contents. In addition, comparative transcriptomic analyses identified GalAK-like and Stv1, two genes involved in pectin biosynthesis, which potentially caused the difference in firmness between 'Zhuji Duanbing' and 'Heizhenzhu'. Transient transformations of PpsGalAK-like and PpsStv1 increase protopectin content and thereby enhance fruit firmness. Our study lays a solid foundation for functional genomic studies and the enhancement of important horticultural traits in Chinese cherries.

Multiomics integrated with sensory evaluations to identify characteristic aromas and key genes in a novel brown navel orange (Citrus sinensis)

'Zong Cheng' navel orange (ZC) is a brown mutant of Lane Late navel orange (LL) and emits a more pleasant odor than that of LL. However, the key volatile compound of this aroma and underlying mechanism remains unclear. In this study, sensory evaluations and volatile profiling were performed throughout fruit development to identify significant differences in sensory perception and metabolites between LL and ZC. It revealed that the sesquiterpene content varied significantly between ZC and LL. Based on aroma extract dilution and gas chromatography-olfactometry analyses, the volatile compound leading to the background aroma of LL and ZC is d-limonene, the orange note in LL was mainly attributed to octanal, whilst valencene, β-myrcene, and (E)-β-ocimene presented balsamic, sweet, and herb notes in ZC. Furthermore, Cs5g12900 and six potential transcription factors were identified as responsible for valencene accumulation in ZC, which is important for enhancing the aroma of ZC.

Haplotype-resolved chromosome-level genome assembly of Huyou (Citrus changshanensis)

Huyou (Citrus changshanensis) is a significant citrus species that originated in Zhejiang Province, China, where it is also primarily cultivated. It is valued for its distinctive flavor and notable health benefits, owing to its high content of bioactive compounds like naringin and limonin. However, the absence of a high quality reference genome has limited the exploration of these health-promoting compounds in Huyou and hindered research into the mechanisms behind its medicinal properties. In this study, we present a phased chromosome-level genome assembly of Huyou. By combining PacBio and Hi-C sequencing, we generated a primary genome assembly and two haplotypes, comprising nine pseudo-chromosomes, with sizes of 339.91 Mb, 323.51 Mb, and 311.89 Mb, respectively. By integrating transcriptome data and annotations of homologous species, we identified a total of 29,775 protein-coding genes in the genome of Huyou. Additionally, we detected lots of structural variants between the two haplotypes. This represents the first reference genome of Huyou, providing a valuable resource for future studies on its agricultural characteristics and medicinal applications.

Genetic Identification of Medicinal Citrus Cultivar 'Local Juhong' Using Molecular Markers and Genomics

The citrus cultivar 'Local Juhong', which has historically been used as a traditional Chinese medicinal material, originated in Yuanjiang County, Hunan Province.Its parental type and genetic background are indistinct as of yet. Morphological observation shows that 'Local Juhong' has a slight oblateness in fruit shape, a relatively smooth pericarp, a fine and slightly raised oil vacuole, and an inward concave at the blossom end. The tree form and fruit and leaf morphology of 'Local Juhong' are similar to those of 'Huangpi' sour orange. To reveal the genetic background of 'Local Juhong', 21 citrus accessions were evaluated using nuclear and chloroplast SSR markers and whole-genome SNP information. 'Local Juhong' was grouped with mandarins and sub-grouped with 'Miyagawa Wase' and 'Yanxi Wanlu' in a nuclear SSR analysis, which indicated that its pollen parent might be mandarins. It was closely clustered with orange and pummelo in the chloroplast SSR analysis. The genomic sequence similarity rate of 'Local Juhong' with mandarin and pummelo heterozygosity was 70.88%; the main part was the heterozygosity, except for the unknown (19.66%), mandarin (8.73%), and pummelo (3.9%) parts. Thus, 'Local Juhong' may be an F1 hybrid with pummelo as the female parent and mandarin as the male parent, sharing sisterhood with 'Huangpi' sour orange.

Ethnobotany and diversity of Citrus spp. (Rutaceae) as a source of "Kem-kem" traditional medicine used among the Karo sub-ethnic in North Sumatra, Indonesia

The Karo sub-ethnic is one of five Batak sub-ethnicities in the Karo Regency, North Sumatra Province, Indonesia. They are famous for their local knowledge about the traditional use of medicinal plants to treat various diseases. The "Kem-kem" traditional medicine is one of the traditional healing practices that involve using plants passed down through generations from their ancestors. One of the plant genus group in the Rutaceae family utilized in the traditional "Kem-kem" healing practice is a citrus known as "Rimo". This study aims to document the local knowledge about the diversity of Citrus spp. as Kem-kem's herbal medicinal plant. This study was conducted from April to July 2023 in the Kabanjahe and Berastagi districts, Karo Regency, North Sumatra. Data was collected using interviews with traditional healers, herbal medicine vendors, and direct observations at traditional markets, involving a total of 8 Citrus spp. The Karo uses "Rimo" with different local names as sources of traditional medicinal ingredients in practicing "Kem-kem". There are 15 local names comprising eight species of Citrus. Four are hybrids, i.e., Citrus x aurantiifolia (Christm.) Swingle, Citrus × aurantium L, Citrus × junos Siebold ex Yu.Tanaka, and Citrus × taitensis Risso. Two of the remaining species are recognized in infraspecific rank, one variety (Citrus medica var. sarcodactylis (Hoola van Nooten) Swingle) and one form (Citrus × aurantium f. deliciosa (Ten.) M.Hiroe). They were used as material sources for Kem-kem traditional medicine to treat at least nine health problems. There are two species with six local names included in the Least Concern (LC) category, namely C. medica (Rimo Gawang, Rimo Hantuantu, Rimo Kayu), C. medica var. sarcodactylis (Rimo Kuku Harimau), and C. medica (Rimo Telur Buaya), C. maxima (Burm.) Merr. (Rimo Malem). Nine local names are included in the Not Evaluated (NE) category, namely C. × junos (Rimo Kejaren), C. × taitensis (Rimo Jungga), C. × aurantium f. deliciosa (Rimo Keling), C. × aurantium (Rimo Kersik), Citrus hystrix DC. (Rimo Mukur), C. × taitensis (Rimo Puraga), C. × aurantium (Rimo Kalele), Citrus swinglei Burkill ex Harms (Rimo Pagar), and C. x aurantiifolia (Rimo Bunga). Rimo Kejaren (C. × junos) is a species that has the most benefits.

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.

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.

Alternative pathways leading to ascorbate biosynthesis in plants: lessons from the last 25 years

l-Ascorbic acid (AsA) is an antioxidant with important roles in plant stress physiology, growth, and development. AsA also plays an essential role in human health, preventing scurvy. Humans do not synthesize AsA, which needs to be supplied via a diet rich in fresh produce. Research efforts have provided progress in the elucidation of a complex metabolic network with at least four routes leading to AsA formation in plants. In this review, three alternative pathways, namely the d-galacturonate, the l-gulose, and the myo-inositol pathways, are presented with the supporting evidence of their operation in multiple plant species. We critically discuss feeding studies using precursors and their conversion to AsA in plant organs, and research where the expression of key genes encoding enzymes involved in the alternative pathways showed >100% AsA content increase in the transgenics and in many cases accompanied by enhanced tolerance to multiple stresses. We propose that the alternative pathways are vital in AsA production in response to stressful conditions and to compensate in cases where the flux through the d-mannose/l-galactose pathway is reduced. The genes and enzymes that have been characterized so far in these alternative pathways represent important tools that are being used to develop more climate-tolerant crops.

Genome-wide identification and characterization of the sweet orange (Citrus sinensis) basic helix-loop-helix (bHLH) family reveals a role for CsbHLH085 as a regulator of citrus bacterial canker resistance

Citrus bacterial canker (CBC) is a harmful bacterial disease caused by Xanthomonas citri subsp. citri (Xcc), negatively impacting citrus production worldwide. The basic helix-loop-helix (bHLH) transcription factor family plays crucial roles in plant development and stress responses. This study aimed to identify and annotate bHLH proteins encoded in the Citrus sinensis genome and explore their involvement and functional importance in regulating CBC resistance. A total of 135 putative CsbHLHs TFs were identified and categorized into 16 subfamilies. Their chromosomal locations, collinearity, and phylogenetic relationships were comprehensively analyzed. Upon Xcc strain YN1 infection, certain CsbHLHs were differentially regulated in CBC-resistant and CBC-sensitive citrus varieties. Among these, CsbHLH085 was selected for further functional characterization. CsbHLH085 was upregulated in the CBC-resistant citrus variety, was localized in the nucleus, and had a transcriptional activation activity. CsbHLH085 overexpression in Citrus significantly enhanced CBC resistance, accompanied by increased levels of salicylic acid (SA), jasmonic acid (JA), reactive oxygen species (ROS), and decreased levels of abscisic acid (ABA) and antioxidant enzymes. Conversely, CsbHLH085 virus-induced gene silencing resulted in opposite phenotypic and biochemical responses. CsbHLH085 silencing also affected the expression of phytohormone biosynthesis and signaling genes involved in SA, JA, and ABA signaling. These findings highlight the crucial role of CsbHLH085 in regulating CBC resistance, suggesting its potential as a target for biotechnological-assisted breeding citrus varieties with improved resistance against phytopathogens.

ACC SYNTHASE4 inhibits gibberellin biosynthesis and FLOWERING LOCUS T expression during citrus flowering

Flowering is an essential process in fruit trees. Flower number and timing have a substantial impact on the yield and maturity of fruit. Ethylene and gibberellin (GA) play vital roles in flowering, but the mechanism of coordinated regulation of flowering in woody plants by GA and ethylene is still unclear. In this study, a lemon (Citrus limon L. Burm) 1-aminocyclopropane-1-carboxylic acid synthase gene (CiACS4) was overexpressed in Nicotiana tabacum and resulted in late flowering and increased flower number. Further transformation of citrus revealed that ethylene and starch content increased, and soluble sugar content decreased in 35S:CiACS4 lemon. Inhibition of CiACS4 in lemon resulted in effects opposite to that of 35S:CiACS4 in transgenic plants. Overexpression of the CiACS4-interacting protein ETHYLENE RESPONSE FACTOR3 (CiERF3) in N. tabacum resulted in delayed flowering and more flowers. Further experiments revealed that the CiACS4-CiERF3 complex can bind the promoters of FLOWERING LOCUS T (CiFT) and GOLDEN2-LIKE (CiFE) and suppress their expression. Moreover, overexpression of CiFE in N. tabacum led to early flowering and decreased flowers, and ethylene, starch, and soluble sugar contents were opposite to those in 35S:CiACS4 transgenic plants. Interestingly, CiFE also bound the promoter of CiFT. Additionally, GA3 and 1-aminocyclopropanecarboxylic acid (ACC) treatments delayed flowering in adult citrus, and treatment with GA and ethylene inhibitors increased flower number. ACC treatment also inhibited the expression of CiFT and CiFE. This study provides a theoretical basis for the application of ethylene to regulate flower number and mitigate the impacts of extreme weather on citrus yield due to delayed flowering.

Identification and characterization of two P450 enzymes from Citrus sinensis involved in TMTT and DMNT biosyntheses and Asian citrus psyllid defense

The homoterpenes (3E)-4,8-dimethyl-1,3,7-nonatriene (DMNT) and (E,E)-4,8,12-trimethyl-1,3,7,11-tridecatetraene (TMTT) are the major herbivore-induced plant volatiles that help in defense directly by acting as repellants and indirectly by recruiting insects' natural enemies. In this study, DMNT and TMTT were confirmed to be emitted from citrus (Citrus sinensis) leaves infested with Asian citrus psyllid (Diaphorina citri Kuwayama; ACP), and two cytochrome P450 (CYP) genes (CsCYP82L1 and CsCYP82L2) were newly identified and characterized. Understanding the functions of these genes in citrus defense will help plan strategies to manage huanglongbing caused by Candidatus Liberibacter asiaticus (CLas) and spread by ACP. Quantitative real-time PCR (qPCR) analysis showed that CsCYP82L1 and CsCYP82L2 were significantly upregulated in citrus leaves after ACP infestation. Yeast recombinant expression and enzyme assays indicated that CsCYP82L1 and CsCYP82L2 convert (E)-nerolidol to DMNT and (E,E)-geranyllinalool to TMTT. However, citrus calluses stably overexpressing CsCYP82L1 generated only DMNT, whereas those overexpressing CsCYP82L2 produced DMNT and TMTT. Furthermore, ACPs preferred wild-type lemon (Citrus limon) over the CsCYP82L1-overexpressing line in dual-choice feeding assays and mineral oil over TMTT or DMNT in behavioral bioassays. Finally, yeast one-hybrid, electrophoretic mobility shift, and dual luciferase assays demonstrated that CsERF017, an AP2/ERF transcription factor, directly bound to the CCGAC motif and activated CsCYP82L1. Moreover, the transient overexpression of CsERF017 in lemon leaves upregulated CsCYP82L1 in the absence and presence of ACP infestation. These results provide novel insights into homoterpene biosynthesis in C. sinensis and demonstrate the effect of homoterpenes on ACP behavior, laying a foundation to genetically manipulate homoterpene biosynthesis for application in huanglongbing and ACP control.

Switching action modes of miR408-5p mediates auxin signaling in rice

MicroRNAs (miRNAs) play fundamental roles in many developmental and physiological processes in eukaryotes. MiRNAs in plants generally regulate their targets via either mRNA cleavage or translation repression; however, which approach plays a major role and whether these two function modes can shift remains elusive. Here, we identify a miRNA, miR408-5p that regulates AUXIN/INDOLE ACETIC ACID 30 (IAA30), a critical repressor in the auxin pathway via switching action modes in rice. We find that miR408-5p usually inhibits IAA30 protein translation, but in a high auxin environment, it promotes the decay of IAA30 mRNA when it is overproduced. We further demonstrate that IDEAL PLANT ARCHITECTURE1 (IPA1), an SPL transcription factor regulated by miR156, mediates leaf inclination through association with miR408-5p precursor promoter. We finally show that the miR156-IPA1-miR408-5p-IAA30 module could be controlled by miR393, which silences auxin receptors. Together, our results define an alternative auxin transduction signaling pathway in rice that involves the switching of function modes by miR408-5p, which contributes to a better understanding of the action machinery as well as the cooperative network of miRNAs in plants.

Genome-Wide Identification and Characterization of the Sweet Orange (Citrus sinensis) GATA Family Reveals a Role for CsGATA12 as a Regulator of Citrus Bacterial Canker Resistance

Citrus bacterial canker (CBC) is a severe bacterial infection caused by Xanthomonas citri subsp. citri (Xcc), which continues to adversely impact citrus production worldwide. Members of the GATA family are important regulators of plant development and regulate plant responses to particular stressors. This report aimed to systematically elucidate the Citrus sinensis genome to identify and annotate genes that encode GATAs and evaluate the functional importance of these CsGATAs as regulators of CBC resistance. In total, 24 CsGATAs were identified and classified into four subfamilies. Furthermore, the phylogenetic relationships, chromosomal locations, collinear relationships, gene structures, and conserved domains for each of these GATA family members were also evaluated. It was observed that Xcc infection induced some CsGATAs, among which CsGATA12 was chosen for further functional validation. CsGATA12 was found to be localized in the nucleus and was differentially upregulated in the CBC-resistant and CBC-sensitive Kumquat and Wanjincheng citrus varieties. When transiently overexpressed, CsGATA12 significantly reduced CBC resistance with a corresponding increase in abscisic acid, jasmonic acid, and antioxidant enzyme levels. These alterations were consistent with lower levels of salicylic acid, ethylene, and reactive oxygen species. Moreover, the bacteria-induced CsGATA12 gene silencing yielded the opposite phenotypic outcomes. This investigation highlights the important role of CsGATA12 in regulating CBC resistance, underscoring its potential utility as a target for breeding citrus varieties with superior phytopathogen resistance.

T2T reference genome assembly and genome-wide association study reveal the genetic basis of Chinese bayberry fruit quality

Chinese bayberry (Myrica rubra or Morella rubra; 2n = 16) produces fruit with a distinctive flavor, high nutritional, and economic value. However, previous versions of the bayberry genome lack sequence continuity. Moreover, to date, no large-scale germplasm resource association analysis has examined the allelic and genetic variations determining fruit quality traits. Therefore, in this study, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar 'Zaojia' using PacBio HiFi long reads. The resulting 292.60 Mb T2T genome, revealed 8 centromeric regions, 15 telomeres, and 28 345 genes. This represents a substantial improvement in the genome continuity and integrity of Chinese bayberry. Subsequently, we re-sequenced 173 accessions, identifying 6 649 674 single nucleotide polymorphisms (SNPs). Further, the phenotypic analyses of 29 fruit quality-related traits enabled a genome-wide association study (GWAS), which identified 1937 SNPs and 1039 genes significantly associated with 28 traits. An SNP cluster pertinent to fruit color was identified on Chr6: 3407532 to 5 153 151 bp region, harboring two MYB genes (MrChr6G07650 and MrChr6G07660), exhibiting differential expression in extreme phenotype transcriptomes, linked to anthocyanin synthesis. An adjacent, closely linked gene, MrChr6G07670 (MLP-like protein), harbored an exonic missense variant and was shown to increase anthocyanin production in tobacco leaves tenfold. This SNP cluster, potentially a quantitative trait locus (QTL), collectively regulates bayberry fruit color. In conclusion, our study presented a complete reference genome, uncovered a suite of allelic variations related to fruit-quality traits, and identified functional genes that could be harnessed to enhance fruit quality and breeding efficiency of bayberries.

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.

The Genome of Vitis zhejiang-adstricta Strengthens the Protection and Utilization of the Endangered Ancient Grape Endemic to China

Vitis zhejiang-adstricta (V. zhejiang-adstricta) is one of the most important and endangered wild grapes. It is a national key protected wild, rare and endangered ancient grape endemic to China and used as a candidate material for resistance breeding owing to its excellent significant disease resistance. Here, we present a high-quality chromosome-level assembly of V. zhejiang-adstricta (IB-VB-01), comprising 506.66 Mb assembled into 19 pseudo-chromosomes. The contig N50 length is 3.91 Mb with 31,196 annotated protein-coding genes. Comparative genome and evolutionary analyses illustrated that V. zhejiang-adstricta has a specific position in the evolution of East Asian Vitis and shared a common ancestor with Vitis vinifera during the divergence of the two species about 10.42 (between 9.34 and 11.12) Mya. The expanded gene families compared with those in plants were related to disease resistance, and constructed gene families were related to plant growth and primary metabolism. With the analysis of gene family expansion and contraction, the evolution of environmental adaptability and especially the NBS-LRR gene family of V. zhejiang-adstricta was elucidated based on the pathways of resistance genes (R genes), unique genes and structural variations. The near-complete and accurate diploid V. zhejiang-adstricta reference genome obtained herein serves as an important complement to wild grape genomes and will provide valuable genomic resources for investigating the genomic architecture of V. zhejiang-adstricta as well as for improving disease resistance breeding strategies in grape.

Chromosome-level genome assembly of navel orange cv. Gannanzao (Citrus sinensis Osbeck cv. Gannanzao)

Navel orange cv. Gannanzao is a variant of the navel orange cv. Newhall (Citrus sinensis Osbeck cv. Newhall) that exhibits an earlier maturation, making it commercially valuable. However, the mechanisms underlying its early maturation remain obscure. To address this question, we conducted genome sequencing and de novo assembly of navel orange cv. Gannanzao. The assembled genome sequence is 334.57 Mb in length with a GC content of 31.48%. It comprises 318 contigs (N50 = 3.23 Mb) and 187 scaffolds (N50 = 31.86 Mb). The Benchmarking Universal Single-Copy Orthologs test demonstrates 94.6% completeness. The annotation revealed 23,037 gene models, 164.95 Mb of repetitive sequences, and 2,554 noncoding RNAs. A comparative analysis identified 323 fruit ripening-related genes in navel orange cv. Gannanzao genome, while navel orange cv. Newhall genome contained 345 such genes. These genes were organized into 320 orthologous gene families, with 30.3% of them exhibiting differences in gene copy numbers between the 2 genomes. Additionally, we identified 15 fruit ripening-related genes that have undergone adaptive evolution, suggesting their potential role in advancing fruit maturation in navel orange cv. Gannanzao. Whole-genome sequencing and annotation of navel orange cv. Gannanzao provides a valuable resource to unravel the early maturation mechanism of citrus and enriches the genomic resources for citrus research.

Multi-omics analyses reveal the importance of chromoplast plastoglobules in carotenoid accumulation in citrus fruit

Chromoplasts act as a metabolic sink for carotenoids, in which plastoglobules serve as versatile lipoprotein particles. PGs in chloroplasts have been characterized. However, the features of PGs from non-photosynthetic plastids are poorly understood. We found that the development of chromoplast plastoglobules (CPGs) in globular and crystalloid chromoplasts of citrus is associated with alterations in carotenoid storage. Using Nycodenz density gradient ultracentrifugation, an efficient protocol for isolating highly purified CPGs from sweet orange (Citrus sinensis) pulp was established. Forty-four proteins were defined as likely comprise the core proteome of CPGs using comparative proteomics analysis. Lipidome analysis of different chromoplast microcompartments revealed that the nonpolar microenvironment within CPGs was modified by 35 triacylglycerides, two sitosterol esters, and one stigmasterol ester. Manipulation of the CPG-localized gene CsELT1 (esterase/lipase/thioesterase) in citrus calli resulted in increased lipids and carotenoids, which is further evidence that the nonpolar microenvironment of CPGs contributes to carotenoid accumulation and storage in the chromoplasts. This multi-feature analysis of CPGs sheds new light on the role of chromoplasts in carotenoid metabolism, paving the way for manipulating carotenoid content in citrus fruit and other crops.

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.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Genome-wide identification of the GRF family in sweet orange (Citrus sinensis) and functional analysis of the CsGRF04 in response to multiple abiotic stresses

BACKGROUND

Citrus is one of the most valuable fruits worldwide and an economic pillar industry in southern China. Nevertheless, it frequently suffers from undesirable environmental stresses during the growth cycle, which severely restricts the growth, development and yield of citrus. In plants, the growth-regulating factor (GRF) family of transcription factors (TF) is extensively distributed and plays an vital part in plant growth and development, hormone response, as well as stress adaptation. However, the systematic identification and functional analysis of GRF TFs in citrus have not been reported.

RESULTS

Here, a genome-wide identification of GRF TFs was performed in Citrus sinensis, 9 members of CsGRFs were systematically identified and discovered to be scattered throughout 5 chromosomes. Subsequently, physical and chemical properties, phylogenetic relationships, structural characteristics, gene duplication events, collinearity and cis-elements of promoter were elaborately analyzed. In particular, the expression patterns of the CsGRF genes in response to multiple phytohormone and abiotic stress treatments were investigated. Predicated on this result, CsGRF04, which exhibited the most differential expression pattern under multiple phytohormone and abiotic stress treatments was screened out. Virus-induced gene silencing (VIGS) technology was utilized to obtain gene silenced plants for CsGRF04 successfully. After the three stress treatments of high salinity, low temperature and drought, the CsGRF04-VIGS lines showed significantly reduced resistance to high salinity and low temperature stresses, but extremely increased resistance to drought stress.

CONCLUSIONS

Taken together, our findings systematically analyzed the genomic characterization of GRF family in Citrus sinensis, and excavated a CsGRF04 with potential functions under multiple abiotic stresses. Our study lay a foundation for further study on the function of CsGRFs in abiotic stress and hormone signaling response.

A metabolomics study in citrus provides insight into bioactive phenylpropanoid metabolism

Citrus fruits are widely consumed worldwide in juices or as fresh and provide a broad range of phytonutrients that are important for human health. Here, a citrus multi-omics resource is presented: comprehensive metabolic profiling of various citrus species was performed and metabolic profiles were compared among species, with a focus on the phenylpropanoid metabolic pathway. A metabolite-based genome-wide association analysis (mGWAS) of 154 pummelo accessions was performed using factored spectrally transformed linear mixed models (FaST-LMM) and efficient mixed-model association eXpedited (EMMAX), and the genetic and biochemical basis of metabolomic variation was comprehensively analysed. A metabolite-single nucleotide polymorphism-gene (metabolite-SNP-gene) interaction network was constructed for pummelo, and many candidate loci controlling the synthesis and regulation of bioactive compounds were identified; among these loci, three BAHD malonyltransferases were involved in the malonylation of flavonoid glycosides. Further investigation revealed that an R2R3-MYB transcription factor CgMYB1 positively controls the metabolism of phenylpropanoid molecules, particularly the flavonoid derivatives. This study provides valuable data resources on the metabolic regulatory networks of bioactive components in citrus, in addition to demonstrating an efficient method for metabolic pathway dissection and providing targets for future breeding work with the aim of improving nutritional value.

High-quality maple genome reveals duplication-facilitated leaf color diversity

Acer truncatum is a horticultural tree species with individuals that display either yellow or red leaves in autumn, giving it high ornamental and economic value. 'Lihong' of A. truncatum is an excellent cultivar due to its characteristic of having autumn leaves that turn a bright and beautiful shade of red, while its closely related cultivar 'Bunge' does not. However, the molecular mechanism underlying the color change in the cultivar 'Lihong' is still unclear. Here, we assembled a high-quality genome sequence of Acer truncatum 'Lihong' (genome size = 688 Mb, scaffold N50 = 9.14 Mb) with 28,438 protein-coding genes predicted. Through comparative genomic analysis, we found that 'Lihong' had experienced more tandem duplication events although it's a high degree of collinearity with 'Bunge'. Especially, the expansion of key enzymes in the anthocyanin synthesis pathway was significantly uneven between the two varieties, with 'Lihong' genome containing a significantly higher number of tandem/dispersed duplication key genes. Further transcriptomic, metabolomic, and molecular functional analyses demonstrated that several UFGT genes, mainly resulting from tandem/dispersed duplication, followed by the promoter sequence variation, may contribute greatly to the leaf color phenotype, which provides new insights into the mechanism of divergent anthocyanin accumulation process in the 'Lihong' and 'Bunge' with yellow leaves in autumn. Further, constitutive expression of two UFGT genes, which showed higher expression in 'Lihong', elevated the anthocyanin content. We proposed that the small-scale duplication events could contribute to phenotype innovation.

Direct Organogenesis of Citrus Cultivars from Shoot Tip Nodal Segments

BACKGROUND

Citrus cultivar improvement via conventional breeding strategies is impeded by factors related to its reproductive biology. The orange is a hybrid between pomelo (Citrus maxima) and mandarin (Citrus reticulata). Among various orange cultivars, Valencia oranges have a bit of bitter tang mixed in with their sweetness, as Navel oranges are, the most widely cultivated citrus species, quite sweeter, and also don't contain any seeds. Tangelo mandarin orange cultivar is a hybrid of C. reticulata × C. maxima or × C. paradisi.

OBJECTIVE

The present study was undertaken to optimize the hormonal composition of the media with regard to plant growth regulators for in vitro propagation of sweet orange cultivars from nodal segment explants.

METHODS

The nodal segment explants were collected from three citrus cultivars, Washington Navel, Valencia and Tangelo. Murashige and Skoog (MS) medium supplemented with sucrose and different concentrations of growth regulators was used for shoot proliferation and root induction, and the optimum medium composition was assessed. The patent for Citrus Tissue Culture was obtained from the Office of Research Affairs, Haramaya University.

RESULTS

The results indicate that the highest shoot response was recorded for Washington's navel with maximum shoot proliferation rate (99.75%), shoot number per explant (1.76), shoot length (10.70 cm), leaf number per explants (3.54) after three weeks of culture. In all experiments, no growth was observed for the basal MS medium. Phytohormone combinations of IAA (1.2 mg/L) and kinetin (2.0 mg/L) were found to be the best for shoot proliferation. Among the cultivars, there were significant differences for the highest rooting rate (81.255), root number (2.22), and root length (2.95 cm) variables for Washington Navel. The lowest rooting rate (48.45%), root number (1.47) and root length (2.26 cm) were observed for Valencia. The highest rooting rate (84.90%), root number per microshoot (2.22) and root length (3.05 cm) was on MS medium supplemented with 1.5 mg/L NAA.

CONCLUSION

A comparison of different concentrations of IAA and NAA on root induction of microshoots from nodal segments of citrus cultivars demonstrated NAA was a more effective hormone than IAA.

Identification and high-throughput genotyping of single nucleotide polymorphism markers in a non-model conifer (Abies nordmanniana (Steven) Spach)

Single nucleotide polymorphism (SNP) markers are powerful tools for investigating population structures, linkage analysis, and genome-wide association studies, as well as for breeding and population management. The availability of SNP markers has been limited to the most commercially important timber species, primarily due to the cost of genome sequencing required for SNP discovery. In this study, a combination of reference-based and reference-free approaches were used to identify SNPs in Nordmann fir (Abies nordmanniana), a species previously lacking genomic sequence information. Using a combination of a genome assembly of the closely related Silver fir (Abies alba) species and a de novo assembly of low-copy regions of the Nordmann fir genome, we identified a high density of reliable SNPs. Reference-based approaches identified two million SNPs in common between the Silver fir genome and low-copy regions of Nordmann fir. A combination of one reference-free and two reference-based approaches identified 250 shared SNPs. A subset of 200 SNPs were used to genotype 342 individuals and thereby tested and validated in the context of identity analysis and/or clone identification. The tested SNPs successfully identified all ramets per clone and five mislabeled individuals via identity and genomic relatedness analysis. The identified SNPs will be used in ad hoc breeding of Nordmann fir in Denmark.

Evolution and diversification of the O-methyltransferase (OMT) gene family in Solanaceae

O-methyltransferases (OMTs) are a group of enzymes involved in several fundamental biological processes in plants, including lignin biosynthesis, pigmentation, and aroma production. Despite the intensive investigation of the role of OMTs in plant secondary metabolism, the evolution and diversification of this gene family in Solanaceae remain poorly understood. Here, we conducted a genome-wide survey of OMT genes in six Solanaceae species, reconstructing gene phylogenetic trees, predicting the potential involvement in biological processes, and investigating the exon/intron structure and chromosomal location. We identified 57 caffeoyl-CoA OMTs (CCoAOMTs) and 196 caffeic acid OMTs (COMTs) in the studied species. We observed a conserved gene block on chromosome 2 that consisted of tandem duplicated copies of OMT genes. Our results suggest that the expansion of the OMT gene family in Solanaceae was driven by whole genome duplication, segmental duplication, and tandem duplication, with multiple genes being retained by neofunctionalization and subfunctionalization. This study represents an essential first step in unraveling the evolutionary history of OMTs in Solanaceae. Our findings deepen our understanding of the crucial role of OMTs in several biological processes and highlight their significance as potential biotechnological targets.

Elucidation of the melitidin biosynthesis pathway in pummelo

Specialized plant metabolism is a rich resource of compounds for drug discovery. The acylated flavonoid glycoside melitidin is being developed as an anti-cholesterol statin drug candidate, but its biosynthetic route in plants has not yet been fully characterized. Here, we describe the gene discovery and functional characterization of a new flavonoid gene cluster (UDP-glucuronosyltransferases (CgUGTs), 1,2 rhamnosyltransferase (Cg1,2RhaT), acyltransferases (CgATs)) that is responsible for melitidin biosynthesis in pummelo (Citrus grandis (L.) Osbeck). Population variation analysis indicated that the tailoring of acyltransferases, specific for bitter substrates, mainly determine the natural abundance of melitidin. Moreover, 3-hydroxy-3-methylglutaryl-CoA reductase enzyme inhibition assays showed that the product from this metabolic gene cluster, melitidin, may be an effective anti-cholesterol statin drug candidate. Co-expression of these clustered genes in Nicotiana benthamiana resulted in the formation of melitidin, demonstrating the potential for metabolic engineering of melitidin in a heterologous plant system. This study establishes a biosynthetic pathway for melitidin, which provides genetic resources for the breeding and genetic improvement of pummelo aimed at fortifying the content of biologically active metabolites.

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.

Highly efficient Agrobacterium rhizogenes-mediated hairy root transformation in citrus seeds and its application in gene functional analysis

Highly efficient genetic transformation technology is beneficial for plant gene functional research and molecular improvement breeding. However, the most commonly used Agrobacterium tumefaciens-mediated genetic transformation technology is time-consuming and recalcitrant for some woody plants such as citrus, hampering the high-throughput functional analysis of citrus genes. Thus, we dedicated to develop a rapid, simple, and highly efficient hairy root transformation system induced by Agrobacterium rhizogenes to analyze citrus gene function. In this report, a rapid, universal, and highly efficient hairy root transformation system in citrus seeds was described. Only 15 days were required for the entire workflow and the system was applicable for various citrus genotypes, with a maximum transformation frequency of 96.1%. After optimization, the transformation frequency of Citrus sinensis, which shows the lowest transformation frequency of 52.3% among four citrus genotypes initially, was increased to 71.4% successfully. To test the applicability of the hairy roots transformation system for gene functional analysis of citrus genes, we evaluated the subcellular localization, gene overexpression and gene editing in transformed hairy roots. Compared with the traditional transient transformation system performed in tobacco leaves, the transgenic citrus hairy roots displayed a more clear and specific subcellular fluorescence localization. Transcript levels of genes were significantly increased in overexpressing transgenic citrus hairy roots as compared with wild-type (WT). Additionally, hairy root transformation system in citrus seeds was successful in obtaining transformants with knocked out targets, indicating that the Agrobacterium rhizogenes-mediated transformation enables the CRISPR/Cas9-mediated gene editing. In summary, we established a highly efficient genetic transformation technology with non-tissue-culture in citrus that can be used for functional analysis such as protein subcellular localization, gene overexpression and gene editing. Since the material used for genetic transformation are roots protruding out of citrus seeds, the process of planting seedlings prior to transformation of conventional tissue culture or non-tissue-culture was eliminated, and the experimental time was greatly reduced. We anticipate that this genetic transformation technology will be a valuable tool for routine research of citrus genes in the future.

Prioritizing strategies for wheat biofortification: Inspiration from underutilized species

The relationship between malnutrition and climate change is still poorly understood but a comprehensive knowledge of their interactions is needed to address the global public health agenda. Limited studies have been conducted to propose robust and economic-friendly strategies to augment the food basket with underutilized species and biofortify the staples for nutritional security. Sea-buckthorn is a known "superfood" rich in vitamin C and iron content. It is found naturally in northern hemispherical temperate Eurasia and can be utilized as a model species for genetic biofortification in cash crops like wheat. This review focuses on the impacts of climate change on inorganic (iron, zinc) and organic (vitamin C) micronutrient malnutrition employing wheat as highly domesticated crop and processed food commodity. As iron and zinc are particularly stored in the outer aleurone and endosperm layers, they are prone to processing losses. Moreover, only 5% Fe and 25% Zn are bioavailable once consumed calling to enhance the bioavailability of these micronutrients. Vitamin C converts non-available iron (Fe3+) to available form (Fe2+) and helps in the synthesis of ferritin while protecting it from degradation at the same time. Similarly, reduced phytic acid content also enhances its bioavailability. This relation urges scientists to look for a common mechanism and genes underlying biosynthesis of vitamin C and uptake of Fe/Zn to biofortify these micronutrients concurrently. The study proposes to scale up the biofortification breeding strategies by focusing on all dimensions i.e., increasing micronutrient content and boosters (vitamin C) and simultaneously reducing anti-nutritional compounds (phytic acid). Mutually, this review identified that genes from the Aldo-keto reductase family are involved both in Fe/Zn uptake and vitamin C biosynthesis and can potentially be targeted for genetic biofortification in crop plants.

Inheritance and Quantitative Trait Loci Mapping of Aromatic Compounds from Clementine (Citrus × clementina Hort. ex Tan.) and Sweet Orange (C. × sinensis (L.) Osb.) Fruit Essential Oils

Despite their importance in food processing, perfumery and cosmetics, the inheritance of sweet orange aromatic compounds, as well as their yield in the fruit peel, has been little analyzed. In the present study, the segregation of aromatic compounds was studied in an F1 population of 77 hybrids resulting from crosses between clementine and blood sweet orange. Fruit-peel essential oils (PEOs) extracted by hydrodistillation were analyzed by gas chromatography coupled with flame ionization detection. Genotyping by sequencing was performed on the parents and the hybrids. The resulting "clementine × sweet blood orange" genetic map consists of 710 SNP markers distributed in nine linkage groups (LGs), representing the nine citrus chromosomes, and spanning 1054 centimorgans. Twenty quantitative trait loci (QTLs) were identified, explaining between 20.5 and 55.0% of the variance of the major aromatic compounds and PEO yield. The QTLs for monoterpenes and aliphatic aldehydes predominantly colocalized on LGs 5 and 8, as did the two QTLs for PEO yield. The sesquiterpene QTLs were located on LGs 1, 3, 6 and 8. The detection of major QTLs associated with the synthesis of aliphatic aldehydes, known for their strong aromatic properties, open the way for marker-assisted selection.

Genome of Phyllanthus emblica: the medicinal plant Amla with super antioxidant properties

Phyllanthus emblica or Indian gooseberry, commonly known as amla, is an important medicinal horticultural plant used in traditional and modern medicines. It bears stone fruits with immense antioxidant properties due to being one of the richest natural sources of vitamin C and numerous flavonoids. This study presents the first genome sequencing of this species performed using 10x Genomics and Oxford Nanopore Technology. The draft genome assembly was 519 Mbp in size and consisted of 4,384 contigs, N50 of 597 Kbp, 98.4% BUSCO score, and 37,858 coding sequences. This study also reports the genome-wide phylogeny of this species with 26 other plant species that resolved the phylogenetic position of P. emblica. The presence of three ascorbate biosynthesis pathways including L-galactose, galacturonate, and myo-inositol pathways was confirmed in this genome. A comprehensive comparative evolutionary genomic analysis including gene family expansion/contraction and identification of multiple signatures of adaptive evolution provided evolutionary insights into ascorbate and flavonoid biosynthesis pathways and stone fruit formation through lignin biosynthesis. The availability of this genome will be beneficial for its horticultural, medicinal, dietary, and cosmetic applications and will also help in comparative genomics analysis studies.

Autotetraploid Origin of Chinese Cherry Revealed by Chromosomal Karyotype and In Situ Hybridization of Seedling Progenies

Polyploidy is considered a driving force in plant evolution and diversification. Chinese cherry [Cerasus pseudocerasus (Lindl.) G.Don], an economically important fruit crop native to China, has evolved at the tetraploid level, with a few pentaploid and hexaploid populations. However, its auto- or allo-polyploid origin remains unclear. To address this issue, we analyzed the ploidy levels and rDNA chromosomal distribution in self- and open-pollinated seedling progenies of tetraploid and hexaploid Chinese cherry. Genomic in situ hybridization (GISH) analysis was conducted to reveal the genomic relationships between Chinese cherry and diploid relatives from the genus Cerasus. Both self- and open-pollinated progenies of tetraploid Chinese cherry exhibited tetraploids, pentaploids, and hexaploids, with tetraploids being the most predominant. In the seedling progenies of hexaploid Chinese cherry, the majority of hexaploids and a few pentaploids were observed. A small number of aneuploids were also observed in the seedling progenies. Chromosome 1, characterized by distinct length characteristics, could be considered the representative chromosome of Chinese cherry. The basic Chinese cherry genome carried two 5S rDNA signals with similar intensity, and polyploids had the expected multiples of this copy number. The 5S rDNA sites were located at the per-centromeric regions of the short arm on chromosomes 4 and 5. Three 45S rDNA sites were detected on chr. 3, 4 and 7 in the haploid complement of Chinese cherry. Tetraploids exhibited 12 signals, while pentaploids and hexaploids showed fewer numbers than expected multiples. Based on the GISH signals, Chinese cherry demonstrated relatively close relationships with C. campanulata and C. conradinae, while being distantly related to another fruiting cherry, C. avium. In combination with the above results, our findings suggested that Chinese cherry likely originated from autotetraploidy.

Evolutionary assessment of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes in citrus relatives with a specific focus on flowering

Phase transition and floral induction in citrus requires several years of juvenility after germination. Such a long period of juvenility has been a major hindrance to its genetic improvement program. Studies have shown that miR156 along with its downstream genes SQUAMOSA PROMOTER BINDING PROTEINS (SBP) and SBP-LIKE (SPL) mediate the phase transition and floral induction process in plants. Our current study has systematically analyzed SPLs in 15 different citrus-related species, systematically annotated them based on their close homology to their respective Arabidopsis orthologs, and confirmed the functional attributes of the selected members in floral precocity. The majority of the species harbored 15 SPLs. Their cis-element assessment suggested the involvement of the SPLs in diverse developmental and physiological processes in response to different biotic and abiotic cues. Among all, SPL5, SPL9, and SPL11 stood out as consistently differentially expressed SPLs in the adult and young tissues of different citrus-related species. Independent overexpression of their F. hindsii orthologs (FhSPL5, FhSPL9, and FhSPL11) brought an enhanced expression of endogenous FLOWERING LOCUS T leading to the significantly precocious flowering in transgenic Arabidopsis lines. Future study of the genes in the citrus plant itself is expected to conclude the assessments made in the current study.

Karyotype and Phylogenetic Relationship Analysis of Five Varieties and Cultivars of Zanthoxylum armatum Based on Oligo-FISH

Green prickly ash (Zanthoxylum armatum) has edible and medicinal value and is an economically significant plant in many countries. Z. armatum has many cultivars and varieties with similar phenotypes that are difficult to distinguish via traditional methods. In this study, we utilized oligo-FISH to distinguish five varieties and cultivars of Z. armatum on the basis of three oligonucleotide probes of 5S rDNA, (AG₃T₃)₃, and (GAA)₆. Karyotype analysis of the five varieties and cultivars of Z. armatum showed that the Z. armatum 'Tengjiao' karyotype formula was 2n = 2x = 98m with karyotype type 1C and an arm ratio of 4.3237, including two pairs of 5S rDNA signals and five pairs of (GAA)₆ signals. The karyotype formula of Z. armatum 'Youkangtengjiao' was 2n = 2x = 128m + 8sm with karyotype type 2B and an arm ratio of 3.5336, including three pairs of 5S rDNA signals and 17 pairs of (GAA)₆ signals. The karyotype formula of Z. armatum var. novemfolius was 2n = 2x = 134m + 2sm with karyotype type 1C and an arm ratio of 5.5224, including two pairs of 5S rDNA signals and eight pairs of (GAA)₆ signals. The karyotype formula of Z. armatum 'YT-03' was 2n = 2x = 2M + 128m + 4sm + 2st with karyotype type 2C and an arm ratio of 4.1829, including three pairs of 5S rDNA signals and nine pairs of (GAA)₆ signals. The karyotype formula of Z. armatum 'YT-06' was 2n = 2x = 126m + 10sm with cytotype 2B and an arm ratio of 3.3011, including three pairs of 5S rDNA signals and two pairs of (GAA)₆ signals. The five varieties and cultivars of Z. armatum had (AG₃T₃)₃ signals on all chromosomes. The chromosomal symmetry of Z. armatum 'Tengjiao' was high, whereas the chromosomal symmetry of Z. armatum 'YT-03' was low, with the karyotypes of the five materials showing a trend toward polyploid evolution. The phylogenetic relationship between Z. armatum 'Tengjiao' and Z. armatum var. novemfolius was the closest, while that between Z. armatum 'YT-03' and Z. armatum 'YT-06' was closer than with Z. armatum 'Youkangtengjiao' according to oligo-FISH. The results provided a karyotype profile and a physical map that contributes to the distinction of varieties and cultivars of Z. armatum and provides strategies for distinguishing other cultivated species.

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.

Citrus ACC synthase CiACS4 regulates plant height by inhibiting gibberellin biosynthesis

Dwarfism is an agronomic trait that has substantial effects on crop yield, lodging resistance, planting density, and a high harvest index. Ethylene plays an important role in plant growth and development, including the determination of plant height. However, the mechanism by which ethylene regulates plant height, especially in woody plants, remains unclear. In this study, a 1-aminocyclopropane-1-carboxylic acid synthase (ACC) gene (ACS), which is involved in ethylene biosynthesis, was isolated from lemon (Citrus limon L. Burm) and named CiACS4. Overexpression of CiACS4 resulted in a dwarf phenotype in Nicotiana tabacum and lemon and increased ethylene release and decreased gibberellin (GA) content in transgenic plants. Inhibition of CiACS4 expression in transgenic citrus significantly increased plant height compared with the controls. Yeast two-hybrid assays revealed that CiACS4 interacted with an ethylene response factor (ERF), CiERF3. Further experiments revealed that the CiACS4-CiERF3 complex can bind to the promoters of 2 citrus GA20-oxidase genes, CiGA20ox1 and CiGA20ox2, and suppress their expression. In addition, another ERF transcription factor, CiERF023, identified using yeast one-hybrid assays, promoted CiACS4 expression by binding to its promoter. Overexpression of CiERF023 in N. tabacum caused a dwarfing phenotype. CiACS4, CiERF3, and CiERF023 expression was inhibited and induced by GA3 and ACC treatments, respectively. These results suggest that the CiACS4-CiERF3 complex may be involved in the regulation of plant height by regulating CiGA20ox1 and CiGA20ox2 expression levels in citrus.

L-Ascorbic acid metabolism and regulation in fruit crops

L-Ascorbic acid (AsA) is more commonly known as vitamin C and is an indispensable compound for human health. As a major antioxidant, AsA not only maintains redox balance and resists biological and abiotic stress but also regulates plant growth, induces flowering, and delays senescence through complex signal transduction networks. However, AsA content varies greatly in horticultural crops, especially in fruit crops. The AsA content of the highest species is approximately 1,800 times higher than that of the lowest species. There have been significant advancements in the understanding of AsA accumulation in the past 20 years. The most noteworthy accomplishment was the identification of the critical rate-limiting genes for the 2 major AsA synthesis pathways (L-galactose pathway and D-galacturonic acid pathway) in fruit crops. The rate-limiting genes of the former are GMP, GME, GGP, and GPP, and the rate-limiting gene of the latter is GalUR. Moreover, APX, MDHAR, and DHAR are also regarded as key genes in degradation and regeneration pathways. Interestingly, some of these key genes are sensitive to environmental factors, such as GGP being induced by light. The efficiency of enhancing AsA content is high by editing upstream open reading frames (uORF) of the key genes and constructing multi-gene expression vectors. In summary, the AsA metabolism has been well understood in fruit crops, but the transport mechanism of AsA and the synergistic improvement of AsA and other traits is less known, which will be the focus of AsA research in fruit crops.

Methylome and transcriptome analysis of flowering branches building of Citrus plants induced by drought stress

Citrus plants exhibit positive floral response under water stress conditions, however, the mechanistic understanding of floral induction remains largely unexplored in water deficit. In this study, DNA methylomic and transcriptomic analyses were integrated to explore the flowering bud formation as well as branches building after light drought stress. While comparing with the conventional watering group (CK), the light drought group treated with five months (LD) showed a significant increase in the flowering branches, whereas an apparent decrease in vegetative branches. Global DNA methylation analysis showed that the LD Group acquired DNA methylation in more than 70090 genomic regions and lost DNA methylation in about 18421 genomic regions compared with normal watering group, this indicates that water deficiency leads to a global increase in the expression of DNA methylation in citrus. In the same time, we verified that the increase of DNA methylation level in LD group was correlated with the decrease of DNA demethylase related gene expression. Interestingly, in transcription analysis, it was found that the promoting flower genes of the LD group did not increase but decreased similarly with repressing genes, which is contrary to the intended result. Thus, we thought the lower expression of suppressors FLC and BFT were the key influencing factor to stimulate the flowering branches formation after LD treatment. Moreover, there was a strong negative correlation between the genes expression level and methylation level of the flowering induction/flower development genes. In general, we thought high global DNA methylation level induced by water deficit regulate the flowering branches building by reducing FLC and BFT genes expression.