Co-Overexpression of OsNAR2.1 and OsNRT2.3a Increased Agronomic Nitrogen Use Efficiency in Transgenic Rice Plants

The NO3 - transporter plays an important role in rice nitrogen acquisition and nitrogen-use efficiency. Our previous studies have shown that the high affinity systems for nitrate uptake in rice is mediated by a two-component NRT2/NAR2 transport system. In this study, transgenic plants were successful developed by overexpression of OsNAR2.1 alone, OsNRT2.3a alone and co-overexpression of OsNAR2.1 and OsNRT2.3a. Our field experiments indicated that transgenic lines expressing p35S:OsNAR2.1 or p35S:OsNAR2.1-p35S:OsNRT2.3a constructs exhibited increased grain yields of approximately 14.1% and 24.6% compared with wild-type (cv. Wuyunjing 7, WT) plants, and the agricultural nitrogen use efficiency increased by 15.8% and 28.6%, respectively. Compared with WT, the 15N influx in roots of p35S:OsNAR2.1 and p35S: OsNAR2.1-p35S:OsNRT2.3a lines increased 18.9%‑27.8% in response to 0.2 mM, 2.5 mM 15NO3 –, and 1.25 mM 15NH4 15NO3, while there was no significant difference between p35S:OsNAR2.1 and p35S:OsNAR2.1-p35S:OsNRT2.3a lines; only the 15N distribution ratio of shoot to root for p35S:OsNAR2.1-p35S:OsNRT2.3a lines increased significantly. However, there were no significant differences in nitrogen use efficiency, 15N influx in roots and the yield between the p35S:NRT2.3a transgenic lines and WT. This study indicated that co-overexpression of OsNAR2.1 and OsNRT2.3a could increase rice yield and nitrogen use efficiency.


INTRODUCTION
Nitrogen (N) is an essential macronutrient for plant growth and crop productivity, it affects all levels of plant function from metabolism to resource allocation, growth, and development (Crawford, 1995;Scheible et al., 2004). Nitrate nitrogen (NO 3 -) and ammonium nitrogen (NH 4 + ) are two main inorganic nitrogen sources in plant growth. Conventionally cultivated rice is submerged, and nitrification is inhibited, NH 4 + is the main inorganic nitrogen in rhizosphere soil, therefore, rice is generally considered as ammonium-preferring plant (Arth et al., 1998;Kronzucker et al., 1999). However, for its well-developed aerial tissue, rice can transport and secrete oxygen from photosynthesis above ground to the rhizosphere (Aurelio et al., 2003). Oxygen can stimulate the growth and reproduction of nitrifying bacteria in the rhizosphere, so the part of NH 4 + in the rhizosphere can be nitrified into NO 3 - (Li et al., 2008). Therefore, in the actual growth of rice, its root system has been in the mixed nutrition of NH 4 + and NO 3 - (Kirk and Kronzucker, 2005). Aurelio et al. (2003) found that most rice varieties could absorb the same amount of NH 4 + and NO 3 -, and the amount of NO 3 absorbed by rice would be much larger than that of NH 4 + because of the strong nitrification in rhizosphere. During the late stage of rice growth and development, when it has been going through the process of long-term alternate wetting and drying irrigation or the whole growth stage of upland rice, the uptake of NO 3 was also higher than that of NH 4 + by rice (Arth et al., 1998;Wang et al., 2004).
In addition to be a nutrient, NO 3 also serves as a signaling molecule, which induces changes in plant growth and gene expression (Kirk and Kronzucker, 2005;Liu et al., 2015). For example, NO 3 could break seed dormancy (Alboresi et al., 2005), and induced multiple genes for plant growth and development (O'Brien et al., 2016). Furthermore, the development of lateral root, leaf growth and altering flowering time could be regulated by NO 3 - (Zhang et al., 1999;Rahayu et al., 2005;Castro et al., 2011;Hsu and Tsay, 2013;Naz et al., 2019). NO 3 could be sensed by the plant, and the uptake of NO 3 -, especially in rice, can increase the uptake of NH 4 + to a certain extent (Duan et al., 2006;Chen et al., 2017).
Two different NO 3 uptake systems in plants, the high-affinity NO 3 uptake systems (HATS) and low-affinity NO 3 uptake systems (LATS) are regulated by NO 3 supply and enable plants to cope with low or high NO 3 concentrations in soils (Miller et al., 2007;Xu et al., 2012). As we known, the NRT2 and NPF families contribute to HATS and LATS responding the NO 3 uptake and translocation in plants (Miller et al., 2007;Leŕan et al., 2014). Some high-affinity NO 3 transporters belonging to the NRT2 family have been reported to require a partner protein, NAR2, for their function (Tong et al., 2005;Yong et al., 2010;Xu et al., 2012). In rice, twocomponent NRT2-NAR2 system also exists in the NO 3 transport process, and OsNAR2.1 is the essential partner of OsNRT2 nitrate transporters for nitrate uptake over low and high concentration range (Miller et al., 2007;Feng et al., 2011;Liu et al., 2015). OsNRT2.1, OsNRT2.2, and OsNRT2.3a were similarly shown to require OsNAR2.1 for NO 3 uptake, and their interaction at the protein level was demonstrated by using a yeast two hybrid assay, an oocyte expression system and western blotting Yan et al., 2011;Liu et al., 2014).

qRT-PCR and Southern Blot Analysis
The total RNA extraction of plants and genes expression analysis were performed as described previously (Chen et al., 2016;Chen et al., 2017). The primers for PCR are shown in Table S2.
The Southern blot was used to identify the copy number of inserted T-DNA. The genomic DNA extraction of the T2 plant, DNA digestion and hybridization were following the reported paper (Chen et al., 2017).

Field Experiments for Harvest Yield
We selected the positive T1 transgenic lines by PCR. The sequence of the resistance screening gene (792 bp) on the vector was amplified by forward primer (5' ATGATTGAACAAGATGGAT TGCA 3') and reverse primer (5' GAAGAACTCGTCAAGA AGGCGAT 3'). The Southern blot was used to identify the copy number of inserted T-DNA of the T2 plant, we chose the singlecopy T2 lines and planted them to screen T3 and T4 generation plants with stable inheritance. T1, T2, and T4 generation plants were grown in plots at the Nanjing Agricultural University from June to October in 2013October in , 2014October in , and 2016 in Nanjing, Jiangsu. The rice plants of T3 generation plants were cultivated in plots at the Experimental Site of Zhejiang University June to October 2015 in Changxing, Zhejiang. Soil properties in Nanjing field experiment was described as before (Chen et al., 2016;Chen et al., 2017). Basal applications of 30 kg P/ha as Ca(H2PO4)2 and 60 kg/K ha (KCl) were made to all plots 3 days before transplanting. N fertilizer as urea accounted for 40%, 30% and 40% of the total N fertilizer was applied prior to transplanting, at tillering, just before the heading stage, respectively. T1-T2 generation transgenic plants and wild-type (cv. Wuyunjing 7, WT) plants were planted in the plots with 300 kg N/ha. T3 generation transgenic plants and WT plants were planted in 3 plots with 150 kg N/ha and the plots without nitrogen fertilizer as blank control. T4 generation transgenic plants and WT plants were planted in 3 plots with 300 kg N/ha and the plots without nitrogen fertilizer as blank control. T2 generation transgenic plants of cv. Nipponbare and wild-type (cv. Nipponbare, NP) plants were planted in 3 plots with 150 kg N/ha and the plots without nitrogen fertilizer as blank control. The plots size was 2 × 2 m, and the seedlings were planted in a 10 × 10 array. Compared with WT, the flowering period of p35S: OsNAR2.1 lines were 4-5 days earlier, that of p35S:OsNRT2.3a lines had no significant change, and that of p35S:OsNAR2.1-p35S: OsNRT2.3a lines were 2-3 days earlier. During rice flowering and mature stages, the samples of T4 generation transgenic plants and WT were collected for further analysis.

Dry Weight, Total Nitrogen Measurement, and Calculation of Nitrogen Use Efficiency
We conducted biomass and nitrogen analysis of the T4 generation shoot samples grown under 300 kg N/ha fertilizer condition according to our previous reported method (Chen et al., 2016). Dry matter at anthesis stage (DMA), grain yield (GY), dry matter at maturity stage (DMM), total nitrogen accumulation at anthesis stage (TNAA), grain nitrogen accumulation at maturity stage (GNAM), total nitrogen accumulation at maturity stage (TNAM), harvest index (HI), the contribution of pre-anthesis assimilates to grain yield (CPAGY), dry matter translocation (DMT), dry matter translocation efficiency (DMTE), contribution of pre-anthesis nitrogen to grain nitrogen accumulation (CPNGN), nitrogen translocation (NT), nitrogen translocation efficiency (NTE), nitrogen harvest index (NHI), agronomic nitrogen use efficiency (ANUE), physiological nitrogen use efficiency (PNUE), and nitrogen recovery efficiency (NRE) were calculated according to Chen et al. (2016;. CPAGY For root 15 N uptake experiment, rice seedlings were grown in 1mM NH 4 + for 3 weeks and then under nitrogen starvation condition for 1 week before 15 N uptake experiment. The plants were rinsed in 0.1 mM CaSO 4 for 1 min, then transferred to a solution containing either 0.2 mM 15 NO 3 -, 2.5 mM 15 NO 3 or 1.25 mM 15 NH 4 15 NO 3 (atom % 15 N: 99%) for 5 min, and finally rinsed again in 0.1 mM CaSO 4 for 1 min.
For root 15 N distribution ratio experiment, WT and transgenic seedlings were grown in 1 mM NH 4 + for 3 weeks, nitrogen starvation was then carried out for 1 week. The plants were rinsed in 0.1 mM CaSO 4 for 1 min, then transferred to the nutrient solution containing 0.2 mM 15 NO 3 -, or 1.25 mM 15 NH 4 15 NO 3 for 12 h, and finally rinsed again in 0.1 mM CaSO 4 for 1 min. The 15 N influx rate and 15 N distribution ratio were calculated following the method reported in Chen et al. (2017) and Tang et al. (2012).

Statistical Analysis
For our data statistical analysis, the single factor analysis of variance (ANOVA) and Tukey's test data analysis was applied. All statistical evaluations were conducted using the IBM SPSS Statistics ver. 20 software. (SPSS Inc., Chicago, IL, United States).

Construction and Field Traits Analysis
In this study, we further introduced the p35S:OsNAR2.1 ( Figure  S1A), p35S:OsNRT2.3a ( Figure S1B) and p35S:OsNAR2.1-p35S: OsNRT2.3a ( Figure S1C) expression construction into Wuyunjing7 (WT) which is a high yield rice cultivar used in Jiangsu, China. We generated 18 lines, including 6 p35S:OsNAR2.1 lines (N lines), 6 p35S:OsNRT2.3a lines (A lines), and 6 p35S:OsNAR2.1-p35S:OsNRT2.3a lines (NA lines). Then we conducted field experiments and analyzed the grain yield and biomass of these transgenic lines. T1-T2 generation transgenic plants and WT plants were planted in the plots with 300 kg N/ha in Nanjing, Jiangsu. Compared to WT plants, the yield and biomass for N lines of T1 plants increased by approximately 17.4% and 7.1%, respectively ( Figure S2A). However, there was no significant difference between A lines and WT ( Figure S2B). The yield and biomass of NA lines of T1 plants increased by approximately 30.5% and 19.5%, respectively ( Figure S2C). Based on RNA expression data for the T1 generation (Figures S2D-F) and T2 generation ( Figures 1D-F), we selected three N T2 lines named N1, N2, N3 ( Figure 1A), three A T2 lines named A1, A2, A3 ( Figure 1B), and three NA T2 lines named NA1, NA2, NA3 ( Figure 1C) for further analysis. Southern blot analysis showed that N1 and N2 were two independent transgenic lines with one transgenic copy respectively ( Figure 1G); A1 and A2 were two independent transgenic lines with one transgenic copy respectively ( Figure 1H); NA1 and NA2 were two independent transgenic lines with one transgenic copy respectively ( Figure 1I).
For the T3 transgenic plants were planted in the plots with 150 kg N/ha in Changxing, Zhejiang. With WT as control, the grain yield of N (N1 and N2) lines was increased by 14.5%, and agronomic nitrogen use efficiency increased by 15.1% ( Figure  S3). The grain yield and agronomic nitrogen use efficiency of NA (NA1 and NA2) lines increased by 27.1% and 31.4% compared with WT ( Figure S3). There was no significant difference in grain yield and agronomic nitrogen use efficiency between A lines and WT ( Figure S3). The data of the field experiment for T4 generation which planted in the plots with 300 kg N/ha in Nanjing, Jiangsu Province were analyzed in detail. Compared with WT, the total tiller number per plant and seed setting rate increased by 26.2% and 16.3% respectively, while the panicle length, grain weight and grain number per panicle was reduced by 4.9%, 11.4%, and 11.6% respectively, and the grain yield increased by 14.1% (Table 1). For the NA lines, there was no significant difference in seed setting rate and tiller number between and N lines, but the grain yield increased by 24.6% (Table 1) compared to the control. However, there was no significant difference between all agronomic traits of A lines and WT (Table 1).

Dry Matter Accumulation and N Analysis in Transgenic Lines
Anthesis stage is a particular time from vegetative growth to reproductive growth of rice, and dry matter and nitrogen accumulation during anthesis stage have positive influence on the formation of final yield of rice (Good et al., 2004). We sampled shoot tissues at the anthesis stage and at the mature stage separately to determine the biomass and the total N content. Analysis using WT as control, we found that at the anthesis stage, the biomass of panicles, leaves and culms in N lines increased by 18.7%, 18.3%, and 20.0% respectively; the biomass of panicles, leaves and culms in A lines had no significant difference; the biomass of panicles, leaves and culms in NA lines increased by 21.1%, 19.3%, and 26.6% ( Figure 2B). At the mature stage, the biomass of panicles in N lines increased by about 14.0%, but the biomass of leaves and culms for N was not significant. There was no significant difference between the biomass of panicles, leaves and culms between A lines and WT. For the NA lines, the biomass of panicles, leaves and culms was increased by 26.4%, 15.5%, and 12.8% respectively compared with WT ( Figure 2C).
At anthesis stage, there was no significant difference of the total nitrogen concentration in leaves between N lines and NA lines, but increased by about 8.8% compared with WT. For the total nitrogen concentration in panicles and culms, there was no significant difference between N lines, NA lines and WT. There was no significant difference in total nitrogen concentration between different parts of A lines and WT ( Figure 2D). At maturity stage, there was no significant difference in total nitrogen concentration between N, A and NA lines and WT ( Figure 2E).
The variation of biomass and total nitrogen concentration eventually led to the difference of total nitrogen content per plant. At anthesis stage, there was no significant difference in total nitrogen concentration in panicles, leaves and culms between N and NA lines, but increased by about 18.3%, 29.2%, and 24.7% compared to WT ( Figure 2F). At maturity stage, the total nitrogen content in panicles of N lines increased by about 14.2%, there was no significant difference in the total nitrogen content in leaves and culms between N lines and WT. The total nitrogen content in panicles, leaves and culms of NA lines increased by about 23.0%, 17.7%, and 15.7% respectively, however, there was no significant difference in total nitrogen content between A lines and WT in different plant parts ( Figure 2G).
Based on a method described by Chen et al. (2016), we calculated the physiological nitrogen use efficiency (PNUE),  Figures 3A, B).
The Expression of OsNAR2.1 and OsNRT2s in Transgenic Lines in Different Growth Stages of WT and Transgenic Lines Chen et al. (2016Chen et al. ( , 2017 reported that the expression of OsNAR2.1 and OsNRT2.1 in culms significantly affected rice yield and nitrogen use efficiency. In this study, we analyzed the expression of OsNAR2.1 and OsNRT2s in culms of WT and transgenic lines at 30 days (the vegetative growth period), 60 days (the anthesis stage), and 75 days (the grain filling stage) separately. During these three periods, the expressions of  Dry matter at anthesis (DMA), dry matter at maturity (DMM), grain yield (GY), total nitrogen accumulation at anthesis (TNAA), total nitrogen accumulation at maturity (TNAM), and grain nitrogen accumulation at maturity (GNAM). Statistical analysis of data from T4 generation; n = 3 plots for each mean. The different letters indicate a significant difference between the transgenic line and the WT. (P < 0.05, one-way ANOVA).
OsNAR2.1 and OsNRT2.3a of N (N1, N2, and N3) or NA (N1, NA2, and N3) lines were significantly higher compared with WT. But there was no significant difference of the expressions of OsNRT2.1and OsNRT2.2 between N and NA lines ( Figure  4). Compared with WT, the expression of OsNAR2.1 increased by about 14 times in N lines, while OsNRT2.3a increased by about 51% in those lines. The expression of OsNAR2.1 and OsNRT2.3a in NA lines increased about 12-fold and 10-fold respectively ( Figure 4). Compared with WT, the expression of OsNRT2.3a in A (A1, A2, and A3) lines increased about 10   times, while the expression of OsNAR2.1, OsNRT2.1, and OsNRT2.2 did not change significantly (Figure 4).

Agronomic Nitrogen Use Efficiency Estimation for Co-Overexpression of
OsNAR2.1 and OsNRT2.3a in the cv.

The Plant Seedling Growth and Total Nitrogen Content Evaluation in Transgenic Lines
We further analyzed the growth and nitrogen uptake of transgenic lines at seedling stage. WT and transgenic rice seedlings were grown in IRRI solution containing 1 mM NH 4 + for 2 weeks and then transferred into 0.2 mM NO 3 -, 2.5 mM NO 3 -, or 1.25 mM NH 4 NO 3 for 3 additional weeks ( Figure S4). Compared with WT, under 0.2 mM NO 3 treatment, the biomass of root and shoot in N lines increased 151.2% and 102.7% respectively, it increased 204.0% and 150.7% respectively in NA lines ( Figure 6A). Under 2.5 mM NO 3 treatment, the biomass of root and shoot in N lines increased 136.8% and 142.8% respectively, it increased 170.1% and 185.7% respectively in NA lines ( Figure 6B). Under 1.25 mM NH 4 NO 3 treatment, the biomass of root shoot in N lines increased 58.9% and 56.2% respectively, it increased 85.5% and 102.1% respectively in NA lines ( Figure 6C).
Moreover, we measured the total nitrogen concentration, and found that there was no significant difference in total nitrogen concentration between N lines and NA lines under different treatments. Compared with WT, the total nitrogen concentration in root and shoot of N lines or NA lines increased by 30.5% and 26.6% respectively under 0.2 mM NO 3 treatment ( Figure 6D), increased by 37.4% and 27.4% respectively under 2.5 mM NO 3 treatment ( Figure 6E), increased by 21.1% and 18.0% respectively under 1.25 mM NH 4 NO 3 treatment ( Figure 6F). The total nitrogen content per plant was also calculated. After 0.2 mM NO 3 treatment, the total nitrogen content of the root and shoot in N lines increased 223.5% and 151.2% respectively, it increased 302.0% and 224.0% respectively in NA lines ( Figure  6G). After 2.5 mM NO 3 treatment, the total nitrogen content increased 221.7% and 201.2% respectively in N lines, it increased 275.0% and 273.9% respectively NA lines ( Figure 6H); under 1.25 mM NH 4 NO 3 treatment, the total nitrogen content increased 89.1% and 81.3% respectively in N lines, it increased 128.7% and 142.7% respectively in NA lines ( Figure 6I).
We further analyzed the expression of OsNAR2s and OsNRT2s in transgenic lines treated with 0.2 mM NO 3 and 1.25 mM NH 4 NO 3 . The expression patterns of the transgenic lines were similar under the treatment of 0.2 mM NO 3 and 1.25 mM NH 4 NO 3 . In the N lines, the expression of OsNAR2.1 increased by about 12.4 times compared with WT, and the expression of OsNRT2.1, OsNRT2.2, and OsNRT2.3a increased by 55%-98% ( Figures S5A, D). In the A lines, the expression of OsNRT2.3a in increased about 11 times, but the expression of OsNAR2.1, OsNRT2.2, and OsNRT2.2 did not change significantly ( Figures S5B, E). However, for the NA lines, the expression of OsNAR2.1 increased by 13.8 times, the expression of OsNRT2.3a increased by 10.1 times, and the expression of OsNRT2.1 and OsNRT2.2 increased by 53%-120% ( Figures S5C, F). The expression of OsNAR2.2, OsNRT2.3b and OsNRT2.4 did not change significantly in N, A, and NA lines ( Figure S6).

N Influx Rates and N Distribution Ratio Influx Determination
We analyzed short-term 15 N uptake in same-size seedlings of the transgenic lines and WT which exposed to 0.2 mM 15 NO 3 -, 2.5 mM 15 NO 3 -, or 1.25 mM 15 NH 4 15 NO 3 for 5 min. There was no significant difference of 15 N influx rate between N lines and A lines. Compared with WT, the influx rate of 15 NO 3 increased 27.6% and 20.1% in response to 0.2 mM 15 NO 3 and 2.5 mM   (Figures 7A-C). The influx rate of 15 N did not change compared with that of WT in the A lines ( Figures 7A-C).
We also detected 15 N distribution ratio experiment in root. WT and transgenic seedlings were grown in 1 mM NH 4 + for 3 weeks and treated in nitrogen starvation condition for 1 week. 15 N concentration experiment was conducted under the nutrient solution containing 0.2 mM 15 NO 3 or 1.25 mM 15 NH 4 15 NO 3 for 12 h. Under 0.2mM 15 NO 3 treatment, there was no significant difference of the 15 N concentration in roots between N and NA lines, but the 15 N concentration increased by 19.7% compared with WT. 15 N concentration in roots of N and NA lines increased by 25.5% and 48.1%, respectively compared with WT ( Figure  8A). Under 1.25 mM 15 NH 4 15 NO 3 treatment, there was no significant difference in 15 N concentration in roots between N and NA lines, but it increased by 20.0% compared with WT. 15 N concentration in roots of N and NA lines increased by 21.2% and 33.2%, respectively ( Figure 8B). Under the treatment of 0.

DISCUSSION
Increasing NO 3 uptake is an effective way to improve yield and nitrogen use efficiency in rice. Fang et al. (2013) found that overexpression of nitrate transporter OsPTR9 could promote lateral root formation, increase grain yield of rice. Overexpression of OsNRT1.1B, a low affinity NO 3 transporter protein gene, increased rice nitrogen use efficiency by about 30% (Hu et al., 2015). Overexpression of high affinity nitrate transporter gene OsNRT2.3b increased the buffer capacity of cell pH and significantly increased grain yield and nitrogen use efficiency of rice (Fan et al., 2016). Overexpression of nitrate transporter gene OsNRT1.1A in rice greatly improved nitrogen utilization and grain yield, and maturation time was also significantly shortened (Wang et al., 2018). Gao et al. (2019) reported that overexpression of OsNR2 increased the activity of nitrate reductase and the uptake of nitrate by rice, thus increasing the yield and nitrogen use efficiency of rice. Overexpression of OsNAR2.1 or OsNRT2.1 by the native promoter of OsNAR2.1 can improve grain yield and nitrogen use efficiency of rice (Chen et al., 2016;Chen et al., 2017;Luo et al., 2018;Chen et al., 2019). OsNAC42 as a transcription factor enhances nitrate uptake in rice by regulating nitrate transporter gene OsNPF6.1, thereby enhancing the regulation of rice nitrogen use efficiency (Tang et al., 2019). In this study, we investigated the effect of co-overexpression of OsNAR2.1 and OsNRT2.3a on NO 3 uptake, yield and nitrogen use efficiency uptake by rice.
Compared with WT, the dry matter and total nitrogen accumulation of NA lines increased significantly at anthesis and maturity ( Figure 2, Table 2), the DMTE, NTE, CPAGY, and CPNGN increased by 28.9%, 9.8%, 29.7%, and 15.7% respectively ( Table 3). During rice filling, 70%-90% of nitrogen was transported from vegetative organs to panicles (Yoneyama et al., 2016). The   Figure 3). Compared with WT, the GY, ANUE, and NRE of N lines increased by 12.1%, 15.8%, and 11.1% respectively ( Table  3, Figure 3), which was significantly lower than that of NA lines. Similarly, co-overexpression of OsNAR2.1 and OsNRT2.3a increased agronomic nitrogen use efficiency of cv. Nipponbare rice. Compared with Nipponbare rice, the grain yields of p35S: OsNAR2.1 and p35S:NAR2.1-p35S:NRT2.3a lines increased by 16.0% and 37.3% ( Figure 5H), and the agronomic nitrogen use efficiency increased by 14.9% and 35.2% ( Figure 5I), respectively. This means the co-overexpression of OsNAR2.1 and OsNRT2.3a transgenic plants have functions in different rice varieties and it can provide a physiological basis for rice breeding. Kronzucker et al. (2000) used 13 N to show that the presence of NO 3 promotes NH 4 + uptake, accumulation, and metabolism in rice. Duan et al. (2006) found that increasing NO 3 uptake promotes dry weight and NO 3 accumulation and assimilation of NH 4 + and NO 3 by 'Nanguang', which is an N-efficient rice cultivar, during the entire growth period. Li et al. (2006) showed that supplying NH 4 + and NO 3 enhances OsAMT1;3, OsAMT1;2, and OsAMT1;1 expression compared with supplying only NH 4 + or NO 3 -, thereby enhancing NH 4 + uptake by rice. The influx rates of 15 NH 4 + and 15 NO 3 in pOsNAR2.1:OsNAR2.1 transgenic lines increased by 21 and 22% in 1.25 mM 15 NH 4 NO 3 and 1.25 mM NH 4 15 NO 3 , respectively (Chen et al., 2017). Further studies showed that when 0.2 mM NO 3 -, 2.5 mM NO 3 or 1.25 mM NH 4 NO 3 were provided as the sole nitrogen source, the dry weight and total nitrogen content of NA lines were higher than that of N lines ( Figure 6). In the 5-min 15 N absorption experiment, we found that there was no significant difference in the influx rate of 15 N between N and NA strains at 0.  Figures  8C, D). OsNRT2.1 and OsNRT2.2 are responsible for NO 3 uptake at the roots, while OsNRT2.3a is responsible for NO 3 transport from roots to shoots Yan et al., 2011;Tang et al., 2012;Chen et al., 2016). Overexpression of OsNAR2.1 can increase the NO 3 and NH 4 NO 3 uptake from roots ( Figures 8A,  B), but it cannot increase the transport ratio from roots to shoots ( Figures 8C, D), this because more OsNRT2.3a needs to be expressed in order to increase the roots-to-shoots transport ratio. In this study, we also found that overexpression of OsNRT2.3a alone did not change rice growth. There was no significant difference in NO 3 uptake, agronomic traits, biomass and nitrogen accumulation between A lines and WT (Tables 1 and 2, Figure 2, Figure S3). The reason may be that the expression of OsNRT2.3a alone did not affect the expression of OsNAR2.1 (Figures 1 and 4, Figure S5).
In conclusion, co-overexpression of OsNAR2.1 and OsNRT2.3a could increase the uptake as well as the transport rate of NO 3 and NH 4 NO 3 from roots to shoots, eventually leading to increasing the total nitrogen content and biomass of rice seedlings under 0.2 mM NO 3 -, 2.5 mM NO 3 -, and 1.25 mM NH 4 NO 3 conditions. Field experiments also showed that co-overexpression of OsNAR2.1 and OsNRT2.3a could increase rice biomass and total nitrogen accumulation, as well as improve rice dry matter transport efficiency and nitrogen transport efficiency, inducing improving rice yield, agronomic nitrogen use efficiency and nitrogen recovery efficiency in rice. Through co-overexpression of OsNAR2.1 and OsNRT2.3a have the same advantages in different rice varieties, it may also function in different crops. This approach provides an effective way to improve grain yield and nitrogen use efficiency in plant.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.