Low iron ameliorates the salinity-induced growth cessation of seminal roots in wheat seedlings.

Wheat plants are ubiquitously simultaneously exposed to salinity and limited iron availability caused by soil saline-alkalisation. Through this study, we found that both low Fe and NaCl severely inhibited the growth of seminal roots in wheat seedlings; however, sufficient Fe caused greater growth cessation of seminal roots than low Fe under salt stress. Low Fe improved the root meristematic division activity, not altering the mature cell sizes compared with sufficient Fe under salt stress. Foliar Fe spray and split-root experiments showed that low Fe-alleviating the salinity-induced growth cessation of seminal roots was dependent on local low Fe signals in the roots. Ionomics combined with TEM/X-ray few differences in the root Na+ uptake and vacuolar Na+ sequestration between two Fe levels under salt stress. Phytohormone profiling and metabolomics revealed salinity-induced overaccumulation of ACC/ethylene and tryptophan/auxin in the roots under sufficient Fe than under low Fe. Differential gene expression, pharmacological inhibitor addition and the root growth performance of transgenic wheat plants revealed that the rootward auxin efflux and was responsible for the low Fe-mediated amelioration of the salinity-induced growth cessation of seminal roots. Our findings will provide novel insights into the modulation of crop root growth under salt stress.

Combined morpho-physiological, ionomic and transcriptomic analyses reveal adaptive responses of allohexaploid wheat (Triticum aestivum L.) to iron deficiency.

BACKGROUND: Plants worldwide are often stressed by low Fe availability around the world, especially in aerobic soils. Therefore, the plant growth, seed yield, and quality of crop species are severely inhibited under Fe deficiency. Fe metabolism in plants is controlled by a series of complex transport, storage, and regulatory mechanisms in cells. Allohexaploid wheat (Triticum aestivum L.) is a staple upland crop species that is highly sensitive to low Fe stresses. Although some studies have been previously conducted on the responses of wheat plants to Fe deficiency, the key mechanisms underlying adaptive responses are still unclear in wheat due to its large and complex genome. RESULTS: Transmission electron microscopy showed that the chloroplast structure was severely damaged under Fe deficiency. Paraffin sectioning revealed that the division rates of meristematic cells were reduced, and the sizes of elongated cells were diminished. ICP-MS-assisted ionmics analysis showed that low-Fe stress significantly limited the absorption of nutrients, including N, P, K, Ca, Mg, Fe, Mn, Cu, Zn, and B nutrients. High-throughput transcriptome sequencing identified 378 and 2,619 genome-wide differentially expressed genes (DEGs) were identified in the shoots and roots between high-Fe and low-Fe conditions, respectively. These DEGs were mainly involved in the Fe chelator biosynthesis, ion transport, photosynthesis, amino acid metabolism, and protein synthesis. Gene coexpression network diagrams indicated that TaIRT1b-4A, TaNAS2-6D, TaNAS1a-6A, TaNAS1-6B, and TaNAAT1b-1D might function as key regulators in the adaptive responses of wheat plants to Fe deficiency. CONCLUSIONS: These results might help us fully understand the morpho-physiological and molecular responses of wheat plants to low-Fe stress, and provide elite genetic resources for the genetic modification of efficient Fe use.

Multiomics reveals an essential role of long-distance translocation in regulating plant cadmium resistance and grain accumulation in allohexaploid wheat (Triticum aestivum).

Cadmium (Cd) is a highly toxic heavy metal that readily enters cereals, such as wheat, via the roots and is translocated to the shoots and grains, thereby posing high risks to human health. However, the vast and complex genome of allohexaploid wheat makes it challenging to understand Cd resistance and accumulation. In this study, a Cd-resistant cultivar of wheat, 'ZM1860', and a Cd-sensitive cultivar, 'ZM32', selected from a panel of 442 accessions, exhibited significantly different plant resistance and grain accumulation. We performed an integrated comparative analysis of the morpho-physiological traits, ionomic and phytohormone profiles, genomic variations, transcriptomic landscapes, and gene functionality in order to identify the mechanisms underlying these differences. Under Cd toxicity, 'ZM1860' outperformed 'ZM32', which showed more severe leaf chlorosis, poorer root architecture, higher accumulation of reactive oxygen species, and disordered phytohormone homeostasis. Ionomics showed that 'ZM32' had a higher root-to-shoot translocation coefficient of Cd and accumulated more Cd in the grains than 'ZM1860'. Whole-genome re-sequencing (WGS) and transcriptome sequencing identified numerous DNA variants and differentially expressed genes involved in abiotic stress responses and ion transport between the two genotypes. Combined ionomics, transcriptomics, and functional gene analysis identified the plasma membrane-localized heavy metal ATPase TaHMA2b-7A as a crucial Cd exporter regulating long-distance Cd translocation in wheat. WGS- and PCR-based analysis of sequence polymorphisms revealed a 25-bp InDel site in the promoter region of TaHMA2b-7A, and this was probably responsible for the differential expression. Our multiomics approach thus enabled the identification of a core transporter involved in long-distance Cd translocation in wheat, and it may provide an elite genetic resource for improving plant Cd resistance and reducing grain Cd accumulation in wheat and other cereal crops.

Genome-Scale Investigation of GARP Family Genes Reveals Their Pivotal Roles in Nutrient Stress Resistance in Allotetraploid Rapeseed.

The GARP genes are plant-specific transcription factors (TFs) and play key roles in regulating plant development and abiotic stress resistance. However, few systematic analyses of GARPs have been reported in allotetraploid rapeseed (Brassica napus L.) yet. In the present study, a total of 146 BnaGARP members were identified from the rapeseed genome based on the sequence signature. The BnaGARP TFs were divided into five subfamilies: ARR, GLK, NIGT1/HRS1/HHO, KAN, and PHL subfamilies, and the members within the same subfamilies shared similar exon-intron structures and conserved motif configuration. Analyses of the Ka/Ks ratios indicated that the GARP family principally underwent purifying selection. Several cis-acting regulatory elements, essential for plant growth and diverse biotic and abiotic stresses, were identified in the promoter regions of BnaGARPs. Further, 29 putative miRNAs were identified to be targeting BnaGARPs. Differential expression of BnaGARPs under low nitrate, ammonium toxicity, limited phosphate, deficient boron, salt stress, and cadmium toxicity conditions indicated their potential involvement in diverse nutrient stress responses. Notably, BnaA9.HHO1 and BnaA1.HHO5 were simultaneously transcriptionally responsive to these nutrient stresses in both hoots and roots, which indicated that BnaA9.HHO1 and BnaA1.HHO5 might play a core role in regulating rapeseed resistance to nutrient stresses. Therefore, this study would enrich our understanding of molecular characteristics of the rapeseed GARPs and will provide valuable candidate genes for further in-depth study of the GARP-mediated nutrient stress resistance in rapeseed.

Integrated ionomic and transcriptomic dissection reveals the core transporter genes responsive to varying cadmium abundances in allotetraploid rapeseed.

BACKGROUND: Oilseed rape (B. napus L.) has great potential for phytoremediation of cadmium (Cd)-polluted soils due to its large plant biomass production and strong metal accumulation. Soil properties and the presence of other soluble compounds or ions, cause a heterogeneous distribution of Cd. RESULTS: The aim of our study was to reveal the differential responses of B. napus to different Cd abundances. Herein, we found that high Cd (50 µM) severely inhibited the growth of B. napus, which was not repressed by low Cd (0.50 µM) under hydroponic culture system. ICP-MS assays showed that the Cd2+ concentrations in both shoots and roots under 50 µM Cd were over 10 times higher than those under 0.50 µM Cd. Under low Cd, the concentrations of only shoot Ca2+/Mn2+ and root Mn2+ were obviously changed (both reduced); under high Cd, the concentrations of most cations assayed were significantly altered in both shoots and roots except root Ca2+ and Mg2+. High-throughput transcriptomic profiling revealed a total of 18,021 and 1408 differentially expressed genes under high Cd and low Cd conditions, respectively. The biological categories related to the biosynthesis of plant cell wall components and response to external stimulus were over-accumulated under low Cd, whereas the terms involving photosynthesis, nitrogen transport and response, and cellular metal ion homeostasis were highly enriched under high Cd. Differential expression of the transporters responsible for Cd uptake (NRAMPs), transport (IRTs and ZIPs), sequestration (HMAs, ABCs, and CAXs), and detoxification (MTPs, PCR, MTs, and PCSs), and some other essential nutrient transporters were investigated, and gene co-expression network analysis revealed the core members of these Cd transporters. Some Cd transporter genes, especially NRAMPs and IRTs, showed opposite responsive patterns between high Cd and low Cd conditions. CONCLUSIONS: Our findings would enrich our understanding of the interaction between essential nutrients and Cd, and might also provide suitable gene resources and important implications for the genetic improvement of plant Cd accumulation and resistance through molecular engineering of these core genes under varying Cd abundances in soils.

Genome-wide identification of Brassicaceae B-BOX genes and molecular characterization of their transcriptional responses to various nutrient stresses in allotetraploid rapeseed.

BACKGROUND: B-box (BBX) genes play important roles in plant growth regulation and responses to abiotic stresses. The plant growth and yield production of allotetraploid rapeseed is usually hindered by diverse nutrient stresses. However, no systematic analysis of Brassicaceae BBXs and the roles of BBXs in the regulation of nutrient stress responses have not been identified and characterized previously. RESULTS: In this study, a total of 536 BBXs were identified from nine brassicaceae species, including 32 AtBBXs, 66 BnaBBXs, 41 BoBBXs, 43 BrBBXs, 26 CrBBXs, 81 CsBBXs, 52 BnBBXs, 93 BjBBXs, and 102 BcBBXs. Syntenic analysis showed that great differences in the gene number of Brassicaceae BBXs might be caused by genome duplication. The BBXs were respectively divided into five subclasses according to their phylogenetic relationships and conserved domains, indicating their diversified functions. Promoter cis-element analysis showed that BBXs probably participated in diverse stress responses. Protein-protein interactions between BnaBBXs indicated their functions in flower induction. The expression profiles of BnaBBXs were investigated in rapeseed plants under boron deficiency, boron toxicity, nitrate limitation, phosphate shortage, potassium starvation, ammonium excess, cadmium toxicity, and salt stress conditions using RNA-seq data. The results showed that different BnaBBXs showed differential transcriptional responses to nutrient stresses, and some of them were simultaneously responsive to diverse nutrient stresses. CONCLUSIONS: Taken together, the findings investigated in this study provided rich resources for studying Brassicaceae BBX gene family and enriched potential clues in the genetic improvement of crop stress resistance.

Multiomics reveal pivotal roles of sodium translocation and compartmentation in regulating salinity resistance in allotetraploid rapeseed.

The large size and complexity of the allotetraploid rapeseed (Brassica napus) genome present huge challenges for understanding salinity resistance in this important crop. In this study, we identified two rapeseed genotypes with significantly different degrees of salinity resistance and examined the underlying mechanisms using an integrated analysis of phenomics, ionomics, genomics, and transcriptomics. Under salinity, a higher accumulation of osmoregulation substances and better root-system architecture was observed in the resistant genotype, H159, than in the sensitive one, L339. A lower shoot Na+ concentration and a higher root vacuolar Na+ concentration indicated lower root-to-shoot translocation and higher compartmentation in H159 than in L339. Whole-genome re-sequencing (WGRS) and transcriptome sequencing identified numerous DNA variants and differentially expressed genes involved in abiotic stress responses and ion transport. Combining ionomics with transcriptomics identified plasma membrane-localized BnaC2.HKT1;1 and tonoplast-localized BnaC5.NHX2 as the central factors regulating differential root xylem unloading and vacuolar sequestration of Na+ between the two genotypes. Identification of polymorphisms by WGRS and PCR revealed two polymorphic MYB-binding sites in the promoter regions that might determine the differential gene expression of BnaC2.HKT1;1 and BnaC5.NHX2. Our multiomics approach thus identified core transporters involved in Na+ translocation and compartmentation that regulate salinity resistance in rapeseed. Our results may provide elite gene resources for the improvement of salinity resistance in this crop, and our multiomics approach can be applied to other similar studies.

Genomic identification of nitrogen assimilation-related genes and transcriptional characterization of their responses to nitrogen in allotetraploid rapeseed.

BACKGROUND: Nitrogen (N) is an essential macronutrient to maintain plant growth and development. Plants absorb nitrate-N or ammonium-N in the environment and undergo reduction reactions catalyzed by nitrate reductase (NR), nitrite reductase (NIR), glutamine synthetase (GS), and glutamine oxoglutarate aminotransferase (GOGAT) within plants. METHODS AND RESULTS: A total of 42 N assimilation-related genes (NAG) members were identified in rapeseed. Darwin's evolutionary pressure analysis showed that rapeseed NAGs underwent purification selection. Cis-element analysis revealed differences in the transcriptional regulation of NAGs between Arabidopsis and rapeseed. Expression analyses revealed that NRs were expressed mainly in old leaves, NIRs were expressed mainly in old leaves and lower stem peels, while the expression situation between different subfamilies of GSs and GOGATs was more complicated. CONCLUSIONS: Differential expression of NAGs suggested that they might be involved in abiotic stresses. The above results greatly enriched our understanding of NAGs' molecular characteristics and provided central gene resources for NAGs-mediated NUE improvement in rapeseed.

Comprehensive dissection into morpho-physiologic responses, ionomic homeostasis, and transcriptomic profiling reveals the systematic resistance of allotetraploid rapeseed to salinity.

BACKGROUND: Salinity severely inhibit crop growth, yield, and quality worldwide. Allotetraploid rapeseed (Brassica napus L.), a major glycophyte oil crop, is susceptible to salinity. Understanding the physiological and molecular strategies of rapeseed salinity resistance is a promising and cost-effective strategy for developing highly resistant cultivars. RESULTS: First, early leaf senescence was identified and root system growth was inhibited in rapeseed plants under severe salinity conditions. Electron microscopic analysis revealed that 200 mM NaCl induced fewer leaf trichomes and stoma, cell plasmolysis, and chloroplast degradation. Primary and secondary metabolite assays showed that salinity led to an obviously increased anthocyanin, osmoregulatory substances, abscisic acid, jasmonic acid, pectin, cellulose, reactive oxygen species, and antioxidant activity, and resulted in markedly decreased photosynthetic pigments, indoleacetic acid, cytokinin, gibberellin, and lignin. ICP-MS assisted ionomics showed that salinity significantly constrained the absorption of essential elements, including the nitrogen, phosphorus, potassium, calcium, magnesium, iron, mangnese, copper, zinc, and boron nutrients, and induced the increase in the sodium/potassium ratio. Genome-wide transcriptomics revealed that the differentially expressed genes were involved mainly in photosynthesis, stimulus response, hormone signal biosynthesis/transduction, and nutrient transport under salinity. CONCLUSIONS: The high-resolution salt-responsive gene expression profiling helped the efficient characterization of central members regulating plant salinity resistance. These findings might enhance integrated comprehensive understanding of the morpho-physiologic and molecular responses to salinity and provide elite genetic resources for the genetic modification of salinity-resistant crop species.

Genome-wide identification of the amino acid permease genes and molecular characterization of their transcriptional responses to various nutrient stresses in allotetraploid rapeseed.

BACKGROUND: Nitrogen (N), referred to as a "life element", is a macronutrient essential for optimal plant growth and yield production. Amino acid (AA) permease (AAP) genes play pivotal roles in root import, long-distance translocation, remobilization of organic amide-N from source organs to sinks, and other environmental stress responses. However, few systematic analyses of AAPs have been reported in Brassica napus so far. RESULTS: In this study, we identified a total of 34 full-length AAP genes representing eight subgroups (AAP1-8) from the allotetraploid rapeseed genome (AnAnCnCn, 2n = 4x = 38). Great differences in the homolog number among the BnaAAP subgroups might indicate their significant differential roles in the growth and development of rapeseed plants. The BnaAAPs were phylogenetically divided into three evolutionary clades, and the members in the same subgroups had similar physiochemical characteristics, gene/protein structures, and conserved AA transport motifs. Darwin's evolutionary analysis suggested that BnaAAPs were subjected to strong purifying selection pressure. Cis-element analysis showed potential differential transcriptional regulation of AAPs between the model Arabidopsis and B. napus. Differential expression of BnaAAPs under nitrate limitation, ammonium excess, phosphate shortage, boron deficiency, cadmium toxicity, and salt stress conditions indicated their potential involvement in diverse nutrient stress responses. CONCLUSIONS: The genome-wide identification of BnaAAPs will provide a comprehensive insight into their family evolution and AAP-mediated AA transport under diverse abiotic stresses. The molecular characterization of core AAPs can provide elite gene resources and contribute to the genetic improvement of crop stress resistance through the modulation of AA transport.

Global Landscapes of the Na+/H+ Antiporter (NHX) Family Members Uncover their Potential Roles in Regulating the Rapeseed Resistance to Salt Stress.

Soil salinity is a main abiotic stress in agriculture worldwide. The Na+/H+ antiporters (NHXs) play pivotal roles in intracellular Na+ excretion and vacuolar Na+ compartmentalization, which are important for plant salt stress resistance (SSR). However, few systematic analyses of NHXs has been reported in allotetraploid rapeseed so far. Here, a total of 18 full-length NHX homologs, representing seven subgroups (NHX1-NHX8 without NHX5), were identified in the rapeseed genome (AnAnCnCn). Number variations of BnaNHXs might indicate their significantly differential roles in the regulation of rapeseed SSR. BnaNHXs were phylogenetically divided into three evolutionary clades, and the members in the same subgroups had similar physiochemical characteristics, gene/protein structures, and conserved Na+ transport motifs. Darwin´s evolutionary pressure analysis suggested that BnaNHXs suffered from strong purifying selection. The cis-element analysis revealed the differential transcriptional regulation of NHXs between the model Arabidopsis and B. napus. Differential expression of BnaNHXs under salt stress, different nitrogen forms (ammonium and nitrate), and low phosphate indicated their potential involvement in the regulation of rapeseed SSR. Global landscapes of BnaNHXs will give an integrated understanding of their family evolution and molecular features, which will provide elite gene resources for the genetic improvement of plant SSR through regulating the NHX-mediated Na+ transport.

Genome-Wide Differential DNA Methylation and miRNA Expression Profiling Reveals Epigenetic Regulatory Mechanisms Underlying Nitrogen-Limitation-Triggered Adaptation and Use Efficiency Enhancement in Allotetraploid Rapeseed.

Improving crop nitrogen (N) limitation adaptation (NLA) is a core approach to enhance N use efficiency (NUE) and reduce N fertilizer application. Rapeseed has a high demand for N nutrients for optimal plant growth and seed production, but it exhibits low NUE. Epigenetic modification, such as DNA methylation and modification from small RNAs, is key to plant adaptive responses to various stresses. However, epigenetic regulatory mechanisms underlying NLA and NUE remain elusive in allotetraploid B. napus. In this study, we identified overaccumulated carbohydrate, and improved primary and lateral roots in rapeseed plants under N limitation, which resulted in decreased plant nitrate concentrations, enhanced root-to-shoot N translocation, and increased NUE. Transcriptomics and RT-qPCR assays revealed that N limitation induced the expression of NRT1.1, NRT1.5, NRT1.7, NRT2.1/NAR2.1, and Gln1;1, and repressed the transcriptional levels of CLCa, NRT1.8, and NIA1. High-resolution whole genome bisulfite sequencing characterized 5094 differentially methylated genes involving ubiquitin-mediated proteolysis, N recycling, and phytohormone metabolism under N limitation. Hypermethylation/hypomethylation in promoter regions or gene bodies of some key N-metabolism genes might be involved in their transcriptional regulation by N limitation. Genome-wide miRNA sequencing identified 224 N limitation-responsive differentially expressed miRNAs regulating leaf development, amino acid metabolism, and plant hormone signal transduction. Furthermore, degradome sequencing and RT-qPCR assays revealed the miR827-NLA pathway regulating limited N-induced leaf senescence as well as the miR171-SCL6 and miR160-ARF17 pathways regulating root growth under N deficiency. Our study provides a comprehensive insight into the epigenetic regulatory mechanisms underlying rapeseed NLA, and it will be helpful for genetic engineering of NUE in crop species through epigenetic modification of some N metabolism-associated genes.

Genome-wide identification of the cation/proton antiporter (CPA) gene family and functional characterization of the key member BnaA05.NHX2 in allotetraploid rapeseed.

Rapeseed (Brassica napus L.) is susceptible to nutrient stresses during growth and development; however, the CPA (cation proton antiporter) family genes have not been identified in B. napus and their biological functions remain unclear. This study was aimed to identify the molecular characteristics of rapeseed CPAs and their transcriptional responses to multiple nutrient stresses. Through bioinformatics analysis, 117 BnaCPAs, consisting of three subfamilies: Na+/H+ antiporter (NHX), K+ efflux antiporter (KEA), and cation/H+ antiporter (CHX), were identified in the rapeseed genome. Transcriptomic profiling showed that BnaCPAs, particularly BnaNHXs, were transcriptionally responsive to diverse nutrient stresses, including Cd toxicity, K starvation, salt stress, NH4+ toxicity, and low Pi. We found that the salt tolerance of the transgenic rapeseed lines overexpressing BnaA05.NHX2 was significantly higher than that of wild type. Subcellular localization showed that BnaA05.NHX2 was localized on the tonoplast, and TEM combined with X-ray energy spectrum analysis revealed that the vacuolar Na+ concentrations of the BnaA05.NHX2-overexpressing rapeseed plants were significantly higher than those of wild type. The findings of this study will provide insights into the complexity of the BnaCPA family and a valuable resource to explore the in-depth functions of CPAs in B. napus.

Morpho-physiological, Genomic, and Transcriptional Diversities in Response to Potassium Deficiency in Rapeseed (Brassica napus L.) Genotypes.

Variations in the resistance to potassium (K) deficiency among rapeseed genotypes emphasize complicated regulatory mechanisms. In this study, a low-K-sensitivity accession (L49) responded to K deficiency with smaller biomasses, severe leaf chlorosis, weaker photosynthesis ability, and deformed stomata morphology compared to a low-K resistant accession (H280). H280 accumulated more K+ than L49 under low K. Whole-genome resequencing (WGS) revealed a total of 5,538,622 single nucleotide polymorphisms (SNPs) and 859,184 insertions/deletions (InDels) between H280 and L49. RNA-seq identified more differentially expressed K+ transporter genes with higher expression in H280 than in L49 under K deficiency. Based on the K+ profiles, differential expression profiling, weighted gene coexpression network analysis, and WGS data between H280 and L49, BnaC4.AKT1 was proposed to be mainly responsible for root K absorption-mediated low K resistance. BnaC4.AKT1 was expressed preferentially in the roots and localized on the plasma membrane. An SNP and an InDel found in the promoter region of BnaC4.AKT1 were proposed to be responsible for its differential expression between rapeseed genotypes. This study identified a gene resource for improving low-K resistance. It also facilitates an integrated knowledge of the differential physiological and transcriptional responses to K deficiency in rapeseed genotypes.

Multiomic Analysis Reveals Core Regulatory Mechanisms underlying Steroidal Glycoalkaloid Metabolism in Potato Tubers.

Steroidal glycoalkaloids (SGAs) present in germinated potato tubers are toxic; however, the mechanisms underlying SGA metabolism are poorly understood. Therefore, integrated transcriptome, metabolome, and hormone analyses were performed in this study to identify and characterize the key regulatory genes, metabolites, and phytohormones related to glycoalkaloid regulation. Based on transcriptome sequencing of bud eyes of germinated and dormant potato tubers, a total of 6260 differentially expressed genes were identified, which were mainly responsible for phytohormone signal transduction, carbohydrate metabolism, and secondary metabolite biosynthesis. Two TCP14 genes were identified as the core transcription factors that potentially regulate SGA synthesis. Metabolite analysis indicated that 149 significantly different metabolites were detected, and they were enriched in metabolic and biosynthetic pathways of secondary metabolites. In these pathways, the α-solanine content was increased and the expression of genes related to glycoalkaloid biosynthesis was upregulated. Levels of gibberellin and jasmonic acid were increased, whereas that of abscisic acid was decreased. This study lays a foundation for investigating the biosynthesis and regulation of SGAs and provides the reference for the production and consumption of potato tubers.

Genome-scale identification of plant defensin (PDF) family genes and molecular characterization of their responses to diverse nutrient stresses in allotetraploid rapeseed.

Plant defensins (PDFs), short peptides with strong antibacterial activity, play important roles in plant growth, development, and stress resistance. However, there are few systematic analyses on PDFs in Brassica napus. Here, bioinformatics methods were used to identify genome-wide PDFs in Brassica napus, and systematically analyze physicochemical properties, expansion pattern, phylogeny, and expression profiling of BnaPDFs under diverse nutrient stresses. A total of 37 full-length PDF homologs, divided into two subgroups (PDF1s and PDF2s), were identified in the rapeseed genome. A total of two distinct clades were identified in the BnaPDF phylogeny. Clade specific conserved motifs were identified within each clade respectively. Most BnaPDFs were proved to undergo powerful purified selection. The PDF members had enriched cis-elements related to growth and development, hormone response, environmental stress response in their promoter regions. GO annotations indicate that the functional pathways of BnaPDFs are mainly involved in cells killing and plant defense responses. In addition, bna-miRNA164 and bna-miRNA172 respectively regulate the expression of their targets BnaA2.PDF2.5 and BnaC7.PDF2.6. The expression patterns of BnaPDFs were analyzed in different tissues. BnaPDF1.2bs was mainly expressed in the roots, whereas BnaPDF2.2s and BnaPDF2.3s were both expressed in stamen, pericarp, silique, and stem. However, the other BnaPDF members showed low expression levels in various tissues. Differential expression of BnaPDFs under nitrate limitation, ammonium excess, phosphorus starvation, potassium deficiency, cadmium toxicity, and salt stress indicated that they might participate in different nutrient stress resistance. The genome-wide identification and characterization of BnaPDFs will enrich understanding of their molecular characteristics and provide elite gene resources for genetic improvement of rapeseed resistance to nutrient stresses.