Virus-induced genome editing (VIGE) leverages viral vectors to deliver CRISPR-Cas components into plants for robust and flexible trait engineering. We describe here a VIGE approach applying an RNA viral vector based on potato virus X (PVX) for genome editing of tomato, a mayor horticultural crop. Viral delivery of single-guide RNA into Cas9-expressing lines resulted in efficient somatic editing with indel frequencies up to 58%. By proof-of-concept VIGE of PHYTOENE DESATURASE (PDS) and plant regeneration from edited somatic tissue, we recovered loss-of-function pds mutant progeny displaying an albino phenotype. VIGE of STAYGREEN 1 (SGR1), a gene involved in fruit color variation, generated sgr1 mutant lines with recolored red-brown fruits and high chlorophyll levels. The obtained editing events were heritable, overall confirming the successful breeding of fruit color. Altogether, our VIGE approach offers great potential for accelerated functional genomics of tomato variation, as well as for precision breeding of novel tomato traits.

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.

CRISPR-Cas technologies allow for precise modifications in plant genomes and promise to revolutionize agriculture. These technologies depend on the delivery of editing components into plant cells and the regeneration of fully edited plants. In vegetatively propagated plants, such as grape, protoplast culture provides one of the best avenues for producing non-chimeric and transgene-free genome-edited plants. However, poor regeneration of plants from protoplasts has hindered their implementation for genome editing. Here, we report an efficient protocol for regenerating plants from protoplasts from multiple grape varieties. By encapsulating the protoplasts in calcium alginate beads and co-culturing them with feeder cultures, the protoplasts divide to form callus colonies that regenerate into embryos and ultimately plants. This protocol worked successfully in wine and table grape (Vitis vinifera) varieties, as well as grape rootstocks and the grapevine wild relative Vitis arizonica. Moreover, by transfecting protoplasts with CRISPR-plasmid or ribonucleoprotein (RNP) complexes, we regenerated albino plants with edits in VvPHYTOENE DESATURASE gene in three varieties and in V. arizonica. The results reveal the potential of this platform to facilitate genome editing in Vitis species.

Peach (Prunus persica) landrace has typical regional characteristics, strong environmental adaptability, and contains many valuable genes that provide the foundation for breeding excellent varieties. Therefore, it is necessary to assemble the genomes of specific landraces to facilitate the localization and utilization of these genes. Here, we de novo assembled a high-quality genome from an ancient blood-fleshed Chinese landrace Tianjin ShuiMi (TJSM) that originated from the China North Plain. The assembled genome size was 243.5 Mb with a contig N50 of 23.7 Mb and a scaffold N50 of 28.6 Mb. Compared with the reported peach genomes, our assembled TJSM genome had the largest number of specific structural variants (SVs) and long terminal repeat-retrotransposons (LTR-RTs). Among the LTR-RTs with the potential to regulate their host genes, we identified a 6688 bp LTR-RT (named it blood TE) in the promoter of NAC transcription factor-encoding PpBL, a gene regulating peach blood-flesh formation. The blood TE was not only co-separated with the blood-flesh phenotype but also associated with fruit maturity date advancement and different intensities of blood-flesh color formation. Our findings provide new insights into the mechanism underlying the development of the blood-flesh color and determination of fruit maturity date and highlight the potential of the TJSM genome to mine more variations related to agronomic traits in peach fruit.

Free amino acids (FAAs) positively determine the tea quality, notably theanine (Thea), endowing umami taste of tea infusion, which is the profoundly prevalent research in albino tea genetic resources. Therefore, 339 tea accessions were collected to study FAAs level for deciphering its variation and accumulation mechanism. Interestingly, alanine (Ala) and Thea which had the highest diversity index (H′) value among three varieties of Camellia sinensis (L.) O. Kuntze were significantly higher than wild relatives (P < 0.05). The intraspecific arginine (Arg) and glutamine (Gln) contents in C. sinensis var. assamica were significantly lower than sinensis and pubilimba varieties. Moreover, the importance of interdependencies operating across FAAs and chlorophyll levels were highlighted via the cell ultrastructure, metabolomics, and transcriptome analysis. We then determined that the association between phytochrome interacting factor 1 (CsPIF1) identified by weighted gene co-expression network analysis (WGCNA) and Thea content. Intriguingly, transient knock-down CsPIF1 expression increased Thea content in tea plant, and the function verification of CsPIF1 in Arabidopsis also indicated that CsPIF1 acts as a negative regulator of Thea content by mainly effecting the genes expression related to Thea biosynthesis, transport, and hydrolysis, especially glutamate synthase (CsGOGAT), which was validated to be associated with Thea content with a nonsynonymous SNP by Kompetitive Allele-Specific PCR (KASP). We also investigated the interspecific and geographical distribution of this SNP. Taken together, these results help us to understand and clarify the variation and profile of major FAAs in tea germplasms and promote efficient utilization in tea genetic improvement and breeding.

Due to the protracted transgenic timeline and low efficiency in stable genetic transformation of woody plants, there has been limited exploration of real-time organelle imaging within stable transgenic woody plant cells. Here, we established an efficient in vivo genetic transformation system for woody plants using an Agrobacterium rhizogenes-mediated approach. This system was successfully validated in multiple perennial woody species. Using citrus as a model, we introduced organelle-targeted fluorescent reporters via genetic transformation and investigated their subcellular localization and dynamics using advanced imaging techniques, such as confocal microscopy and live-cell imaging. Moreover, we subjected transgenic MT-GFP-labeled mitochondria in root cells to stress conditions simulating agricultural adversities faced by fruit crops. The stress-induced experiments revealed notable alterations in mitochondrial morphology. Our study contributes novel insights into membrane trafficking processes, protein localization dynamics, and cellular physiology in woody plants, while also providing stable and efficient genetic transformation methods for perennial woody species.

Clubroot disease caused by Plasmodiophora brassicae (P. brassicae) severely threatens the cultivation of Cruciferous plants, especially Chinese cabbage. Recently, resistance genes in plants have been reported to encode for a Ca²⁺-permeable channel in the plasma membrane, which can mediate the cytosolic Ca²⁺ increase in plant cells upon pathogen attack. However, the downstream Ca²⁺ sensor and decoder are still unknown. In this study, we identified the virulent and avirulent P. brassicae isolates (Pbs) of two near isogenic lines, CR 3-2 and CS 3-2, with CR 3-2 harboring clubroot resistant gene BraCRa. The transcriptomic analysis was then conducted with CR 3-2 after inoculating with virulent isolate PbE and avirulent isolate Pb4. From the differentially expressed genes of transcriptomic data, we identified a Ca²⁺-sensor encoding gene, BraCBL1.2, that was highly induced in CR 3-2 during infection by Pb4 but not by PbE. Moreover, GUS histochemical staining and subcellular localization analysis revealed that BraCBL1.2 was specifically expressed in the root hair cells of Arabidopsis and encoded a putative Ca²⁺ sensor localized in the plasma membrane. We also developed an assay to investigate the BraCRa-mediated hypersensitive response (HR) in tobacco leaves. The results suggest that BraCBL1.2 is involved in the BraCRa-mediated plant ETI immune response against P. brassicae. In addition, we verified that overexpression of BraCBL1.2 enhanced clubroot resistance in Arabidopsis. Collectively, our data identified the involvement of a Ca²⁺ sensor in BraCRa-mediated clubroot resistance in Chinese cabbage, providing a theoretical basis for further research on the resistance of Chinese cabbage to P. brassicae.

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.

Protoberberine alkaloids are a group of tetracyclic isoquinoline compounds known for their well-established antimicrobial and anti-inflammatory properties. The richness and diversity of protoberberine alkaloids accumulated in the Coptis genus necessitate a comprehensive examination of the biosynthetic machinery to understand their ecological significance. Here, from Coptis chinensis we identified CcCYP719A1, which could install a methylenedioxy bridge on either ring A or ring D of the protoberberine backbone, thus diverging metabolite flux towards the biosynthesis of various protoberberine components. We also obtained CcCYP719A2 and CcCYP719A3, which underwent positive selection after diverging from CcCYP719A1 and maintained specific catalytic activity on ring D. Further, we resolved the biosynthetic pathway of jatrorrhizine by identifying two demethylases, which could also modulate protoberberine composition by removing the C-3 methyl group and methylenedioxy bridge of ring D, allowing demethylated metabolites to be redirected into different routes. Moreover, we characterized 2-O-methyltransferase CcOMT1 and flavin-dependent oxidase CcTHBO, respectively responsible for the commonly observed 2-O-methylation and aromatic ring-C assembly in protoberberine alkaloids. Overall, this study reveals an interconnected metabolite network from which diverse protoberberine alkaloids originate. It provides valuable insights into the existence of undiscovered protoberberine components, and paves the way for the targeted production of desired protoberberine components for potential therapeutic development.

Scutellaria baicalensis Georgi, also known as huang-qin in traditional Chinese medicine, is a widely used herbal remedy due to its anticancer, antivirus, and hepatoprotective properties. The S. baicalensis genome was sequenced many years ago; by contrast, the proteome as the executer of most biological processes of S. baicalensis in the aerial parts, as well as the secondary structure of the roots (xylem, phloem, and periderm), is far less comprehensively characterized. Here we attempt to depict the molecular landscape of the non-model plant S. baicalensis through a multi-omics approach, with the goal of constructing a highly informative and valuable reference dataset. Furthermore, we provide an in-depth characterization dissection to explain the two distinct flavonoid biosynthesis pathways that exist in the aerial parts and root, at the protein and phosphorylated protein levels. Our study provides detailed spatial proteomic and phosphoproteomic information in the context of secondary structures, with implications for the molecular profiling of secondary metabolite biosynthesis in non-model medicinal plants.

A high-quality reference genome is indispensable for resolving biologically essential traits. Ficus hispida is a dioecious plant. A complete Ficus reference genome will be crucial for understanding their sex evolution and important biological characteristics, such as aerial roots, mutualistic symbiosis with ficus-wasps, and fruiting from old stems. Here, we generated a telomere-to-telomere (T2T) genome for F. hispida using PacBio HiFi and Oxford Nanopore Ultra-long sequencing technologies. The genome contiguity and completeness has shown improvement compared with the previously released genome, with the annotation of six centromeres and 28 telomeres. We have refined our previously reported 2-Mb male-specific region into a 7.2-Mb genomic region containing 51 newly predicted genes and candidate sex-determination genes AG2 and AG3. Many of these genes showed extremely low expression, likely attributed to hypermethylation in the gene body and promoter regions. Gene regulatory networks (GRNs) revealed that AG2 and AG3 are related to the regulation of stamen development in male flowers, while the AG1 gene is responsible for regulating female flowers’ defense responses and secondary metabolite processes. Comparative analysis of GRNs showed that the NAC, WRKY, and MYB transcription factor families dominate the female GRN, whereas the MADS and MYB transcription factor families are prevalent in the male GRN.

Podosphaera xanthii is the main causal agent of powdery mildew (PM) on Cucurbitaceae. In Cucumis melo, the Pm-w resistance gene, which confers resistance to P. xanthii, is located on chromosome 5 in a cluster of nucleotide-binding leucine-rich repeat receptors (NLRs). We used positional cloning and transgenesis, to isolate the Pm-w^{WMR 29} gene encoding a coiled-coil NLR (CC-NLR). Pm-w^{WMR 29} conferred high level of resistance to race 1 of PM and intermediate level of resistance to race 3 of PM. Pm-w^{WMR 29} turned out to be a homolog of the Aphis gossypii resistance gene Vat-1^{PI 161375}. We confirmed that Pm-w^{WMR 29} did not confer resistance to aphids, while Vat-1^{PI 161375} did not confer resistance to PM. We showed that both homologs were included in a highly diversified cluster of NLRs, the Vat cluster. Specific Vat-1^{PI 161375} and Pm-w^{WMR 29} markers were present in 10% to 13% of 678 accessions representative of wild and cultivated melon types worldwide. Phylogenic reconstruction of 34 protein homologs of Vat-1^{PI 161375} and Pm-w^WMR 29 identified in 24 melon accessions revealed an ancestor with four R65aa—a specific motif in the LRR domain, evolved towards aphid and virus resistance, while an ancestor with five R65aa evolved towards PM resistance. The complexity of the cluster comprising the Vat/Pm-w genes and its diversity in melon suggest that Vat homologs may contribute to the recognition of a broad range of yet to be identified pests and pathogens.

Populus cathayana Rehder, an indigenous poplar species of ecological and economic importance, is widely distributed in a high-elevation range from southwest to northeast China. Further development of this species as a sustainable poplar resource has been hindered by a lack of genome information the at the population level. Here, we produced a chromosome-level genome assembly of P. cathayana, covering 406.55 Mb (scaffold N50 = 20.86 Mb) and consisting of 19 chromosomes, with 35 977 protein-coding genes. Subsequently, we made a genomic variation atlas of 438 wild individuals covering 36 representative geographic areas of P. cathayana, which were divided into four geographic groups. It was inferred that the Northwest China regions served as the genetic diversity centers and a population bottleneck happened during the history of P. cathayana. By genotype-environment association analysis, 947 environment-association loci were significantly associated with temperature, solar radiation, precipitation, and altitude variables. We identified local adaptation genes involved in DNA repair and UV radiation response, among which UVR8, HY5, and CUL4 had key roles in high-altitude adaptation of P. cathayana. Predictions of adaptive potential under future climate conditions showed that P. cathayana populations in areas with drastic climate change were anticipated to have greater maladaptation risk. These results provide comprehensive insights for understanding wild poplar evolution and optimizing adaptive potential in molecular breeding.

Gray mold caused by Botrytis cinerea is one of the major threats in lily production. However, limited information is available about the underlying defense mechanism against B. cinerea in lily. Here, we characterized a nuclear-localized class A heat stress transcription factor (HSF)-LlHSFA4 from lily (Lilium longiflorum), which positively regulated the response to B. cinerea infection. LlHSFA4 transcript and its promoter activity were increased by B. cinerea infection in lily, indicating its involvement in the response to B. cinerea. Virus-induced gene silencing (VIGS) of LlHSFA4 impaired the resistance of lily to B. cinerea. Consistent with its role in lily, overexpression of LlHSFA4 in Arabidopsis (Arabidopsis thaliana) enhanced the resistance of transgenic Arabidopsis to B. cinerea infection. Further analysis showed that LlWRKY33 directly activated LlHSFA4 expression. We also found that both LlHSFA4 and LlWRKY33 positively regulated plant response to B. cinerea through reducing cell death and H₂O₂ accumulation and activating the expression of the reactive oxygen species (ROS) scavenging enzyme gene LlCAT2 (Catalase 2) by binding its prompter, which might contribute to reducing H₂O₂ accumulation in the infected area. Taken together, our data suggested that there may be a LlWRKY33-LlHSFA4-LlCAT2 regulatory module which confers B. cinerea resistance via reducing cell death and the ROS accumulation.