Multi-omics strategies uncover the molecular mechanisms of nitrogen, phosphorus and potassium deficiency responses in Brassica napus.

BACKGROUND: Nitrogen (N), phosphorus (P) and potassium (K) are critical macronutrients in crops, such that deficiency in any of N, P or K has substantial effects on crop growth. However, the specific commonalities of plant responses to different macronutrient deficiencies remain largely unknown. METHODS: Here, we assessed the phenotypic and physiological performances along with whole transcriptome and metabolomic profiles of rapeseed seedlings exposed to N, P and K deficiency stresses. RESULTS: Quantities of reactive oxygen species were significantly increased by all macronutrient deficiencies. N and K deficiencies resulted in more severe root development responses than P deficiency, as well as greater chlorophyll content reduction in leaves (associated with disrupted chloroplast structure). Transcriptome and metabolome analyses validated the macronutrient-specific responses, with more pronounced effects of N and P deficiencies on mRNAs, microRNAs (miRNAs), circular RNAs (circRNAs) and metabolites relative to K deficiency. Tissue-specific responses also occurred, with greater effects of macronutrient deficiencies on roots compared with shoots. We further uncovered a set of common responders with simultaneous roles in all three macronutrient deficiencies, including 112 mRNAs and 10 miRNAs involved in hormonal signaling, ion transport and oxidative stress in the root, and 33 mRNAs and 6 miRNAs with roles in abiotic stress response and photosynthesis in the shoot. 27 and seven common miRNA-mRNA pairs with role in miRNA-mediated regulation of oxidoreduction processes and ion transmembrane transport were identified in all three macronutrient deficiencies. No circRNA was responsive to three macronutrient deficiency stresses, but two common circRNAs were identified for two macronutrient deficiencies. Combined analysis of circRNAs, miRNAs and mRNAs suggested that two circRNAs act as decoys for miR156 and participate in oxidoreduction processes and transmembrane transport in both N- and P-deprived roots. Simultaneously, dramatic alterations of metabolites also occurred. Associations of RNAs with metabolites were observed, and suggested potential positive regulatory roles for tricarboxylic acids, azoles, carbohydrates, sterols and auxins, and negative regulatory roles for aromatic and aspartate amino acids, glucosamine-containing compounds, cinnamic acid, and nicotianamine in plant adaptation to macronutrient deficiency. CONCLUSIONS: Our findings revealed strategies to rescue rapeseed from macronutrient deficiency stress, including reducing the expression of non-essential genes and activating or enhancing the expression of anti-stress genes, aided by plant hormones, ion transporters and stress responders. The common responders to different macronutrient deficiencies identified could be targeted to enhance nutrient use efficiency in rapeseed.

LMI1-like genes involved in leaf margin development of Brassica napus.

In rapeseed (Brassica napus L.), leaf margins are variable and can be entire, serrate, or lobed. In our previous study, the lobed-leaf gene (LOBED-LEAF 1, BnLL1) was mapped to a 32.1 kb section of B. napus A10. Two LMI1-like genes, BnaA10g26320D and BnaA10g26330D, were considered the potential genes that controlled the lobed-leaf trait in rapeseed. In the present study, these two genes and another homologous gene (BnaC04g00850D) were transformed into Arabidopsis thaliana (L.) Heynh. plants to identify their functions. All three LMI1-like genes of B. napus produced serrate leaf margins. The expression analysis indicated that the expression level of BnaA10g26320D determined the difference between lobed- and entire-leaved lines in rapeseed. Therefore, it is likely that BnaA10g26320D corresponds to BnLL1.

Characterization of the quantitative trait locus OilA1 for oil content in Brassica napus.

Increasing seed oil content has become one of the most important breeding criteria in rapeseed (Brassica napus). However, oil content is a complex quantitative trait. QTL mapping in a double haploid population (SG population) emerging from a cross between a German (Sollux) and Chinese (Gaoyou) cultivars revealed one QTL for oil content on linkage group A1 (OilA1), which was mapped to a 17 cM genetic interval. To further validate and characterize the OilA1, we constructed a high-resolution map using B. rapa sequence resources and developed a set of near-isogenic lines (NILs) by employing a DH line SG-DH267 as donor and Chinese parent Gaoyou as recurrent background. The results showed highly conserved synteny order between B. rapa and B. napus within the linkage group A1 and revealed a possible centromere region between two markers ZAASA1-38 and NTP3 (2.5 cM). OilA1 was firstly validated by 250 BC5F2 plants and was confirmed in a 10.6 cM interval between the markers ZAASA1-47 and ZAASA1-77. Further substitution mapping was conducted by using two generations of QTL-NILs, 283 lines from eight BC5F3:4 families and 428 plants from six BC5F4 sub-NILs and thus narrowed the OilA1 interval to 6.9 cM and 4.3 cM (1.4 Mb), respectively. Field investigations with two replications using homozygous BC5F3:4 sister sub-NILs indicated that NILs, which carry a Sollux chromosome segment across the target region showed significant higher oil content (1.26 %, p < 0.001) than their sister NILs containing Gaoyou chromosome. The OilA1 locus is of particular interest for breeding purpose in China because 80 % of Chinese cultivars do not carry this desirable allele.

Development and validation of functional kompetitive allele-specific PCR markers for herbicide resistance in Brassica napus.

Effective weed control in the field is essential for maintaining favorable growing conditions and rapeseed yields. Sulfonylurea herbicides are one kind of most widely used herbicides worldwide, which control weeds by inhibiting acetolactate synthase (ALS). Molecular markers have been designed from polymorphic sites within the sequences of ALS genes, aiding marker-assisted selection in breeding herbicide-resistant rapeseed cultivars. However, most of them are not breeder friendly and have relatively limited application due to higher costs and lower throughput in the breeding projects. The aims of this study were to develop high throughput kompetitive allele-specific PCR (KASP) assays for herbicide resistance. We first cloned and sequenced BnALS1 and BnALS3 genes from susceptible cultivars and resistant 5N (als1als1/als3als3 double mutant). Sequence alignments of BnALS1 and BnALS3 genes for cultivars and 5N showed single nucleotide polymorphisms (SNPs) at positions 1676 and 1667 respectively. These two SNPs for BnALS1 and BnALS3 resulted in amino acid substitutions and were used to develop a KASP assay. These functional markers were validated in three distinct BC1F2 populations. The KASP assay developed in this study will be valuable for the high-throughput selection of elite materials with high herbicide resistance in rapeseed breeding programs.

Molecular mapping of Arabidopsis thaliana lipid-related orthologous genes in Brassica napus.

Quantitative Trait Loci (QTL) for oil content has been previously analyzed in a SG-DH population from a cross between a Chinese cultivar and a European cultivar of Brassica napus. Eight QTL with additive and epistatic effects, and with environmental interactions were evaluated. Here we present an integrated linkage map of this population predominantly based on informative markers derived from Brassica sequences, including 249 orthologous A. thaliana genes, where nearly half (112) are acyl lipid metabolism related genes. Comparative genomic analysis between B. napus and A. thaliana revealed 33 colinearity regions. Each of the conserved A. thaliana segments is present two to six times in the B. napus genome. Approximately half of the mapped lipid-related orthologous gene loci (76/137) were assigned in these conserved colinearity regions. QTL analysis for seed oil content was performed using the new map and phenotypic data from 11 different field trials. Nine significant QTL were identified on linkage groups A1, A5, A7, A9, C2, C3, C6 and C8, together explaining 57.79% of the total phenotypic variation. A total of 14 lipid related candidate gene loci were located in the confidence intervals of six of these QTL, of which ten were assigned in the conserved colinearity regions and felled in the most frequently overlapped QTL intervals. The information obtained from this study demonstrates the potential role of the suggested candidate genes in rapeseed kernel oil accumulation.

CRISPR/Cas9-Mediated Gene Editing of BnFAD2 and BnFAE1 Modifies Fatty Acid Profiles in Brassica napus.

Fatty acid (FA) composition determines the quality of oil from oilseed crops, and thus is a major target for genetic improvement. FAD2 (Fatty acid dehydrogenase 2) and FAE1 (fatty acid elongase 1) are critical FA synthetic genes, and have been the focus of genetic manipulation to alter fatty acid composition in oilseed plants. In this study, to improve the nutritional quality of rapeseed cultivar CY2 (about 50% oil content; of which 40% erucic acid), we generated novel knockout plants by CRISPR/Cas9 mediated genome editing of BnFAD2 and BnFAE1 genes. Two guide RNAs were designed to target one copy of the BnFAD2 gene and two copies of the BnFAE1 gene, respectively. A number of lines with mutations at three target sites of BnFAD2 and BnFAE1 genes were identified by sequence analysis. Three of these lines showed mutations in all three target sites of the BnFAD2 and BnFAE1 genes. Fatty acid composition analysis of seeds revealed that mutations at all three sites resulted in significantly increased oleic acid (70-80%) content compared with that of CY2 (20%), greatly reduced erucic acid levels and slightly decreased polyunsaturated fatty acids content. Our results confirmed that the CRISPR/Cas9 system is an effective tool for improving this important trait.