Two aquaporins, PIP1;1 and PIP2;1, mediate the uptake of neonicotinoid pesticides in plants.

Neonicotinoids (NEOs), a large class of organic compounds, are a type of commonly used pesticide for crop protection. Their uptake and accumulation in plants are prerequisites for their intra- and intercellular movements, transformation, and function. Understanding the molecular mechanisms that underpin NEO uptake by plants is crucial for effective application, which remains elusive. Here, we demonstrate that NEOs enter plant cells primarily through the transmembrane symplastic pathway and accumulate mainly in the cytosol. Two plasma membrane intrinsic proteins discovered in Brassica rapa, BraPIP1;1 and BraPIP2;1, were found to encode aquaporins (AQPs) that are highly permeable to NEOs in different plant species and facilitate NEO subcellular diffusion and accumulation. Their conserved transport function was further demonstrated in Xenopus laevis oocyte and yeast assays. BraPIP1;1 and BraPIP2;1 gene knockouts and interaction assays suggested that their proteins can form functional heterotetramers. Assessment of the potential of mean force indicated a negative correlation between NEO uptake and the energy barrier of BraPIP1;1 channels. This study shows that AQPs transport organic compounds with greater osmolarity than previously thought, providing new insight into the molecular mechanisms of organic compound uptake and facilitating innovations in systemic pesticides.

Targeted mutagenesis of BnTTG1 homologues generated yellow-seeded rapeseed with increased oil content and seed germination under abiotic stress.

Yellow seed is one desirable trait with great potential to improve seed oil quality and yield. The present study surveys the redundant role of BnTTG1 genes in the proanthocyanidins (PA) biosynthesis, oil content and abiotic stress resistance. Stable yellow seed mutants were generated after mutating BnTTG1 by CRISPR/Cas9 genome editing system. Yellow seed phenotype could be obtained only when both functional homologues of BnTTG1 were simultaneously knocked out. Homozygous mutants of BnTTG1 homologues showed decreased thickness and PA accumulation in seed coat. Transcriptome and qRT-PCR analysis indicated that BnTTG1 mutation inhibited the expression of genes involved in phenylpropanoid and flavonoid biosynthetic pathways. Increased seed oil content and alteration of fatty acid (FA) composition were observed in homozygous mutants of BnTTG1 with enriched expression of genes involved in FA biosynthesis pathway. In addition, target mutation of BnTTG1 accelerated seed germination rate under salt and cold stresses. Enhanced seed germination capacity in BnTTG1 mutants was correlated with the change of expression level of ABA responsive genes. Overall, this study elucidated the redundant role of BnTTG1 in regulating seed coat color and established an efficient approach for generating yellow-seeded oilseed rape genetic resources with increase oil content, modified FA composition and resistance to multiple abiotic stresses.

Genomics of predictive radiation mutagenesis in oilseed rape: modifying seed oil composition.

Rapeseed is a crop of global importance but there is a need to broaden the genetic diversity available to address breeding objectives. Radiation mutagenesis, supported by genomics, has the potential to supersede genome editing for both gene knockout and copy number increase, but detailed knowledge of the molecular outcomes of radiation treatment is lacking. To address this, we produced a genome re-sequenced panel of 1133 M2 generation rapeseed plants and analysed large-scale deletions, single nucleotide variants and small insertion-deletion variants affecting gene open reading frames. We show that high radiation doses (2000 Gy) are tolerated, gamma radiation and fast neutron radiation have similar impacts and that segments deleted from the genomes of some plants are inherited as additional copies by their siblings, enabling gene dosage decrease. Of relevance for species with larger genomes, we showed that these large-scale impacts can also be detected using transcriptome re-sequencing. To test the utility of the approach for predictive alteration of oil fatty acid composition, we produced lines with both decreased and increased copy numbers of Bna.FAE1 and confirmed the anticipated impacts on erucic acid content. We detected and tested a 21-base deletion expected to abolish function of Bna.FAD2.A5, for which we confirmed the predicted reduction in seed oil polyunsaturated fatty acid content. Our improved understanding of the molecular effects of radiation mutagenesis will underpin genomics-led approaches to more efficient introduction of novel genetic variation into the breeding of this crop and provides an exemplar for the predictive improvement of other crops.

Two-Component System Genes in Brassica napus: Identification, Analysis, and Expression Patterns in Response to Abiotic and Biotic Stresses.

The two-component system (TCS), consisting of histidine kinases (HKs), histidine phosphotransfer proteins (HPs) and response regulators (RRs) in eukaryotes, plays pivotal roles in regulating plant growth, development, and responses to environment stimuli. However, the TCS genes were poorly characterized in rapeseed, which is an important tetraploid crop in Brassicaceae. In this work, a total of 182 BnaTCS genes were identified, including 43 HKs, 16 HPs, and 123 RRs, which was more than that in other crops due to segmental duplications during the process of polyploidization. It was significantly different in genetic diversity between the three subfamilies, and some members showed substantial genetic differentiation among the three rapeseed ecotypes. Several hormone- and stress-responsive cis-elements were identified in the putative promoter regions of BnaTCS genes. Furthermore, the expression of BnaTCS genes under abiotic stresses, exogenous phytohormone, and biotic stresses was analyzed, and numerous candidate stress-responsive genes were screened out. Meanwhile, using a natural population with 505 B. napus accessions, we explored the genetic effects of BnaTCS genes on salt tolerance by association mapping analysis and detected some significant association SNPs/genes. The result will help to further understand the functions of TCS genes in the developmental and stress tolerance improvement in B. napus.

Downregulation of the expression of subgenomic chromosome A7 genes promotes plant height in resynthesized allopolyploid Brassica napus.

KEY MESSAGE: Homoeolog expression bias and the gene dosage effect induce downregulation of genes on chromosome A7, causing a significant increase in the plant height of resynthesized allopolyploid Brassica napus. Gene expression levels in allopolyploid plants are not equivalent to the simple average of the expression levels in the parents and are associated with several non-additive expression phenomena, including homoeolog expression bias. However, hardly any information is available on the effect of homoeolog expression bias on traits. Here, we studied the effects of gene expression-related characteristics on agronomic traits using six isogenic resynthesized Brassica napus lines across the first ten generations. We found a group of genes located on chromosome A7 whose expression levels were significantly negatively correlated with plant height. They were expressed at significantly lower levels than their homoeologous genes, owing to allopolyploidy rather than inheritance from parents. Homoeolog expression bias resulted in resynthesized allopolyploids with a plant height similar to their female Brassica oleracea parent, but significantly higher than that of the male Brassica rapa parent. Notably, aneuploid lines carrying monosomic and trisomic chromosome A7 had the highest and lowest plant heights, respectively, due to changes in the expression bias of homoeologous genes because of alterations in the gene dosage. These findings suggest that the downregulation of the expression of homoeologous genes on a single chromosome can result in the partial improvement of traits to a significant extent in the nascent allopolyploid B. napus.

Transcriptomic Analysis of the Reduction in Seed Oil Content through Increased Nitrogen Application Rate in Rapeseed (Brassica napus L.).

Nitrogen is essential for improving the seed oil yield of rapeseed (Brassica napus L.). However, the molecular mechanism by which increased nitrogen rates impact seed oil content is largely unknown. Therefore, a field experiment was conducted to determine how three nitrogen application rates (120, 240, and 360 kg ha-1) regulated seed oil content via transcriptomic analysis. The results showed that the seed yield and the protein and total N contents increased from N1 to N3, with average increases of 57.2%, 16.9%, and 79.5%, respectively. However, the seed oil content significantly decreased from N1 to N3, with an average decrease of 8.6%. These results were repeated over a number of years. The quantity of oil protein bodies observed under a transmission electron microscope was in accordance with the ultimate seed oil and protein contents. As the nitrogen application rate increased, a substantial number of genes involved in the photosynthesis, glycolysis, and phenylpropanoid biosynthesis pathways were up-regulated, as were TF families, such as AP2/ERF, MYB, and NAC. The newly identified genes were mainly involved in carbohydrate, lipid, and amino acid metabolism. Metabolic flux analysis showed that most of the genes involved in glycolysis and fatty acid biosynthesis had higher transcript levels in the early development stages. Our results provide new insights into the molecular regulation of rapeseed seed oil content through increased nitrogen application rates.

The lack of negative association between TE load and subgenome dominance in synthesized Brassica allotetraploids.

Polyploidization is important to the evolution of plants. Subgenome dominance is a distinct phenomenon associated with most allopolyploids. A gene on the dominant subgenome tends to express to higher RNA levels in all organs as compared to the expression of its syntenic paralogue (homoeolog). The mechanism that underlies the formation of subgenome dominance remains unknown, but there is evidence for the involvement of transposon/DNA methylation density differences nearby the genes of parents as being causal. The subgenome with lower density of transposon and methylation near genes is positively associated with subgenome dominance. Here, we generated eight generations of allotetraploid progenies from the merging of parental genomes Brassica rapa and Brassica oleracea. We found that transposon/methylation density differ near genes between the parental (rapa:oleracea) existed in the wide hybrid, persisted in the neotetraploids (the synthetic Brassica napus), but these neotetraploids expressed no expected subgenome dominance. This absence of B. rapa vs. B. oleracea subgenome dominance is particularly significant because, while there is no negative relationship between transposon/methylation level and subgenome dominance in the neotetraploids, the more ancient parental subgenomes for all Brassica did show differences in transposon/methylation densities near genes and did express, in the same samples of cells, biased gene expression diagnostic of subgenome dominance. We conclude that subgenome differences in methylated transposon near genes are not sufficient to initiate the biased gene expressions defining subgenome dominance. Our result was unexpected, and we suggest a "nuclear chimera" model to explain our data.

A novel type of Brassica napus with higher stearic acid in seeds developed through genome editing of BnaSAD2 family.

KEY MESSAGE: Modifications of multiple copies of the BnaSAD2 gene family with genomic editing technology result in higher stearic acid content in the seed of polyploidy rapeseed. Solid fats from vegetable oils are widely used in food processing industry. Accumulating data showed that stearic acid is more favorite as the major composite among the saturate fatty acids in solid fats in considerations of its effects on human health. Rapeseed is the third largest oil crop worldwide, and has potential to be manipulated to produce higher saturated fatty acids as raw materials of solid fats. Toward that end, we identified four SAD2 gene family members in B. napus genome and established spatiotemporal expression pattern of the BnaSAD2 members. Genomic editing technology was applied to mutate all the copies of BnaSAD2 in this allopolyploid species and mutants at multiple alleles were generated and characterized to understand the effect of each BnaSAD2 member on blocking desaturation of stearic acid. Mutations occurred at BnaSAD2.A3 resulted in more dramatic changes of fatty acid profile than ones on BnaSAD2.C3, BnaSAD2.A5 and BnaSAD2.C4. The content of stearic acid in mutant seeds with single locus increased dramatically with a range of 3.1-8.2%. Furthermore, combination of different mutated alleles of BnaSAD2 resulted in more dramatic changes in fatty acid profiles and the double mutant at BnaSAD2.A3 and BnaSAD2.C3 showed the most dramatic phenotypic changes compared with its single mutants and other double mutants, leading to 11.1% of stearic acid in the seeds. Our results demonstrated that the members of BnaSAD2 have differentiated in their efficacy as a Δ9-Stearoyl-ACP-Desaturase and provided valuable rapeseed germplasm for breeding high stearic rapeseed oil.

Combining transcriptomics and metabolomics to identify key response genes for aluminum toxicity in the root system of Brassica napus L. seedlings.

KEY MESSAGE: By integrating QTL mapping, transcriptomics and metabolomics, 138 hub genes were identified in rapeseed root response to aluminum stress and mainly involved in metabolism of lipids, carbohydrates and secondary metabolites. Aluminum (Al) toxicity has become one of the important abiotic stress factors in areas with acid soil, which hinders the absorption of water and nutrients by roots, and consequently retards the growth of crops. A deeper understanding of the stress-response mechanism of Brassica napus may allow us to identify the tolerance gene(s) and use this information in breeding-resistant crop varieties. In this study, a population of 138 recombinant inbred lines (RILs) was subjected to aluminum stress, and QTL (quantitative trait locus) mapping was used to preliminarily locate quantitative trait loci related to aluminum stress. Root tissues from seedlings of an aluminum-resistant (R) line and an aluminum-sensitive (S) line from the RIL population were harvested for transcriptome sequencing and metabolome determination. By combining the data on quantitative trait genes (QTGs), differentially expressed genes (DEGs), and differentially accumulated metabolites (DAMs), key candidate genes related to aluminum tolerance in rapeseed were determined. The results showed that there were 3186 QTGs in the RIL population, 14,232 DEGs and 457 DAMs in the comparison between R and S lines. Lastly, 138 hub genes were selected to have a strong positive or negative correlation with 30 important metabolites (|R|≥ 0.95). These genes were mainly involved in the metabolism of lipids, carbohydrates and secondary metabolites in response to Al toxicity stress. In summary, this study provides an effective method for screening key genes by combining QTLs, transcriptome sequencing and metabolomic analysis, but also lists key genes for exploring the molecular mechanism of Al tolerance in rapeseed seedling roots.

Genetic factors inherited from both diploid parents interact to affect genome stability and fertility in resynthesized allotetraploid Brassica napus.

Established allopolyploids are known to be genomically stable and fertile. However, in contrast, most newly resynthesized allopolyploids are infertile and meiotically unstable. Identifying the genetic factors responsible for genome stability in newly formed allopolyploid is key to understanding how 2 genomes come together to form a species. One hypothesis is that established allopolyploids may have inherited specific alleles from their diploid progenitors which conferred meiotic stability. Resynthesized Brassica napus lines are often unstable and infertile, unlike B. napus cultivars. We tested this hypothesis by characterizing 41 resynthesized B. napus lines produced by crosses between 8 Brassica rapa and 8 Brassica oleracea lines for copy number variation resulting from nonhomologous recombination events and fertility. We resequenced 8 B. rapa and 5 B. oleracea parent accessions and analyzed 19 resynthesized lines for allelic variation in a list of meiosis gene homologs. SNP genotyping was performed using the Illumina Infinium Brassica 60K array for 3 individuals per line. Self-pollinated seed set and genome stability (number of copy number variants) were significantly affected by the interaction between both B. rapa and B. oleracea parental genotypes. We identified 13 putative meiosis gene candidates which were significantly associated with frequency of copy number variants and which contained putatively harmful mutations in meiosis gene haplotypes for further investigation. Our results support the hypothesis that allelic variants inherited from parental genotypes affect genome stability and fertility in resynthesized rapeseed.

Engineering the Staple Oil Crop Brassica napus Enriched with α-Linolenic Acid Using the Perilla FAD2-FAD3 Fusion Gene.

Modern people generally suffer from α-linolenic acid (ALA) deficiency, since most staple food oils are low in ALA content. Thus, the enhancement of ALA in staple oil crops is of importance. In this study, the FAD2 and FAD3 coding regions from the ALA-king species Perilla frutescens were fused using a newly designed double linker LP4-2A, driven by a seed-specific promoter PNAP, and engineered into a rapeseed elite cultivar ZS10 with canola quality background. The mean ALA content in the seed oil of PNAP:PfFAD2-PfFAD3 (N23) T5 lines was 3.34-fold that of the control (32.08 vs 9.59%), with the best line being up to 37.47%. There are no significant side effects of the engineered constructs on the background traits including oil content. In fatty acid biosynthesis pathways, the expression levels of structural genes as well as regulatory genes were significantly upregulated in N23 lines. On the other hand, the expression levels of genes encoding the positive regulators of flavonoid-proanthocyanidin biosynthesis but negative regulators of oil accumulation were significantly downregulated. Surprisingly, the ALA level in PfFAD2-PfFAD3 transgenic rapeseed lines driven by the constitutive promoter PD35S was not increased or even showed a slight decrease due to the lower level of foreign gene expression and downregulation of the endogenous orthologous genes BnFAD2 and BnFAD3.

Transcriptome analysis reveals different mechanisms of selenite and selenate regulation of cadmium translocation in Brassica rapa.

Selenium (Se) inhibits cadmium (Cd) root-to-shoot translocation and accumulation in the shoots of pak choi; however, the mechanism by which Se regulates Cd retention in roots is still poorly understood. A time-dependent hydroponic experiment was conducted to compare the effects of selenite and selenate on Cd translocation and retention in the roots. The underlying mechanisms were investigated regarding Se biotransformation and metal transportation in roots using HPLC and transcriptome analyses. Selenite showed reducing effects on Cd translocation and accumulation in shoots earlier than selenate. Selenite is mainly biotransformed into selenomethionine (80% of total Se in roots) at 72 h, while SeO42- was the dominant species in the selenate treatments (68% in shoots). Selenite up-regulated genes involved in the biosynthesis of lignin, suberin, and phytochelatins and those involved in stress signaling, thereby helping to retain Cd in the roots, whereas essentially, selenate had opposite effects and impaired the symplastic and apoplastic retention of Cd. These results suggest that cell-wall reinforcement and Cd retention in roots may be the key processes by which Se regulates Cd accumulation, and faster biotransformation into organic seleno-compounds could lead to earlier effects.

Subgenome-dominant expression and alternative splicing in response to Sclerotinia infection in polyploid Brassica napus and progenitors.

Polyploidy has played an extensive role in the evolution of flowering plants. Allopolyploids, with subgenomes containing duplicated gene pairs called homeologs, can show rapid transcriptome changes including novel alternative splicing (AS) patterns. The extent to which abiotic stress modulates AS of homeologs is a nascent topic in polyploidy research. We subjected both resynthesized and natural lines of polyploid Brassica napus, along with the progenitors Brassica rapa and Brassica oleracea, to infection with the fungal pathogen Sclerotinia sclerotiorum. RNA-sequencing analyses revealed widespread divergence between polyploid subgenomes in both gene expression and AS patterns. Resynthesized B. napus displayed significantly more A and C subgenome biased homeologs under pathogen infection than during uninfected growth. Differential AS (DAS) in response to infection was highest in natural B. napus (12 709 DAS events) and lower in resynthesized B. napus (8863 DAS events). Natural B. napus had more upregulated events and fewer downregulated events. There was a global expression bias towards the B. oleracea-derived (C) subgenome in both resynthesized and natural B. napus, enhanced by widespread non-parental downregulation of the B. rapa-derived (A) homeolog. In the resynthesized B. napus, this resulted in a disproportionate C subgenome contribution to the pathogen defense response, characterized by biases in both transcript expression levels and the proportion of induced genes. Our results elucidate the complex ways in which Sclerotinia infection affects expression and AS of homeologous genes in resynthesized and natural B. napus.

Genome-wide survey of catalase genes in Brassica rapa, Brassica oleracea, and Brassica napus: identification, characterization, molecular evolution, and expression profiling of BnCATs in response to salt and cadmium stress.

Catalase (CAT, EC 1.11.1.6), one of the most important antioxidant enzymes, can control excess levels of H2O2 produced under oxidative stress in plants. In this study, 16, 8, and 7 CAT genes in the genome of Brassica napus, B. rapa, and B. oleracea were identified, respectively. Phylogenetic studies showed that CATs could be divided into two main groups, each containing specific monocotyledon and dicotyledon subgroups. Motifs, gene structure, and intron phase of CATs in B. napus, Brassica rapa, and Brassica oleracea are highly conserved. Analysis of codon usage bias showed the mutation pressure and natural selection of the codon usage of CATs. Segmental duplication and polyploid were major factors in the expansion of this gene family in B. napus, and genes have experienced  negative selection during evolution. Existence of hormones and stress-responsive cis-elements and identifying miRNA molecules affecting CATs showed that these genes are complexly regulated at the transcriptional and posttranscriptional levels. Based on RNA-seq data, CATs are divided into two groups; the first group has moderate and specific expression in flowers, leaves, stems, and roots, while the second group shows expression in most tissues. qRT-PCR analysis showed that the expression of these genes is dynamic and has a specific expression consistent with other CAT genes in response to salinity and cadmium (Cd) stresses. These results provide information for further investigation of the function of CAT genes in response to stresses and the development of tolerant  plants.

Kinase CIPK9 integrates glucose and abscisic acid signaling to regulate seed oil metabolism in rapeseed.

Rapeseed (Brassica napus), an important oil crop worldwide, provides large amounts of lipids for human requirements. Calcineurin B-like (CBL)-interacting protein kinase 9 (CIPK9) was reported to regulate seed oil content in the plant. Here, we generated gene-silenced lines through RNA interference biotechnology and loss-of-function mutant bnacipk9 using CRISPR/Cas9 to further study BnaCIPK9 functions in the seed oil metabolism of rapeseeds. We discovered that compared with wild-type (WT) lines, gene-silenced and bnacipk9 lines had substantially different oil contents and fatty acid compositions: seed oil content was improved by 3%-5% and 1%-6% in bnacipk9 lines and gene-silenced lines, respectively; both lines were with increased levels of monounsaturated fatty acids and decreased levels of polyunsaturated fatty acids. Additionally, hormone and glucose content analyses revealed that compared with WT lines the bnacipk9 lines showed significant differences: in bnacipk9 seeds, indoleacetic acid and abscisic acid (ABA) levels were higher; glucose and sucrose contents were higher with a higher hexose-to-sucrose ratio in bnacipk9 mid-to-late maturation development seeds. Furthermore, the bnacipk9 was less sensitive to glucose and ABA than the WT according to stomatal aperture regulation assays and the expression levels of genes involved in glucose and ABA regulating pathways in rapeseeds. Notably, in Arabidopsis (Arabidopsis thaliana), exogenous ABA and glucose imposed on developing seeds revealed the effects of ABA and glucose signaling on seed oil accumulation. Altogether, our results strongly suggest a role of CIPK9 in mediating the interaction between glucose flux and ABA hormone signaling to regulate seed oil metabolism in rapeseed.

Differences in pseudogene evolution contributed to the contrasting flavors of turnip and Chiifu, two Brassica rapa subspecies.

Pseudogenes are important resources for investigation of genome evolution and genomic diversity because they are nonfunctional but have regulatory effects that influence plant adaptation and diversification. However, few systematic comparative analyses of pseudogenes in closely related species have been conducted. Here, we present a turnip (Brassica rapa ssp. rapa) genome sequence and characterize pseudogenes among diploid Brassica species/subspecies. The results revealed that the number of pseudogenes was greatest in Brassica oleracea (CC genome), followed by B. rapa (AA genome) and then Brassica nigra (BB genome), implying that pseudogene differences emerged after species differentiation. In Brassica AA genomes, pseudogenes were distributed asymmetrically on chromosomes because of numerous chromosomal insertions/rearrangements, which contributed to the diversity among subspecies. Pseudogene differences among subspecies were reflected in the flavor-related glucosinolate (GSL) pathway. Specifically, turnip had the highest content of pungent substances, probably because of expansion of the methylthioalkylmalate synthase-encoding gene family in turnips; these genes were converted into pseudogenes in B. rapa ssp. pekinensis (Chiifu). RNA interference-based silencing of the gene encoding 2-oxoglutarate-dependent dioxygenase 2, which is also associated with flavor and anticancer substances in the GSL pathway, resulted in increased abundance of anticancer compounds and decreased pungency of turnip and Chiifu. These findings revealed that pseudogene differences between turnip and Chiifu influenced the evolution of flavor-associated GSL metabolism-related genes, ultimately resulting in the different flavors of turnip and Chiifu.

Glutamic Acid and Poly-γ-glutamic Acid Enhanced the Heat Resistance of Chinese Cabbage (Brassica rapa L. ssp. pekinensis) by Improving Carotenoid Biosynthesis, Photosynthesis, and ROS Signaling.

Heat stress is one of the most common agrometeorological risks in crop production in the middle and lower reaches of the Yangtze River in China. This study aimed to investigate whether glutamic acid (Glu) or poly-γ-glutamic acid (γ-PGA) biostimulants can improve the thermotolerance of a cool-season Chinese cabbage (Brassica rapa L. ssp. pekinensis) crop. Priming with Glu (2.0 mM) or γ-PGA (20 mg·L-1) was conducted at the third leaf stage by applying as daily foliar sprays for 5 days before 5 days of heat stress (45 °C in 16-h light/35 °C in 8-h dark). Coupled with morpho-physiological and biochemical analyses, transcriptomes of Glu or γ-PGA-primed Chinese cabbage under heat stress were examined by RNA-seq analysis. The results showed that the thermotolerance conferred by Glu and γ-PGA priming was associated with the increased parameters of vegetative growth, gas exchange, and chlorophyll fluorescence. Compared with the control, the dry weights of plants treated with Glu and γ-PGA increased by 51.52% and 39.39%, respectively. Glu and γ-PGA application also significantly increased the contents of total chlorophyll by 42.21% and 23.12%, and carotenoid by 32.00% and 24.00%, respectively. In addition, Glu- and γ-PGA-primed plants markedly inhibited the levels of malondialdehyde, electrolyte leakage, and super-oxide anion radical, which was accompanied by enhanced activity levels of superoxide dismutase (SOD), catalase (CAT), ascorbate peroxidase (APX), and peroxidase (POD). Enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) categories within the differentially expressed genes (DEGs) functional clusters of RNA-seq data indicated that the expression levels of the genes for DNA replication, DNA repair system, linoleic acid metabolism, cysteine and methionine metabolism, glutathione metabolism, purine and pyrimidine metabolism, carotenoid biosynthesis, and plant-pathogen interaction were commonly up-regulated by both Glu and γ-PGA priming. Glu treatment enhanced the expression levels of the genes involved in aliphatic glucosinolate and 2-oxocarboxylic acid, while γ-PGA treatment activated carotenoid cleavage reaction to synthesize abscisic acid. Taken together, both Glu and γ-PGA have great potential for the preadaptation of Chinese cabbage seedlings to heat stress, with Glu being more effective than γ-PGA.

Genome-Wide Identification and Expression Analysis of nsLTP Gene Family in Rapeseed (Brassica napus) Reveals Their Critical Roles in Biotic and Abiotic Stress Responses.

Non-specific lipid transfer proteins (nsLTPs) are small cysteine-rich basic proteins which play essential roles in plant growth, development and abiotic/biotic stress response. However, there is limited information about the nsLTP gene (BnLTP) family in rapeseed (Brassica napus). In this study, 283 BnLTP genes were identified in rapeseed, which were distributed randomly in 19 chromosomes of rapeseed. Phylogenetic analysis showed that BnLTP proteins were divided into seven groups. Exon/intron structure and MEME motifs both remained highly conserved in each BnLTP group. Segmental duplication and hybridization of rapeseed's two sub-genomes mainly contributed to the expansion of the BnLTP gene family. Various potential cis-elements that respond to plant growth, development, biotic/abiotic stresses, and phytohormone signals existed in BnLTP gene promoters. Transcriptome analysis showed that BnLTP genes were expressed in various tissues/organs with different levels and were also involved in the response to heat, drought, NaCl, cold, IAA and ABA stresses, as well as the treatment of fungal pathogens (Sclerotinia sclerotiorum and Leptosphaeria maculans). The qRT-PCR assay validated the results of RNA-seq expression analysis of two top Sclerotinia-responsive BnLTP genes, BnLTP129 and BnLTP161. Moreover, batches of BnLTPs might be regulated by BnTT1 and BnbZIP67 to play roles in the development, metabolism or adaptability of the seed coat and embryo in rapeseed. This work provides an important basis for further functional study of the BnLTP genes in rapeseed quality improvement and stress resistance.

Genome-Wide Analysis of the Type-B Authentic Response Regulator Gene Family in Brassica napus.

The type-B authentic response regulators (type-B ARRs) are positive regulators of cytokinin signaling and involved in plant growth and stress responses. In this study, we used bioinformatics, RNA-seq, and qPCR to study the phylogenetic and expression pattern of 35 type-B ARRs in Brassica napus. The BnARRs experienced gene expansion and loss during genome polyploidization and were classified into seven groups. Whole-genome duplication (WGD) and segmental duplication were the main forces driving type-B ARR expansion in B. napus. Several BnARRs with specific expression patterns during rapeseed development were identified, including BnARR12/14/18/23/33. Moreover, we found the type-B BnARRs were involved in rapeseed development and stress responses, through participating in cytokinin and ABA signaling pathways. This study revealed the origin, evolutionary history, and expression pattern of type-B ARRs in B. napus and will be helpful to the functional characterization of BnARRs.

Transcriptomics Integrated with Metabolomics Unveil Carotenoids Accumulation and Correlated Gene Regulation in White and Yellow-Fleshed Turnip (Brassica rapa ssp. rapa).

Turnip (Brassica rapa ssp. rapa) is considered to be a highly nutritious and health-promoting vegetable crop, whose flesh color can be divided into yellow and white. It is widely accepted that yellow-fleshed turnips have higher nutritional value. However, reports about flesh color formation is lacking. Here, the white-fleshed inbred line, W21, and yellow-fleshed inbred line, W25, were profiled from the swollen root of the turnip at three developmental periods to elucidate the yellow color formation. Transcriptomics integrated with metabolomics analysis showed that the PSY gene was the key gene affecting the carotenoids formation in W25. The coding sequence of BrrPSY-W25 was 1278 bp and that of BrrPSY-W21 was 1275 bp, and BrrPSY was more highly expressed in swollen roots in W25 than in W21. Transient transgenic tobacco leaf over-expressing BrrPSY-W and BrrPSY-Y showed higher transcript levels and carotenoids contents. Results revealed that yellow turnip formation is due to high expression of the PSY gene rather than mutations in the PSY gene, indicating that a post-transcriptional regulatory mechanism may affect carotenoids formation. Results obtained in this study will be helpful for explaining the carotenoids accumulation of turnips.