Unraveling the mechanism of flower color variation in Brassica napus by integrated metabolome and transcriptome analyses.

Brassica napus is one of the most important oil crops in the world. Breeding oilseed rape with colorful flowers can greatly enhance the ornamental value of B. napus and thus improve the economic benefits of planting. As water-soluble flavonoid secondary metabolites, anthocyanins are very important for the synthesis and accumulation of pigments in the petals of plants, giving them a wide range of bright colors. Despite the documentation of over 60 distinct flower shades in B. napus, the intricacies underlying flower color variation remain elusive. Particularly, the mechanisms driving color development across varying flower color backgrounds necessitate further comprehensive investigation. This research undertook a comprehensive exploration through the integration of transcriptome and metabolome analyses to pinpoint pivotal genes and metabolites underpinning an array of flower colors, including beige, beige-red, yellow, orange-red, deep orange-red, white, light-purple, and purple. First, we used a two-way BLAST search to find 275 genes in the reference genome of B. napus Darmor v10 that were involved in making anthocyanins. The subsequent scrutiny of RNA-seq outcomes underscored notable upregulation in the structural genes F3H and UGT, alongside the MYB75, GL3, and TTG1 transcriptional regulators within petals, showing anthocyanin accumulation. By synergizing this data with a weighted gene co-expression network analysis, we identified CHS, F3H, MYB75, MYB12, and MYB111 as the key players driving anthocyanin synthesis in beige-red, orange-red, deep orange-red, light-purple, and purple petals. By integrating transcriptome and weighted gene co-expression network analysis findings with anthocyanin metabolism data, it is hypothesized that the upregulation of MYB75, which, in turn, enhances F3H expression, plays a pivotal role in the development of pigmented oilseed rape flowers. These findings help to understand the transcriptional regulation of anthocyanin biosynthesis in B. napus and provide valuable genetic resources for breeding B. napus varieties with novel flower colors.

Analysis of Altered Flowering Related Genes in a Multi-Silique Rapeseed (Brassica napus L.) Line zws-ms Based on Combination of Genome, Transcriptome and Proteome Data.

Based on previous researches, we further investigated the multi-silique trait in rapeseed (Brassica napus L.) line zws-ms. In this study, we used a relatively comprehensive list of flowering related genes in rapeseed and compared them between zws-ms and its near-isogenic line (NIL) zws-217. Genes were studied on genome, transcriptome and proteome levels and then we focused on genes with non-synonymous single nucleotide polymorphism (SNP) or frame-shift insertion-deletion (InDel), finding some genes on the list which changes their sequences. Then, combined with their annotation and the information of their orthologs, certain genes such as BnaA09g05900D, ortholog of AGAMOUS-LIKE 42 (AGL42), which encodes an MADS-box protein, were assumed as probably responsible for the multi-silique trait. Also, we analyzed the Differentially Accumulated Proteins (DAPs) between zws-ms and zws-217, revealing some genes involved in homologous recombination and mismatch repair pathways. Since the development of flowers/siliques is crucial to crops and it influences the yield of rapeseed, this study paved a way to deeply understand the mechanism of the multi-pistil flower formation, which may facilitate researches on rapeseed production in future.

Investigation of Thermomorphogenesis-Related Genes for a Multi-Silique Trait in Brassica napus by Comparative Transcriptome Analysis.

In higher plants, the structure of a flower is precisely controlled by a series of genes. An aberrance flower results in abnormal fruit morphology. Previously, we reported multi-silique rapeseed (Brassica napus) line zws-ms. We identified two associated regions and investigated differentially expressed genes (DEGs); thus, some candidate genes underlying the multi-silique phenotype in warm area Xindu were selected. However, this phenotype was switched off by lower temperature, and the responsive genes, known as thermomorphogenesis-related genes, remained elusive. So, based on that, in this study, we further investigated the transcriptome data from buds of zws-ms and its near-isogenic line zws-217 grown in colder area Ma'erkang, where both lines showed normal siliques only, and the DEGs between them analyzed. We compared the 129 DEGs from Xindu to the 117 ones from Ma'erkang and found that 33 of them represented the same or similar expression trends, whereas the other 96 DEGs showed different expression trends, which were defined as environment-specific. Furthermore, we combined this with the gene annotations and ortholog information and then selected BnaA09g45320D (chaperonin gene CPN10-homologous) and BnaC08g41780D [Seryl-tRNA synthetase gene OVULE ABORTION 7 (OVA7)-homologous] the possible thermomorphogenesis-related genes, which probably switched off the multi-silique under lower temperature. This study paves a way to a new perspective into flower/fruit development in Brassica plants.

Detection of Blackleg Resistance Gene Rlm1 in Double-Low Rapeseed Accessions from Sichuan Province, by Kompetitive Allele-Specific PCR.

Blackleg is a serious disease in Brassica plants, causing moderate to severe yield losses in rapeseed worldwide. Although China has not suffered from this disease yet (more aggressive Leptosphaeria maculans is not present yet), it is crucial to take provisions in breeding for disease resistance to have excellent blackleg-resistant cultivars already in the fields or in the breeding pipeline. The most efficient strategy for controlling this disease is breeding plants with identified resistance genes. We selected 135 rapeseed accessions in Sichuan, including 30 parental materials and 105 hybrids, and we determined their glucosinolate and erucic acid content and confirmed 17 double-low materials. A recently developed single-nucleotide polymorphism (SNP) marker, SNP_208, was used to genotype allelic Rlm1/rlm1 on chromosome A07, and 87 AvrLm1-resistant materials. Combined with the above-mentioned seed quality data, we identified 11 AvrLm1-resistant double-low rapeseed accessions, including nine parental materials and two hybrids. This study lays the foundation of specific R gene-oriented breeding, in the case that the aggressive Leptosphaeria maculans invades and establishes in China in the future and a robust and less labor consuming method to identify resistance in canola germplasm.

Investigation for a multi-silique trait in Brassica napus by alternative splicing analysis.

BACKGROUND: Flower and fruit development are vital stages of the angiosperm lifecycle. We previously investigated the multi-silique trait in the rapeseed (Brassica napus) line zws-ms on a genomic and transcriptomic level, leading to the identification of two genomic regions and several candidate genes associated with this trait. However, some events on the transcriptome level, like alternative splicing, were poorly understood. METHODS: Plants from zws-ms and its near-isogenic line (NIL) zws-217 were both grown in Xindu with normal conditions and a colder area Ma'erkang. Buds from the two lines were sampled and RNA was isolated to perform the transcriptomic sequencing. The numbers and types of alternative splicing (AS) events from the two lines were counted and classified. Genes with AS events and expressed differentially between the two lines, as well as genes with AS events which occurred in only one line were emphasized. Their annotations were further studied. RESULTS: From the plants in Xindu District, an average of 205,496 AS events, which could be sorted into five AS types, were identified. zws-ms and zws-217 shared highly similar ratios of each AS type: The alternative 5' and 3' splice site types were the most common, while the exon skipping type was observed least often. Eleven differentially expressed AS genes were identified, of which four were upregulated and seven were downregulated in zws-ms. Their annotations implied that five of these genes were directly associated with the multi-silique trait. While samples from colder area Ma'erkang generated generally reduced number of each type of AS events except for Intron Retention; but the number of differentially expressed AS genes increased significantly. Further analysis found that among the 11 differentially expressed AS genes from Xindu, three of them maintained the same expression models, while the other eight genes did not show significant difference between the two lines in expression level. Additionally, the 205 line-specific expressed AS genes were analyzed, of which 187 could be annotated, and two were considered to be important. DISCUSSION: This study provides new insights into the molecular mechanism of the agronomically important multi-silique trait in rapeseed on the transcriptome level and screens out some environment-responding candidate genes.

Fine Mapping of a Locus Underlying the Ectopic Blade-Like Outgrowths on Leaf and Screening Its Candidate Genes in Rapeseed (Brassica napus L.).

Leaf is an important organ for higher plants, and the shape of it is one of the crucial traits of crops. In this study, we investigated a unique aberrant leaf morphology trait in a mutational rapeseed material, which displayed ectopic blade-like outgrowths on the adaxial side of leaf. The abnormal line 132000B-3 was crossed with the normal line 827-3. Based on the F2 : 3 family, we constructed two DNA pools (normal pool and abnormal pool) by the bulked segregant analysis (BSA) method and performed whole genome re-sequencing (WGR), obtaining the single-nucleotide polymorphism (SNP) and insertion/deletion (InDel) data. The SNP-index method was used to calculate the Δ(SNP/InDel-index), and then an association region was identified on chromosome A10 with a length of 5.5 Mbp, harboring 1048 genes totally. Subsequently, the fine mapping was conducted by using the penta-primer amplification refractory mutation system (PARMS), and the associated region was narrowed down to a 35.1-kbp segment, containing only seven genes. These seven genes were then analyzed according to their annotations and finally, BnA10g0422620 and BnA10g0422610, orthologs of LATE MERISTEM IDENTITY1 (LMI1) gene from Arabidopsis and REDUCED COMPLEXITY (RCO) gene from its relative Cardamine hirsuta, respectively, were identified as the candidate genes responding to this blade-like outgrowth trait in rapeseed. This study provides a novel perspective into the leaf formation in Brassica plants.

Identification of genomic regions associated with multi-silique trait in Brassica napus.

BACKGROUND: Although rapeseed (Brassica napus L.) mutant forming multiple siliques was morphologically described and considered to increase the silique number per plant, an important agronomic trait in this crop, the molecular mechanism underlying this beneficial trait remains unclear. Here, we combined bulked-segregant analysis (BSA) and whole genome re-sequencing (WGR) to map the genomic regions responsible for the multi-silique trait using two pools of DNA from the near-isogenic lines (NILs) zws-ms (multi-silique) and zws-217 (single-silique). We used the Euclidean Distance (ED) to identify genomic regions associated with this trait based on both SNPs and InDels. We also conducted transcriptome sequencing to identify differentially expressed genes (DEGs) between zws-ms and zws-217. RESULTS: Genetic analysis using the ED algorithm identified three SNP- and two InDel-associated regions for the multi-silique trait. Two highly overlapped parts of the SNP- and InDel-associated regions were identified as important intersecting regions, which are located on chromosomes A09 and C08, respectively, including 2044 genes in 10.20-MB length totally. Transcriptome sequencing revealed 129 DEGs between zws-ms and zws-217 in buds, including 39 DEGs located in the two abovementioned associated regions. We identified candidate genes involved in multi-silique formation in rapeseed based on the results of functional annotation. CONCLUSIONS: This study identified the genomic regions and candidate genes related to the multi-silique trait in rapeseed.

Development of a Fertility Restorer for inap CMS (Isatis indigotica) Brassica napus Through Genetic Introgression of One Alien Addition.

Novel Brassica napus cytoplasmic male sterility (CMS) with carpelloid stamens (inap CMS) was produced by intertribal somatic hybridization with Isatis indigotica (Chinese woad), but its RF (restorer of fertility) gene(s) existed in one particular woad chromosome that was carried by one fertile monosomic alien addition line (MAAL) of rapeseed. Herein, the selfed progenies of this MAAL were extensively selected and analyzed to screen the rapeseed-type plants (2n = 38) with good male fertility and to produce their doubled haploid (DH) lines by microspore culture. From the investigation of fertility restoration in the F1 hybrids with inap CMS, one DH line (RF 39) was identified to adequately restore male fertility and likely carried one dominant RF gene. Specifically, this restorer produced brown pollen grains, similar to the woad and the MAAL, suggesting that this trait is closely linked with the RF gene(s) and serves as one phenotypic marker for the restorer. This restorer contained 38 chromosomes of rapeseed and no intact chromosomes of woad, but some DNA fragments of woad origin were detected at low frequency. This restorer was much improved for pollen and seed fertility and for low glucosinolate content. The successful breeding of the restorer for inap CMS rendered this new pollination control system feasible for rapeseed hybrid production.

Genome-Wide Identification of Flowering-Time Genes in Brassica Species and Reveals a Correlation between Selective Pressure and Expression Patterns of Vernalization-Pathway Genes in Brassica napus.

Flowering time is a key agronomic trait, directly influencing crop yield and quality. Many flowering-time genes have been identified and characterized in the model plant Arabidopsis thaliana; however, these genes remain uncharacterized in many agronomically important Brassica crops. In this study, we identified 1064, 510, and 524 putative orthologs of A. thaliana flowering-time genes from Brassica napus, Brassica rapa, and Brassica oleracea, respectively, and found that genes involved in the aging and ambient temperature pathways were fewer than those in other flowering pathways. Flowering-time genes were distributed mostly on chromosome C03 in B. napus and B. oleracea, and on chromosome A09 in B. rapa. Calculation of non-synonymous (Ka)/synonymous substitution (Ks) ratios suggested that flowering-time genes in vernalization pathways experienced higher selection pressure than those in other pathways. Expression analysis showed that most vernalization-pathway genes were expressed in flowering organs. Approximately 40% of these genes were highly expressed in the anther, whereas flowering-time integrator genes were expressed in a highly organ-specific manner. Evolutionary selection pressures were negatively correlated with the breadth and expression levels of vernalization-pathway genes. These findings provide an integrated framework of flowering-time genes in these three Brassica crops and provide a foundation for deciphering the relationship between gene expression patterns and their evolutionary selection pressures in Brassica napus.

Genome-wide identification, putative functionality and interactions between lncRNAs and miRNAs in Brassica species.

Non-coding RNA (ncRNA) is abundant in plant genomes, but is poorly described with unknown functionality in most species. Using whole genome RNA sequencing, we identified 1885, 1910 and 1299 lncRNAs and 186, 157 and 161 miRNAs at the whole genome level in the three Brassica species B. napus, B. oleracea and B. rapa, respectively. The lncRNA sequences were divergent between the three Brassica species. One quarter of lncRNAs were located in tandem repeat (TR) region. The expression of both lncRNAs and miRNAs was strongly biased towards the A rather than the C subgenome in B. napus, unlike mRNA expression. miRNAs in genic regions had higher average expression than miRNAs in non-genic regions in B. napus and B. oleracea. We provide a comprehensive reference for the distribution, functionality and interactions of lncRNAs and miRNAs in Brassica.

The high-quality genome of Brassica napus cultivar 'ZS11' reveals the introgression history in semi-winter morphotype.

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi-winter oilseed rape cultivar 'ZS11' and its comprehensive genomic comparison with the genomes of the winter-type cultivar 'Darmor-bzh' as well as two progenitors. The integrated BAC-to-BAC and whole-genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high-quality genome assembly of B. napus 'ZS11'. Within a short evolutionary period (~6700 years ago), semi-winter-type 'ZS11' and the winter-type 'Darmor-bzh' maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to 'Darmor-bzh', both two subgenomes of 'ZS11' are closely related to its progenitors, and the 'ZS11' genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi-winter-type 'ZS11' underwent potential genomic introgressions with B. rapa (Ar ). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization-responsive flowering time in 'ZS11' was first experienced HE, and then underwent genomic introgression event with Ar , which potentially has led to genetic differences in controlling vernalization in the semi-winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi-winter oilseed rape in Asia.

Transcriptome Analysis of Interspecific Hybrid between Brassica napus and B. rapa Reveals Heterosis for Oil Rape Improvement.

The hybrid between Brassica napus and B. rapa displays obvious heterosis in both growth performance and stress tolerances. A comparative transcriptome analysis for B. napus (A(n)A(n)CC genome), B. rapa (A(r)A(r) genome), and its hybrid F1 (A(n)A(r)C genome) was carried out to reveal the possible molecular mechanisms of heterosis at the gene expression level. A total of 40,320 nonredundant unigenes were identified using B. rapa (AA genome) and B. oleracea (CC genome) as reference genomes. A total of 6,816 differentially expressed genes (DEGs) were mapped in the A and C genomes with 4,946 DEGs displayed nonadditively by comparing the gene expression patterns among the three samples. The coexistence of nonadditive DEGs including high-parent dominance, low-parent dominance, overdominance, and underdominance was observed in the gene action modes of F1 hybrid, which were potentially related to the heterosis. The coexistence of multiple gene actions in the hybrid was observed and provided a list of candidate genes and pathways for heterosis. The expression bias of transposable element-associated genes was also observed in the hybrid compared to their parents. The present study could be helpful for the better understanding of the determination and regulation of mechanisms of heterosis to aid Brassica improvement.

The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes.

Polyploidization has provided much genetic variation for plant adaptive evolution, but the mechanisms by which the molecular evolution of polyploid genomes establishes genetic architecture underlying species differentiation are unclear. Brassica is an ideal model to increase knowledge of polyploid evolution. Here we describe a draft genome sequence of Brassica oleracea, comparing it with that of its sister species B. rapa to reveal numerous chromosome rearrangements and asymmetrical gene loss in duplicated genomic blocks, asymmetrical amplification of transposable elements, differential gene co-retention for specific pathways and variation in gene expression, including alternative splicing, among a large number of paralogous and orthologous genes. Genes related to the production of anticancer phytochemicals and morphological variations illustrate consequences of genome duplication and gene divergence, imparting biochemical and morphological variation to B. oleracea. This study provides insights into Brassica genome evolution and will underpin research into the many important crops in this genus.