Search

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.

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.

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.

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. Electronic supplementary material The online version of this article (10.1186/s12864-019-5675-4) contains supplementary material, which is available to authorized users.

Plant height is a vital agronomic trait for crops, including oilseed crops such as rapeseed (Brassica napus L.). It affects the crop yield, oil content, and lodging resistance in rapeseed. In this study, we investigated a dwarf trait controlled by a semi-dominant allele in rapeseed. A dwarf line, YA2016-12, was crossed with a tall line, G184-189, and an F2 population was established. Forty of the tallest plants and 40 of the shortest plants from the F2 population were selected and two DNA pools (tall and dwarf) were constructed by the bulked segregant analysis (BSA) method. The two DNA pools and two parental DNAs were then re-sequenced. A sliding window analysis was used to calculate the &Delta, (SNP-index) and discover an association region on chromosome A03 with a length of 12.4 Mb. Within this region, we found 1225 genes, including 811 genes with non-synonymous or frameshift mutations between YA2016-12 and G184-189. Alignment to known plant height-related orthologs in Arabidopsis thaliana, as well as KEGG pathway and gene ontology annotations, was used to identify nine candidate genes (BnaA03g31770D, BnaA03g37960D, BnaA03g24740D, BnaA03g40550D, BnaA03g26120D, BnaA03g35130D, BnaA03g42350D, BnaA03g25610D, and BnaA03g39850D) involved in gibberellin or cytokinin signaling. Identification of the causal gene for this trait, and of genetic markers linked to favorable alleles, has potential utility for marker-assisted selection to breed rapeseed varieties with improved height.

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.

Table S2. Information about the mapped reads based on the transcriptome sequencing data. (DOCX 15 kb)

Table S7. GO annotations of the 104 genes with non-synonymous mutation SNPs in the associated regions. (DOCX 34 kb)

Table S8. GO annotations of the 19 genes with frame-shift mutation InDels in the associated regions. (DOCX 17 kb)

Figure S1. Whole genome with associated regions. (DOCX 209 kb)

Table S1. Summary of the transcriptome sequencing data. (DOCX 15 kb)

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.

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.

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.

The hybrid betweenBrassica napusandB. rapadisplays obvious heterosis in both growth performance and stress tolerances. A comparative transcriptome analysis forB. napus(AnAnCCgenome),B. rapa(ArArgenome), and its hybrid F1 (AnArCgenome) 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 usingB. rapa(AA genome) andB. 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 aidBrassicaimprovement.

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., Brassica oleracea is plant species comprising economically important vegetable crops. Here, the authors report the draft genome sequence of B. oleracea and, through a comparative analysis with the closely related B. rapa, reveal insights into Brassica evolution and divergence of interspecific genomes and intraspecific subgenomes.