Genomics of predictive radiation mutagenesis in oilseed rape: modifying seed oil composition

Rapeseed is a crop of global importance but there is a need to broaden the genetic diversity available to address breeding objectives. Radiation mutagenesis, supported by genomics, has the potential to supersede genome editing for both gene knockout and copy number increase, but detailed knowledge of the molecular outcomes of radiation treatment is lacking. To address this, we produced a genome re-sequenced panel of 1133 M2 generation rapeseed plants and analysed large-scale deletions, single nucleotide variants and small insertion-deletion variants affecting gene open reading frames. We show that high radiation doses (2000 Gy) are tolerated, gamma radiation and fast neutron radiation have similar impacts and that segments deleted from the genomes of some plants are inherited as additional copies by their siblings, enabling gene dosage decrease. Of relevance for species with larger genomes, we showed that these large-scale impacts can also be detected using transcriptome re-sequencing. To test the utility of the approach for predictive alteration of oil fatty acid composition, we produced lines with both decreased and increased copy numbers of Bna.FAE1 and confirmed the anticipated impacts on erucic acid content. We detected and tested a 21-base deletion expected to abolish function of Bna.FAD2.A5, for which we confirmed the predicted reduction in seed oil polyunsaturated fatty acid content. Our improved understanding of the molecular effects of radiation mutagenesis will underpin genomics-led approaches to more efficient introduction of novel genetic variation into the breeding of this crop and provides an exemplar for the predictive improvement of other crops.

Co-expression network analysis of diverse wheat landraces reveals markers of early thermotolerance and a candidate master regulator of thermotolerance genes

Triticum aestivum L. (bread wheat) is a crop relied upon by billions of people around the world, as a major source of both income and calories. Rising global temperatures, however, pose a genuine threat to the livelihood of these people, as wheat growth and yields are extremely vulnerable to damage by heat stress. Here we present the YoGI wheat landrace panel, comprising 342 accessions that show remarkable phenotypic and genetic diversity thanks to their adaptation to different climates. We quantified the abundance of 110 790 transcripts from the panel and used these data to conduct weighted co-expression network analysis and to identify hub genes in modules associated with abiotic stress tolerance. We found that the expression of three hub genes, all heat-shock proteins (HSPs), were significantly correlated with early thermotolerance in a validation panel of landraces. These hub genes belong to the same module, with one (TraesCS4D01G207500.1) being a candidate master-regulator potentially controlling the expression of the other two hub genes, as well as a suite of other HSPs and heat-stress transcription factors (HSFs). In this work, therefore, we identify three validated hub genes, the expression of which can serve as markers of thermotolerance during early development, and suggest that TraesCS4D01G207500.1 is a potential master regulator of HSP and HSF expression - presenting the YoGI landrace panel as an invaluable tool for breeders wishing to determine and introduce novel alleles into modern varieties, for the production of climate-resilient crops.

Natural Variation in OASC Gene for Mitochondrial O-Acetylserine Thiollyase Affects Sulfate Levels in Arabidopsis

Sulfur plays a vital role in the primary and secondary metabolism of plants, and carries an important function in a large number of different compounds. Despite this importance, compared to other mineral nutrients, relatively little is known about sulfur sensing and signalling, as well as about the mechanisms controlling sulfur metabolism and homeostasis. Sulfur contents in plants vary largely not only among different species, but also among accessions of the same species. We previously used associative transcriptomics to identify several genes potentially controlling variation in sulfate content in the leaves of Brassica napus, including an OASC gene for mitochondrial O-acetylserine thiollyase (OAS-TL), an enzyme involved in cysteine synthesis. Here, we show that loss of OASC in Arabidopsis thaliana lowers not only sulfate, but also glutathione levels in the leaves. The reduced accumulation is caused by lower sulfate uptake and translocation to the shoots; however, the flux through the pathway is not affected. In addition, we identified a single nucleotide polymorphism in the OASC gene among A. thaliana accessions that is linked to variation in sulfate content. Both genetic and transgenic complementation confirmed that the exchange of arginine at position 81 for lysine in numerous accessions resulted in a less active OASC and a lower sulfate content in the leaves. The mitochondrial isoform of OAS-TL is, thus, after the ATPS1 isoform of sulfurylase and the APR2 form of APS reductase 2, the next metabolic enzyme with a role in regulation of sulfate content in Arabidopsis.

Genomic signatures of vegetable and oilseed allopolyploid Brassica juncea and genetic loci controlling the accumulation of glucosinolates

Allopolyploid Brassica juncea crops in Brassicaceae are becoming increasingly revitalized as vegetables and oilseeds owing to wide adaptability and significant economic values. However, the genomic differentiation of diversified vegetables and oilseed B. juncea and the genetic basis underlying glucosinolates accumulation have yet to be elucidated. To address this knowledge gap, we report the sequencing of pairwise genomes of vegetable and oilseed B. juncea at chromosome scale. Comparative genomics analysis unveils panoramic structural variation footprints, particularly the genetic loci of HSP20 and TGA1 associated with abiotic and biotic stresses responses between oilseed and vegetable subgroups. We anchored two major loci of MYB28 (HAG1) orthologues caused by copy number variations on A02 and A09 chromosomes using scored genomic SNPs-based GWAS that are responsible for seed oil quality-determining glucosinolates biosynthesis. These findings will provide valuable repertories of polyploidy genomic information enabling polyploidy genome evolution studies and precise genomic selections for crucial traits like functional components of glucosinolates in B. juncea crops and beyond.

Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids

Plant genomes demonstrate significant presence/absence variation (PAV) within a species; however, the factors that lead to this variation have not been studied systematically in Brassica across diploids and polyploids. Here, we developed pangenomes of polyploid Brassica napus and its two diploid progenitor genomes B. rapa and B. oleracea to infer how PAV may differ between diploids and polyploids. Modelling of gene loss suggests that loss propensity is primarily associated with transposable elements in the diploids while in B. napus, gene loss propensity is associated with homoeologous recombination. We use these results to gain insights into the different causes of gene loss, both in diploids and following polyploidization, and pave the way for the application of machine learning methods to understanding the underlying biological and physical causes of gene presence/absence.

Validation of a novel associative transcriptomics pipeline in Brassica oleracea: identifying candidates for vernalisation response

BACKGROUND

Associative transcriptomics has been used extensively in Brassica napus to enable the rapid identification of markers correlated with traits of interest. However, within the important vegetable crop species, Brassica oleracea, the use of associative transcriptomics has been limited due to a lack of fixed genetic resources and the difficulties in generating material due to self-incompatibility. Within Brassica vegetables, the harvestable product can be vegetative or floral tissues and therefore synchronisation of the floral transition is an important goal for growers and breeders. Vernalisation is known to be a key determinant of the floral transition, yet how different vernalisation treatments influence flowering in B. oleracea is not well understood.

RESULTS

Here, we present results from phenotyping a diverse set of 69 B. oleracea accessions for heading and flowering traits under different environmental conditions. We developed a new associative transcriptomics pipeline, and inferred and validated a population structure, for the phenotyped accessions. A genome-wide association study identified miR172D as a candidate for the vernalisation response. Gene expression marker association identified variation in expression of BoFLC.C2 as a further candidate for vernalisation response.

CONCLUSIONS

This study describes a new pipeline for performing associative transcriptomics studies in B. oleracea. Using flowering time as an example trait, it provides insights into the genetic basis of vernalisation response in B. oleracea through associative transcriptomics and confirms its characterisation as a complex G x E trait. Candidate leads were identified in miR172D and BoFLC.C2. These results could facilitate marker-based breeding efforts to produce B. oleracea lines with more synchronous heading dates, potentially leading to improved yields.

Genome structural evolution in Brassica crops

The cultivated Brassica species include numerous vegetable and oil crops of global importance. Three genomes (designated A, B and C) share mesohexapolyploid ancestry and occur both singly and in each pairwise combination to define the Brassica species. With organizational errors (such as misplaced genome segments) corrected, we showed that the fundamental structure of each of the genomes is the same, irrespective of the species in which it occurs. This enabled us to clarify genome evolutionary pathways, including updating the Ancestral Crucifer Karyotype (ACK) block organization and providing support for the Brassica mesohexaploidy having occurred via a two-step process. We then constructed genus-wide pan-genomes, drawing from genes present in any species in which the respective genome occurs, which enabled us to provide a global gene nomenclature system for the cultivated Brassica species and develop a methodology to cost-effectively elucidate the genomic impacts of alien introgressions. Our advances not only underpin knowledge-based approaches to the more efficient breeding of Brassica crops but also provide an exemplar for the study of other polyploids.

QTL-seq identifies BnaFT.A02 and BnaFLC.A02 as candidates for variation in vernalization requirement and response in winter oilseed rape (Brassica napus)

Winter, spring and biennial varieties of Brassica napus that vary in vernalization requirement are grown for vegetable and oil production. Here, we show that the obligate or facultative nature of the vernalization requirement in European winter oilseed rape is determined by allelic variation at a 10 Mbp region on chromosome A02. This region includes orthologues of the key floral regulators FLOWERING LOCUS C (BnaFLC.A02) and FLOWERING LOCUS T (BnaFT.A02). Polymorphism at BnaFLC.A02 and BnaFT.A02, mostly in cis-regulatory regions, results in distinct gene expression dynamics in response to vernalization treatment. Our data suggest allelic variation at BnaFT.A02 is associated with flowering time in the absence of vernalization, while variation at BnaFLC.A02 is associated with flowering time under vernalizing conditions. We hypothesize selection for BnaFLC.A02 and BnaFT.A02 gene expression variation has facilitated the generation of European winter oilseed rape varieties that are adapted to different winter climates. This knowledge will allow for the selection of alleles of flowering time regulators that alter the vernalization requirement of oilseed rape, informing the generation of new varieties with adapted flowering times and improved yields.

The WRKY6 transcription factor affects seed oil accumulation and alters fatty acid compositions in Arabidopsis thaliana

In rapeseed, the oil content of the seed not only supplies energy for seed germination and seedling development but also provides essential dietary nutrients for humans and livestock. Recent studies have revealed that many transcription factors (TFs) regulate the accumulation of fatty acids (FAs) during seed development. WRKY6, a WRKY6 family TF, was reported to serve a function in the plant senescence processes, pathogen defense mechanisms and abiotic stress responses. However, the precise role of WRKY6 in influencing FA accumulation in seeds is still unknown. In this study, we demonstrate that WRKY6 has a high expression level in developing seeds and plays an essential role in regulating the accumulation of FAs in developing seeds of Arabidopsis. Mutation of WRKY6 resulted in significant increase in seed size, accompanied by an increase in FA content and changes in FA composition. Ultrastructure analyses showed that the absence of WRKY6 resulted in more and higher percentage of oil body in the cell of mature seeds. Quantitative real-time PCR analysis revealed changes in the expression of several genes related to photosynthesis and FA biosynthesis in wrky6 mutants at 10 or 16 days after pollination. These results reveal a novel function of WRKY6 influencing seed oil content and FAs compositions. This gene could be used as a promising gene resource to improve FA accumulation and seed yield in Brassica napus through genetic manipulation.

Validation of an Associative Transcriptomics platform in the polyploid crop species Brassica juncea by dissection of the genetic architecture of agronomic and quality traits

The development of more productive crops will be key to addressing the challenges that climate change, population growth and diminishing resources pose to global food security. Advanced 'omics techniques can help to accelerate breeding by facilitating the identification of genetic markers for use in marker-assisted selection. Here, we present the validation of a new Associative Transcriptomics platform in the important oilseed crop Brassica juncea. To develop this platform, we established a pan-transcriptome reference for B. juncea, to which we mapped transcriptome data from a diverse panel of B. juncea accessions. From this panel, we identified 355 050 single nucleotide polymorphism variants and quantified the abundance of 93 963 transcripts. Subsequent association analysis of functional genotypes against a number of important agronomic and quality traits revealed a promising candidate gene for seed weight, BjA.TTL, as well as additional markers linked to seed colour and vitamin E content. The establishment of the first full-scale Associative Transcriptomics platform for B. juncea enables rapid progress to be made towards an understanding of the genetic architecture of trait variation in this important species, and provides an exemplar for other crops.

A high-contiguity Brassica nigra genome localizes active centromeres and defines the ancestral Brassica genome

It is only recently, with the advent of long-read sequencing technologies, that we are beginning to uncover previously uncharted regions of complex and inherently recursive plant genomes. To comprehensively study and exploit the genome of the neglected oilseed Brassica nigra, we generated two high-quality nanopore de novo genome assemblies. The N50 contig lengths for the two assemblies were 17.1 Mb (12 contigs), one of the best among 324 sequenced plant genomes, and 0.29 Mb (424 contigs), respectively, reflecting recent improvements in the technology. Comparison with a de novo short-read assembly corroborated genome integrity and quantified sequence-related error rates (0.2%). The contiguity and coverage allowed unprecedented access to low-complexity regions of the genome. Pericentromeric regions and coincidence of hypomethylation enabled localization of active centromeres and identified centromere-associated ALE family retro-elements that appear to have proliferated through relatively recent nested transposition events (<1 Ma). Genomic distances calculated based on synteny relationships were used to define a post-triplication Brassica-specific ancestral genome, and to calculate the extensive rearrangements that define the evolutionary distance separating B. nigra from its diploid relatives.

A Brassica napus Reductase Gene Dissected by Associative Transcriptomics Enhances Plant Adaption to Freezing Stress

Cold treatment (vernalization) is required for winter crops such as rapeseed (Brassica napus L.). However, excessive exposure to low temperature (LT) in winter is also a stress for the semi-winter, early-flowering rapeseed varieties widely cultivated in China. Photosynthetic efficiency is one of the key determinants, and thus a good indicator for LT tolerance in plants. So far, the genetic basis underlying photosynthetic efficiency is poorly understood in rapeseed. Here the current study used Associative Transcriptomics to identify genetic loci controlling photosynthetic gas exchange parameters in a diversity panel comprising 123 accessions. A total of 201 significant Single Nucleotide Polymorphisms (SNPs) and 147 Gene Expression Markers (GEMs) were detected, leading to the identification of 22 candidate genes. Of these, Cab026133.1, an ortholog of the Arabidopsis gene AT2G29300.2 encoding a tropinone reductase (BnTR1), was further confirmed to be closely linked to transpiration rate. Ectopic expressing BnTR1 in Arabidopsis plants significantly increased the transpiration rate and enhanced LT tolerance under freezing conditions. Also, a much higher level of alkaloids content was observed in the transgenic Arabidopsis plants, which could help protect against LT stress. Together, the current study showed that AT is an effective approach for dissecting LT tolerance trait in rapeseed and that BnTR1 is a good target gene for the genetic improvement of LT tolerance in plant.

The impact of reducing fatty acid desaturation on the composition and thermal stability of rapeseed oil

Oilseed rape (Brassica napus) is the third largest source of vegetable oil globally. In addition to food uses, there are industrial applications that exploit the ability of the species to accumulate the very-long-chain fatty acid (VLCFA) erucic acid in its seed oil, controlled by orthologues of FATTY ACID ELONGASE 1 (Bna.FAE1.A8 and Bna.FAE1.C3). The proportion of polyunsaturated fatty acids (PUFAs) in rapeseed oil is predicted to affect its thermal stability and is controlled by orthologues of FATTY ACID DESATURASE 2, particularly Bna.FAD2.C5. Our aim was to develop rapeseed lines combining high erucic and low PUFA characters and to assess the impact on thermal stability of the oil they produce. The new type of rapeseed oil (high erucic low polyunsaturate; HELP) contained a substantially greater proportion of erucic acid (54%) compared with high erucic rapeseed oil (46%). Although the total VLCFA content was greater in oil from HELP lines (64%) than from high erucic rapeseed (57%), analysis of triacylglycerol composition showed negligible incorporation of VLCFAs into the sn-2 position. Rancimat analysis showed that the thermal stability of rapeseed oil was improved greatly as a consequence of reduction of PUFA content, from 3.8 and 4.2 h in conventional low erucic and high erucic rapeseed oils, respectively, to 11.3 and 16.4 h in high oleic low PUFA (HOLP) and HELP oils, respectively. Our results demonstrate that engineering of the lipid biosynthetic pathway of rapeseed, using traditional approaches, enables the production of renewable industrial oils with novel composition and properties.

Genome-wide selection footprints and deleterious variations in young Asian allotetraploid rapeseed

Brassica napus (AACC, 2n = 38) is an important oilseed crop grown worldwide. However, little is known about the population evolution of this species, the genomic difference between its major genetic groups, such as European and Asian rapeseed, and the impacts of historical large-scale introgression events on this young tetraploid. In this study, we reported the de novo assembly of the genome sequences of an Asian rapeseed (B. napus), Ningyou 7, and its four progenitors and compared these genomes with other available genomic data from diverse European and Asian cultivars. Our results showed that Asian rapeseed originally derived from European rapeseed but subsequently significantly diverged, with rapid genome differentiation after hybridization and intensive local selective breeding. The first historical introgression of B. rapa dramatically broadened the allelic pool but decreased the deleterious variations of Asian rapeseed. The second historical introgression of the double-low traits of European rapeseed (canola) has reshaped Asian rapeseed into two groups (double-low and double-high), accompanied by an increase in genetic load in the double-low group. This study demonstrates distinctive genomic footprints and deleterious SNP (single nucleotide polymorphism) variants for local adaptation by recent intra- and interspecies introgression events and provides novel insights for understanding the rapid genome evolution of a young allopolyploid crop.

Genetic architecture of glucosinolate variation in Brassica napus

The diverse biological activities of glucosinolate (GSL) hydrolysis products play significant biological and economical roles in the defense system and nutritional qualities of Brassica napus (oilseed rape). Yet, genomic-based study of the B. napus GSL regulatory mechanisms are scarce due to the complexity of working with polyploid species. To address these challenges, we used transcriptome-based GWAS approach, Associative Transcriptomics (AT), across a diversity panel of 288 B. napus genotypes to uncover the underlying genetic basis controlling quantitative variation of GSLs in B. napus vegetative tissues. Single nucleotide polymorphism (SNP) markers and gene expression markers (GEMs) associations identify orthologues of MYB28/HAG1 (AT5G61420), specifically the copies on chromosome A9 and C2, to be the key regulators of aliphatic GSL variation in leaves. We show that the positive correlation observed between aliphatic GSLs in seed and leaf is due to the amount synthesized, as controlled by Bna.HAG1.A9 and Bna.HAG1.C2, rather than by variation in the transport processes. In addition, AT and differential expression analysis in root tissues implicate an orthologue of MYB29/HAG3 (AT5G07690), Bna.HAG3.A3, as controlling root aromatic GSL variation. Based on the root expression data we also propose Bna.MAM3.A3 to have a role in controlling phenylalanine chain elongation for aromatic GSL biosynthesis. This work uncovers a regulator of homophenylalanine-derived aromatic GSLs and implicates the shared biosynthetic pathways between aliphatic and aromatic GSLs.

Data in support of genetic architecture of glucosinolate variations in Brassica napus

The transcriptome-based GWAS approach, Associative Transcriptomics (AT), which was employed to uncover the genetic basis controlling quantitative variation of glucosinolates in Brassica napus vegetative tissues is described. This article includes the phenotypic data of leaf and root glucosinolate (GSL) profiles across a diversity panel of 288 B. napus genotypes, as well as information on population structure and levels of GSLs grouped by crop types. Moreover, data on genetic associations of single nucleotide polymorphism (SNP) markers and gene expression markers (GEMs) for the major GSL types are presented in detail, while Manhattan plots and QQ plots for the associations of individual GSLs are also included. Root genetic association are supported by differential expression analysis generated from root RNA-seq. For further interpretation and details, please see the related research article entitled 'Genetic architecture of glucosinolate variation in Brassica napus' (Kittipol et al., 2019).

Analysing the genetic architecture of clubroot resistance variation in Brassica napus by associative transcriptomics

Clubroot is a destructive soil-borne pathogen of Brassicaceae that causes significant recurrent reductions in yield of cruciferous crops. Although there is some resistance in oilseed rape (a crop type of the species Brassica napus), the genetic basis of that resistance is poorly understood. In this study, we used an associative transcriptomics approach to elucidate the genetic basis of resistance to clubroot pathotype ECD 17/31/31 across a genetic diversity panel of 245 accessions of B. napus. A single nucleotide polymorphism (SNP) association analysis was performed with 256,397 SNPs distributed across the genome of B. napus and combined with transcript abundance data of 53,889 coding DNA sequence (CDS) gene models. The SNP association analysis identified two major loci (on chromosomes A2 and A3) controlling resistance and seven minor loci. Within these were a total of 86 SNP markers. Altogether, 392 genes were found in these regions. Another 21 genes were implicated as potentially involved in resistance using gene expression marker (GEM) analysis. After GO enrichment analysis and InterPro functional analysis of the identified genes, 82 candidate genes were identified as having roles in clubroot resistance. These results provide useful information for marker-assisted breeding which could lead to acceleration of pyramiding of multiple clubroot resistance genes in new varieties.

Cytonuclear interactions remain stable during allopolyploid evolution despite repeated whole-genome duplications in Brassica

Several plastid macromolecular protein complexes are encoded by both nuclear and plastid genes. Therefore, cytonuclear interactions are held in place to prevent genomic conflicts that may lead to incompatibilities. Allopolyploidy resulting from hybridization and genome doubling of two divergent species can disrupt these fine-tuned interactions, as newly formed allopolyploid species confront biparental nuclear chromosomes with a uniparentally inherited plastid genome. To avoid any deleterious effects of unequal genome inheritance, preferential transcription of the plastid donor over the other donor has been hypothesized to occur in allopolyploids. We used Brassica as a model to study the effects of paleopolyploidy in diploid parental species, as well as the effects of recent and ancient allopolyploidy in Brassica napus, on genes implicated in plastid protein complexes. We first identified redundant nuclear copies involved in those complexes. Compared with cytosolic protein complexes and with genome-wide retention rates, genes involved in plastid protein complexes show a higher retention of genes in duplicated and triplicated copies. Those redundant copies are functional and are undergoing strong purifying selection. We then compared transcription patterns and sequences of those redundant gene copies between resynthesized allopolyploids and their diploid parents. The neopolyploids showed no biased subgenome expression or maternal homogenization via gene conversion, despite the presence of some non-synonymous substitutions between plastid genomes of parental progenitors. Instead, subgenome dominance was observed regardless of the maternal progenitor. Our results provide new insights on the evolution of plastid protein complexes that could be tested and generalized in other allopolyploid species.

Species-Wide Variation in Shoot Nitrate Concentration, and Genetic Loci Controlling Nitrate, Phosphorus and Potassium Accumulation in Brassica napus L

Large nitrogen, phosphorus and potassium fertilizer inputs are used in many crop systems. Identifying genetic loci controlling nutrient accumulation may be useful in crop breeding strategies to increase fertilizer use efficiency and reduce financial and environmental costs. Here, variation in leaf nitrate concentration across a diversity population of 383 genotypes of Brassica napus was characterized. Genetic loci controlling variation in leaf nitrate, phosphorus and potassium concentration were then identified through Associative Transcriptomics using single nucleotide polymorphism (SNP) markers and gene expression markers (GEMs). Leaf nitrate concentration varied over 8-fold across the diversity population. A total of 455 SNP markers were associated with leaf nitrate concentration after false-discovery-rate (FDR) correction. In linkage disequilibrium of highly associated markers are a number of known nitrate transporters and sensors, including a gene thought to mediate expression of the major nitrate transporter NRT1.1. Several genes influencing root and root-hair development co-localize with chromosomal regions associated with leaf P concentration. Orthologs of three ABC-transporters involved in suberin synthesis in roots also co-localize with association peaks for both leaf nitrate and phosphorus. Allelic variation at nearby, highly associated SNPs confers large variation in leaf nitrate and phosphorus concentration. A total of five GEMs associated with leaf K concentration after FDR correction including a GEM that corresponds to an auxin-response family protein. Candidate loci, genes and favorable alleles identified here may prove useful in marker-assisted selection strategies to improve fertilizer use efficiency in B. napus.

Validation of an updated Associative Transcriptomics platform for the polyploid crop species Brassica napus by dissection of the genetic architecture of erucic acid and tocopherol isoform variation in seeds

An updated platform was developed to underpin association genetics studies in the polyploid crop species Brassica napus (oilseed rape). Based on 1.92 × 10¹² bases of leaf mRNAseq data, functional genotypes, comprising 355 536 single-nucleotide polymorphism markers and transcript abundance were scored across a genetic diversity panel of 383 accessions using a transcriptome reference comprising 116 098 ordered coding DNA sequence (CDS) gene models. The use of the platform for Associative Transcriptomics was first tested by analysing the genetic architecture of variation in seed erucic acid content, as high-erucic rapeseed oil is highly valued for a variety of applications in industry. Known loci were identified, along with a previously undetected minor-effect locus. The platform was then used to analyse variation for the relative proportions of tocopherol (vitamin E) forms in seeds, and the validity of the most significant markers was assessed using a take-one-out approach. Furthermore, the analysis implicated expression variation of the gene Bo2g050970.1, an orthologue of VTE4 (which encodes a γ-tocopherol methyl transferase converting γ-tocopherol into α-tocopherol) associated with the observed trait variation. The establishment of the first full-scale Associative Transcriptomics platform for B. napus enables rapid progress to be made towards an understanding of the genetic architecture of trait variation in this important species, and provides an exemplar for other crops.

Dissecting the complex regulation of lodging resistance in Brassica napus

Lodging continues to be a major cause of yield loss in important crop species such as Brassica napus. Understanding the genetic regulation of lodging resistance is therefore of key interest to breeders worldwide. Current strategies aimed at minimising lodging risk involve the incorporation of dwarfing genes or the application of plant growth regulators. However, despite these efforts, lodging continues to be a persistent problem and it is therefore of high interest that novel, complimentary strategies for lodging control are implemented. One approach would be to focus on understanding the genetic properties underlying stem mechanical strength. With this in mind, we screened a training genetic diversity panel of B. napus accession for variation in stem mechanical strength and related traits. Using Associative Transcriptomics, we identified molecular markers for a suite of valuable traits. Using an independent test genetic diversity panel, we show that the methods employed are robust for identification of predictive markers. Furthermore, based on conserved synteny with Arabidopsis thaliana, we are able to provide a biological context to the marker associations detected and provide evidence for a role in pectin methylesterification in contributing to stem mechanical strength in Brassicaceae.

Assembly and comparison of two closely related Brassica napus genomes

As an increasing number of plant genome sequences become available, it is clear that gene content varies between individuals, and the challenge arises to predict the gene content of a species. However, genome comparison is often confounded by variation in assembly and annotation. Differentiating between true gene absence and variation in assembly or annotation is essential for the accurate identification of conserved and variable genes in a species. Here, we present the de novo assembly of the B. napus cultivar Tapidor and comparison with an improved assembly of the Brassica napus cultivar Darmor-bzh. Both cultivars were annotated using the same method to allow comparison of gene content. We identified genes unique to each cultivar and differentiate these from artefacts due to variation in the assembly and annotation. We demonstrate that using a common annotation pipeline can result in different gene predictions, even for closely related cultivars, and repeat regions which collapse during assembly impact whole genome comparison. After accounting for differences in assembly and annotation, we demonstrate that the genome of Darmor-bzh contains a greater number of genes than the genome of Tapidor. Our results are the first step towards comparison of the true differences between B. napus genomes and highlight the potential sources of error in future production of a B. napus pangenome.

Variation across a wheat genetic diversity panel for saccharification of hydrothermally pretreated straw

BACKGROUND

Wheat straw forms an important, reliable source of lignocellulosic biomass for use in second-generation ethanol production. However, there is limited understanding of the variation in quality of straw from current breeding cultivars, and studies on such variation have generally employed suboptimal pretreatments. There is also a degree of confusion regarding phenotypic characteristics relevant to optimising the enzymatic saccharification of cellulose after suitable pretreatments for biorefining compared with those which determine good ruminant digestibility. The aim of this study has been to (a) evaluate and compare the levels of glucose enzymatically released from straw obtained from 89 cultivars of winter wheat after optimised hydrothermal pretreatments and (b) identify the underlying phenotypic characteristics relevant to enhanced glucose production with special reference to the ratios of constituent tissue types.

RESULTS

Optimised pretreatment involved hydrothermal extraction at 210 °C for 10 min. Using excess cellulases, quantitative saccharification was achieved within 24 h. The amount of glucose released ranged from 192 to 275 mg/g. The extent of glucose release was correlated with (a) the level of internode tissue (R = 0.498; p = 6.84 × 10^-7), (b) stem height (R = 0.491; p = 1.03 × 10^-6), and (c) chemical characteristics particular to stem tissues including higher levels of cellulose (R = 0.552; p = 2.06 × 10^-8) and higher levels of lignin R = 0.494; p = 8.67 × 10^-7.

CONCLUSIONS

In order to achieve maximum yields of cellulosic glucose for second-generation ethanol production, a predisposition for wheat to produce cellulose-enriched internode stem tissue, particularly of longer length, would be beneficial. This contrasts with the ideotype for ruminant nutrition, in which an increased proportion of leaf tissue is preferable.

Extensive homoeologous genome exchanges in allopolyploid crops revealed by mRNAseq-based visualization

Polyploidy, the possession of multiple sets of chromosomes, has been a predominant factor in the evolution and success of the angiosperms. Although artificially formed allopolyploids show a high rate of genome rearrangement, the genomes of cultivars and germplasm used for crop breeding were assumed stable and genome structural variation under the artificial selection process of commercial breeding has remained little studied. Here, we show, using a repurposed visualization method based on transcriptome sequence data, that genome structural rearrangement occurs frequently in varieties of three polyploid crops (oilseed rape, mustard rape and bread wheat), meaning that the extent of genome structural variation present in commercial crops is much higher than expected. Exchanges were found to occur most frequently where homoeologous chromosome segments are collinear to telomeres and in material produced as doubled haploids. The new insights into genome structural evolution enable us to reinterpret the results of recent studies and implicate homoeologous exchanges, not deletions, as being responsible for variation controlling important seed quality traits in rapeseed. Having begun to identify the extent of genome structural variation in polyploid crops, we can envisage new strategies for the global challenge of broadening crop genetic diversity and accelerating adaptation, such as the molecular identification and selection of genome deletions or duplications encompassing genes with trait-controlling dosage effects.

Development of an efficient glucosinolate extraction method

BACKGROUND

Glucosinolates, anionic sulfur rich secondary metabolites, have been extensively studied because of their occurrence in the agriculturally important brassicaceae and their impact on human and animal health. There is also increasing interest in the biofumigant properties of toxic glucosinolate hydrolysis products as a method to control agricultural pests. Evaluating biofumigation potential requires rapid and accurate quantification of glucosinolates, but current commonly used methods of extraction prior to analysis involve a number of time consuming and hazardous steps; this study aimed to develop an improved method for glucosinolate extraction.

RESULTS

Three methods previously used to extract glucosinolates from brassicaceae tissues, namely extraction in cold methanol, extraction in boiling methanol, and extraction in boiling water were compared across tissue type (root, stem leaf) and four brassicaceae species (B. juncea, S. alba, R. sativus, and E. sativa). Cold methanol extraction was shown to perform as well or better than all other tested methods for extraction of glucosinolates with the exception of glucoraphasatin in R. sativus shoots. It was also demonstrated that lyophilisation methods, routinely used during extraction to allow tissue disruption, can reduce final glucosinolate concentrations and that extracting from frozen wet tissue samples in cold 80% methanol is more effective.

CONCLUSIONS

We present a simplified method for extracting glucosinolates from plant tissues which does not require the use of a freeze drier or boiling methanol, and is therefore less hazardous, and more time and cost effective. The presented method has been shown to have comparable or improved glucosinolate extraction efficiency relative to the commonly used ISO method for major glucosinolates in the Brassicaceae species studied: sinigrin and gluconasturtiin in B. juncea; sinalbin, glucotropaeolin, and gluconasturtiin in S. alba; glucoraphenin and glucoraphasatin in R. sativus; and glucosatavin, glucoerucin and glucoraphanin in E. sativa.

Development of an efficient glucosinolate extraction method

BACKGROUND

Glucosinolates, anionic sulfur rich secondary metabolites, have been extensively studied because of their occurrence in the agriculturally important brassicaceae and their impact on human and animal health. There is also increasing interest in the biofumigant properties of toxic glucosinolate hydrolysis products as a method to control agricultural pests. Evaluating biofumigation potential requires rapid and accurate quantification of glucosinolates, but current commonly used methods of extraction prior to analysis involve a number of time consuming and hazardous steps; this study aimed to develop an improved method for glucosinolate extraction.

RESULTS

Three methods previously used to extract glucosinolates from brassicaceae tissues, namely extraction in cold methanol, extraction in boiling methanol, and extraction in boiling water were compared across tissue type (root, stem leaf) and four brassicaceae species (B. juncea, S. alba, R. sativus, and E. sativa). Cold methanol extraction was shown to perform as well or better than all other tested methods for extraction of glucosinolates with the exception of glucoraphasatin in R. sativus shoots. It was also demonstrated that lyophilisation methods, routinely used during extraction to allow tissue disruption, can reduce final glucosinolate concentrations and that extracting from frozen wet tissue samples in cold 80% methanol is more effective.

CONCLUSIONS

We present a simplified method for extracting glucosinolates from plant tissues which does not require the use of a freeze drier or boiling methanol, and is therefore less hazardous, and more time and cost effective. The presented method has been shown to have comparable or improved glucosinolate extraction efficiency relative to the commonly used ISO method for major glucosinolates in the Brassicaceae species studied: sinigrin and gluconasturtiin in B. juncea; sinalbin, glucotropaeolin, and gluconasturtiin in S. alba; glucoraphenin and glucoraphasatin in R. sativus; and glucosatavin, glucoerucin and glucoraphanin in E. sativa.

Identification of Candidate Genes for Calcium and Magnesium Accumulation in Brassica napus L. by Association Genetics

Calcium (Ca) and magnesium (Mg) are essential plant nutrients and vital for human and animal nutrition. Biofortification of crops has previously been suggested to alleviate widespread human Ca and Mg deficiencies. In this study, new candidate genes influencing the leaf accumulation of Ca and Mg were identified in young Brassica napus plants using associative transcriptomics of ionomics datasets. A total of 247 and 166 SNP markers were associated with leaf Ca and Mg concentration, respectively, after false discovery rate correction and removal of SNPs with low second allele frequency. Gene expression markers at similar positions were also associated with leaf Ca and Mg concentration, including loci on chromosomes A10 and C2, within which lie previously identified transporter genes ACA8 and MGT7. Further candidate genes were selected from seven loci and the mineral composition of whole Arabidopsis thaliana shoots were characterized from lines mutated in orthologous genes. Four and two mutant lines had reduced shoot Ca and Mg concentration, respectively, compared to wild type plants. Three of these mutations were found to have tissue specific effects; notably reduced silique Ca in all three such mutant lines. This knowledge could be applied in targeted breeding, with the possibility of increasing Ca and Mg in plant tissue for improving human and livestock nutrition.

Root morphology and seed and leaf ionomic traits in a Brassica napus L. diversity panel show wide phenotypic variation and are characteristic of crop habit

BACKGROUND

Mineral nutrient uptake and utilisation by plants are controlled by many traits relating to root morphology, ion transport, sequestration and translocation. The aims of this study were to determine the phenotypic diversity in root morphology and leaf and seed mineral composition of a polyploid crop species, Brassica napus L., and how these traits relate to crop habit. Traits were quantified in a diversity panel of up to 387 genotypes: 163 winter, 127 spring, and seven semiwinter oilseed rape (OSR) habits, 35 swede, 15 winter fodder, and 40 exotic/unspecified habits. Root traits of 14 d old seedlings were measured in a 'pouch and wick' system (n = ~24 replicates per genotype). The mineral composition of 3-6 rosette-stage leaves, and mature seeds, was determined on compost-grown plants from a designed experiment (n = 5) by inductively coupled plasma-mass spectrometry (ICP-MS).

RESULTS

Seed size explained a large proportion of the variation in root length. Winter OSR and fodder habits had longer primary and lateral roots than spring OSR habits, with generally lower mineral concentrations. A comparison of the ratios of elements in leaf and seed parts revealed differences in translocation processes between crop habits, including those likely to be associated with crop-selection for OSR seeds with lower sulphur-containing glucosinolates. Combining root, leaf and seed traits in a discriminant analysis provided the most accurate characterisation of crop habit, illustrating the interdependence of plant tissues.

CONCLUSIONS

High-throughput morphological and composition phenotyping reveals complex interrelationships between mineral acquisition and accumulation linked to genetic control within and between crop types (habits) in B. napus. Despite its recent genetic ancestry (<10 ky), root morphology, and leaf and seed composition traits could potentially be used in crop improvement, if suitable markers can be identified and if these correspond with suitable agronomy and quality traits.

Elucidation of the genetic basis of variation for stem strength characteristics in bread wheat by Associative Transcriptomics

Background

The current approach to reducing the tendency for wheat grown under high fertilizer conditions to collapse (lodge) under the weight of its grain is based on reducing stem height via the introduction of Rht genes. However, these reduce the yield of straw (itself an important commodity) and introduce other undesirable characteristics. Identification of alternative height-control loci is therefore of key interest. In addition, the improvement of stem mechanical strength provides a further way through which lodging can be reduced.

Results

To investigate the prospects for genetic alternatives to Rht, we assessed variation for plant height and stem strength properties in a training genetic diversity panel of 100 wheat accessions fixed for Rht. Using mRNAseq data derived from RNA purified from leaves, functional genotypes were developed for the panel comprising 42,066 Single Nucleotide Polymorphism (SNP) markers and 94,060 Gene Expression Markers (GEMs). In the first application in wheat of the recently-developed method of Associative Transcriptomics, we identified associations between trait variation and both SNPs and GEMs. Analysis of marker-trait associations revealed candidates for the causative genes underlying the trait variation, implicating xylan acetylation and the COP9 signalosome as contributing to stem strength and auxin in the control of the observed variation for plant height. Predictive capabilities of key markers for stem strength were validated using a test genetic diversity panel of 30 further wheat accessions.

Conclusions

This work illustrates the power of Associative Transcriptomics for the exploration of complex traits of high agronomic importance in wheat. The careful selection of genotypes included in the analysis, allowed for high resolution mapping of novel trait-controlling loci in this staple crop. The use of Gene Expression markers coupled with the more traditional sequence-based markers, provides the power required to understand the biological context of the marker-trait associations observed. This not only adds to the wealth of knowledge that we strive to accumulate regarding gene function and plant adaptation, but also provides breeders with the information required to make more informed decisions regarding the potential consequences of incorporating the use of particular markers into future breeding programmes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2775-2) contains supplementary material, which is available to authorized users.

Genome distribution of differential homoeologue contributions to leaf gene expression in bread wheat

Using a combination of de novo transcriptome assembly, a newly developed 9495-marker transcriptome SNP genetic linkage map and comparative genomics approaches, we developed an ordered set of nonredundant transcripts for each of the subgenomes of hexaploid wheat: A (47 160 unigenes), B (59 663 unigenes) and D (40 588 unigenes). We used these as reference sequences against which to map Illumina mRNA-Seq reads derived from young leaf tissue. Transcript abundance was quantified for each unigene. Using a three-way reciprocal BLAST approach, 15 527 triplet sets of homoeologues (one from each genome) were identified. Differential expression (P < 0.05) was identified for 5248 unigenes, with 2906 represented at greater abundance than their two homoeologues and 2342 represented at lower abundance than their two homoeologues. Analysis of gene ontology terms revealed no biases between homoeologues. There was no evidence of genomewide dominance effects, rather the more highly transcribed individual genes were distributed throughout all three genomes. Transcriptome display tile plot, a visualization approach based on CMYK colour space, was developed and used to assess the genome for regions of skewed homoeologue transcript abundance. Extensive striation was revealed, indicative of many small regions of genome dominance (transcripts of homoeologues from one genome more abundant than the others) and many larger regions of genome repression (transcripts of homoeologues from one genome less abundant than the others).

Molecular markers for tolerance of European ash (Fraxinus excelsior) to dieback disease identified using Associative Transcriptomics

Tree disease epidemics are a global problem, impacting food security, biodiversity and national economies. The potential for conservation and breeding in trees is hampered by complex genomes and long lifecycles, with most species lacking genomic resources. The European Ash tree Fraxinus excelsior is being devastated by the fungal pathogen Hymenoscyphus fraxineus, which causes ash dieback disease. Taking this system as an example and utilizing Associative Transcriptomics for the first time in a plant pathology study, we discovered gene sequence and gene expression variants across a genetic diversity panel scored for disease symptoms and identified markers strongly associated with canopy damage in infected trees. Using these markers we predicted phenotypes in a test panel of unrelated trees, successfully identifying individuals with a low level of susceptibility to the disease. Co-expression analysis suggested that pre-priming of defence responses may underlie reduced susceptibility to ash dieback.

Identification, evolution, and expression partitioning of miRNAs in allopolyploid Brassica napus

The recently published genome of Brassica napus offers for the first time the opportunity to gain insights into the genomic organization and the evolution of miRNAs in oilseed rape. In this study, 12 small RNA libraries from two B. napus cultivars (Tapidor and Ningyou7) and their four double-haploid lines were sequenced, employing the newly sequenced B. napus genome, together with genomes of its progenitors Brassica rapa and Brassica oleracea. A total of 645 miRNAs including 280 conserved and 365 novel miRNAs were identified. Comparative analysis revealed a high level of genomic conservation of MIRNAs (75.9%) between the subgenomes of B. napus and its two progenitors' genomes, and MIRNA lost/gain events (133) occurred in B. napus after its speciation. Furthermore, significant partitioning of miRNA expressions between the two subgenomes in B. napus was detected. The data of degradome sequencing, miRNA-mediated cleavage, and expression analyses support specific interactions between miRNAs and their targets in the modulation of diverse physiological processes in roots and leaves, as well as in biosynthesis of, for example, glucosinolates and lipids in oilseed rape. These data provide a first genome-wide view on the origin, evolution, and genomic organization of B. napus MIRNAs.

Collinearity analysis of Brassica A and C genomes based on an updated inferred unigene order

This data article includes SNP scoring across lines of the Brassica napus TNDH population based on Illumina sequencing of mRNA, expanded to 75 lines. The 21, 323 mapped markers defined 887 recombination bins, representing an updated genetic linkage map for the species. Based on this new map, 5 genome sequence scaffolds were split and the order and orientation of scaffolds updated to establish a new pseudomolecule specification. The order of unigenes and SNP array probes within these pseudomolecules was determined. Unigenes were assessed for sequence similarity to the A and C genomes. The 57, 246 that mapped to both enabled the collinearity of the A and C genomes to be illustrated graphically. Although the great majority was in collinear positions, some were not. Analyses of 60 such instances are presented, suggesting that the breakdown in collinearity was largely due to either the absence of the homoeologue on one genome (resulting in sequence match to a paralogue) or multiple similar sequences being present. The mRNAseq datasets for the TNDH lines are available from the SRA repository (ERA283648); the remaining datasets are supplied with this article.

Construction of Brassica A and C genome-based ordered pan-transcriptomes for use in rapeseed genomic research

This data article reports the establishment of the first pan-transcriptome resources for the Brassica A and C genomes. These were developed using existing coding DNA sequence (CDS) gene models from the now-published Brassica oleracea TO1000 and Brassica napus Darmor-bzh genome sequence assemblies representing the chromosomes of these species, along with preliminary CDS models from an updated Brassica rapa Chiifu genome sequence assembly. The B. rapa genome sequence scaffolds required splitting and re-ordering to match the expected genome organisation based on a high density SNP linkage map, but the B. oleracea assembly was used unchanged. The resulting B. rapa (A genome) pseudomolecules contained 47,656 ordered CDS models and the B. oleracea (C genome) pseudomolecules contained 54,766 ordered CDS models. Interpolation of B. napus CDS models not already represented by orthologues resulted in 52,790 and 63,308 ordered CDS models in the A and C pan-transcriptomes, an increase of 13,676 overall. Comparison of the organisation of this resource with publicly available genome sequences for B. napus showed excellent consistency for the B. napus Darmor-bzh resource, but more breakdown of collinearity for the B. napus ZS11 resource. CDS datasets comprising the pan-transcriptomes are available with this article (B. rapa) or from public repositories (B. oleracea and B. napus).

Identification of a gene controlling variation in the salt tolerance of rapeseed (Brassica napus L.)

By genome-wide association study, QTLs for salt tolerance in rapeseed were detected, and a TSN1 ortholog was identified as a candidate gene responsible for genetic variation in cultivars. Dissecting the genomic regions governing abiotic stress tolerance is necessary for marker-assisted breeding to produce elite breeding lines. In this study, a world-wide collection of rapeseed was evaluated for salt tolerance. These rapeseed accessions showed a large variation for salt tolerance index ranging from 0.311 to 0.999. Although no significant correlation between salt tolerance and Na(+) content was observed, there was a significant negative correlation between shoot biomass production under a control condition and salt tolerance. These rapeseed accessions were genotyped by DArTseq for a total of 51,109 genetic markers, which were aligned with 'pseudomolecules' representative of the genome of rapeseed to locate their hypothetical order for association mapping. A total of 62 QTLs for salt tolerance, shoot biomass, and ion-homeostasis-related traits were identified by association mapping using both the P and Q+K models. Candidate genes located within the QTL regions were also shortlisted. Sequence analysis showed many polymorphisms for BnaaTSN1. Three of them in the coding region resulting in a premature stop codon or frameshift were found in most of the sensitive lines. Loss-of-function mutations showed a significant association with salt tolerance in B. napus.

Seedling development traits in Brassica napus examined by gene expression analysis and association mapping

BACKGROUND

An optimal seedling development of Brassica napus plants leads to a higher yield stability even under suboptimal growing conditions and has therefore a high importance for plant breeders. The objectives of our study were to (i) examine the expression levels of candidate genes in seedling leaves of B. napus and correlate these with seedling development as well as (ii) detect genome regions associated with gene expression levels and seedling development traits in B. napus by genome-wide association mapping.

RESULTS

The expression levels of the 15 candidate genes examined in the 509 B. napus inbreds showed an averaged standard deviation of 5.6 across all inbreds and ranged from 3.2 to 8.8. The gene expression differences between the 509 B. napus inbreds were more than adequate for the correlation with phenotypic variation of seedling development. The average of the absolute value correlations of the correlation coefficients of 0.11 were observed with a range from 0.00 to 0.39. The candidate genes GER1, AILP1, PECT, and FBP were strongly correlated with the seedling development traits. In a genome-wide association study, we detected a total of 63 associations between single nucleotide polymorphisms (SNPs) and the seedling development traits and 31 SNP-gene associations for the candidate genes with a P-value < 0.0001. For the projected leaf area traits we identified five different association hot spots on the chromosomes A2, A7, C3, C6, and C7.

CONCLUSION

A total of 99.4% of the adjacent SNPs on the A genome and 93.0% of the adjacent SNPs on the C genome had a distance smaller than the average range of linkage disequilibrium. Therefore, this genome-wide association study is expected to result on average in 14.7% of the possible power. Compared to previous studies in B. napus, the SNP marker density of our study is expected to provide a higher power to detect SNP-trait/-gene associations in the B. napus diversity set. The large number of associations detected for the examined 14 seedling development traits indicated that these are genetically complex inherited. The results of our analyses suggested that the studied genes ribulose 1,5-bisphosphate carboxylase/oxygenase small subunit (RBC) on the chromosomes A4 and C4 and fructose-1,6-bisphosphatase precursor (FBP) on the chromosomes A9 and C8 are cis-regulated.

Effect of Brassica napus cultivar on cellulosic ethanol yield

BACKGROUND

Intraspecific variations in biomass composition are likely to influence their suitability for biorefining. This may be particularly important in species such as Brassica napus, which contain many different crop types bred for different purposes. Here, straw derived from 17 B. napus cultivars, of varying crop types, were steam exploded, saccharified and fermented to establish differences in biomass composition relevant to cellulosic ethanol production.

RESULTS

Despite being grown and processed in the same manner, straw from the various cultivars produced different saccharification and fermentation yields after processing. Fermentation inhibitor abundances released by steam explosion also varied between genotypes. Cultivars with glucan-rich straw did not necessarily produce higher saccharification or ethanol yields after processing. Instead, the compositions of non-cellulosic components were more reliable indicators of substrate quality. The abundance of pectins and arabinogalactans had the greatest influence on saccharification efficiency between straw genotypes.

CONCLUSIONS

In dicotyledonous species, such as B. napus, variations in the abundance of pectins between crop cultivars are likely to influence processing efficiency for bioethanol production. Knowledge of these genotypic variants provides targets for plant breeding and could aid in the development of improved cellulase cocktails.

Associative transcriptomics study dissects the genetic architecture of seed glucosinolate content in Brassica napus

Breeding new varieties with low seed glucosinolate (GS) concentrations has long been a prime target in Brassica napus. In this study, a novel association mapping methodology termed 'associative transcriptomics' (AT) was applied to a panel of 101 B. napus lines to define genetic regions and also candidate genes controlling total seed GS contents. Over 100,000 informative single-nucleotide polymorphisms (SNPs) and gene expression markers (GEMs) were developed for AT analysis, which led to the identification of 10 SNP and 7 GEM association peaks. Within these peaks, 26 genes were inferred to be involved in GS biosynthesis. A weighted gene co-expression network analysis provided additional 40 candidate genes. The transcript abundance in leaves of two candidate genes, BnaA.GTR2a located on chromosome A2 and BnaC.HAG3b on C9, was correlated with seed GS content, explaining 18.8 and 16.8% of phenotypic variation, respectively. Resequencing of genomic regions revealed six new SNPs in BnaA.GTR2a and four insertions or deletions in BnaC.HAG3b. These deletion polymorphisms were then successfully converted into polymerase chain reaction-based diagnostic markers that can, due to high linkage disequilibrium observed in these regions of the genome, be used for marker-assisted breeding for low seed GS lines.

Dissection of the control of anion homeostasis by associative transcriptomics in Brassica napus

To assess the variation in nutrient homeostasis in oilseed rape and to identify the genes responsible for this variation, we determined foliar anion levels in a diversity panel of Brassica napus accessions, 84 of which had been genotyped previously using messenger RNA sequencing. We applied associative transcriptomics to identify sequence polymorphisms linked to variation in nitrate, phosphate, or sulfate in these accessions. The analysis identified several hundred significant associations for each anion. Using functional annotation of Arabidopsis (Arabidopsis thaliana) homologs and available microarray data, we identified 60 candidate genes for controlling variation in the anion contents. To verify that these genes function in the control of nutrient homeostasis, we obtained Arabidopsis transfer DNA insertion lines for these candidates and tested them for the accumulation of nitrate, phosphate, and sulfate. Fourteen lines differed significantly in levels of the corresponding anions. Several of these genes have been shown previously to affect the accumulation of the corresponding anions in Arabidopsis mutants. These results thus confirm the power of associative transcriptomics in dissection of the genetic control of complex traits and present a set of candidate genes for use in the improvement of efficiency of B. napus mineral nutrition.

Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome

Oilseed rape (Brassica napus L.) was formed ~7500 years ago by hybridization between B. rapa and B. oleracea, followed by chromosome doubling, a process known as allopolyploidy. Together with more ancient polyploidizations, this conferred an aggregate 72× genome multiplication since the origin of angiosperms and high gene content. We examined the B. napus genome and the consequences of its recent duplication. The constituent An and Cn subgenomes are engaged in subtle structural, functional, and epigenetic cross-talk, with abundant homeologous exchanges. Incipient gene loss and expression divergence have begun. Selection in B. napus oilseed types has accelerated the loss of glucosinolate genes, while preserving expansion of oil biosynthesis genes. These processes provide insights into allopolyploid evolution and its relationship with crop domestication and improvement.

Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome

Oilseed rape (Brassica napus L.) was formed ~7500 years ago by hybridization between B. rapa and B. oleracea, followed by chromosome doubling, a process known as allopolyploidy. Together with more ancient polyploidizations, this conferred an aggregate 72× genome multiplication since the origin of angiosperms and high gene content. We examined the B. napus genome and the consequences of its recent duplication. The constituent An and Cn subgenomes are engaged in subtle structural, functional, and epigenetic cross-talk, with abundant homeologous exchanges. Incipient gene loss and expression divergence have begun. Selection in B. napus oilseed types has accelerated the loss of glucosinolate genes, while preserving expansion of oil biosynthesis genes. These processes provide insights into allopolyploid evolution and its relationship with crop domestication and improvement.

Variation in the chemical composition of wheat straw: the role of tissue ratio and composition

BACKGROUND

Wheat straw is an attractive substrate for second generation ethanol production because it will complement and augment wheat production rather than competing with food production. However, like other sources of lignocellulosic biomass, even from a single species, it is heterogeneous in nature due to the different tissues and cell types, and this has implications for saccharification efficiency. The aim of this study has been to use Fourier transform infrared (FTIR) spectroscopy and Partial least squares (PLS) modelling to rapidly screen wheat cultivars for the levels of component tissues, the carbohydrate composition and lignin content, and the levels of simple cross-linking phenolics such as ferulic and diferulic acids.

RESULTS

FTIR spectroscopy and PLS modelling was used to analyze the tissue and chemical composition of wheat straw biomass. Predictive models were developed to evaluate the variability in the concentrations of the cell wall sugars, cell wall phenolics and acid-insoluble lignin. Models for the main sugars, phenolics and lignin were validated and then used to evaluate the variation in total biomass composition across 90 cultivars of wheat grown over two seasons.

CONCLUSIONS

Whilst carbohydrate and lignin components varied across the varieties, this mainly reflected differences in the ratios of the component tissues rather than differences in the composition of those tissues. Further analysis indicated that on a mol% basis, relative levels of sugars within the tissues varied to only a small degree. There were no clear associations between simple phenolics and tissues. The results provide a basis for improving biomass quality for biofuels production through selection of cultivars with appropriate tissue ratios.

Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (Brassica napus L.)

Association mapping can quickly and efficiently dissect complex agronomic traits. Rapeseed is one of the most economically important polyploid oil crops, although its genome sequence is not yet published. In this study, a recently developed 60K Brassica Infinium(®) SNP array was used to analyse an association panel with 472 accessions. The single-nucleotide polymorphisms (SNPs) of the array were in silico mapped using 'pseudomolecules' representative of the genome of rapeseed to establish their hypothetical order and to perform association mapping of seed weight and seed quality. As a result, two significant associations on A8 and C3 of Brassica napus were detected for erucic acid content, and the peak SNPs were found to be only 233 and 128 kb away from the key genes BnaA.FAE1 and BnaC.FAE1. BnaA.FAE1 was also identified to be significantly associated with the oil content. Orthologues of Arabidopsis thaliana HAG1 were identified close to four clusters of SNPs associated with glucosinolate content on A9, C2, C7 and C9. For seed weight, we detected two association signals on A7 and A9, which were consistent with previous studies of quantitative trait loci mapping. The results indicate that our association mapping approach is suitable for fine mapping of the complex traits in rapeseed.

Genome-wide investigation of genetic changes during modern breeding of Brassica napus

Considerable genome variation had been incorporated within rapeseed breeding programs over past decades. In past decades, there have been substantial changes in phenotypic properties of rapeseed as a result of extensive breeding effort. Uncovering the underlying patterns of allelic variation in the context of genome organisation would provide knowledge to guide future genetic improvement. We assessed genome-wide genetic changes, including population structure, genetic relatedness, the extent of linkage disequilibrium, nucleotide diversity and genetic differentiation based on F ST outlier detection, for a panel of 472 Brassica napus inbred accessions using a 60 k Brassica Infinium® SNP array. We found genetic diversity varied in different sub-groups. Moreover, the genetic diversity increased from 1950 to 1980 and then remained at a similar level in China and Europe. We also found ~6-10 % genomic regions revealed high F ST values. Some QTLs previously associated with important agronomic traits overlapped with these regions. Overall, the B. napus C genome was found to have more high F ST signals than the A genome, and we concluded that the C genome may contribute more valuable alleles to generate elite traits. The results of this study indicate that considerable genome variation had been incorporated within rapeseed breeding programs over past decades. These results also contribute to understanding the impact of rapeseed improvement on available genome variation and the potential for dissecting complex agronomic traits.

Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea

Background

Brassica oleracea is a valuable vegetable species that has contributed to human health and nutrition for hundreds of years and comprises multiple distinct cultivar groups with diverse morphological and phytochemical attributes. In addition to this phenotypic wealth, B. oleracea offers unique insights into polyploid evolution, as it results from multiple ancestral polyploidy events and a final Brassiceae-specific triplication event. Further, B. oleracea represents one of the diploid genomes that formed the economically important allopolyploid oilseed, Brassica napus. A deeper understanding of B. oleracea genome architecture provides a foundation for crop improvement strategies throughout the Brassica genus.

Results

We generate an assembly representing 75% of the predicted B. oleracea genome using a hybrid Illumina/Roche 454 approach. Two dense genetic maps are generated to anchor almost 92% of the assembled scaffolds to nine pseudo-chromosomes. Over 50,000 genes are annotated and 40% of the genome predicted to be repetitive, thus contributing to the increased genome size of B. oleracea compared to its close relative B. rapa. A snapshot of both the leaf transcriptome and methylome allows comparisons to be made across the triplicated sub-genomes, which resulted from the most recent Brassiceae-specific polyploidy event.

Conclusions

Differential expression of the triplicated syntelogs and cytosine methylation levels across the sub-genomes suggest residual marks of the genome dominance that led to the current genome architecture. Although cytosine methylation does not correlate with individual gene dominance, the independent methylation patterns of triplicated copies suggest epigenetic mechanisms play a role in the functional diversification of duplicate genes.

Steam explosion of oilseed rape straw: establishing key determinants of saccharification efficiency

Oilseed rape straw was steam exploded into hot water at a range of severities. The residues were fractionated into solid and liquid phases and chemically characterised. The effect of steam explosion on enzymatic hydrolysis of the water-insoluble fractions was investigated by studying initial cellulase binding and hydrolysis yields for different cellulase doses. Time-course data was modelled to establish rate-dependent differences in saccharification as a function of pretreatment severity and associated chemical composition. The study concluded: (1) the initial hydrolysis rate was limited by the amount of (pectic) uronic acid remaining in the substrate; (2) the proportion of rapidly hydrolysable carbohydrate was most closely and positively related to lignin abundance and (3) the final sugar yield most closely related to xylan removal from the substrate. Comparisons between milled and un-milled steam exploded straw highlighted the influence that physical structure has on hydrolysis rates and yields, particularly at low severities.

The control of seed oil polyunsaturate content in the polyploid crop species Brassica napus

Many important plant species have polyploidy in their recent ancestry, complicating inferences about the genetic basis of trait variation. Although the principal locus controlling the proportion of polyunsaturated fatty acids (PUFAs) in seeds of Arabidopsis thaliana is known (fatty acid desaturase 2; FAD2), commercial cultivars of a related crop, oilseed rape (Brassica napus), with very low PUFA content have yet to be developed. We showed that a cultivar of oilseed rape with lower than usual PUFA content has non-functional alleles at three of the four orthologous FAD2 loci. To explore the genetic basis further, we developed an ethyl methanesulphonate mutagenised population, JBnaCAB_E, and used it to identify lines that also carried mutations in the remaining functional copy. This confirmed the hypothesised basis of variation, resulting in an allelic series of mutant lines showing a spectrum of PUFA contents of seed oil. Several lines had PUFA content of ~6 % and oleic acid content of ~84 %, achieving a long-standing industry objective: very high oleic, very low PUFA rapeseed without the use of genetic modification technology. The population contains a high rate of mutations and represents an important resource for research in B. napus.

Sequencing-based variant detection in the polyploid crop oilseed rape

BACKGROUND

The detection and exploitation of genetic variation underpins crop improvement. However, the polyploid nature of the genomes of many of our most important crops represents a barrier, particularly for the analysis of variation within genes. To overcome this, we aimed to develop methodologies based on amplicon sequencing that involve the incorporation of barcoded amplification tags (BATs) into PCR products.

RESULTS

A protocol was developed to tag PCR products with 5' 6-base oligonucleotide barcode extensions before pooling for sequencing library production using standard Illumina adapters. A computational method was developed for the de-convolution of products and the robust detection and scoring of sequence variants. Using this methodology, amplicons targeted to gene sequences were screened across a B. napus mapping population and the resulting allele scoring strings for 24 markers linkage mapped to the expected regions of the genome. Furthermore, using one-dimensional 8-fold pooling, 4608 lines of a B. napus mutation population were screened for induced mutations in a locus-specific amplicon (an orthologue of GL2.b) and mixed product of three co-amplified loci (orthologues of FAD2), identifying 10 and 41 mutants respectively.

CONCLUSIONS

The utilisation of barcode tags to de-convolute pooled PCR products in multiplexed, variation screening via Illumina sequencing provides a cost effective method for SNP genotyping and mutation detection and, potentially, markers for causative changes, even in polyploid species. Combining this approach with existing Illumina multiplexing workflows allows the analysis of thousands of lines cheaply and efficiently in a single sequencing run with minimal library production costs.