The story of a decade: Genomics, functional genomics, and molecular breeding in Brassica napus

Rapeseed (Brassica napus L.) is one of the major global sources of edible vegetable oil and is also used as a feed and pioneer crop and for sightseeing and industrial purposes. Improvements in genome sequencing and molecular marker technology have fueled a boom in functional genomic studies of major agronomic characters such as yield, quality, flowering time, and stress resistance. Moreover, introgression and pyramiding of key functional genes have greatly accelerated the genetic improvement of important traits. Here we summarize recent progress in rapeseed genomics and genetics, and we discuss effective molecular breeding strategies by exploring these findings in rapeseed. These insights will extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture throughout the world.

Cytokinin regulates the activity of the inflorescence meristem and components of seed yield in oilseed rape

The number of flowers and seed-bearing structures formed by the inflorescence meristem and the formation of ovules in the female reproductive part of the flowers are important yield-related traits of crop plants. It has been shown that cytokinin is a pivotal factor regulating these traits. Here, we explore the impact of mutation of CYTOKININ OXIDASE/DEHYDROGENASE (CKX) genes encoding cytokinin-degrading enzymes on these yield-related traits in oilseed rape (Brassica napus L.). We describe the identification of four BnCKX3 and two BnCKX5 genes as regulators of reproductive development in the allotetraploid B. napus. RNA-seq analysis and in situ hybridization showed expression of these genes in reproductive organs. Loss-of-function mutants for each of these CKX gene copies were identified by targeting induced local lesions in genomes (TILLING) and combined by crossing. Sextuple ckx3 ckx5 mutants showed an increased cytokinin concentration and larger and more active inflorescence meristems. They also produced up to 72% more flowers with gynoecia containing 32% more ovules and up to 54% more pods on the main stem. The weight of seeds harvested from the main stem of plants grown in the greenhouse or in the field was increased by 20-32%. Our results show that cytokinin regulates inflorescence meristem and placenta activity in oilseed rape. The work demonstrates the potential to achieve yield enhancement in a dicot crop plant by modulating the cytokinin status through mutagenesis of specific CKX genes.

Quantitative Trait Locus Mapping and Identification of Candidate Genes Controlling Flowering Time in Brassica napus L

Flowering time plays a vital role in determining the life-cycle period, yield, and seed quality of rapeseed (Brassica napus L.) in certain environments. Quantitative trait locus (QTL) mapping to identify the genetic architecture of genes controlling flowering time helps accelerate the early maturity breeding process. In this study, simple sequence repeats (SSR) and specific-locus amplified fragment sequencing (SLAF-seq) technologies were adopted to map the QTLs for flowering time in four environments. As a result, three target intervals, FTA09, FTA10, and FTC05 were identified. Among this, FTA09 was considered as a novel interval, FTA10 and FTC05 as stable regions. Based on the parental re-sequencing data, 7,022 single nucleotide polymorphisms (SNPs) and 2,195 insertion-deletions (InDels) between the two parents were identified in these three target regions. A total of 186 genes possessed genetic variations in these intervals, 14 of which were related to flowering time involved in photoperiod, circadian clock, vernalization, and gibberellin pathways. Six InDel markers linked to flowering time were developed in the three target intervals, indicating that the results were credible in this study. These results laid a good foundation for further genetic studies on flowering-time regulation in B. napus L.

The vernalisation regulator FLOWERING LOCUS C is differentially expressed in biennial and annual Brassica napus

Plants in temperate areas evolved vernalisation requirement to avoid pre-winter flowering. In Brassicaceae, a period of extended cold reduces the expression of the flowering inhibitor FLOWERING LOCUS C (FLC) and paves the way for the expression of downstream flowering regulators. As with all polyploid species of the Brassicaceae, the model allotetraploid Brassica napus (rapeseed, canola) is highly duplicated and carries 9 annotated copies of Bna.FLC. To investigate whether these multiple homeologs and paralogs have retained their original function in vernalisation or undergone subfunctionalisation, we compared the expression patterns of all 9 copies between vernalisation-dependent (biennial, winter type) and vernalisation-independent (annual, spring type) accessions, using RT-qPCR with copy-specific primers and RNAseq data from a diversity set. Our results show that only 3 copies – Bna.FLC.A03b, Bna.FLC.A10 and to some extent Bna.FLC.C02 – are differentially expressed between the two growth types, showing that expression of the other 6 copies does not correlate with growth type. One of those 6 copies, Bna.FLC.C03b, was not expressed at all, indicating a pseudogene, while three further copies, Bna.FLC.C03a and Bna.FLC.C09ab, did not respond to cold treatment. Sequence variation at the COOLAIR binding site of Bna.FLC.A10 was found to explain most of the variation in gene expression. However, we also found that Bna.FLC.A10 expression is not fully predictive of growth type.

Genetic dissection of the mechanism of flowering time based on an environmentally stable and specific QTL in Brassica napus

Flowering time is an important agronomic trait that is highly influenced by the environment. To elucidate the genetic mechanism of flowering time in rapeseed (Brassica napus L.), a genome-wide QTL analysis was performed in a doubled haploid population grown in winter, semi-winter and spring ecological conditions. Fifty-five consensus QTLs were identified after combining phenotype and genomic data, including 12 environment-stable QTLs and 43 environment-specific QTLs. Importantly, six major QTLs for flowering time were identified, of which two were considered environment-specific QTLs in spring ecological condition and four were considered environment-stable QTLs in winter and semi-winter ecological conditions. Through QTL comparison, 18 QTLs were colocalized with QTLs from six other published studies. Combining the candidate genes with their functional annotation, in 49 of 55 consensus QTLs, 151 candidate genes in B. napus corresponding to 95 homologous genes in Arabidopsis thaliana related to flowering were identified, including BnaC03g32910D (CO), BnaA02g12130D (FT) and BnaA03g13630D (FLC). Most of the candidate genes were involved in different flowering regulatory pathways. Based on re-sequencing and differences in sequence annotation between the two parents, we found that regions containing some candidate genes have numerous non-frameshift InDels and many non- synonymous mutations, which might directly lead to gene functional variation. Flowering time was negativly correlated with seed yield and thousand seed weight based on a QTL comparison of flowering time and seed yield traits, which has implications in breeding new early-maturing varieties of B. napus. Moreover, a putative flowering regulatory network was constructed, including the photoperiod, circadian clock, vernalization, autonomous and gibberellin pathways. Multiple copies of genes led to functional difference among the different copies of homologous genes, which also increased the complexity of the flowering regulatory networks. Taken together, the present results not only provide new insights into the genetic regulatory network underlying the control of flowering time but also improve our understanding of flowering time regulatory pathways in rapeseed.

Current status and future possibilities of molecular genetics techniques in Brassica napus

As PCR methods have improved over the last 15 years, there has been an upsurge in the number of new DNA marker tools, which has allowed the generation of high-density molecular maps for all the key Brassica crop types. Biotechnology and molecular plant breeding have emerged as a significant tool for molecular understanding that led to a significant crop improvement in the Brassica napus species. Brassica napus possess a very complicated polyploidy-based genomics. The quantitative trait locus (QTL) is not sufficient to develop effective markers for trait introgression. In the coming years, the molecular marker techniques will be more effective to determine the whole genome impairing desired traits. Available genetic markers using the single-nucleotide sequence (SNP) technique and high-throughput sequencing are effective in determining the maps and genome polymorphisms amongst candidate genes and allele interactions. High-throughput sequencing and gene mapping techniques are involved in discovering new alleles and gene pairs, serving as a bridge between the gene map and genome evaluation. The decreasing cost for DNA sequencing will help in discovering full genome sequences with less resources and time. This review describes (1) the current use of integrated approaches, such as molecular marker technologies, to determine genome arrangements and interspecific outcomes combined with cost-effective genomes to increase the efficiency in prognostic breeding efforts. (2) It also focused on functional genomics, proteomics and field-based breeding practices to achieve insight into the genetics underlying both simple and complex traits in canola.

Diversity and Genome Analysis of Australian and Global Oilseed Brassica napus L. Germplasm Using Transcriptomics and Whole Genome Re-sequencing

Intensive breeding of Brassica napus has resulted in relatively low diversity, such that B. napus would benefit from germplasm improvement schemes that sustain diversity. As such, samples representative of global germplasm pools need to be assessed for existing population structure, diversity and linkage disequilibrium (LD). Complexity reduction genotyping-by-sequencing (GBS) methods, including GBS-transcriptomics (GBS-t), enable cost-effective screening of a large number of samples, while whole genome re-sequencing (WGR) delivers the ability to generate large numbers of unbiased genomic single nucleotide polymorphisms (SNPs), and identify structural variants (SVs). Furthermore, the development of genomic tools based on whole genomes representative of global oilseed diversity and orientated by the reference genome has substantial industry relevance and will be highly beneficial for canola breeding. As recent studies have focused on European and Chinese varieties, a global diversity panel as well as a substantial number of Australian spring types were included in this study. Focusing on industry relevance, 633 varieties were initially genotyped using GBS-t to examine population structure using 61,037 SNPs. Subsequently, 149 samples representative of global diversity were selected for WGR and both data sets used for a side-by-side evaluation of diversity and LD. The WGR data was further used to develop genomic resources consisting of a list of 4,029,750 high-confidence SNPs annotated using SnpEff, and SVs in the form of 10,976 deletions and 2,556 insertions. These resources form the basis of a reliable and repeatable system allowing greater integration between canola genomics studies, with a strong focus on breeding germplasm and industry applicability.

Oilseed rape: learning about ancient and recent polyploid evolution from a recent crop species

Oilseed rape (Brassica napus) is one of our youngest crop species, arising several times under cultivation in the last few thousand years and completely unknown in the wild. Oilseed rape originated from hybridisation events between progenitor diploid species B. rapa and B. oleracea, both important vegetable species. The diploid progenitors are also ancient polyploids, with remnants of two previous polyploidisation events evident in the triplicated genome structure. This history of polyploid evolution and human agricultural selection makes B. napus an excellent model with which to investigate processes of genomic evolution and selection in polyploid crops. The ease of de novo interspecific hybridisation, responsiveness to tissue culture, and the close relationship of oilseed rape to the model plant Arabidopsis thaliana, coupled with the recent availability of reference genome sequences and suites of molecular cytogenetic and high-throughput genotyping tools, allow detailed dissection of genetic, genomic and phenotypic interactions in this crop. In this review we discuss the past and present uses of B. napus as a model for polyploid speciation and evolution in crop species, along with current and developing analysis tools and resources. We further outline unanswered questions that may now be tractable to investigation.

A high-density SNP genotyping array for Brassica napus and its ancestral diploid species based on optimised selection of single-locus markers in the allotetraploid genome

The Brassica napus Illumina array provides genome-wide markers linked to the available genome sequence, a significant tool for genetic analyses of the allotetraploid B. napus and its progenitor diploid genomes. A high-density single nucleotide polymorphism (SNP) Illumina Infinium array, containing 52,157 markers, was developed for the allotetraploid Brassica napus. A stringent selection process employing the short probe sequence for each SNP assay was used to limit the majority of the selected markers to those represented a minimum number of times across the highly replicated genome. As a result approximately 60 % of the SNP assays display genome-specificity, resolving as three clearly separated clusters (AA, AB, and BB) when tested with a diverse range of B. napus material. This genome specificity was supported by the analysis of the diploid ancestors of B. napus, whereby 26,504 and 29,720 markers were scorable in B. oleracea and B. rapa, respectively. Forty-four percent of the assayed loci on the array were genetically mapped in a single doubled-haploid B. napus population allowing alignment of their physical and genetic coordinates. Although strong conservation of the two positions was shown, at least 3 % of the loci were genetically mapped to a homoeologous position compared to their presumed physical position in the respective genome, underlying the importance of genetic corroboration of locus identity. In addition, the alignments identified multiple rearrangements between the diploid and tetraploid Brassica genomes. Although mostly attributed to genome assembly errors, some are likely evidence of rearrangements that occurred since the hybridisation of the progenitor genomes in the B. napus nucleus. Based on estimates for linkage disequilibrium decay, the array is a valuable tool for genetic fine mapping and genome-wide association studies in B. napus and its progenitor genomes.

Transcriptome Analysis of Stem and Globally Comparison with Other Tissues in Brassica napus

Brassica napus is one of the most important oilseed crops in the world. However, there is currently no enough stem transcriptome information and comparative transcriptome analysis of different tissues, which impedes further functional genomics research on B. napus. In this study, the stem transcriptome of B. napus was characterized by RNA-seq technology. Approximately 13.4 Gb high-quality clean reads with an average length of 100 bp were generated and used for comparative transcriptome analysis with the existing transcriptome sequencing data of roots, leaves, flower buds, and immature embryos of B. napus. All the transcripts were annotated against GO and KEGG databases. The common genes in five tissues, differentially expressed genes (DEGs) of the common genes between stems and other tissues, and tissue-specific genes were detected, and the main biochemical activities and pathways implying the common genes, DEGs and tissue-specific genes were investigated. Accordingly, the common transcription factors (TFs) in the five tissues and tissue-specific TFs were identified, and a TFs-based regulation network between TFs and the target genes involved in 'Phenylpropanoid biosynthesis' pathway were constructed to show several important TFs and key nodes in the regulation process. Collectively, this study not only provided an available stem transcriptome resource in B. napus, but also revealed valuable comparative transcriptome information of five tissues of B. napus for future investigation on specific processes, functions and pathways.

Hairy Canola (Brasssica napus) re-visited: Down-regulating TTG1 in an AtGL3-enhanced hairy leaf background improves growth, leaf trichome coverage, and metabolite gene expression diversity

BACKGROUND

Through evolution, some plants have developed natural resistance to insects by having hairs (trichomes) on leaves and other tissues. The hairy trait has been neglected in Brassica breeding programs, which mainly focus on disease resistance, yield, and overall crop productivity. In Arabidopsis, a network of three classes of proteins consisting of TTG1 (a WD40 repeat protein), GL3 (a bHLH factor) and GL1 (a MYB transcription factor), activates trichome initiation and patterning. Introduction of a trichome regulatory gene AtGL3 from Arabidopsis into semi-glabrous Brassica napus resulted in hairy canola plants which showed tolerance to flea beetles and diamondback moths; however plant growth was negatively affected. In addition, the role of BnTTG1 transcription in the new germplasm was not understood.

RESULTS

Here, we show that two ultra-hairy lines (K-5-8 and K-6-3) with BnTTG1 knock-down in the hairy AtGL3+ B. napus background showed stable enhancement of trichome coverage, density, and length and restored wild type growth similar to growth of the semi-glabrous Westar plant. In contrast, over-expression of BnTTG1 in the hairy AtGL3+ B. napus background gave consistently glabrous plants of very low fertility and poor stability, with only one glabrous plant (O-3-7) surviving to the T3 generation. Q-PCR trichome gene expression data in leaf samples combining several leaf stages for these lines suggested that BnGL2 controlled B. napus trichome length and out-growth and that strong BnTTG1 transcription together with strong GL3 expression inhibited this process. Weak expression of BnTRY in both glabrous and trichome-bearing leaves of B. napus in the latter Q-PCR experiment suggested that TRY may have functions other than as an inhibitor of trichome initiation in the Brassicas. A role for BnTTG1 in the lateral inhibition of trichome formation in neighbouring cells was also proposed for B. napus. RNA sequencing of first leaves identified a much larger array of genes with altered expression patterns in the K-5-8 line compared to the hairy AtGL3(+) B. napus background (relative to the Westar control plant). These genes particularly included transcription factors, protein degradation and modification genes, but also included pathways that coded for anthocyanins, flavonols, terpenes, glucosinolates, alkaloids, shikimates, cell wall biosynthesis, and hormones. A 2nd Q-PCR experiment was conducted on redox, cell wall carbohydrate, lignin, and trichome genes using young first leaves, including T4 O-3-7-5 plants that had partially reverted to yield two linked growth and trichome phenotypes. Most of the trichome genes tested showed to be consistant with leaf trichome phenotypes and with RNA sequencing data in three of the lines. Two redox genes showed highest overall expression in K-5-8 leaves and lowest in O-3-7-5 leaves, while one redox gene and three cell wall genes were consistently higher in the two less robust lines compared with the two robust lines.

CONCLUSION

The data support the strong impact of BnTTG1 knockdown (in the presence of strong AtGL3 expression) at restoring growth, enhancing trichome coverage and length, and enhancing expression and diversity of growth, metabolic, and anti-oxidant genes important for stress tolerance and plant health in B. napus. Our data also suggests that the combination of strong (up-regulated) BnTTG1 expression in concert with strong AtGL3 expression is unstable and lethal to the plant.

The fusion of genomes leads to more options: A comparative investigation on the desulfo-glucosinolate sulfotransferases of Brassica napus and homologous proteins of Arabidopsis thaliana

Sulfotransferases (SOTs) (EC 2.8.2.-) play a crucial role in the glucosinolate (Gl) biosynthesis, by catalyzing the final step of the core glucosinolate formation. In Arabidopsis thaliana the three desulfo (ds)-Gl SOTs AtSOT16, AtSOT17 and AtSOT18 were previously characterized, showing different affinities to ds-Gls. But can the knowledge about these SOTs be generally transferred to other Gl-synthesizing plants? It was investigated how many SOTs are present in the economically relevant crop plant Brassica napus L., and if it is possible to predict their characteristics by sequence analysis. The recently sequenced B. napus is a hybrid of Brassica rapa and Brassica oleracea. By database research, 71 putative functional BnSOT family members were identified and at least eleven of those are putative ds-Gl SOTs. Besides the homologs of AtSOT16 - 18, phylogenetic analyses revealed new subfamilies of ds-Gl SOTs, which are not present in A. thaliana. Three of the B. napus ds-Gl SOT proteins were expressed and purified, and characterized by determining the substrate affinities to different ds-Gls. Two of them, BnSOT16-a and BnSOT16-b, showed a significantly higher affinity to an indolic ds-Gl, similarly to AtSOT16. Additionally, BnSOT17-a was characterized and showed a higher affinity to long chained aliphatic Gls, similarly to AtSOT17. Identification of homologs to AtSOT18 was less reliable, because putative SOT18 sequences are more heterogeneous and confirmation of similar characteristics was not possible.

Single-nucleotide polymorphism identification and genotyping in Camelina sativa

Camelina sativa, a largely relict crop, has recently returned to interest due to its potential as an industrial oilseed. Molecular markers are key tools that will allow C. sativa to benefit from modern breeding approaches. Two complementary methodologies, capture of 3' cDNA tags and genomic reduced-representation libraries, both of which exploited second generation sequencing platforms, were used to develop a low density (768) Illumina GoldenGate single nucleotide polymorphism (SNP) array. The array allowed 533 SNP loci to be genetically mapped in a recombinant inbred population of C. sativa. Alignment of the SNP loci to the C. sativa genome identified the underlying sequenced regions that would delimit potential candidate genes in any mapping project. In addition, the SNP array was used to assess genetic variation among a collection of 175 accessions of C. sativa, identifying two sub-populations, yet low overall gene diversity. The SNP loci will provide useful tools for future crop improvement of C. sativa.

Gene-based SNP discovery and genetic mapping in pea

KEY MESSAGE

Gene-based SNPs were identified and mapped in pea using five recombinant inbred line populations segregating for traits of agronomic importance. Pea (Pisum sativum L.) is one of the world's oldest domesticated crops and has been a model system in plant biology and genetics since the work of Gregor Mendel. Pea is the second most widely grown pulse crop in the world following common bean. The importance of pea as a food crop is growing due to its combination of moderate protein concentration, slowly digestible starch, high dietary fiber concentration, and its richness in micronutrients; however, pea has lagged behind other major crops in harnessing recent advances in molecular biology, genomics and bioinformatics, partly due to its large genome size with a large proportion of repetitive sequence, and to the relatively limited investment in research in this crop globally. The objective of this research was the development of a genome-wide transcriptome-based pea single-nucleotide polymorphism (SNP) marker platform using next-generation sequencing technology. A total of 1,536 polymorphic SNP loci selected from over 20,000 non-redundant SNPs identified using deep transcriptome sequencing of eight diverse Pisum accessions were used for genotyping in five RIL populations using an Illumina GoldenGate assay. The first high-density pea SNP map defining all seven linkage groups was generated by integrating with previously published anchor markers. Syntenic relationships of this map with the model legume Medicago truncatula and lentil (Lens culinaris Medik.) maps were established. The genic SNP map establishes a foundation for future molecular breeding efforts by enabling both the identification and tracking of introgression of genomic regions harbouring QTLs related to agronomic and seed quality traits.

Genome wide SNP identification in chickpea for use in development of a high density genetic map and improvement of chickpea reference genome assembly

BACKGROUND

In the whole genome sequencing, genetic map provides an essential framework for accurate and efficient genome assembly and validation. The main objectives of this study were to develop a high-density genetic map using RAD-Seq (Restriction-site Associated DNA Sequencing) genotyping-by-sequencing (RAD-Seq GBS) and Illumina GoldenGate assays, and to examine the alignment of the current map with the kabuli chickpea genome assembly.

RESULTS

Genic single nucleotide polymorphisms (SNPs) totaling 51,632 SNPs were identified by 454 transcriptome sequencing of Cicer arietinum and Cicer reticulatum genotypes. Subsequently, an Illumina GoldenGate assay for 1,536 SNPs was developed. A total of 1,519 SNPs were successfully assayed across 92 recombinant inbred lines (RILs), of which 761 SNPs were polymorphic between the two parents. In addition, the next generation sequencing (NGS)-based GBS was applied to the same population generating 29,464 high quality SNPs. These SNPs were clustered into 626 recombination bins based on common segregation patterns. Data from the two approaches were used for the construction of a genetic map using a population derived from an intraspecific cross. The map consisted of 1,336 SNPs including 604 RAD recombination bins and 732 SNPs from Illumina GoldenGate assay. The map covered 653 cM of the chickpea genome with an average distance between adjacent markers of 0.5 cM. To date, this is the most extensive genetic map of chickpea using an intraspecific population. The alignment of the map with the CDC Frontier genome assembly revealed an overall conserved marker order; however, a few local inconsistencies within the Cicer arietinum pseudochromosome 1 (Ca1), Ca5 and Ca8 were detected. The map enabled the alignment of 215 unplaced scaffolds from the CDC Frontier draft genome assembly. The alignment also revealed varying degrees of recombination rates and hotspots across the chickpea genome.

CONCLUSIONS

A high-density genetic map using RAD-Seq GBS and Illumina GoldenGate assay was developed and aligned with the existing kabuli chickpea draft genome sequence. The analysis revealed an overall conserved marker order, although some localized inversions between draft genome assembly and the genetic map were detected. The current analysis provides an insight of the recombination rates and hotspots across the chickpea genome.

Capturing sequence variation among flowering-time regulatory gene homologs in the allopolyploid crop species Brassica napus

Flowering, the transition from the vegetative to the generative phase, is a decisive time point in the lifecycle of a plant. Flowering is controlled by a complex network of transcription factors, photoreceptors, enzymes and miRNAs. In recent years, several studies gave rise to the hypothesis that this network is also strongly involved in the regulation of other important lifecycle processes ranging from germination and seed development through to fundamental developmental and yield-related traits. In the allopolyploid crop species Brassica napus, (genome AACC), homoeologous copies of flowering time regulatory genes are implicated in major phenological variation within the species, however the extent and control of intraspecific and intergenomic variation among flowering-time regulators is still unclear. To investigate differences among B. napus morphotypes in relation to flowering-time gene variation, we performed targeted deep sequencing of 29 regulatory flowering-time genes in four genetically and phenologically diverse B. napus accessions. The genotype panel included a winter-type oilseed rape, a winter fodder rape, a spring-type oilseed rape (all B. napus ssp. napus) and a swede (B. napus ssp. napobrassica), which show extreme differences in winter-hardiness, vernalization requirement and flowering behavior. A broad range of genetic variation was detected in the targeted genes for the different morphotypes, including non-synonymous SNPs, copy number variation and presence-absence variation. The results suggest that this broad variation in vernalization, clock and signaling genes could be a key driver of morphological differentiation for flowering-related traits in this recent allopolyploid crop species.

Genomic DNA enrichment using sequence capture microarrays: a novel approach to discover sequence nucleotide polymorphisms (SNP) in Brassica napus L

Targeted genomic selection methodologies, or sequence capture, allow for DNA enrichment and large-scale resequencing and characterization of natural genetic variation in species with complex genomes, such as rapeseed canola (Brassica napus L., AACC, 2n=38). The main goal of this project was to combine sequence capture with next generation sequencing (NGS) to discover single nucleotide polymorphisms (SNPs) in specific areas of the B. napus genome historically associated (via quantitative trait loci –QTL– analysis) to traits of agronomical and nutritional importance. A 2.1 million feature sequence capture platform was designed to interrogate DNA sequence variation across 47 specific genomic regions, representing 51.2 Mb of the Brassica A and C genomes, in ten diverse rapeseed genotypes. All ten genotypes were sequenced using the 454 Life Sciences chemistry and to assess the effect of increased sequence depth, two genotypes were also sequenced using Illumina HiSeq chemistry. As a result, 589,367 potentially useful SNPs were identified. Analysis of sequence coverage indicated a four-fold increased representation of target regions, with 57% of the filtered SNPs falling within these regions. Sixty percent of discovered SNPs corresponded to transitions while 40% were transversions. Interestingly, fifty eight percent of the SNPs were found in genic regions while 42% were found in intergenic regions. Further, a high percentage of genic SNPs was found in exons (65% and 64% for the A and C genomes, respectively). Two different genotyping assays were used to validate the discovered SNPs. Validation rates ranged from 61.5% to 84% of tested SNPs, underpinning the effectiveness of this SNP discovery approach. Most importantly, the discovered SNPs were associated with agronomically important regions of the B. napus genome generating a novel data resource for research and breeding this crop species.

Applications and challenges of next-generation sequencing in Brassica species

Next-generation sequencing (NGS) produces numerous (often millions) short DNA sequence reads, typically varying between 25 and 400 bp in length, at a relatively low cost and in a short time. This revolutionary technology is being increasingly applied in whole-genome, transcriptome, epigenome and small RNA sequencing, molecular marker and gene discovery, comparative and evolutionary genomics, and association studies. The Brassica genus comprises some of the most agro-economically important crops, providing abundant vegetables, condiments, fodder, oil and medicinal products. Many Brassica species have undergone the process of polyploidization, which makes their genomes exceptionally complex and can create difficulties in genomics research. NGS injects new vigor into Brassica research, yet also faces specific challenges in the analysis of complex crop genomes and traits. In this article, we review the advantages and limitations of different NGS technologies and their applications and challenges, using Brassica as an advanced model system for agronomically important, polyploid crops. Specifically, we focus on the use of NGS for genome resequencing, transcriptome sequencing, development of single-nucleotide polymorphism markers, and identification of novel microRNAs and their targets. We present trends and advances in NGS technology in relation to Brassica crop improvement, with wide application for sophisticated genomics research into agronomically important polyploid crops.

Ancient orphan crop joins modern era: gene-based SNP discovery and mapping in lentil

Background

The genus Lens comprises a range of closely related species within the galegoid clade of the Papilionoideae family. The clade includes other important crops (e.g. chickpea and pea) as well as a sequenced model legume (Medicago truncatula). Lentil is a global food crop increasing in importance in the Indian sub-continent and elsewhere due to its nutritional value and quick cooking time. Despite this importance there has been a dearth of genetic and genomic resources for the crop and this has limited the application of marker-assisted selection strategies in breeding.

Results

We describe here the development of a deep and diverse transcriptome resource for lentil using next generation sequencing technology. The generation of data in multiple cultivated (L. culinaris) and wild (L. ervoides) genotypes together with the utilization of a bioinformatics workflow enabled the identification of a large collection of SNPs and the subsequent development of a genotyping platform that was used to establish the first comprehensive genetic map of the L. culinaris genome. Extensive collinearity with M. truncatula was evident on the basis of sequence homology between mapped markers and the model genome and large translocations and inversions relative to M. truncatula were identified. An estimate for the time divergence of L. culinaris from L. ervoides and of both from M. truncatula was also calculated.

Conclusions

The availability of the genomic and derived molecular marker resources presented here will help change lentil breeding strategies and lead to increased genetic gain in the future.

A tailed PCR procedure for cost-effective, two-order multiplex sequencing of candidate genes in polyploid plants

Complex polyploid crop genomes can be recalcitrant towards conventional DNA sequencing approaches for allele mining in candidate genes for valuable traits. In the past, this has greatly complicated the transfer of knowledge on promising candidate genes from model plants to even closely related polyploid crops. Next-generation sequencing offers diverse solutions to overcome such difficulties. Here, we present a method for multiplexed 454 sequencing in gene-specific PCR amplicons that can simultaneously address multiple homologues of given target genes. We devised a simple two-step PCR procedure employing a set of barcoded M13/T7 universal fusion primers that enable a cost-effective and efficient amplification of large numbers of target gene amplicons. Sequencing-ready amplicons are generated that can be simultaneously sequenced in pools comprising multiple amplicons from multiple genotypes. High-depth sequencing allows resolution of the resulting sequence reads into contigs representing multiple homologous loci, with only insignificant off-target capture of paralogues or PCR artefacts. In a case study, the procedure was tested in the complex polyploid genome of Brassica napus for a set of nine genes identified in Arabidopsis as candidates for regulation of seed development and oil content. Up to six copies of these genes were expected in B. napus. SNP discovery was performed by pooled multiplex sequencing of 30 amplicons in 20 diverse B. napus accessions with interesting trait variation for oil content, providing a basis for comparative mapping to relevant quantitative trait loci and for subsequent marker-assisted breeding.

Analysis of expressed sequence tags in Brassica napus cotyledons damaged by crucifer flea beetle feeding

The molecular basis of canola ( Brassica napus L.) susceptibility to the crucifer flea beetle (FB, Phyllotreta cruciferae Goeze) was investigated by comparing transcript representation in FB-damaged and undamaged cotyledons. The B. napus cotyledon transcriptome increased and diversified substantially after FB feeding damage. Twenty-two genes encoding proteins with unknown function, six encoding proteins involved in signaling, and a gene encoding a B-box zinc finger transcription factor were moderately or strongly changed in representation with FB feeding damage. Zinc finger and calcium-dependent genes formed the largest portion of transcription factors and signaling factors with changes in representation. Six genes with unknown function, one transcription factor, and one signaling gene specific to the FB-damaged library were co-represented in a FB-damaged leaf library. Out of 188 transcription factor and signaling gene families screened for "early" expression changes, 16 showed changes in expression within 8 h. Four of these early factors were zinc finger genes with representation only in the FB-damaged cotyledon. These genes are now available to test their potential at initiating or specifying cotyledon responses to crucifer FB feeding.

Diversity array technology markers: genetic diversity analyses and linkage map construction in rapeseed (Brassica napus L.)

We developed Diversity Array Technology (DArT) markers for application in genetic studies of Brassica napus and other Brassica species with A or C genomes. Genomic representation from 107 diverse genotypes of B. napus L. var. oleifera (rapeseed, AACC genomes) and B. rapa (AA genome) was used to develop a DArT array comprising 11 520 clones generated using PstI/BanII and PstI/BstN1 complexity reduction methods. In total, 1547 polymorphic DArT markers of high technical quality were identified and used to assess molecular diversity among 89 accessions of B. napus, B. rapa, B. juncea, and B. carinata collected from different parts of the world. Hierarchical cluster and principal component analyses based on genetic distance matrices identified distinct populations clustering mainly according to their origin/pedigrees. DArT markers were also mapped in a new doubled haploid population comprising 131 lines from a cross between spring rapeseed lines ‘Lynx-037DH’ and ‘Monty-028DH’. Linkage groups were assigned on the basis of previously mapped simple sequence repeat (SSRs), intron polymorphism (IP), and gene-based markers. The map consisted of 437 DArT, 135 SSR, 6 IP, and 6 gene-based markers and spanned 2288 cM. Our results demonstrate that DArT markers are suitable for genetic diversity analysis and linkage map construction in rapeseed.