NPF genes excavation and their expression response to vernalization and nitrogen deficiency in allotetraploid rapeseed

doi:10.21203/rs.3.rs-236072/v1

Background

The NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY (NPF) genes, initially characterized as nitrate or peptide transporters in plants, involve in the transport of a large variety of substrates including amino acids, nitrate, auxin (IAA), jasmonates (JAs), abscisic acid (ABA) and gibberellins (GAs) and glucosinolates. The evolution and expression diversification of genes determine their functional differentiation in polyploid species.

Results

Among 169 NPF genes excavated in Brassica napus, 97 B. napus NPF (BnaNPF) genes evolved from B. rapa, and 72 BnaNPF genes from B. olereaca. They unevenly distributed on B. napus genome and exhibited obvious synteny with NPF genes in Arabidopsis thaliana, B. rapa and B. olereaca. BnaNPF genes were identified to show diversified expression patterns in 90 different organs or tissues based on transcriptome profile data. Besides, they exhibited complex expression changes in the development process of leaves, silique wall and seeds, which indicated that the expression of BnaNPF genes maybe respond to altered phytohormone and secondary metabolite content through combining with promoter elements enrichment analysis. Furthermore, many BnaNPF genes were detected to response to vernalization with two different patterns and 20 BnaNPF genes responded to nitrate deficiency.

Conclusion

The evolution of BnaNPF genes and their expression pattern including response to vernalization and nitrogen deficiency were characterized and provide valuable information for further functional characterization in rapeseed.

Plant Molecular Biology and Genetics

Brassica napus

NPF gene

transfer protein

evolution

gene expression

Plant NPF (NRT1/PTR FAMIL) proteins, display sequence homology with proton-coupled oligopeptide transporter (POT) family of peptide transporters, belong to the large peptide transporter (PTR) family [1] and are involved in dietary nitrogen absorption in the form of di- and tripeptides [2, 3]. In plants, NPF members are initially characterized as nitrate or peptide transporters. AtNPF6.3, known as AtNRT1.1, is the first plant member discovered as a nitrate transporter in Arabidopsis [4], and the crystal structure of AtNPF6.3 recently determined showed the similarity with PTRs in the bacterial [5, 6]. Subsequently, some NPF members are identified in different plants and demonstrated that they behave as dipeptide transporters [7, 8]. So, nitrate/peptide transport function is believed to be a specific feature of this family in plants over a period of time. However, in recent years, several studies demonstrated that some NPF members could transport an even wider range of substrates, including nitrite, chloride, auxin (IAA), abscisic acid (ABA), jasmonates (JAs), and gibberellins (GAs) and glucosinolates [9–12]. Additionally, some of NPF members are even able to transport more than one different substrate: nitrate/IAA, nitrate/ABA, nitrate/glucosinolates, peptides/amino acids, or GA/JA.

After identifying the first nitrate transporter NPF6.3/NRT1.1 [4], Tsay et al. (1993) subsequently characterized more than half of the NPF nitrate transporters [8]. So far, more than 21 NPF members have been demonstrated to be able to transport nitrate in Arabidopsis, many of which are able to transport other substrates, such as ABA, GA, JA and glucosinolates [8, 13]. Screening 45 out of 53 AtNPF members using modified yeast two-hybrid system with ABA, GA and JA-Ile specific receptor complexes, Chiba et al. (2015) confirmed that nine NPF members had ABA transport function, and 18 NPF members that were able to transport GA and 13 NPF members able to transport the bioactive JA/JA-Ile [14]. NRT1.1/NPF6.3 was confirmed to be able to transport IAA and 2,4-D [15, 16]. In addition, PTR3/NPF5.2, PTR1/NPF8.1, PTR5/NPF8.2 and PTR2/NTR1/NPF8.3 were identified as dipeptide and tripeptide transporters using complementation of yeast strains deficient for peptide uptake [17, 18].

Based on the phylogenetic relationship, NPF members were divided into eight subfamilies (from NPF1 to NPF8) [1]. In NPF1 subfamily, MtNPF1.7, from Medicago truncatula has been characterized to involve in nodulation and root architecture and behaves as a high-affinity nitrate transporter [19–21]. NPF1.1 and NPF1.2, from Arabidopsis, are confirmed to be important for redistributing xylem-borne nitrate to enhance plant growth [22]. NPF2 subfamily contains nitrate, phytohormone and glucosinolates transporters, and NPF2 members show a wide range of tissue and developmental specificity in Arabidopsis. AtNPF2.3, as a root stele transporter, contributes to nitrate translocation to shoots under salt stress [23], and low-affinity nitrate transporter AtNPF2.12 is responsible for nitrate-dependent early embryo development [24]. Besides, five AtNPF members that are capable of transporting glucosinolates belong to NPF2 subfamily [11]. AtNPF3.1 is the only member of the NPF3 subfamily in Arabidopsis and plays the role of nitrite transport [25]. ABA could be transported by NPF4 subfamily members have been demonstrated [10], and some members of the NPF4 subfamily are able to transport GA, such as NPF4.1 and NPF4.2. NPF5 is the largest NPF subfamily and contains numerous members involved in nitrate, ABA, GA, JA and di-peptides transport [14, 26, 27]. In the NPF6 subfamily, most of the members have been demonstrated to be nitrate transporters, including the first identified NPF member AtNPF6.3 [4, 28–30]. At least two different substrates were able to be transported by NPF7 subfamily members based on three well functionally characterized proteins: AtNPF7.2/NTR1.8 and AtNPF7.3/NTR1.5 for nitrate transport [31, 32], as well as OsNPF7.3 for dipeptide transport [33]. Three Arabidopsis NPF8 subfamily members, AtNPF8.1/PTR1, AtNPF8.2/PTR5 and AtNPF8.3/PTR2 have been proved to be dipeptide transporters, while AtNPF8.1 and AtNPF8.2 are able to transport JA-Ile simultaneously [34, 35]. It is noteworthy that 38 and 17 NPF genes have been characterized at the functional level in Arabidopsis and rice respectively, and more than half of them function on nitrate transport and are distributed on eight NPF subfamilies, and play important roles in nitrate absorption, translocation and utilization [13, 31, 36, 37]. Generally, NPF transporters with a very broad substrate specificity have an important function in plant growth and development, and genome-wide identification has been implemented in poplar [38], Medicago truncatula [36], apple [39] and wheat [13], but a few systematic analyses have been conducted for NPF genes in Brassica species.

B. napus, an allotetraploid evolved from the hybridization between two diploid progenitors (Brassica rapa and Brassica oleracea), is an important oil crop in the world. Compared with Arabidopsis, B. napus have experienced a whole-genome triplication, occurred between 7.9 and 14.6 million years ago [40], and a hybridization event via the natural crossing of B. rapa and B. oleracea (both of which diverged from a common ancestor ~\n4 million years ago) about 7500 years ago [41–43]. Recently, the genomes of ‘Darmor-bzh’ and ‘ZS11’ have been successfully sequenced and assembled [42, 44], and Pan-genomes have been constructed based on next-generation sequencing technologies for B. napus, which facilitate to systematically excavate NPF genes in rapeseed. In this study, 169 NPF genes were excavated in the B. napus, and their evolution and orthologous duplication with NPF genes in Arabidopsis, B. rapa and B. oleracea genome were analyzed. Moreover, the expression profiles of BnaNPF genes in diverse tissues, as well as expression changes at different development stages and under nitrate deficiency were determined. These results will provide useful information for further investigation of the biological function of BnaNPF genes for growth and development in B. napus.

Distribution and synteny analysis of NPF genes in four Brassica species

Based on BLASTP using 53 Arabidopsis NPF protein and phylogenetic analysis (Additional file 1: Figure S1), a total of 169 NPF genes encoding 186 proteins were identified in B. napus genome. To investigate the evolution of BnaNPF genes, the synteny of NPF gene pairs between B. rapa and Arabidopsis genome, B. oleracea and Arabidopsis genome, B. napus and B. rapa genome, B. napus and B. oleracea was performed to further understand the expansion mechanism of NPF genes in B. napus (Fig. 1). The result showed that most of the BnaNPF genes exhibited evolutionary and syntenic relationships with NPF genes in Arabidopsis, B. rapa, and B. olereaca (Additional file 2: Figure S2), suggesting the contribution to the evolution of BnaNPF gene family. Furthermore, Ka, Ks and Ka/Ks of orthologous pairs on BnaNPF and AtNPF genes were calculated to test the evolutionary selection pressure (Additional file 6: Table S2). The majority of orthologous BnaNPF gene pairs had Ka/Ks < 1, which suggested that most of BnaNPF genes have undergone purifying selection to preserve gene function. The mean value of NPF3 (Ka/Ks = 0.10), NPF6 (Ka/Ks = 0.11) and NPF7 (Ka/Ks = 0.13) gene pairs was lower than other subfamilies, showing that these three subfamilies may have suffered robust purifying selective pressure during evolution. However, some of BnaNPF genes had Ka/Ks > 1, including BnaA01NPF2.8, BnaC01NPF2.9, BnaA06NPF2.10, BnaC03NPF2.12 and BnaC01NPF2.25 in NPF2 subfamily, BnaA09NPF4.15 in NPF4 subfamily, BnaA05NPF5.1, BnaC04NPF5.3, BnaC03NPF5.7, BnaA03NPF5.8, BnaA02NPF5.15, BnaA02NPF5.40, BnaC02NPF5.41 and BnaC06NPF5.42 in NPF5 subfamily, and BnaC09NPF8.19 in NPF8 subfamily, suggesting that these BnaNPF genes are subjected positive selection during the evolution from Arabidopsis to rapeseed.

The distribution and synteny of NPF genes were marked on the chromosomes of B. rapa, B. oleracea and B. napus (Fig. 1b). NPF genes are unevenly distributed on every chromosome, and often organized as clusters in the genome of B. rapa, B. oleracea and B. napus. In B. napus genome, the chromosomes A09 and C06 possess the most BnaNPF genes (15 respectively), and A08 possess only 4 BnaNPF genes which were clustered on the chromosome terminal. NPF genes distributed on B. rapa and B. oleracea genome keep good collinearity with NPF genes on A and C sub-genome of B. napus, respectively. B. rapa genome contains 82 NPF genes, and the corresponding A sub-genome of B. napus contains only 76 NPF genes; B. oleracea genome contains 70 NPF genes, and the corresponding C sub-genome of B. napus contains 93 NPF genes, which indicates parts of NPF genes from B. rapa genome were lost or recombined to the C genome of B. napus in evolution process. For example, BraNPF5.21 on the terminal of the chromosome BraA05 was replicated and recombined to BnaC05 chromosome (BnaC05NPF5.37 and BnaC05NPF5.38). According to the synteny analysis, 97 BnaNPF genes evolved from B. rapa genome, and 72 BnaNPF genes from B. oleracea genome. Furthermore, 73 BraNPF genes retained synteny with NPF genes in B. napus genome, including 55 BraNPF genes with 1:1 synteny relationship, and 16 BraNPF genes with a 1:2 relationship (duplication in B. napus genome) and even two BraNPF genes with more 1:2 relationship (1:3 and 1:5) (Table 1 and Additional file 7: Table S3). Nine BraNPF orthologs were not identified in B. napus genome (1:0 relationship) and two BnaNPF orthologs were not identified in B. rapa genome (0:1 relationship), which suggesting loss of the gene during evolutionary. Sixty-one BolNPF genes retained synteny with NPF genes in B. napus genome, including 54 with a 1:1 relationship and 7 with a 1:2 relationship. Twenty-six and five translocations were identified for NPF genes when comparing B. napus genome with B. rapa and B. oleracea genome, respectively. Besides, because the genomic data of B. napus has not yet been fully mapped to the chromosome, the chromosomal location and evolution of three BnaNPF genes (BnaNPF2.26, BnaNPF2.29 and BnaNPF2.30) is still unclear.

Table 1. The synteny relationship of NPF genes between B. rapa and B. napus, and between B. oleracea and B. napus

Ratio^a	0:1	1:0	1:1	1:2	1:3	1:5
*B. rape*	2	9	55	16	1	1
*B.oleracea*	4	9	54	7

^a Orthologous NPF gene ratio by comparing B. rapa and B. oleracea with B. napus genome. 0:1 represents NPF orthologs lost in B. rapa or B. oleracea genome, 1:0 represent lost in B. napus genome, 1:2, 1:3 and 1:5 represent different replication multiples in B. napus genome.

B. napus genome possessed the most NPF genes

Using the sequences of 53 Arabidopsis NPF family protein as queries to perform BLASTp and the information from the article Leran et al (2014) reported [1], NPF proteins from sequenced 36 species were retrieved including B. rapa, B. oleracea and B. napus. The information of genome size and number of NPF genes were shown in Table 2. The genome sizes of these 34 plant species ranged from 127.42 Mb (Arabidopsis) to 2,271.03 Mb (Zea mays), and the NPF gene number was varied from 23 (Physcomitrella patens) to 169 (B. napus). The B. napus genome possessed the most NPF genes (167) through its genome size smaller than that of Malus domestica and Zea mays, which indicated the copy number variations of B. napusNPF genes might be attributed to their requirement for (un)specific substrates as a result of evolutionary selection, such as some NPF2 members for transporting glucosinolates [8]. All NPF genes were grouped into eight clades with known 53 NPF members from Arabidopsis. Most plants have more NPF2, NPF4 and NPF5 subfamily members. NPF1 and NPF2 subfamilies are absent from the two lower plants Physcomitrella patens and Selaginella moellendorffii. In addition, based on BnPIR database that provides more detailed annotation for B. napus genes, 11 BnaNFP genes were identified to encode two proteins derived from two different transcripts, and three BnaNFP genes encode three proteins translated from three different transcripts. Therefore, a total of 186 BnaNPF proteins were identified in B. napus including 17 proteins from different transcripts. Based on the phylogenetic tree (Additional file 1: Figure S1), the evolutionary relationship of NPF proteins between B. napus and Arabidopsis was easy to compare and provided a good guide for studying the function of NPF genes in B. napus. According to known ArabidopsisNPF proteins subfamily information and phylogenetic tree branch, eight unambiguous clades that represented eight B. napusNPF subfamilies were identified. The BnaNPF5 subfamily was the largest because of a larger number of Arabidopsis NPF5 members and possessed 63 members (more than a third of the total number of BnaNPF genes), followed by NPF2 (30), NPF8 (19), NPF4 (16), NPF6 (15), NPF7 (10), NPF1 (10), and NPF3 (6). Additionally, BnaA05NPF5.1 and BnaC04NPF5.2, located in the same branch with AtNPF5.1, were grouped into NPF2 clade, which suggested that the two NPF genes might be more closely related to B. napus NPF2 in evolution. Similarly, BnaA02NPF6.14 and BnaC02NPF6.15 seemed to be more closely related to NPF7. Most of the phylogenetic branches within the same clade showed high bootstrap value (> 0.80), which reflected the low genetic differentiation of Arabidopsis and B. napus NPF genes within the subfamily.

Table 2. Copy number variations (CNVs) of the NPF genes in 36 plant species

Organism Name	*NPF*1	*NPF*2	*NPF*3	*NPF*4	*NPF*5	*NPF*6	*NPF*7	*NPF*8	Total	Genome Size (Mb)
Arabidopsis lyrata (D)	3	14	1	9	17	4	3	5	56	202.97
*Arabidopsis*thaliana (D)	3	14	1	7	16	4	3	5	53	127.42
Aquilaria agallochum (D)	6	7	3	12	13	5	3	6	55	726.71
Brachypodium distachyon (M)	2	6	4	13	21	8	11	17	82	271.3
*Brassica rapa* (D)	4	23	3	9	23	7	5	8	82	401.93
*Brassica oleracea* (D)	4	15	2	8	26	6	5	4	70	554.98
*Brassica napus* (D)	10	30	6	16	63	15	10	19	169	976.19
Carica papaya (D)	4	14	3	8	12	8	6	4	59	370.42
Capsella rubella (D)	3	12	1	6	17	4	3	5	51	133.06
Citrus clementina (D)	9	7	3	9	17	6	4	4	59	301.37
Citrus sinensis (D)	8	7	3	10	17	6	4	4	59	319.23
Cuscuta campestris (D)	4	9	3	8	19	6	5	5	59	476.79
Eucalyptus grandis (D)	6	12	4	11	19	6	4	6	68	691.43
Fragaria vesca (D)	0	13	2	8	23	3	5	6	60	214.37
Glycine max (D)	13	14	6	22	41	11	14	13	134	927.71
Gossypium raimondii (D)	7	10	4	14	14	11	7	8	75	773.77
Linum usitatissimum (D)	12	7	4	14	25	9	11	10	92	316.17
Malus domestica (D)	2	34	4	21	44	17	8	9	139	1,874.77
Manihot esculenta (D)	7	12	6	10	23	7	5	5	75	292.1
Medicago truncatula (D)	8	12	3	14	25	8	9	1	80	412.92
Oryza sativa (M)	3	6	5	12	29	6	11	21	93	389.75
Phaseolus vulgaris (D)	8	11	3	12	22	5	7	6	74	521.08
Populus trichocarpa (D)	15	9	5	12	26	6	5	7	85	434.29
Prunus persica (D)	2	15	1	8	16	5	5	5	57	214.22
Ricinus communis (D)	5	20	3	7	13	5	4	3	60	350.62
Setaria italica (M)	4	11	8	16	19	7	12	21	98	405.87
Solanum tuberosum (D)	17	10	2	15	8	9	4	8	73	772.25
Solanum lycopersicum (D)	19	16	2	12	11	12	7	11	90	760.07
Sorghum bicolor (M)	4	8	7	16	22	6	9	19	91	709.35
Theobroma cacao (D)	4	14	3	10	19	7	4	5	66	345.99
Vitis vinifera (D)	4	7	2	6	21	5	4	3	52	486.2
Zea mays (M)	4	4	6	12	17	8	12	16	79	2,271.03
Amborella trichopoda (D)	1	5	2	7	15	4	3	7	45	706.60
Physcomitrella patens (L)	0	0	1	1	8	6	3	4	23	472.081
Selaginella moellendorffii (L)	0	0	4	4	11	6	5	16	46	212.315
Selaginella moellendorffii (L)	0	0	4	4	11	6	5	16	46	212.315

D dicots, M monocots, L lower plants

BnaNPF gene owning PTR2 functional domain and might be regulated by multiple phytohormones

The gene structures (number and organization exon-intron) are typical evolutionary imprints within certain gene families and are closely related to their function. The exon/intron arrangements of 169 BnaNPF genes were analyzed together with 53 AtNPF by comparing CDS and the corresponding genomic DNA sequences within and between subgroups based on the phylogenetic tree (Additional file 3: Figure S3). The BnaNPF genes have a higher degree of divergence among gene structure than NPF genes in Arabidopsis and contained the numbers of exons varying from 2 to 18. BnaC02NPF1.8 and BnaC09NPF1.9 in NPF1 subfamily, BnaC05NPF2.6, BnaA06NPF2.7 and BnaA06NPF6.8 in the BnaNPF2 subfamily were significantly longer than other genes and contained the most exons (16, 16, 18, 18 and 8, respectively), however, most of the BnaNPF genes contained no more than 6 exons. BnaNPF genes in different branches exhibited different gene structural features, while the genes in the same branch generally had similar intron/exon distribution patterns. For instance, BnaA05NPF1.4, BnaC05NPF1.5 and BnaA06NPF1.7 in BnaNPF1 subgroup, BnaC05NPF5.56, BnaA09NPF5.57, BnaA07NPF5.58 and BnaC07NPF5.59 in BnaNPF5 subgroup, and BnaC07NPF7.3, BnaA03NPF7.4, BnaC01NPF7.5 and BnaA01NPF7.6 in BnaNPF7 subgroup had almost the same exon/intron distribution characteristics within subgroup, and different distribution patterns between subgroups. To further explore the specific and conserved regions of 186 BnaNPF proteins, four conserved domains (PTR2, MFS_1, Chorismate_bind and PDDEXK_6) were identified by the HMMER website (Additional file 3: Figure S3). PTR2 domain, responsible for proton-dependent transport, is the signature domain of NPF protein and could be found in each BnaNPF member, which suggesting functional conservation. Major facilitator superfamily MFS_1 domain feature was detected to partially overlaps or within the PTR2 domain in some BnaNPF members (45/186). The chorismite_bind domain involved in chorismate-utilizing was found in BnaC05NPF2.6 and BnaA06NPF2.7, and BnaC03NPF4.4 contained an unknown function PDDEXK_6 domain.

Transcription factors bind to CREs in the promoter and regulate the expression of the target genes (Wittkopp and Kalay 2011). Generally, genes with similar CREs show the same expression patterns. The 2.0-kb upstream regulatory regions of the BnaNPF genes were used to explore the CREs (Fig. 2 and Additional file 8: Table S4). The result showed that 157 BnaNPF genes contained at least one type of CREs in the promoter regions, which indicated that complex transcriptional regulatory might be implicated for BnaNPF genes. Apart from the common CREs, such as the CAAT-box, TATA-box and some light-responsive elements (G-box, Box 4, GT1-motif and TCT-motif), some phytohormone responsive elements, such as the auxin-responsive elements (TGA-element, AuxRR-core, GATA-box, TGA-box and AuxRE), the ABA-responsive element (ABRE) and the JA-responsive elements (CGTCA-motif and TGACG-motif), and some abiotic stress-responsive elements, such as the low-temperature responsive element (LTR), the salicylic acid responsive element (TCA-element) and the anaerobic responsive element (ARE) were identified. Some over-presented CREs, including ARE, ABRE, CGTCA-motif, TGACG-motif, LTR and TC-rich repeats, were involved in the molecular response of plants to phytohormone, defense and stress responsiveness (Fig. 2a). Among these, the MYB recognition site was most enriched, implying that the MYB transcript factors may play crucial roles in the transcriptional regulation of the BnaNPF genes. Besides, RY-element, the CRE involved in seed-specific regulation, was identified in the promoters of the 15 BnaNPF genes, which indicated that these BnaNPF genes might function in the process of seed development and matter storage.

Gene expression pattern analysis of NPF genes in diverse tissues of B. napus

In order to explore the potential tissues in which NPF genes function in B. napus, the expression profiles were characterized in 90 different organs or tissues, including cotyledon, root, vegetative rosette, stem peel (peel of upper, middle and lower stem), leaf (23 parts or periods), sepal, petal, filament, pollen, bud, silique wall (30 development periods) and seed (24 development periods) based on transcriptome information from BnTIR (http://yanglab.hzau.edu.cn/BnTIR/eFP). Except for half of the genes in the BnaNPF2, BnaNPF5 and BnaNPF8 subfamily that has relatively low expression values (FPKM <1) or no expression, most of the BnaNPF genes had preferential expression profiles in the 90 tissues (Fig. 3). For instance, half of BnaNPF1 genes showed high expression levels in the silique wall at the early and middle development stages and in leaves of all parts; One-third of BnaNPF2 genes (10/30) showed specific expression in the seeds at early and middle development stages; Most of BnaNPF7 genes (8/10) were highly expressed in the bud, petal, pollen and seeds. In general, expression patterns were conserved in each clade within a subfamily, but were quite different across different subfamilies, suggesting the expression differentiation trend of this gene family. For instance, expression patterns of BnaNPF2 and BnaNPF4 subfamilies were classified into three conserved patterns that were consistent with the three major clades in these two subfamilies, and while the expression profile of the BnaNPF3 genes was similar in this subfamily.

Based on the expression profiles in seeds, silique wall and leaves from multiple development periods or plant parts, the expression patterns of BnaNPF genes in leaves, silique wall and seeds could be clarified clearly (Fig. 4). Although some members of both BnaNPF1 (4/10) and BnaNPF2 (5/30) were highly expressed in the silique wall of the developing silique, BnaNPF1 genes showed higher expression levels at the middle development stages, and BnaNPF2 genes were higher expressed at the later development stages. The members of BnaNPF3 with high expression levels in the silique wall, BnaA07NPF3.1 and BnaC06NPF3.2, were higher expressed at the early than later development stage of the silique. However, some members of BnaNPF4 and BnaNPF5, such as BnaA06NPF4.10, BnaC03NPF4.11, BnaC04NPF5.3, BnaA04NPF5.4, BnaA08NPF5.54 and BnaC08NPF5.55 were higher expressed in silique wall at the later than early development stages of the silique. BnaC09NPF5.9 and BnaA09NPF5.10 were found preferential high expression in aged leaves and silique walls and nearly mature seeds. BnaC01NPF7.5 and BnaA01NPF7.6 showed higher expression at later development stages of seed, and BnaC07NPF8.7, BnaA06NPF8.8 and BnaC09NPF8.9 were preferential higher expressed in aged leaves and silique wall.

Expression dynamic of NPF genes during the growth of B. napus under vernalization

There are differences in nutrition utilization and phytohormone distribution at different stages of plant growth. In order to explore the function and expression variation of BnaNPF genes, the expression trend of BnaNPF genes in leaves was analyzed during the growth of B. napus based on ZS11 transcription data from online data resources BnPIR (http://cbi.hzau.edu.cn/bnapus/). Although the expression levels were quite different, members of the same subfamily usually have the same expression trend in leaves during the growth of B. napus (Fig. 5). The members in subfamily NPF1, NPF4, NPF5 and NPF7 seem to be the same expression trend, decline at beginning of vernalization or in the early stage of vernalization and rise after vernalization. For example, BnaC06NPF1.6, as an ortholog of AtNPF1.2 that was able to transport GA and JA, has the exactly same expression trend and high expression level with BnaC04NPF5.3 (homologous with AtNPF5.1 that was able to transport GA, JA and ABA), which indicated that they might play important roles in phytohormone transport for a developmental phase transition. Some other members in NPF2, NPF3 and NPF6 shared this similar expression trend, the expression level raise during vernalization and declined after vernalization. In typical cases, the expression level of BnaC02NPF2.6, BnaC06NPF3.2 and BnaC05NPF6.10 are dramatically raised from T1 to T2, and then begin to decline, which indicated that these members played an important role in the development stage during vernalization. Many BnaNPF genes showed diverse expression levels in the leaves of different cultivars at certain development stages (Additional file 4: Figure S4). For example, BnaA05NPF1.4 and BnaC05NPF1.5 have no expression or lower expression levels in Shengli than other cultivars at T3 and T4 stages. At the T2 stage, BnaC02NPF2.16 showed obviously a higher expression level in cultivars Quinta, Shengli and Tapidor than others. BnaA06NPF8.8 has almost no expression during the whole development process in the three cultivars Shengli, Tapidor and Westar in comparison to other cultivars. These expression variations might lead to differences in nitrogen utilization efficiency, peptide transport and polar transport of phytohormone among the cultivars.

Transcriptional analysis of BnaNPF genes under nitrate deficiency

Nitrate is the main substrate that NPF proteins transport and more than one-third NPF members have been reported to have nitrate transport function in Arabidopsis [8]. Here, we analyzed the expression changes of BnaNPF genes under the condition of nitrogen suitability and deficiency. A total of 20 BnaNPF genes were detected to have relatively high expression and showed significant expression changes in shoot and/or root (Fig. 6). Among them, six BnaNPF genes (BnaC06NPF4.16, BnaC04NPF6.1, BnaA07NPF6.2, BnaA08NPF6.6, BnaC05NPF7.7 and BnaA09NPF7.8) were expressed at a high level in both shoot and root, and the expression levels were significantly elevated in both shoot and root after treated with low nitrogen. Ten BnaNPF genes were specifically expressed in root, of which seven (BnaA06NPF2.7, BnaC06NPF2.20, BnaC08NPF6.4, BnaA09NPF6.5, BnaA06NPF6.8, BnaC05NPF6.9 and BnaA05NPF7.10) were induced to highly express after low nitrogen treatment, which suggested they have a positive function for nitrogen absorption by roots. However, the expression of the other three BnaNPF genes (BnaC06NPF4.8, BnaC09NPF4.14 and BnaC07NPF7.3) that specifically expressed in roots were declined under low nitrogen treatment. In addition, four BnaNPF genes that were specifically expressed in shoots also showed different expression changes under low nitrogen treatment: BnaC02NPF2.16 and BnaA06NPF4.10 were upregulated, and the other two (BnaC06NPF3.2 and BnaC07NPF6.3) were declined after treated by low nitrogen.

Although its genome is not the largest in comparison to the genomes of 33 plant species displayed in our study, B. napus contained the most NPF genes (Table 1). A total of 169 BnaNPF genes coding 186 proteins were identified in the B. napus genome in this study and designated as BnaNPF1.1 to BnaNPF8.19 in eight subfamilies based on phylogenetic analysis, and they exhibited evolutionary and syntenic relationships with NPF genes in Arabidopsis, B. rapa, and B. olereaca. Furthermore, the expression profiles of BnaNPF genes in 90 diverse tissues, as well as expression changes at different development stages under vernalization between 8 rapeseed cultivars and under nitrate deficiency were determined. This study provides a piece of basic information for further functional characterization of BnaNPF genes in the growth and development of B. napus. Recently, in apple and wheat [13, 39], the NPF protein family was characterized and also defined as eight subfamilies, and according to eight subfamilies of NPF defined in these species and phylogenetic analysis, the 169 BnaNPF genes were identified and classified into eight unambiguous subfamilies, and BnaNPF subfamilies showed similar member proportions with these in Arabidopsis and wheat [1, 13].

Brassicaceae species experienced a common whole genome soon after divergence from the Arabidopsis lineage approximately 17 to 20 million years ago [40, 45]. B. napus is an allotetraploid (AnAnCnCn) that evolved from a spontaneous hybridization event between B. rape (ArAr) and B. oleracea (CoCo) about 7500 years ago [42], and have suffered a whole-genome triplication and a hybridization event compared with Arabidopsis. In theory, there should be three-time as much NPF genes in B. rapa and B. oleracea (53 × 3), and six times as much NPF genes in B. napus (53 × 3 × 2) as in Arabidopsis. Fifty-three NPF members were identified in the Arabidopsis genome, it was expected that NPF genes may be expanded to about 160 genes in B. rapa or B. oleracea, and about 320 genes in B. napus genomes, respectively. However, only 82, 70 and 169 BnaNPF genes were identified in these three species respectively in this study (Table 1), which revealed genome replication was accompanied by large scale loss of genes during evolution that identical to the previous reports [46, 47]. The A sub-genome and C sub-genome of B. napus (AACC) were originated from B. rapa (AA) and B. oleracea (CC), respectively. Comparison to both ancestral species, fewer NPF genes (76) were identified in the A sub-genome and more NPF genes (93) were discovered in the C sub-genome of B. napus. NPF genes in the C sub-genome were amplified obviously, which happened probably due to chromosome rearrangement or gene replication when B. napus formed by hybridization between B. rapa and B. oleracea about 7500 years ago [42]. Besides, NPF genes distributed on B. rapa (88.16%) and B. oleracea (55.91%) genome keep good collinearity with NPF genes on A and C sub-genome of B. napus, respectively (Fig. 1). These results indicated that BnaNPF genes have undergone not only chromosome segment replication, but complex recombination and gene loss in evolution processes.

The function and expression level of a gene is usually closely related to its gene characteristic and CREs [48]. Therefore, BnaNPF genes were furtherly characterized for gene structures, protein conserved domains and CREs in this study. Most of the BnaNPF genes exhibited relatively concentrated distributional property in gene length (246-1568bp) and amino acid number (400-600). In gene structure, most of the BnaNPF genes contained at most six extron, and different BnaNPF subfamily genes exhibited significant exon-intron structural divergences, but BnaNPF within the same branches share similar gene structures, motifs, and localization patterns. Besides, CREs involved in hormone responses and MYB recognition site were detected in the promoter region of BnaNPF genes except for the common CREs, which indicated the expression of BnaNPF genes regulated by phytohormones and secondary metabolites (Fig. 2).

Gene expression patterns provided imperative clues to map out gene functionality. In this study, the expression level of BnaNPF genes was investigated in diverse tissues or organs of B. napus using the released transcriptome information resource (http://yanglab.hzau.edu.cn/BnTIR/eFP). Gene expression analysis showed BnaNPF genes have significantly different and complex expression patterns across different subfamilies, but the expression pattern was conserved in each clade within a subfamily, which reflected structure and function uniform (Fig. 3). Some BnaNPF genes showed obvious tissue preferential expression. Half of BnaNPF1 genes having ultrahigh expression levels in silique wall at the early and middle development stages of silique indicated efficient nitrogen transport for nutrient synthesis in the seed. BnaA07NPF2.18, orthologous to AtNPF2.13/NRT1.7 that was able to transport nitrate, glucosinolates, JAs and GAs [14, 49], was expressed at ultrahigh expression level in the bud, pollen, filament and petal, which contribute to the previous reports that nitrate and nitrogen regulated flowering, and high nitrate/nitrogen helped promote flowering [50], and phytohormone played an important role in the regulation of flower organogenesis [51, 52].

NPF proteins can transport a huge variety of substrates, including dipeptides, nitrate, glucosinolates, amino-acids and several plant hormones [8], the complex gene expression pattern would endow BnaNPF versatile roles in the growth and development of B. napus. Many BnaNPF genes were found to have a changing expression in the development of organs or tissues (Fig. 4). Such as, in the BnaNPF2 subfamily, BnaA06NPF2.10, BnaA06NPF2.11 and BnaC03NPF2.12, were orthologs to AtNPF2.10; and BnaC02NPF2.16, BnaC06NPF2.17 and BnaA07NPF2.18, were orthologs to AtNPF2.13, which two Arabidopsis orthologs had the function of transporting glucosinolates [11, 53], showed up-regulated expression in later stages of silique and seed development of B. napus. Many CREs involved in hormone responses were detected in the promoter region of BnaNPF genes, including IAA- (103/169 genes), ABA- (123/169 genes), and MeJA-responsive CRE (119/169 genes) (Fig. 2 and Additional file 8: Table S4), which suggested their potential hormone-inducing characteristics. The process of plant growth and development was regulated by phytohormone directly, which might be why the transcription of BnaNPF genes was regulated responding to growth and development. Besides, with the development of plant organs, secondary metabolite accumulation level maybe also played a part in the expression changes of some BnaNPF genes because of the existence of MYB recognition site in their promoters [54].

The expression changes of BnaNPF genes during the growth of B. napus under natural vernalization were analyzed in this study. Vernalization is an important process that regulates the transition from vegetative growth to reproduction in B. napus [55, 56], and involved in the regulation of various environmental factors and hormones [57]. The BnaNPF genes that expressed in leaves exhibited two expression trends: the first one, decline at beginning of vernalization or in the early stage of vernalization and rise after vernalization (the most of members of NPF1, NPF4, NPF5 and NPF7 subfamily); the second one, the expression level was raised during vernalization and declined after vernalization (the most of members of NPF2, NPF3 and NPF6 subfamily) (Fig. 5). These results indicated that many BnaNPF genes might participate in floral transition and play different roles in the reproduction and development of B. napus. Based on the transcriptome data of eight cultivars from BnPIR database, significant expression variation was found for some BnaNPF genes in different cultivars (Additional file 4: Figure S4). These expression variations might lead to differences in transport of the corresponding substrates among the cultivars, which is expected for further functional research in the future.

Nitrate and phytohormone signaling pathways crosstalk to modulate growth and developmental programs in a multifactorial manner [58]. So far, more than half of functionally characterized NPF genes have been demonstrated to be able to transport nitrate in Arabidopsis [13]. Here, twenty BnaNPF genes, ortholog to 11 AtNPF genes, were detected to respond to low nitrate treatment (Fig. 6). The six members of the BnaNPF6 subfamily, BnaNPF6.4-6.9 orthologous with AtNPF6.3/NRT1.1, were predominantly expressed in roots and were significantly up-regulated under low nitrogen treatment, suggesting their functional importance in nitrogen utilization efficiency. AtNPF6.3 is the first plant NPF member that characterized for functioning in nitrate uptake in the root, root-to-shoot transport and transceptor in sensing/signaling, and govern many molecular, physiological, and morphological responses to nitrate [6, 15, 16]. The gene expansion and consistent expression patterns in B. napus indicated the function uniformity of NPF6.3 orthologs in nitrogen utilization efficiency as previously reported [59, 60]. As the ortholog of AtNPF2.9/NRT1.9 has been reported to facilitate the loading of nitrate into the root phloem and enhance downward nitrate transport in roots [61], BnaA06NPF2.7 were up-regulated significantly in root under low nitrate condition. While its homologs in rice, OsNPF2.4, was discovered by a genome-wide association study (GWAS) on nitrogen utilization efficiency-related agronomic traits [62]. So BnaA06NPF2.7 might also play an important role in nitrate transport in roots in B. napus, which needs to be characterized in the future. In addition to the role as a nutrient, nitrate acts as a signal molecular, and N nutrition and plant hormone signaling pathways are closely interconnected [58]. BnaC02NPF2.16 and BnaC06NPF3.2, orthologues with AtNPF2.13 and AtNPF3.1 respectively, were predominantly expressed in leaves. According to previous reports about AtNPF2.13 and AtNPF3.1 that remobilizing nitrate from old to young leaves and involved in GA accumulation and responses [63-65]. BnaC02NPF2.16 and BnaC06NPF3.2 might function in nitrite accumulation in leaves coupling hormone signal, which may be possible but needs to be verified in the future.

In this study, we provided a comprehensive knowledge of the evolution and expression characteristics of BnaNPF genes in B. napus. It gives an important implication for further understanding the biological functions of individual BnaNPF genes. However, the study only provided a preliminary characterization of BnaNPF genes and large functional validation work needs to be done in further work to understanding the roles of BnaNPF genes.

A total of 169 NPF gene members were identified in the B. napus genome and classified into eight subfamilies in this study. The BnaNPF genes was unevenly distributed in B. napus genome and exhibited obvious synteny and orthologous duplication with NPF genes in Arabidopsis, B. rapa and B. olereaca. Moreover, the complex expression patterns of NPF genes in various tissues and periods were investigated, and the expression changes at different development stages under nature vernalization and response to nitrate deficiency were determined in B. napus. The evolution and expression pattern analysis of NPF genes will provide valuable information for further functional characterization in rapeseed.、

Data resource related to NPF genes acquisition

The 53 NPF protein sequences from Arabidopsis were used as query, BLASTp search (E-value < e-10) was performed to identify NPF genes in B. rapa, B. oleracea, B. napus and other plant species through the online BLAST tool in the databases, included The Arabidopsis Information Resource (TAIR, https://www.Arabidopsis.org/, Genome Version Araport11), the Brassica Database (BRAD, http://brassicadb.org/brad/, var. ‘Chiifu’ and ‘TO1000’), Brassica napus Pan-Genome Information Resource (BnPIR, http://cbi.hzau.edu.cn/bnapus/, var. ZS11), National Center for Biotechnology Information (NCBI, https://blast.ncbi.nlm.nih.gov/Blast.cgi). Then, the potential NPF proteins were confirmed through the hidden Markov model (HMM) search program (HMMER v3.0, http://hmmer.janelia.org/) with E-value below e-200, and the conserved domain database (CDD) (http://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi) based on the presence of the PTR2 domains (PF00854). Furthermore, the candidate NPF protein of B. rapa, B. oleracea and B. napus with less than 200 amino acid residues and 20% of PTR2 domains missing were removed, and the rests were thought to be considered to be functional and used for further analysis.

Multiple sequence alignment and BnaNPF genes nomenclature

Full-length sequences of the NPF proteins from Arabidopsis and three Brassica crops were aligned with ClustalW, and then these alignments were used to construct the phylogenetic trees by the software MEGA Version 10.1.0 [66] with the neighbor-joining method. P-distance, pairwise deletion, and bootstrapping (1000 replicates) were set as the required parameters.

NPF genes were named according to the nomenclature Leran et al. (2014) recommended [1]. According to eight subfamilies of NPF in Arabidopsis and phylogenetic relationships, the clade number of NPF members would be ensured. Then NPF members were named with two or three letters to identify the species, followed by “NPF + clade number (followed by a point) + the order number (which they are identified in phylogenetic tree)”, for instance, “BraNPF2.3”. Consequently, this second number is used to differentiate genes within the species and does not reflect orthologous relationships. The NPF members from B. napus obeyed the nomenclature convention but modified with adding chromosome between species name and “NPF”. If multiple NPF proteins were translated from the transcripts of the same gene, they were distinguished by the English letters “a”, “b” and “c”.

Chromosomal location and syntenic analysis

The genomic locations of all BnaNPF genes were mapped to chromosomes of the B. napus genome according to the reference genome information of ZS11 in BnPIR database. The synteny orthologous gene pairs were identified based on BLASTP (identity > 75%, and E-value < e-20) and phylogenetic relationship. The chromosomal regions within 200 kb containing a string of two or more genes were defined as tandem duplication [67]. The nonsynonymous rate (Ka), synonymous rate (Ks), and Ka/Ks between the orthologous gene pairs were calculated using the NY method implemented in the Ka/Ks calculator program [68] according to gene CDS pairwise alignment performed with Clustal W (https://www.genome.jp/tools-bin/clustalw).

Functional domain validation and cis-acting regulatory elements (CREs) prediction

The protein sequence and full-length coding sequences (CDS) information of the AtNPFs and BnaNPFs were retrieved and extracted from the Arabidopsis Information Resource (TAIR: https://www.Arabidopsis.org/index.jsp) and Brassica napus pan-genome information resource (BnPIR: http://cbi.hzau.edu.cn/bnapus/index.php). To examine the structural divergence among the NPF proteins in Arabidopsis and B. napus, the protein sequences were subjected to the HMMER software (http://www.ebi.ac.uk/Tools/hmmer) to predict and characterize the conserved domains with default parameters. A 2.0-kb genomic sequence upstream from the start codon was downloaded for each gene from the BnPIR website (http://cbi.hzau.edu.cn/bnapus/index.php). These sequences were subjected to plantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) to identify putative CREs and CREs distribution in the promoter region were displayed by TBtools software [69]. The Gene Structure Display Server (GSDS Version 2.0, http://gsds.cbi.pku.edu.cn/index.php) was used to display the exon-intron structures of the NPFs in Arabidopsis and B. napus.

Identification of expression pattern of BnaNPF genes in B. napus

The fragments per kilobase of exon model per million mapped fragments (FPKM) value of 169 BnaNPF genes in different organs or tissues, including cotyledon, root, vegetative rosette, stem peel (peel of upper, middle and lower stem), leaves (23 parts or periods), sepal, petal, filament, pollen, bud, silique wall (30 development periods) and seeds (24 development periods) based on transcriptome information were retrieved from BnTIR (http://yanglab.hzau.edu.cn/BnTIR/eFP). RNA-seq for leaves of all eight rapeseed accessions at five stages with one-month interval derived from BnPIR (http://cbi.hzau.edu.cn/bnapus/expression/), including one stage before vernalization (T0), three stages during vernalization (T1, T2 and T3), and one stage post vernalization (T4), were used to analyze BnaNPF gene expression patterns at different development stages.

Expression analysis of BnaNPF genes under low nitrate stress

To further investigate the transcriptional responses of BnaNPF genes under low nitrate stress, the uniform B. napus seedlings (var. ZS11) were hydroponically cultured in Hoagland nutrient solution for 10 days at 7 days after seed germination, and then parts of them were transferred to Hoagland nutrient solution modified with low nitrate (0.3 mM NO^3-) for 3 days. The rapeseed seedlings were cultivated in the culture room as Cui et al. (2020) described [70]. The shoots and roots under low nitrogen treatment for 72 h and control were individually harvested and immediately stored at -80 ℃ for RNA isolation, and each sample contained 3 independent biological replicates. Total RNA was isolated from the frozen samples using a RNAprep Pure Plant Kit (Tiangen), and first-strand cDNA was synthesized from the total RNA using a PrimeScriptTM RT Master Mix Kit (TaKaRa). The cDNA was subjected to quantitative analysis using SYBR^® Premix Ex Taq^TM (Takara Bio) on the Applied Biosystems StepOne^TM Plus Real-time PCR System (Thermo Fisher Scientific, Waltham, MA, USA), as previously described [70]. The BnaNPF primer sequences were obtained from the qPCR Primer Database [71] and are listed in Additional file 5: Table S1. The housekeeping gene BnaACTIN7 was used as a reference gene for normalization and to analyze the BnaNPF gene expression levels via the 2^-ΔΔCt method. Three independent technical replicates were performed for each sample.

Acknowledgments

Not applicable

Authors’ contributions

HC participated in the experiments designing, data analysis, and drafting the manuscript. JH participated in gene and protein sequence retrieval and partial data analysis. WZ and HF did plant culturing and gene expression analysis. YH collects and provides some transcriptome data and gives some proposals. ML and JH took part in designing experiments and supervising the study. All authors have read and approved the final manuscript.

Funding

The work was supported by the Key Research Plan Project of Shaanxi province (2020ZDLNY04-01), Chinese Postdoctoral Science Foundation (2019M662539), Key Scientific Research Projects of Colleges and Universities in Henan Province (K20A180031) and Postdoctoral Research Grant of Henan Province (201901004).

Availability of data and materials

The datasets generated and analyzed during the present study are available from the corresponding author on reasonable request.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that the work was performed without any commercial or financial relationship that could be understood as a potential conflict of interest.

NPF: NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY; IAA: nitrate, auxin; JA: jasmonate; ABA: abscisic acid; GA: gibberellin; BnaNPF: Brassica napus NPF ; Brassica rapa: B. rapa; B. olereaca: Brassica oleracea; Arabidopsis: Arabidopsis thaliana; CRE: cis-acting regulatory element; CDS: coding sequences; Ka: nonsynonymous rate; Ks: synonymous rate; FPKM: fragments per kilobase of exon model per million mapped fragments; T0: 24 days post sowing and before vernalization; T1: 54 days post sowing and during vernalization; T2: 82 days post sowing and during vernalization; T3: 115 days post sowing and during vernalization; T4: 147 days post sowing and post vernalization; DAF: days after flowering.

Leran S, Varala K, Boyer JC, Chiurazzi M, Crawford N, Daniel-Vedele F, David L, Dickstein R, Fernandez E, Forde B et al. A unified nomenclature of NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER family members in plants. Trends Plant Sci. 2014;19(1):5-9.
Daniel H, Spanier B, Kottra G, Weitz D. From bacteria to man: archaic proton-dependent peptide transporters at work. Physiology (Bethesda). 2006;21:93-102.
Meredith D. Review. The mammalian proton-coupled peptide cotransporter PepT1: sitting on the transporter-channel fence? Philos Trans R Soc Lond B Biol Sci. 2009;364(1514):203-7.
Tsay YF, Schroeder JI, Feldmann KA, Crawford NM. The herbicide sensitivity gene CHL1 of Arabidopsis encodes a nitrate-inducible nitrate transporter. Cell. 1993;72(5):705-13.
Tsay YF. PLANT SCIENCE How to switch affinity. Nature. 2014;507(7490):44-45.
Sun J, Zheng N. Molecular Mechanism Underlying the Plant NRT1.1 Dual-Affinity Nitrate Transporter. Front Physiol. 2015;6:386.
Tegeder M, Rentsch D. Uptake and Partitioning of Amino Acids and Peptides. Mol Plant. 2010;3(6):997-1011.
Corratge-Faillie C, Lacombe B. Substrate (un)specificity of Arabidopsis NRT1/PTR FAMILY (NPF) proteins. J Exp Bot. 2017;68(12):3107-13.
Krouk G, Lacombe B, Bielach A, Perrine-Walker F, Malinska K, Mounier E, Hoyerova K, Tillard P, Leon S, Ljung K et al. Nitrate-regulated auxin transport by NRT1.1 defines a mechanism for nutrient sensing in plants. Dev Cell. 2010;18(6):927-37.
Kanno Y, Hanada A, Chiba Y, Ichikawa T, Nakazawa M, Matsui M, Koshiba T, Kamiya Y, Seo M. Identification of an abscisic acid transporter by functional screening using the receptor complex as a sensor. Proc Natl Acad Sci U S A. 2012;109(24):9653-58.
Nour-Eldin HH, Andersen TG, Burow M, Madsen SR, Jorgensen ME, Olsen CE, Dreyer I, Hedrich R, Geiger D, Halkier BA. NRT/PTR transporters are essential for translocation of glucosinolate defence compounds to seeds. Nature. 2012;488(7412):531-34.
Boursiac Y, Leran S, Corratge-Faillie C, Gojon A, Krouk G, Lacombe B. ABA transport and transporters. Trends Plant Sci. 2013;18(6):325-33.
Wang H, Wan Y, Buchner P, King R, Ma H, Hawkesford MJ. Phylogeny and gene expression of the complete NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY in Triticum aestivum. J Exp Bot. 2020;71(15):4531-46.
Chiba Y, Shimizu T, Miyakawa S, Kanno Y, Koshiba T, Kamiya Y, Seo M. Identification of Arabidopsis thaliana NRT1/PTR FAMILY (NPF) proteins capable of transporting plant hormones. J Plant Res. 2015;128(4):679-86.
Krouk G, Crawford NM, Coruzzi GM, Tsay YF. Nitrate signaling: adaptation to fluctuating environments. Curr Opin Plant Biol. 2010;13(3):266-73.
Bouguyon E, Brun F, Meynard D, Kubes M, Pervent M, Leran S, Lacombe B, Krouk G, Guiderdoni E, Zazimalova E et al. Multiple mechanisms of nitrate sensing by Arabidopsis nitrate transceptor NRT1.1. Nat Plants. 2015;1:15015.
Chiang CS, Stacey G, Tsay YF. Mechanisms and functional properties of two peptide transporters, AtPTR2 and fPTR2. J Biol Chem. 2004;279(29):30150-7.
Hammes UZ, Meier S, Dietrich D, Ward JM, Rentsch D. Functional properties of the Arabidopsis peptide transporters AtPTR1 and AtPTR5. J Biol Chem. 2010;285(51):39710-7.
Yendrek CR, Lee YC, Morris V, Liang Y, Pislariu CI, Burkart G, Meckfessel MH, Salehin M, Kessler H, Wessler H et al. A putative transporter is essential for integrating nutrient and hormone signaling with lateral root growth and nodule development in Medicago truncatula. Plant J. 2010;62(1):100-12.
Bagchi R, Salehin M, Adeyemo OS, Salazar C, Shulaev V, Sherrier DJ, Dickstein R. Functional assessment of the Medicago truncatula NIP/LATD protein demonstrates that it is a high-affinity nitrate transporter. Plant Physiol. 2012;160(2):906-16.
Salehin M, Huang YS, Bagchi R, Sherrier DJ, Dickstein R. Allelic differences in Medicago truncatula NIP/LATD mutants correlate with their encoded proteins' transport activities in planta. Plant Signal Behav. 2013;8(2):e22813.
Hsu PK, Tsay YF. Two phloem nitrate transporters, NRT1.11 and NRT1.12, are important for redistributing xylem-borne nitrate to enhance plant growth. Plant Physiol. 2013;163(2):844-56.
Taochy C, Gaillard I, Ipotesi E, Oomen R, Leonhardt N, Zimmermann S, Peltier JB, Szponarski W, Simonneau T, Sentenac H et al. The Arabidopsis root stele transporter NPF2.3 contributes to nitrate translocation to shoots under salt stress. Plant J. 2015;83(3):466-79.
Almagro A, Lin SH, Tsay YF. Characterization of the Arabidopsis nitrate transporter NRT1.6 reveals a role of nitrate in early embryo development. Plant Cell. 2008;20(12):3289-99.
Sugiura M, Georgescu MN, Takahashi M. A nitrite transporter associated with nitrite uptake by higher plant chloroplasts. Plant Cell Physiol. 2007;48(7):1022-35.
Karim S, Holmstrom KO, Mandal A, Dahl P, Hohmann S, Brader G, Palva ET, Pirhonen M. AtPTR3, a wound-induced peptide transporter needed for defence against virulent bacterial pathogens in Arabidopsis. Planta. 2007;225(6):1431-45.
Leran S, Edel KH, Pervent M, Hashimoto K, Corratge-Faillie C, Offenborn JN, Tillard P, Gojon A, Kudla J, Lacombe B. Nitrate sensing and uptake in Arabidopsis are enhanced by ABI2, a phosphatase inactivated by the stress hormone abscisic acid. Sci Signal. 2015;8(375):ra43.
Okamoto M, Vidmar JJ, Glass AD. Regulation of NRT1 and NRT2 gene families of Arabidopsis thaliana: responses to nitrate provision. Plant Cell Physiol. 2003;44(3):304-17.
Chiu CC, Lin CS, Hsia AP, Su RC, Lin HL, Tsay YF. Mutation of a nitrate transporter, AtNRT1:4, results in a reduced petiole nitrate content and altered leaf development. Plant Cell Physiol. 2004;45(9):1139-48.
Morere-Le Paven MC, Viau L, Hamon A, Vandecasteele C, Pellizzaro A, Bourdin C, Laffont C, Lapied B, Lepetit M, Frugier F et al. Characterization of a dual-affinity nitrate transporter MtNRT1.3 in the model legume Medicago truncatula. J Exp Bot. 2011;62(15):5595-605.
Li JY, Fu YL, Pike SM, Bao J, Tian W, Zhang Y, Chen CZ, Zhang Y, Li HM, Huang J et al. The Arabidopsis nitrate transporter NRT1.8 functions in nitrate removal from the xylem sap and mediates cadmium tolerance. Plant Cell. 2010;22(5):1633-46.
Chen CZ, Lv XF, Li JY, Yi HY, Gong JM. Arabidopsis NRT1.5 is another essential component in the regulation of nitrate reallocation and stress tolerance. Plant Physiol. 2012;159(4):1582-90.
Ouyang J, Cai ZY, Xia KF, Wang YQ, Duan J, Zhang MY. Identification and analysis of eight peptide transporter homologs in rice. Plant Sci. 2010;179(4):374-82.
Komarova NY, Thor K, Gubler A, Meier S, Dietrich D, Weichert A, Suter Grotemeyer M, Tegeder M, Rentsch D. AtPTR1 and AtPTR5 transport dipeptides in planta. Plant Physiol. 2008;148(2):856-69.
Choi MG, Kim EJ, Song JY, Choi SB, Cho SW, Park CS, Kang CS, Park YI. Peptide transporter2 (PTR2) enhances water uptake during early seed germination in Arabidopsis thaliana. Plant Mol Biol. 2020;102(6):615-24.
Pellizzaro A, Alibert B, Planchet E, Limami AM, Morere-Le Paven MC. Nitrate transporters: an overview in legumes. Planta. 2017;246(4):585-95.
Babst BA, Gao F, Acosta-Gamboa LM, Karve A, Schueller MJ, Lorence A. Three NPF genes in Arabidopsis are necessary for normal nitrogen cycling under low nitrogen stress. Plant Physiol Biochem. 2019;143:1-10.
Bai H, Euring D, Volmer K, Janz D, Polle A. The nitrate transporter (NRT) gene family in poplar. PLoS One. 2013;8(8):e72126.
Wang Q, Liu C, Dong Q, Huang D, Li C, Li P, Ma F. Genome-Wide Identification and Analysis of Apple NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER Family (NPF) Genes reveals MdNPF6.5 Confers high capacity for nitrogen uptake under low-nitrogen conditions. Int J Mol Sci. 2018;19(9).
Lysak MA, Koch MA, Pecinka A, Schubert I. Chromosome triplication found across the tribe Brassiceae. Genome Res. 2005;15(4):516-25.
Wang X, Wang H, Wang J, Sun R, Wu J, Liu S, Bai Y, Mun JH, Bancroft I, Cheng F et al. The genome of the mesopolyploid crop species Brassica rapa. Nat Genet. 2011;43(10):1035-9.
Chalhoub B, Denoeud F, Liu S, Parkin IA, Tang H, Wang X, Chiquet J, Belcram H, Tong C, Samans B et al. Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science. 2014;345(6199):950-3.
Liu S, Liu Y, Yang X, Tong C, Edwards D, Parkin IA, Zhao M, Ma J, Yu J, Huang S et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat Commun. 2014;5:3930.
Sun F, Fan G, Hu Q, Zhou Y, Guan M, Tong C, Li J, Du D, Qi C, Jiang L et al. The high-quality genome of Brassica napus cultivar 'ZS11' reveals the introgression history in semi-winter morphotype. Plant J. 2017;92(3):452-68.
Town CD, Cheung F, Maiti R, Crabtree J, Haas BJ, Wortman JR, Hine EE, Althoff R, Arbogast TS, Tallon LJ et al. Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy. Plant Cell. 2006;18(6):1348-59.
De Bodt S, Maere S, Van de Peer Y. Genome duplication and the origin of angiosperms. Trends Ecol Evol. 2005;20(11):591-7.
Mun JH, Kwon SJ, Yang TJ, Seol YJ, Jin M, Kim JA, Lim MH, Kim JS, Baek S, Choi BS et al. Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication. Genome Biol. 2009;10(10):R111.
Jeffares DC, Penkett CJ, Bahler J. Rapidly regulated genes are intron poor. Trends Genet. 2008;24(8):375-8.
Fan SC, Lin CS, Hsu PK, Lin SH, Tsay YF. The Arabidopsis nitrate transporter NRT1.7, expressed in phloem, is responsible for source-to-sink remobilization of nitrate. Plant Cell. 2009;21(9):2750-61.
Lin YL, Tsay YF. Influence of differing nitrate and nitrogen availability on flowering control in Arabidopsis. J Exp Bot. 2017;68(10):2603-09.
Li Q, Liu B. Genetic regulation of maize flower development and sex determination. Planta. 2017;245(1):1-14.
Fredes I, Moreno S, Diaz FP, Gutierrez RA. Nitrate signaling and the control of Arabidopsis growth and development. Curr Opin Plant Biol. 2019;47:112-18.
Andersen TG, Nour-Eldin HH, Fuller VL, Olsen CE, Burow M, Halkier BA. Integration of biosynthesis and long-distance transport establish organ-specific glucosinolate profiles in vegetative Arabidopsis. Plant Cell. 2013;25(8):3133-45.
Zhou M, Memelink J. Jasmonate-responsive transcription factors regulating plant secondary metabolism. Biotechnol Adv. 2016;34(4):441-49.
Li B, Zhao W, Li D, Chao H, Zhao X, Ta N, Li Y, Guan Z, Guo L, Zhang L et al. Genetic dissection of the mechanism of flowering time based on an environmentally stable and specific QTL in Brassica napus. Plant Sci. 2018;277:296-310.
Li H, Fan Y, Yu J, Chai L, Zhang J, Jiang J, Cui C, Zheng B, Jiang L, Lu K. Genome-wide identification of flowering-time genes in brassica species and reveals a correlation between selective pressure and expression patterns of vernalization-pathway genes in Brassica napus. Int J Mol Sci. 2018;19(11):3632.
Bernier G, Perilleux C. A physiological overview of the genetics of flowering time control. Plant Biotechnol J. 2005;3(1):3-16.
Vega A, O'Brien JA, Gutierrez RA. Nitrate and hormonal signaling crosstalk for plant growth and development. Curr Opin Plant Biol. 2019;52:155-63.
Wen ZY, Tyerman SD, Dechorgnat J, Ovchinnikova E, Dhugga KS, Kaiser BN. Maize NPF6 Proteins are homologs of Arabidopsis CHL1 that are selective for both nitrate and chloride. Plant Cell. 2017;29(10):2581-96.
Wang W, Hu B, Yuan D, Liu Y, Che R, Hu Y, Ou S, Liu Y, Zhang Z, Wang H et al. Expression of the nitrate transporter gene OsNRT1.1A/OsNPF6.3 confers high yield and early maturation in rice. Plant Cell. 2018;30(3):638-51.
Wang YY, Tsay YF. Arabidopsis Nitrate Transporter NRT1.9 Is Important in Phloem Nitrate Transport. Plant Cell. 2011;23(5):1945-57.
Tang WJ, Ye J, Yao XM, Zhao PZ, Xuan W, Tian YL, Zhang YY, Xu S, An HZ, Chen GM et al. Genome-wide associated study identifies NAC42-activated nitrate transporter conferring high nitrogen use efficiency in rice. Nat Commun. 2019;10(1):5279.
Pike S, Gao F, Kim MJ, Kim SH, Schachtman DP, Gassmann W. Members of the NPF3 transporter subfamily encode pathogen-inducible nitrate/nitrite transporters in grapevine and Arabidopsis. Plant Cell Physiol. 2014;55(1):162-70.
David LC, Berquin P, Kanno Y, Seo M, Daniel-Vedele F, Ferrario-Mery S. N availability modulates the role of NPF3.1, a gibberellin transporter, in GA-mediated phenotypes in Arabidopsis. Planta. 2016;244(6):1315-28.
Tal I, Zhang Y, Jorgensen ME, Pisanty O, Barbosa ICR, Zourelidou M, Regnault T, Crocoll C, Olsen CE, Weinstain R et al. The Arabidopsis NPF3 protein is a GA transporter. Nat Commun. 2016;7.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725-29.
Holub EB. The arms race is ancient history in Arabidopsis, the wildflower. Nat Rev Genet. 2001;2(7):516-27.
Zhang Z, Li J, Zhao XQ, Wang J, Wong GK, Yu J. KaKs_Calculator: calculating Ka and Ks through model selection and model averaging. Genomics Proteomics Bioinformatics. 2006;4(4):259-63.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194-202.
Cui JQ, Hua YP, Zhou T, Liu Y, Huang JY, Yue CP. Global landscapes of the Na(+)/H(+) antiporter (NHX) family members uncover their potential roles in regulating the rapeseed resistance to salt stress. Int J Mol Sci. 2020;21(10).
Lu K, Li T, He J, Chang W, Zhang R, Liu M, Yu M, Fan Y, Ma J, Sun W et al. qPrimerDB: a thermodynamics-based gene-specific qPCR primer database for 147 organisms. Nucleic Acids Res. 2018;46(1):1229-36.

Tables

Table 1. The synteny relationship of NPF genes between B. rapa and B. napus, and between B. oleracea and B. napus

Ratio^a	0:1	1:0	1:1	1:2	1:3	1:5
*B. rape*	2	9	55	16	1	1
*B.oleracea*	4	9	54	7

^a Orthologous NPF gene ratio by comparing B. rapa and B. oleracea with B. napus genome. 0:1 represents NPF orthologs lost in B. rapa or B. oleracea genome, 1:0 represent lost in B. napus genome, 1:2, 1:3 and 1:5 represent different replication multiples in B. napus genome.

Table 2. Copy number variations (CNVs) of the NPF genes in 36 plant species

Organism Name	*NPF*1	*NPF*2	*NPF*3	*NPF*4	*NPF*5	*NPF*6	*NPF*7	*NPF*8	Total	Genome Size (Mb)
Arabidopsis lyrata (D)	3	14	1	9	17	4	3	5	56	202.97
*Arabidopsis* thaliana (D)	3	14	1	7	16	4	3	5	53	127.42
Aquilaria agallochum (D)	6	7	3	12	13	5	3	6	55	726.71
Brachypodium distachyon (M)	2	6	4	13	21	8	11	17	82	271.3
*Brassica rapa* (D)	4	23	3	9	23	7	5	8	82	401.93
*Brassica oleracea* (D)	4	15	2	8	26	6	5	4	70	554.98
*Brassica napus* (D)	10	30	6	16	63	15	10	19	169	976.19
Carica papaya (D)	4	14	3	8	12	8	6	4	59	370.42
Capsella rubella (D)	3	12	1	6	17	4	3	5	51	133.06
Citrus clementina (D)	9	7	3	9	17	6	4	4	59	301.37
Citrus sinensis (D)	8	7	3	10	17	6	4	4	59	319.23
Cuscuta campestris (D)	4	9	3	8	19	6	5	5	59	476.79
Eucalyptus grandis (D)	6	12	4	11	19	6	4	6	68	691.43
Fragaria vesca (D)	0	13	2	8	23	3	5	6	60	214.37
Glycine max (D)	13	14	6	22	41	11	14	13	134	927.71
Gossypium raimondii (D)	7	10	4	14	14	11	7	8	75	773.77
Linum usitatissimum (D)	12	7	4	14	25	9	11	10	92	316.17
Malus domestica (D)	2	34	4	21	44	17	8	9	139	1,874.77
Manihot esculenta (D)	7	12	6	10	23	7	5	5	75	292.1
Medicago truncatula (D)	8	12	3	14	25	8	9	1	80	412.92
Oryza sativa (M)	3	6	5	12	29	6	11	21	93	389.75
Phaseolus vulgaris (D)	8	11	3	12	22	5	7	6	74	521.08
Populus trichocarpa (D)	15	9	5	12	26	6	5	7	85	434.29
Prunus persica (D)	2	15	1	8	16	5	5	5	57	214.22
Ricinus communis (D)	5	20	3	7	13	5	4	3	60	350.62
Setaria italica (M)	4	11	8	16	19	7	12	21	98	405.87
Solanum tuberosum (D)	17	10	2	15	8	9	4	8	73	772.25
Solanum lycopersicum (D)	19	16	2	12	11	12	7	11	90	760.07
Sorghum bicolor (M)	4	8	7	16	22	6	9	19	91	709.35
Theobroma cacao (D)	4	14	3	10	19	7	4	5	66	345.99
Vitis vinifera (D)	4	7	2	6	21	5	4	3	52	486.2
Zea mays (M)	4	4	6	12	17	8	12	16	79	2,271.03
Amborella trichopoda (D)	1	5	2	7	15	4	3	7	45	706.60
Physcomitrella patens (L)	0	0	1	1	8	6	3	4	23	472.081
Selaginella moellendorffii (L)	0	0	4	4	11	6	5	16	46	212.315
Selaginella moellendorffii (L)	0	0	4	4	11	6	5	16	46	212.315

D dicots, M monocots, L lower plants

No competing interests reported.

NPF genes excavation and their expression response to vernalization and nitrogen deficiency in allotetraploid rapeseed

Status:

Version 1

Abstract

Background

Results

Conclusion

Figures

Background

Result

Discussion

Conclusion

Methods

Declarations

Abbreviations

References

Tables

Tables

Additional Declarations

Supplementary Files

Status:

Version 1