Overexpression of rice gene OsATG8b confers tolerance to nitrogen starvation and increases yield and nitrogen use efficiency (NUE) in Arabidopsis

Rice breeding for low input agriculture

A low-input-based farming system can reduce the adverse effects of modern agriculture through proper utilization of natural resources. Modern varieties often need to improve in low-input settings since they are not adapted to these systems. In addition, rice is one of the most widely cultivated crops worldwide. Enhancing rice performance under a low input system will significantly reduce the environmental concerns related to rice cultivation. Traits that help rice to maintain yield performance under minimum inputs like seedling vigor, appropriate root architecture for nutrient use efficiency should be incorporated into varieties for low input systems through integrated breeding approaches. Genes or QTLs controlling nutrient uptake, nutrient assimilation, nutrient remobilization, and root morphology need to be properly incorporated into the rice breeding pipeline. Also, genes/QTLs controlling suitable rice cultivars for sustainable farming. Since several variables influence performance under low input conditions, conventional breeding techniques make it challenging to work on many traits. However, recent advances in omics technologies have created enormous opportunities for rapidly improving multiple characteristics. This review highlights current research on features pertinent to low-input agriculture and provides an overview of alternative genomics-based breeding strategies for enhancing genetic gain in rice suitable for low-input farming practices.

ScATG8 Gene Cloned from Desert Moss Syntrichia caninervis Exhibits Multiple Stress Tolerance

Syntrichia caninervis is the dominant species of biological soil crust in the desert, including the Gurbantunggut Desert in China. It is widely distributed in drylands and considered to be a new model of vegetative desiccation tolerance moss. Here, we cloned an ATG8 gene from S. caninervis and confirmed its function under multiple abiotic stresses, both in situ and in Physcomitrium patens. The results showed that the ScATG8 gene encoded a protein with a highly conserved ATG8 functional domain. ScATG8 gene was increasingly expressed under different abiotic stresses. Under desiccation stress, the overexpression of ScATG8 enhanced the tolerance of S. caninervis and its ability to scavenge ROS. In addition, ScATG8 overexpression promoted the growth of P. patens under multiple stress conditions. Thus, ScATG8 may be a multifunctional gene, and it plays a critical role in the survival of S. caninervis under various abiotic stresses. Our results provide new insights into the function of ATG8 in enabling desiccation tolerance and open up more possibilities for subsequent plant molecular breeding and the mining of the resistance genes of S. caninervis and other moss species.

Soil bacterial community composition in rice-turtle coculture systems with different planting years

The rice-turtle coculture system is the most special rice-fish integrated farming system. In this study, we selected four paddy fields, including a rice monoculture paddy and three rice-turtle paddies with different planting years, to investigate the soil bacterial community composition with Illumina MiSeq sequencing technology. The results indicated that the contents of soil available nitrogen (AN), soil available phosphorus (AP) and soil organic matter (OM) in 9th year of rice-turtle paddy (RT9) were increased by 5.40%, 51.11% and 23.33% compared with rice monoculture paddy (CK), respectively. Significant differences of Acidobacteria, Desulfobacteria, Crenarchaeota were observed among the different rice farming systems. The relative abundance of Methylomonadaceae, Methylococcaceae and Methylophilaceae in RT9 was significantly higher than that in other treatments. RT9 had significantly lower relative abundance of Acidobacteria, but significantly higher relative abundance of Proteobacteria than other treatments. Redundancy analysis showed that soil AN and AP contents were the major factors influencing the abundance of the dominant microbes, wherein Methylomonadaceae, Methylococcaceae and Methylophilaceae were positively correlated with OM. The findings revealed the rice-turtle coculture system in the 9th year had higher soil nutrients and soil bacterial diversity, but there was also a risk of increasing methane emissions.

Mineral nutrients in plants under changing environments: A road to future food and nutrition security

Plant nutrition is an important aspect that contributes significantly to sustainable agriculture, whereas minerals enrichment in edible source implies global human health; hence, both strategies need to be bridged to ensure "One Health" strategies. Abiotic stress-induced nutritional imbalance impairs plant growth. In this context, we discuss the molecular mechanisms related to the readjustment of nutrient pools for sustained plant growth under harsh conditions, and channeling the minerals to edible source (seeds) to address future nutritional security. This review particularly highlights interventions on (i) the physiological and molecular responses of mineral nutrients in crop plants under stressful environments; (ii) the deployment of breeding and biotechnological strategies for the optimization of nutrient acquisition, their transport, and distribution in plants under changing environments. Furthermore, the present review also infers the recent advancements in breeding and biotechnology-based biofortification approaches for nutrient enhancement in crop plants to optimize yield and grain mineral concentrations under control and stress-prone environments to address food and nutritional security.

Regulation of Tomato Fruit Autophagic Flux and Promotion of Fruit Ripening by the Autophagy-Related Gene SlATG8f

Autophagy is a highly conserved self-degradation process that involves the degradation and recycling of cellular components and organelles. Although the involvement of autophagy in metabolic changes during fruit ripening has been preliminarily demonstrated, the variations in autophagic flux and specific functional roles in tomato fruit ripening remain to be elucidated. In this study, we analyzed the variations in autophagic flux during tomato fruit ripening. The results revealed differential expression of the SlATG8 family members during tomato fruit ripening. Transmission electron microscopy observations and dansylcadaverine (MDC) staining confirmed the presence of autophagy at the cellular level in tomato fruits. Furthermore, the overexpression of SlATG8f induced the formation of autophagosomes, increased autophagic flux within tomato fruits, and effectively enhanced the expression of ATG8 proteins during the color-transition phase of fruit ripening, thus promoting tomato fruit maturation. SlATG8f overexpression also led to the accumulation of vitamin C (VC) and soluble solids while reducing acidity in the fruit. Collectively, our findings highlight the pivotal role of SlATG8f in enhancing tomato fruit ripening, providing insights into the mechanistic involvement of autophagy in this process. This research contributes to a better understanding of the key factors that regulate tomato fruit quality and offers a theoretical basis for tomato variety improvement.

Overexpression of CsATG3a improves tolerance to nitrogen deficiency and increases nitrogen use efficiency in arabidopsis

Nitrogen (N) is a major nutrition element for tea plant. However, application of high levels of N negatively causes environmental problems. Therefore, improved N use efficiency (NUE) of tea plant will be highly desirable and crucial for sustainable tea cultivation. Autophagy plays a central role in N recycling and holds potential to improve N utilization, and many AuTophaGy-related genes (ATGs) are involved in the autophagy process. Here, CsATG3a was identified from Camellia sinensis, and the functions involved in N utilization was characterized in arabidopsis (Arabidopsis thaliana). The transcript level of CsATG3a in tea leaves increases with their maturity. Relative to the wild type (WT) arabidopsis, two CsATG3a-overexpressing (CsATG3a-OE) lines exhibited improved vegetative growth, delayed reproductive stage, and upregulated expression of AtATGs (AtATG3, AtATG5 and AtATG8b) in a low N (LN) hydroponic condition. The expression levels of AtNRT1.1, AtNRT2.1, AtNRT2.2, AtAMT1.1 and AtAMT1.3 for N uptake and transport in roots were all significantly higher in CsATG3a-OE lines compared with those in the WT under LN. Meanwhile, the overexpression of CsATG3a in arabidopsis also increased N and dry matter allocation into both rosette leaves and roots under LN. Additionally, compared with WT, improved HI (harvest index), NHI (N harvest index), NUtE (N utilization efficiency) and NUE (N use efficiency) of CsATG3a-OE lines were further confirmed in a low-N soil cultured experiment. Together, these results concluded that CsATG3a is involved in N recycling and enhances tolerance to LN, indicating that CsATG3a holds potential promise to improve NUE in tea plant.

A dominant gene Ihrl1 is tightly linked to and inhibits the gene Ndhrl1 mediating nitrogen-dependent hypersensitive reaction-like phenotype in wheat

Identification and mapping of an inhibitor of Ndhrl1 mediating nitrogen-dependent hypersensitive reaction-like phenotype in wheat. Hypersensitive reaction-like (HRL) traits are characteristic of spontaneous lesions including yellowish spots, brown spots or white-stripe that appeared randomly and dispersedly on all the leaves in the absence of plant pathogens. Our previous studies have shown that the wheat line P7001 showed an HRL trait at low nitrogen supply, and this trait was controlled by the gene Ndhrl1 (Nitrogen-dependent hypersensitive reaction-like 1). In order to investigate the robustness of the trait expression mediated by Ndhrl1 under different genetic backgrounds, seven genetic populations, with P7001 being the common female parent, were constructed and analyzed. F₁ plants from six of the seven combinations showed HRL trait and Ndhrl1 segregated in a dominant way of HRL: non-HRL = 3:1 in the six populations (F₂). Exceptionally, the F₁ plants of P7001/Fielder combination showed non-HRL trait and HRL trait in the F₂ population showed a contrasting recessive segregation ratio of HRL: non-HRL = 1:3, suggesting Fielder may have another HRL-related gene. Using 55 K SNP array and PCR-based markers, the HRL-related gene in Fielder was mapped to an interval of 5.63-12.91 Mb on the short arm of chromosome 2B with the flanking markers Yzu660R075552 and Yzu660F075941. A recombinant with genomic region of Fielder at Ndhrl1 locus showing HRL trait demonstrated that Fielder also harbored Ndhrl1 same as P7001. Thus, Fielder carries a single dominant suppressor of Ndhrl1, designated as Ihrl1 (Inhibitor of hypersensitive reaction-like). Interestingly, Ihrl1 is tightly linked to Ndhrl1 and may be also involved in nitrogen metabolic and (or) signaling pathways.

Quantitative trait loci mapping reveals important genomic regions controlling root architecture and shoot biomass under nitrogen, phosphorus, and potassium stress in rapeseed (Brassica napus L.)

Plants rely on root systems for nutrient uptake from soils. Marker-assisted selection helps breeders to select desirable root traits for effective nutrient uptake. Here, 12 root and biomass traits were investigated at the seedling stage under low nitrogen (LN), low phosphorus (LP), and low potassium (LK) conditions, respectively, in a recombinant inbred line (RIL) population, which was generated from Brassica napus L. Zhongshuang11 and 4D122 with significant differences in root traits and nutrient efficiency. Significant differences for all the investigated traits were observed among RILs, with high heritabilities (0.43-0.74) and high correlations between the different treatments. Quantitative trait loci (QTL) mapping identified 57, 27, and 36 loci, explaining 4.1-10.9, 4.6-10.8, and 4.9-17.4% phenotypic variances under LN, LP, and LK, respectively. Through QTL-meta analysis, these loci were integrated into 18 significant QTL clusters. Four major QTL clusters involved 25 QTLs that could be repeatedly detected and explained more than 10% phenotypic variances, including two NPK-common and two specific QTL clusters (K and NK-specific), indicating their critical role in cooperative nutrients uptake of N, P, and K. Moreover, 264 genes within the four major QTL clusters having high expressions in roots and SNP/InDel variations between two parents were identified as potential candidate genes. Thirty-eight of them have been reported to be associated with root growth and development and/or nutrient stress tolerance. These key loci and candidate genes lay the foundation for deeper dissection of the NPK starvation response mechanisms in B. napus.

CsATG101 Delays Growth and Accelerates Senescence Response to Low Nitrogen Stress in Arabidopsis thaliana

For tea plants, nitrogen (N) is a foundational element and large quantities of N are required during periods of roundly vigorous growth. However, the fluctuation of N in the tea garden could not always meet the dynamic demand of the tea plants. Autophagy, an intracellular degradation process for materials recycling in eukaryotes, plays an important role in nutrient remobilization upon stressful conditions and leaf senescence. Studies have proven that numerous autophagy-related genes (ATGs) are involved in N utilization efficiency in Arabidopsis thaliana and other species. Here, we identified an ATG gene, CsATG101, and characterized the potential functions in response to N in A. thaliana. The expression patterns of CsATG101 in four categories of aging gradient leaves among 24 tea cultivars indicated that autophagy mainly occurred in mature leaves at a relatively high level. Further, the in planta heterologous expression of CsATG101 in A. thaliana was employed to investigate the response of CsATG101 to low N stress. The results illustrated a delayed transition from vegetative to reproductive growth under normal N conditions, while premature senescence under N deficient conditions in transgenic plants vs. the wild type. The expression profiles of 12 AtATGs confirmed the autophagy process, especially in mature leaves of transgenic plants. Also, the relatively high expression levels for AtAAP1, AtLHT1, AtGLN1;1, and AtNIA1 in mature leaves illustrated that the mature leaves act as the source leaves in transgenic plants. Altogether, the findings demonstrated that CsATG101 is a candidate gene for improving annual fresh tea leaves yield under both deficient and sufficient N conditions via the autophagy process.

Recent Advances in Agronomic and Physio-Molecular Approaches for Improving Nitrogen Use Efficiency in Crop Plants

The efficiency with which plants use nutrients to create biomass and/or grain is determined by the interaction of environmental and plant intrinsic factors. The major macronutrients, especially nitrogen (N), limit plant growth and development (1.5-2% of dry biomass) and have a direct impact on global food supply, fertilizer demand, and concern with environmental health. In the present time, the global consumption of N fertilizer is nearly 120 MT (million tons), and the N efficiency ranges from 25 to 50% of applied N. The dynamic range of ideal internal N concentrations is extremely large, necessitating stringent management to ensure that its requirements are met across various categories of developmental and environmental situations. Furthermore, approximately 60 percent of arable land is mineral deficient and/or mineral toxic around the world. The use of chemical fertilizers adds to the cost of production for the farmers and also increases environmental pollution. Therefore, the present study focused on the advancement in fertilizer approaches, comprising the use of biochar, zeolite, and customized nano and bio-fertilizers which had shown to be effective in improving nitrogen use efficiency (NUE) with lower soil degradation. Consequently, adopting precision farming, crop modeling, and the use of remote sensing technologies such as chlorophyll meters, leaf color charts, etc. assist in reducing the application of N fertilizer. This study also discussed the role of crucial plant attributes such as root structure architecture in improving the uptake and transport of N efficiency. The crosstalk of N with other soil nutrients plays a crucial role in nutrient homeostasis, which is also discussed thoroughly in this analysis. At the end, this review highlights the more efficient and accurate molecular strategies and techniques such as N transporters, transgenes, and omics, which are opening up intriguing possibilities for the detailed investigation of the molecular components that contribute to nitrogen utilization efficiency, thus expanding our knowledge of plant nutrition for future global food security.

Metabolism and autophagy in plants - A perfect match

Autophagy is a eukaryotic cellular transport mechanism that delivers intracellular macromolecules, proteins, and even organelles to a lytic organelle (vacuole in yeast and plants/lysosome in animals) for degradation and nutrient recycling. The process is mediated by highly conserved Autophagy-Related (ATG) proteins. In plants, autophagy maintains cellular homeostasis under favorable conditions, guaranteeing normal plant growth and fitness. Severe stress such as nutrient starvation and plant senescence further induce it, thus ensuring plant survival under unfavorable conditions by providing nutrients through the removal of damaged or aged proteins, or organelles. In this article, we examine the interplay between metabolism and autophagy, focusing on the different aspects of this reciprocal relationship. We show that autophagy has a strong influence on a range of metabolic processes, whereas, at the same time, even single metabolites can activate autophagy. We highlight the involvement of ATG genes in metabolism, examine the role of the macronutrients carbon and nitrogen, as well as various micronutrients, and take a closer look at how the interaction between autophagy and metabolism impacts on plant phenotypes and yield.

Molecular Regulatory Networks for Improving Nitrogen Use Efficiency in Rice

Nitrogen is an important factor limiting the growth and yield of rice. However, the excessive application of nitrogen will lead to water eutrophication and economic costs. To create rice varieties with high nitrogen use efficiency (NUE) has always been an arduous task in rice breeding. The processes for improving NUE include nitrogen uptake, nitrogen transport from root to shoot, nitrogen assimilation, and nitrogen redistribution, with each step being indispensable to the improvement of NUE. Here, we summarize the effects of absorption, transport, and metabolism of nitrate, ammonium, and amino acids on NUE, as well as the role of hormones in improving rice NUE. Our discussion provide insight for further research in the future.

Targeting Nitrogen Metabolism and Transport Processes to Improve Plant Nitrogen Use Efficiency

In agricultural cropping systems, relatively large amounts of nitrogen (N) are applied for plant growth and development, and to achieve high yields. However, with increasing N application, plant N use efficiency generally decreases, which results in losses of N into the environment and subsequently detrimental consequences for both ecosystems and human health. A strategy for reducing N input and environmental losses while maintaining or increasing plant performance is the development of crops that effectively obtain, distribute, and utilize the available N. Generally, N is acquired from the soil in the inorganic forms of nitrate or ammonium and assimilated in roots or leaves as amino acids. The amino acids may be used within the source organs, but they are also the principal N compounds transported from source to sink in support of metabolism and growth. N uptake, synthesis of amino acids, and their partitioning within sources and toward sinks, as well as N utilization within sinks represent potential bottlenecks in the effective use of N for vegetative and reproductive growth. This review addresses recent discoveries in N metabolism and transport and their relevance for improving N use efficiency under high and low N conditions.

Genetic regulation of the traits contributing to wheat nitrogen use efficiency

High nitrogen application aimed at increasing crop yield is offset by higher production costs and negative environmental consequences. For wheat, only one third of the applied nitrogen is utilized, which indicates there is scope for increasing Nitrogen Use Efficiency (NUE). However, achieving greater NUE is challenged by the complexity of the trait, which comprises processes associated with nitrogen uptake, transport, reduction, assimilation, translocation and remobilization. Thus, knowledge of the genetic regulation of these processes is critical in increasing NUE. Although primary nitrogen uptake and metabolism-related genes have been well studied, the relative influence of each towards NUE is not fully understood. Recent attention has focused on engineering transcription factors and identification of miRNAs acting on expression of specific genes related to NUE. Knowledge obtained from model species needs to be translated into wheat using recently-released whole genome sequences, and by exploring genetic variations of NUE-related traits in wild relatives and ancient germplasm. Recent findings indicate the genetic basis of NUE is complex. Pyramiding various genes will be the most effective approach to achieve a satisfactory level of NUE in the field.

Autophagy mediates grain yield and nitrogen stress resistance by modulating nitrogen remobilization in rice

Autophagy, a conserved cellular process in eukaryotes, has evolved to a sophisticated process to dispose of intracellular constituents and plays important roles in plant development, metabolism, and efficient nutrients remobilization under suboptimal nutrients conditions. Here, we show that OsATG8b, an AUTOPHAGY-RELATED8 (ATG8) gene in rice, was highly induced by nitrogen (N) starvation. Elevated expression of OsATG8b significantly increased ATG8 lipidation, autophagic flux, and grain yield in rice under both sufficient and deficient N conditions. Overexpressing of OsATG8b could greatly increase the activities of enzymes related to N metabolism. Intriguingly, the 15N-labeling assay further revealed that more N was remobilized to seeds in OsATG8b-overexpressing rice, which significantly increased the N remobilization efficiency (NRE), N harvest index, N utilization efficiency (NUE), and N uptake efficiency (NUpE). Conversely, the osatg8b knock-out mutants had the opposite results on these characters. The substantial transcriptional changes of the overexpressed transgenic lines indicated the presence of complex signaling to developmental, metabolic process, and hormone, etc. Excitingly, the transgenic rice under different backgrounds all similarly be boosted in yield and NUE with OsATG8b overexpression. This work provides an excellent candidate gene for improving N remobilization, utilization, and yield in crops simultaneously.

Signaling Responses to N Starvation: Focusing on Wheat and Filling the Putative Gaps With Findings Obtained in Other Plants. A Review

Wheat is one of the most important food crops worldwide. In recent decades, fertilizers, especially nitrogen (N), have been increasingly utilized to maximize wheat productivity. However, a large proportion of N is not used by plants and is in fact lost into the environment and causes serious environmental pollution. Therefore, achieving a low N optimum via efficient physiological and biochemical processes in wheat grown under low-N conditions is highly important for agricultural sustainability. Although N stress-related N capture in wheat has become a heavily researched subject, how this plant adapts and responds to N starvation has not been fully elucidated. This review summarizes the current knowledge on the signaling mechanisms activated in wheat plants in response to N starvation. Furthermore, we filled the putative gaps on this subject with findings obtained in other plants, primarily rice, maize, and Arabidopsis. Phytohormones have been determined to play essential roles in sensing environmental N starvation and transducing this signal into an adjustment of N transporters and phenotypic adaptation. The critical roles played by protein kinases and critical kinases and phosphatases, such as MAPK and PP2C, as well as the multifaceted functions of transcription factors, such as NF-Y, MYB, DOF, and WRKY, in regulating the expression levels of their target genes (proteins) for low-N tolerance are also discussed. Optimization of root system architecture (RSA) via root branching and thinning, improvement of N acquisition and assimilation, and fine-tuned autophagy are pivotal strategies by which plants respond to N starvation. In light of these findings, we attempted to construct regulatory networks for RSA modification and N uptake, transport, assimilation, and remobilization.

Genome-Wide Identification of CsATGs in Tea Plant and the Involvement of CsATG8e in Nitrogen Utilization

Nitrogen (N) is a macroelement with an indispensable role in the growth and development of plants, and tea plant (Camellia sinensis) is an evergreen perennial woody species with young shoots for harvest. During senescence or upon N stress, autophagy has been shown to be induced in leaves, involving a variety of autophagy-related genes (ATGs), which have not been characterized in tea plant yet. In this study, a genome-wide survey in tea plant genome identified a total of 80 Camellia Sinensis autophagy-related genes, CsATGs. The expression of CsATG8s in the tea plant showed an obvious increase from S1 (stage 1) to S4 (stage 4), especially for CsATG8e. The expression levels of AtATGs (Arabidopsis thaliana) and genes involved in N transport and assimilation were greatly improved in CsATG8e-overexpressed Arabidopsis. Compared with wild type, the overexpression plants showed earlier bolting, an increase in amino N content, as well as a decrease in biomass and the levels of N, phosphorus and potassium. However, the N level was found significantly higher in APER (aerial part excluding rosette) in the overexpression plants relative to wild type. All these results demonstrated a convincing function of CsATG8e in N remobilization and plant development, indicating CsATG8e as a potential gene for modifying plant nutrient utilization.

Autophagy Dances with Phytohormones upon Multiple Stresses

Autophagy is an evolutionarily conserved process for turning over unwanted cellular components, thus promoting nutrient recycling and maintaining cellular homeostasis, which eventually enables plants to survive unfavorable growth conditions. In addition to plant growth and development, previous studies have demonstrated that autophagy is involved in the responses to various environmental challenges through interplaying with multiple phytohormones, including abscisic acid (ABA), jasmonic acid (JA), and salicylic acid (SA). In this review, we summarize the advances made in their synergistic interactions in response to multiple abiotic and biotic stresses; we also discuss the remaining issues and perspectives regarding their crosstalk.

A Rice Autophagy Gene OsATG8b Is Involved in Nitrogen Remobilization and Control of Grain Quality

Enhancing nitrogen (N) use efficiency is a potential way to reduce excessive nitrogen application and increase yield. Autophagy is a conserved degradation system in the evolution of eukaryotic cells and plays an important role in plant development and stress response. Autophagic cores have two conjugation pathways that attach the product of autophagy-related gene 8 (ATG8) to phosphatidylethanolamine (PE) and ATG5 to ATG12, respectively, which then help with vesicle elongation and enclosure. Rice has six ATG8 genes, which have not been functionally confirmed so far. We identified the rice gene OsATG8b and characterized its role in N remobilization to affect grain quality by generating transgenic plants with its over-expression and knockdown. Our study confirmed the autophagy activity of OsATG8b through the complementation of the yeast autophagy-defective mutant scatg8 and by observation of autophagosome formation in rice. The autophagy activity is higher in OsATG8b-OE lines and lower in OsATG8b-RNAi than that in wild type (ZH11). ¹⁵N pulse-chase analysis revealed that OsATG8b-OE plants conferred higher N recycling efficiency to grains, while OsATG8b-RNAi transgenic plants exhibited lower N recycling efficiency and poorer grain quality. The autophagic role of OsATG8b was experimentally confirmed, and it was concluded that OsATG8b-mediated autophagy is involved in N recycling to grains and contributes to the grain quality, indicating that OsATG8b may be a potential gene for molecular breeding and cultivation of rice.

Impact of Autophagy on Gene Expression and Tapetal Programmed Cell Death During Pollen Development in Rice

Autophagy has recently been shown to be required for tapetal programmed cell death (PCD) and pollen maturation in rice. A transcriptional regulatory network is also known to play a key role in the progression of tapetal PCD. However, the relationship between the gene regulatory network and autophagy in rice anther development is mostly unknown. Here, we comprehensively analyzed the effect of autophagy disruption on gene expression profile during the tapetal PCD in rice anther development using high-throughput RNA sequencing. Expression of thousands of genes, including specific transcription factors and several proteases required for tapetal degradation, fluctuated synchronously at specific stages during tapetal PCD progression in the wild-type anthers, while this fluctuation showed significant delay in the autophagy-deficient mutant Osatg7-1. Moreover, gene ontology enrichment analysis in combination with self-organizing map clustering as well as pathway analysis revealed that the expression patterns of a variety of organelle-related genes as well as genes involved in carbohydrate/lipid metabolism were affected in the Osatg7-1 mutant during pollen maturation. These results suggest that autophagy is required for proper regulation of gene expression and quality control of organelles and timely progression of tapetal PCD during rice pollen development.