Integration of linkage maps for the Amphidiploid Brassica napus and comparative mapping with Arabidopsis and Brassica rapa

Global transcriptome analysis reveals potential genes associated with genic male sterility of rapeseed (Brassica napus L.)

Rapeseed is the third leading source of edible oil in the world. Genic male sterility (GMS) lines provide crucial material for harnessing heterosis for rapeseed. GMS lines have been widely used successfully for rapeseed hybrid production. The physiological and molecular mechanism of pollen development in GMS lines of rapeseed (Brassica napus L.) need to be determined for the creation of hybrids and cultivation of new varieties. However, limited studies have focused on systematically mining genes that regulate the pollen development of GMS lines in B. napus. In the present study, to determine the stage at which pollen development begins to show abnormality in the GMS lines, we performed semi-thin section analysis of the anthers with five pollen development stages. The results indicated that the abnormal pollen development in DGMS lines might start at the meiotic stage, and abnormal pollen development in RGMS lines probably occurred before the tetrad stage. To investigate the critical genes and pathways involved in pollen development in GMS lines, we constructed and sequenced 24 transcriptome libraries for the flower buds from the fertile and sterile lines of two recessive GMS (RGMS) lines (6251AB and 6284AB) and two dominant GMS (DGMS) lines (4001AB and 4006AB). A total of 23,554 redundant DEGs with over two-fold change between sterile and fertile lines were obtained. A total of 346 DEGs were specifically related to DGMS, while 1,553 DEGs were specifically related to RGMS. A total of 1,545 DEGs were shared between DGMS and RGMS. And 253 transcription factors were found to be differentially expressed between the sterile and fertile lines of GMS. In addition, 6,099 DEGs possibly related to anther, pollen, and microspore development processes were identified. Many of these genes have been reported to be involved in anther and microspore developmental processes. Several DEGs were speculated to be key genes involved in the regulation of fertility. Three differentially expressed genes were randomly selected and their expression levels were verified by quantitative PCR (qRT-PCR). The results of qRT-PCR largely agreed with the transcriptome sequencing results. Our findings provide a global view of genes that are potentially involved in GMS occurrence. The expression profiles and function analysis of these DEGs were provided to expand our understanding of the complex molecular mechanism in pollen and sterility development in B. napus.

Mapping and identification of a new potential dominant resistance gene to turnip mosaic virus in Brassica rapa

By constructing an F₂ population, a new potential dominant resistance gene to TuMV in Brassica rapa was mapped and identified. Brassica rapa is the most widely grown vegetable crop in China, and turnip mosaic virus (TuMV) is a great threat to its production. Hence, it is a very important work to excavate more and novel resistance genes in B. rapa. In this study, the resistant line B80124 and the susceptible line B80450 were used to construct the F₂ populations, and through genetic analysis, the resistance to TuMV was found to be controlled by a dominant gene. Bulked segregant analysis sequence (BSA-seq) was used for the primary mapping, and an intersection (22.25-25.03 Mb) was obtained. After fine mapping using single nucleotide polymorphisms (SNP) markers, the candidate region was narrowed to 330 kb between the SNP markers A06S11 and A06S14, including eight genes relating to disease resistance. Using the transcriptome analysis and sequence identification, BraA06g035130.3C was screened as the final candidate gene, and it contained two deletion mutations, leading to frameshift in the susceptible line B80450. In addition, the phylogenetic analysis, hydrophilia and hydrophobicity analysis, subcellular location prediction analysis, amino acid bias analysis, and 3D modeling structures of BraA06g035130.3C were conducted to predict its functions. This study was conducive to the identification of a new TuMV resistance gene in B. rapa, which is of important scientific significance and application value for the improvement of TuMV resistance traits and molecular design breeding for Brassica crops.

Meiosis and beyond - understanding the mechanistic and evolutionary processes shaping the germline genome

The separation of germ cell populations from the soma is part of the evolutionary transition to multicellularity. Only genetic information present in the germ cells will be inherited by future generations, and any molecular processes affecting the germline genome are therefore likely to be passed on. Despite its prevalence across taxonomic kingdoms, we are only starting to understand details of the underlying micro-evolutionary processes occurring at the germline genome level. These include segregation, recombination, mutation and selection and can occur at any stage during germline differentiation and mitotic germline proliferation to meiosis and post-meiotic gamete maturation. Selection acting on germ cells at any stage from the diploid germ cell to the haploid gametes may cause significant deviations from Mendelian inheritance and may be more widespread than previously assumed. The mechanisms that affect and potentially alter the genomic sequence and allele frequencies in the germline are pivotal to our understanding of heritability. With the rise of new sequencing technologies, we are now able to address some of these unanswered questions. In this review, we comment on the most recent developments in this field and identify current gaps in our knowledge.

Identification of Chromosomes and Chromosome Rearrangements in Crop Brassicas and Raphanus sativus: A Cytogenetic Toolkit Using Synthesized Massive Oligonucleotide Libraries

Crop brassicas include three diploid [Brassica rapa (AA; 2n = 2x = 16), B. nigra (BB; 2n = 2x = 18), and B. oleracea (CC; 2n = 2x = 20)] and three derived allotetraploid species. It is difficult to distinguish Brassica chromosomes as they are small and morphologically similar. We aimed to develop a genome-sequence based cytogenetic toolkit for reproducible identification of Brassica chromosomes and their structural variations. A bioinformatic pipeline was used to extract repeat-free sequences from the whole genome assembly of B. rapa. Identified sequences were subsequently used to develop four c. 47-mer oligonucleotide libraries comprising 27,100, 11,084, 9,291, and 16,312 oligonucleotides. We selected these oligonucleotides after removing repeats from 18 identified sites (500-1,000 kb) with 1,997-5,420 oligonucleotides localized at each site in B. rapa. For one set of probes, a new method for amplification or immortalization of the library is described. oligonucleotide probes produced specific and reproducible in situ hybridization patterns for all chromosomes belonging to A, B, C, and R (Raphanus sativus) genomes. The probes were able to identify structural changes between the genomes, including translocations, fusions, and deletions. Furthermore, the probes were able to identify a structural translocation between a pak choi and turnip cultivar of B. rapa. Overall, the comparative chromosomal mapping helps understand the role of chromosome structural changes during genome evolution and speciation in the family Brassicaceae. The probes can also be used to identify chromosomes in aneuploids such as addition lines used for gene mapping, and to track transfer of chromosomes in hybridization and breeding programs.

Engineering Multiple Abiotic Stress Tolerance in Canola, Brassica napus

Impacts of climate change like global warming, drought, flooding, and other extreme events are posing severe challenges to global crop production. Contribution of Brassica napus towards the oilseed industry makes it an essential component of international trade and agroeconomics. Consequences from increasing occurrences of multiple abiotic stresses on this crop are leading to agroeconomic losses making it vital to endow B. napus crop with an ability to survive and maintain yield when faced with simultaneous exposure to multiple abiotic stresses. For an improved understanding of the stress sensing machinery, there is a need for analyzing regulatory pathways of multiple stress-responsive genes and other regulatory elements such as non-coding RNAs. However, our understanding of these pathways and their interactions in B. napus is far from complete. This review outlines the current knowledge of stress-responsive genes and their role in imparting multiple stress tolerance in B. napus. Analysis of network cross-talk through omics data mining is now making it possible to unravel the underlying complexity required for stress sensing and signaling in plants. Novel biotechnological approaches such as transgene-free genome editing and utilization of nanoparticles as gene delivery tools are also discussed. These can contribute to providing solutions for developing climate change resilient B. napus varieties with reduced regulatory limitations. The potential ability of synthetic biology to engineer and modify networks through fine-tuning of stress regulatory elements for plant responses to stress adaption is also highlighted.

Identification of Flower-Specific Promoters through Comparative Transcriptome Analysis in Brassica napus

Brassica napus (oilseed rape) is an economically important oil crop worldwide. Sclerotinia stem rot (SSR) caused by Sclerotinia sclerotiorum is a threat to oilseed rape production. Because the flower petals play pivotal roles in the SSR disease cycle, it is useful to express the resistance-related genes specifically in flowers to hinder further infection with S. sclerotiorum. To screen flower-specific promoters, we first analyzed the transcriptome data from 12 different tissues of the B. napus line ZS11. In total, 249 flower-specific candidate genes with high expression in petals were identified, and the expression patterns of 30 candidate genes were verified by quantitative real-time transcription-PCR (qRT-PCR) analysis. Furthermore, two novel flower-specific promoters (FSP046 and FSP061 promoter) were identified, and the tissue specificity and continuous expression in petals were determined in transgenic Arabidopsis thaliana with fusing the promoters to β-glucuronidase (GUS)-reporter gene. GUS staining, transcript expression pattern, and GUS activity analysis indicated that FSP046 and FSP061 promoter were strictly flower-specific promoters, and FSP046 promoter had a stronger activity. The two promoters were further confirmed to be able to direct GUS expression in B. napus flowers using transient expression system. The transcriptome data and the flower-specific promoters screened in the present study will benefit fundamental research for improving the agronomic traits as well as disease and pest control in a tissue-specific manner.

QTL and Candidate Gene Identification for Silique Length Based on High-Dense Genetic Map in Brassica napus L

Silique length (SL) is an important yield trait and positively correlates with seeds per silique and seed weight. In the present study, two double haploid (DH) populations, established from crosses Zhongshuang11 × R11 (ZR) and R1 × R2 (RR), containing 280 and 95 DH lines, respectively, were used to map quantitative trait loci (QTL) for SL. A high-dense genetic map from ZR population was constructed comprising 14,658 bins on 19 linkage groups, with map length of 2,198.85 cM and an average marker distance of 0.15 cM. Genetic linkage map from RR population was constructed by using 2,046 mapped markers anchored to 19 chromosomes with 2,217-cM map length and an average marker distance of 1.08 cM. Major QTL qSL_ZR_A09 and qSL_RR_A09b on A09 were identified from ZR and RR populations, respectively. Both QTL could be stably detected in four environments. QTL qSL_RR_A09b and qSL_ZR_A09 were located on 68.5-70.8 cM and 91.33-91.94 cM interval with R² values of 14.99-39.07% and 15.00-20.36% in RR and ZR populations, respectively. Based on the physical positions of single nucleotide polymorphism (SNP) markers flanking qSL_ZR_A09 and gene annotation in Arabidopsis, 26 genes were identified with SNP/Indel variation between parents and two genes (BnaA09g41180D and BnaA09g41380D) were selected as the candidate genes. Expression analysis further revealed BnaA09g41180D, encoding homologs of Arabidopsis fasciclin-like arabinogalactan proteins (FLA3), as the most promising candidate gene for qSL_ZR_A09. The QTL identification and candidate gene analysis will provide new insight into the genomic regions controlling SL in Brassica napus as well as candidate genes underlying the QTL.

Syntenic quantitative trait loci and genomic divergence for Sclerotinia resistance and flowering time in Brassica napus

Oilseed rape (Brassica napus) is an allotetraploid with two subgenomes descended from a common ancestor. Accordingly, its genome contains syntenic regions with many duplicate genes, some of which may have retained their original functions, whereas others may have diverged. Here, we mapped quantitative trait loci (QTL) for stem rot resistance (SRR), a disease caused by the fungus Sclerotinia sclerotiorum, and flowering time (FT) in a recombinant inbred line population. The population was genotyped using B. napus 60K single nucleotide polymorphism arrays and phenotyped in six (FT) and nine (SSR) experimental conditions or environments. In total, we detected 30 SRR QTL and 22 FT QTL and show that some of the major QTL associated with these two traits were co-localized, suggesting a genetic linkage between them. Two SRR QTL on chromosome A2 and two on chromosome C2 were shown to be syntenic, suggesting the functional conservation of these regions. We used the syntenic properties of the genomic regions to exclude genes for selection candidates responsible for QTL-associated traits. For example, 152 of the 185 genes could be excluded from a syntenic A2-C2 region. These findings will help to elucidate polyploid genomics in future studies, in addition to providing useful information for B. napus breeding programs.

Construction of a High-Density Genetic Map and Analysis of Seed-Related Traits Using Specific Length Amplified Fragment Sequencing for Cucurbita maxima

Seed traits are agronomically important for Cucurbita breeding, but the genes controlling seed size, seed weight and seed number have not been mapped in Cucurbita maxima (C. maxima). In this study, 100 F₂ individual derived from two parental lines, "2013-12" and "9-6", were applied to construct a 3,376.87-cM genetic map containing 20 linkage groups (LGs) with an average genetic distance of 0.47 cM using a total of 8,406 specific length amplified fragment (SLAF) markers in C. maxima. Ten quantitative trait loci (QTLs) of seed width (SW), seed length (SL) and hundred-seed weight (HSW) were identified using the composite interval mapping (CIM) method. The QTLs affecting SW, SL and HSW explained a maximum of 38.6%, 28.9% and 17.2% of the phenotypic variation and were detected in LG6, LG6 and LG17, respectively. To validate these results, an additional 150 F₂ individuals were used for QTL mapping of SW and SL with cleaved amplified polymorphic sequence (CAPS) markers. We found that two major QTLs, SL6-1 and SW6-1, could be detected in both SLAF-seq and CAPS markers in an overlapped region. Based on gene annotation and non-synonymous single-nucleotide polymorphisms (SNPs) in the major SWand SL-associated regions, we found that two genes encoding a VQ motif and an E3 ubiquitin-protein ligase may be candidate genes influencing SL, while an F-box and leucinerich repeat (LRR) domain-containing protein is the potential regulator for SW in C. maxima. This study provides the first high-density linkage map of C. maxima using SNPs developed by SLAF-seq technology, which is a powerful tool for associated mapping of important agronomic traits, map-based gene cloning and marker-assisted selection (MAS)-based breeding in C. maxima.

Construction of a high-density genetic linkage map and QTL mapping of oleic acid content and three agronomic traits in sunflower (Helianthus annuus L.) using specific-locus amplified fragment sequencing (SLAF-seq)

High-density genetic linkage maps are particularly important for quantitative trait loci (QTL) mapping, genome assembly, and marker-assisted selection (MAS) in plants. In this study, a high-density genetic linkage map of sunflower (Helianthus annuus L.) was constructed using an F₂ population generated from a cross between Helianthus annuus L. '86-1' and 'L-1-OL-1' via specific-locus amplified fragment sequencing (SLAF-seq). After sequence preprocessing, 530.50 M reads (105.60 Gb) were obtained that contained a total of 343,197 SLAFs, of which 39,589 were polymorphic. Of the polymorphic SLAFs, 6,136 were organized into a linkage map consisting of 17 linkage groups (LGs) spanning 2,221.86 cM, with an average genetic distance of 0.36 cM between SLAFs. Based on this high-density genetic map, QTL analysis was performed that focused on four sunflower phenotypic traits: oleic acid content (OAC), plant height (PH), head diameter (HD), and stem diameter (SD). Subsequently, for these four traits eight QTLs were detected that will likely be useful for increasing our understanding of genetic factors underlying these traits and for use in marker-assisted selection (MAS) for future sunflower breeding.

High-density multi-population consensus genetic linkage map for peach

Highly saturated genetic linkage maps are extremely helpful to breeders and are an essential prerequisite for many biological applications such as the identification of marker-trait associations, mapping quantitative trait loci (QTL), candidate gene identification, development of molecular markers for marker-assisted selection (MAS) and comparative genetic studies. Several high-density genetic maps, constructed using the 9K SNP peach array, are available for peach. However, each of these maps is based on a single mapping population and has limited use for QTL discovery and comparative studies. A consensus genetic linkage map developed from multiple populations provides not only a higher marker density and a greater genome coverage when compared to the individual maps, but also serves as a valuable tool for estimating genetic positions of unmapped markers. In this study, a previously developed linkage map from the cross between two peach cultivars 'Zin Dai' and 'Crimson Lady' (ZC2) was improved by genotyping additional progenies. In addition, a peach consensus map was developed based on the combination of the improved ZC2 genetic linkage map with three existing high-density genetic maps of peach and a reference map of Prunus. A total of 1,476 SNPs representing 351 unique marker positions were mapped across eight linkage groups on the ZC2 genetic map. The ZC2 linkage map spans 483.3 cM with an average distance between markers of 1.38 cM/marker. The MergeMap and LPmerge tools were used for the construction of a consensus map based on markers shared across five genetic linkage maps. The consensus linkage map contains a total of 3,092 molecular markers, consisting of 2,975 SNPs, 116 SSRs and 1 morphological marker associated with slow ripening in peach (SR). The consensus map provides valuable information on marker order and genetic position for QTL identification in peach and other genetic studies within Prunus and Rosaceae.

Highly dense linkage maps from 31 full-sibling families of turbot (Scophthalmus maximus) provide insights into recombination patterns and chromosome rearrangements throughout a newly refined genome assembly

Highly dense linkage maps enable positioning thousands of landmarks useful for anchoring the whole genome and for analysing genome properties. Turbot is the most important cultured flatfish worldwide and breeding programs in the fifth generation of selection are targeted to improve growth rate, obtain disease resistant broodstock and understand sex determination to control sex ratio. Using a Restriction-site Associated DNA approach, we genotyped 18,214 single nucleotide polymorphism in 1,268 turbot individuals from 31 full-sibling families. Individual linkage maps were combined to obtain a male, female and species consensus maps. The turbot consensus map contained 11,845 markers distributed across 22 linkage groups representing a total normalised length of 3,753.9 cM. The turbot genome was anchored to this map, and scaffolds representing 96% of the assembly were ordered and oriented to obtain the expected 22 megascaffolds according to its karyotype. Recombination rate was lower in males, especially around centromeres, and pairwise comparison of 44 individual maps suggested chromosome polymorphism at specific genomic regions. Genome comparison across flatfish provided new evidence on karyotype reorganisations occurring across the evolution of this fish group.

Interactions of WRKY15 and WRKY33 transcription factors and their roles in the resistance of oilseed rape to Sclerotinia infection

WRKY transcription factors are known to participate in the defence responses of higher plants. However, little is known about the roles of such proteins, especially regarding their functions in the resistance of oilseed rape (Brassica napus) to Sclerotinia sclerotiorum, a necrotrophic fungal pathogen that causes stem rot. In this study, we identified BnWRKY33 as a S. sclerotiorum-responsive gene that positively regulates resistance to this pathogen by enhancing the expression of genes involved in camalexin synthesis and genes regulated by salicylic acid (SA) and jasmonic acid (JA). We also identified a S. sclerotiorum-responsive region in the promoter of BnWRKY33, which we revealed to be a relatively conserved W-box region in the promoters of homologous genes in different species. Using this S. sclerotiorum-responsive region as bait in a yeast one-hybrid assay, we identified another WRKY transcription factor, BnWRKY15, and observed that both BnWRKY15 and BnWRKY33 could bind to this region. In addition, BnWRKY15 overexpression simultaneously increased the susceptibility of B. napus to S. sclerotiorum and down-regulated BnWRKY33 after different durations of infection. Furthermore, BnWRKY15, which contains a transcriptional repression domain, exhibited reduced transactivation ability and could reduce the transactivation ability of BnWRKY33 in Arabidopsis protoplast assays. Therefore, we suggest that the increased susceptibility of BnWRKY15-overexpressing plants results from reduced BnWRKY33 expression, which is due to the inhibition of BnWRKY33 transcriptional activation by BnWRKY15.

Genome-Wide Identification and Structural Analysis of bZIP Transcription Factor Genes in Brassica napus

The basic region/leucine zipper motif (bZIP) transcription factor family is one of the largest families of transcriptional regulators in plants. bZIP genes have been systematically characterized in some plants, but not in rapeseed (Brassica napus). In this study, we identified 247 BnbZIP genes in the rapeseed genome, which we classified into 10 subfamilies based on phylogenetic analysis of their deduced protein sequences. The BnbZIP genes were grouped into functional clades with Arabidopsis genes with similar putative functions, indicating functional conservation. Genome mapping analysis revealed that the BnbZIPs are distributed unevenly across all 19 chromosomes, and that some of these genes arose through whole-genome duplication and dispersed duplication events. All expression profiles of 247 bZIP genes were extracted from RNA-sequencing data obtained from 17 different B. napus ZS11 tissues with 42 various developmental stages. These genes exhibited different expression patterns in various tissues, revealing that these genes are differentially regulated. Our results provide a valuable foundation for functional dissection of the different BnbZIP homologs in B. napus and its parental lines and for molecular breeding studies of bZIP genes in B. napus.

Gene Introgression in Weeds Depends on Initial Gene Location in the Crop: Brassica napus-Raphanus raphanistrum Model

The effect of gene location within a crop genome on its transfer to a weed genome remains an open question for gene flow assessment. To elucidate this question, we analyzed advanced generations of intergeneric hybrids, derived from an initial pollination of known oilseed rape varieties (Brassica napus, AACC, 2n = 38) by a local population of wild radish (Raphanus raphanistrum, RrRr, 2n = 18). After five generations of recurrent pollination, 307 G5 plants with a chromosome number similar to wild radish were genotyped using 105 B. napus specific markers well distributed along the chromosomes. They revealed that 49.8% of G5 plants carried at least one B. napus genomic region. According to the frequency of B. napus markers (0-28%), four classes were defined: Class 1 (near zero frequency), with 75 markers covering ∼70% of oilseed rape genome; Class 2 (low frequency), with 20 markers located on 11 genomic regions; Class 3 (high frequency), with eight markers on three genomic regions; and Class 4 (higher frequency), with two adjacent markers detected on A10. Therefore, some regions of the oilseed rape genome are more prone than others to be introgressed into wild radish. Inheritance and growth of plant progeny revealed that genomic regions of oilseed rape could be stably introduced into wild radish and variably impact the plant fitness (plant height and seed number). Our results pinpoint that novel technologies enabling the targeted insertion of transgenes should select genomic regions that are less likely to be introgressed into the weed genome, thereby reducing gene flow.

Introducing the Brassica Information Portal: Towards integrating genotypic and phenotypic Brassica crop data

The Brassica Information Portal (BIP) is a centralised repository for brassica phenotypic data. The site hosts trait data associated with brassica research and breeding experiments conducted on brassica crops, that are used as oilseeds, vegetables, livestock forage and fodder and for biofuels. A key feature is the explicit management of meta-data describing the provenance and relationships between experimental plant materials, as well as trial design and trait descriptors. BIP is an open access and open source project, built on the schema of CropStoreDB, and as such can provide trait data management strategies for any crop data. A new user interface and programmatic submission/retrieval system helps to simplify data access for researchers, breeders and other end-users. BIP opens up the opportunity to apply integrative, cross-project analyses to data generated by the Brassica Research Community. Here, we present a short description of the current status of the repository.

Genetic dissection of seed oil and protein content and identification of networks associated with oil content in Brassica napus

High-density linkage maps can improve the precision of QTL localization. A high-density SNP-based linkage map containing 3207 markers covering 3072.7 cM of the Brassica napus genome was constructed in the KenC-8 × N53-2 (KNDH) population. A total of 67 and 38 QTLs for seed oil and protein content were identified with an average confidence interval of 5.26 and 4.38 cM, which could explain up to 22.24% and 27.48% of the phenotypic variation, respectively. Thirty-eight associated genomic regions from BSA overlapped with and/or narrowed the SOC-QTLs, further confirming the QTL mapping results based on the high-density linkage map. Potential candidates related to acyl-lipid and seed storage underlying SOC and SPC, respectively, were identified and analyzed, among which six were checked and showed expression differences between the two parents during different embryonic developmental periods. A large primary carbohydrate pathway based on potential candidates underlying SOC- and SPC-QTLs, and interaction networks based on potential candidates underlying SOC-QTLs, was constructed to dissect the complex mechanism based on metabolic and gene regulatory features, respectively. Accurate QTL mapping and potential candidates identified based on high-density linkage map and BSA analyses provide new insights into the complex genetic mechanism of oil and protein accumulation in the seeds of rapeseed.

Brassica napus DS-3, encoding a DELLA protein, negatively regulates stem elongation through gibberellin signaling pathway

Identification and characterization of a semi-dwarfing gene ds-3 encoding a mutant DELLA protein regulating plant height through gibberellin signaling pathway. Lodging is one of the most important factors causing severe yield loss in oilseed rape. Utilization of semi-dwarf varieties has been proved the most effective way to increase lodging resistance and yield in many crops. To develop semi-dwarf germplasm in oilseed rape, we identified a semi-dwarf mutant ds-3 which showed a reduced response to phytohormones gibberellins (GAs). Genetic analysis indicated the dwarfism was controlled by a single semi-dominant gene, ds-3. The DS-3 gene was mapped to a genomic region on chromosome C07, which is syntenic to the region of a previously identified semi-dwarf gene ds-1 (BnaA06.RGA). In this region, DS-3 (BnaC07.RGA) gene was identified to encode a DELLA protein that functions as a repressor in GA signaling pathway. A substitution of proline to leucine was identified in ds-3 in the conserved VHYNP motif, which is essential for GA-dependent interaction between gibberellin receptor GID1 and DELLA proteins. Segregation analysis in the F₂ population derived from the cross between ds-1 and ds-3 demonstrated that BnaA06.RGA displayed a stronger effect on plant height than BnaC07.RGA, indicating that different RGA genes may play different roles in stem elongation. In addition to BnaA06.RGA and BnaC07.RGA, two more RGA genes (BnaA09.RGA and BnaC09.RGA) were identified in the Brassica napus (B. napus) genome. Reverse-transcription polymerase chain reaction (RT-PCR) and yeast two-hybrid (Y2H) assays suggest that both BnaA09.RGA and BnaC09.RGA are transcribed in leaves and stems and can mediate GA signaling in vivo. These genes represent potential targets for screening ideal semi-dwarfing alleles for oilseed rape breeding.

Mapping X-Disease Phytoplasma Resistance in Prunus virginiana

Phytoplasmas such as "Candidatus Phytoplasma pruni," the causal agent of X-disease of stone fruits, lack detailed biological analysis. This has limited the understanding of plant resistance mechanisms. Chokecherry (Prunus virginiana L.) is a promising model to be used for the plant-phytoplasma interaction due to its documented ability to resist X-disease infection. A consensus chokecherry genetic map "Cho" was developed with JoinMap 4.0 by joining two parental maps. The new map contains a complete set of 16 linkage groups, spanning a genetic distance of 2,172 cM with an average marker density of 3.97 cM. Three significant quantitative trait loci (QTL) associated with X-disease resistance were identified contributing to a total of 45.9% of the phenotypic variation. This updated genetic linkage map and the identified QTL will provide the framework needed to facilitate molecular genetics, genomics, breeding, and biotechnology research concerning X-disease in chokecherry and other Prunus species.

A Major Locus for Manganese Tolerance Maps on Chromosome A09 in a Doubled Haploid Population of Brassica napus L

Soil acidity poses a major threat to productivity of several crops; mainly due to the prevalence of toxic levels of Al³⁺ and Mn²⁺. Crop productivity could be harnessed on acid soils via the development of plant varieties tolerant to phytotoxic levels of these cations. In this study, we investigated the extent of natural variation for Mn²⁺ tolerance among ten parental lines of the Australian and International canola mapping populations. Response to Mn²⁺ toxicity was measured on the bases of cotyledon chlorosis, shoot biomass, and leaf area in nutrient solution under control (9 μM of MnCl₂⋅4H₂O) and Mn treatment (125 μM of MnCl₂⋅4H₂O). Among parental lines, we selected Darmor-bzh and Yudal that showed significant and contrasting variation in Mn²⁺ tolerance to understand genetic control and identify the quantitative trait loci (QTL) underlying Mn²⁺ tolerance. We evaluated parental lines and their doubled haploid (DH) progenies (196 lines) derived from an F₁ cross, Darmor-bzh/Yudal for Mn²⁺ tolerance. Mn²⁺-tolerant genotypes had significantly higher shoot biomass and leaf area compared to Mn²⁺-sensitive genotypes. A genetic linkage map based on 7,805 DArTseq markers corresponding to 2,094 unique loci was constructed and further utilized for QTL identification. A major locus, BnMn²⁺.A09 was further mapped with a SNP marker, Bn-A09-p29012402 (LOD score of 34.6) accounting for most of the variation in Mn²⁺ tolerance on chromosome A09. This is the first report on the genomic localization of a Mn²⁺ tolerance locus in B. napus. Additionally, an ortholog of A. thaliana encoding for cation efflux facilitator transporter was located within 3,991 bp from significant SNP marker associated with BnMn²⁺.A09. A suite of genome sequence based markers (DArTseq and Illumina Infinium SNPs) flanking the BnMn²⁺.A09 locus would provide an invaluable tool for various molecular breeding applications to improve canola production and profitability on Mn²⁺ toxic soils.

Extraction of the Constituent Subgenomes of the Natural Allopolyploid Rapeseed (Brassica napus L.)

As the dynamic nature of progenitor genomes accompanies the speciation by interspecific hybridization, the extraction of the constituent subgenome(s) from a natural allopolyploid species of long history and then restitution of the progenitor(s) provides the unique opportunity to study the genome evolution and interplay. Herein, the A subgenome from the allotetraploid oilseed rape (Brassica napus L., AACC) was extracted through inducing the preferential elimination of C-subgenome chromosomes in intertribal crosses and the progenitor B. rapa was restituted (RBR). Then by crossing and backcrossing RBR with B. napus donor, the C subgenome was in situ dissected by adding each of its nine chromosomes to the extracted A subgenome and establishing the whole set of monosonic alien addition lines (MAALs). RBR from spring-type B. napus genotype "Oro" expressed a phenotype resembling some type of B. rapa never observed before, but showed a winter-type flowering habit. This RBR had weaker growth vigor and suffered more seriously from biotic and abiotic stresses compared with Oro. The phenotypes specific for these MAALs showed the location of the related genes on the particular C-subgenome chromosomes. These MAALs exhibited obviously different frequencies in homeologous pairing and transmission of additional C-subgenome chromosomes, which were associated with the distinct degrees of their relatedness, and even with the possible genetic regulation for meiotic pairing evolved in B. napus Finally, large scaffolds undetermined for sequence assembly of B. napus were anchored to specific C-subgenome chromosomes using MAALs.

Molecular and systems approaches towards drought-tolerant canola crops

1169 I. 1170 II. 1170 III. 1172 IV. 1176 V. 1181 VI. 1182 1183 References 1183 SUMMARY: Modern agriculture is facing multiple challenges including the necessity for a substantial increase in production to meet the needs of a burgeoning human population. Water shortage is a deleterious consequence of both population growth and climate change and is one of the most severe factors limiting global crop productivity. Brassica species, particularly canola varieties, are cultivated worldwide for edible oil, animal feed, and biodiesel, and suffer dramatic yield loss upon drought stress. The recent release of the Brassica napus genome supplies essential genetic information to facilitate identification of drought-related genes and provides new information for agricultural improvement in this species. Here we summarize current knowledge regarding drought responses of canola, including physiological and -omics effects of drought. We further discuss knowledge gained through translational biology based on discoveries in the closely related reference species Arabidopsis thaliana and through genetic strategies such as genome-wide association studies and analysis of natural variation. Knowledge of drought tolerance/resistance responses in canola together with research outcomes arising from new technologies and methodologies will inform novel strategies for improvement of drought tolerance and yield in this and other important crop species.

Microsatellite markers used for genome-wide association mapping of partial resistance to Sclerotinia sclerotiorum in a world collection of Brassica napus

The fungal pathogen Sclerotinia sclerotiorum causes stem rot of oilseed rape (Brassica napus) worldwide. In preparation for genome-wide association mapping (GWAM) of sclerotinia resistance in B. napus, 152 accessions from diverse geographical regions were screened with a single Canadian isolate, #321. Plants were inoculated by attaching mycelium plugs to the main stem at full flower. Lesion lengths measured 7, 14 and 21 days after inoculation were used to calculate the area under the disease progress curve (AUDPC). Depth of penetration was noted and used to calculate percent soft and collapsed lesions (% s + c). The two disease traits were highly correlated (r = 0.93). Partially resistant accessions (AUDPC <7 and % s + c <2) were identified primarily from South Korea and Japan with a few from Pakistan, China and Europe. Genotyping of accessions with 84 simple sequence repeat markers provided 690 polymorphic loci for GWAM. The general linear model in TASSEL best fitted the data when adjusted for population structure (STRUCTURE), GLM + Q. After correction for positive false discovery rate, 34 loci were significantly associated with both disease traits of which 21 alleles contributed to resistance, while the remaining enhanced susceptibility. The phenotypic variation explained by the loci ranged from 6 to 25 %. Five loci mapped to published quantitative trait loci conferring sclerotinia resistance in Chinese lines.

Population Genomic Analysis Reveals Differential Evolutionary Histories and Patterns of Diversity across Subgenomes and Subpopulations of Brassica napus L

The allotetraploid species Brassica napus L. is a global crop of major economic importance, providing canola oil (seed) and vegetables for human consumption and fodder and meal for livestock feed. Characterizing the genetic diversity present in the extant germplasm pool of B. napus is fundamental to better conserve, manage and utilize the genetic resources of this species. We used sequence-based genotyping to identify and genotype 30,881 SNPs in a diversity panel of 782 B. napus accessions, representing samples of winter and spring growth habits originating from 33 countries across Europe, Asia, and America. We detected strong population structure broadly concordant with growth habit and geography, and identified three major genetic groups: spring (SP), winter Europe (WE), and winter Asia (WA). Subpopulation-specific polymorphism patterns suggest enriched genetic diversity within the WA group and a smaller effective breeding population for the SP group compared to WE. Interestingly, the two subgenomes of B. napus appear to have different geographic origins, with phylogenetic analysis placing WE and WA as basal clades for the other subpopulations in the C and A subgenomes, respectively. Finally, we identified 16 genomic regions where the patterns of diversity differed markedly from the genome-wide average, several of which are suggestive of genomic inversions. The results obtained in this study constitute a valuable resource for worldwide breeding efforts and the genetic dissection and prediction of complex B. napus traits.

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.

Molecular Mapping and QTL for Expression Profiles of Flavonoid Genes in Brassica napus

Flavonoids are secondary metabolites that are extensively distributed in the plant kingdom and contribute to seed coat color formation in rapeseed. To decipher the genetic networks underlying flavonoid biosynthesis in rapeseed, we constructed a high-density genetic linkage map with 1089 polymorphic loci (including 464 SSR loci, 97 RAPD loci, 451 SRAP loci, and 75 IBP loci) using recombinant inbred lines (RILs). The map consists of 19 linkage groups and covers 2775 cM of the B. napus genome with an average distance of 2.54 cM between adjacent markers. We then performed expression quantitative trait locus (eQTL) analysis to detect transcript-level variation of 18 flavonoid biosynthesis pathway genes in the seeds of the 94 RILs. In total, 72 eQTLs were detected and found to be distributed among 15 different linkage groups that account for 4.11% to 52.70% of the phenotypic variance atrributed to each eQTL. Using a genetical genomics approach, four eQTL hotspots together harboring 28 eQTLs associated with 18 genes were found on chromosomes A03, A09, and C08 and had high levels of synteny with genome sequences of A. thaliana and Brassica species. Associated with the trans-eQTL hotspots on chromosomes A03, A09, and C08 were 5, 17, and 1 genes encoding transcription factors, suggesting that these genes have essential roles in the flavonoid biosynthesis pathway. Importantly, bZIP25, which is expressed specifically in seeds, MYC1, which controls flavonoid biosynthesis, and the R2R3-type gene MYB51, which is involved in the synthesis of secondary metabolites, were associated with the eQTL hotspots, and these genes might thus be involved in different flavonoid biosynthesis pathways in rapeseed. Hence, further studies of the functions of these genes will provide insight into the regulatory mechanism underlying flavonoid biosynthesis, and lay the foundation for elaborating the molecular mechanism of seed coat color formation in B. napus.

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.

High density linkage mapping of genomic and transcriptomic SNPs for synteny analysis and anchoring the genome sequence of chickpea

This study presents genome-wide discovery of SNPs through next generation sequencing of the genome of Cicer reticulatum. Mapping of the C. reticulatum sequenced reads onto the draft genome assembly of C. arietinum (desi chickpea) resulted in identification of 842,104 genomic SNPs which were utilized along with an additional 36,446 genic SNPs identified from transcriptome sequences of the aforementioned varieties. Two new chickpea Oligo Pool All (OPAs) each having 3,072 SNPs were designed and utilized for SNP genotyping of 129 Recombinant Inbred Lines (RILs). Using Illumina GoldenGate Technology genotyping data of 5,041 SNPs were generated and combined with the 1,673 marker data from previously published studies, to generate a high resolution linkage map. The map comprised of 6698 markers distributed on eight linkage groups spanning 1083.93 cM with an average inter-marker distance of 0.16 cM. Utility of the present map was demonstrated for improving the anchoring of the earlier reported draft genome sequence of desi chickpea by ~30% and that of kabuli chickpea by 18%. The genetic map reported in this study represents the most dense linkage map of chickpea , with the potential to facilitate efficient anchoring of the draft genome sequences of desi as well as kabuli chickpea varieties.

A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding

Consensus linkage maps are important tools in crop genomics. We have assembled a high-density tetraploid wheat consensus map by integrating 13 data sets from independent biparental populations involving durum wheat cultivars (Triticum turgidum ssp. durum), cultivated emmer (T. turgidum ssp. dicoccum) and their ancestor (wild emmer, T. turgidum ssp. dicoccoides). The consensus map harboured 30 144 markers (including 26 626 SNPs and 791 SSRs) half of which were present in at least two component maps. The final map spanned 2631 cM of all 14 durum wheat chromosomes and, differently from the individual component maps, all markers fell within the 14 linkage groups. Marker density per genetic distance unit peaked at centromeric regions, likely due to a combination of low recombination rate in the centromeric regions and even gene distribution along the chromosomes. Comparisons with bread wheat indicated fewer regions with recombination suppression, making this consensus map valuable for mapping in the A and B genomes of both durum and bread wheat. Sequence similarity analysis allowed us to relate mapped gene-derived SNPs to chromosome-specific transcripts. Dense patterns of homeologous relationships have been established between the A- and B-genome maps and between nonsyntenic homeologous chromosome regions as well, the latter tracing to ancient translocation events. The gene-based homeologous relationships are valuable to infer the map location of homeologs of target loci/QTLs. Because most SNP and SSR markers were previously mapped in bread wheat, this consensus map will facilitate a more effective integration and exploitation of genes and QTL for wheat breeding purposes.

Genome-wide identification of SSR and SNP markers from the non-heading Chinese cabbage for comparative genomic analyses

BACKGROUND

Non-heading Chinese cabbage (NHCC), belonging to Brassica, is an important leaf vegetable in Asia. Although genetic analyses have been performed through conventional selection and breeding efforts, the domestication history of NHCC and the genetics underlying its morphological diversity remain unclear. Thus, the reliable molecular markers representative of the whole genome are required for molecular-assisted selection in NHCC.

RESULTS

A total of 20,836 simple sequence repeats (SSRs) were detected in NHCC, containing repeat types from mononucleotide to nonanucleotide. The average density was 62.93 SSRs/Mb. In gene regions, 5,435 SSRs were identified in 4,569 genes. A total of 5,008 primer pairs were designed, and 74 were randomly selected for validation. Among these, 60 (81.08%) were polymorphic in 18 Cruciferae. The number of polymorphic bands ranged from two to five, with an average of 2.70 for each primer. The average values of the polymorphism information content, observed heterozygosity, Hardy-Weinberg equilibrium, and Shannon's information index were 0.2970, 0.4136, 0.5706, and 0.5885, respectively. Four clusters were classified according to the unweighted pair-group method with arithmetic average cluster analysis of 18 genotypes. In addition, a total of 1,228,979 single nucleotide polymorphisms (SNPs) were identified in the NHCC through a comparison with the genome of Chinese cabbage, and the average SNP density in the whole genome was 4.33/Kb. The number of SNPs ranged from 341,939 to 591,586 in the 10 accessions, and the average heterozygous SNPs ratio was ~42.53%. All analyses showed these markers were high quality and reliable. Therefore, they could be used in the construction of a linkage map and for genetic diversity studies for NHCC in future.

CONCLUSIONS

This is the first systematic and comprehensive analysis and identification of SSRs in NHCC and 17 species. The development of a large number of SNP and SSR markers was successfully achieved for NHCC. These novel markers are valuable for constructing genetic linkage maps, comparative genome analysis, quantitative trait locus (QTL) mapping, genome-wide association studies, and marker-assisted selection in NHCC breeding system research.

Digital gene expression analysis of gene expression differences within Brassica diploids and allopolyploids

BACKGROUND

Brassica includes many successfully cultivated crop species of polyploid origin, either by ancestral genome triplication or by hybridization between two diploid progenitors, displaying complex repetitive sequences and transposons. The U's triangle, which consists of three diploids and three amphidiploids, is optimal for the analysis of complicated genomes after polyploidization. Next-generation sequencing enables the transcriptome profiling of polyploids on a global scale.

RESULTS

We examined the gene expression patterns of three diploids (Brassica rapa, B. nigra, and B. oleracea) and three amphidiploids (B. napus, B. juncea, and B. carinata) via digital gene expression analysis. In total, the libraries generated between 5.7 and 6.1 million raw reads, and the clean tags of each library were mapped to 18547-21995 genes of B. rapa genome. The unambiguous tag-mapped genes in the libraries were compared. Moreover, the majority of differentially expressed genes (DEGs) were explored among diploids as well as between diploids and amphidiploids. Gene ontological analysis was performed to functionally categorize these DEGs into different classes. The Kyoto Encyclopedia of Genes and Genomes analysis was performed to assign these DEGs into approximately 120 pathways, among which the metabolic pathway, biosynthesis of secondary metabolites, and peroxisomal pathway were enriched. The non-additive genes in Brassica amphidiploids were analyzed, and the results indicated that orthologous genes in polyploids are frequently expressed in a non-additive pattern. Methyltransferase genes showed differential expression pattern in Brassica species.

CONCLUSION

Our results provided an understanding of the transcriptome complexity of natural Brassica species. The gene expression changes in diploids and allopolyploids may help elucidate the morphological and physiological differences among Brassica species.

Construction of a high-density integrated genetic linkage map of rubber tree (Hevea brasiliensis) using genotyping-by-sequencing (GBS)

Construction of linkage maps is crucial for genetic studies and marker-assisted breeding programs. Recent advances in next generation sequencing technologies allow for the generation of high-density linkage maps, especially in non-model species lacking extensive genomic resources. Here, we constructed a high-density integrated genetic linkage map of rubber tree (Hevea brasiliensis), the sole commercial producer of high-quality natural rubber. We applied a genotyping-by-sequencing (GBS) technique to simultaneously discover and genotype single nucleotide polymorphism (SNP) markers in two rubber tree populations. A total of 21,353 single nucleotide substitutions were identified, 55% of which represented transition events. GBS-based genetic maps of populations P and C comprised 1704 and 1719 markers and encompassed 2041 cM and 1874 cM, respectively. The average marker densities of these two maps were one SNP in 1.23-1.25 cM. A total of 1114 shared SNP markers were used to merge the two component maps. An integrated linkage map consisted of 2321 markers and spanned the cumulative length of 2052 cM. The composite map showed a substantial improvement in marker density, with one SNP marker in every 0.89 cM. To our knowledge, this is the most saturated genetic map in rubber tree to date. This integrated map allowed us to anchor 28,965 contigs, covering 135 Mb or 12% of the published rubber tree genome. We demonstrated that GBS is a robust and cost-effective approach for generating a common set of genome-wide SNP data suitable for constructing integrated linkage maps from multiple populations in a highly heterozygous agricultural species.

Assessing and broadening genetic diversity of a rapeseed germplasm collection

Assessing the level of genetic diversity within a germplasm collection contributes to evaluating the potential for its utilization as a gene pool to improve the performance of cultivars. In this study, 45 high-quality simple sequence repeat (SSR) markers were screened and used to estimate the genetic base of a world-wide collection of 248 rapeseed (Brassica napus) inbred lines. For the whole collection, the genetic diversity of A genome was higher than that of C genome. The genetic diversity of C genome for the semi-winter type was the lowest among the different germplasm types. Because B. oleracea is usually used to broaden the genetic diversity of C genome in rapeseed, we evaluated the potential of 25 wild B. oleracea lines. More allelic variations and a higher genetic diversity were observed in B. oleracea than in rapeseed. One B. oleracea line and one oilseed B. rapa line were used to generate a resynthesized Brassica napus line, which was then crossed with six semi-winter rapeseed cultivars to produce 7 F1 hybrids. Not only the allele introgression but also mutations were observed in the hybrids, resulting in significant improvement of the genetic base.

A complex recombination pattern in the genome of allotetraploid Brassica napus as revealed by a high-density genetic map

Polyploidy plays a crucial role in plant evolution. Brassica napus (2n = 38, AACC), the most important oil crop in the Brassica genus, is an allotetraploid that originated through natural doubling of chromosomes after the hybridization of its progenitor species, B. rapa (2n = 20, AA) and B. oleracea (2n = 18, CC). A better understanding of the evolutionary relationship between B. napus and B. rapa, B. oleracea, as well as Arabidopsis, which has a common ancestor with these three species, will provide valuable information about the generation and evolution of allopolyploidy. Based on a high-density genetic map with single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers, we performed a comparative genomic analysis of B. napus with Arabidopsis and its progenitor species B. rapa and B. oleracea. Based on the collinear relationship of B. rapa and B. oleracea in the B. napus genetic map, the B. napus genome was found to consist of 70.1% of the skeleton components of the chromosomes of B. rapa and B. oleracea, with 17.7% of sequences derived from reciprocal translocation between homoeologous chromosomes between the A- and C-genome and 3.6% of sequences derived from reciprocal translocation between non-homologous chromosomes at both intra- and inter-genomic levels. The current study thus provides insights into the formation and evolution of the allotetraploid B. napus genome, which will allow for more accurate transfer of genomic information from B. rapa, B. oleracea and Arabidopsis to B. napus.

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

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

Sixteen cytosolic glutamine synthetase genes identified in the Brassica napus L. genome are differentially regulated depending on nitrogen regimes and leaf senescence

A total of 16 BnaGLN1 genes coding for cytosolic glutamine synthetase isoforms (EC 6.3.1.2.) were found in the Brassica napus genome. The total number of BnaGLN1 genes, their phylogenetic relationships, and genetic locations are in agreement with the evolutionary history of Brassica species. Two BnaGLN1.1, two BnaGLN1.2, six BnaGLN1.3, four BnaGLN1.4, and two BnaGLN1.5 genes were found and named according to the standardized nomenclature for the Brassica genus. Gene expression showed conserved responses to nitrogen availability and leaf senescence among the Brassiceae tribe. The BnaGLN1.1 and BnaGLN1.4 families are overexpressed during leaf senescence and in response to nitrogen limitation. The BnaGLN1.2 family is up-regulated under high nitrogen regimes. The members of the BnaGLN1.3 family are not affected by nitrogen availability and are more expressed in stems than in leaves. Expression of the two BnaGLN1.5 genes is almost undetectable in vegetative tissues. Regulations arising from plant interactions with their environment (such as nitrogen resources), final architecture, and therefore sink-source relations in planta, seem to be globally conserved between Arabidopsis and B. napus. Similarities of the coding sequence (CDS) and protein sequences, expression profiles, response to nitrogen availability, and ageing suggest that the roles of the different GLN1 families have been conserved among the Brassiceae tribe. These findings are encouraging the transfer of knowledge from the Arabidopsis model plant to the B. napus crop plant. They are of special interest when considering the role of glutamine synthetase in crop yield and grain quality in maize and wheat.

Comparative genomics of Brassicaceae crops

The family Brassicaceae is one of the major groups of the plant kingdom and comprises diverse species of great economic, agronomic and scientific importance, including the model plant Arabidopsis. The sequencing of the Arabidopsis genome has revolutionized our knowledge in the field of plant biology and provides a foundation in genomics and comparative biology. Genomic resources have been utilized in Brassica for diversity analyses, construction of genetic maps and identification of agronomic traits. In Brassicaceae, comparative sequence analysis across the species has been utilized to understand genome structure, evolution and the detection of conserved genomic segments. In this review, we focus on the progress made in genetic resource development, genome sequencing and comparative mapping in Brassica and related species. The utilization of genomic resources and next-generation sequencing approaches in improvement of Brassica crops is also discussed.

A combined linkage and regional association mapping validation and fine mapping of two major pleiotropic QTLs for seed weight and silique length in rapeseed (Brassica napus L.)

BACKGROUND

Seed weight (SW) and silique length (SL) are important determinants of the yield potential in rapeseed (Brassica napus L.). However, the genetic basis of both traits is poorly understood. The main objectives of this study were to dissect the genetic basis of SW and SL in rapeseed through the preliminary mapping of quantitative trait locus (QTL) by linkage analysis and fine mapping of the target major QTL by regional association analysis.

RESULTS

Preliminary linkage mapping identified thirteen and nine consensus QTLs for SW and SL, respectively. These QTLs explained 0.7-67.1% and 2.1-54.4% of the phenotypic variance for SW and SL, respectively. Of these QTLs, three pairs of SW and SL QTLs were co-localized and integrated into three unique QTLs. In addition, the significance level and genetic effect of the three co-localized QTLs for both SW and SL showed great variation before and after the conditional analysis. Moreover, the allelic effects of the three QTLs for SW were highly consistent with those for SL. Two of the three co-localized QTLs, uq.A09-1 (mean R(2) = 20.1% and 19.0% for SW and SL, respectively) and uq.A09-3 (mean R(2) = 13.5% and 13.2% for SW and SL, respectively), were detected in all four environments and showed the opposite additive-effect direction. These QTLs were validated and fine mapped (their confidence intervals were narrowed down from 5.3 cM to 1 cM for uq.A09-1 and 13.2 cM to 2.5 cM for uq.A09-3) by regional association analysis with a panel of 576 inbred lines, which has a relatively rapid linkage disequilibrium decay (0.3 Mb) in the target QTL region.

CONCLUSIONS

A few QTLs with major effects and several QTLs with moderate effects might contribute to the natural variation of SW and SL in rapeseed. The meta-, conditional and allelic effect analyses suggested that pleiotropy, rather than tight linkage, was the genetic basis of the three pairs of co-localized of SW and SL QTLs. Regional association analysis was an effective and highly efficient strategy for the direct fine mapping of target major QTL identified by preliminary linkage mapping.

Cross-family translational genomics of abiotic stress-responsive genes between Arabidopsis and Medicago truncatula

Cross-species translation of genomic information may play a pivotal role in applying biological knowledge gained from relatively simple model system to other less studied, but related, genomes. The information of abiotic stress (ABS)-responsive genes in Arabidopsis was identified and translated into the legume model system, Medicago truncatula. Various data resources, such as TAIR/AtGI DB, expression profiles and literatures, were used to build a genome-wide list of ABS genes. tBlastX/BlastP similarity search tools and manual inspection of alignments were used to identify orthologous genes between the two genomes. A total of 1,377 genes were finally collected and classified into 18 functional criteria of gene ontology (GO). The data analysis according to the expression cues showed that there was substantial level of interaction among three major types (i.e., drought, salinity and cold stress) of abiotic stresses. In an attempt to translate the ABS genes between these two species, genomic locations for each gene were mapped using an in-house-developed comparative analysis platform. The comparative analysis revealed that fragmental colinearity, represented by only 37 synteny blocks, existed between Arabidopsis and M. truncatula. Based on the combination of E-value and alignment remarks, estimated translation rate was 60.2% for this cross-family translation. As a prelude of the functional comparative genomic approaches, in-silico gene network/interactome analyses were conducted to predict key components in the ABS responses, and one of the sub-networks was integrated with corresponding comparative map. The results demonstrated that core members of the sub-network were well aligned with previously reported ABS regulatory networks. Taken together, the results indicate that network-based integrative approaches of comparative and functional genomics are important to interpret and translate genomic information for complex traits such as abiotic stresses.

The fate of chromosomes and alleles in an allohexaploid Brassica population

Production of allohexaploid Brassica (2n = AABBCC) is a promising goal for plant breeders due to the potential for hybrid heterosis and useful allelic contributions from all three of the Brassica genomes present in the cultivated diploid (2n = AA, 2n = BB, 2n = CC) and allotetraploid (2n = AABB, 2n = AACC, and 2n = BBCC) crop species (canola, cabbages, mustards). We used high-throughput SNP molecular marker assays, flow cytometry, and fluorescent in situ hybridization (FISH) to characterize a population of putative allohexaploids derived from self-pollination of a hybrid from the novel cross (B. napus × B. carinata) × B. juncea to investigate whether fertile, stable allohexaploid Brassica can be produced. Allelic segregation in the A and C genomes generally followed Mendelian expectations for an F2 population, with minimal nonhomologous chromosome pairing. However, we detected no strong selection for complete 2n = AABBCC chromosome complements, with weak correlations between DNA content and fertility (r(2) = 0.11) and no correlation between missing chromosomes or chromosome segments and fertility. Investigation of next-generation progeny resulting from one highly fertile F2 plant using FISH revealed general maintenance of high chromosome numbers but severe distortions in karyotype, as evidenced by recombinant chromosomes and putative loss/duplication of A- and C-genome chromosome pairs. Our results show promise for the development of meiotically stable allohexaploid lines, but highlight the necessity of selection for 2n = AABBCC karyotypes.

The widespread crucifer species Cardamine flexuosa is an allotetraploid with a conserved subgenomic structure

The origin of Cardamine flexuosa (Wavy Bittercress) has been a conundrum for more than six decades. Here we identify its parental species, analyse its genome structure in comparison to parental genomes and describe intergenomic structural variations in C. flexuosa. Genomic in situ hybridization (GISH) and comparative chromosome painting (CCP) uncovered the parental genomes and the chromosome composition of C. flexuosa and its presumed diploid progenitors. Cardamine flexuosa is an allotetraploid (2n = 4x = 32), originating from two diploid species, Cardamine amara and Cardamine hirsuta (2n = 2x = 16). The two parental species display almost perfectly conserved chromosomal collinearity for seven out of the eight chromosomes. A 13 Mb pericentric inversion distinguishes chromosome CA1 from CH1. A comparative cytomolecular map was established for C. flexuosa by CCP/GISH. Whereas conserved chromosome collinearity between the C. amara and C. hirsuta subgenomes might have promoted intergenomic rearrangements through homeologous recombination, only one reciprocal translocation between two homeologues has occurred since the origin of C. flexuosa. The genome of C. flexuosa demonstrates that allopolyploids can maintain remarkably stable subgenomes over 10(4) -10(5)  yr throughout a wide distribution range. By contrast, the rRNA genes underwent genome-specific elimination towards a diploid-like number of loci.

Identification and mapping of a novel dominant resistance gene, TuRB07 to Turnip mosaic virus in Brassica rapa

A novel dominant resistance gene, TuRB07, was found to confer resistance to an isolate of TuMV strain C4 in B. rapa line VC1 and mapped on the top of chromosome A06. The inheritance of resistance to Turnip mosaic virus in Brassica rapa was investigated by crossing the resistant line, VC1 with the susceptible line, SR5, and genotyping and phenotyping diverse progenies derived from this cross. Both a doubled haploid population, VCS3M-DH, an F2 and two BC1 (F1 × VC1 and F1 × SR5) populations were created. Population tests revealed that the resistance to the TuMV C4 isolate in B. rapa is controlled by a single dominant gene. This resistance gene, TuRB07 was positioned on the top of linkage group A06 of the B. rapa genome through bulk segregation analysis and fine mapping recombinants in three doubled haploid- and one backcross population using microsatellite markers developed from BAC end sequences. Within the region between the two closely linked markers flanking TuRB07, H132A24-s1, and KS10960, in the Chiifu reference genome, two genes encoding nucleotide-binding site and leucine-rich repeat proteins with a coiled-coil motif (CC-NBS-LRR), Bra018862 and Bra018863 were identified as candidate resistance genes. The gene Bra018862 is truncated, but the gene Bra018863 has all the domains to function. Furthermore, the analysis of structural variation using resequencing data of VC1 and SR5 revealed that Bra018863 might be a functional gene because the gene has no structural variation in the resistant line VC1 when compared with Chiifu, whereas at the other NBS-LRR genes large deletions were identified in the resistant line. Allelic differences of Bra018863 were found between VC1 and SR5, supporting the notion that this gene is a putative candidate gene for the virus resistance.

A reference consensus genetic map for molecular markers and economically important traits in faba bean (Vicia faba L.)

BACKGROUND

Faba bean (Vicia faba L.) is among the earliest domesticated crops from the Near East. Today this legume is a key protein feed and food worldwide and continues to serve an important role in culinary traditions throughout Middle East, Mediterranean region, China and Ethiopia. Adapted to a wide range of soil types, the main faba bean breeding objectives are to improve yield, resistance to biotic and abiotic stresses, seed quality and other agronomic traits. Genomic approaches aimed at enhancing faba bean breeding programs require high-quality genetic linkage maps to facilitate quantitative trait locus analysis and gene tagging for use in a marker-assisted selection. The objective of this study was to construct a reference consensus map in faba bean by joining the information from the most relevant maps reported so far in this crop.

RESULTS

A combination of two approaches, increasing the number of anchor loci in diverse mapping populations and joining the corresponding genetic maps, was used to develop a reference consensus map in faba bean. The map was constructed from three main recombinant inbreed populations derived from four parental lines, incorporates 729 markers and is based on 69 common loci. It spans 4,602 cM with a range from 323 to 1041 loci in six main linkage groups or chromosomes, and an average marker density of one locus every 6 cM. Locus order is generally well maintained between the consensus map and the individual maps.

CONCLUSION

We have constructed a reliable and fairly dense consensus genetic linkage map that will serve as a basis for genomic approaches in faba bean research and breeding. The core map contains a larger number of markers than any previous individual map, covers existing gaps and achieves a wider coverage of the large faba bean genome as a whole. This tool can be used as a reference resource for studies in different genetic backgrounds, and provides a framework for transferring genetic information when using different marker technologies. Combined with syntenic approaches, the consensus map will increase marker density in selected genomic regions and will be useful for future faba bean molecular breeding applications.

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.

Identification of QTLs associated with oil content in a high-oil Brassica napus cultivar and construction of a high-density consensus map for QTLs comparison in B. napus

Increasing seed oil content is one of the most important goals in breeding of rapeseed (B. napus L.). To dissect the genetic basis of oil content in B. napus, a large and new double haploid (DH) population containing 348 lines was obtained from a cross between 'KenC-8' and 'N53-2', two varieties with >10% difference in seed oil content, and this population was named the KN DH population. A genetic linkage map consisting of 403 markers was constructed, which covered a total length of 1783.9 cM with an average marker interval of 4.4 cM. The KN DH population was phenotyped in eight natural environments and subjected to quantitative trait loci (QTL) analysis for oil content. A total of 63 identified QTLs explaining 2.64-17.88% of the phenotypic variation were identified, and these QTLs were further integrated into 24 consensus QTLs located on 11 chromosomes using meta-analysis. A high-density consensus map with 1335 marker loci was constructed by combining the KN DH map with seven other published maps based on the common markers. Of the 24 consensus QTLs in the KN DH population, 14 were new QTLs including five new QTLs in A genome and nine in C genome. The analysis revealed that a larger population with significant differences in oil content gave a higher power detecting new QTLs for oil content, and the construction of the consensus map provided a new clue for comparing the QTLs detected in different populations. These findings enriched our knowledge of QTLs for oil content and should be a potential in marker-assisted breeding of B. napus.

Characterization of the quantitative trait locus OilA1 for oil content in Brassica napus

Increasing seed oil content has become one of the most important breeding criteria in rapeseed (Brassica napus). However, oil content is a complex quantitative trait. QTL mapping in a double haploid population (SG population) emerging from a cross between a German (Sollux) and Chinese (Gaoyou) cultivars revealed one QTL for oil content on linkage group A1 (OilA1), which was mapped to a 17 cM genetic interval. To further validate and characterize the OilA1, we constructed a high-resolution map using B. rapa sequence resources and developed a set of near-isogenic lines (NILs) by employing a DH line SG-DH267 as donor and Chinese parent Gaoyou as recurrent background. The results showed highly conserved synteny order between B. rapa and B. napus within the linkage group A1 and revealed a possible centromere region between two markers ZAASA1-38 and NTP3 (2.5 cM). OilA1 was firstly validated by 250 BC5F2 plants and was confirmed in a 10.6 cM interval between the markers ZAASA1-47 and ZAASA1-77. Further substitution mapping was conducted by using two generations of QTL-NILs, 283 lines from eight BC5F3:4 families and 428 plants from six BC5F4 sub-NILs and thus narrowed the OilA1 interval to 6.9 cM and 4.3 cM (1.4 Mb), respectively. Field investigations with two replications using homozygous BC5F3:4 sister sub-NILs indicated that NILs, which carry a Sollux chromosome segment across the target region showed significant higher oil content (1.26 %, p < 0.001) than their sister NILs containing Gaoyou chromosome. The OilA1 locus is of particular interest for breeding purpose in China because 80 % of Chinese cultivars do not carry this desirable allele.

Construction of a genetic linkage map and QTL analysis of erucic acid content and glucosinolate components in yellow mustard (Sinapis alba L.)

BACKGROUND

Yellow mustard (Sinapis alba L.) is an important condiment crop for the spice trade in the world. It has lagged behind oilseed Brassica species in molecular marker development and application. Intron length polymorphism (ILP) markers are highly polymorphic, co-dominant and cost-effective. The cross-species applicability of ILP markers from Brassica species and Arabidopsis makes them possible to be used for genetic linkage mapping and further QTL analysis of agronomic traits in yellow mustard.

RESULTS

A total of 250 ILP and 14 SSR markers were mapped on 12 linkage groups and designated as Sal01-12 in yellow mustard. The constructed map covered a total genetic length of 890.4 cM with an average marker interval of 3.3 cM. The QTL for erucic content co-localized with the fatty acid elongase 1 (FAE1) gene on Sal03. The self-(in)compatibility gene was assigned to Sal08. The 4-hydroxybenzyl, 3-indolylmethyl and 4-hydroxy-3-indolylmethyl glucosinolate contents were each controlled by one major QTL, all of which were located on Sal02. Two QTLs, accounting for the respective 20.4% and 19.2% of the total variation of 2-hydroxy-3-butenyl glucosinolate content, were identified and mapped to Sal02 and Sal11. Comparative synteny analysis revealed that yellow mustard was phylogenetically related to Arabidopsis thaliana and had undergone extensive chromosomal rearrangements during speciation.

CONCLUSION

The linkage map based on ILP and SSR markers was constructed and used for QTL analysis of seed quality traits in yellow mustard. The markers tightly linked with the genes for different glucosinolate components will be used for marker-assisted selection and map-based cloning. The ILP markers and linkage map provide useful molecular tools for yellow mustard breeding.

An integrated genetic linkage map for white clover (Trifolium repens L.) with alignment to Medicago

BACKGROUND

White clover (Trifolium repens L.) is a temperate forage legume with an allotetraploid genome (2n=4×=32) estimated at 1093 Mb. Several linkage maps of various sizes, marker sources and completeness are available, however, no integrated map and marker set has explored consistency of linkage analysis among unrelated mapping populations. Such integrative analysis requires tools for homoeologue matching among populations. Development of these tools provides for a consistent framework map of the white clover genome, and facilitates in silico alignment with the model forage legume, Medicago truncatula.

RESULTS

This is the first report of integration of independent linkage maps in white clover, and adds to the literature on methyl filtered GeneThresher®-derived microsatellite (simple sequence repeat; SSR) markers for linkage mapping. Gene-targeted SSR markers were discovered in a GeneThresher® (TrGT) methyl-filtered database of 364,539 sequences, which yielded 15,647 SSR arrays. Primers were designed for 4,038 arrays and of these, 465 TrGT-SSR markers were used for parental consensus genetic linkage analysis in an F1 mapping population (MP2). This was merged with an EST-SSR consensus genetic map of an independent population (MP1), using markers to match homoeologues and develop a multi-population integrated map of the white clover genome. This integrated map (IM) includes 1109 loci based on 804 SSRs over 1274 cM, covering 97% of the genome at a moderate density of one locus per 1.2 cM. Eighteen candidate genes and one morphological marker were also placed on the IM. Despite being derived from disparate populations and marker sources, the component maps and the derived IM had consistent representations of the white clover genome for marker order and genetic length. In silico analysis at an E-value threshold of 1e-20 revealed substantial co-linearity with the Medicago truncatula genome, and indicates a translocation between T. repens groups 2 and 6 relative to M. truncatula.

CONCLUSIONS

This integrated genetic linkage analysis provides a consistent and comprehensive linkage analysis of the white clover genome, with alignment to a model forage legume. Associated marker locus information, particularly the homoeologue-specific markers, offers a new resource for forage legume research to enable genetic analysis and improvement of this forage and grassland species.

A consensus map of rapeseed (Brassica napus L.) based on diversity array technology markers: applications in genetic dissection of qualitative and quantitative traits

Background

Dense consensus genetic maps based on high-throughput genotyping platforms are valuable for making genetic gains in Brassica napus through quantitative trait locus identification, efficient predictive molecular breeding, and map-based gene cloning. This report describes the construction of the first B. napus consensus map consisting of a 1,359 anchored array based genotyping platform; Diversity Arrays Technology (DArT), and non-DArT markers from six populations originating from Australia, Canada, China and Europe. We aligned the B. napus DArT sequences with genomic scaffolds from Brassica rapa and Brassica oleracea, and identified DArT loci that showed linkage with qualitative and quantitative loci associated with agronomic traits.

Results

The integrated consensus map covered a total of 1,987.2 cM and represented all 19 chromosomes of the A and C genomes, with an average map density of one marker per 1.46 cM, corresponding to approximately 0.88 Mbp of the haploid genome. Through in silico physical mapping 2,457 out of 3,072 (80%) DArT clones were assigned to the genomic scaffolds of B. rapa (A genome) and B. oleracea (C genome). These were used to orientate the genetic consensus map with the chromosomal sequences. The DArT markers showed linkage with previously identified non-DArT markers associated with qualitative and quantitative trait loci for plant architecture, phenological components, seed and oil quality attributes, boron efficiency, sucrose transport, male sterility, and race-specific resistance to blackleg disease.

Conclusions

The DArT markers provide increased marker density across the B. napus genome. Most of the DArT markers represented on the current array were sequenced and aligned with the B. rapa and B. oleracea genomes, providing insight into the Brassica A and C genomes. This information can be utilised for comparative genomics and genomic evolution studies. In summary, this consensus map can be used to (i) integrate new generation markers such as SNP arrays and next generation sequencing data; (ii) anchor physical maps to facilitate assembly of B. napus genome sequences; and (iii) identify candidate genes underlying natural genetic variation for traits of interest.