Development of B. carinata with super-high erucic acid content through interspecific hybridization.

KEY MESSAGE: Disomic alien chromosome addition Brassica carinata lines with super-high erucic acid content were developed through interspecific hybridization with B. juncea and characterized using molecular, cytological and biochemical techniques. Brassica carinata [A.] Braun (BBCC, 2n = 34) is a climate-resilient oilseed. Its seed oil is high in erucic acid (> 40%), rendering it well suited for the production of biofuel and other bio-based applications. To enhance the competitiveness of B. carinata with high erucic B. napus (HEAR), lines with super-high erucic acid content were developed through interspecific hybridization. To this end, a fad2B null allele from Brassica juncea (AABB, 2n = 36) was introgressed into B. carinata, resulting in a B. carinata fad2B mutant with erucic acid levels of over 50%. Subsequently, the FAE allele from B. rapa spp. yellow sarson (AA, 2n = 20) was transferred to the fad2B B. carinata line, yielding lines with erucic acid contents of up to 57.9%. Molecular analysis using the Brassica 90 K Illumina Infinium™ SNP genotyping array identified these lines as disomic alien chromosome addition lines, with two extra A08 chromosomes containing the BrFAE gene. The alien chromosomes from B. rapa were clearly distinguished by molecular cytogenetics in one of the addition lines. Analysis of microspore-derived offspring and hybrids from crosses with a CMS B. carinata line showed that the transfer rate of the A08 chromosome into male gametes was over 98%, resulting in almost completely stable transmission of an A08 chromosome copy into the progeny. The increase in erucic acid levels was accompanied by changes in the proportions of other fatty acids depending on the genetic changes that were introduced in the interspecific hybrids, providing valuable insights into erucic acid metabolism in Brassica.

Fine mapping of Rcr1 and analyses of its effect on transcriptome patterns during infection by Plasmodiophora brassicae.

BACKGROUND: The protist Plasmodiophora brassicae is a biotrophic soil-borne pathogen that causes clubroot on Brassica crops worldwide. Clubroot disease is a serious threat to the 8 M ha of canola (Brassica napus) grown annually in western Canada. While host resistance is the key to clubroot management, sources of resistance are limited. RESULTS: To identify new sources of clubroot resistance (CR), we fine mapped a CR gene (Rcr1) from B. rapa ssp. chinensis to the region between 24.26 Mb and 24.50 Mb on the linkage group A03, with several closely linked markers identified. Transcriptome analysis was conducted using RNA sequencing on a segregating F1 population inoculated with P. brassicae, with 2,212 differentially expressed genes (DEGs) identified between plants carrying and not carrying Rcr1. Functional annotation of these DEGs showed that several defense-related biological processes, including signaling and metabolism of jasmonate and ethylene, defensive deposition of callose and biosynthesis of indole-containing compounds, were up-regulated significantly in plants carrying Rcr1 while genes involved in salicylic acid metabolic and signaling pathways were generally not elevated. Several DEGs involved in metabolism potentially related to clubroot symptom development, including auxin biosynthesis and cell growth/development, showed significantly lower expression in plants carrying Rcr1. CONCLUSION: The CR gene Rcr1 and closely linked markers will be highly useful for breeding new resistant canola cultivars. The identification of DEGs between inoculated plants carrying and not carrying Rcr1 is an important step towards understanding of specific metabolic/signaling pathways in clubroot resistance mediated by Rcr1. This information may help judicious use of CR genes with complementary resistance mechanisms for durable clubroot resistance.

Analysis of genome-wide variants through bulked segregant RNA sequencing reveals a major gene for resistance to Plasmodiophora brassicae in Brassica oleracea.

Two cabbage (Brassica oleracea) cultivars 'Tekila' and 'Kilaherb' were identified as resistant to several pathotypes of Plasmodiophora brassicae. In this study, we identified a clubroot resistance gene (Rcr7) in 'Tekila' for resistance to pathotype 3 of P. brassicae from a segregating population derived from 'Tekila' crossed with the susceptible line T010000DH3. Genetic mapping was performed by identifying the percentage of polymorphic variants (PPV), a new method proposed in this study, through bulked segregant RNA sequencing. Chromosome C7 carried the highest PPV (42%) compared to the 30-34% in the remaining chromosomes. A peak with PPV (56-73%) was found within the physical interval 41-44 Mb, which indicated that Rcr7 might be located in this region. Kompetitive Allele-Specific PCR was used to confirm the association of Rcr7 with SNPs in the region. Rcr7 was flanked by two SNP markers and co-segregated with three SNP markers in the segregating population of 465 plants. Seven genes encoding TIR-NBS-LRR disease resistance proteins were identified in the target region, but only two genes, Bo7g108760 and Bo7g109000, were expressed. Resistance to pathotype 5X was also mapped to the same region as Rcr7. B. oleracea lines including 'Kilaherb' were tested with five SNP markers for Rcr7 and for resistance to pathotype 3; 11 of 25 lines were resistant, but 'Kilaherb' was the only line that carried the SNP alleles associated with Rcr7. The presence of Rcr7 in 'Kilaherb' for resistance to both pathotypes 3 and 5X was confirmed through linkage analysis.

Genotyping-by-sequencing reveals three QTL for clubroot resistance to six pathotypes of Plasmodiophora brassicae in Brassica rapa.

Clubroot, caused by Plasmodiophora brassicae, is an important disease of Brassica crops worldwide. F1 progeny from the Brassica rapa lines T19 (resistant) × ACDC (susceptible) were backcrossed with ACDC, then self-pollinated to produce BC1S1 lines, From genotyping-by-sequencing (GBS) of the parental lines and BC1 plants, about 1.32 M sequences from T19 were aligned into the reference genome of B. rapa with 0.4-fold coverage, and 1.77 M sequences with 0.5-fold coverage in ACDC. The number of aligned short reads per plant in the BC1 ranged from 0.07 to 1.41 M sequences with 0.1-fold coverage. A total of 1584 high quality SNP loci were obtained, distributed on 10 chromosomes. A single co-localized QTL, designated as Rcr4 on chromosome A03, conferred resistance to pathotypes 2, 3, 5, 6 and 8. The peak was at SNP locus A03_23710236, where LOD values were 30.3 to 38.8, with phenotypic variation explained (PVE) of 85-95%. Two QTLs for resistance to a novel P. brassicae pathotype 5x, designated Rcr8 on chromosome A02 and Rcr9 on A08, were detected with 15.0 LOD and 15.8 LOD, and PVE of 36% and 39%, respectively. Bulked segregant analysis was performed to examine TIR-NBS-LRR proteins in the regions harboring the QTL.

Fine Mapping of a Clubroot Resistance Gene in Chinese Cabbage Using SNP Markers Identified from Bulked Segregant RNA Sequencing.

Clubroot, caused by Plasmodiophora brassicae, is an important disease of canola (Brassica napus) in western Canada and worldwide. In this study, a clubroot resistance gene (Rcr2) was identified and fine mapped in Chinese cabbage cv. "Jazz" using single-nucleotide polymorphisms (SNP) markers identified from bulked segregant RNA sequencing (BSR-Seq) and molecular markers were developed for use in marker assisted selection. In total, 203.9 million raw reads were generated from one pooled resistant (R) and one pooled susceptible (S) sample, and >173,000 polymorphic SNP sites were identified between the R and S samples. One significant peak was observed between 22 and 26 Mb of chromosome A03, which had been predicted by BSR-Seq to contain the causal gene Rcr2. There were 490 polymorphic SNP sites identified in the region. A segregating population consisting of 675 plants was analyzed with 15 SNP sites in the region using the Kompetitive Allele Specific PCR method, and Rcr2 was fine mapped between two SNP markers, SNP_A03_32 and SNP_A03_67 with 0.1 and 0.3 cM from Rcr2, respectively. Five SNP markers co-segregated with Rcr2 in this region. Variants were identified in 14 of 36 genes annotated in the Rcr2 target region. The numbers of poly variants differed among the genes. Four genes encode TIR-NBS-LRR proteins and two of them Bra019410 and Bra019413, had high numbers of polymorphic variants and so are the most likely candidates of Rcr2.

Identification of Genome-Wide Variants and Discovery of Variants Associated with Brassica rapa Clubroot Resistance Gene Rcr1 through Bulked Segregant RNA Sequencing.

Clubroot, caused by Plasmodiophora brassicae, is an important disease on Brassica species worldwide. A clubroot resistance gene, Rcr1, with efficacy against pathotype 3 of P. brassicae, was previously mapped to chromosome A03 of B. rapa in pak choy cultivar "Flower Nabana". In the current study, resistance to pathotypes 2, 5 and 6 was shown to be associated with Rcr1 region on chromosome A03. Bulked segregant RNA sequencing was performed and short read sequences were assembled into 10 chromosomes of the B. rapa reference genome v1.5. For the resistant (R) bulks, a total of 351.8 million (M) sequences, 30,836.5 million bases (Mb) in length, produced 120-fold coverage of the reference genome. For the susceptible (S) bulks, 322.9 M sequences, 28,216.6 Mb in length, produced 109-fold coverage. In total, 776.2 K single nucleotide polymorphisms (SNPs) and 122.2 K insertion / deletion (InDels) in R bulks and 762.8 K SNPs and 118.7 K InDels in S bulks were identified; each chromosome had about 87% SNPs and 13% InDels, with 78% monomorphic and 22% polymorphic variants between the R and S bulks. Polymorphic variants on each chromosome were usually below 23%, but made up 34% of the variants on chromosome A03. There were 35 genes annotated in the Rcr1 target region and variants were identified in 21 genes. The numbers of poly variants differed significantly among the genes. Four out of them encode Toll-Interleukin-1 receptor / nucleotide-binding site / leucine-rich-repeat proteins; Bra019409 and Bra019410 harbored the higher numbers of polymorphic variants, which indicates that they are more likely candidates of Rcr1. Fourteen SNP markers in the target region were genotyped using the Kompetitive Allele Specific PCR method and were confirmed to associate with Rcr1. Selected SNP markers were analyzed with 26 recombinants obtained from a segregating population consisting of 1587 plants, indicating that they were completely linked to Rcr1. Nine SNP markers were used for marker-assisted introgression of Rcr1 into B. napus canola from B. rapa, with 100% accuracy in this study.

Mid-infrared spectral characteristics of lipid molecular structures in Brassica carinata seeds: relationship to oil content, fatty acid and glucosinolate profiles, polyphenols, and condensed tannins.

The objectives of this study were to quantify lipid-related inherent molecular structures using a Fourier transform infrared spectroscopy (FT-IR) technique and determine their relationship to oil content, fatty acid and glucosinolate profile, total polyphenols, and condensed tannins in seeds from newly developed yellow-seeded and brown-seeded Brassica carinata lines. Canola seeds were used as a reference. The lipid-related molecular spectral band intensities were strongly correlated to the contents of oil, fatty acids, glucosinolates, and polyphenols. The regression equations gave relatively high predictive power for the estimation of oil (R² = 0.99); all measured fatty acids (R² > 0.80), except C14:0, C20:3n-3, C22:2n-9, and C22:2n-6; 3-butenyl, 2-OH-3-butenyl, 4-OH-3-CH3-indolyl, and total glucosinolates (R² > 0.686); and total polyphenols (R² = 0.935). However, further study is required to obtain predictive equations based on large numbers of samples from diverse sources to illustrate the general applicability of these regression equations.

Studies on Brassica carinata seed. 2. Carbohydrate molecular structure in relation to carbohydrate chemical profile, energy values, and biodegradation characteristics.

The objectives of this study were to investigate (1) the carbohydrate chemical profile, (2) the energy values, (3) the rumen neutral detergent fiber (NDF) degradation kinetics, (4) the carbohydrate-related functional group structural features using a Fourier transform infrared (FTIR) spectroscopic technique with attenuated total reflectance (ATR), and (5) the correlations between carbohydrate intrinsic structural features and nutritional profiles in three strains of Brassica carinata in yellow and brown seed coats, with comparison to canola seed as a reference. The results showed that yellow B. carinata strains 111000EM and AAC A100 were lower for contents of neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), and carbohydrate (CHO) and higher for contents of total digestible nutrients (TDN), energy values, and effective degradable NDF (EDNDF) than brown-seeded 110915EM. In comparison, brown canola seed (Brassica napus L.) had more fiber content and less EDNDF. Also, carinata strains showed significantly different IR intensities in structural carbohydrate (SCHO), cellulosic compounds (CELC), and total CHO profiles. These structural variations might be one of the possible reasons for various fiber profile and biodegradation characteristics for ruminants in oilseeds. However, multivariate analyses within carbohydrate regions indicated there were still some structural relationships among the four oilseed samples. Moreover, the correlation study showed that the changes of CELC and CHO peak intensities were highly related with some changes in CHO chemical profile, energy values, and in situ NDF degradation kinetics in B. carinata and canola seeds. Further study with a large sample size is still necessary to figure out whether CHO molecular spectral information could be used to predict nutrient values and biological behavior in oilseeds.

Studies on Brassica carinata seed. 1. Protein molecular structure in relation to protein nutritive values and metabolic characteristics.

The objectives of this study were to investigate (1) the protein chemical profile, (2) the protein subfractions partitioned by the Cornell Net Carbohydrate and Protein System (CNCPS), (3) the rumen crude protein (CP) degradation kinetics, (4) the protein supply predicted by the DVE/OEB system, (5) the protein structural features using a Fourier transform infrared (FTIR) spectroscopic technique with attenuated total reflectance (ATR), and (6) the correlations between protein intrinsic structural features and nutritional profiles in three strains of Brassica carinata in yellow and brown seed coats, with comparison to canola seed as a reference. The results showed that carinata seed strains were different in both nutritional values and IR absorbance within the protein spectral region (ca. 1720-1482 cm(-1)). The comparison between yellow and brown B. carinata seeds indicated that the former was lower in acid detergent insoluble crude protein (ADICP; P = 0.002) and undegradable protein fraction (PC; P = 0.002) and greater in the degradable (D) fraction (P = 0.004) and true absorbed protein in the small intestine (DVE; P = 0.02) as well as feed milk value (FMV; P = 0.02) than the latter. The brown canola seed (Brassica napus L.) was also not in full accordance with B. carinata seed on these parameters. The FTIR studies showed significant differences in protein amide II peak height, amide I peak area, and ß-sheet height among different B. carinata strains. However, multivariate spectral analyses indicated a similarity in protein structural makeup in these four kinds of oilseed. The not very strong correlations shown in this study implied that the limited sample size and narrow range in biological and spectral variation might be responses for the weak relationships between chemical profile and mid-IR spectral data. Further studies using sufficient samples with wide and diverse range in nutritional properties are needed to illustrate the actual relationship between spectroscopic data and nutritional profiles in oilseeds.

Differential metabolite profiles and salinity tolerance between two genetically related brown-seeded and yellow-seeded Brassica carinata lines.

Brassica carinata (Ethiopian mustard) has previously been identified as a potential crop species suitable for marginal land in the North American prairies due to its relatively high salt tolerance. Two genetically related B. carinata lines with brown-seeded (BS) and yellow-seeded (YS) phenotypes were assessed for their tolerance to sodium sulfate. Specifically, each line was greenhouse-grown under 0, 50 and 100mM of salt, and analyzed after four weeks and eight weeks of treatment. Generally, the height of the BS line was greater than the YS line under both salt treatments, indicating enhanced salt tolerance of the BS line. NMR-based metabolite profiling and PCA analyses indicated a more pronounced shift in key stem metabolites after four weeks of treatment with the YS line compared to the BS line. For example, tryptophan and formate levels increased in the YS line after four weeks of 100mM salt treatment, while proline and threonine levels varied uniquely compared to other metabolites of the lines. Together, the data indicate that the brown-seeded line has greater sodium tolerance than the yellow-seed line, provide clues to the biochemical underpinnings for the phenotypic variation, and highlight the utility of B. carinata as a biorefinery crop for saline land.