Genome-Wide Identification and Characterization of the CCT Gene Family in Rapeseed (Brassica napus L.).

The CCT gene family is present in plants and is involved in biological processes such as flowering, circadian rhythm regulation, plant growth and development, and stress resistance. We identified 87, 62, 46, and 40 CCTs at the whole-genome level in B. napus, B. rapa, B. oleracea, and A. thaliana, respectively. The CCTs can be classified into five groups based on evolutionary relationships, and each of these groups can be further subdivided into three subfamilies (COL, CMF, and PRR) based on function. Our analysis of chromosome localization, gene structure, collinearity, cis-acting elements, and expression patterns in B. napus revealed that the distribution of the 87 BnaCCTs on the chromosomes of B. napus was uneven. Analysis of gene structure and conserved motifs revealed that, with the exception of a few genes that may have lost structural domains, the majority of genes within the same group exhibited similar structures and conserved domains. The gene collinearity analysis identified 72 orthologous genes, indicating gene duplication and expansion during the evolution of BnaCCTs. Analysis of cis-acting elements identified several elements related to abiotic and biotic stress, plant hormone response, and plant growth and development in the promoter regions of BnaCCTs. Expression pattern and protein interaction network analysis showed that BnaCCTs are differentially expressed in various tissues and under stress conditions. The PRR subfamily genes have the highest number of interacting proteins, indicating their significant role in the growth, development, and response to abiotic stress of B. napus.

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.

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.

Joint QTL mapping and transcriptome sequencing analysis reveal candidate flowering time genes in Brassica napus L.

BACKGROUND: Optimum flowering time is a key agronomic trait in Brassica napus. To investigate the genetic architecture and genetic regulation of flowering time in this important crop, we conducted quantitative trait loci (QTL) analysis of flowering time in a recombinant inbred line (RIL) population, including lines with extreme differences in flowering time, in six environments, along with RNA-Seq analysis. RESULTS: We detected 27 QTLs distributed on eight chromosomes among six environments, including one major QTL on chromosome C02 that explained 11-25% of the phenotypic variation and was stably detected in all six environments. RNA-Seq analysis revealed 105 flowering time-related differentially expressed genes (DEGs) that play roles in the circadian clock/photoperiod, autonomous pathway, and hormone and vernalization pathways. We focused on DEGs related to the regulation of flowering time, especially DEGs in QTL regions. CONCLUSIONS: We identified 45 flowering time-related genes in these QTL regions, eight of which are DEGs, including key flowering time genes PSEUDO RESPONSE REGULATOR 7 (PRR7) and FY (located in a major QTL region on C02). These findings provide insights into the genetic architecture of flowering time in B. napus.

Genome-wide exploration and characterization of miR172/euAP2 genes in Brassica napus L. for likely role in flower organ development.

BACKGROUND: APETALA2-like genes encode plant-specific transcription factors, some of which possess one microRNA172 (miR172) binding site. The miR172 and its target euAP2 genes are involved in the process of phase transformation and flower organ development in many plants. However, the roles of miR172 and its target AP2 genes remain largely unknown in Brassica napus (B. napus). RESULTS: In this study, 19 euAP2 and four miR172 genes were identified in the B. napus genome. A sequence analysis suggested that 17 euAP2 genes were targeted by Bna-miR172 in the 3' coding region. EuAP2s were classified into five major groups in B.napus. This classification was consistent with the exon-intron structure and motif organization. An analysis of the nonsynonymous and synonymous substitution rates revealed that the euAP2 genes had gone through purifying selection. Whole genome duplication (WGD) or segmental duplication events played a major role in the expansion of the euAP2 gene family. A cis-regulatory element (CRE) analysis suggested that the euAP2s were involved in the response to light, hormones, stress, and developmental processes including circadian control, endosperm and meristem expression. Expression analysis of the miR172-targeted euAP2s in nine different tissues showed diverse spatiotemporal expression patterns. Most euAP2 genes were highly expressed in the floral organs, suggesting their specific functions in flower development. BnaAP2-1, BnaAP2-5 and BnaTOE1-2 had higher expression levels in late-flowering material than early-flowering material based on RNA-seq and qRT-PCR, indicating that they may act as floral suppressors. CONCLUSIONS: Overall, analyses of the evolution, structure, tissue specificity and expression of the euAP2 genes were peformed in B.napus. Based on the RNA-seq and experimental data, euAP2 may be involved in flower development. Three euAP2 genes (BnaAP2-1, BnaAP2-5 and BnaTOE1-2) might be regarded as floral suppressors. The results of this study provide insights for further functional characterization of the miR172 /euAP2 module in B.napus.

Genome-wide identification and characterization of Gretchen Hagen3 (GH3) family genes in Brassica napus.

The hormone auxin is involved in many biological processes throughout a plant's lifecycle. However, genes in the GH3 (Gretchen Hagen3) family, one of the three major auxin-responsive gene families, have not yet been identified in oilseed rape (Brassica napus). In this study, we identified 63 BnaGH3 genes in oilseed rape using homology searches. We analyzed the chromosome locations, gene structures, and phylogenetic relationships of the BnaGH3 genes, as well as the cis-elements in their promoters. Most BnaGH3 genes are located on chromosomes A03, A09, C02, C03, and C09, each with 4-7 members. In addition, we analyzed the expression patterns of BnaGH3 genes in seven tissues by transcriptome sequencing and quantitative RT-PCR analysis of plants under exogenous IAA treatment. The BnaGH3 genes showed different expression patterns in various tissues. BnaA.GH3.2-1 and BnaC.GH3.2-1 were expressed in the seed and seed coat during development and in response to IAA treatment. These results shed light on the possible roles of the GH3 gene family in oilseed rape.

Genome-Wide Identification of the LAC Gene Family and Its Expression Analysis Under Stress in Brassica napus.

Lignin is an important biological polymer in plants that is necessary for plant secondary cell wall ontogenesis. The laccase (LAC) gene family catalyzes lignification and has been suggested to play a vital role in the plant kingdom. In this study, we identified 45 LAC genes from the Brassica napus genome (BnLACs), 25 LAC genes from the Brassica rapa genome (BrLACs) and 8 LAC genes from the Brassica oleracea genome (BoLACs). These LAC genes could be divided into five groups in a cladogram and members in same group had similar structures and conserved motifs. All BnLACs contained hormone- and stress- related elements determined by cis-element analysis. The expression of BnLACs was relatively higher in the root, seed coat and stem than in other tissues. Furthermore, BnLAC4 and its predicted downstream genes showed earlier expression in the silique pericarps of short silique lines than long silique lines. Three miRNAs (miR397a, miR397b and miR6034) target 11 BnLACs were also predicted. The expression changes of BnLACs under series of stresses were further investigated by RNA sequencing (RNA-seq) and quantitative real-time polymerase chain reaction (qRT-PCR). The study will give a deeper understanding of the LAC gene family evolution and functions in B. napus.

Genome-Wide Identification of MicroRNAs in Response to Cadmium Stress in Oilseed Rape (Brassica napus L.) Using High-Throughput Sequencing.

MicroRNAs (miRNAs) have important roles in regulating stress-response genes in plants. However, identification of miRNAs and the corresponding target genes that are induced in response to cadmium (Cd) stress in Brassica napus remains limited. In the current study, we sequenced three small-RNA libraries from B. napus after 0 days, 1 days, and 3 days of Cd treatment. In total, 44 known miRNAs (belonging to 27 families) and 103 novel miRNAs were identified. A comprehensive analysis of miRNA expression profiles found 39 differentially expressed miRNAs between control and Cd-treated plants; 13 differentially expressed miRNAs were confirmed by qRT-PCR. Characterization of the corresponding target genes indicated functions in processes including transcription factor regulation, biotic stress response, ion homeostasis, and secondary metabolism. Furthermore, we propose a hypothetical model of the Cd-response mechanism in B. napus. Combined with qRT-PCR confirmation, our data suggested that miRNAs were involved in the regulations of TFs, biotic stress defense, ion homeostasis and secondary metabolism synthesis to respond Cd stress 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 association mapping and Identification of candidate genes for fatty acid composition in Brassica napus L. using SNP markers.

BACKGROUND: B. napus (oilseed) is an important source of edible vegetable oil, and its nutritional and economic value is determined by its fatty acid composition and content. RESULTS: Using the Brassica 60 K SNP array, we performed a genome-wide association study of fatty acid composition in a population of 520 genetically diverse oilseed accessions. Using the PCA + K model in TASSEL 5.2.1, we identified 62 genomic regions that were significantly associated with the composition of seven fatty acids, and five consensus regions that mapped to the A2, A8, A9, C1, and C3 chromosomes, respectively, of the Brassica napus Darmor-bzh genome. We then identified 24 orthologs of the functional candidate genes involved in fatty acid biosynthesis, excluding BnaA.FAE1 and BnaC.FAE1 on the A8 and C3 homologous genome blocks, which are known to have critical roles in the fatty acid biosynthesis pathway, and potential orthologs of these genes (e.g., LACS9, KCR1, FAB1, LPAT4, KCS17, CER4, TT16, and ACBP5). CONCLUSIONS: Our results demonstrate the power of association mapping in identifying genes of interest in B. napus and provide insight into the genetic basis of fatty acid biosynthesis in B. napus. Furthermore, our findings may facilitate marker-based breeding efforts aimed at improving fatty acid composition and quality in B. napus.

Genetic and transcriptomic analyses of lignin- and lodging-related traits in Brassica napus.

KEY MESSAGE: Candidate genes associated with lignin and lodging traits were identified by combining phenotypic, genotypic, and gene expression data in B. napus. Brassica napus is one of the world's most important oilseed crops, but its yield can be dramatically reduced by lodging, bending, and falling of its vertical stems. Lignin has been shown to contribute to stem mechanical strength. In this study, we found that the syringyl/guaiacyl (S/G) monolignol ratio exhibits a significant negative correlation with disease and lodging resistance. A total of 92 and 50 SNP and SSR loci, respectively, were found to be significantly associated with five traits, breaking force, breaking strength, lodging coefficient, acid detergent lignin content, and the S/G monolignol ratio using GWAS. To identify novel genes involved in lignin biosynthesis, transcriptome sequencing of high- (H) and low (L)-ADL content accessions was performed. The up-regulated genes were mainly involved in glycoside catabolic processes (especially glucosinolate catabolism) and cell wall biogenesis, while down-regulated genes were involved in glucosinolate biosynthesis, indicating that crosstalk exists between glucosinolate metabolic processes and lignin biosynthesis. Integrating this differential expression with the GWAS analysis, we identified four candidate genes regulating lignin, including glycosyl hydrolase (BnaA01g00480D), CYT1 (BnaA04g22820D), and two encoding transcription factors, SHINE1 (ERF family) and DAR6 (LIM family). This study provides insight into the genetic control of lodging and lignin in B. napus.

Development of genic cleavage markers in association with seed glucosinolate content in canola.

KEY MESSAGE: The orthologues of Arabidopsis involved in seed glucosinolates metabolism within QTL confidence intervals were identified, and functional markers were developed to facilitate breeding for ultra-low glucosinolates in canola. Further reducing the content of seed glucosinolates will have a positive impact on the seed quality of canola (Brassica napus). In this study 43 quantitative trait loci (QTL) for seed glucosinolate (GSL) content in a low-GSL genetic background were mapped over seven environments in Germany and China in a doubled haploid population from a cross between two low-GSL oilseed rape parents with transgressive segregation. By anchoring these QTL to the reference genomes of B. rapa and B. oleracea, we identified 23 orthologues of Arabidopsis involved in GSL metabolism within the QTL confidence intervals. Sequence polymorphisms between the corresponding coding regions of the parental lines were used to develop cleaved amplified polymorphic site markers for two QTL-linked genes, ISOPROPYLMALATE DEHYDROGENASE1 and ADENOSINE 5'-PHOSPHOSULFATE REDUCTASE 3. The genic cleavage markers were mapped in the DH population into the corresponding intervals of QTL explaining 3.36-6.88 and 4.55-8.67 % of the phenotypic variation for seed GSL, respectively. The markers will facilitate breeding for ultra-low seed GSL content in canola.

Differential accumulation of phenolic compounds and expression of related genes in black- and yellow-seeded Brassica napus.

Developing yellow-seeded Brassica napus (rapeseed) with improved qualities is a major breeding goal. The intermediate and final metabolites of the phenylpropanoid and flavonoid pathways affect not only oil quality but also seed coat colour of B. napus. Here, the accumulation of phenolic compounds was analysed in the seed coats of black-seeded (ZY821) and yellow-seeded (GH06) B. napus. Using toluidine blue O staining and liquid chromatography-mass spectrometry, histochemical and biochemical differences were identified in the accumulation of phenolic compounds between ZY821 and GH06. Two and 13 unique flavonol derivatives were detected in ZY821 and GH06, respectively. Quantitative real-time PCR analysis revealed significant differences between ZY821 and GH06 in the expression of common phenylpropanoid biosynthetic genes (BnPAL and BnC4H), common flavonoid biosynthetic genes (BnTT4 and BnTT6), anthocyanin- and proanthocyandin-specific genes (BnTT3 and BnTT18), proanthocyandin-specific genes (BnTT12, BnTT10, and BnUGT2) and three transcription factor genes (BnTTG1, BnTTG2, and BnTT8) that function in the flavonoid biosynthetic pathway. These data provide insight into pigment accumulation in B. napus, and serve as a useful resource for researchers analysing the formation of seed coat colour and the underlying regulatory mechanisms in B. napus.

Comparative genomic analyses reveal the genetic basis of the yellow-seed trait in Brassica napus.

Yellow-seed trait is a desirable breeding characteristic of rapeseed (Brassica napus) that could greatly improve seed oil yield and quality. However, the underlying mechanisms controlling this phenotype in B. napus plants are difficult to discern because of their complexity. Here, we assemble high-quality genomes of yellow-seeded (GH06) and black-seeded (ZY821). Combining in-depth fine mapping of a quantitative trait locus (QTL) for seed color with other omics data reveal BnA09MYB47a, encoding an R2R3-MYB-type transcription factor, as the causal gene of a major QTL controlling the yellow-seed trait. Functional studies show that sequence variation of BnA09MYB47a underlies the functional divergence between the yellow- and black-seeded B. napus. The black-seed allele BnA09MYB47aZY821, but not the yellow-seed allele BnA09MYB47aGH06, promotes flavonoid biosynthesis by directly activating the expression of BnTT18. Our discovery suggests a possible approach to breeding B. napus for improved commercial value and facilitates flavonoid biosynthesis studies in Brassica crops.

Metabolite Characteristics Analysis of Siliques and Effects of Lights on the Accumulation of Glucosinolates in Siliques of Rapeseed.

Glucosinolates (GSLs) are naturally occurring secondary metabolites found in the Brassicaceae family, which mainly synthesize in the siliques with a wide range of functions. In this study, we investigated the effects of lights on metabolites in siliques of rapeseed through ultra high-performance liquid chromatography (UPLC)-heated electrospray ionization (HESI)-tandem mass spectrometry (MS/MS). A total of 249 metabolites, including 29 phenolic acids, 38 flavonoids, 22 GSLs, 93 uncalculated and 67 unknown compounds, were identified in siliques of rapeseed. Meanwhile, 62 metabolites showed significant differences after shading treatment, which were mainly GSLs and unknown compounds. Interestingly, the amounts of 10 GSLs had high accumulation levels in siliques, while the expression levels of their corresponding biosynthetic genes (AOP, GSL-OH, IGMT, and ST5a) were obviously reduced after shading treatment. Further evidence showed that the amounts of GSLs were significantly reduced in seeds, in accordance with the expression profiles of transporter genes (BnaGTRs). Our findings indicated that lights could affect the accumulation and transportation of GSLs from siliques to seeds in rapeseed. Therefore, this study facilitates a better understanding of metabolic characteristics of siliques and provides insight into the importance of light for GSLs accumulation and transportation in siliques and seeds of rapeseed.

Genome-Wide Analysis of the YABBY Transcription Factor Family in Rapeseed (Brassica napus L.).

The YABBY family of plant-specific transcription factors play important regulatory roles during the development of leaves and floral organs, but their functions in Brassica species are incompletely understood. Here, we identified 79 YABBY genes from Arabidopsis thaliana and five Brassica species (B. rapa, B. nigra, B. oleracea, B. juncea, and B. napus). A phylogenetic analysis of YABBY proteins separated them into five clusters (YAB1-YAB5) with representatives from all five Brassica species, suggesting a high degree of conservation and similar functions within each subfamily. We determined the gene structure, chromosomal location, and expression patterns of the 21 BnaYAB genes identified, revealing extensive duplication events and gene loss following polyploidization. Changes in exon-intron structure during evolution may have driven differentiation in expression patterns and functions, combined with purifying selection, as evidenced by Ka/Ks values below 1. Based on transcriptome sequencing data, we selected nine genes with high expression at the flowering stage. qRT-PCR analysis further indicated that most BnaYAB family members are tissue-specific and exhibit different expression patterns in various tissues and organs of B. napus. This preliminary study of the characteristics of the YABBY gene family in the Brassica napus genome provides theoretical support and reference for the later functional identification of the family genes.

Characterization of cold stress responses in different rapeseed ecotypes based on metabolomics and transcriptomics analyses.

The winter oilseed ecotype is more tolerant to low temperature than the spring ecotype. Transcriptome and metabolome analyses of leaf samples of five spring Brassica napus L. (B. napus) ecotype lines and five winter B. napus ecotype lines treated at 4 °C and 28 °C were performed. A total of 25,460 differentially expressed genes (DEGs) of the spring oilseed ecotype and 28,512 DEGs of the winter oilseed ecotype were identified after cold stress; there were 41 differentially expressed metabolites (DEMs) in the spring and 47 in the winter oilseed ecotypes. Moreover, more than 46.2% DEGs were commonly detected in both ecotypes, and the extent of the changes were much more pronounced in the winter than spring ecotype. By contrast, only six DEMs were detected in both the spring and winter oilseed ecotypes. Eighty-one DEMs mainly belonged to primary metabolites, including amino acids, organic acids and sugars. The large number of specific genes and metabolites emphasizes the complex regulatory mechanisms involved in the cold stress response in oilseed rape. Furthermore, these data suggest that lipid, ABA, secondary metabolism, signal transduction and transcription factors may play distinct roles in the spring and winter ecotypes in response to cold stress. Differences in gene expression and metabolite levels after cold stress treatment may have contributed to the cold tolerance of the different oilseed ecotypes.

Genome-wide identification AINTEGUMENTA-like (AIL) genes in Brassica species and expression patterns during reproductive development in Brassica napus L.

The AINTEGUMENTA-like (AIL) proteins, which belong to the AP2 family, play important roles in regulating the growth and development of plant organs. The AIL family has not yet been comprehensively studied in rapeseed (Brassica napus), an allotetraploid and model organism for the study of polyploid evolution. In the present study, 99 AIL family genes were identified and characterized from B. rapa, B. oleracea, B. napus, B. juncea, and B. nigra using a comprehensive genome-wide study, including analyses of phylogeny, gene structure, chromosomal localization, and expression pattern. Using a phylogenetic analysis, the AIL genes were divided into eight groups, which were closely related to the eight AtAIL genes, and which shared highly conserved structural features within the same subfamily. The non-synonymous/synonymous substitution ratios of the paralogs and orthologs were less than 1, suggesting that the AIL genes mainly experienced purifying selection during evolution. In addition, the RNA sequencing data and qRT-PCR analysis revealed that the B. napus AIL genes exhibited organ- and developmental stage-specific expression patterns. Certain genes were highly expressed in the developing seeds (BnaAIL1, BnaAIL2, BnaAIL5, and BnaAIL6), the roots (BnaANT, BnaAIL5, and BnaAIL6), and the stem (BnaAIL7B). Our results provide valuable information for further functional analysis of the AIL family in B. napus and related Brassica species.

Comparative Analysis of the Metabolic Profiles of Yellow- versus Black-Seeded Rapeseed Using UPLC-HESI-MS/MS and Transcriptome Analysis.

The high levels of secondary metabolites in rapeseed play important roles in determining the oil quality and feeding value. Here, we characterized the metabolic profiles in seeds of various yellow- and black-seeded rapeseed accessions. Two hundred and forty-eight features were characterized, including 31 phenolic acids, 54 flavonoids, 24 glucosinolates, 65 lipid compounds, and 74 other polar compounds. The most abundant phenolic acids and various flavonoids (epicatechin, isorhamnetin, kaempferol, quercetin, and their derivatives) were widely detected and showed significant differences in distribution between the yellow- and black-seeded rapeseed. Furthermore, the related genes (e.g., BnTT3, BnTT18, BnTT10, BnTT12, and BnBAN) involved in the proanthocyanidin pathway had lower expression levels in yellow-seeded rapeseed, strongly suggesting that the seed coat color could be mainly determined by the levels of epicatechin and their derivatives. These results improve our understanding of the primary constituents of rapeseed and lay the foundation for breeding novel varieties with a high nutritional value.