Citrus Huanglongbing is a pathogen-triggered immune disease that can be mitigated with antioxidants and gibberellin

Ma, Wenxiu; Pang, Zhiqian; Huang, Xiaoen; Xu, Jin; Pandey, Sheo Shankar; Li, Jinyun; Achor, Diann S.; Vasconcelos, Fernanda N. C.; Hendrich, Connor; Huang, Yixiao; Wang, Wenting; Lee, Donghwan; Stanton, Daniel; Wang, Nian

doi:10.1038/s41467-022-28189-9

Download PDF

Article
Open access
Published: 26 January 2022

Citrus Huanglongbing is a pathogen-triggered immune disease that can be mitigated with antioxidants and gibberellin

Nature Communications volume 13, Article number: 529 (2022) Cite this article

25k Accesses
58 Citations
39 Altmetric
Metrics details

Subjects

Abstract

Huanglongbing (HLB) is a devastating disease of citrus, caused by the phloem-colonizing bacterium Candidatus Liberibacter asiaticus (CLas). Here, we present evidence that HLB is an immune-mediated disease. We show that CLas infection of Citrus sinensis stimulates systemic and chronic immune responses in phloem tissue, including callose deposition, production of reactive oxygen species (ROS) such as H₂O₂, and induction of immunity-related genes. The infection also upregulates genes encoding ROS-producing NADPH oxidases, and downregulates antioxidant enzyme genes, supporting that CLas causes oxidative stress. CLas-triggered ROS production localizes in phloem-enriched bark tissue and is followed by systemic cell death of companion and sieve element cells. Inhibition of ROS levels in CLas-positive stems by NADPH oxidase inhibitor diphenyleneiodonium (DPI) indicates that NADPH oxidases contribute to CLas-triggered ROS production. To investigate potential treatments, we show that addition of the growth hormone gibberellin (known to have immunoregulatory activities) upregulates genes encoding H₂O₂-scavenging enzymes and downregulates NADPH oxidases. Furthermore, foliar spray of HLB-affected citrus with gibberellin or antioxidants (uric acid, rutin) reduces H₂O₂ concentrations and cell death in phloem tissues and reduces HLB symptoms. Thus, our results indicate that HLB is an immune-mediated disease that can be mitigated with antioxidants and gibberellin.

Legume rhizodeposition promotes nitrogen fixation by soil microbiota under crop diversification

Article Open access 04 April 2024

Mengjie Qiao, Ruibo Sun, … Yan Chen

A plant NLR receptor employs ABA central regulator PP2C-SnRK2 to activate antiviral immunity

Article Open access 13 April 2024

Shen Huang, Chunli Wang, … Xiaorong Tao

Mechanisms of calcium homeostasis orchestrate plant growth and immunity

Article 28 February 2024

Chao Wang, Ren-Jie Tang, … Sheng Luan

Introduction

Plants and animals utilize their immune system to fight pathogens. Both plants and animals have innate immunity, whereas animals also have adapted immunity. The plant innate immune system consists of pattern-triggered immunity (PTI), which is triggered by pathogen-associated molecular patterns (PAMPs) via cell surface-localized pattern-recognition receptors, and effector-triggered immunity (ETI), which is instigated by pathogen effector proteins via intracellular receptors called nucleotide-binding, leucine-rich repeat receptors (NLRs)^1,2,3. However, some human diseases are mediated by the immune response itself, including autoimmune diseases such as inflammatory bowel disease, and non-autoimmune diseases, like asthma, sepsis induced by various microbes⁴, and cryptococcal meningitis caused by Cryptococcus neoformans⁵. Host immune responses have been known to be an important factor in addition to microbial pathogenicity factors for human diseases caused by microbial pathogens, which has resulted in the proposal of a damage–response framework of microbial pathogenesis⁶. Immune-mediated diseases have not been reported in the Plantae Kingdom^7,8.

The damaging effect of plant diseases has been assumed to directly result from the impact of pathogenicity factors of the causal pathogens⁹. Common pathogenicity factors include effectors, toxins, cell wall degrading enzymes, and biofilm that are directly responsible for causing disease symptoms. For instance, the transcriptional activator-like effector PthA4 is responsible for the hypertrophy and hyperplasia symptoms of citrus canker caused by Xanthomonas citri subsp. citri¹⁰. Xylem blockage caused by biofilm of Xylella fastidiosa leads to the wilting of grapevine plants with Pierce’s disease¹¹. Citrus Huanglongbing (HLB, also known as citrus greening) is currently the most devastating citrus disease, causing billions of dollars of economic losses worldwide annually^12,13. HLB presents an unprecedented challenge for the citrus industry despite some promising progress in research^14,15,16. HLB is caused by the phloem-colonizing Candidatus Liberibacter asiaticus (CLas), Ca. L. americanus and Ca. L. africanus with CLas being the most prevalent in the world¹⁷. CLas is a biotroph, vectored by the Asian citrus psyllid (Diaphorina citri), and in the Rhizobiaceae family. Its ~1.23 Mb genome is significantly reduced compared to other free-living members of the Rhizobiaceae family, resulting from the reductive evolution within the nutrient-rich phloem of citrus and psyllid tissues¹⁸.

Despite its economic importance, how Ca. Liberibacter damages infected citrus plants remains poorly understood, partly because HLB pathogens have not been cultured in artificial media. No pathogenicity factors have been confirmed to be responsible for the HLB symptoms of characteristic blotchy mottle on leaves, hardened leaves, small and upright leaves, leaves showing zinc or manganese deficiency, corky veins, twig dieback, stunted growth of seedlings, thin canopy, small and lopsided fruit and root decay¹⁷.

In this work, we present evidence that citrus HLB is an immune-mediated disease. This hypothesis explains most HLB phenomena, is consistent with the genetic information of Ca. Liberibacter spp., and provides guidance regarding HLB management.

Results

CLas does not contain pathogenicity factors that directly cause HLB symptoms

We conducted a comprehensive analysis of CLas proteins and did not identify any homologs with known pathogenicity factors that are directly responsible for causing plant disease symptoms. Ca. Liberibacter spp. do not contain type II, III, or IV secretion systems that secrete such pathogenicity factors. The pathogenicity genes of closely related Agrobacterium and Rhizobium pathogens^19,20 were not identified in CLas. To test whether CLas contains any novel pathogenicity factors responsible for causing HLB symptoms, 47 putative virulence factors including serralysin and hemolysin (substrates of type I secretion system) and proteins containing Sec secretion signals (Table S1)²¹ were overexpressed in Arabidopsis thaliana, Nicotiana tabacum or Citrus paradisi. These virulence factors refer to genes that contribute to bacterial growth in plants but are not directly responsible for disease symptoms, and genes that contribute to virulence in non-plant hosts, such as serralysin. None of the overexpressed CLas proteins caused HLB-like symptoms, consistent with the bioinformatic analyses that CLas does not contain pathogenicity factors that directly cause HLB symptoms. Intriguingly, it was reported that CLas triggers immune responses^{22,23,24,25,26}, which we hypothesize is responsible for causing the devastating damages of HLB, similar to immune-mediated diseases of human.

CLas infection triggers immune response and cell death in the phloem tissue

Next, we tested if and how CLas triggers an immune response and cell death. Newly emerged citrus flush from HLB-positive citrus trees is initially free of CLas. We trimmed HLB-positive and healthy 2-year-old C. sinensis “Valencia” trees in a greenhouse to trigger young flush and conducted dynamic analyses of the relationship between CLas infection, immune response, cell death, and HLB symptom development. CLas was detected in young leaves of HLB-positive trees at ~15 days post-bud initiation based on quantitative PCR (qPCR) (Fig. 1). H₂O₂ (an indicator of reactive oxygen species (ROS)) content in CLas-positive flush was significantly higher than that of healthy plants at 15- and 21-day post-bud initiation (Fig. 1a). Significantly more callose deposition, an indicator of immune response²⁷, was observed in CLas-positive flush than that of the healthy plants starting at 18 days post-bud initiation and thereafter (Fig. 1b). The observation that CLas triggers ROS production earlier than callose deposition is consistent with the fact that ROS positively regulate callose deposition²⁸. On the other hand, significantly more starch accumulation, which was suggested to cause disruption of chloroplast inner grana structure and contribute to the yellowing symptoms²⁹, was observed in CLas-positive samples than in healthy samples beginning 18 days after bud initiation (Fig. 1c). Symptoms began to appear at ~40 days after bud initiation. It appears that CLas infection triggered plant immune responses, such as ROS (e.g., H₂O₂) production and callose deposition, that in turn, incited symptom development.

**Fig. 1: CLas infection causes ROS production, callose deposition, and starch accumulation in young citrus leaves.**

To better understand the nature of immune responses induced by CLas, we conducted temporal expression analyses of immune marker genes (PR1, PR2, PR3, and PR5) in young leaves at 15-, 18-, 21-, 24-, 27-, and 60-days post-bud initiation for the CLas infected and healthy C. sinensis plants. PR genes were consistently induced by CLas despite some fluctuations between 15- and 60-day-post-bud initiation (Fig. S1).

We observed cell death of sieve element and companion cells via transmission electron microscopy (TEM) of asymptomatic young leaves of HLB-positive C. sinensis “Valencia” trees (Fig. 2a–d), indicating cell death of phloem tissues occurs prior to the appearance of HLB symptoms. CLas was observed in phloem tissue with intact sieve element and companion cells, but not in sieve element cells undergoing the cell death process (Fig. 2c–e). It was noteworthy that some sieve element and companion cells underwent cell death while others remained intact in the same field (Fig. 2d), explaining the reduced function of the phloem, rather than its complete loss of function. This is consistent with the observation that CLas moves primarily vertically, but not laterally during infection³⁰. TEM observation of CLas-infected leaves confirmed that cell death was limited to sieve element and companion cells³¹, but not in surrounding parenchyma cells (Fig. 2c).

**Fig. 2: Transmission Electron Microscopy (TEM) analysis of asymptomatic young leaves of CLas-infected *Citrus sinensis*.**

We confirmed cell death in mature C. sinensis leaves exhibiting various symptoms using trypan-blue staining. No cell death was observed in healthy leaves collected from CLas-free plants. Cell death was observed in CLas-positive leaves and was more severe in leaves with severe symptoms (Fig. 3a), congruent with the TEM observation of the death of sieve element and companion cells in the midribs of CLas-infected mature leaves and CLas-positive stem tissues (Figs. S2, S3). Higher incidence of cell death was observed in leaves with higher CLas titers, suggesting CLas infection is responsible for cell death in the phloem tissue (Figs. S2, S3). TEM demonstrated that death of both companion cells and sieve element cells (Fig. 2c, d) or death of companion cells alone (Figs. 2b, S3c), but not death of sieve element cells alone occurred, indicating that death of companion cells occurred before that of sieve elements. We observed significantly higher H₂O₂ concentrations in CLas-infected mature leaves than CLas-free leaves (Fig. 3b), consistent with a previous study by Pitino et al.²⁶. In addition, significantly higher H₂O₂ concentrations were detected in the exudates of phloem-enriched tissues of CLas-positive stems than in stems of CLas-free trees (Fig. 3c).

**Fig. 3: *Candidatus* Liberibacter asiaticus (CLas) induces ROS production and cell death in phloem tissue of *C. sinensis*.**

Cell death is usually accompanied by ion leakage. Surprisingly, no difference in ion leakage between leaf blades or midribs of healthy, asymptomatic, mildly symptomatic, and severely symptomatic leaves was observed (Fig. S4), contrary to the results of TEM and trypan-blue staining (Figs. 2, 3a, S2). However, ion concentrations in exudates extracted from phloem-enriched bark tissues of CLas-infected samples were significantly higher than that of healthy samples (Fig. 3d), consistent with a model in which cell death occurs in phloem tissues, but not in surrounding parenchyma (Fig. 2c). Any change in ion leakage in the leaf blades and midribs of CLas-infected samples (Fig. S4) is probably masked by parenchyma cells because companion and sieve element cells make up only ~1% of the total cell population³².

Next, we used callose deposition as an indicator to investigate the localization of the immune response of citrus leaves in response to CLas infection. For this test, we investigated callose deposition in different sections of asymptomatic and symptomatic leaves of HLB-positive C. sinensis trees (Figs. 4a–g, S5). Callose deposition in the petiole, midrib, and lamina of asymptomatic leaves was significantly lower than that in symptomatic leaves (Figs. 4a–g, S5). No callose deposition was observed in the CLas-free lamina of asymptomatic leaves (Figs. 4c, S5). Together, the correlation between callose deposition and CLas titers suggests that CLas is responsible for inducing callose deposition (Figs. 4a–f, S5) as observed in the past²⁴. Callose deposition was observed only in phloem tissues as observed previously³³, but not in mesophyll cells. In contrast, pathogens infecting the apoplast, such as Xanthomonas, induce callose that is deposited between the plasma membrane and the cell wall in the vicinity of the pathogen^34,35. Hence, CLas induces a systemic immune response in phloem tissues following its systemic infection of the tree.

**Fig. 4: CLas infection induces callose deposition, H₂O₂ production and cell death in the phloem tissue of *C. sinensis*.**

To further verify that CLas induces immune responses in phloem tissues, we monitored ROS formation and localization in phloem-enriched bark tissues using the fluorescent probe 2′,7′-dichlorodihydrofluorescein diacetate (H₂DCFDA) and confocal laser microscopy. H₂DCFDA is a commonly used cell-permeable probe for measuring cellular H₂O₂³⁶. Significantly higher levels of H₂O₂ were detected in the phloem-enriched bark tissues of CLas-infected citrus plants than that of CLas-free plants (Fig. 4h–j).

Next, we aimed to further establish the causal relationship between CLas infection and ROS induction and cell death in phloem tissues. CLas-positive 5-year-old C. sinensis trees were treated with streptomycin via trunk injection to kill the pathogen³⁷. At 7 days after treatment, streptomycin significantly reduced CLas titers, H₂O₂ content, and ion leakage in the phloem tissue (Fig. 4k–l). Collectively, this establishes the causative relationship between CLas infection and ROS induction and cell death in phloem tissues.

HLB-mediated cell death is triggered by ROS

Cell death can be instigated by ROS³⁸. At high concentrations, ROS trigger necrotic cell death, but also induce programmed cell death at lower concentrations³⁹. We thus analyzed H₂O₂ contents in C. sinensis leaves infected by CLas. The H₂O₂ concentration in young CLas-infected leaves was ~6 µmol g^–1 FW, but reached 10–15 µmol g^–1 FW in mature leaves (Fig. 3b). Moreover, this method probably underestimated the H₂O₂ concentration in the phloem tissue because it could not differentiate levels in phloem cells, where H₂O₂ concentrates (Fig. 4h–j), from that in the more abundant parenchyma cells. Intriguingly, Xanthomonas citri subsp. citri, another bacterial pathogen of citrus, triggers H₂O₂ production by 2 days after inoculation of kumquat (Citrus japonica, syn: Fortunella crassifolia) which peaks (9.86 µmol g^–1 FW) 8 days after inoculation, eventually leading to cell death^40,41. X. citri subsp. citri-induced cell death is a slow process, occurring 6–8 days after inoculation⁴¹. Hence, we conjectured that ROS induced by CLas can reach a threshold necessary to trigger death of companion and sieve element cells in mature leaves. ROS production triggered by CLas is distinct from that triggered by incompatible pathogens, which is typified by a biphasic oxidative burst⁴². Instead, high levels of ROS were consistently detected in both young leaves during early stages of infection (Fig. 1a) as well as in mature infected leaves and stems (Fig. 3b, c), probably triggered by CLas colonizing and multiplying in the previously unoccupied phloem tissue as well as that triggered by damage-associated molecular patterns resulting from dying companion and sieve element cells.

H₂O₂ concentrations in the phloem-enriched bark tissues from symptomatic (1.80 ± 0.13 mmol/L) branches were much higher than that (0.59 ± 0.01 mmol/L) of healthy trees (Fig. 3c). H₂O₂ induces necrosis of immortalized rat embryo fibroblasts at a concentration of 700 μmol/L⁴³. H₂O₂ at a concentration of 1.8 mmol/L but not 0.6 mmol/L or lower induced cell death of C. sinensis protoplast cells (Figs. 5a, b, S6). Similar results were observed for suspension culture cells (Fig. S7). Addition of uric acid (0.2 mM), a ROS scavenger, reduced both H₂O₂ concentrations (Fig. 5c) and cell death (Fig. 5a, b), indicating that H₂O₂ induced by CLas is sufficient to kill phloem tissue. It is important to note that the growth of plants (e.g., Arabidopsis) is inhibited by 1 mM H₂O₂⁴⁴, partly explaining the growth stunting phenotype of CLas-infected young citrus trees.

**Fig. 5: ROS are responsible for cell death of phloem tissues of CLas-infected citrus.**

In addition to H₂O₂, ROS induced by pathogens include hydroxyl radicals, superoxide anions, and singlet oxygen. To further corroborate that HLB caused cell death is instigated by ROS, we sprayed HLB-positive C. sinensis ‘Valencia’ trees weekly with the antioxidants uric acid (1.8 mM) and rutin (0.6 mM). After treatment for 6 weeks, ROS levels in phloem-enriched bark tissues were much lower than control plants (Fig. 5c), and the uric acid and rutin treatments reduced cell death (Fig. 5d). Taken together, CLas infection of citrus phloem tissue induces ROS production, which subsequently causes cell death of phloem tissue.

CLas infection significantly altered expression of pathways related to oxidative stress and immune responses

Next, we investigated 20 previously published gene expression profiles of C. sinensis response to CLas infection that studied different tissues (leaf, stem, and fruit), different environments (greenhouse or groves), and different infection stages (Table S2). Analyses of differentially expressed genes (DEGs) clearly demonstrated that the expression of genes related to ROS and immune response are significantly altered by CLas infection (Table S3). The combined analysis showed an overall downregulation of antioxidant enzymes and upregulation of transmembrane localized NADPH oxidases, known as RBOHs, explaining the oxidative stress response in response to CLas infection (Fig. S8). Critically, expression of respiratory burst oxidative homolog D (RBOHD) gene, which encodes NADPH oxidase implicated in the generation of ROS during defense responses, was induced by CLas infection in several studies. RBOHD is primarily responsible for ROS produced upon PAMP recognition and is required for cell death that is initiated after pathogen detection⁴⁵. In addition, RBOHB and RBOHF have also been reported to be involved in ROS production in response to pathogen infection^45,46. qRT-PCR assays revealed that both RBOHB and RBOHD were induced by CLas infection of leaf samples collected in both field and greenhouse conditions, whereas RBOHF was induced only under field conditions (Fig. S9). Intriguingly, ROS levels in CLas-positive stems were reduced by NADPH oxidase inhibitor diphenyleneiodonium (DPI) (Fig. S10), supporting the notion that RBOHs contribute to the ROS accumulation triggered by CLas, which probably also results from contributions from other components such as peroxidases.

Taken together, our analysis implies that CLas infection causes oxidative stress to citrus in most conditions, consistent with the H₂O₂ production it triggers (Figs. 1a, 3b, c, 4h–k, 5c). Intriguingly, the involvement of oxidative stress in HLB disease has been suggested by multiple previous studies^26,47,48. Comparative transcriptome analysis revealed complex expression changes associated with the immune response pathways in response to CLas infection (Fig. S11). Approximately 66 NLR genes showed overall induction by CLas in many studies (Fig. S11B). Transcriptome analyses of sweet orange response to CLas infection therefore support our hypothesis that HLB is an immune-mediated disease that results from ROS induced cell death of phloem tissue triggered by the pathogen.

Suppressing ROS-mediated cell death mitigates HLB symptoms

Antioxidants, and immunoregulators are commonly used to treat human immune-mediated diseases by halting or reducing ROS-mediated cell death^49,50,51. Correspondingly, we tested whether growth hormone gibberellin (GA), and antioxidants (uric acid and rutin) could mitigate the ROS-mediated cell death triggered by CLas infection, and thus block or reduce HLB symptoms. The plant growth hormone GA was selected because it is known to modulate PAMP-triggered immunity and PAMP-induced plant growth inhibition⁵². Both uric acid and rutin were also assessed since they are well-known ROS scavengers^53,54.

Foliar sprays of HLB-positive C. sinensis trees with GA at both 5 mg/L and 25 mg/L reduced tissue H₂O₂ levels and ion leakage (Fig. 6a, b). GA (5 mg/L) also suppressed death of C. sinensis protoplast cells treated with 1.8 mM H₂O₂ (Fig. 6c, d). RNA-seq analyses of GA treated vs. non-GA treated C. sinensis protoplast cells in the presence of 1.8 mM H₂O₂ demonstrated that GA induced the expression of genes encoding H₂O₂ scavenging enzymes catalases, ascorbate peroxidases, and glutathione peroxidases. GA also inhibited the expression of RBOHD (Fig. 7, Supplementary Data 1). GA treatment therefore clearly alleviates the oxidative stress caused by H₂O₂ in citrus.

**Fig. 6: The immunoregulator gibberellin (GA) suppresses HLB development.**

**Fig. 7: The expression profiling of ROS related genes between GA and non-GA (nGA) treated *Citrus sinensis* protoplast cells in the presence of 1.8 mM H₂O_2.**

Six weeks after initiation of foliar sprays, treated plants exhibited reduced HLB symptoms (i.e., less blotchy mottle) compared with that of plants before treatment, whereas the blotchy mottling symptoms of plants treated only with water became more severe in this period (Fig. S12). In addition, GA promoted plant growth as indicated by new leaf flushes that were not observed on water-treated plants (Fig. S12). GA is registered for use on citrus in the US, enabling its use in large-scale field trials of its efficacy for control of HLB. Foliar spray of C. sinensis with GA significantly reduced HLB disease symptoms 8 months after application in such large-scale trials (Figs. 6e–g, S13). The treated trees appeared much healthier than the non-treated control trees even though all were infected (Fig. 6g). Because symptomatic leaves demonstrated significantly more cell death than asymptomatic leaves based on trypan-blue staining (Fig. 3a) and TEM observation (Fig. S2), we used the proportion of symptomatic leaves on a tree as an indicator of cell death caused by HLB. GA treatment significantly reduced the incidence of symptomatic leaves (i.e., leaves showing blotchy mottle, yellowing, and nutrient deficiency) (Figs. 6e–g, S13), suggesting reduced death of sieve element and companion cells in treated leaves. In addition, foliar sprays of GA on HLB-positive 6-year-old C. sinensis var. ‘Valencia’ and var. ‘Vernia’ significantly promoted plant growth as measured by tree height, trunk diameter and canopy volume when assessed 8 months after application (Fig. S14). We presume GA reduces HLB symptoms via its direct effect on both mitigating ROS (Figs. 6a, c, d, 7) and promoting plant growth (Fig. S14), thus alleviating the growth inhibition caused by CLas.

To determine whether mitigating ROS can directly halt or reduce HLB symptoms, we sprayed HLB-positive symptomatic C. sinensis trees weekly with the antioxidants uric acid (1.8 mM) and rutin (0.6 mM). Remarkably, 6 weeks after initial treatment, both antioxidants significantly reduced HLB symptoms compared to that observed before treatment, whereas the plants treated with water became more symptomatic in the same duration (Fig. S12). Taken together, our data suggest that suppression of ROS-mediated cell death can mitigate HLB symptoms. Consequently, we have established the causal relationship that CLas triggers ROS production in the phloem tissue, which subsequently causes cell death of phloem tissue, leading to HLB symptoms.

Discussion

In this study, we provide evidence that citrus HLB is an immune-mediated disease. CLas induces a systemic chronic immune response, mimicking systemic chronic inflammation diseases of humans⁵⁵. Systemic chronic inflammation diseases may result from collateral damage to tissues and organs over time by oxidative stress⁵⁶. ROS concentrations triggered by CLas infection are above the threshold needed to induce cell death, probably resulting from the combined effect of programmed cell death induced at low ROS concentrations and necrotic cell death stimulated at high ROS concentrations. In addition, ROS positively regulate callose deposition^28,57 and inhibit plant growth including roots^58,59. Persistent induction of ROS by systemic CLas infection leads to systemic cell death of phloem tissue and other effects owing to diverse roles of ROS, which subsequently affect phloem function, hormone synthesis and transport, metabolic transport, and rerouting energy to immune responses rather than to growth. This model can explain most HLB phenomena. For instance, phloem dysfunction resulting from the death of companion and sieve element cells can lead to starch accumulation, and the resultant blotchy mottle symptoms. Hardened leaves likely result from the action of ROS which are known to directly cause strengthening of host cell walls⁴². Cell death of the phloem tissue, reduced transport of photosynthates, and ROS inhibition of root growth may be responsible for root decay. Stunted growth probably results from the direct effect of ROS and reduced transport of carbohydrates and hormones. The detailed molecular mechanism of how CLas activates immune responses remains unknown. We anticipate that cytoplasmic receptors, such as nucleotide-binding leucine-rich repeat (NLR) proteins are primarily responsible for intracellular detection of CLas through recognition of PAMPs inside companion and sieve element cells. Although they have not been previously known for plants, immune-mediated diseases are likely prevalent in the Kingdom Plantae, particularly for diseases caused by phloem-colonizing bacteria (e.g., Ca. Liberibacter, Spiroplasma, and Ca. Phytoplasma), viruses, and fungi as well as non-phloem-colonizing pathogens.

The finding of citrus HLB as an immune-mediated plant disease helps guide the battle against this notorious disease. It seems likely that horticultural approaches that suppress oxidative stress can alleviate the immune-mediated damage caused by CLas in HLB endemic citrus production areas. These approaches include optimized usage of plant growth hormones, such as GA and brassinosteroids⁶⁰, and antioxidant treatments. Even though we did not test the effect of nutritional modulation of immune function on HLB in this study, citrus growers in Florida have observed that modulation of macronutrients (N, P, and K) and micronutrients (e.g., B, Cu, Fe, Mn, Mo, Ni, Se and Zn) can reduce HLB symptoms. This is consistent with the observation that a deficiency of macronutrients can lead to oxidative stress⁶¹, whereas micronutrients (B, Cu, Fe, Mn, Mo, Ni, Se and Zn) at low concentrations can activate antioxidative enzymes⁶². Antioxidants (e.g., uric acid), growth hormones (e.g., GA) and nutritional modulation (e.g., micronutrients) could be used to directly alleviate oxidative stress to reduce cell death of the phloem tissue to mitigate HLB symptoms. Moreover, growth hormones and nutrients, by promoting new growth, decrease the proportion of dead cells in phloem tissue, further mitigating HLB symptoms. The horticultural measures used to mitigate ROS and cell death are unable to reduce or eliminate CLas inoculum, thus are not recommended for citrus production areas with low HLB incidence. For those areas, region-wide comprehensive implementation of roguing infected trees, tree replacement, and insecticide applications has been shown to successfully control citrus HLB^17,63. Genetic improvements that enhance plant tolerance of oxidative stress, prevent overproduction of ROS, or evade recognition of CLas may generate HLB resistant/tolerant citrus varieties. In particular, enhanced companion cell- or phloem-specific overexpression of antioxidant enzymes (such as superoxide dismutase, catalases, glutathione peroxidases, ascorbate peroxidase, and glutathione reductase) could be achieved using CRISPR-mediated gene editing, transgenic, or cisgenic approaches or citrus tristeza virus vectors). Additionally, gene editing of the NLRs responsible for recognizing CLas might also render citrus tolerant to HLB. It is plausible that some of the HLB-tolerant citrus genotypes might not recognize PAMPs or proteins of CLas or the recognition is too weak to reach the ROS threshold to cause cell death, thus avoiding the immune-mediated damages. On the other hand, some HLB-tolerant citrus genotypes might be superior in antioxidant mechanism or phloem regeneration. Intriguingly, it has been reported that HLB-tolerant varieties contain higher levels of antioxidants than more susceptible varieties⁶⁴. “Sugar Belle” mandarin, which is HLB-tolerant, has more vigorous phloem regeneration than the HLB-susceptible C. sinensis varieties¹⁵. In summary, citrus HLB is an immune-mediated disease and mitigating ROS via antioxidant mechanisms and promoting new growth both can reduce cell death of phloem tissues, thus controlling HLB.

Methods

Transgenic expression analysis of CLas proteins containing Sec secretion signals and other putative virulence factors

Transgenic expression of CLas genes was conducted⁶⁵. For the citrus transformation, CLas genes were amplified without signal peptide sequence and cloned into the binary vector RCsVMV-erGFP-pCAMBIA-1380N-35S-BXKES-3xHA, which has a C-terminal 3×HA tag, to generate the CLas gene overexpression vectors. The resulting binary vectors were transferred into Agrobacterium tumefaciens strain EHA105 for citrus transformation. The empty vector (EV) was used in citrus transformation as a negative control. Agrobacterium-mediated transformation of epicotyl segments of Duncan grapefruit (Citrus paradisi) was carried out to generate transgenic citrus plants⁶⁶. Transgenic lines showing kanamycin-resistance and erGFP-specific fluorescence were selected and then micro-grafted in vitro onto 1-month-old Carrizo citrange rootstock seedlings. After 1 month of growth in vitro, the grafted shoots were potted into a peat-based commercial potting medium and acclimated under greenhouse conditions for the phenotype evaluation. Transgenic plants were confirmed by PCR, qRT-PCR at the RNA level, or western blot using HA Tag Antibodies (Cat. #11867423001, Roche).

For the tobacco transformation, Agrobacterium-mediated transformation of leaf discs of Nicotiana tabacum was carried out to generate the transgenic tobacco⁶⁷. A. tumefaciens strain EHA105 containing the vectors was used for transformation. Transgenic positive shoots showing kanamycin-resistance and erGFP-specific fluorescence were selected and transferred to the rooting medium. Evaluation of the transgenic N. tabacum was conducted in a growth chamber. Transgenic plants were confirmed by PCR, qRT-PCR at the RNA level, or western blot using HA Tag Antibodies (Cat. #11867423001, Roche).

For the gene overexpression in Arabidopsis thaliana. CLas genes without signal peptides were PCR amplified and cloned into the binary vector pCambia1380-35S-EYFP, which has a C-terminal EYFP protein tag, and transferred into A. tumefaciens strain GV2260. Agrobacterium-mediated floral dip method was used for the Arabidopsis transformation⁶⁸. The T1 generation transgenic plants were screened on the Hygromycin B selection medium. Positive plants were further confirmed by PCR and western blot using GFP antibodies (Cat. #SAB4301138, Sigma-Aldrich). The positive T2 generation transgenic plants were evaluated in a growth chamber. Both an-HA antibody and GFP antibody were confirmed to be functional (Fig. S15).

Plant materials used for investigation of the relationship between CLas infection, immune response, phloem blockage, cell death and HLB symptom development

We used 2-year-old CLas-infected and healthy Valencia sweet orange (Citrus sinensis) plants maintained in a greenhouse (28 ± 2 °C, relative humidity of 50 ± 5%, natural light period). We selected young flushes from the twigs of HLB-positive sweet orange trees. Both “Valencia” and “Hamlin” sweet orange are C. sinensis and susceptible to HLB without observable differences in symptoms. The plant was infected with CLas by graft inoculation and maintained in a greenhouse. The sweet orange trees from the groves were naturally infected by CLas. Healthy plants were maintained in a glasshouse with natural light and without temperature control.

Quantification of H₂O₂ concentrations

H₂O₂ concentrations were quantified following the procedure described elsewhere⁴⁰. CLas-positive asymptomatic mature leaves, mature leaves with mild or severe symptoms were collected from HLB-positive C. sinensis “Valencia” trees in citrus groves of Citrus Research and Education Center, University of Florida/Institute of Food and Agricultural Sciences. CLas negative mature leaves were collected from healthy C. sinensis “Valencia” trees in glasshouse. Briefly, leaf samples (0.5 g) were grinded in 0.1% (w/v) trichloroacetic acid (TCA) and centrifuged at 12,000 × g for 15 min at 4 °C. The supernatant (0.3 mL) was mixed with 1.7 mL 1.0 M potassium phosphate buffer (pH 7.0) and 1.0 mL of 1.0 M potassium iodide solution, then incubated for 5 min before measuring the absorbance of the oxidation product at 390 nm. H₂O₂ concentrations were calculated using a standard curve prepared with known concentrations of H₂O₂ and expressed in µmol/g fresh weight. For measuring H₂O₂ concentrations in the exudates of phloem-enriched bark tissues, the same procedure was used except the TCA step and H₂O₂ concentrations were expressed in mmol/L. CLas-positive symptomatic branches and CLas-negative branches from the spring flush were used for collection of bark tissues from the C. sinensis “Valencia” trees mentioned above.

Ion leakage

Conductivity of the exudates extracted from phloem-enriched bark tissues was measured using a CON 700 conductivity/°C/°F bench meter (Cat. #WD-35411-00, OAKTON Instruments, Vernon Hills, IL, USA).

Callose deposition assay

Leaf samples were fixed with FAA (37% formaldehyde/glacial acetic acid/95% ethanol/deionized water at a volume ratio of 50:5:10:35) solution overnight. Samples were embedded in the Tissue Plus O.C.T compound (Cat. #23-730-571, Thermo-Fisher, Waltham, MA, USA), sectioned with a Harris Cryostat Microtome (International Equipment, Boston, MA, USA) and stained with 0.005% aniline blue solution prior to analysis. Samples were observed in an Olympus BX61 epifluorescence microscope (Olympus Corporation, Center Valley, PA, USA). Callose spots were counted per slide area for all sample types.

CLas quantification using qPCR

Tissues (100 mg) were homogenized into powders using a TissueLyser II (Cat. #85300, Qiagen, Valencia, CA, USA). DNA was extracted using the DNeasy Plant kit (Cat. #69204, Qiagen), following the manufacturer’s instructions, and eluted in 100 µL nuclease free water. DNA concentration was measured using a Synergy LX plate reader (BioTek, Winooski, VT, USA). Quantification of CLas in plant tissues was performed as described elsewhere⁶⁹. Briefly, qPCR was carried out with primers and probe for CLas⁷⁰. qPCR assays were performed with QuantiStudio3 (Cat. #A28567, Thermo-Fisher, Waltham, MA) using the Quantitec Probe PCR Master Mix (Cat. #204343, Qiagen) in a 25-µl reaction. The standard amplification protocol was 95 °C for 10 min followed by 40 cycles at 95 °C for 15 s and 60 °C for 60 s. All reactions were conducted in triplicate with CLas-positive and water controls. Quantification of CLas was conducted using the following equation:⁷¹

$${{{{{\rm{Y}}}}}}\,=\,-0.288\,\times\, ({{{{{\rm{CLas}}}}}}\,{{{{{\rm{Ct}}}}}})\,+\,11.607$$

Starch assay

The samples (100 mg) were powdered using a TissueLyser II (Qiagen, Hilden, Germany). The powdered samples were used to quantify the starch. The starch estimation was performed using the Total Starch Assay Kit (AA/AMG) (Cat. #K-TSTA-100A, Megazyme, Bray, Ireland) as instructed by the manufacturer. The experiments were repeated thrice with similar result.

Statistical analyses

All statistical analyses were performed using SAS statistical software (Version 9.4, SAS Institute, Cary, NC, USA).

Gene expression assays using reverse transcription quantitative PCR (RT-qPCR)

Total RNA was extracted using the RNeasy Plant Mini Kit (Cat. #74904, Qiagen), according to manufacturer’s instructions. cDNA was synthesized with Quantitect Reverse Transcription Kit (Cat. #205311, Qiagen) according to manufacturer’s instructions and diluted 10 times for RT-qPCR. Reactions were carried out by adding 1 µL of cDNA, 1 µL of each specific primer, 7 µL water and 10 µL Fast SYBR Green Master Mix (Cat. #4385617, Thermo-Fisher Scientific, Waltham, MA, USA) performed with QuantiStudio3 (Thermo-Fisher) using the standard fast protocol of 95 °C for 20 s followed by 40 cycles of 95 °C for 1 s and 60 °C for 20 s. Denaturation protocol consisted of 95 °C for 1 s, 60 °C for 20 s and a final dissociation step of 95 °C. Relative gene expression was calculated using the $2^{-\Delta\Delta {C_{T}}}$method⁷². CsGAPDH was used as an endogenous control.

TEM analysis

Small sections of the leaf and stem samples were collected under a stereomicroscope (Swift Table Stereo Zoom Microscope, Carlsbad, CA, USA). The leaf samples were transferred to 3% glutaraldehyde overnight at 4 °C for fixation. Then, the samples were postfixed in 2% osmium tetroxide prepared in 3% glutaraldehyde for 4 h at room temperature in a fume hood. The samples were dehydrated by sequential treatment with 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100% (thrice) acetone for 10 min each. The leaf samples were incubated sequentially in 50, 75, and 100% (twice) Spurr’s low-viscosity epoxy resin prepared in acetone for 8 h each. One-micrometer sections were cut with glass knives using an ultramicrotome followed by staining with methylene blue/azure A for 30 s and basic fuchsin (0.1 g in 10 mL of 50% ethanol) for 30 s. The sections were observed under a Leitz Laborlux S compound microscope (Leica Microsystems, Wetzlar, Germany) for the right spot with a vascular system. The same blocks were trimmed with a surgical blade and then sectioned to 0.1 µm using a diamond knife under an ultramicrotome. The thin sections were collected on 200-mesh copper grids. The samples were stained with 2% aqueous uranyl acetate for 5 min, washed in water, and again stained with lead citrate followed by water wash. The micrographs were prepared and analyzed using a Morgagni 268 (FEI Company, Hillsboro, OR, USA) transmission electron microscope equipped with an AMT digital camera (Advanced Microscopy Techniques Corp., Danvers, MA, USA).

Trypan-blue staining

Trypan-blue staining was conducted as described by Fernández-Bautista et al.⁷³.

Monitoring ROS formation and localization in phloem tissues by use of the fluorescent probe 2′,7′-dichlorodihydrofluorescein diacetate (H₂DCFDA) and confocal laser microscopy

Four CLas-infected branches collected from HLB-positive C. sinensis “Valencia” trees and four CLas-free branches collected from healthy C. sinensis “Valencia” trees were collected and placed in glass test tube with 30 mL water containing 10 μM 2′,7′-dichlorodihydrofluorescein diacetate (H₂DCFDA) (Cat. #D399, Thermo-Fisher Scientific). Four leaves/branch at the top were kept to facilitate transpiration. Glass tubes were wrapped with aluminum foil and kept in room temperature for 24 h. Bark was peeled from the stem section that was submerged in water and placed on slide with inner side upwards. HLB-positive and healthy branches were also incubated in water without H₂DCFDA as controls. 2′,7′-dichlorofluorescein (DCF) fluorescence was visualized by confocal laser scanning microscopy (CLSM) (Leica TCS‐SP5, Mannheim, Germany) with excitation/emission at 495 nm/525 nm.

Trunk injection of HLB-positive 5-year-old C. sinensis trees

Trunk injection was conducted as described elsewhere³⁷. For each tree, approximately 0.4 g streptomycin sulfate (Cat.#AC612240500, laboratory grade; Thermo-Fisher Scientific) at 5 g/L was injected into the trunk of three trees. The amount of streptomycin injected was calculated to reach the concentration needed to kill CLas in planta based on the canopy volume³⁷.

Exudates of phloem-enriched bark tissues for H₂O₂ and ion leakage assays

Exudates of phloem-enriched bark tissues were extracted from stems following the procedure described elsewhere⁷⁴. Stems were collected from small branches of spring spouts with mildly symptomatic leaves.

Protoplast

Protoplast cells of C. sinensis “Hamlin” were prepared as described by⁷⁵. Embryogenic calli were subcultured on solid MT (Murashige and Tucker) media (Cat. #M5525, Phytotech) every 2 weeks. From the maintained calli suspension, cells were prepared and maintained in DOG liquid media as described elsewhere⁷⁶. The final isolated protoplast cells were suspended in W5 solution (154 mM NaCl,125 mM CaCl2, 5 mM KCl, 2 mM MES at pH 5.7) at 1 × 10⁷ cells/mL for different treatments.

Treatment of protoplast cells with H₂O₂, antioxidants, and Gibberellic acid (GA)

For H₂O₂ treatment, H₂O₂ was freshly prepared as 1 M stock solution with sterile double stilled water. For protoplast treatments with different concentrations of H₂O₂, H₂O₂ was further prepared as 100× stock with protoplast buffer (W5 solution). Protoplast cells were then treated with different concentrations of H₂O₂ for 24 h. For each treatment, at least three biological replicates were conducted. After 24 h treatment, samples were stained with Fluorescein Diacetate (FDA) (Cat. #F1303, Invitrogen) for viability observation. For 50 µL of sample, 2 µL of FDA was added. Immediately after staining, the samples were observed under Olympus BX53 Epi-fluorescence microscope with green channel. The ratio of green cells (living cells) to total cells was calculated as viability rate. All steps and chemical treatments were performed at room temperature.

Protoplast cells co-treated with H₂O₂, antioxidants or Gibberellic acid (Alfa Aesar) were conducted as described above. Uric acid (Cat. # A13346, Thermo-Fisher Scientific) was dissolved in protoplast buffer (W5 solution). Rutin hydrate (Cat. # R5143, Sigma-Aldrich) was dissolved in DMSO as stock, Gibberellic acid (Cat. #A17843, Alfa Aesar) was dissolved in ethanol as stock.

Treatment of citrus suspension culture cells with H₂O₂

Sweet orang Hamlin suspension culture 7 days after subculture was used. Five mL of suspension culture cells were aliquoted into a 50-mL Falcon tube. Freshly prepared H₂O₂ was added into each tube at a concentration of 0 (water control), 0.6, 1.5, 1.8, or 3.6 mM. The tubes were incubated at room temperature with gentle shaking (100 rpm). Twenty-four hours after treatment, 50 µL of cells were pipetted into a 1.5-mL tube from each treatment. Each sample was stained with 2 µL of fluorescein diacetate (stain only living cells, green color) and 2 µL of propidium iodide (stain only dead cells, orange to red color). One minute after staining, the stained samples were observed under a fluorescent microscope with green and red dual channels.

Foliar spray with antioxidants and GA

Five-year-old Valencia sweet orange trees with similar symptoms were used for foliar spray treatments. All trees in the grove were HLB-positive. The experiment was a completely randomized design with 5 treatments. Each treatment consisted of four trees. The treatments were applied by foliar spray with 2.5 mL/L of Induce non-ionic surfactant (Helena Ag, Collier, TN, USA). One liter of solution per plant were applied at approximately 400 kPa using a handheld pump sprayer to apply on the whole tree. This pressure resulted a fine mist and was sufficient to produce runoff from the leaves to ensure complete coverage. The five treatments were applied to the various trees as follows: uric acid (1.8 mM), rutin (0.6 mM), GA (5 mg/L), GA (25 mg/L) with water as the negative control. Foliar spray was conducted in the evening to facilitate absorption. The chemicals GA, and uric acid were purchased from Fisher Scientific. Rutin was purchased from Sigma-Aldrich (St. Louis, MO, USA).

GA treatment via foliar spray

GA foliar spay was conducted in the first week of November, 2020. For the GA application, 20 ounces of Pro Gibb LV (Valent U.S.A. LLC, Walnut Creek, CA, USA) was mixed with water in a 100 Gal tank. 64 ounces of WIDESPREAD MAX (A.I. organosilicone) was included as the surfactant for leaf spray with airblast. Applications were conducted during night. One block of Valencia sweet orange on rootstock 942 was treated with GA, whereas the nearby Valencia/942 block was not treated with GA as a negative control. In addition, one block of Vernia sweet orange on X639 rootstock was treated with GA with one nearby Vernia/X639 block as a negative control. All blocks are 10 acres or more with approximately 140 trees/acre.

RNA-seq analyses of GA treatment on citrus protoplast cells in the presence of H₂O₂

Protoplast cells were prepared as described above and suspended in W5 solution (154 mM NaCl,125 mM CaCl₂, 5 mM KCl, 2 mM MES at pH 5.7) at 1 × 10⁷ cells/mL. All steps were performed at room temperature. The following two treatments were conducted: (1) Protoplast + H₂O₂ (1.8 μmol/mL) + Gibberellin (5 mg/L), and (2) Protoplast + H₂O₂ (1.8 μmol/mL). Protoplast cells were maintained at room temperature without shaking. RNA was collected at 6 h after treatment. Four biological replicates were included for each treatment.

Total RNA was extracted using the RNeasy plant kit (Qiagen, Valencia, CA), followed by treatment with RQ1 RNase-Free DNase (Cat. #M6101, Promega, Madison, WI). RNA concentration and quality were measured by a Nanodrop One Microvolume UV-Vis Spectrophotometer (Thermo-Fisher Scientific, Waltham, MA). Samples meeting the following requirements (Concentration ≥20 ng/μL, OD260/280 > 2.0) were sent to Novogene (Novogene, Davis CA) for cDNA libraries construction and RNA-seq analyses. Libraries were constructed with the NEBNext Ultra II RNA Library Prep Kit for Illumina (Illumina, San Diego, CA). Samples were sequenced to generate 150 bp paired-end reads using the Illumina NovaSeq 6000 platform (Illumina). Raw RNA-seq data were filtered by removing low-quality reads and adapters, and then aligned to the Citrus sinensis reference genome⁷⁷ using HISAT2 version 2.2.1⁷⁸ and SAMtools version 1.2⁷⁹ each gene was quantified using HTSeq-count version 0.11.2⁸⁰. Different expressed genes (DEGs) analysis was performed using DESeq2 packages version 1.30.1 in R version 4.0.5⁸¹. Genes were considered significantly expressed with adjusted p value <0.05 (FDR method). The heatmap plots of expression profiling of DEGs were drawn using the pheatmap package version 1.0.12 in R program version 4.0.5⁸².

NADPH oxidase inhibitor diphenyleneiodonium (DPI) treatment on ROS levels in CLas-positive stems

HLB-positive branches from the summer flush of Valencia sweet orange in the field were collected and then soaked in DPI (Cat. #AAJ64838MC, Thermo-Fisher) solution (25 μM) or water (control). After 48 h treatment, phloem-enriched bark tissues were collected for H₂O₂ concentration measurement. Experiments were repeated two times and representative result is shown.

Evaluation of citrus tree growth and HLB symptoms in response to GA treatment

Tree growth was evaluated by estimating trunk diameter, tree height (TH), and canopy volume (CV) on both GA treated and untreated control trees. For each treatment group, a total of 10 trees (n = 10) were randomly selected for the evaluation. A digital caliper (Fowler, Newton, MA) was used to take two measurements of trunk diameter (north-south and east-west orientation) at ~20 cm above the ground. A tape measure was used to measure the TH above the ground from the soil surface to the apical point of the plant. CV was estimated by taking the average of two independent measurements of the diameter of the canopy at different directions (north-south and east-west). The CV was estimated using the equation: V = (2/3) × p × h × (d/2)², where h is the TH and d is the average diameter of the tree canopy⁸³. All statistical analyses were performed using SAS V9.4 (SAS Institute Inc., Cary, NC). The data were first tested for normality and homogeneity of variance using the Shapiro–Wilk’s test and Levene’s test, respectively. A Student’s two-tailed t test was performed to explore differences between GA treated and untreated control trees in growth performance traits.

HLB disease incidence in different treatments was evaluated by randomly checking 200 trees/treatment. Ratio of symptomatic leaves vs. total leaves in different treatments were investigated by evaluating three groups of branches/treatment with each group containing 16 branches that were selected randomly from 8 trees (2 branches/tree).

Data analyses of previously published RNA-seq data between HLB-positive and healthy samples

To generate the comprehensive expression pattern of sweet orange in response to CLas infection, we collected 15 microarray and 5 RNA-seq data sets from NCBI SRA and GEO databases (Table S2). The differentially expressed genes (DEGs) were determined using Limma version 3.46.0⁸⁴ and DESeq2 version 1.30.1⁸¹ packages in R version 4.0.5 for microarray and RNA-seq data, respectively (adjusted p value <0.05 and |log2 fold change|>1). Gene ontology (GO) term enrichment of DEGs was conducted using agriGO v2.0: a GO analysis toolkit for the agricultural community⁸⁵ using the singular enrichment analysis tool. The heatmap plots were drawn using the gplots package version 3.1.1 in R program version 4.0.5⁸⁶.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

The raw RNA-seq data generated in this study have been deposited in the NCBI Bioproject database under accession code PRJNA780217. Public transcriptomic data were obtained from NCBI GEO and SRA databases under accession numbers GSE33003, GSE33004, GSE33459, GSE29633, GSE33373, GSE101381, and SRP022979. Source data are provided with this paper.

References

Yu, X., Feng, B., He, P. & Shan, L. From chaos to harmony: responses and signaling upon microbial pattern recognition. Annu. Rev. Phytopathol. 55, 109–137 (2017).
Article CAS PubMed PubMed Central Google Scholar
Couto, D. et al. The arabidopsis protein phosphatase PP2C38 negatively regulates the central immune kinase BIK1. PLoS Pathog. 12, e1005811 (2016).
Article PubMed PubMed Central Google Scholar
Cui, H., Tsuda, K. & Parker, J. E. Effector-triggered immunity: from pathogen perception to robust defense. Annu. Rev. Plant Biol. 66, 487–511 (2015).
Article CAS PubMed Google Scholar
Warner, E. A. & Moldawer, L. L. Using innate immunity to characterize the host response to microbial invasion in severe sepsis. Future Microbiol. 3, 177–189 (2008).
Article CAS PubMed Google Scholar
Neal, L. M. et al. CD4(+) T cells orchestrate lethal immune pathology despite fungal clearance during cryptococcus neoformans meningoencephalitis. mBio 8, e01415-17 (2017).
Article PubMed PubMed Central Google Scholar
Casadevall, A. & Pirofski, L. A. The damage-response framework of microbial pathogenesis. Nat. Rev. Microbiol. 1, 17–24 (2003).
Article CAS PubMed PubMed Central Google Scholar
Adapa, D., Sai, Y., Anand, S., Mehaboobi, S. & Aramalla, E. A brief review on immune mediated diseases. J. Clin. Cell. Immunol. S11, 001 (2011).
Google Scholar
Barth, K., Remick, D. G. & Genco, C. A. Disruption of immune regulation by microbial pathogens and resulting chronic inflammation. J. Cell Physiol. 228, 1413–1422 (2013).
Article CAS PubMed PubMed Central Google Scholar
Speth, E. B., Lee, Y. N. & He, S. Y. Pathogen virulence factors as molecular probes of basic plant cellular functions. Curr. Opin. Plant Biol. 10, 580–586 (2007).
Article PubMed PubMed Central Google Scholar
Hu, Y. et al. Lateral organ boundaries 1 is a disease susceptibility gene for citrus bacterial canker disease. Proc. Natl Acad. Sci. U.S.A. 111, E521–E529 (2014).
Article CAS PubMed PubMed Central Google Scholar
Chatterjee, S., Almeida, R. P. & Lindow, S. Living in two worlds: the plant and insect lifestyles of Xylella fastidiosa. Annu. Rev. Phytopathol. 46, 243–271 (2008).
Article CAS PubMed Google Scholar
Wang, N. The citrus huanglongbing crisis and potential solutions. Mol. Plant 12, 607–609 (2019).
Article CAS PubMed Google Scholar
Singerman, A. & Rogers, M. E. The Economic Challenges of Dealing with Citrus Greening: The Case of Florida. Journal of Integrated Pest Management 11, https://doi.org/10.1093/jipm/pmz037 (2020).
Yuan, X. et al. Region-wide comprehensive implementation of roguing infected trees, tree replacement, and insecticide applications successfully controls citrus huanglongbing. Phytopathology 111, 1361–1368 (2021).
Article PubMed Google Scholar
Deng, H. et al. Phloem regeneration is a mechanism for huanglongbing-tolerance of “Bearss” lemon and “LB8-9” sugar belle. Front. Plant. Sci. 10, 277 (2019).
Article PubMed PubMed Central Google Scholar
Huang, C.-Y. et al. A stable antimicrobial peptide with dual functions of treating and preventing citrus Huanglongbing. Proc. Natl Acad. Sci. U.S.A. 118, e2019628118 (2021).
Article CAS PubMed PubMed Central Google Scholar
Bové, J. M. Huanglongbing: a destructive, newly-emerging, century-old disease of citrus. J. Plant Path. 88, 7–37 (2006).
Google Scholar
Wang, N. et al. The candidatus liberibacter-host interface: insights into pathogenesis mechanisms and disease control. Annu. Rev. Phytopathol. https://doi.org/10.1146/annurev-phyto-080516-035513 (2017).
Kuzmanović, N., Puławska, J., Hao, L. & Burr, T. J. The ecology of agrobacterium vitis and management of crown gall disease in vineyards. Curr. Top Microbiol. Immunol. 418, 15–53 (2018).
PubMed Google Scholar
Chilton, M.-D. et al. Agrobacterium rhizogenes inserts T-DNA into the genomes of the host plant root cells. Nature 295, 432–434 (1982).
Article CAS ADS Google Scholar
Thapa, S. P. et al. Genome-wide analyses of Liberibacter species provides insights into evolution, phylogenetic relationships, and virulence factors. Mol. Plant Pathol. https://doi.org/10.1111/mpp.12925 (2020).
Clark, K. et al. An effector from the Huanglongbing-associated pathogen targets citrus proteases. Nat. Commun. 9, 1718 (2018).
Article PubMed PubMed Central ADS Google Scholar
Li, J. et al. ‘Candidatus Liberibacter asiaticus’ encodes a functional salicylic acid (SA) hydroxylase that degrades SA to suppress plant defenses. Mol. Plant Microbe Interact. 30, 620–630 (2017).
Article CAS PubMed Google Scholar
Kim, J. S., Sagaram, U. S., Burns, J. K., Li, J. L. & Wang, N. Response of sweet orange (Citrus sinensis) to ‘Candidatus Liberibacter asiaticus’ infection: microscopy and microarray analyses. Phytopathology 99, 50–57 (2009).
Article PubMed Google Scholar
Zou, H. et al. The destructive citrus pathogen, ‘Candidatus Liberibacter asiaticus’ encodes a functional flagellin characteristic of a pathogen-associated molecular pattern. PLoS ONE 7, e46447 (2012).
Article CAS PubMed PubMed Central ADS Google Scholar
Pitino, M., Armstrong, C. M. & Duan, Y. Molecular mechanisms behind the accumulation of ATP and H2O2 in citrus plants in response to ‘Candidatus Liberibacter asiaticus’ infection. Hortic. Res. 4, 17040 (2017).
Article PubMed Central PubMed Google Scholar
Luna, E. et al. Callose deposition: a multifaceted plant defense response. Mol. Plant Microbe Interact. 24, 183–193 (2011).
Article CAS PubMed Google Scholar
Galletti, R. et al. The AtrbohD-mediated oxidative burst elicited by oligogalacturonides in Arabidopsis is dispensable for the activation of defense responses effective against Botrytis cinerea. Plant Physiol. 148, 1695–1706 (2008).
Article CAS PubMed PubMed Central Google Scholar
Achor, D. S. et al. Sequence of anatomical symptom observations in citrus affected with Huanglongbing disease. Plant Pathol. J. 9, 56–64 (2010).
Article Google Scholar
Raiol-Junior, L. L., Cifuentes-Arenas, J. C., de Carvalho, E. V., Girardi, E. A. & Lopes, S. A. Evidence that ‘Candidatus Liberibacter asiaticus’ moves predominantly toward new tissue growth in citrus plants. Plant Dis. 105, 34–42 (2021).
Article PubMed Google Scholar
Otero, S. & Helariutta, Y. Companion cells: a diamond in the rough. J. Exp. Bot. 68, 71–78 (2017).
Article CAS PubMed Google Scholar
Kaur, P., Gonzalez, P., Dutt, M. & Etxeberria, E. Identification of sieve elements and companion cell protoplasts by a combination of brightfield and fluorescence microscopy. Appl. Plant Sci. 6, e01179 (2018).
Article PubMed PubMed Central Google Scholar
Koh, E. J. et al. Callose deposition in the phloem plasmodesmata and inhibition of phloem transport in citrus leaves infected with “Candidatus Liberibacter asiaticus”. Protoplasma 249, 687–697 (2012).
Article PubMed Google Scholar
Rong, W., Feng, F., Zhou, J. & He, C. Effector-triggered innate immunity contributes Arabidopsis resistance to Xanthomonas campestris. Mol. Plant Pathol. 11, 783–793 (2010).
CAS PubMed PubMed Central Google Scholar
Wang, Y., Li, X., Fan, B., Zhu, C. & Chen, Z. Regulation and function of defense-related callose deposition in plants. Int. J. Mol. Sci. 22, https://doi.org/10.3390/ijms22052393 (2021).
McLennan, H. R. & Degli Esposti, M. The contribution of mitochondrial respiratory complexes to the production of reactive oxygen species. J. Bioenerg. Biomembr. 32, 153–162 (2000).
Article CAS PubMed Google Scholar
Li, J. et al. Residue dynamics of streptomycin in citrus delivered by foliar spray and trunk injection and effect on ‘Candidatus Liberibacter asiaticus’ titer. Phytopathology, https://doi.org/10.1094/phyto-09-20-0427-r (2020).
Gechev, T. S., Van Breusegem, F., Stone, J. M., Denev, I. & Laloi, C. Reactive oxygen species as signals that modulate plant stress responses and programmed cell death. Bioessays 28, 1091–1101 (2006).
Article CAS PubMed Google Scholar
Montillet, J. L. et al. Fatty acid hydroperoxides and H2O2 in the execution of hypersensitive cell death in tobacco leaves. Plant Physiol. 138, 1516–1526 (2005).
Article CAS PubMed PubMed Central Google Scholar
Kumar, N., Ebel, R. C. & Roberts, P. D. Superoxide dismutase activity in kumquat leaves infected with Xanthomonas axonopodis pv. citri. J. Hortic. Sci. Biotechnol. 86, 62–68 (2011).
Article CAS Google Scholar
Teper, D., Xu, J., Li, J. & Wang, N. The immunity of Meiwa kumquat against Xanthomonas citri is associated with a known susceptibility gene induced by a transcription activator-like effector. PLoS Pathog. 16, e1008886 (2020).
Article CAS PubMed PubMed Central Google Scholar
Torres, M. A., Jones, J. D. & Dangl, J. L. Reactive oxygen species signaling in response to pathogens. Plant Physiol. 141, 373–378 (2006).
Article CAS PubMed PubMed Central Google Scholar
Guénal, I., Sidoti-de Fraisse, C., Gaumer, S. & Mignotte, B. Bcl-2 and Hsp27 act at different levels to suppress programmed cell death. Oncogene 15, 347–360 (1997).
Article PubMed Google Scholar
Claeys, H., Van Landeghem, S., Dubois, M., Maleux, K. & Inzé, D. What is stress? Dose-response effects in commonly used in vitro stress assays. Plant Physiol. 165, 519–527 (2014).
Article CAS PubMed PubMed Central Google Scholar
Torres, M. A., Dangl, J. L. & Jones, J. D. Arabidopsis gp91phox homologues AtrbohD and AtrbohF are required for accumulation of reactive oxygen intermediates in the plant defense response. Proc. Natl Acad. Sci U.S.A. 99, 517–522 (2002).
Article CAS PubMed ADS Google Scholar
Kobayashi, M. et al. Calcium-dependent protein kinases regulate the production of reactive oxygen species by potato NADPH oxidase. Plant Cell 19, 1065–1080 (2007).
Article CAS PubMed PubMed Central Google Scholar
Tang, L. & Vashisth, T. New insight in Huanglongbing-associated mature fruit drop in citrus and its link to oxidative stress. Sci. Hortic. 265, 109246 (2020).
Article CAS Google Scholar
Martinelli, F. & Dandekar, A. M. Genetic mechanisms of the devious intruder candidatus liberibacter in citrus. Front. Plant Sci. 8, 904–904 (2017).
Article PubMed PubMed Central Google Scholar
Levine, D. S. Immune modulating therapies for idiopathic inflammatory bowel diseases. Adv. Pharmacol. 25, 171–234 (1994).
Article CAS PubMed Google Scholar
Wu, D., Lewis, E. D., Pae, M. & Meydani, S. N. Nutritional modulation of immune function: analysis of evidence, mechanisms, and clinical relevance. Front. Immunol. 9, 3160 (2018).
Article CAS PubMed Google Scholar
Campbell, N. K., Fitzgerald, H. K. & Dunne, A. Regulation of inflammation by the antioxidant haem oxygenase 1. Nat. Rev. Immunol. 21, 411–425 (2021).
Article CAS PubMed Google Scholar
Huot, B., Yao, J., Montgomery, B. L. & He, S. Y. Growth-defense tradeoffs in plants: a balancing act to optimize fitness. Mol. Plant 7, 1267–1287 (2014).
Article CAS PubMed PubMed Central Google Scholar
Ma, X. et al. Dual and opposing roles of xanthine dehydrogenase in defense-associated reactive oxygen species metabolism in arabidopsis. Plant Cell. 28, 1108–1126 (2016).
Article CAS PubMed PubMed Central Google Scholar
Chen, S. et al. NtMYB4 and NtCHS1 are critical factors in the regulation of flavonoid biosynthesis and are involved in salinity responsiveness. Front. Plant Sci. 10, 178 (2019).
Article CAS PubMed PubMed Central Google Scholar
Furman, D. et al. Chronic inflammation in the etiology of disease across the life span. Nat. Med. 25, 1822–1832 (2019).
Article CAS PubMed PubMed Central Google Scholar
Kotas, M. E. & Medzhitov, R. Homeostasis, inflammation, and disease susceptibility. Cell 160, 816–827 (2015).
Article CAS PubMed PubMed Central Google Scholar
Zhang, J. et al. A Pseudomonas syringae effector inactivates MAPKs to suppress PAMP-induced immunity in plants. Cell Host Microbe 1, 175–185 (2007).
Article CAS PubMed Google Scholar
Dunand, C., Crèvecoeur, M. & Penel, C. Distribution of superoxide and hydrogen peroxide in Arabidopsis root and their influence on root development: possible interaction with peroxidases. New Phytol. 174, 332–341 (2007).
Article CAS PubMed Google Scholar
Jing, Y. et al. Danger-associated peptide regulates root immune responses and root growth by affecting ROS formation in. Int. J. Mol. Sci. 21, https://doi.org/10.3390/ijms21134590 (2020).
Canales, E. et al. ‘Candidatus Liberibacter asiaticus’,causal agent of citrus huanglongbing, is reduced by treatment with brassinosteroids. PLoS ONE 11, e0146223 (2016).
Article PubMed PubMed Central Google Scholar
Zhang, Z., Lynch, J. P., Zhang, B. & Wang, Q. in Plant Macronutrient Use Efficiency (ed Takehiro Kamiya Mohammad Anwar Hossain, David J. Burritt, Lam-Son Phan Tran, Toru Fujiwara) 245-265 (Academic Press, 2017).
Tavanti, T. R. et al. Micronutrient fertilization enhances ROS scavenging system for alleviation of abiotic stresses in plants. Plant Physiol. Biochem. 160, 386–396 (2021).
Article CAS PubMed Google Scholar
Yuan, X. et al. Region-wide comprehensive implementation of roguing infected trees, tree replacement, and insecticide applications successfully controls citrus HLB. Phytopathology, https://doi.org/10.1094/phyto-09-20-0436-r (2020).
Hijaz, F., Al-Rimawi, F., Manthey, J. A. & Killiny, N. Phenolics, flavonoids and antioxidant capacities in in Citrus species with different degree of tolerance to Huanglongbing. Plant Signal. Behav. 15, 1752447 (2020).
Article PubMed PubMed Central Google Scholar
Pang, Z. et al. Citrus CsACD2 is a target of candidatus liberibacter asiaticus in huanglongbing disease. Plant Physiol. 184, 792–805 (2020).
Article CAS PubMed PubMed Central Google Scholar
Orbović, V. & Grosser, J. W. Citrus transformation using juvenile tissue explants. Methods Mol. Biol. 1224, 245–257 (2015).
Article PubMed Google Scholar
Curtis, I. S., Davey, M. R. & Power, J. B. Leaf disk transformation. Methods Mol. Biol. 44, 59–70 (1995).
CAS PubMed Google Scholar
Clough, S. J. & Bent, A. F. Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 16, 735–743 (1998).
Article CAS PubMed Google Scholar
Vasconcelos, N. C. F., Li, J., Pang, Z., Vincent, C. & Wang, N. The Total Population Size of ‘Candidatus Liberibacter asiaticus’ Inside the Phloem of Citrus Trees and the Corresponding Metabolic Burden Related to Huanglongbing Disease Development. Phytopathology, https://doi.org/10.1094/PHYTO-09-20-0388-R (2020).
Wang, Z. et al. Development and application of molecular-based diagnosis for ‘Candidatus Liberibacter asiaticus’, the causal pathogen of citrus huanglongbing. Plant Pathol. 55, 630–638 (2006).
Article CAS Google Scholar
Trivedi, P. et al. Quantification of viable Candidatus Liberibacter asiaticus in hosts using quantitative PCR with the aid of ethidium monoazide (EMA). Eur. J. Plant Pathol. 124, 553–563 (2009).
Article CAS Google Scholar
Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408 (2001).
Article CAS PubMed Google Scholar
Fernández-Bautista, N., Domínguez-Núñez, J., Moreno, M. & Berrocal-Lobo, M. Plant tissue trypan blue staining during phytopathogen infection. Bio Protoc. 6, e2078 (2016).
Article Google Scholar
Franco, J. Y. et al. Citrus vascular proteomics highlights the role of peroxidases and serine proteases during huanglongbing disease progression. Mol. Cell. Proteomics 19, 1936–1952 (2020).
Article CAS PubMed PubMed Central Google Scholar
Huang, X., Wang, Y., Xu, J. & Wang, N. Development of multiplex genome editing toolkits for citrus with high efficacy in biallelic and homozygous mutations. Plant Mol. Biol. 104, 297–307 (2020).
Article CAS PubMed Google Scholar
Omar, A. A., Dutt, M., Gmitter, F. G. & Grosser, J. W. Somatic embryogenesis: still a relevant technique in citrus improvement. Methods Mol. Biol. 1359, 289–327 (2016).
Article CAS PubMed Google Scholar
Xu, Q. et al. The draft genome of sweet orange (Citrus sinensis). Nat. Genet. 45, 59–66 (2013).
Article CAS PubMed Google Scholar
Kim, D., Paggi, J. M., Park, C., Bennett, C. & Salzberg, S. L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915 (2019).
Article CAS PubMed PubMed Central Google Scholar
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
Article PubMed PubMed Central Google Scholar
Anders, S., Pyl, P. T. & Huber, W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166–169 (2015).
Article CAS PubMed Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Article PubMed PubMed Central Google Scholar
Xu, J. et al. The structure and function of the global citrus rhizosphere microbiome. Nat. Commun. 9, 4894 (2018).
Article PubMed PubMed Central ADS Google Scholar
Pons, E., Peris, J. E. & Peña, L. Field performance of transgenic citrus trees: assessment of the long-term expression of uidA and nptII transgenes and its impact on relevant agronomic and phenotypic characteristics. BMC Biotechnol. 12, 41 (2012).
Article CAS PubMed PubMed Central Google Scholar
Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43, e47 (2015).
Article PubMed PubMed Central Google Scholar
Tian, T. et al. agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 45, W122–w129 (2017).
Article CAS PubMed PubMed Central Google Scholar
Warnes, G. et al. Gplots: various R programming tools for plotting data. R Package. version 2.17.0. http://cran.r-project.org/package=gplots (2015).

Download references

Acknowledgements

The research has been supported by Florida Citrus Initiative, Florida Citrus Research and Development Foundation, USDA National Institute of Food and Agriculture grant # 2018-70016-27412, #2016-70016-24833, and #2019-70016-29796. We thank Dr. Steven E. Lindow for critical reading of the paper.

Author information

These authors contributed equally: Wenxiu Ma, Zhiqian Pang, Xiaoen Huang, Jin Xu, Sheo Shankar Pandey, Jinyun Li, Diann S. Achor.

Authors and Affiliations

Citrus Research and Education Center, Department of Microbiology and Cell Science, IFAS, University of Florida, Lake Alfred, FL, USA
Wenxiu Ma, Zhiqian Pang, Xiaoen Huang, Jin Xu, Sheo Shankar Pandey, Jinyun Li, Diann S. Achor, Fernanda N. C. Vasconcelos, Connor Hendrich, Yixiao Huang, Wenting Wang, Donghwan Lee, Daniel Stanton & Nian Wang

Authors

Wenxiu Ma
View author publications
You can also search for this author in PubMed Google Scholar
Zhiqian Pang
View author publications
You can also search for this author in PubMed Google Scholar
Xiaoen Huang
View author publications
You can also search for this author in PubMed Google Scholar
Jin Xu
View author publications
You can also search for this author in PubMed Google Scholar
Sheo Shankar Pandey
View author publications
You can also search for this author in PubMed Google Scholar
Jinyun Li
View author publications
You can also search for this author in PubMed Google Scholar
Diann S. Achor
View author publications
You can also search for this author in PubMed Google Scholar
Fernanda N. C. Vasconcelos
View author publications
You can also search for this author in PubMed Google Scholar
Connor Hendrich
View author publications
You can also search for this author in PubMed Google Scholar
Yixiao Huang
View author publications
You can also search for this author in PubMed Google Scholar
Wenting Wang
View author publications
You can also search for this author in PubMed Google Scholar
Donghwan Lee
View author publications
You can also search for this author in PubMed Google Scholar
Daniel Stanton
View author publications
You can also search for this author in PubMed Google Scholar
Nian Wang
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

N.W. conceptualized, designed the experiments, and supervised the project. W.X.M. performed ROS related experiments; Z.Q.P. conducted ion leakage experiments; W.X.M. and Z.Q.P. conducted trypan-blue staining together; X.E.H. conducted research related to protoplast; S.S.P. tested H₂O₂ concentrations, callose deposition, and gene expression of young leaves; J.Y.L. conducted foliar spray and trunk injection; J.X. conducted bioinformatic analysis of transcriptomic data; D.A. conducted TEM analyses; F.V. conducted callose deposition analyses of mature leaves; Z.Q.P., F.V., and Y.X.H. conducted overexpression of CLas genes in citrus, tobacco and Arabidopsis; C.H. was involved in gene expression using RT-qPCR; W.T.W., Y.X.H., and J.Y.L. participated in ROS and ion leakage experiments; D.H.L. participated in field trials. D.S. conducted confocal microscopy analysis. N.W. wrote the paper with input from all co-authors.

Corresponding author

Correspondence to Nian Wang.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks Daniel Klessig, Everlyne Muleke and the other, anonymous, reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Data 1

Reporting Summary

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Ma, W., Pang, Z., Huang, X. et al. Citrus Huanglongbing is a pathogen-triggered immune disease that can be mitigated with antioxidants and gibberellin. Nat Commun 13, 529 (2022). https://doi.org/10.1038/s41467-022-28189-9

Download citation

Received: 07 September 2021
Accepted: 04 January 2022
Published: 26 January 2022
DOI: https://doi.org/10.1038/s41467-022-28189-9

This article is cited by

Transcriptome analysis of Citrus sinensis reveals potential responsive events triggered by Candidatus Liberibacter asiaticus
- Chang Liu
- Xiaopeng Chang
- Ruimin Li
Protoplasma (2024)
Controlling citrus Huanglongbing based on soil remediation and biocontrol
- Yang Zhang
- Yi Chen
- Yujian He
European Journal of Plant Pathology (2024)
Humic substances trigger plant immune responses
- Rakiely Martins da Silva
- Natália Aguiar Canellas
- Luciano Pasqualoto Canellas
Chemical and Biological Technologies in Agriculture (2023)
Polyethyleneimine-coated MXene quantum dots improve cotton tolerance to Verticillium dahliae by maintaining ROS homeostasis
- Ping Qiu
- Jiayue Li
- Longfu Zhu
Nature Communications (2023)
The protein interactome of the citrus Huanglongbing pathogen Candidatus Liberibacter asiaticus
- Erica W. Carter
- Orlene Guerra Peraza
- Nian Wang
Nature Communications (2023)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.