Genetic map construction and QTL mapping of resistance to blackleg (Leptosphaeria maculans) disease in Australian canola (Brassica napus L.) cultivars

Brassica napus Haploid and Double Haploid Production and Its Latest Applications

Rapeseed is one of the most important oil crops in the world. Increasing demand for oil and limited agronomic capabilities of present-day rapeseed result in the need for rapid development of new, superior cultivars. Double haploid (DH) technology is a fast and convenient approach in plant breeding as well as genetic research. Brassica napus is considered a model species for DH production based on microspore embryogenesis; however, the molecular mechanisms underlying microspore reprogramming are still vague. It is known that morphological changes are accompanied by gene and protein expression patterns, alongside carbohydrate and lipid metabolism. Novel, more efficient methods for DH rapeseed production have been reported. This review covers new findings and advances in Brassica napus DH production as well as the latest reports related to agronomically important traits in molecular studies employing the double haploid rapeseed lines.

Construction of a Quantitative Genomic Map, Identification and Expression Analysis of Candidate Genes for Agronomic and Disease-Related Traits in Brassica napus

Rapeseed is the second most important oil crop in the world. Improving seed yield and seed oil content are the two main highlights of the research. Unfortunately, rapeseed development is frequently affected by different diseases. Extensive research has been made through many years to develop elite cultivars with high oil, high yield, and/or disease resistance. Quantitative trait locus (QTL) analysis has been one of the most important strategies in the genetic deciphering of agronomic characteristics. To comprehend the distribution of these QTLs and to uncover the key regions that could simultaneously control multiple traits, 4,555 QTLs that have been identified during the last 25 years were aligned in one unique map, and a quantitative genomic map which involved 128 traits from 79 populations developed in 12 countries was constructed. The present study revealed 517 regions of overlapping QTLs which harbored 2,744 candidate genes and might affect multiple traits, simultaneously. They could be selected to customize super-rapeseed cultivars. The gene ontology and the interaction network of those candidates revealed genes that highly interacted with the other genes and might have a strong influence on them. The expression and structure of these candidate genes were compared in eight rapeseed accessions and revealed genes of similar structures which were expressed differently. The present study enriches our knowledge of rapeseed genome characteristics and diversity, and it also provided indications for rapeseed molecular breeding improvement in the future.

Status and advances in mining for blackleg (Leptosphaeria maculans) quantitative resistance (QR) in oilseed rape (Brassica napus)

KEY MESSAGE

Quantitative resistance (QR) loci discovered through genetic and genomic analyses are abundant in the Brassica napus genome, providing an opportunity for their utilization in enhancing blackleg resistance. Quantitative resistance (QR) has long been utilized to manage blackleg in Brassica napus (canola, oilseed rape), even before major resistance genes (R-genes) were extensively explored in breeding programmes. In contrast to R-gene-mediated qualitative resistance, QR reduces blackleg symptoms rather than completely eliminating the disease. As a polygenic trait, QR is controlled by numerous genes with modest effects, which exerts less pressure on the pathogen to evolve; hence, its effectiveness is more durable compared to R-gene-mediated resistance. Furthermore, combining QR with major R-genes has been shown to enhance resistance against diseases in important crops, including oilseed rape. For these reasons, there has been a renewed interest among breeders in utilizing QR in crop improvement. However, the mechanisms governing QR are largely unknown, limiting its deployment. Advances in genomics are facilitating the dissection of the genetic and molecular underpinnings of QR, resulting in the discovery of several loci and genes that can be potentially deployed to enhance blackleg resistance. Here, we summarize the efforts undertaken to identify blackleg QR loci in oilseed rape using linkage and association analysis. We update the knowledge on the possible mechanisms governing QR and the advances in searching for the underlying genes. Lastly, we lay out strategies to accelerate the genetic improvement of blackleg QR in oilseed rape using improved phenotyping approaches and genomic prediction tools.

A gene-for-gene interaction involving a 'late' effector contributes to quantitative resistance to the stem canker disease in Brassica napus

The control of stem canker disease of Brassica napus (rapeseed), caused by the fungus Leptosphaeria maculans is based largely on plant genetic resistance: single-gene specific resistance (Rlm genes) or quantitative, polygenic, adult-stage resistance. Our working hypothesis was that quantitative resistance partly obeys the gene-for-gene model, with resistance genes 'recognizing' fungal effectors expressed during late systemic colonization. Five LmSTEE (stem-expressed effector) genes were selected and placed under the control of the AvrLm4-7 promoter, an effector gene highly expressed at the cotyledon stage of infection, for miniaturized cotyledon inoculation test screening of a gene pool of 204 rapeseed genotypes. We identified a rapeseed genotype, 'Yudal', expressing hypersensitive response to LmSTEE98. The LmSTEE98-RlmSTEE98 interaction was further validated by inactivation of the LmSTEE98 gene with a CRISPR-Cas9 approach. Isolates with mutated versions of LmSTEE98 induced more severe stem symptoms than the wild-type isolate in 'Yudal'. This single-gene resistance was mapped in a 0.6 cM interval of the 'Darmor_bzh' × 'Yudal' genetic map. One typical gene-for-gene interaction contributes partly to quantitative resistance when L. maculans colonizes the stems of rapeseed. With numerous other effectors specific to stem colonization, our study provides a new route for resistance gene discovery, elucidation of quantitative resistance mechanisms and selection for durable resistance.

Genetic Analysis of a Horizontal Resistance Locus BLMR2 in Brassica napus

Leptosphaeria maculans causes blackleg disease in Brassica napus. The blackleg disease is mainly controlled by resistance genes in B. napus. Previous studies have shown that the blackleg resistant BLMR2 locus that conferred horizontal resistance under field conditions, is located on chromosome A10 of B. napus. The purpose of this study is to fine map this locus and hence identify a candidate gene underlying horizontal resistance. The spectrum of resistance to L. maculans isolates of the resistance locus BLMR2 was analyzed using near isogenic lines, resistant, and susceptible cultivars. The results showed that this locus was horizontally resistant to all isolates tested. Sequence characterized amplified regions (SCAR), simple sequence repeats (SSR), and single nucleotide polymorphism (SNP) markers were developed in the chromosome region of BLMR2 and a fine genetic map was constructed. Two molecular markers narrowed BLMR2 in a 53.37 kb region where six genes were annotated. Among the six annotated genes, BnaA10g11280D/BnaA10g11290D encoding a cytochrome P450 protein were predicted as the candidate of BLMR2. Based on the profiling of pathogen induced transcriptome, three expressed genes in the six annotated genes were identified while only cytochrome P450 showed upregulation. The candidate corresponds to the gene involved in the indole glucosinolate biosynthesis pathway and plant basal defense in Arabidopsis thaliana. The molecular markers identified in this study will allow the quick incorporation of the BLMR2 allele in rapeseed cultivars to enhance blackleg resistance.

An extended KASP-SNP resource for molecular breeding in Chinese cabbage(Brassica rapa L. ssp. pekinensis)

Kompetitive allele-specific PCR (KASP) is a cost-effective single-step SNP genotyping technology, With an objective to enhance the marker repertoire and develop high efficient KASP-SNP markers in Chinese cabbage, we re-sequenced four Chinese cabbage doubled haploid (DH) lines, Y177-47, Y635-10, Y510-1 and Y510-9, and generated a total of more than 38.5 billion clean base pairs. A total of 827,720 SNP loci were identified with an estimated density of 3,217 SNPs/Mb. Further, a total of 387,354 SNPs with at least 30 bp to the next most adjacent SNPs on either side were selected as resource for KASP markers. From this resource, 258 (96.27%) of 268 SNP loci were successfully transformed into KASP-SNP markers using a Roche LightCycler 480-II instrument. Among these markers, 221 (85.66%) were co-dominant markers, 220 (85.27%) were non-synonymous SNPs, and 257 (99.6%) were newly developed markers. In addition, 53 markers were applied for genotyping of 34 Brassica rapa accessions. Cluster analysis separated these 34 accessions into three clusters based on heading types. The millions of SNP loci, a large set of resource for KASP markers, as well as the newly developed KASP markers in this study may facilitate further genetic and molecular breeding studies in Brassica rapa.

Insights on SNP types, detection methods and their utilization in Brassica species: Recent progress and future perspectives

The genus Brassica, family Brassicaceae (Cruciferae), comprises many important species of oil crops, vegetables and medicinal plants including B. rapa, B. oleracea, B. nigra, B. napus, B. juncea, B. carinata. Genomic researches in Brassica species is constrained by polyploidization, mainly due to its complicated genomic structure. However, rapid development of methods for detecting single nucleotide polymorphisms (SNP), such as next generation sequencing and SNP microarray, has accelerated release of reference Brassica species genomes as well as discovery of large numbers and genome-wide SNPs, thus intensifying forward genetics in this genus. In this review, we summarize biological characteristics, classification and various methods for detecting SNPs, focusing on high-throughput techniques. Moreover, we describe the pivotal roles of SNPs in genetic diversity, linkage map construction and QTL mapping, comparative genomics, linkage disequilibrium and genome-wide association studies. These insights are expected to deepen our understanding and guide further advancements in Brassica species research.

Meta-analysis of GWAS in canola blackleg (Leptosphaeria maculans) disease traits demonstrates increased power from imputed whole-genome sequence

Blackleg disease causes yield losses in canola (Brassica napus L.). To identify resistance genes and genomic regions, genome-wide association studies (GWAS) of 585 diverse winter and spring canola accessions were performed using imputed whole-genome sequence (WGS) and transcriptome genotype-by-sequencing (GBSt). Blackleg disease phenotypes were collected across three years in six trials. GWAS were performed in several ways and their respective power was judged by the number of significant single nucleotide polymorphisms (SNP), the false discovery rate (FDR), and the percentage of SNP that validated in additional field trials in two subsequent years. WGS GWAS with 1,234,708 million SNP detected a larger number of significant SNP, achieved a lower FDR and a higher validation rate than GBSt with 64,072 SNP. A meta-analysis combining survival and average internal infection resulted in lower FDR but also lower validation rates. The meta-analysis GWAS identified 79 genomic regions (674 SNP) conferring potential resistance to L. maculans. While several GWAS signals localised in regions of known Rlm genes, fifty-three new potential resistance regions were detected. Seventeen regions had underlying genes with putative functions related to disease defence or stress response in Arabidopsis thaliana. This study provides insight into the genetic architecture and potential molecular mechanisms underlying canola L. maculans resistance.

Genome-Wide Association Mapping Identifies Novel Loci for Quantitative Resistance to Blackleg Disease in Canola

Blackleg disease, caused by the fungal pathogen Leptosphaeria maculans, continues to be a major concern for sustainable production of canola (Brassica napus L.) in many parts of the world. The deployment of effective quantitative resistance (QR) is recognized as a durable strategy in providing natural defense to pathogens. Herein, we uncover loci for resistance to blackleg in a genetically diverse panel of canola accessions by exploiting historic recombination events which occurred during domestication and selective breeding by genome-wide association analysis (GWAS). We found extensive variation in resistance to blackleg at the adult plant stage, including for upper canopy infection. Using the linkage disequilibrium and genetic relationship estimates from 12,414 high quality SNPs, GWAS identified 59 statistically significant and "suggestive" SNPs on 17 chromosomes of B. napus genome that underlie variation in resistance to blackleg, evaluated under field and shade-house conditions. Each of the SNP association accounted for up to 25.1% of additive genetic variance in resistance among diverse panel of accessions. To understand the homology of QR genomic regions with Arabidopsis thaliana genome, we searched the synteny between QR regions with 22 ancestral blocks of Brassicaceae. Comparative analyses revealed that 25 SNP associations for QR were localized in nine ancestral blocks, as a result of genomic rearrangements. We further showed that phenological traits such as flowering time, plant height, and maturity confound the genetic variation in resistance. Altogether, these findings provided new insights on the complex genetic control of the blackleg resistance and further expanded our understanding of its genetic architecture.

Identification of resistance loci in Chinese and Canadian canola/rapeseed varieties against Leptosphaeria maculans based on genome-wide association studies

Background

The fungal pathogen Leptosphaeria maculans (Lm). causes blackleg disease on canola/rapeseed in many parts of the world. It is important to use resistant cultivars to manage the disease and minimize yield losses. In this study, twenty-two Lm isolates were used to identify resistance genes in a collection of 243 canola/rapeseed (Brassica napus L.) accessions from Canada and China. These Lm isolates carry different compliments of avirulence genes, and the investigation was based on a genome-wide association study (GWAS) and genotype-by-sequencing (GBS).

Results

Using the CROP-SNP pipeline, a total of 81,471 variants, including 78,632 SNPs and 2839 InDels, were identified. The GWAS was performed using TASSEL 5.0 with GLM + Q model. Thirty-two and 13 SNPs were identified from the Canadian and Chinese accessions, respectively, tightly associated with blackleg resistance with P values < 1 × 10^− 4. These SNP loci were distributed on chromosomes A03, A05, A08, A09, C01, C04, C05, and C07, with the majority of them on A08 followed by A09 and A03. The significant SNPs identified on A08 were all located in a 2010-kb region and associated with resistance to 12 of the 22 Lm isolates. Furthermore, 25 resistance gene analogues (RGAs) were identified in these regions, including two nucleotide binding site (NBS) domain proteins, fourteen RLKs, three RLPs and six TM-CCs. These RGAs can be the potential candidate genes for blackleg resistance.

Conclusion

This study provides insights into potentially new genomic regions for discovery of additional blackleg resistance genes. The identified regions associated with blackleg resistance in the germplasm collection may also contribute directly to the development of canola varieties with novel resistance genes against blackleg of canola.

Genetic and physical mapping of loci for resistance to blackleg disease in canola (Brassica napus L.)

Sustainable canola production is essential to meet growing human demands for vegetable oil, biodiesel, and meal for stock feed markets. Blackleg, caused by the fungal pathogen, Leptosphaeria maculans is a devastating disease that can lead to significant yield loss in many canola production regions worldwide. Breakdown of race-specific resistance to L. maculans in commercial cultivars poses a constant threat to the canola industry. To identify new alleles, especially for quantitative resistance (QR), we analyzed 177 doubled haploid (DH) lines derived from an RP04/Ag-Outback cross. DH lines were evaluated for QR under field conditions in three experiments conducted at Wagga Wagga (2013, 2014) and Lake Green (2015), and under shade house conditions using the 'ascospore shower' test. DH lines were also characterized for qualitative R gene-mediated resistance via cotyledon tests with two differential single spore isolates, IBCN17 and IBCN76, under glasshouse conditions. Based on 18,851 DArTseq markers, a linkage map representing 2,019 unique marker bins was constructed and then utilized for QTL detection. Marker regression analysis identified 22 significant marker associations for resistance, allowing identification of two race-specific resistance R genes, Rlm3 and Rlm4, and 21 marker associations for QR loci. At least three SNP associations for QR were repeatedly detected on chromosomes A03, A07 and C04 across phenotyping environments. Physical mapping of markers linked with these consistent QR loci on the B. napus genome assembly revealed their localization in close proximity of the candidate genes of B. napus BnaA03g26760D (A03), BnaA07g20240D (A07) and BnaC04g02040D (C04). Annotation of these candidate genes revealed their association with protein kinase and jumonji proteins implicated in defense resistance. Both Rlm3 and Rlm4 genes identified in this DH population did not show any association with resistance loci detected under either field and/or shade house conditions (ascospore shower) suggesting that both genes are ineffective in conferring resistance to L. maculans in Australian field conditions. Taken together, our study identified sequence-based molecular markers for dissecting R and QR loci to L. maculans in a canola DH population from the RP04/Ag-Outback cross.

A multiomics approach reveals the pivotal role of subcellular reallocation in determining rapeseed resistance to cadmium toxicity

Oilseed rape (Brassica napus) has great potential for phytoremediation of cadmium (Cd)-polluted soils due to its large plant biomass production and strong metal accumulation. Enhanced plant Cd resistance (PCR) is a crucial prerequisite for phytoremediation through hyper-accumulation of excess Cd. However, the complexity of the allotetraploid genome of rapeseed hinders our understanding of PCR. To explore rapeseed Cd-resistance mechanisms, we examined two genotypes, 'ZS11' (Cd-resistant) and 'W10' (Cd-sensitive), that exhibit contrasting PCR while having similar tissue Cd concentrations, and characterized their different fingerprints in terms of plant morphophysiology (electron microscopy), ion abundance (inductively coupled plasma mass spectrometry), DNA variation (whole-genome resequencing), transcriptomics (high-throughput mRNA sequencing), and metabolomics (ultra-high performance liquid chromatography-mass spectrometry). Fine isolation of cell components combined with ionomics revealed that more Cd accumulated in the shoot vacuoles and root pectins of the resistant genotype than in the sensitive one. Genome and transcriptome sequencing identified numerous DNA variants and differentially expressed genes involved in pectin modification, ion binding, and compartmentalization. Transcriptomics-assisted gene co-expression networks characterized BnaCn.ABCC3 and BnaA8.PME3 as the central members involved in the determination of rapeseed PCR. High-resolution metabolic profiles revealed greater accumulation of shoot Cd chelates, and stronger biosynthesis and higher demethylation of root pectins in the resistant genotype than in the sensitive one. Our comprehensive examination using a multiomics approach has greatly improved our understanding of the role of subcellular reallocation of Cd in the determination of PCR.

NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci

Genes containing nucleotide-binding sites (NBS) play an important role in pathogen resistance in plants. However, the evolutionary fate of NBS-encoding genes after formation of allotetraploid Brassica napus (A_nA_nC_nC_n, 2n = 38) is still unknown. We performed a genome-wide comparison of putatively functional NBS-encoding genes in B. napus and its progenitor species Brassica rapa (A_rA_r, 2n = 20) and Brassica oleracea (C_oC_o, 2n = 18), identifying 464, 202, and 146 putatively functional NBS-encoding genes respectively, with genes unevenly distributed in several clusters. The A_n-subgenome of B. napus possessed similar numbers of NBS-encoding genes (191 genes) to the A_r genome of B. rapa (202 genes) and similar clustering patterns. However, the C_n genome of B. napus had many more genes (273) than the B. oleracea C_o genome (146), with different clustering trends. Only 97 NBS-encoding genes (66.4%) in B. oleracea were homologous with NBS-encoding genes in B. napus, while 176 NBS-encoding genes (87.1%) were homologous between B. rapa and B. napus. These results suggest a greater diversification of NBS-encoding genes in the C genome may have occurred after formation of B. napus. Although most NBS-encoding genes in B. napus appeared to derive from the progenitors, the birth and death of several NBS-encoding genes was also putatively mediated by non-homologous recombination. The Ka/Ks values of most homologous pairs between B. napus and the progenitor species were less than 1, suggesting purifying selection during B. napus evolution. The majority of NBS-encoding genes (60% in all species) showed higher expression levels in root tissue (out of root, leaf, stem, seed and flower tissue types). Comparative analysis of NBS-encoding genes with mapped resistance QTL against three major diseases of B. napus (blackleg, clubroot and Sclerotinia stem rot) found 204 NBS-encoding genes in B. napus located within 71 resistance QTL intervals. The majority of NBS-encoding genes were co-located with resistance QTLs against a single disease, while 47 genes were co-located with QTLs against two diseases and 3 genes were co-located with QTLs against all three. Our results revealed significant variation as well as interesting evolutionary trajectories of NBS-encoding genes in the different Brassica subgenomes, while co-localization of NBS-encoding genes and resistance QTL may facilitate resistance breeding in oilseed rape.

Finding invisible quantitative trait loci with missing data

Evolutionary processes during plant polyploidization and speciation have led to extensive presence-absence variation (PAV) in crop genomes, and there is increasing evidence that PAV associates with important traits. Today, high-resolution genetic analysis in major crops frequently implements simple, cost-effective, high-throughput genotyping from single nucleotide polymorphism (SNP) hybridization arrays; however, these are normally not designed to distinguish PAV from failed SNP calls caused by hybridization artefacts. Here, we describe a strategy to recover valuable information from single nucleotide absence polymorphisms (SNaPs) by population-based quality filtering of SNP hybridization data to distinguish patterns associated with genuine deletions from those caused by technical failures. We reveal that including SNaPs in genetic analyses elucidate segregation of small to large-scale structural variants in nested association mapping populations of oilseed rape (Brassica napus), a recent polyploid crop with widespread structural variation. Including SNaP markers in genomewide association studies identified numerous quantitative trait loci, invisible using SNP markers alone, for resistance to two major fungal diseases of oilseed rape, Sclerotinia stem rot and blackleg disease. Our results indicate that PAV has a strong influence on quantitative disease resistance in B. napus and that SNaP analysis using cost-effective SNP array data can provide extensive added value from 'missing data'. This strategy might also be applicable for improving the precision of genetic mapping in many important crop species.

Stable Quantitative Resistance Loci to Blackleg Disease in Canola (Brassica napus L.) Over Continents

The hemibiotrophic fungus, Leptosphaeria maculans is the most devastating pathogen, causing blackleg disease in canola (Brassica napus L). To study the genomic regions involved in quantitative resistance (QR), 259-276 DH lines from Darmor-bzh/Yudal (DYDH) population were assessed for resistance to blackleg under shade house and field conditions across 3 years. In different experiments, the broad sense heritability varied from 43 to 95%. A total of 27 significant quantitative trait loci (QTL) for QR were detected on 12 chromosomes and explained between 2.14 and 10.13% of the genotypic variance. Of the significant QTL, at least seven were repeatedly detected across different experiments on chromosomes A02, A07, A09, A10, C01, and C09. Resistance alleles were mainly contributed by 'Darmor-bzh' but 'Yudal' also contributed few of them. Our results suggest that plant maturity and plant height may have a pleiotropic effect on QR in our conditions. We confirmed that Rlm9 which is present in 'Darmor-bzh' is not effective to confer resistance in our Australian field conditions. Comparative mapping showed that several R genes coding for nucleotide-binding leucine-rich repeat (LRR) receptors map in close proximity (within 200 Kb) of the significant trait-marker associations on the reference 'Darmor-bzh' genome assembly. More importantly, eight significant QTL regions were detected across diverse growing environments: Australia, France, and United Kingdom. These stable QTL identified herein can be utilized for enhancing QR in elite canola germplasm via marker- assisted or genomic selection strategies.

Brassicaceae Mustards: Traditional and Agronomic Uses in Australia and New Zealand

Commonly cultivated Brassicaceae mustards, namely garlic mustard (Alliaria petiolata), white mustard (Brassica alba), Ethiopian mustard (B. carinata), Asian mustard (B. juncea), oilseed rape (B. napus), black mustard (B. nigra), rapeseed (B. rapa), white ball mustard (Calepina irregularis), ball mustard (Neslia paniculata), treacle mustard (Erysimum repandum), hedge mustard (Sisymbrium officinale), Asian hedge mustard (S. orientale), smooth mustard (S. erysimoides) and canola are the major economically important oilseed crops in many countries. Mustards were naturalized to Australia and New Zealand and Australia is currently the second largest exporter of Brassicaceae oilseeds to meet the global demand for a healthy plant-derived oil, high in polyunsaturated fats. Apart from providing edible oil, various parts of these plants and many of their phytochemicals have been used traditionally for both agronomic as well as medicinal purposes, with evidence of their use by early Australian and New Zealand settlers and also the indigenous population. This review provides an overview of the current knowledge of traditional and agronomic uses of Brassicaceae oilseeds and mustards with a focus on their importance in Australia and New Zealand.

Genome-wide Association Study Identifies New Loci for Resistance to Leptosphaeria maculans in Canola

Key message "We identified both quantitative and quantitative resistance loci to Leptosphaeria maculans, a fungal pathogen, causing blackleg disease in canola. Several genome-wide significant associations were detected at known and new loci for blackleg resistance. We further validated statistically significant associations in four genetic mapping populations, demonstrating that GWAS marker loci are indeed associated with resistance to L. maculans. One of the novel loci identified for the first time, Rlm12, conveys adult plant resistance in canola." Blackleg, caused by Leptosphaeria maculans, is a significant disease which affects the sustainable production of canola (Brassica napus). This study reports a genome-wide association study based on 18,804 polymorphic SNPs to identify loci associated with qualitative and quantitative resistance to L. maculans. Genomic regions delimited with 694 significant SNP markers, that are associated with resistance evaluated using 12 single spore isolates and pathotypes from four canola stubble were identified. Several significant associations were detected at known disease resistance loci including in the vicinity of recently cloned Rlm2/LepR3 genes, and at new loci on chromosomes A01/C01, A02/C02, A03/C03, A05/C05, A06, A08, and A09. In addition, we validated statistically significant associations on A01, A07, and A10 in four genetic mapping populations, demonstrating that GWAS marker loci are indeed associated with resistance to L. maculans. One of the novel loci identified for the first time, Rlm12, conveys adult plant resistance and mapped within 13.2 kb from Arabidopsis R gene of TIR-NBS class. We showed that resistance loci are located in the vicinity of R genes of Arabidopsis thaliana and Brassica napus on the sequenced genome of B. napus cv. Darmor-bzh. Significantly associated SNP markers provide a valuable tool to enrich germplasm for favorable alleles in order to improve the level of resistance to L. maculans in canola.

Multi-environment QTL studies suggest a role for cysteine-rich protein kinase genes in quantitative resistance to blackleg disease in Brassica napus

BACKGROUND

Resistance to the blackleg disease of Brassica napus (canola/oilseed rape), caused by the hemibiotrophic fungal pathogen Leptosphaeria maculans, is determined by both race-specific resistance (R) genes and quantitative resistance loci (QTL), or adult-plant resistance (APR). While the introgression of R genes into breeding material is relatively simple, QTL are often detected sporadically, making them harder to capture in breeding programs. For the effective deployment of APR in crop varieties, resistance QTL need to have a reliable influence on phenotype in multiple environments and be well defined genetically to enable marker-assisted selection (MAS).

RESULTS

Doubled-haploid populations produced from the susceptible B. napus variety Topas and APR varieties AG-Castle and AV-Sapphire were analysed for resistance to blackleg in two locations over 3 and 4 years, respectively. Three stable QTL were detected in each population, with two loci appearing to be common to both APR varieties. Physical delineation of three QTL regions was sufficient to identify candidate defense-related genes, including a cluster of cysteine-rich receptor-like kinases contained within a 49 gene QTL interval on chromosome A01. Individual L. maculans isolates were used to define the physical intervals for the race-specific R genes Rlm3 and Rlm4 and to identify QTL common to both field studies and the cotyledon resistance response.

CONCLUSION

Through multi-environment QTL analysis we have identified and delineated four significant and stable QTL suitable for MAS of quantitative blackleg resistance in B. napus, and identified candidate genes which potentially play a role in quantitative defense responses to L. maculans.

Pterostilbene Is a Potential Candidate for Control of Blackleg in Canola

Two stilbenes, resveratrol and pterostilbene, exhibit antifungal activity against Leptosphaeria maculans, the fungal pathogen responsible for blackleg (stem canker) in canola (Brassica napus). In vitro studies on the effect of these stilbenes on L. maculans mycelial growth and conidia germination showed that pterostilbene is a potent fungicide and sporicide, but resveratrol only exerted minor inhibition on L. maculans. Cell viability of hyphae cultures was markedly reduced by pterostilbene and SYTOX green staining showed that cell membrane integrity was compromised. We demonstrate that pterostilbene exerts fungicidal activity across 10 different L. maculans isolates and the compound confers protection to the blackleg-susceptible canola cv. Westar seedlings. The potential of pterostilbene as a control agent against blackleg in canola is discussed.

Identification of environmentally stable QTL for resistance against Leptosphaeria maculans in oilseed rape (Brassica napus)

Six stable QTL for resistance against L. maculans (phoma stem canker) have been identified by QTL × environment interaction analysis using data from five winter oilseed rape field experiments. Phoma stem canker, caused by Leptosphaeria maculans, is a disease of worldwide importance on oilseed rape (Brassica napus). Quantitative trait loci (QTL)-mediated resistance against L. maculans in B. napus is considered to be race non-specific and potentially durable. Identification and evaluation of QTL for resistance to L. maculans is important for breeding oilseed rape cultivars with durable resistance. An oilseed rape mapping population was used to detect QTL for resistance against L. maculans in five winter oilseed rape field experiments under different environments. A total of 17 QTL involved in 'field' quantitative resistance against L. maculans were detected and collectively explained 51% of the phenotypic variation. The number of QTL detected in each experiment ranged from two to nine and individual QTL explained 2-25% of the phenotypic variation. QTL × environment interaction analysis suggested that six of these QTL were less sensitive to environmental factors, so they were considered to be stable QTL. Markers linked to these stable QTL will be valuable for selection to breed for effective resistance against L. maculans in different environments, which will contribute to sustainable management of the disease.

Genetic variation and inheritance of phytosterol and oil content in a doubled haploid population derived from the winter oilseed rape Sansibar × Oase cross

Identification of QTL for phytosterol content, oil content, fatty acids content, protein content of defatted meal, and seed weight by multiple interval mapping in a Brassica napus DH population. Phytosterols are minor seed constituents in oilseed rape which have recently drawn wide-interest from the food and nutrition industry due to their health benefit in lowering LDL cholesterol in humans. To understand the genetic basis of phytosterol content and its relationship with other seed quality traits in oilseed rape, QTL mapping was performed in a segregating DH population derived from the cross of two winter oilseed rape varieties, Sansibar and Oase, termed SODH population. Both parental lines are of canola quality which differ in phytosterol and oil content in seed. A genetic map was constructed for SODH population based on a total of 1638 markers organized in 23 linkage groups and covering a map length of 2350 cM with a mean marker interval of 2.0 cM. The SODH population and the parental lines were cultivated at six environments in Europe and were phenotyped for phytosterol content, oil content, fatty acids content, protein content of the defatted meal, and seed weight. Multiple interval mapping identified between one and six QTL for nine phytosterol traits, between two and six QTL for four fatty acids, five QTL for oil content, four QTL for protein content of defatted meal, and three QTL for seed weight. Colocalizations of QTL for different traits were more frequently observed than individual isolated QTL. Major QTL (R (2) ≥ 25%) were all located in the A genome, and the possible candidate genes were investigated by physical localization of the QTL to the reference genome sequence of Brassica rapa.

Linkage and regional association analysis reveal two new tightly-linked major-QTLs for pod number and seed number per pod in rapeseed (Brassica napus L.)

To facilitate the pseudochromosomes assembly and gene cloning in rapeseed, we developed a reference genetic population/map (named BnaZNF2) from two sequenced cultivars, Zhongshuang11 and No.73290, those exhibit significant differences in many traits, particularly yield components. The BnaZNF2 genetic map exhibited perfect collinearity with the physical map of B. napus, indicating its high quality. Comparative mapping revealed several genomic rearrangements between B. napus and B. rapa or B. oleracea. A total of eight and 16 QTLs were identified for pod number and seed number per pod, respectively, and of which three and five QTLs are identical to previously identified ones, whereas the other five and 11 are novel. Two new major QTL respectively for pod number and seed number per pod, qPN.A06-1 and qSN.A06-1 (R(2 )= 22.8% and 32.1%), were colocalised with opposite effects, and only qPN.A06-1 was confirmed and narrowed by regional association analysis to 180 kb including only 33 annotated genes. Conditional QTL analysis and subsequent NILs test indicated that tight linkage, rather than pleiotropy, was the genetic causation of their colocalisation. Our study demonstrates potential of this reference genetic population/map for precise QTL mapping and as a base for positional gene cloning in rapeseed.

Comparative genomic analysis of duplicated homoeologous regions involved in the resistance of Brassica napus to stem canker

All crop species are current or ancient polyploids. Following whole genome duplication, structural and functional modifications result in differential gene content or regulation in the duplicated regions, which can play a fundamental role in the diversification of genes underlying complex traits. We have investigated this issue in Brassica napus, a species with a highly duplicated genome, with the aim of studying the structural and functional organization of duplicated regions involved in quantitative resistance to stem canker, a disease caused by the fungal pathogen Leptosphaeria maculans. Genome-wide association analysis on two oilseed rape panels confirmed that duplicated regions of ancestral blocks E, J, R, U, and W were involved in resistance to stem canker. The structural analysis of the duplicated genomic regions showed a higher gene density on the A genome than on the C genome and a better collinearity between homoeologous regions than paralogous regions, as overall in the whole B. napus genome. The three ancestral sub-genomes were involved in the resistance to stem canker and the fractionation profile of the duplicated regions corresponded to what was expected from results on the B. napus progenitors. About 60% of the genes identified in these duplicated regions were single-copy genes while less than 5% were retained in all the duplicated copies of a given ancestral block. Genes retained in several copies were mainly involved in response to stress, signaling, or transcription regulation. Genes with resistance-associated markers were mainly retained in more than two copies. These results suggested that some genes underlying quantitative resistance to stem canker might be duplicated genes. Genes with a hydrolase activity that were retained in one copy or R-like genes might also account for resistance in some regions. Further analyses need to be conducted to indicate to what extent duplicated genes contribute to the expression of the resistance phenotype.

Construction of a high-density genetic map using specific length amplified fragment markers and identification of a quantitative trait locus for anthracnose resistance in walnut (Juglans regia L.)

BACKGROUND

Walnut (Juglans regia, 2n = 32, approximately 606 Mb per 1C genome) is an economically important tree crop. Resistance to anthracnose, caused by Colletotrichum gloeosporioides, is a major objective of walnut genetic improvement in China. The recently developed specific length amplified fragment sequencing (SLAF-seq) is an efficient strategy that can obtain large numbers of markers with sufficient sequence information to construct high-density genetic maps and permits detection of quantitative trait loci (QTLs) for molecular breeding.

RESULTS

SLAF-seq generated 161.64 M paired-end reads. 153,820 SLAF markers were obtained, of which 49,174 were polymorphic. 13,635 polymorphic markers were sorted into five segregation types and 2,577 markers of them were used to construct genetic linkage maps: 2,395 of these fell into 16 linkage groups (LGs) for the female map, 448 markers for the male map, and 2,577 markers for the integrated map. Taking into account the size of all LGs, the marker coverage was 2,664.36 cM for the female map, 1,305.58 cM for the male map, and 2,457.82 cM for the integrated map. The average intervals between two adjacent mapped markers were 1.11 cM, 2.91 cM and 0.95 cM for three maps, respectively. 'SNP_only' markers accounted for 89.25% of the markers on the integrated map. Mapping markers contained 5,043 single nucleotide polymorphisms (SNPs) loci, which corresponded to two SNP loci per SLAF marker. According to the integrated map, we used interval mapping (Logarithm of odds, LOD > 3.0) to detect our quantitative trait. One QTL was detected for anthracnose resistance. The interval of this QTL ranged from 165.51 cM to 176.33 cM on LG14, and ten markers in this interval that were above the threshold value were considered to be linked markers to the anthracnose resistance trait. The phenotypic variance explained by each marker ranged from 16.2 to 19.9%, and their LOD scores varied from 3.22 to 4.04.

CONCLUSIONS

High-density genetic maps for walnut containing 16 LGs were constructed using the SLAF-seq method with an F1 population. One QTL for walnut anthracnose resistance was identified based on the map. The results will aid molecular marker-assisted breeding and walnut resistance genes identification.

Molecular marker applications in plants

Individuals within a population of a sexually reproducing species will have some degree of heritable genomic variation caused by mutations, insertion/deletions (INDELS), inversions, duplications, and translocations. Such variation can be detected and screened using molecular, or genetic, markers. By definition, molecular markers are genetic loci that can be easily tracked and quantified in a population and may be associated with a particular gene or trait of interest. This chapter will review the current major applications of molecular markers in plants.

Homoeologous duplicated regions are involved in quantitative resistance of Brassica napus to stem canker

BACKGROUND

Several major crop species are current or ancient polyploids. To better describe the genetic factors controlling traits of agronomic interest (QTL), it is necessary to understand the structural and functional organisation of these QTL regions in relation to genome duplication. We investigated quantitative resistance to the fungal disease stem canker in Brassica napus, a highly duplicated amphidiploid species, to assess the proportion of resistance QTL located at duplicated positions.

RESULTS

Genome-wide association analysis on a panel of 116 oilseed rape varieties genotyped with 3228 SNP indicated that 321 markers, corresponding to 64 genomic regions, are associated with resistance to stem canker. These genomic regions are relatively equally distributed on the A (53%) and C (47%) genomes of B. napus. Overall, 44% of these regions (28/64) are duplicated homoeologous regions. They are located in duplications of six (E, J, R, T, U and W) of the 24 ancestral blocks that constitute the B. napus genome. Overall, these six ancestral blocks have 34 duplicated copies in the B.napus genome. Almost all of the duplicated copies (82% of the 34 regions) harboured resistance associated markers for stem canker resistance, which suggests structural and functional conservation of genetic factors involved in this trait in B. napus.

CONCLUSIONS

Our study provides information on the involvement of duplicated loci in the control of stem canker resistance in B. napus. Further investigation of the similarity/divergence in sequence and gene content of these duplicated regions will provide insight into the conservation and allelic diversity of the underlying genes.

Identifying genetic diversity of avirulence genes in Leptosphaeria maculans using whole genome sequencing

Next generation sequencing technology allows rapid re-sequencing of individuals, as well as the discovery of single nucleotide polymorphisms (SNPs), for genomic diversity and evolutionary analyses. By sequencing two isolates of the fungal plant pathogen Leptosphaeria maculans, the causal agent of blackleg disease in Brassica crops, we have generated a resource of over 76 million sequence reads aligned to the reference genome. We identified over 21,000 SNPs with an overall SNP frequency of one SNP every 2,065 bp. Sequence validation of a selection of these SNPs in additional isolates collected throughout Australia indicates a high degree of polymorphism in the Australian population. In preliminary phylogenetic analysis, isolates from Western Australia clustered together and those collected from Brassica juncea stubble were identical. These SNPs provide a novel marker resource to study the genetic diversity of this pathogen. We demonstrate that re-sequencing provides a method of validating previously characterised SNPs and analysing differences in important genes, such as the disease related avirulence genes of L. maculans. Understanding the genetic characteristics of this devastating pathogen is vital in developing long-term solutions to managing blackleg disease in Brassica crops.

Development of a core set of single-locus SSR markers for allotetraploid rapeseed (Brassica napus L.)

Brassica napus (AACC) is a recent allotetraploid species evolved through hybridization between two diploids, B. rapa (AA) and B. oleracea (CC). Due to extensive genome duplication and homoeology within and between the A and C genomes of B. napus, most SSR markers display multiple fragments or loci, which limit their application in genetics and breeding studies of this economically important crop. In this study, we collected 3,890 SSR markers from previous studies and also developed 5,968 SSR markers from genomic sequences of B. rapa, B. oleracea and B. napus. Of these, 2,701 markers that produced single amplicons were putative single-locus markers in the B. napus genome. Finally, a set of 230 high-quality single-locus SSR markers were established and assigned to the 19 linkage groups of B. napus using a segregating population with 154 DH individuals. A subset of 78 selected single-locus SSR markers was proved to be highly stable and could successfully discriminate each of the 45 inbred lines and hybrids. In addition, most of the 230 SSR markers showed the single-locus nature in at least one of the Brassica species of the U's triangle besides B. napus. These results indicated that this set of single-locus SSR markers has a wide range of coverage with excellent stability and would be useful for gene tagging, sequence scaffold assignment, comparative mapping, diversity analysis, variety identification and association mapping in Brassica species.

High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napus L

Background

High density genetic maps built with SNP markers that are polymorphic in various genetic backgrounds are very useful for studying the genetics of agronomical traits as well as genome organization and evolution. Simultaneous dense SNP genotyping of segregating populations and variety collections was applied to oilseed rape (Brassica napus L.) to obtain a high density genetic map for this species and to study the linkage disequilibrium pattern.

Results

We developed an integrated genetic map for oilseed rape by high throughput SNP genotyping of four segregating doubled haploid populations. A very high level of collinearity was observed between the four individual maps and a large number of markers (>59%) was common to more than two maps. The precise integrated map comprises 5764 SNP and 1603 PCR markers. With a total genetic length of 2250 cM, the integrated map contains a density of 3.27 markers (2.56 SNP) per cM. Genotyping of these mapped SNP markers in oilseed rape collections allowed polymorphism level and linkage disequilibrium (LD) to be studied across the different collections (winter vs spring, different seed quality types) and along the linkage groups. Overall, polymorphism level was higher and LD decayed faster in spring than in “00” winter oilseed rape types but this was shown to vary greatly along the linkage groups.

Conclusions

Our study provides a valuable resource for further genetic studies using linkage or association mapping, for marker assisted breeding and for Brassica napus sequence assembly and genome organization analyses.

Identification and mapping of a novel blackleg resistance locus LepR4 in the progenies from Brassica napus × B. rapa subsp. sylvestris

Blackleg, caused by Leptosphaeria maculans, is one of the most economically important diseases of Brassica napus worldwide. Two blackleg-resistant lines, 16S and 61446, were developed through interspecific hybridization between B. napus and B. rapa subsp. sylvestris and backcrossing to B. napus. Classical genetic analysis demonstrated that a single recessive gene in both lines conferred resistance to L. maculans and that the resistance alleles were allelic. Using BC(1) progeny derived from each resistant plant, this locus was mapped to B. napus linkage group N6 and was flanked by microsatellite markers sN2189b and sORH72a in an interval of about 10 cM, in a region equivalent to about 6 Mb of B. rapa DNA sequence. This new resistance gene locus was designated as LepR4. The two lines were evaluated for resistance to a wide range of L. maculans isolates using cotyledon inoculation tests under controlled environment conditions, and for stem canker resistance in blackleg field nurseries. Results indicated that line 16S, carrying LepR4a, was highly resistant to all isolates tested on cotyledons and had a high level of stem canker resistance under field conditions. Line 61446, carrying LepR4b, was only resistant to some of the isolates tested on cotyledons and was weakly resistant to stem canker under field conditions.

Identification and characterization of candidate Rlm4 blackleg resistance genes in Brassica napus using next-generation sequencing

A thorough understanding of the relationships between plants and pathogens is essential if we are to continue to meet the agricultural needs of the world's growing population. The identification of genes underlying important quantitative trait loci is extremely challenging in complex genomes such as Brassica napus (canola, oilseed rape or rapeseed). However, recent advances in next-generation sequencing (NGS) enable much quicker identification of candidate genes for traits of interest. Here, we demonstrate this with the identification of candidate disease resistance genes from B. napus for its most devastating fungal pathogen, Leptosphaeria maculans (blackleg fungus). These two species are locked in an evolutionary arms race whereby a gene-for-gene interaction confers either resistance or susceptibility in the plant depending on the genotype of the plant and pathogen. Preliminary analysis of the complete genome sequence of Brassica rapa, the diploid progenitor of B. napus, identified numerous candidate genes with disease resistance characteristics, several of which were clustered around a region syntenic with a major locus (Rlm4) for blackleg resistance on A7 of B. napus. Molecular analyses of the candidate genes using B. napus NGS data are presented, and the difficulties associated with identifying functional gene copies within the highly duplicated Brassica genome are discussed.

Molecular mapping of qualitative and quantitative loci for resistance to Leptosphaeria maculans causing blackleg disease in canola (Brassica napus L.)

Blackleg, caused by Leptosphaeria maculans, is one of the most important diseases of oilseed and vegetable crucifiers worldwide. The present study describes (1) the construction of a genetic linkage map, comprising 255 markers, based upon simple sequence repeats (SSR), sequence-related amplified polymorphism, sequence tagged sites, and EST-SSRs and (2) the localization of qualitative (race-specific) and quantitative (race non-specific) trait loci controlling blackleg resistance in a doubled-haploid population derived from the Australian canola (Brassica napus L.) cultivars Skipton and Ag-Spectrum using the whole-genome average interval mapping approach. Marker regression analyses revealed that at least 14 genomic regions with LOD ≥ 2.0 were associated with qualitative and quantitative blackleg resistance, explaining 4.6-88.9 % of genotypic variation. A major qualitative locus, designated RlmSkipton (Rlm4), was mapped on chromosome A7, within 0.8 cM of the SSR marker Xbrms075. Alignment of the molecular markers underlying this QTL region with the genome sequence data of B. rapa L. suggests that RlmSkipton is located approximately 80 kb from the Xbrms075 locus. Molecular marker-RlmSkipton linkage was further validated in an F(2) population from Skipton/Ag-Spectrum. Our results show that SSR markers linked to consistent genomic regions are suitable for enrichment of favourable alleles for blackleg resistance in canola breeding programs.

Genomic advances will herald new insights into the Brassica: Leptosphaeria maculans pathosystem

The study of the relationship between plants and phytopathogenic fungi is one of the most rapidly moving fields in the plant sciences, the findings of which have contributed to the development of new strategies and technologies to protect crops. Plants employ sophisticated mechanisms to perceive and appropriately defend themselves against pathogens. A good example of plant and pathogen evolution is the gene-for-gene interaction between the fungal pathogen Leptosphaeria maculans, the causal agent of blackleg disease, and Brassica crops. This interaction has been studied at the genetic and physiological level due to its agro-economic importance. The newly available genome sequence for Brassica spp. and L. maculans will provide the resources to study the co-evolution of this plant and pathogen. Particularly, an understanding of the co-evolution of genes responsible for virulence and resistance will lead to improved plant protection strategies for Brassica canola and provide a model to understand plant-pathogen interactions in other major crops. This review summarises the research-to-date in the study of the Brassica-L. maculans gene-for-gene interaction, with a focus on the genetics of resistance in Brassica and the wealth of information to be gained from genome sequencing efforts.

Molecular mapping of Arabidopsis thaliana lipid-related orthologous genes in Brassica napus

Quantitative Trait Loci (QTL) for oil content has been previously analyzed in a SG-DH population from a cross between a Chinese cultivar and a European cultivar of Brassica napus. Eight QTL with additive and epistatic effects, and with environmental interactions were evaluated. Here we present an integrated linkage map of this population predominantly based on informative markers derived from Brassica sequences, including 249 orthologous A. thaliana genes, where nearly half (112) are acyl lipid metabolism related genes. Comparative genomic analysis between B. napus and A. thaliana revealed 33 colinearity regions. Each of the conserved A. thaliana segments is present two to six times in the B. napus genome. Approximately half of the mapped lipid-related orthologous gene loci (76/137) were assigned in these conserved colinearity regions. QTL analysis for seed oil content was performed using the new map and phenotypic data from 11 different field trials. Nine significant QTL were identified on linkage groups A1, A5, A7, A9, C2, C3, C6 and C8, together explaining 57.79% of the total phenotypic variation. A total of 14 lipid related candidate gene loci were located in the confidence intervals of six of these QTL, of which ten were assigned in the conserved colinearity regions and felled in the most frequently overlapped QTL intervals. The information obtained from this study demonstrates the potential role of the suggested candidate genes in rapeseed kernel oil accumulation.

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.

Large-scale development of functional markers in Brassica species

Numerous quantitative trait loci (QTL) have been detected in Brassica species, but fine-mapping of major QTL has advanced slowly. The development of functional markers can overcome this barrier. We used publicly available PlantGDB-assembled unique transcripts (PUTs) from Brassica species to design 7836 functional simple sequence repeat (SSR) primer pairs. Functional annotation of the PUTs containing SSRs was done by Blast2GO. The PUTs harbouring SSRs were mainly involved with nucleotide or protein binding and enzyme activity, and preferentially functioned in membranes and cytoplasm. Totally, 210 PUT primer pairs were selected to test their polymorphism, stability, and PCR quality. Approximately 70% (147) of the primer pairs resulted in successful amplification with an average polymorphic information content (PIC) value of 0.49. The highest level of polymorphism was dinucleotide repeat SSRs, followed by tri- and mononucleotide repeats. Approximately 60% of the primer pairs showed good transferability among Brassica species. These results show that the development of markers from PUTs is a feasible and simple approach to develop functional SSR markers on a large scale across Brassica species. In addition, these markers can provide a novel alternative that is a putative approach for rapid determination of candidate genes, genetic mapping, genetic diversity analysis, and comparative mapping in Brassica species.

Genic microsatellite markers in Brassica rapa: development, characterization, mapping, and their utility in other cultivated and wild Brassica relatives

Genic microsatellite markers, also known as functional markers, are preferred over anonymous markers as they reveal the variation in transcribed genes among individuals. In this study, we developed a total of 707 expressed sequence tag-derived simple sequence repeat markers (EST-SSRs) and used for development of a high-density integrated map using four individual mapping populations of B. rapa. This map contains a total of 1426 markers, consisting of 306 EST-SSRs, 153 intron polymorphic markers, 395 bacterial artificial chromosome-derived SSRs (BAC-SSRs), and 572 public SSRs and other markers covering a total distance of 1245.9 cM of the B. rapa genome. Analysis of allelic diversity in 24 B. rapa germplasm using 234 mapped EST-SSR markers showed amplification of 2 alleles by majority of EST-SSRs, although amplification of alleles ranging from 2 to 8 was found. Transferability analysis of 167 EST-SSRs in 35 species belonging to cultivated and wild brassica relatives showed 42.51% (Sysimprium leteum) to 100% (B. carinata, B. juncea, and B. napus) amplification. Our newly developed EST-SSRs and high-density linkage map based on highly transferable genic markers would facilitate the molecular mapping of quantitative trait loci and the positional cloning of specific genes, in addition to marker-assisted selection and comparative genomic studies of B. rapa with other related species.

Transcriptome sequencing of lentil based on second-generation technology permits large-scale unigene assembly and SSR marker discovery

Background

Lentil (Lens culinaris Medik.) is a cool-season grain legume which provides a rich source of protein for human consumption. In terms of genomic resources, lentil is relatively underdeveloped, in comparison to other Fabaceae species, with limited available data. There is hence a significant need to enhance such resources in order to identify novel genes and alleles for molecular breeding to increase crop productivity and quality.

Results

Tissue-specific cDNA samples from six distinct lentil genotypes were sequenced using Roche 454 GS-FLX Titanium technology, generating c. 1.38 × 10⁶ expressed sequence tags (ESTs). De novo assembly generated a total of 15,354 contigs and 68,715 singletons. The complete unigene set was sequence-analysed against genome drafts of the model legume species Medicago truncatula and Arabidopsis thaliana to identify 12,639, and 7,476 unique matches, respectively. When compared to the genome of Glycine max, a total of 20,419 unique hits were observed corresponding to c. 31% of the known gene space. A total of 25,592 lentil unigenes were subsequently annoated from GenBank. Simple sequence repeat (SSR)-containing ESTs were identified from consensus sequences and a total of 2,393 primer pairs were designed. A subset of 192 EST-SSR markers was screened for validation across a panel 12 cultivated lentil genotypes and one wild relative species. A total of 166 primer pairs obtained successful amplification, of which 47.5% detected genetic polymorphism.

Conclusions

A substantial collection of ESTs has been developed from sequence analysis of lentil genotypes using second-generation technology, permitting unigene definition across a broad range of functional categories. As well as providing resources for functional genomics studies, the unigene set has permitted significant enhancement of the number of publicly-available molecular genetic markers as tools for improvement of this species.

Genomics and bioinformatics resources for crop improvement

Recent remarkable innovations in platforms for omics-based research and application development provide crucial resources to promote research in model and applied plant species. A combinatorial approach using multiple omics platforms and integration of their outcomes is now an effective strategy for clarifying molecular systems integral to improving plant productivity. Furthermore, promotion of comparative genomics among model and applied plants allows us to grasp the biological properties of each species and to accelerate gene discovery and functional analyses of genes. Bioinformatics platforms and their associated databases are also essential for the effective design of approaches making the best use of genomic resources, including resource integration. We review recent advances in research platforms and resources in plant omics together with related databases and advances in technology.