Integrating GWAS, RNA-Seq and functional analysis revealed that BnaA02.SE mediates silique elongation by affecting cell proliferation and expansion in Brassica napus.

Rapeseed (Brassica napus) silique is the major carbohydrate source for seed development, and the final silique length has attracted great attention from breeders. However, no studies had focused on the dynamic character of silique elongation length (SEL). Here, the dynamic SEL investigation in a natural population including 588 lines over two years indicate that dynamic SEL during 0-20 days after flowering was the most essential stage associated with seed number per silique (SPS) and thousand seed weight (TSW). Then, nine loci were identified to be associated with SEL based on GWAS analysis, among which five SNPs (over 50%) distributed on the A02 chromosome within 6.08 to 6.48 Mb. Subsequently, we screened 5078 differentially expressed genes between two extreme materials. An unknown protein, BnaA02.SE, was identified combining with GWAS and RNA-Seq analysis. Subcellular localization and expression profiles analysis demonstrated that BnaA02.SE is a chloroplast- and nucleus-localized protein mainly expressed in pericarps and leaves. Furthermore, transgenic verification and dynamic cytological observation reveal that overexpressed BnaA02.SE can promote silique elongation by regulating JA and IAA contents, affecting cell proliferation and expansion, respectively, and finally enhance seed yield by influencing SPS and TSW. Haplotype analysis reveal that the homologs of BnaA02.SE may also be involved in silique elongation regulation. Our findings provided comprehensive insights into a newly SEL trait, and cloned the first gene (BnaA02.SE) controlling silique elongation in B. napus. The identified BnaA02.SE and its homologs can offer a valuable target for improving B. napus yield.

Spatio-temporal transcriptome profiling and subgenome analysis in Brassica napus.

Brassica napus is an important oil crop and an allotetraploid species. However, the detailed analysis of gene function and homoeologous gene expression in all tissues at different developmental stages was not explored. In this study, we performed a global transcriptome analysis of 24 vegetative and reproductive tissues at six developmental stages (totally 111 tissues). These samples were clustered into eight groups. The gene functions of silique pericarp were similar to roots, stems and leaves. In particular, glucosinolate metabolic process was associated with root and silique pericarp. Genes involved in protein phosphorylation were often associated with stamen, anther and the early developmental stage of seeds. Transcription factor (TF) genes were more specific than structural genes. A total of 17 100 genes that were preferentially expressed in one tissue (tissue-preferred genes, TPGs), including 889 TFs (5.2%), were identified in the 24 tissues. Some TPGs were identified as hub genes in the co-expression network analysis, and some TPGs in different tissues were involved in different hormone pathways. About 67.0% of the homoeologs showed balanced expression, whereas biased expression of homoeologs was associated with structural divergence. In addition, the spatiotemporal expression of homoeologs was related to the presence of transposable elements (TEs) and regulatory elements (REs); more TEs and fewer REs in the promoters resulted in divergent expression in different tissues. This study provides a valuable transcriptional map for understanding the growth and development of B. napus, for identifying important genes for future crop improvement, and for exploring gene expression patterns in the B. napus.

LESION MIMIC MUTANT 1 confers basal resistance to Sclerotinia sclerotiorum in rapeseed via a salicylic acid-dependent pathway.

Rapeseed (Brassica napus) is a major edible oilseed crop consumed worldwide. However, its yield is seriously affected by infection from the broad-spectrum non-obligate pathogen Sclerotinia sclerotiorum due to a lack of highly resistant germplasm. Here, we identified a Sclerotinia-resistant and light-dependent lesion mimic mutant from an ethyl methanesulfonate-mutagenized population of the rapeseed inbred Zhongshuang 11 (ZS11) named lesion mimic mutant 1 (lmm1). The phenotype of lmm1 is controlled by a single recessive gene, named LESION MIMIC MUTANT 1 (LMM1), which mapped onto chromosome C04 by bulked segregant analysis within a 2.71-Mb interval. Histochemical analysis indicated that H2O2 strongly accumulated and cell death occurred around the lesion mimic spots. Among 877 differentially expressed genes (DEGs) between ZS11 and lmm1 leaves, 188 DEGs were enriched in the defense response, including 95 DEGs involved in systemic acquired resistance, which is consistent with the higher salicylic acid levels in lmm1. Combining bulked segregant analysis and transcriptome analysis, we identified a significantly up-regulated gene, BnaC4.PR2, which encodes ß-1,3-glucanase, as the candidate gene for LMM1. Overexpression of BnaC4.PR2 may induce a reactive oxygen species burst to trigger partial cell death and systemic acquired resistance. Our study provides a new genetic resource for S. sclerotiorum resistance as well as new insights into disease resistance breeding in B. napus.

Genome-Wide Identification and Posttranscriptional Regulation Analyses Elucidate Roles of Key Argonautes and Their miRNA Triggers in Regulating Complex Yield Traits in Rapeseed.

Argonautes (AGOs) interact with microRNAs (miRNAs) to form the RNA-induced silencing complex (RISC), which can posttranscriptionally regulate the expression of targeted genes. To date, however, the AGOs and their miRNA triggers remain elusive in rapeseed (Brassica napus). Here, we systematically performed a phylogenetic analysis and examined the collinear relationships of the AGOs among four Brassicaceae species. Their physicochemical properties, gene structures, and expression patterns among 81 tissues from multiple materials and developmental stages were further analyzed. Additionally, their posttranscriptional regulation was analyzed using psRNATarget prediction, miRNA-/mRNA-Seq analyses, and a qRT-PCR verification. We finally identified 10 AtAGOs, 13 BolAGOs, 11 BraAGOs, and 24 BnaAGOs. An expression analysis of the BnaAGOs in the B. napus cultivar ZS11, as well as genotypes with extreme phenotypes in various yield-related traits, revealed the conservation and diversity of these genes. Furthermore, we speculated the posttranscriptional regulation of the B. napus miR168a-AGO1s and miR403-AGO2s modules. Combining miRNA-Seq and mRNA-Seq analyses, we found that the B. napus miR168a-AGO1s module may play an essential role in negatively regulating yield traits, whereas the miR403-AGO2s module positively impacts yield. This is the first attempt to comprehensively analyze the AGOs and their miRNA triggers in B. napus and provides a theoretical basis for breeding high-yielding varieties through the manipulation of the miRNA-AGOs modules.

Knockout of the lignin pathway gene BnF5H decreases the S/G lignin compositional ratio and improves Sclerotinia sclerotiorum resistance in Brassica napus.

Ferulate-5-hydroxylase is a key enzyme involved in the conversion of the guaiacyl monolignol to the syringyl monolignol in angiosperms. The monolignol ratio has been proposed to affect biomass recalcitrance and the resistance to plant disease. Stem rot caused by the fungus Sclerotinia sclerotiorum in Brassica napus causes severe losses in its production. To date, there is no information about the effect of the lignin monomer ratio on the resistance to S. sclerotiorum in B. napus. Four dominantly expressed ferulate-5-hydroxylase genes were concertedly knocked out by CRISPR/Cas9 in B. napus, and three mutant lines were generated. The S/G lignin compositional ratio was decreased compared to that of the wild type based on the results of Mӓule staining and 2D-NMR profiling in KO-7. The resistance to S. sclerotiorum in stems and leaves increased for the three f5h mutant lines compared with WT. Furthermore, we found that the stem strength of f5h mutant lines was significantly increased. Overall, we demonstrate for the first time that decreasing the S/G ratio by knocking out of the F5H gene improves S. sclerotiorum resistance in B. napus and increases stem strength.

Multi-omics analysis reveals the mechanism of seed coat color formation in Brassica rapa L.

KEYMESSAGE: Multi-omics analysis of the transcriptome, metabolome and genome identified major and minor loci and candidate genes for seed coat color and explored the mechanism of flavonoid metabolites biosynthesis in Brassica rapa. Yellow seed trait is considered an agronomically desirable trait with great potential for improving seed quality of Brassica crops. Mechanisms of the yellow seed trait are complex and not well understood. In this study, we performed an integrated metabolome, transcriptome and genome-wide association study (GWAS) on different B. rapa varieties to explore the mechanisms underlying the seed coat color formation. A total of 2,499 differentially expressed genes and 116 differential metabolites between yellow and black seeds with strong association with the flavonoid biosynthesis pathway was identified. In addition, 330 hub genes involved in the seed coat color formation, and the most significantly differential flavonoids biosynthesis were detected based on weighted gene co-expression network analysis. Metabolite GWAS analysis using the contents of 42 flavonoids in developing seeds of 159 B. rapa lines resulted in the identification of 1,626 quantitative trait nucleotides (QTNs) and 37 chromosomal intervals, including one major locus on chromosome A09. A combination of QTNs detection, transcriptome and functional analyses led to the identification of 241 candidate genes that were associated with different flavonoid metabolites. The flavonoid biosynthesis pathway in B. rapa was assembled based on the identified flavonoid metabolites and candidate genes. Furthermore, BrMYB111 members (BraA09g004490.3C and BraA06g034790.3C) involved in the biosynthesis of taxifolin were functionally analyzed in vitro. Our findings lay a foundation and provide a reference for systematically investigating the mechanism of seed coat color in B. rapa and in the other plants.

Comprehensive analysis of polygalacturonase genes offers new insights into their origin and functional evolution in land plants.

Polygalacturonase (PG) is a hydrolase that participates in pectin degradation, pod shattering and fruit softening. Here, we identified 2786 PG genes across 54 plants, which could be divided into three groups. Evolutionary analysis suggested that PG family originated from the charophyte green algae, and Subgroups A2-A4 evolved from the Subgroup A1 after the tracheophyte-angiosperm split. Whole-genome duplication was the major force leading to PG gene expansion. Interestingly, the PG genes continuously expanded in eudicots, whereas it contracted in monocots after the eudicot-monocot split. PG genes in Group A are expressed at high levels in floral organs, whereas genes in Groups B and C are expressed at high levels in various tissues. Moreover, three BnaPG15 members were found for their potential possibility in pod shattering in Brassica napus. Our results provide new insight into the evolutionary history of PG family, and their potentially functional role in plants.

Comparative transcriptomic analysis of seed coats with high and low lignin contents reveals lignin and flavonoid biosynthesis in Brassica napus.

BACKGROUND: Brassica napus L. (2n = 38, AACC) is one of the most important oil crops and sources of protein for animal feed worldwide. Lignin is a large molecule aromatic polymer and a major cell wall component. However, lignin in the seed coat reduces the availability and restricts the development of rapeseed cake. Therefore, it is critical to reduce the lignin content of the seed coat. Here, high-lignin (H-lignin) and low-lignin (L-lignin) content recombinant inbred lines (RILs) were selected from an RIL population for analysis. RESULTS: The cross-section results indicated that the seed coat of the H-lignin lines was thicker than that of the L-lignin lines, especially the palisade layer. The seed coats and embryos at 35, 40 and 46 days after flowering (DAF) were subjected to RNA sequencing (RNA-Seq), and the expression of the BnPAL and BnC4H gene families in the lignin pathway was significantly higher in the H-lignin seed coat than in the L-lignin seed coat. The Bn4CL gene family also showed this trend. In addition, among the genes related to plant hormone synthesis, BnaC02g01710D was upregulated and BnaA07g11700D and BnaC09g00190D were downregulated in H-lignin lines. Some transcription factors were upregulated, such as BnNAC080, BnNAC083, BnMYB9, BnMYB9-1, BnMYB60 and BnMYB60-1, while BnMYB91 was downregulated in H-lignin lines. Moreover, most genes of the flavonoid pathway, such as BnCHS and BnDFR, were strongly expressed in H-lignin seed coat. CONCLUSIONS: In Our study, some key genes such as hormone synthesis genes, transcription factors and miRNAs related to lignin and flavonoid biosynthesis were identified. A regulatory model of B. napus seed coat lignin was proposed. These results provide new insight into lignin and flavonoid biosynthesis in B. napus.

Comparative transcriptomic and metabolomic analyses of carotenoid biosynthesis reveal the basis of white petal color in Brassica napus.

MAIN CONCLUSION: The molecular mechanism underlying white petal color in Brassica napus was revealed by transcriptomic and metabolomic analyses. Rapeseed (Brassica napus L.) is one of the most important oilseed crops worldwide, but the mechanisms underlying flower color in this crop are known less. Here, we performed metabolomic and transcriptomic analyses of the yellow-flowered rapeseed cultivar 'Zhongshuang 11' (ZS11) and the white-flowered inbred line 'White Petal' (WP). The total carotenoid contents were 1.778-fold and 1.969-fold higher in ZS11 vs. WP petals at stages S2 and S4, respectively. Our findings suggest that white petal color in WP flowers is primarily due to decreased lutein and zeaxanthin contents. Transcriptome analysis revealed 10,116 differentially expressed genes with a fourfold or greater change in expression (P-value less than 0.001) in WP vs. ZS11 petals, including 1,209 genes that were differentially expressed at four different stages and 20 genes in the carotenoid metabolism pathway. BnNCED4b, encoding a protein involved in carotenoid degradation, was expressed at abnormally high levels in WP petals, suggesting it might play a key role in white petal formation. The results of qRT-PCR were consistent with the transcriptome data. The results of this study provide important insights into the molecular mechanisms of the carotenoid metabolic pathway in rapeseed petals, and the candidate genes identified in this study provide a resource for the creation of new B. napus germplasms with different petal colors.

Genome-wide association study and transcriptome comparison reveal novel QTL and candidate genes that control petal size in rapeseed.

Petal size determines the value of ornamental plants, and thus their economic value. However, the molecular mechanisms controlling petal size remain unclear in most non-model species. To identify quantitative trait loci and candidate genes controlling petal size in rapeseed (Brassica napus), we performed a genome-wide association study (GWAS) using data from 588 accessions over three consecutive years. We detected 16 significant single nucleotide polymorphisms (SNPs) associated with petal size, with the most significant SNPs located on chromosomes A05 and C06. A combination of GWAS and transcriptomic sequencing based on two accessions with contrasting differences in petal size identified 52 differentially expressed genes (DEGs) that may control petal size variation in rapeseed. In particular, the rapeseed gene BnaA05.RAP2.2, homologous to Arabidopsis RAP2.2, may be critical to the negative control of petal size through the ethylene signaling pathway. In addition, a comparison of petal epidermal cells indicated that petal size differences between the two contrasting accessions were determined mainly by differences in cell number. Finally, we propose a model for the control of petal size in rapeseed through ethylene and cytokinin signaling pathways. Our results provide insights into the genetic mechanisms regulating petal size in flowering plants.

Comparative transcriptome and metabolomic profiling reveal the complex mechanisms underlying the developmental dynamics of tobacco leaves.

Although the leaf is the most important photosynthetic organ in most plants, many of the molecular mechanisms underlying leaf developmental dynamics remain to be explored. To better understand the transcriptional regulatory mechanisms involved in leaf development, we conducted comparative transcriptomic and metabolomic analysis of leaves from seven positions on tobacco (Nicotiana tabacum) plants. A total of 35,622 unique differentially expressed genes and 79 metabolites were identified. A time-series expression analysis detected two interesting transcriptional profiles, one comprising 10,197 genes that displayed continual up-regulation during leaf development and another comprising 4696 genes that displayed continual down-regulation. Combining these data with co-expression network results identified four important regulatory networks involved in photorespiration and the tricarboxylic acid cycle; these networks may regulate carbon/nitrogen balance during leaf development. We also found that the transcription factor NtGATA5 acts as a hub associated with C and N metabolism and chloroplast development during leaf development through regulation of phytohormones. Furthermore, we investigated the transcriptional dynamics of genes involved in the auxin, cytokinin, and jasmonic acid biosynthesis and signaling pathways during tobacco leaf development. Overall, our study greatly expands the understanding of the regulatory network controlling developmental dynamics in plant leaves.

Genome-Wide Identification and Analysis of MKK and MAPK Gene Families in Brassica Species and Response to Stress in Brassica napus.

Mitogen-activated protein kinase (MAPK) cascades are common and conserved signal transduction pathways and play important roles in various biotic and abiotic stress responses and growth and developmental processes in plants. With the advancement of sequencing technology, more systematic genetic information is being explored. The work presented here focuses on two protein families in Brassica species: MAPK kinases (MKKs) and their phosphorylation substrates MAPKs. Forty-seven MKKs and ninety-two MAPKs were identified and extensively analyzed from two tetraploid (B. juncea and B. napus) and three diploid (B. nigra, B. oleracea, and B. rapa) Brassica species. Phylogenetic relationships clearly distinguished both MKK and MAPK families into four groups, labeled A-D, which were also supported by gene structure and conserved protein motif analysis. Furthermore, their spatial and temporal expression patterns and response to stresses (cold, drought, heat, and shading) were analyzed, indicating that BnaMKK and BnaMAPK transcript levels were generally modulated by growth, development, and stress signals. In addition, several protein interaction pairs between BnaMKKs and C group BnaMAPKs were detected by yeast two-hybrid assays, in which BnaMKK3 and BnaMKK9 showed strong interactions with BnaMAPK1/2/7, suggesting that interaction between BnaMKKs and C group BnaMAPKs play key roles in the crosstalk between growth and development processes and abiotic stresses. Taken together, our data provide a deeper foundation for the evolutionary and functional characterization of MKK and MAPK gene families in Brassica species, paving the way for unraveling the biological roles of these important signaling molecules in plants.

Metabolite Profiling and Transcriptome Analysis Provide Insight into Seed Coat Color in Brassica juncea.

The allotetraploid species Brassica juncea (mustard) is grown worldwide as oilseed and vegetable crops; the yellow seed-color trait is particularly important for oilseed crops. Here, to examine the factors affecting seed coat color, we performed a metabolic and transcriptomic analysis of yellow- and dark-seeded B. juncea seeds. In this study, we identified 236 compounds, including 31 phenolic acids, 47 flavonoids, 17 glucosinolates, 38 lipids, 69 other hydroxycinnamic acid compounds, and 34 novel unknown compounds. Of these, 36 compounds (especially epicatechin and its derivatives) accumulated significantly different levels during the development of yellow- and dark-seeded B. juncea. In addition, the transcript levels of BjuDFR, BjuANS,BjuBAN, BjuTT8, and BjuTT19 were closely associated with changes to epicatechin and its derivatives during seed development, implicating this pathway in the seed coat color determinant in B. juncea. Furthermore, we found numerous variations of sequences in the TT8A genes that may be associated with the stability of seed coat color in B. rapa, B. napus, and B. juncea, which might have undergone functional differentiation during polyploidization in the Brassica species. The results provide valuable information for understanding the accumulation of metabolites in the seed coat color of B. juncea and lay a foundation for exploring the underlying mechanism.

Genome-wide identification and comparative analysis of diacylglycerol kinase (DGK) gene family and their expression profiling in Brassica napus under abiotic stress.

BACKGROUND: Diacylglycerol kinases (DGKs) are signaling enzymes that play pivotal roles in response to abiotic and biotic stresses by phosphorylating diacylglycerol (DAG) to form phosphatidic acid (PA). However, no comprehensive analysis of the DGK gene family had previously been reported in B. napus and its diploid progenitors (B. rapa and B. oleracea). RESULTS: In present study, we identified 21, 10, and 11 DGK genes from B. napus, B. rapa, and B. oleracea, respectively, which all contained conserved catalytic domain and were further divided into three clusters. Molecular evolutionary analysis showed that speciation and whole-genome triplication (WGT) was critical for the divergence of duplicated DGK genes. RNA-seq transcriptome data revealed that, with the exception of BnaDGK4 and BnaDGK6, BnaDGK genes have divergent expression patterns in most tissues. Furthermore, some DGK genes were upregulated or downregulated in response to hormone treatment and metal ion (arsenic and cadmium) stress. Quantitative real-time PCR analysis revealed that different BnaDGK genes contribute to seed oil content. CONCLUSIONS: Together, our results indicate that DGK genes have diverse roles in plant growth and development, hormone response, and metal ion stress, and in determining seed oil content, and lay a foundation for further elucidating the roles of DGKs in Brassica species.

qPrimerDB: a thermodynamics-based gene-specific qPCR primer database for 147 organisms.

Real-time quantitative polymerase chain reaction (qPCR) is one of the most important methods for analyzing the expression patterns of target genes. However, successful qPCR experiments rely heavily on the use of high-quality primers. Various qPCR primer databases have been developed to address this issue, but these databases target only a few important organisms. Here, we developed the qPrimerDB database, founded on an automatic gene-specific qPCR primer design and thermodynamics-based validation workflow. The qPrimerDB database is the most comprehensive qPCR primer database available to date, with a web front-end providing gene-specific and pre-computed primer pairs across 147 important organisms, including human, mouse, zebrafish, yeast, thale cress, rice and maize. In this database, we provide 3331426 of the best primer pairs for each gene, based on primer pair coverage, as well as 47760359 alternative gene-specific primer pairs, which can be conveniently batch downloaded. The specificity and efficiency was validated for qPCR primer pairs for 66 randomly selected genes, in six different organisms, through qPCR assays and gel electrophoresis. The qPrimerDB database represents a valuable, timesaving resource for gene expression analysis. This resource, which will be routinely updated, is publically accessible at http://biodb.swu.edu.cn/qprimerdb.

BrassicaEDB: A Gene Expression Database for Brassica Crops.

The genus Brassica contains several economically important crops, including rapeseed (Brassica napus, 2n = 38, AACC), the second largest source of seed oil and protein meal worldwide. However, research in rapeseed is hampered because it is complicated and time-consuming for researchers to access different types of expression data. We therefore developed the Brassica Expression Database (BrassicaEDB) for the research community. In the current BrassicaEDB, we only focused on the transcriptome level in rapeseed. We conducted RNA sequencing (RNA-Seq) of 103 tissues from rapeseed cultivar ZhongShuang11 (ZS11) at seven developmental stages (seed germination, seedling, bolting, initial flowering, full-bloom, podding, and maturation). We determined the expression patterns of 101,040 genes via FPKM analysis and displayed the results using the eFP browser. We also analyzed transcriptome data for rapeseed from 70 BioProjects in the SRA database and obtained three types of expression level data (FPKM, TPM, and read counts). We used this information to develop the BrassicaEDB, including "eFP", "Treatment", "Coexpression", and "SRA Project" modules based on gene expression profiles and "Gene Feature", "qPCR Primer", and "BLAST" modules based on gene sequences. The BrassicaEDB provides comprehensive gene expression profile information and a user-friendly visualization interface for rapeseed researchers. Using this database, researchers can quickly retrieve the expression level data for target genes in different tissues and in response to different treatments to elucidate gene functions and explore the biology of rapeseed at the transcriptome level.

Genome-Wide Analysis of Phosphorus Transporter Genes in Brassica and Their Roles in Heavy Metal Stress Tolerance.

Phosphorus transporter (PHT) genes encode H2PO4-/H+ co-transporters that absorb and transport inorganic nutrient elements required for plant development and growth and protect plants from heavy metal stress. However, little is known about the roles of PHTs in Brassica compared to Arabidopsis thaliana. In this study, we identified and extensively analyzed 336 PHTs from three diploid (B. rapa, B. oleracea, and B. nigra) and two allotetraploid (B. juncea and B. napus) Brassica species. We categorized the PHTs into five phylogenetic clusters (PHT1-PHT5), including 201 PHT1 homologs, 15 PHT2 homologs, 40 PHT3 homologs, 54 PHT4 homologs, and 26 PHT5 homologs, which are unevenly distributed on the corresponding chromosomes of the five Brassica species. All PHT family genes from Brassica are more closely related to Arabidopsis PHTs in the same vs. other clusters, suggesting they are highly conserved and have similar functions. Duplication and synteny analysis revealed that segmental and tandem duplications led to the expansion of the PHT gene family during the process of polyploidization and that members of this family have undergone purifying selection during evolution based on Ka/Ks values. Finally, we explored the expression profiles of BnaPHT family genes in specific tissues, at various developmental stages, and under heavy metal stress via RNA-seq analysis and qRT-PCR. BnaPHTs that were induced by heavy metal treatment might mediate the response of rapeseed to this important stress. This study represents the first genome-wide analysis of PHT family genes in Brassica species. Our findings improve our understanding of PHT family genes and provide a basis for further studies of BnaPHTs in plant tolerance to heavy metal stress.

Genome-Wide Identification and Comparative Expression Profile Analysis of the Long-Chain Acyl-CoA synthetase (LACS) Gene Family in Two Different Oil Content Cultivars of Brassica napus.

Long-chain acyl-CoA synthetase (LACS) is one of the key enzymes involved in fatty acid metabolism, including phospholipid biosynthesis, triacylglycerol (TAG) biosynthesis, and fatty acid ß-oxidation in plants. However, the characterization of LACSs family in seed oil biosynthesis of Brassica napus (B. napus) remains unknown. In the present study, we performed a comprehensive genome-wide analysis of this gene family in B. napus, and 34 B. napus LACS genes (BnaLACSs) were identified. Phylogenetic analysis classified the BnaLACS proteins into four groups (A, B, C, and D), which were supported by highly conserved gene structures and consensus motifs. RNA-Sequencing (RNA-Seq) and qRT-PCR combined analysis revealed that 18 BnaLACSs (BnaLACS1-2, 1-3, 1-4, 1-9, 1-10, 2-1, 2-2, 4-1, 4-2, 6-1, 6-2, 6-4, 7-1, 7-2, 8-1, 8-2, 9-3, and 9-4) were highly expressed in developmental seeds. Comparative expression analysis between extremely high oil content (P1-HO) and low oil content (P2-LO) B. napus cultivars revealed that BnaLACS6-4, BnaLACS9-3, and BnaLACS9-4 may be involved in fatty acid synthesis in chloroplast, and BnaLACS1-10 and 4-1 may play a vital role in lipid biosynthesis in B. napus, which is important for further seed oil accumulation in oilseed rape. The present study provides important information for functional characterization of BnaLACSs in seed oil metabolism in B. napus.

Genome-Wide Characterization and Analysis of Metallothionein Family Genes That Function in Metal Stress Tolerance in Brassica napus L.

Brassica plants exhibit both high biomass productivity and high rates of heavy metal absorption. Metallothionein (MT) proteins are low molecular weight, cysteine-rich, metal-binding proteins that play crucial roles in protecting plants from heavy metal toxicity. However, to date, MT proteins have not been systematically characterized in Brassica. In this study, we identified 60 MTs from Arabidopsis thaliana and five Brassica species. All the MT family genes from Brassica are closely related to Arabidopsis MTs, encoding putative proteins that share similar functions within the same clades. Genome mapping analysis revealed high levels of synteny throughout the genome due to whole genome duplication and segmental duplication events. We analyzed the expression levels of 16 Brassica napus MTs (BnaMTs) by RNA-sequencing and real-time RT-PCR (RT-qPCR) analysis in plants under As3+ stress. These genes exhibited different expression patterns in various tissues. Our results suggest that BnaMT3C plays a key role in the response to As3+ stress in B. napus. This study provides insight into the phylogeny, origin, and evolution of MT family members in Brassica, laying the foundation for further studies of the roles of MT proteins in these important crops.

Genome-Wide Identification and Characterization of NODULE-INCEPTION-Like Protein (NLP) Family Genes in Brassica napus.

NODULE-INCEPTION-like proteins (NLPs) are conserved, plant-specific transcription factors that play crucial roles in responses to nitrogen deficiency. However, the evolutionary relationships and characteristics of NLP family genes in Brassica napus are unclear. In this study, we identified 31 NLP genes in B. napus, including 16 genes located in the A subgenome and 15 in the C subgenome. Subcellular localization predictions indicated that most BnaNLP proteins are localized to the nucleus. Phylogenetic analysis suggested that the NLP gene family could be divided into three groups and that at least three ancient copies of NLP genes existed in the ancestor of both monocots and dicots prior to their divergence. The ancestor of group III NLP genes may have experienced duplication more than once in the Brassicaceae species. Three-dimensional structural analysis suggested that 14 amino acids in BnaNLP7-1 protein are involved in DNA binding, whereas no binding sites were identified in the two RWP-RK and PB1 domains conserved in BnaNLP proteins. Expression profile analysis indicated that BnaNLP genes are expressed in most organs but tend to be highly expressed in a single organ. For example, BnaNLP6 subfamily members are primarily expressed in roots, while the four BnaNLP7 subfamily members are highly expressed in leaves. BnaNLP genes also showed different expression patterns in response to nitrogen-deficient conditions. Under nitrogen deficiency, all members of the BnaNLP1/4/5/9 subfamilies were upregulated, all BnaNLP2/6 subfamily members were downregulated, and BnaNLP7/8 subfamily members showed various expression patterns in different organs. These results provide a comprehensive evolutionary history of NLP genes in B. napus, and insight into the biological functions of BnaNLP genes in response to nitrogen deficiency.