Genome-wide survey of flavonoid biosynthesis genes and gene expression analysis between black- and yellow-seeded Brassica napus

Flavonoids, the compounds that impart color to fruits, flowers, and seeds, are the most widespread secondary metabolites in plants. However, a systematic analysis of these loci has not been performed in Brassicaceae. In this study, we isolated 649 nucleotide sequences related to flavonoid biosynthesis, i.e., the Transparent Testa (TT) genes, and their associated amino acid sequences in 17 Brassicaceae species, grouped into Arabidopsis or Brassicaceae subgroups. Moreover, 36 copies of 21 genes of the flavonoid biosynthesis pathway were identified in A. thaliana, 53 were identified in B. rapa, 50 in B. oleracea, and 95 in B. napus, followed the genomic distribution, collinearity analysis and genes triplication of them among Brassicaceae species. The results showed that the extensive gene loss, whole genome triplication, and diploidization that occurred after divergence from the common ancestor. Using qRT-PCR methods, we analyzed the expression of eighteen flavonoid biosynthesis genes in 6 yellow- and black-seeded B. napus inbred lines with different genetic background, found that 12 of which were preferentially expressed during seed development, whereas the remaining genes were expressed in all B. napus tissues examined. Moreover, fourteen of these genes showed significant differences in expression level during seed development, and all but four of these (i.e., BnTT5, BnTT7, BnTT10, and BnTTG1) had similar expression patterns among the yellow- and black-seeded B. napus. Results showed that the structural genes (BnTT3, BnTT18 and BnBAN), regulatory genes (BnTTG2 and BnTT16) and three encoding transfer proteins (BnTT12, BnTT19, and BnAHA10) might play an crucial roles in the formation of different seed coat colors in B. napus. These data will be helpful for illustrating the molecular mechanisms of flavonoid biosynthesis in Brassicaceae species

Fine mapping of the major QTL for seed coat color in Brassica rapa var. Yellow Sarson by use of NIL populations and transcriptome sequencing for identification of the candidate genes.

Yellow seed is a desirable trait in Brassica oilseed crops. The B. rapa var. Yellow Sarson carry unique yellow seed color genes which are not only important for the development of yellow-seeded oilseed B. rapa cultivars but this variant can also be used to develop yellow-seeded B. napus. In this study, we developed near-isogenic lines (NILs) of Yellow Sarson for the major seed coat color QTL SCA9-2 of the chromosome A9 and used the NILs to fine map this QTL region and to identify the candidate genes through linkage mapping and transcriptome sequencing of the developing seeds. From the 18.4 to 22.79 Mb region of SCA9-2, six SSR markers showing 0.63 to 5.65% recombination were developed through linkage analysis and physical mapping. A total of 55 differentially expressed genes (DEGs) were identified in the SCA9-2 region through transcriptome analysis; these included three transcription factors, Bra028039 (NAC), Bra023223 (C2H2 type zinc finger), Bra032362 (TIFY), and several other genes which encode unknown or nucleic acid binding protein; these genes might be the candidates and involved in the regulation of seed coat color in the materials used in this study. Several biosynthetic pathways, including the flavonoid, phenylpropanoid and suberin biosynthetic pathways were significantly enriched through GO and KEGG enrichment analysis of the DEGs. This is the first comprehensive study to understand the yellow seed trait of Yellow Sarson through employing linkage mapping and global transcriptome analysis approaches

Gene silencing of BnTT10 family genes causes retarded pigmentation and lignin reduction in the seed coat of Brassica napus.

Yellow-seed (i.e., yellow seed coat) is one of the most important agronomic traits of Brassica plants, which is correlated with seed oil and meal qualities. Previous studies on the Brassicaceae, including Arabidopsis and Brassica species, proposed that the seed-color trait is correlative to flavonoid and lignin biosynthesis, at the molecular level. In Arabidopsis thaliana, the oxidative polymerization of flavonoid and biosynthesis of lignin has been demonstrated to be catalyzed by laccase 15, a functional enzyme encoded by the AtTT10 gene. In this study, eight Brassica TT10 genes (three from B. napus, three from B. rapa and two from B. oleracea) were isolated and their roles in flavonoid oxidation/polymerization and lignin biosynthesis were investigated. Based on our phylogenetic analysis, these genes could be divided into two groups with obvious structural and functional differentiation. Expression studies showed that Brassica TT10 genes are active in developing seeds, but with differential expression patterns in yellow- and black-seeded near-isogenic lines. For functional analyses, three black-seeded B. napus cultivars were chosen for transgenic studies. Transgenic B. napus plants expressing antisense TT10 constructs exhibited retarded pigmentation in the seed coat. Chemical composition analysis revealed increased levels of soluble proanthocyanidins, and decreased extractable lignin in the seed coats of these transgenic plants compared with that of the controls. These findings indicate a role for the Brassica TT10 genes in proanthocyanidin polymerization and lignin biosynthesis, as well as seed coat pigmentation in B. napus

BnKAT2 Positively Regulates the Main Inflorescence Length and Silique Number in Brassica napus by Regulating the Auxin and Cytokinin Signaling Pathways

Brassica napus is the dominant oil crop cultivated in China for its high quality and high yield. The length of the main inflorescence and the number of siliques produced are important traits contributing to rapeseed yield. Therefore, studying genes related to main inflorescence and silique number is beneficial to increase rapeseed yield. Herein, we focused on the effects of BnKAT2 on the main inflorescence length and silique number in B. napus. We explored the mechanism of BnKAT2 increasing the effective length of main inflorescence and the number of siliques through bioinformatics analysis, transgenic technology, and transcriptome sequencing analysis. The full BnKAT2(BnaA01g09060D) sequence is 3674 bp, while its open reading frame is 2055 bp, and the encoded protein comprises 684 amino acids. BnKAT2 is predicted to possess two structural domains, namely KHA and CNMP-binding domains. The overexpression of BnKAT2 effectively increased the length of the main inflorescence and the number of siliques in B. napus, as well as in transgenic Arabidopsis thaliana. The type-A Arabidopsis response regulator (A-ARR), negative regulators of the cytokinin, are downregulated in the BnKAT2-overexpressing lines. The Aux/IAA, key genes in auxin signaling pathways, are downregulated in the BnKAT2-overexpressing lines. These results indicate that BnKAT2 might regulate the effective length of the main inflorescence and the number of siliques through the auxin and cytokinin signaling pathways. Our study provides a new potential function gene responsible for improvement of main inflorescence length and silique number, as well as a candidate gene for developing markers used in MAS (marker-assisted selection) breeding to improve rapeseed yield

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

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

Molecular mapping and QTL for expression profiles of flavonoid genes in Brassica napus

Flavonoids are secondary metabolites that are extensively distributed in the plant kingdom and contribute to seed coat color formation in rapeseed. To decipher the genetic networks underlying flavonoid biosynthesis in rapeseed, we constructed a high-density genetic linkage map with 1089 polymorphic loci (including 464 SSR loci, 97 RAPD loci, 451 SRAP loci, and 75 IBP loci) using recombinant inbred lines (RILs). The map consists of 19 linkage groups and covers 2,775 cM of the B. napus genome with an average distance of 2.54 cM between adjacent markers. We then performed expression quantitative trait locus (eQTL) analysis to detect transcript-level variation of 18 flavonoid biosynthesis pathway genes in the seeds of the 94 RILs. In total, 72 eQTLs were detected and found to be distributed among 15 different linkage groups that account for 4.11% to 52.70% of the phenotypic variance atrributed to each eQTL. Using a genetical genomics approach, four eQTL hotspots together harboring 28 eQTLs associated with 18 genes were found on chromosomes A03, A09, and C08 and had high levels of synteny with genome sequences of A. thaliana and Brassica species. Associated with the trans-eQTL hotspots on chromosomes A03, A09, and C08 were 5, 17, and 1 genes encoding transcription factors, suggesting that these genes have essential roles in the flavonoid biosynthesis pathway. Importantly, bZIP25, which is expressed specifically in seeds, MYC1, which controls flavonoid biosynthesis, and the R2R3-type gene MYB51, which is involved in the synthesis of secondary metabolites, were associated with the eQTL hotspots, and these genes might thus be involved in different flavonoid biosynthesis pathways in rapeseed. Hence, further studies of the functions of these genes will provide insight into the regulatory mechanism underlying flavonoid biosynthesis, and lay the foundation for elaborating the molecular mechanism of seed coat color formation in B. napus

Identification and characterization of a curly-leaf locus CL1 encoding an IAA2 protein in Brassica napus

The leaf is the main organ for rapeseed photosynthesis, and its morphology influences photosynthetic efficiency and supports increased planting density and yield. However, the molecular regulatory mechanism of leaf morphology in Brassica napus is poorly understood, restricting progress in breeding for the trait. We describe a novel dominant mutation, curly leaf 1 (cl1), which confers uneven dorsal–ventral axis development, irregular cellular structure and influenced gravitropic response in the seedling stage. The CL1 locus was mapped to a 1.573-Mb interval on chromosome A05 using simple sequence repeat (SSR) markers, and co-segregated with the phenotype of plants in the curly F2 population. A substitution (P62S) was identified in the highly conserved degron motif (GWSPV) of the IAA2 protein in the cl1 mutant, and the P62S substitution impaired the interaction between IAA2 and TIR1 in the presence of auxin, influencing auxin signaling. The P62S substitution-induced curly leaf phenotype was verified by ectopic expression of BnaA05.iaa2 in Arabidopsis and B. napus. Our findings explain the function of IAA2 in rapeseed, providing a foundation for future investigation of auxin signaling and the mechanisms underlying leaf development in B. napus

Genome-Wide Identification and Expression Profiling of Monosaccharide Transporter Genes Associated with High Harvest Index Values in Rapeseed (Brassica napus L.)

Sugars are important throughout a plant’s lifecycle. Monosaccharide transporters (MST) are essential sugar transporters that have been identified in many plants, but little is known about the evolution or functions of MST genes in rapeseed (Brassica napus). In this study, we identified 175 MST genes in B. napus, 87 in Brassica oleracea, and 83 in Brassica rapa. These genes were separated into the sugar transport protein (STP), polyol transporter (PLT), vacuolar glucose transporter (VGT), tonoplast monosaccharide transporter (TMT), inositol transporter (INT), plastidic glucose transporter (pGlcT), and ERD6-like subfamilies, respectively. Phylogenetic and syntenic analysis indicated that gene redundancy and gene elimination have commonly occurred in Brassica species during polyploidization. Changes in exon-intron structures during evolution likely resulted in the differences in coding regions, expression patterns, and functions seen among BnMST genes. In total, 31 differentially expressed genes (DEGs) were identified through RNA-seq among materials with high and low harvest index (HI) values, which were divided into two categories based on the qRT-PCR results, expressed more highly in source or sink organs. We finally identified four genes, including BnSTP5, BnSTP13, BnPLT5, and BnERD6-like14, which might be involved in monosaccharide uptake or unloading and further affect the HI of rapeseed. These findings provide fundamental information about MST genes in Brassica and reveal the importance of BnMST genes to high HI in B. napus