Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding.

Hybrid potato breeding will transform the crop from a clonally propagated tetraploid to a seed-reproducing diploid. Historical accumulation of deleterious mutations in potato genomes has hindered the development of elite inbred lines and hybrids. Utilizing a whole-genome phylogeny of 92 Solanaceae and its sister clade species, we employ an evolutionary strategy to identify deleterious mutations. The deep phylogeny reveals the genome-wide landscape of highly constrained sites, comprising ∼2.4% of the genome. Based on a diploid potato diversity panel, we infer 367,499 deleterious variants, of which 50% occur at non-coding and 15% at synonymous sites. Counterintuitively, diploid lines with relatively high homozygous deleterious burden can be better starting material for inbred-line development, despite showing less vigorous growth. Inclusion of inferred deleterious mutations increases genomic-prediction accuracy for yield by 24.7%. Our study generates insights into the genome-wide incidence and properties of deleterious mutations and their far-reaching consequences for breeding.

Genome design of hybrid potato.

Reinventing potato from a clonally propagated tetraploid into a seed-propagated diploid, hybrid potato, is an important innovation in agriculture. Due to deleterious mutations, it has remained a challenge to develop highly homozygous inbred lines, a prerequisite to breed hybrid potato. Here, we employed genome design to develop a generation of pure and fertile potato lines and thereby the uniform, vigorous F1s. The metrics we applied in genome design included the percentage of genome homozygosity and the number of deleterious mutations in the starting material, the number of segregation distortions in the S1 population, the haplotype information to infer the break of tight linkage between beneficial and deleterious alleles, and the genome complementarity of the parental lines. This study transforms potato breeding from a slow, non-accumulative mode into a fast-iterative one, thereby potentiating a broad spectrum of benefits to farmers and consumers.

Genome evolution and diversity of wild and cultivated potatoes.

Potato (Solanum tuberosum L.) is the world's most important non-cereal food crop, and the vast majority of commercially grown cultivars are highly heterozygous tetraploids. Advances in diploid hybrid breeding based on true seeds have the potential to revolutionize future potato breeding and production1-4. So far, relatively few studies have examined the genome evolution and diversity of wild and cultivated landrace potatoes, which limits the application of their diversity in potato breeding. Here we assemble 44 high-quality diploid potato genomes from 24 wild and 20 cultivated accessions that are representative of Solanum section Petota, the tuber-bearing clade, as well as 2 genomes from the neighbouring section, Etuberosum. Extensive discordance of phylogenomic relationships suggests the complexity of potato evolution. We find that the potato genome substantially expanded its repertoire of disease-resistance genes when compared with closely related seed-propagated solanaceous crops, indicative of the effect of tuber-based propagation strategies on the evolution of the potato genome. We discover a transcription factor that determines tuber identity and interacts with the mobile tuberization inductive signal SP6A. We also identify 561,433 high-confidence structural variants and construct a map of large inversions, which provides insights for improving inbred lines and precluding potential linkage drag, as exemplified by a 5.8-Mb inversion that is associated with carotenoid content in tubers. This study will accelerate hybrid potato breeding and enrich our understanding of the evolution and biology of potato as a global staple food crop.

The Paris Climate Agreement and future sea-level rise from Antarctica.

The Paris Agreement aims to limit global mean warming in the twenty-first century to less than 2 degrees Celsius above preindustrial levels, and to promote further efforts to limit warming to 1.5 degrees Celsius1. The amount of greenhouse gas emissions in coming decades will be consequential for global mean sea level (GMSL) on century and longer timescales through a combination of ocean thermal expansion and loss of land ice2. The Antarctic Ice Sheet (AIS) is Earth's largest land ice reservoir (equivalent to 57.9 metres of GMSL)3, and its ice loss is accelerating4. Extensive regions of the AIS are grounded below sea level and susceptible to dynamical instabilities5-8 that are capable of producing very rapid retreat8. Yet the potential for the implementation of the Paris Agreement temperature targets to slow or stop the onset of these instabilities has not been directly tested with physics-based models. Here we use an observationally calibrated ice sheet-shelf model to show that with global warming limited to 2 degrees Celsius or less, Antarctic ice loss will continue at a pace similar to today's throughout the twenty-first century. However, scenarios more consistent with current policies (allowing 3 degrees Celsius of warming) give an abrupt jump in the pace of Antarctic ice loss after around 2060, contributing about 0.5 centimetres GMSL rise per year by 2100-an order of magnitude faster than today4. More fossil-fuel-intensive scenarios9 result in even greater acceleration. Ice-sheet retreat initiated by the thinning and loss of buttressing ice shelves continues for centuries, regardless of bedrock and sea-level feedback mechanisms10-12 or geoengineered carbon dioxide reduction. These results demonstrate the possibility that rapid and unstoppable sea-level rise from Antarctica will be triggered if Paris Agreement targets are exceeded.

Effective green treatment of sewage sludge from Fenton reactions: Utilizing MoS2 for sustainable resource recovery.

Effectively managing sewage sludge from Fenton reactions in an eco-friendly way is vital for Fenton technology's viability in pollution treatment. This study focuses on sewage sludge across various treatment stages, including generation, concentration, dehydration, and landfill, and employs chemical composite MoS2 to facilitate green resource utilization of all types of sludge. MoS2, with exposed Mo4+ and low-coordination sulfur, enhances iron cycling and creates an acidic microenvironment on the sludge surface. The MoS2-modified iron sludge exhibits outstanding (>95%) phenol and pollutant degradation in hydrogen peroxide and peroxymonosulfate-based Fenton systems, unlike unmodified sludge. This modified sludge maintains excellent Fenton activity in various water conditions and with multiple anions, allowing extended phenol degradation for over 14 d. Notably, the generated chemical oxygen demand (COD) in sludge modification process can be efficiently eliminated through the Fenton reaction, ensuring effluent COD compliance and enabling eco-friendly sewage sludge resource utilization.

A viral protein disrupts vacuolar acidification to facilitate virus infection in plants.

Vacuolar acidification is essential for vacuoles in diverse physiological functions. However, its role in plant defense, and whether and how pathogens affect vacuolar acidification to promote infection remain unknown. Here, we show that Barley stripe mosaic virus (BSMV) replicase γa, but not its mutant γaR569A , directly blocks acidification of vacuolar lumen and suppresses autophagic degradation to promote viral infection in plants. These were achieved via molecular interaction between γa and V-ATPase catalytic subunit B2 (VHA-B2), leading to disruption of the interaction between VHA-B2 and V-ATPase catalytic subunit E (VHA-E), which impairs the membrane localization of VHA-B2 and suppresses V-ATPase activity. Furthermore, a mutant virus BSMVR569A with the R569A point mutation possesses less viral pathogenicity. Interestingly, multiple viral infections block vacuolar acidification. These findings reveal that functional vacuolar acidification is required for plant antiviral defense and disruption of vacuolar acidification could be a general viral counter-defense strategy employed by multiple viruses.

Palmitoylation of γb protein directs a dynamic switch between Barley stripe mosaic virus replication and movement.

Viral replication and movement are intimately linked; however, the molecular mechanisms regulating the transition between replication and subsequent movement remain largely unknown. We previously demonstrated that the Barley stripe mosaic virus (BSMV) γb protein promotes viral replication and movement by interacting with the αa replicase and TGB1 movement proteins. Here, we found that γb is palmitoylated at Cys-10, Cys-19, and Cys-60 in Nicotiana benthamiana, which supports BSMV infection. Intriguingly, non-palmitoylated γb is anchored to chloroplast replication sites and enhances BSMV replication, whereas palmitoylated γb protein recruits TGB1 to the chloroplasts and forms viral replication-movement intermediate complexes. At the late stages of replication, γb interacts with NbPAT15 and NbPAT21 and is palmitoylated at the chloroplast periphery, thereby shifting viral replication to intracellular and intercellular movement. We also show that palmitoylated γb promotes virus cell-to-cell movement by interacting with NbREM1 to inhibit callose deposition at the plasmodesmata. Altogether, our experiments reveal a model whereby palmitoylation of γb directs a dynamic switch between BSMV replication and movement events during infection.

Zn2+-dependent association of cysteine-rich protein with virion orchestrates morphogenesis of rod-shaped viruses.

The majority of rod-shaped and some filamentous plant viruses encode a cysteine-rich protein (CRP) that functions in viral virulence; however, the roles of these CRPs in viral infection remain largely unknown. Here, we used barley stripe mosaic virus (BSMV) as a model to investigate the essential role of its CRP in virus morphogenesis. The CRP protein γb directly interacts with BSMV coat protein (CP), the mutations either on the His-85 site in γb predicted to generate a potential CCCH motif or on the His-13 site in CP exposed to the surface of the virions abolish the zinc-binding activity and their interaction. Immunogold-labeling assays show that γb binds to the surface of rod-shaped BSMV virions in a Zn2+-dependent manner, which enhances the RNA binding activity of CP and facilitates virion assembly and stability, suggesting that the Zn2+-dependent physical association of γb with the virion is crucial for BSMV morphogenesis. Intriguingly, the tightly binding of diverse CRPs to their rod-shaped virions is a general feature employed by the members in the families Virgaviridae (excluding the genus Tobamovirus) and Benyviridae. Together, these results reveal a hitherto unknown role of CRPs in the assembly and stability of virus particles, and expand our understanding of the molecular mechanism underlying virus morphogenesis.

RETICULON-LIKE PROTEIN B2 is a proviral factor co-opted for the biogenesis of viral replication organelles in plants.

Endomembrane remodeling to form a viral replication complex (VRC) is crucial for a virus to establish infection in a host. Although the composition and function of VRCs have been intensively studied, host factors involved in the assembly of VRCs for plant RNA viruses have not been fully explored. TurboID-based proximity labeling (PL) has emerged as a robust tool for probing molecular interactions in planta. However, few studies have employed the TurboID-based PL technique for investigating plant virus replication. Here, we used Beet black scorch virus (BBSV), an endoplasmic reticulum (ER)-replicating virus, as a model and systematically investigated the composition of BBSV VRCs in Nicotiana benthamiana by fusing the TurboID enzyme to viral replication protein p23. Among the 185 identified p23-proximal proteins, the reticulon family of proteins showed high reproducibility in the mass spectrometry data sets. We focused on RETICULON-LIKE PROTEIN B2 (RTNLB2) and demonstrated its proviral functions in BBSV replication. We showed that RTNLB2 binds to p23, induces ER membrane curvature, and constricts ER tubules to facilitate the assembly of BBSV VRCs. Our comprehensive proximal interactome analysis of BBSV VRCs provides a resource for understanding plant viral replication and offers additional insights into the formation of membrane scaffolds for viral RNA synthesis.

The MAPK-Alfin-like 7 module negatively regulates ROS scavenging genes to promote NLR-mediated immunity.

Nucleotide-binding leucine-rich repeat (NLR) receptor-mediated immunity includes rapid production of reactive oxygen species (ROS) and transcriptional reprogramming, which is controlled by transcription factors (TFs). Although some TFs have been reported to participate in NLR-mediated immune response, most TFs are transcriptional activators, and whether and how transcriptional repressors regulate NLR-mediated plant defenses remains largely unknown. Here, we show that the Alfin-like 7 (AL7) interacts with N NLR and functions as a transcriptional repressor. Knockdown and knockout of AL7 compromise N NLR-mediated resistance against tobacco mosaic virus, whereas AL7 overexpression enhances defense, indicating a positive regulatory role for AL7 in immunity. AL7 binds to the promoters of ROS scavenging genes to inhibit their transcription during immune responses. Mitogen-activated protein kinases (MAPKs), salicylic acid-induced protein kinase (SIPK), and wound-induced protein kinase (WIPK) directly interact with and phosphorylate AL7, which impairs the AL7-N interaction and enhances its DNA binding activity, which promotes ROS accumulation and enables immune activation. In addition to N, AL7 is also required for the function of other Toll interleukin 1 receptor/nucleotide-binding/leucine-rich repeats (TNLs) including Roq1 and RRS1-R/RPS4. Our findings reveal a hitherto unknown MAPK-AL7 module that negatively regulates ROS scavenging genes to promote NLR-mediated immunity.

Barley stripe mosaic virus γb protein disrupts chloroplast antioxidant defenses to optimize viral replication.

The plant antioxidant system plays important roles in response to diverse abiotic and biotic stresses. However, the effects of virus infection on host redox homeostasis and how antioxidant defense pathway is manipulated by viruses remain poorly understood. We previously demonstrated that the Barley stripe mosaic virus (BSMV) γb protein is recruited to the chloroplast by the viral αa replicase to enhance viral replication. Here, we show that BSMV infection induces chloroplast oxidative stress. The versatile γb protein interacts directly with NADPH-dependent thioredoxin reductase C (NTRC), a core component of chloroplast antioxidant systems. Overexpression of NbNTRC significantly impairs BSMV replication in Nicotiana benthamiana plants, whereas disruption of NbNTRC expression leads to increased viral accumulation and infection severity. To counter NTRC-mediated defenses, BSMV employs the γb protein to competitively interfere with NbNTRC binding to 2-Cys Prx. Altogether, this study indicates that beyond acting as a helicase enhancer, γb also subverts NTRC-mediated chloroplast antioxidant defenses to create an oxidative microenvironment conducive to viral replication.

VirusPredictor: XGBoost-based software to predict virus-related sequences in human data.

MOTIVATION: Discovering disease causative pathogens, particularly viruses without reference genomes, poses a technical challenge as they are often unidentifiable through sequence alignment. Machine learning prediction of patient high-throughput sequences unmappable to human and pathogen genomes may reveal sequences originating from uncharacterized viruses. Currently, there is a lack of software specifically designed for accurately predicting such viral sequences in human data. RESULTS: We developed a fast XGBoost method and software VirusPredictor leveraging an in-house viral genome database. Our two-step XGBoost models first classify each query sequence into one of three groups: infectious virus, endogenous retrovirus (ERV) or non-ERV human. The prediction accuracies increased as the sequences became longer, i.e. 0.76, 0.93, and 0.98 for 150-350 (Illumina short reads), 850-950 (Sanger sequencing data), and 2000-5000 bp sequences, respectively. Then, sequences predicted to be from infectious viruses are further classified into one of six virus taxonomic subgroups, and the accuracies increased from 0.92 to >0.98 when query sequences increased from 150-350 to >850 bp. The results suggest that Illumina short reads should be de novo assembled into contigs (e.g. ∼1000 bp or longer) before prediction whenever possible. We applied VirusPredictor to multiple real genomic and metagenomic datasets and obtained high accuracies. VirusPredictor, a user-friendly open-source Python software, is useful for predicting the origins of patients' unmappable sequences. This study is the first to classify ERVs in infectious viral sequence prediction. This is also the first study combining virus sub-group predictions. AVAILABILITY AND IMPLEMENTATION: www.dllab.org/software/VirusPredictor.html.

Multiplex gene editing reveals cucumber MILDEW RESISTANCE LOCUS O family roles in powdery mildew resistance.

Powdery mildew (PM) is one of the most widespread and prevalent diseases that affects a wide range of crops. In cucumber (Cucumis sativus L.), previous forward genetic studies have identified MILDEW RESISTANCE LOCUS O 8 (CsMLO8) as necessary but alone insufficient for cucumber PM resistance (PMR) and suggested the involvement of other members of the CsMLO family. However, the function of other CsMLO family members in cucumber remains largely unknown. Here, we developed a highly efficient multiplex gene editing system in cucumber to generate a series of Csmlo mutants from all the 13 family members. Systematic analysis of these mutants revealed growth effects of these CsMLO family members on development and PMR. Importantly, we obtained the Csmlo1/8/11 triple mutant with complete resistance to PM. Transcriptome and proteome analysis of PM-resistant Csmlo mutants suggested that the kinesin-like calmodulin-binding protein (KCBP)-interacting Ca2+-binding protein (CsKIC), calmodulin-like protein 28 (CsCML28), and Ca2+-dependent protein kinase 11 (CsCPK11)-mediated calcium signaling pathway is involved in PMR. CsMLO8 interacted directly with CsKIC, and the simultaneous silencing of both genes resulted in a phenotype that resembled the silencing of CsKIC alone. Silencing CsCML28 and CsCPK11 increased susceptibility to PM, whereas overexpressing CsCPK11 through genetic transformation enhanced cucumber's PMR, demonstrating their positive regulatory roles in PMR. Given the importance of PMR for cucurbit crops, this research provides unprecedented insights into the function of the proteins encoded by the CsMLO gene family as well as the plant defense response to PM pathogen.

COLMARppm: A Web Server Tool for the Accurate and Rapid Prediction of 1H and 13C NMR Chemical Shifts of Organic Molecules and Metabolites.

Despite rapid progress in metabolomics research, a major bottleneck is the large number of metabolites whose chemical structures are unknown or whose spectra have not been deposited in metabolomics databases. Nuclear magnetic resonance (NMR) spectroscopy has a long history of elucidating chemical structures from experimentally measured 1H and 13C chemical shifts. One approach to characterizing the chemical structures of an unknown metabolite is to predict the 1H and 13C chemical shifts of candidate compounds (e.g., metabolites from the Human Metabolome Database (HMDB)) and compare them with chemical shifts of the unknown. However, accurate prediction of NMR chemical shifts in aqueous solution is challenging due to limitations of experimental chemical shift libraries and the high computational cost of quantum chemical methods. To improve NMR prediction accuracy and applicability, an empirical prediction strategy is introduced here to provide an accurately predicted chemical shift for organic molecules and metabolites within seconds. Unique features of COLMARppm include (i) the training library exclusively consisting of high quality NMR spectra measured under standard conditions in aqueous solution, (ii) utilization of NMR motif information, and (iii) leveraging of the improved prediction accuracy for the automated assignment of experimental chemical shifts for candidate structures. COLMARppm is demonstrated in terms of accuracy and speed for a set of 20 compounds taken from the HMDB for chemical shift prediction and resonance assignment. COLMARppm is applicable to a wide range of small molecules and can be directly incorporated into metabolomics workflows.

Nanoconfined Electrokinetic Chromatography (NEC): Gradient Separation and Sensing of Short DNA Fragments at the Single-Molecule Level.

Glass nanopipets have been demonstrated to be a powerful tool for the sensing and discrimination of biomolecules, such as DNA strands with different lengths or configurations. Despite progress made in nanopipet-based sensors, it remains challenging to develop effective strategies that separate and sense in one operation. In this study, we demonstrate an agarose gel-filled nanopipet that enables hyphenated length-dependent separation and electrochemical sensing of short DNA fragments based on the electrokinetic flow of DNA molecules in the nanoconfined channel at the tip of the nanopipet. This nanoconfined electrokinetic chromatography (NEC) method is used to distinguish the mixture of DNA strands without labels, and the ionic current signals measured in real time show that the mixed DNA strands pass through the tip hole in order according to the molecular weight. With NEC, gradient separation and electrochemical measurement of biomolecules can be achieved simultaneously at the single-molecule level, which is further applied for programmable gene delivery into single living cells. Overall, NEC provides a multipurpose platform integrating separation, sensing, single-cell delivery, and manipulation, which may bring new insights into advanced bioapplication.

High-Resolution and Dynamic Visualization of Intracellular Redox Potential Using a Metal-Organic Framework-Functionalized Nanopotentiometer.

Redox potential plays a key role in regulating intracellular signaling pathways, with its quantitative analysis in individual cells benefiting our understanding of the underlying mechanism in the pathophysiological events. Here, a metal organic framework (MOF)-functionalized SERS nanopotentiometer has been developed for the dynamic monitoring of intracellular redox potential. The approach is based on the encapsulation of zirconium-based MOF (Uio-66-F4) on a surface of gold-silver nanorods (Au-Ag NRs) that is modified with the newly synthesized redox-sensitive probe ortho-mercaptohydroquinone (HQ). Thanks to size exclusion of MOF as the chemical protector, the nanopotentiometer can be adapted to long-term use and possess high anti-interference ability toward nonredox species. Combining the superior fingerprint identification of SERS with the electrochemical activity of the quinone/hydroquinone, the nanopotentiometer shows a reversible redox responsivity and can quantify redox potential with a relatively wide range of -250-100 mV. Furthermore, the nanopotentiometer allows for dynamic visualization of intracellular redox potential changes induced by drugs' stimulation in a high-resolution manner. The developed approach would be promising for offering new insights into the correlation between redox potential and tumor proliferation-involved processes such as oxidative stress and hypoxia.

Holistic Prediction of AuNP Aggregation in Diverse Aqueous Suspensions Based on Machine Vision and Dark-Field Scattering Imaging.

The localized surface-plasmon resonance of the AuNP in aqueous media is extremely sensitive to environmental changes. By measuring the signal of plasmon scattering light, the dark-field microscopic (DFM) imaging technique has been used to monitor the aggregation of AuNPs, which has attracted great attention because of its simplicity, low cost, high sensitivity, and universal applicability. However, it is still challenging to interpret DFM images of AuNP aggregation due to the heterogeneous characteristics of the isolated and discontinuous color distribution. Herein, we introduce machine vision algorithms for the training of DFM images of AuNPs in different saline aqueous media. A visual deep learning framework based on AlexNet is constructed for studying the aggregation patterns of AuNPs in aqueous suspensions, which allows for rapid and accurate identification of the aggregation extent of AuNPs, with a prediction accuracy higher than 0.96. With the aid of machine learning analysis, we further demonstrate the prediction ability of various aggregation phenomena induced by both cation species and the concentration of the external saline solution. Our results suggest the great potential of machine vision frameworks in the accurate recognition of subtle pattern changes in DFM images, which can help researchers build predictive analytics based on DFM imaging data.

A multi-dimensional approach to unravel the intricacies of lactylation related signature for prognostic and therapeutic insight in colorectal cancer.

BACKGROUND: Lactylation, a novel contributor to post-translational protein modifications, exhibits dysregulation across various tumors. Nevertheless, its intricate involvement in colorectal carcinoma, particularly for non-histone lactylation and its intersection with metabolism and immune evasion, remains enigmatic. METHODS: Employing immunohistochemistry on tissue microarray with clinical information and immunofluorescence on colorectal cell lines, we investigated the presence of global lactylation and its association with development and progression in colorectal cancer as well as its functional location. Leveraging the AUCell algorithm alongside correlation analysis in single-cell RNA sequencing data, as well as cox-regression and lasso-regression analysis in TCGA dataset and confirmed in GEO dataset, we identified a 23-gene signature predicting colorectal cancer prognosis. Subsequently, we analyzed the associations between the lactylation related gene risk and clinical characteristics, mutation landscapes, biological functions, immune cell infiltration, immunotherapy responses, and drug sensitivity. Core genes were further explored for deep biological insights through bioinformatics and in vitro experiments. RESULTS: Our study innovatively reveals a significant elevation of global lactylation in colorectal cancer, particularly in malignant tumors, confirming it as an independent prognostic factor for CRC. Through a comprehensive analysis integrating tumor tissue arrays, TCGA dataset, GEO dataset, combining in silico investigations and in vitro experiments, we identified a 23-gene Lactylation-Related Gene risk model capable of predicting the prognosis of colorectal cancer patients. Noteworthy variations were observed in clinical characteristics, biological functions, immune cell infiltration, immune checkpoint expression, immunotherapy responses and drug sensitivity among distinct risk groups. CONCLUSIONS: The Lactylation-Related Gene risk model exhibits significant potential for improving the management of colorectal cancer patients and enhancing therapeutic outcomes, particularly at the intersection of metabolism and immune evasion. This finding underscores the clinical relevance of global lactylation in CRC and lays the groundwork for mechanism investigation and targeted therapeutic strategies given the high lactate concentration in CRC.

Dynamic m6A mRNA methylation reveals the involvement of AcALKBH10 in ripening-related quality regulation in kiwifruit.

Kiwifruit ripening is a complex and highly coordinated process that occurs in conjunction with the formation of fruit edible quality. The significance of epigenetic changes, particularly the impact of N6-methyladenosine (m6A) RNA modification on fruit ripening and quality formation, has been largely overlooked. We monitored m6A levels and gene expression changes in kiwifruit at four different stages using LC-MS/MS, MeRIP, RNA-seq, and validated the function of AcALKBH10 through heterologous transgenic expression in tomato. Notable m6A modifications occurred predominantly at the stop codons and the 3' UTRs and exhibited a gradual reduction in m6A levels during the fruit ripening process. Moreover, these m6A modifications in the aforementioned sites demonstrated a discernible inverse relationship with the levels of mRNA abundance throughout the ripening process, suggesting a repression effect of m6A modification in the modulation of kiwifruit ripening. We further demonstrated that AcALKBH10 rather than AcECT9 predominantly regulates m6A levels in ripening-related genes, thereby exerting the regulatory control over the ripening process and the accumulation of soluble sugars and organic acids, ultimately influencing fruit ripening and quality formation. In conclusion, our findings illuminate the epi-regulatory mechanism involving m6A in kiwifruit ripening, offering a fresh perspective for cultivating high-quality kiwifruit with enhanced nutritional attributes.

Kiwifruit bZIP transcription factor AcePosF21 elicits ascorbic acid biosynthesis during cold stress.

Cold stress seriously affects plant development, resulting in heavy agricultural losses. L-ascorbic acid (AsA, vitamin C) is an antioxidant implicated in abiotic stress tolerance and metabolism of reactive oxygen species (ROS). Understanding whether and how cold stress elicits AsA biosynthesis to reduce oxidative damage is important for developing cold-resistant plants. Here, we show that the accumulation of AsA in response to cold stress is a common mechanism conserved across the plant kingdom, from single-cell algae to angiosperms. We identified a basic leucine zipper domain (bZIP) transcription factor (TF) of kiwifruit (Actinidia eriantha Benth.), AcePosF21, which was triggered by cold and is involved in the regulation of kiwifruit AsA biosynthesis and defense responses against cold stress. AcePosF21 interacted with the R2R3-MYB TF AceMYB102 and directly bound to the promoter of the gene encoding GDP-L-galactose phosphorylase 3 (AceGGP3), a key conduit for regulating AsA biosynthesis, to up-regulate AceGGP3 expression and produce more AsA, which neutralized the excess ROS induced by cold stress. On the contrary, VIGS or CRISPR-Cas9-mediated editing of AcePosF21 decreased AsA content and increased the generation of ROS in kiwifruit under cold stress. Taken together, we illustrated a model for the regulatory mechanism of AcePosF21-mediated regulation of AceGGP3 expression and AsA biosynthesis to reduce oxidative damage by cold stress, which provides valuable clues for manipulating the cold resistance of kiwifruit.