A Thermostable mRNA Vaccine against COVID-19.

There is an urgent need for vaccines against coronavirus disease 2019 (COVID-19) because of the ongoing SARS-CoV-2 pandemic. Among all approaches, a messenger RNA (mRNA)-based vaccine has emerged as a rapid and versatile platform to quickly respond to this challenge. Here, we developed a lipid nanoparticle-encapsulated mRNA (mRNA-LNP) encoding the receptor binding domain (RBD) of SARS-CoV-2 as a vaccine candidate (called ARCoV). Intramuscular immunization of ARCoV mRNA-LNP elicited robust neutralizing antibodies against SARS-CoV-2 as well as a Th1-biased cellular response in mice and non-human primates. Two doses of ARCoV immunization in mice conferred complete protection against the challenge of a SARS-CoV-2 mouse-adapted strain. Additionally, ARCoV is manufactured as a liquid formulation and can be stored at room temperature for at least 1 week. ARCoV is currently being evaluated in phase 1 clinical trials.

Salivary factor LTRIN from Aedes aegypti facilitates the transmission of Zika virus by interfering with the lymphotoxin-ß receptor.

Pathogens have co-evolved with mosquitoes to optimize transmission to hosts. Mosquito salivary-gland extract is known to modulate host immune responses and facilitate pathogen transmission, but the underlying molecular mechanisms of this have remained unknown. In this study, we identified and characterized a prominent 15-kilodalton protein, LTRIN, obtained from the salivary glands of the mosquito Aedes aegypti. LTRIN expression was upregulated in blood-fed mosquitoes, and LTRIN facilitated the transmission of Zika virus (ZIKV) and exacerbated its pathogenicity by interfering with signaling through the lymphotoxin-ß receptor (LTßR). Mechanically, LTRIN bound to LTßR and 'preferentially' inhibited signaling via the transcription factor NF-κB and the production of inflammatory cytokines by interfering with the dimerization of LTßR during infection with ZIKV. Furthermore, treatment with antibody to LTRIN inhibited mosquito-mediated infection with ZIKV, and abolishing LTßR potentiated the infectivity of ZIKV both in vitro and in vivo. This study provides deeper insight into the transmission of mosquito-borne diseases in nature and supports the therapeutic potential of inhibiting the action of LTRIN to disrupt ZIKV transmission.

Detection of SARS-CoV-2-Specific Humoral and Cellular Immunity in COVID-19 Convalescent Individuals.

The World Health Organization has declared SARS-CoV-2 virus outbreak a worldwide pandemic. However, there is very limited understanding on the immune responses, especially adaptive immune responses to SARS-CoV-2 infection. Here, we collected blood from COVID-19 patients who have recently become virus-free, and therefore were discharged, and detected SARS-CoV-2-specific humoral and cellular immunity in eight newly discharged patients. Follow-up analysis on another cohort of six patients 2 weeks post discharge also revealed high titers of immunoglobulin G (IgG) antibodies. In all 14 patients tested, 13 displayed serum-neutralizing activities in a pseudotype entry assay. Notably, there was a strong correlation between neutralization antibody titers and the numbers of virus-specific T cells. Our work provides a basis for further analysis of protective immunity to SARS-CoV-2, and understanding the pathogenesis of COVID-19, especially in the severe cases. It also has implications in developing an effective vaccine to SARS-CoV-2 infection.

Memory B cell repertoire from triple vaccinees against diverse SARS-CoV-2 variants.

Omicron (B.1.1.529), the most heavily mutated SARS-CoV-2 variant so far, is highly resistant to neutralizing antibodies, raising concerns about the effectiveness of antibody therapies and vaccines1,2. Here we examined whether sera from individuals who received two or three doses of inactivated SARS-CoV-2 vaccine could neutralize authentic Omicron. The seroconversion rates of neutralizing antibodies were 3.3% (2 out of 60) and 95% (57 out of 60) for individuals who had received 2 and 3 doses of vaccine, respectively. For recipients of three vaccine doses, the geometric mean neutralization antibody titre for Omicron was 16.5-fold lower than for the ancestral virus (254). We isolated 323 human monoclonal antibodies derived from memory B cells in triple vaccinees, half of which recognized the receptor-binding domain, and showed that a subset (24 out of 163) potently neutralized all SARS-CoV-2 variants of concern, including Omicron. Therapeutic treatments with representative broadly neutralizing monoclonal antibodies were highly protective against infection of mice with SARS-CoV-2 Beta (B.1.351) and Omicron. Atomic structures of the Omicron spike protein in complex with three classes of antibodies that were active against all five variants of concern defined the binding and neutralizing determinants and revealed a key antibody escape site, G446S, that confers greater resistance to a class of antibodies that bind on the right shoulder of the receptor-binding domain by altering local conformation at the binding interface. Our results rationalize the use of three-dose immunization regimens and suggest that the fundamental epitopes revealed by these broadly ultrapotent antibodies are rational targets for a universal sarbecovirus vaccine.

scGIR: deciphering cellular heterogeneity via gene ranking in single-cell weighted gene correlation networks.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular heterogeneity through high-throughput analysis of individual cells. Nevertheless, challenges arise from prevalent sequencing dropout events and noise effects, impacting subsequent analyses. Here, we introduce a novel algorithm, Single-cell Gene Importance Ranking (scGIR), which utilizes a single-cell gene correlation network to evaluate gene importance. The algorithm transforms single-cell sequencing data into a robust gene correlation network through statistical independence, with correlation edges weighted by gene expression levels. We then constructed a random walk model on the resulting weighted gene correlation network to rank the importance of genes. Our analysis of gene importance using PageRank algorithm across nine authentic scRNA-seq datasets indicates that scGIR can effectively surmount technical noise, enabling the identification of cell types and inference of developmental trajectories. We demonstrated that the edges of gene correlation, weighted by expression, play a critical role in enhancing the algorithm's performance. Our findings emphasize that scGIR outperforms in enhancing the clustering of cell subtypes, reverse identifying differentially expressed marker genes, and uncovering genes with potential differential importance. Overall, we proposed a promising method capable of extracting more information from single-cell RNA sequencing datasets, potentially shedding new lights on cellular processes and disease mechanisms.

Human Virus-Derived Small RNAs Can Confer Antiviral Immunity in Mammals.

RNA interference (RNAi) functions as a potent antiviral immunity in plants and invertebrates; however, whether RNAi plays antiviral roles in mammals remains unclear. Here, using human enterovirus 71 (HEV71) as a model, we showed HEV71 3A protein as an authentic viral suppressor of RNAi during viral infection. When the 3A-mediated RNAi suppression was impaired, the mutant HEV71 readily triggered the production of abundant HEV71-derived small RNAs with canonical siRNA properties in cells and mice. These virus-derived siRNAs were produced from viral dsRNA replicative intermediates in a Dicer-dependent manner and loaded into AGO, and they were fully active in degrading cognate viral RNAs. Recombinant HEV71 deficient in 3A-mediated RNAi suppression was significantly restricted in human somatic cells and mice, whereas Dicer deficiency rescued HEV71 infection independently of type I interferon response. Thus, RNAi can function as an antiviral immunity, which is induced and suppressed by a human virus, in mammals.

25-Hydroxycholesterol Protects Host against Zika Virus Infection and Its Associated Microcephaly in a Mouse Model.

Zika virus (ZIKV) has become a public health threat due to its global transmission and link to severe congenital disorders. The host immune responses to ZIKV infection have not been fully elucidated, and effective therapeutics are not currently available. Herein, we demonstrated that cholesterol-25-hydroxylase (CH25H) was induced in response to ZIKV infection and that its enzymatic product, 25-hydroxycholesterol (25HC), was a critical mediator of host protection against ZIKV. Synthetic 25HC addition inhibited ZIKV infection in vitro by blocking viral entry, and treatment with 25HC reduced viremia and conferred protection against ZIKV in mice and rhesus macaques. 25HC suppressed ZIKV infection and reduced tissue damage in human cortical organoids and the embryonic brain of the ZIKV-induced mouse microcephaly model. Our findings highlight the protective role of CH25H during ZIKV infection and the potential use of 25HC as a natural antiviral agent to combat ZIKV infection and prevent ZIKV-associated outcomes, such as microcephaly.

The lack of negative association between TE load and subgenome dominance in synthesized Brassica allotetraploids.

Polyploidization is important to the evolution of plants. Subgenome dominance is a distinct phenomenon associated with most allopolyploids. A gene on the dominant subgenome tends to express to higher RNA levels in all organs as compared to the expression of its syntenic paralogue (homoeolog). The mechanism that underlies the formation of subgenome dominance remains unknown, but there is evidence for the involvement of transposon/DNA methylation density differences nearby the genes of parents as being causal. The subgenome with lower density of transposon and methylation near genes is positively associated with subgenome dominance. Here, we generated eight generations of allotetraploid progenies from the merging of parental genomes Brassica rapa and Brassica oleracea. We found that transposon/methylation density differ near genes between the parental (rapa:oleracea) existed in the wide hybrid, persisted in the neotetraploids (the synthetic Brassica napus), but these neotetraploids expressed no expected subgenome dominance. This absence of B. rapa vs. B. oleracea subgenome dominance is particularly significant because, while there is no negative relationship between transposon/methylation level and subgenome dominance in the neotetraploids, the more ancient parental subgenomes for all Brassica did show differences in transposon/methylation densities near genes and did express, in the same samples of cells, biased gene expression diagnostic of subgenome dominance. We conclude that subgenome differences in methylated transposon near genes are not sufficient to initiate the biased gene expressions defining subgenome dominance. Our result was unexpected, and we suggest a "nuclear chimera" model to explain our data.

Deep learning can predict subgenome dominance in ancient but not in neo/synthetic polyploidized genomes.

Deep learning offers new approaches to investigate the mechanisms underlying complex biological phenomena, such as subgenome dominance. Subgenome dominance refers to the dominant expression and/or biased fractionation of genes in one subgenome of allopolyploids, which has shaped the evolution of a large group of plants. However, the underlying cause of subgenome dominance remains elusive. Here, we adopt deep learning to construct two convolutional neural network (CNN) models, binary expression model (BEM) and homoeolog contrast model (HCM), to investigate the mechanism underlying subgenome dominance using DNA sequence and methylation sites. We apply these CNN models to analyze three representative polyploidization systems, Brassica, Gossypium, and Cucurbitaceae, each with available ancient and neo/synthetic polyploidized genomes. The BEM shows that DNA sequence of the promoter region can accurately predict whether a gene is expressed or not. More importantly, the HCM shows that the DNA sequence of the promoter region predicts dominant expression status between homoeologous gene pairs retained from ancient polyploidizations, thus predicting subgenome dominance associated with these events. However, HCM fails to predict gene expression dominance between new homoeologous gene pairs arising from the neo/synthetic polyploidizations. These results are consistent across the three plant polyploidization systems, indicating broad applicability of our models. Furthermore, the two models based on methylation sites produce similar results. These results show that subgenome dominance is associated with long-term sequence differentiation between the promoters of homoeologs, suggesting that subgenome expression dominance precedes and is the driving force or even the determining factor for sequence divergence between subgenomes following polyploidization.

Chromosome-level genome of putative autohexaploid Actinidia deliciosa provides insights into polyploidisation and evolution.

Actinidia ('Mihoutao' in Chinese) includes species with complex ploidy, among which diploid Actinidia chinensis and hexaploid Actinidia deliciosa are economically and nutritionally important fruit crops. Actinidia deliciosa has been proposed to be an autohexaploid (2n = 174) with diploid A. chinensis (2n = 58) as the putative parent. A CCS-based assembly anchored to a high-resolution linkage map provided a chromosome-resolved genome for hexaploid A. deliciosa yielded a 3.91-Gb assembly of 174 pseudochromosomes comprising 29 homologous groups with 6 members each, which contain 39 854 genes with an average of 4.57 alleles per gene. Here we provide evidence that much of the hexaploid genome matches diploid A. chinensis; 95.5% of homologous gene pairs exhibited >90% similarity. However, intragenome and intergenome comparisons of synteny indicate chromosomal changes. Our data, therefore, indicate that if A. deliciosa is an autoploid, chromosomal rearrangement occurred following autohexaploidy. A highly diversified pattern of gene expression and a history of rapid population expansion after polyploidisation likely facilitated the adaptation and niche differentiation of A. deliciosa in nature. The allele-defined hexaploid genome of A. deliciosa provides new genomic resources to accelerate crop improvement and to understand polyploid genome evolution.

Gene and Allele-Specific Expression Underlying the Electric Signal Divergence in African Weakly Electric Fish.

In the African weakly electric fish genus Campylomormyrus, electric organ discharge signals are strikingly different in shape and duration among closely related species, contribute to prezygotic isolation, and may have triggered an adaptive radiation. We performed mRNA sequencing on electric organs and skeletal muscles (from which the electric organs derive) from 3 species with short (0.4 ms), medium (5 ms), and long (40 ms) electric organ discharges and 2 different cross-species hybrids. We identified 1,444 upregulated genes in electric organ shared by all 5 species/hybrid cohorts, rendering them candidate genes for electric organ-specific properties in Campylomormyrus. We further identified several candidate genes, including KCNJ2 and KLF5, and their upregulation may contribute to increased electric organ discharge duration. Hybrids between a short (Campylomormyrus compressirostris) and a long (Campylomormyrus rhynchophorus) discharging species exhibit electric organ discharges of intermediate duration and showed imbalanced expression of KCNJ2 alleles, pointing toward a cis-regulatory difference at this locus, relative to electric organ discharge duration. KLF5 is a transcription factor potentially balancing potassium channel gene expression, a crucial process for the formation of an electric organ discharge. Unraveling the genetic basis of the species-specific modulation of the electric organ discharge in Campylomormyrus is crucial for understanding the adaptive radiation of this emerging model taxon of ecological (perhaps even sympatric) speciation.

The subgenomic flaviviral RNA suppresses RNA interference through competing with siRNAs for binding RISC components.

During the life cycle of mosquito-borne flaviviruses, substantial subgenomic flaviviral RNA (sfRNA) is produced via incomplete degradation of viral genomic RNA by host XRN1. Zika virus (ZIKV) sfRNA has been detected in mosquito and mammalian somatic cells. Human neural progenitor cells (hNPCs) in the developing brain are the major target cells of ZIKV, and antiviral RNA interference (RNAi) plays a critical role in hNPCs. However, whether ZIKV sfRNA was produced in ZIKV-infected hNPCs as well as its function remains not known. In this study, we demonstrate that abundant sfRNA was produced in ZIKV-infected hNPCs. RNA pulldown and mass spectrum assays showed ZIKV sfRNA interacted with host proteins RHA and PACT, both of which are RNA-induced silencing complex (RISC) components. Functionally, ZIKV sfRNA can antagonize RNAi by outcompeting small interfering RNAs (siRNAs) in binding to RHA and PACT. Furthermore, the 3' stem loop (3'SL) of sfRNA was responsible for RISC components binding and RNAi inhibition, and 3'SL can enhance the replication of a viral suppressor of RNAi (VSR)-deficient virus in a RHA- and PACT-dependent manner. More importantly, the ability of binding to RISC components is conversed among multiple flaviviral 3'SLs. Together, our results identified flavivirus 3'SL as a potent VSR in RNA format, highlighting the complexity in virus-host interaction during flavivirus infection.IMPORTANCEZika virus (ZIKV) infection mainly targets human neural progenitor cells (hNPCs) and induces cell death and dysregulated cell-cycle progression, leading to microcephaly and other central nervous system abnormalities. RNA interference (RNAi) plays critical roles during ZIKV infections in hNPCs, and ZIKV has evolved to encode specific viral proteins to antagonize RNAi. Herein, we first show that abundant sfRNA was produced in ZIKV-infected hNPCs in a similar pattern to that in other cells. Importantly, ZIKV sfRNA acts as a potent viral suppressor of RNAi (VSR) by competing with siRNAs for binding RISC components, RHA and PACT. The 3'SL of sfRNA is responsible for binding RISC components, which is a conserved feature among mosquito-borne flaviviruses. As most known VSRs are viral proteins, our findings highlight the importance of viral non-coding RNAs during the antagonism of host RNAi-based antiviral innate immunity.

Harnessing ZIKV NS2A RNA for alleviating acute hepatitis and cytokine release storm by targeting translation machinery.

BACKGROUND AND AIMS: Hyperactivated inflammatory responses induced by cytokine release syndrome (CRS) are the primary causes of tissue damage and even death. The translation process is precisely regulated to control the production of proinflammatory cytokines. However, it is largely unknown whether targeting translation can effectively limit the hyperactivated inflammatory responses during acute hepatitis and graft-versus-host disease (GVHD). APPROACH AND RESULTS: By using in vitro translation and cellular overexpression systems, we have found that the non-structural protein gene NS2A of Zika virus (ZIKV) functions as RNA molecules to suppress the translation of both ectopic genes and endogenous proinflammatory cytokines. Mechanistically, results from RNA pulldown and co-immunoprecipitation (Co-IP) assays have demonstrated that NS2A RNA interacts with the translation initiation factor eIF2α to disrupt the dynamic balance of the eIF2/eIF2B complex and translation initiation, which is the rate-limiting step during the translation process. In the acetaminophen (APAP)-, LPS/D-galactosamine (D-GalN)-, viral infection-induced acute hepatitis, and GVHD mouse models, mice with myeloid cell-specific knock-in of NS2A show decreased levels of serum proinflammatory cytokines and reduced tissue damage. CONCLUSIONS: ZIKV NS2A dampens the production of proinflammatory cytokines and alleviates inflammatory injuries by interfering translation process as RNA molecules, which suggests that NS2A RNA is potentially used to treat numerous acute inflammatory diseases characterized by CRS.

Regional Active Transcription Associates with Homoeologous Exchange Breakpoints in Synthetic Brassica Tetraploids.

Polyploidization plays a crucial role in plant evolution and is becoming increasingly important in breeding. Structural variations and epigenomic repatterning have been observed in synthetic polyploidizations. However, the mechanisms underlying the occurrence and their effects on gene expression and phenotype remain unknown. Here, we investigated genome-wide large deletion/duplication regions (DelDups) and genomic methylation dynamics in leaf organs of progeny from the first eight generations of synthetic tetraploids derived from Chinese cabbage (Brassica rapa L. ssp. pekinensis) and cabbage (Brassica oleracea L. var. capitata). One- or two-copy DelDups, with a mean size of 5.70 Mb (400 kb - 65.85 Mb), occurred from the first generation of selfing and thereafter. The duplication of a fragment in one subgenome consistently coincided with the deletion of its syntenic fragment in the other subgenome, and vice versa, indicating that these DelDups were generated by homoeologous exchanges (HEs). Interestingly, the larger the genomic syntenic region, the higher the frequency of DelDups, further suggesting that the pairing of large homoeologous fragments is crucial for HEs. Moreover, we found that the active transcription of continuously distributed genes in local regions is positively associated with the occurrence of HE breakpoints. In addition, the expression of genes within DelDups exhibited a dosage effect, and plants with extra parental genomic fragments generally displayed phenotypes biased towards the corresponding parent. Genome-wide methylation fluctuated remarkably, which did not clearly affect gene expression on a large scale. Our findings provide insights into the early evolution of polyploid genomes, offering valuable knowledge for polyploidization-based breeding.

Dedifferentiation-like reprogramming of degenerative nucleus pulposus cells into notochordal-like cells by defined factors.

The extensive degeneration of functional somatic cells and the depletion of endogenous stem/progenitor populations present significant challenges to tissue regeneration in degenerative diseases. Currently, a cellular reprogramming approach enabling directly generating corresponding progenitor populations from degenerative somatic cells remains elusive. The present study focused on intervertebral disc degeneration (IVDD) and identified a three-factor combination (OCT4, FOXA2, TBXT [OFT]) that could induce the dedifferentiation-like reprogramming of degenerative nucleus pulposus cells (dNPCs) toward induced notochordal-like cells (iNCs). Single-cell transcriptomics dissected the transitions of cell identity during reprogramming. Further, OCT4 was found to directly interact with bromodomain PHD-finger transcription factor to remodel the chromatin during the early phases, which was crucial for initiating this dedifferentiation-like reprogramming. In rat models, intradiscal injection of adeno-associated virus carrying OFT generated iNCs from in situ dNPCs and reversed IVDD. These results collectively present a proof-of-concept for dedifferentiation-like reprogramming of degenerated somatic cells into corresponding progenitors through the development of a factor-based strategy, providing a promising approach for regeneration in degenerative disc diseases.

The meso-octoploid Heliophila variabilis genome sheds a new light on the impact of polyploidization and diploidization on the diversity of the Cape flora.

Although the South African Cape flora is one of the most remarkable biodiversity hotspots, its high diversity has not been associated with polyploidy. Here, we report the chromosome-scale genome assembly of an ephemeral cruciferous species Heliophila variabilis (~334 Mb, n = 11) adapted to South African semiarid biomes. Two pairs of differently fractionated subgenomes suggest an allo-octoploid origin of the genome at least 12 million years ago. The ancestral octoploid Heliophila genome (2n = 8x = ~60) has probably originated through hybridization between two allotetraploids (2n = 4x = ~30) formed by distant, intertribal, hybridization. Rediploidization of the ancestral genome was marked by extensive reorganization of parental subgenomes, genome downsizing, and speciation events in the genus Heliophila. We found evidence for loss-of-function changes in genes associated with leaf development and early flowering, and over-retention and sub/neofunctionalization of genes involved in pathogen response and chemical defense. The genomic resources of H. variabilis will help elucidate the role of polyploidization and genome diploidization in plant adaptation to hot arid environments and origin of the Cape flora. The sequenced H. variabilis represents the first chromosome-scale genome assembly of a meso-octoploid representative of the mustard family.

ZCCHC8 p.P410A disrupts nucleocytoplasmic localization, promoting idiopathic pulmonary fibrosis and chronic obstructive pulmonary disease.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is a special kind of chronic interstitial lung disease with insidious onset. Previous studies have revealed that mutations in ZCCHC8 may lead to IPF. The aim of this study is to explore the ZCCHC8 mutations in Chinese IPF patients. METHODS: Here, we enrolled 124 patients with interstitial lung disease from 2017 to 2023 in our hospital. Whole exome sequencing and Sanger sequencing were employed to explore the genetic lesions of these patients. RESULTS: Among these 124 patients, a novel mutation (NM_017612: c.1228 C > G/p.P410A) of Zinc Finger CCHC-Type Containing 8 (ZCCHC8)was identified in a family with IPF and chronic obstructive lung disease. As a component of the nuclear exosome-targeting complex that regulates the turnover of human telomerase RNA, ZCCHC8 mutations have been reported may lead to IPF in European population and American population. Functional study confirmed that the novel mutation can disrupt the nucleocytoplasmic localization of ZCCHC8, which further decreased the expression of DKC1 and RTEL1, and finally reduced the length of telomere and led to IPF and related disorders. CONCLUSIONS: We may first report the ZCCHC8 mutation in Asian population with IPF. Our study broadens the mutation, phenotype, and population spectrum of ZCCHC8 deficiency.

Evaluation of reverse transcription yield of RNA standards and forensic samples based on droplet digital PCR.

RNA analysis has shown great value in forensic science, such as body fluids and tissue identification, postmortem interval estimation, biological age prediction, etc. Currently, most RNA follow-up experiments involve reverse transcription (RT) procedures. It has been shown that the RT step is variable and has a greater impact on subsequent data analysis, especially for forensic trace samples. However, the pattern of variation between different RNA template inputs and complementary DNA (cDNA) yield is unclear. In this study, a series of 2-fold gradient dilutions of RNA standards (1 µg/µL - 0.24 ng/µL) and forensic samples (including blood samples, saliva samples, bloodstains, and saliva stains) were reverse-transcribed using EasyQuick RT MasterMix. The obtained cDNA was quantified by droplet digital PCR (ddPCR) to assess the RT yield of the ACTB gene. The results showed that the 125 ng RNA template had the highest RT yield in a 10 µL RT reaction system with the selected kit. For all stain samples, the RT yield improved as the amount of RNA template input increased since RNA quantities were below 125 ng. As many commercialized reverse transcription kits using different kinds of enzymes are available for forensic RNA research, we recommend that systematic experiments should be performed in advance to determine the amount of RNA input at the optimum RT yield when using any kit for reverse transcription experiments.

Nitric oxide mediates positive regulation of Nostoc flagelliforme polysaccharide yield via potential S-nitrosylation of G6PDH and UGDH.

Based on our previous findings that salicylic acid and jasmonic acid increased Nostoc flagelliforme polysaccharide yield by regulating intracellular nitric oxide (NO) levels, the mechanism through which NO affects polysaccharide biosynthesis in Nostoc flagelliforme was explored from the perspective of S-nitrosylation (SNO). The addition of NO donor and scavenger showed that intracellular NO had a significant positive effect on the polysaccharide yield of N. flagelliforme. To explore the mechanism, we investigated the relationship between NO levels and the activity of several key enzymes involved in polysaccharide biosynthesis, including fructose 1,6-bisphosphate aldolase (FBA), glucokinase (GK), glucose 6-phosphate dehydrogenase (G6PDH), mitochondrial isocitrate dehydrogenase (ICDH), and UDP-glucose dehydrogenase (UGDH). The enzymatic activities of G6PDH, ICDH, and UGDH were shown to be significantly correlated with the shifts in intracellular NO levels. For further validation, G6PDH, ICDH, and UGDH were heterologously expressed in Escherichia coli and purified via Ni+-NAT affinity chromatography, and subjected to a biotin switch assay and western blot analysis, which revealed that UGDH and G6PDH were susceptible to SNO. Furthermore, mass spectrometry analysis of proteins treated with S-nitrosoglutathione (GSNO) identified the SNO modification sites for UGDH and G6PDH as cysteine 423 and cysteine 249, respectively. These findings suggest that NO modulates polysaccharide biosynthesis in N. flagelliforme through SNO of UGDH and G6PDH. This reveals a potential mechanism through which NO promotes polysaccharide synthesis in N. flagelliforme, while also providing a new strategy for improving the industrial production of polysaccharides.

Pan-genome analysis of 13 Spinacia accessions reveals structural variations associated with sex chromosome evolution and domestication traits in spinach.

Structural variations (SVs) are major genetic variants that can be involved in the origin, adaptation and domestication of species. However, the identification and characterization of SVs in Spinacia species are rare due to the lack of a pan-genome. Here, we report eight chromosome-scale assemblies of cultivated spinach and its two wild species. After integration with five existing assemblies, we constructed a comprehensive Spinacia pan-genome and identified 193 661 pan-SVs, which were genotyped in 452 Spinacia accessions. Our pan-SVs enabled genome-wide association study identified signals associated with sex and clarified the evolutionary direction of spinach. Most sex-linked SVs (86%) were biased to occur on the Y chromosome during the evolution of the sex-linked region, resulting in reduced Y-linked gene expression. The frequency of pan-SVs among Spinacia accessions further illustrated the contribution of these SVs to domestication, such as bolting time and seed dormancy. Furthermore, compared with SNPs, pan-SVs act as efficient variants in genomic selection (GS) because of their ability to capture missing heritability information and higher prediction accuracy. Overall, this study provides a valuable resource for spinach genomics and highlights the potential utility of pan-SV in crop improvement and breeding programmes.