Duplicated OsATG9 Genes Antagonize Autophagy to Balance Growth and Drought Tolerance in RiceYiming Li1,†, Yuantai Liu1,†, Mengzhao Shi1,†, Xiaoyun Luo1, Yanshu Huang1, Hao Zeng1, Yunfeng Liu2, Yifeng Huang3, Peng Xu4, Yangwen Qian5, Xixian Li1, Jieying Wang1, Qingjun Xie1,*, Qianying Yang1,*1State Key Laboratory for Conservation and Utilization of Subtropical Agro-Bioresources, Guangdong Basic Research Center of Excellence for Precise Breeding of Future Crops, Key Laboratory for Enhancing Resource Use Efficiency of Crops in South China, Ministry of Agriculture and Rural Affairs, Guangdong Provincial Key Laboratory of Plant Molecular Breeding, College of Agriculture, South China Agricultural University, Guangzhou 510642, China.2State Key Laboratory for Conservation and Utilization of Subtropical Agro-Bioresources, College of Life Sciences and Technology, Guangxi University, Nanning 530004, China.3Institute of Crop and Nuclear Technology Utilization, Zhejiang Academy of Agricultural Science, Hangzhou, 310001, China.4CAS Key Laboratory of Tropical Plant Resources and Sustainable Use, The Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Menglun, Mengla, Yunnan 666303, China.5WIMI Biotechnology Co. Ltd., Changzhou, 213000, China.*Corresponding author: Dr. Qianying Yang (yangqianying@scau.edu.cn) and Dr. Qingjun Xie (qjxie@scau.edu.cn)†These authors contributed equally to this work.Gene duplication drives functional innovation in evolution, enabling organisms to develop new mechanisms to cope with environmental challenges. One such mechanism is autophagy—a conserved intracellular degradation pathway that not only supports normal development but is also strongly induced by abiotic stresses. Indeed, when plants face drought, autophagy is activated to recycle cellular components (Yagyu et al., 2023). Autophagy initiation involves conserved ATG genes (AuTophaGy related genes), with ATG9 facilitating autophagosome formation (Zhuang et al., 2017). While most plants possess a single ATG9 gene, rice experienced a segmental duplication that yielded two paralogs, OsATG9a and OsATG9b, with positive selection (Ka/Ks >1) (Xia et al., 2011). This raises the question of whether these duplicates have evolved distinct functions, particularly in the context of drought-induced autophagy, which remains poorly understood in rice. In our previous study, we observed that OsATG9b plays a role in determining grain size and quality in rice. Mutations of OsATG9b result in smaller grains and increased chalkiness, whereas its overexpression enhances grain size and quality (Liu et al., 2023). Considering OsATG9b is the duplicated gene from OsATG9a and they experienced positive selection, we hypothesized potential divergent functions of these genes during rice development and in response to environmental stressors.Firstly, we constructed the evolution analysis of ATG9s from 21 species, which revealed that OsATG9a and OsATG9b belong to different clusters (Supporting Information S1: Figure 1). The superfamily domain of ATG9 is highly conserved across different species (Supporting Information S1: Figure 2). Using fluorescent‐tagged constructs, we showed that both OsATG9a and OsATG9b localize to the Golgi apparatus and colocalize with each other, as well as with the core autophagy marker OsATG8b (Supporting Information S1: Figure 3-4). Furthermore, OsATG9a share the similar expression pattern with OsATG9b (Supporting Information S1: Figure 5). To assess their functional roles, we then analyzed two newosatg9a single mutants, two osatg9a osatg9b double mutants, and previous published osatg9b line to confirm that, OsATG9a also positively regulates autophagy and affects rice grain development by degrading TGW6 and other agronomic traits (Supporting Information S1: Figure 6-9), mirroring OsATG9b’s roles in rice development (Liu et al., 2023).Given that OsATG9b originated from a duplication of OsATG9a with a Ka/Ks ratio > 1 (Xia et al. , 2011), it is plausible that OsATG9b has evolved distinct functions compared to OsATG9a under abiotic stress. To explore this hypothesis, we subjected WT, osatg9a-1 ,osatg9a-2 , osatg9b , osatg9a osatg9b-1 andosatg9a osatg9b-2 plants to 20% PEG (simulating drought) and soil-drying conditions. Surprisingly, osatg9a-1 andosatg9a-2 exhibited significantly lower survival rates than WT, whereas osatg9b , osatg9a osatg9b-1 and osatg9a osatg9b-2 mutants displayed higher survival rates than WT (Figure 1A, B, D, E, Supporting Information S1: Figure 10). These results suggested that OsATG9a positively regulates drought stress response, while OsATG9b negatively regulated it, potentially epistatic to OsATG9a. We then assessed the stomatal movement, water-loss rates, drought-related physiological parameters, and transcript levels of drought-resistance marker genes in the leaves after 20% PEG treatment (Supporting Information S1: Figure 11-13). Consistent with survival rate,osatg9a mutants were more drought-sensitive, whereasosatg9b and osatg9a osatg9b mutants were more drought-resistant than WT.We further investigate whether the divergent drought‐resistance functions between OsATG9a and OsATG9b are linked to autophagy. Hence, we firstly checked the autophagy function in response to drought stress. Under 20% PEG treatment, the drought marker gene OsERD1 (Early-Responsive to Dehydration stress ) peaked early and then declined by day 7 under 20% PEG (Supporting Information S1: Figure 14A, B). OsATG9a expression similarly rose and then fell, whereasOsATG9b peaked early but remained unchanged after 24 h (Supporting Information S1: Figure 14C, D). NBR1 (Next to BRCA1 gene 1 protein) is a selective autophagy substrate and it also act as cargo receptors for degradation of other substrates, thus NBR1 can be used as the autophagy markers (Marshall & Vierstra, 2018). In WT plants, OsNBR1 protein levels dropped in 24 h under PEG and soared above baseline after 3 days (Supporting Information S1: Figure 14E-H). AllOsATG genes were upregulated at 6 h post-PEG (Supporting Information S1: Figure 15), suggesting an initial induction of autophagy followed by suppression after one day. By employing the same approaches using previously characterized autophagy mutants and overexpression lines (osatg5 , OsATG5-OE , osatg7 , andOsATG8b-OE ) (Gou et al., 2019), we confirmed that autophagy contributes positively to drought tolerance (Supporting Information S1: Figure 16-18).We further assessed NBR1 protein levels and ATG8-PE levels (phosphatidylethanolamine), which is another stringent criterion for assessing autophagic flux (Marshall & Vierstra, 2018) in osatg9 mutants. The NBR1 PEG:CK ratios of osatg9a-1 and osatg9a-2 were higher than in WT (Figure 1G, H), consistent with the observations in osatg5 and osatg7 (Supporting Information S1: Figure 16E-H), indicating a reduction in autophagy activity in osatg9a in response to drought stress. However, the PEG:CK ratios inosatg9b , osatg9a osatg9b-1 and osatg9a osatg9b-2 were significantly lower than in WT (Figure 1K, L, Supporting Information S1: Figure 19), suggesting a much higher induction of autophagy flux in osatg9b and osatg9a osatg9b compared to WT. Using the ATG8 antibody, our immunoblotting analyses revealed that in osatg9 mutants, reduced ATG8-PE/ATG8 ratios under non-stress conditions indicated compromised basal autophagic flux (Figure 1I, J, M, N). Paradoxically, osatg9b mutants exhibited significantly stronger induction of ATG8-PE conjugation after PEG treatment (Figure 1I, J, M, N), demonstrating enhanced autophagic activation.This functional antagonism was further validated through overexpression studies. The survival rate of OsATG9b-OE was significantly lower than WT under PEG treatment (Figure 1C, F), and soil-drying conditions (Supporting Information S1: Figure 20A, B). Under both normal and PEG conditions, the percentages of completely open stomata and partially open stomata in OsATG9b-OE were higher than in WT, opposite to the trend observed in osatg9b (Supporting Information S1: Figure 20C), which led to higher water loss in OsATG9b-OE compared to WT (Supporting Information S1: Figure 20D). Concordantly, OsATG9b overexpression suppressed drought-induced autophagic flux, as evidenced by attenuated NBR1 degradation (Figure 1K, L). The OsATG8-PE: OsATG8 ratio revealed that under normal conditions, the ratio was higher inOsATG9b-OE line compared to WT (Figure 1M, N), indicating increased autophagic flux in the OsATG9b-OE line under normal condition. However, after PEG treatment, the ratio in theOsATG9b-OE line did not significantly differ from that in WT, but was lower than under normal condition (Figure 1M, N), suggesting a suppression of autophagic flux in the OsATG9b-OE line following drought stress. The reciprocal regulation patterns establish that OsATG9a promotes drought resilience through autophagy potentiation like canonical ATG5/7/8b components, while OsATG9b acts as a novel negative regulator by constraining autophagic flux during stress adaptation.Gene duplicates typically undergo three evolutionary trajectories: non-functionalization, neofunctionalization/subfunctionalization (retained for functional innovation) (Lynch and Conery, 2000). Evolutionarily significant duplicates often develop specialized roles, as exemplified by legume CHI1B ’s exclusive nodulation control versus nonfunctional CHI1A (Liu et al., 2024b). Such divergence enables biological systems to expand adaptive capacity while maintaining regulatory equilibrium.Abscisic acid (ABA) is a key phytohormone that regulates plant drought resistance and response (Haverroth et al., 2023; Liu et al., 2024a). Hence, the ABA-mediated drought responses were investigated to determine functional links between OsATG9a/9b and ABA signaling. The expression levels of ABA pathway components, including OsDREB2A (AP2/EREBP transcription factor ), OsLEA3 (late embryogenesis abundant 3 ), OsbZIP23 (bZIP transcription factor ), MYB2 (MYB Domain Protein 2 ) and OsPYL7 (Pyrabactin Resistance-like Abscisic Acid Receptor ) showed PEG-induced upregulation in WT, osatg9a , and osatg9b mutants. Strikingly, osatg9a exhibited attenuated induction of these genes, while osatg9b displayed amplified responses (Figure 1O), demonstrating divergent ABA pathway regulation. In plants, ABA and autophagy form a feedback loop via ATG8‐interacting proteins, transcriptional regulation, and TOR‐mediated phosphorylation to balance stress responses and growth (Gou et al., 2019). OsATG9b may affect the ABA pathway under drought conditions; it is possible that a key factor in this pathway, in turn, regulates autophagic activity or that OsATG9b’s suppression of autophagy leads to the accumulation of undegraded proteins, thereby disrupting the ABA pathway and compromising drought resistance. We further proposed that OsATG9b has evolved as a critical regulator of autophagy that balances growth and drought tolerance in concert with OsATG9a (Figure 1P).