High-density single nucleotide polymorphism (SNP) array mapping in Brassica oleracea: identification of QTL associated with carotenoid variation in broccoli florets

Genetic Diversity and Population Structure in Ethiopian Mustard (Brassica carinata A. Braun) as Revealed by Single Nucleotide Polymorphism Markers

Ethiopian mustard (Brassica carinata A. Braun) is currently one of the potential oilseeds dedicated to the production for biofuel and other bio-industrial applications. The crop is assumed to be native to Ethiopia where a number of diversified B. carinata germplasms are found and conserved ex situ. However, there is very limited information on the genetic diversity and population structure of the species. This study aimed to investigate the genetic diversity and population structure of B. carinata genotypes of different origins using high-throughput single nucleotide polymorphism (SNP) markers. We used Brassica 90K Illumina InfiniumTM SNP array for genotyping 90 B. carinata genotypes, and a total of 11,499 informative SNP markers were used for investigating the population structure and genetic diversity. The structure analysis, principal coordinate analysis (PcoA) and neighbor-joining tree analysis clustered the 90 B. carinata genotypes into two distinct subpopulations (Pop1 and Pop2). The majority of accessions (65%) were clustered in Pop1, mainly obtained from Oromia and South West Ethiopian People (SWEP) regions. Pop2 constituted dominantly of breeding lines and varieties, implying target selection contributed to the formation of distinct populations. Analysis of molecular variance (AMOVA) revealed a higher genetic variation (93%) within populations than between populations (7%), with low genetic differentiation (PhiPT = 0.07) and poor correlation between genetic and geographical distance (R = 0.02). This implies the presence of gene flow (Nm > 1) and weak geographical structure of accessions. Genetic diversity indices showed the presence of moderate genetic diversity in B. carinata populations with an average genetic diversity value (HE = 0.31) and polymorphism information content (PIC = 0.26). The findings of this study provide important and relevant information for future breeding and conservation efforts of B. carinata.

Mapping and Validation of BrGOLDEN: A Dominant Gene Regulating Carotenoid Accumulation in Brassica rapa

In plants, the accumulation of carotenoids can maintain the balance of the photosystem and improve crop nutritional quality. Therefore, the molecular mechanisms underlying carotenoid synthesis and accumulation should be further explored. In this study, carotenoid accumulation differed significantly among parental Brassica rapa. Genetic analysis was carried out using the golden inner leaf ‘1900264′ line and the light−yellow inner leaf ‘1900262′ line, showing that the golden inner leaf phenotype was controlled by a single dominant gene. Using bulked−segregant analysis sequencing, BraA09g007080.3C encoding the ORANGE protein was selected as a candidate gene. Sequence alignment revealed that a 4.67 kb long terminal repeat insertion in the third exon of the BrGOLDEN resulted in three alternatively spliced transcripts. The spatiotemporal expression results indicated that BrGOLDEN might regulate the expression levels of carotenoid−synthesis−related genes. After transforming BrGOLDEN into Arabidopsis thaliana, the seed−derived callus showed that BrGOLDEN_Ins and BrGOLDEN_Del lines presented a yellow color and the BrGOLDEN_Ldel line presented a transparent phenotype. In addition, using the yeast two−hybrid assay, BrGOLDEN_Ins, BrGOLDEN_Ldel, and Brgolden_wt exhibited strong interactions with BrPSY1, but BrGOLDEN_Del did not interact with BrPSY1 in the split−ubiquitin membrane system. In the secondary and 3D structure analysis, BrGOLDEN_Del was shown to have lost the PNFPSFIPFLPPL sequences at the 125 amino acid position, which resulted in the α−helices of BrGOLDEN_Del being disrupted, restricting the formation of the 3D structure and affecting the functions of the protein. These findings may provide new insights into the regulation of carotenoid synthesis in B. rapa.

Different Shades of Kale-Approaches to Analyze Kale Variety Interrelations

Brassica oleracea is a vegetable crop with an amazing morphological diversity. Among the various crops derived from B. oleracea, kale has been in the spotlight globally due to its various health-benefitting compounds and many different varieties. Knowledge of the existing genetic diversity is essential for the improved breeding of kale. Here, we analyze the interrelationships, population structures, and genetic diversity of 72 kale and cabbage varieties by extending our previous diversity analysis and evaluating the use of summed potential lengths of shared haplotypes (SPLoSH) as a new method for such analyses. To this end, we made use of the high-density Brassica 60K SNP array, analyzed SNPs included in an available Brassica genetic map, and used these resources to generate and evaluate the information from SPLoSH data. With our results we could consistently differentiate four groups of kale across all analyses: the curly kale varieties, Italian, American, and Russian varieties, as well as wild and cultivated types. The best results were achieved by using SPLoSH information, thus validating the use of this information in improving analyses of interrelations in kale. In conclusion, our definition of kale includes the curly varieties as the kales in a strict sense, regardless of their origin. These results contribute to a better understanding of the huge diversity of kale and its interrelations.

Development and application of the Faba_bean_130K targeted next-generation sequencing SNP genotyping platform based on transcriptome sequencing

KEY MESSAGE

Large-scale faba bean transcriptome data are available, and the first genotyping platform based on liquid-phase probe targeted capture technology was developed for genetic and molecular breeding studies. Faba bean (Vicia faba L., 2n = 12) is an important food legume crop that is widely grown for multiple uses worldwide. However, no reference genome is currently available due to its very large genome size (approximately 13 Gb) and limited single nucleotide polymorphism (SNP) markers as well as highly efficient genotyping tools have been reported for faba bean. In this study, 16.7 billion clean reads were obtained from transcriptome libraries of flowers and leaves of 102 global faba bean accessions. A total of 243,120 unigenes were de novo assembled and functionally annotated. Moreover, a total of 1,579,411 SNPs were identified and further filtered according to a selection pipeline to develop a high-throughput, flexible, low-cost Faba_bean_130K targeted next-generation sequencing (TNGS) genotyping platform. A set of 69 Chinese faba bean accessions were genotyped with the TNGS genotyping platform, and the average mapping rate of captured reads to reference transcripts was 93.14%, of which 53.23% were located in the targeted regions. The TNGS genotyping results were validated by Sanger sequencing and the average consistency rate reached 93.6%. Comprehensive population genetic analysis was performed on the 69 Chinese faba bean accessions and identified four genetic subgroups correlated with the geographic distribution. This study provides valuable genomic resources and a reliable genotyping tool that could be implemented in genetic and molecular breeding studies to accelerate new cultivar development and improvement in faba bean.

Advances in Genetics and Molecular Breeding of Broccoli

Broccoli (Brassica oleracea L. var. italica) is one of the most important vegetable crops cultivated worldwide. The market demand for broccoli is still increasing due to its richness in vitamins, anthocyanins, mineral substances, fiber, secondary metabolites and other nutrients. The famous secondary metabolites, glucosinolates, sulforaphane and selenium have protective effects against cancer. Significant progress has been made in fine-mapping and cloning genes that are responsible for important traits; this progress provides a foundation for marker-assisted selection (MAS) in broccoli breeding. Genetic engineering by the well-developed Agrobacterium tumefaciens-mediated transformation in broccoli has contributed to the improvement of quality; postharvest life; glucosinolate and sulforaphane content; and resistance to insects, pathogens and abiotic stresses. Here, we review recent progress in the genetics and molecular breeding of broccoli. Future perspectives for improving broccoli are also briefly discussed.

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.

SLAF-Based Construction of a High-Density Genetic Map and Its Application in QTL Mapping of Carotenoids Content in Citrus Fruit

Carotenoids are important antioxidant components in the human diet. To develop carotenoid-rich agricultural products by genetic intervention, understanding the genetic basis of carotenoids variation is essential. In this study, we constructed a high-density integrated genetic map with 3817 molecular markers using specific locus amplified fragment (SLAF) sequencing from a C. reticulata × P. trifoliata F₁ pseudotestcross population. A total of 17 significant quantitative trait loci (QTLs) distributed on Chromosomes (Chr) 2, 3, 5, 6, and 9 were detected to determine the carotenoid variation in the population. In particular, three QTL colocalizations for multiple carotenoid constituents were observed on Chr 2, 3, and 9, one of which was located on Chr2:34,654,608-35430715 accounted for 20.1-25.4% of the variation of luteoxanthin, auroxanthin, lutein, violaxanthin, and total carotenoid content. Overall, this study provides a genetic foundation for marker-assisted selection (MAS) breeding of nutritionally enhanced citrus fruit.

Complex Horticultural Quality Traits in Broccoli Are Illuminated by Evaluation of the Immortal BolTBDH Mapping Population

Improving horticultural quality in regionally adapted broccoli (Brassica oleracea var. italica) and other B. oleracea crops is challenging due to complex genetic control of traits affecting morphology, development, and yield. Mapping horticultural quality traits to genomic loci is an essential step in these improvement efforts. Understanding the mechanisms underlying horticultural quality enables multi-trait marker-assisted selection for improved, resilient, and regionally adapted B. oleracea germplasm. The publicly-available biparental double-haploid BolTBDH mapping population (Chinese kale × broccoli; N = 175) was evaluated for 25 horticultural traits in six trait classes (architecture, biomass, phenology, leaf morphology, floral morphology, and head quality) by multiple quantitative trait loci mapping using 1,881 genotype-by-sequencing derived single nucleotide polymorphisms. The physical locations of 56 single and 41 epistatic quantitative trait locus (QTL) were identified. Four head quality QTL (OQ_C03@57.0, OQ_C04@33.3, OQ_CC08@25.5, and OQ_C09@49.7) explain a cumulative 81.9% of phenotypic variance in the broccoli heading phenotype, contain the FLOWERING LOCUS C (FLC) homologs Bo9g173400 and Bo9g173370, and exhibit epistatic effects. Three key genomic hotspots associated with pleiotropic control of the broccoli heading phenotype were identified. One phenology hotspot reduces days to flowering by 7.0 days and includes an additional FLC homolog Bo3g024250 that does not exhibit epistatic effects with the three horticultural quality hotspots. Strong candidates for other horticultural traits were identified: BoLMI1 (Bo3g002560) associated with serrated leaf margins and leaf apex shape, BoCCD4 (Bo3g158650) implicated in flower color, and BoAP2 (Bo1g004960) implicated in the hooked sepal horticultural trait. The BolTBDH population provides a framework for B. oleracea improvement by targeting key genomic loci contributing to high horticultural quality broccoli and enabling de novo mapping of currently unexplored traits.

Bioactive Compounds in Brassicaceae Vegetables with a Role in the Prevention of Chronic Diseases

The beneficial role of the Mediterranean diet in the prevention of chronic diseases, including cardiovascular diseases, diabetes, and obesity, is well-recognized. In this context, Brassicaceae are considered important vegetables due to several evidences of their health promoting effects that are associated to bioactive compounds present in the edible parts of the plants. In this review, the mechanisms of action and the factors regulating the levels of the bioactive compounds in Brassicaceae have been discussed. In addition, the impact of industrial and domestic processing on the amount of these compounds have been considered, in order to identify the best conditions that are able to preserve the functional properties of the Brassicaceae products before consumption. Finally, the main strategies used to increase the content of health-promoting metabolites in Brassica plants through biofortification have been analyzed.

Quantitative trait loci mapping of heat tolerance in broccoli (Brassica oleracea var. italica) using genotyping-by-sequencing

Five quantitative trait loci and one epistatic interaction were associated with heat tolerance in a doubled haploid population of broccoli evaluated in three summer field trials. Predicted rising global temperatures due to climate change have generated a demand for crops that are resistant to yield and quality losses from heat stress. Broccoli (Brassica oleracea var. italica) is a cool weather crop with high temperatures during production decreasing both head quality and yield. Breeding for heat tolerance in broccoli has potential to both expand viable production areas and extend the growing season but breeding efficiency is constrained by limited genetic information. A doubled haploid (DH) broccoli population segregating for heat tolerance was evaluated for head quality in three summer fields in Charleston, SC, USA. Multiple quantitative trait loci (QTL) mapping of 1,423 single nucleotide polymorphisms developed through genotyping-by-sequencing identified five QTL and one positive epistatic interaction that explained 62.1% of variation in heat tolerance. The QTL identified here can be used to develop markers for marker-assisted selection and to increase our understanding of the molecular mechanisms underlying plant response to heat stress.

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

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

Multigenic Control of Pod Shattering Resistance in Chinese Rapeseed Germplasm Revealed by Genome-Wide Association and Linkage Analyses

The majority of rapeseed cultivars shatter seeds upon maturity especially under hot-dry and windy conditions, reducing yield and gross margin return to growers. Here, we identified quantitative trait loci (QTL) for resistance to pod shatter in an unstructured diverse panel of 143 rapeseed accessions, and two structured populations derived from bi-parental doubled haploid (DH) and inter-mated (IF2) crosses derived from R1 (resistant to pod shattering) and R2 (prone to pod shattering) accessions. Genome-wide association analysis identified six significant QTL for resistance to pod shatter located on chromosomes A01, A06, A07, A09, C02, and C05. Two of the QTL, qSRI.A09 delimited with the SNP marker Bn-A09-p30171993 (A09) and qSRI.A06 delimited with the SNP marker Bn-A06-p115948 (A06) could be repeatedly detected across environments in a diversity panel, DH and IF2 populations, suggesting that at least two loci on chromosomes A06 and A09 were the main contributors to pod shatter resistance in Chinese germplasm. Significant SNP markers identified in this study especially those that appeared repeatedly across environments provide a cost-effective and an efficient method for introgression and pyramiding of favorable alleles for pod shatter resistance via marker-assisted selection in rapeseed improvement programs.

QTL Analysis of Head Splitting Resistance in Cabbage (Brassica oleracea L. var. capitata) Using SSR and InDel Makers Based on Whole-Genome Re-Sequencing

Head splitting resistance (HSR) in cabbage is an important trait closely related to both quality and yield of head. However, the genetic control of this trait remains unclear. In this study, a doubled haploid (DH) population derived from an intra-cross between head splitting-susceptible inbred cabbage line 79-156 and resistant line 96-100 was obtained and used to analyze inheritance and detect quantitative trait loci (QTLs) for HSR using a mixed major gene/polygene inheritance analysis and QTL mapping. HSR can be attributed to additive-epistatic effects of three major gene pairs combined with those of polygenes. Negative and significant correlations were also detected between head Hsr and head vertical diameter (Hvd), head transverse diameter (Htd) and head weight (Hw). Using the DH population, a genetic map was constructed with simple sequence repeat (SSR) and insertion-deletion (InDel) markers, with a total length of 1065.9 cM and average interval length of 4.4 cM between adjacent markers. Nine QTLs for HSR were located on chromosomes C3, C4, C7, and C9 based on 2 years of phenotypic data using both multiple-QTL mapping and inclusive composite interval mapping. The identified QTLs collectively explained 39.4 to 59.1% of phenotypic variation. Three major QTLs (Hsr 3.2, 4.2, 9.2) showing a relatively larger effect were robustly detected in different years or with different mapping methods. The HSR trait was shown to have complex genetic mechanisms. Results from QTL mapping and classical genetic analysis were consistent. The QTLs obtained in this study should be useful for molecular marker-assisted selection in cabbage breeding and provide a foundation for further research on HSR genetic regulation.

Genetic analysis of glucosinolate variability in broccoli florets using genome-anchored single nucleotide polymorphisms

The identification of genetic factors influencing the accumulation of individual glucosinolates in broccoli florets provides novel insight into the regulation of glucosinolate levels in Brassica vegetables and will accelerate the development of vegetables with glucosinolate profiles tailored to promote human health. Quantitative trait loci analysis of glucosinolate (GSL) variability was conducted with a B. oleracea (broccoli) mapping population, saturated with single nucleotide polymorphism markers from a high-density array designed for rapeseed (Brassica napus). In 4 years of analysis, 14 QTLs were associated with the accumulation of aliphatic, indolic, or aromatic GSLs in floret tissue. The accumulation of 3-carbon aliphatic GSLs (2-propenyl and 3-methylsulfinylpropyl) was primarily associated with a single QTL on C05, but common regulation of 4-carbon aliphatic GSLs was not observed. A single locus on C09, associated with up to 40 % of the phenotypic variability of 2-hydroxy-3-butenyl GSL over multiple years, was not associated with the variability of precursor compounds. Similarly, QTLs on C02, C04, and C09 were associated with 4-methylsulfinylbutyl GSL concentration over multiple years but were not significantly associated with downstream compounds. Genome-specific SNP markers were used to identify candidate genes that co-localized to marker intervals and previously sequenced Brassica oleracea BAC clones containing known GSL genes (GSL-ALK, GSL-PRO, and GSL-ELONG) were aligned to the genomic sequence, providing support that at least three of our 14 QTLs likely correspond to previously identified GSL loci. The results demonstrate that previously identified loci do not fully explain GSL variation in broccoli. The identification of additional genetic factors influencing the accumulation of GSL in broccoli florets provides novel insight into the regulation of GSL levels in Brassicaceae and will accelerate development of vegetables with modified or enhanced GSL profiles.

Recent progress in the use of 'omics technologies in brassicaceous vegetables

Continuing advances in 'omics methodologies and instrumentation is enhancing the understanding of how plants cope with the dynamic nature of their growing environment. 'Omics platforms have been only recently extended to cover horticultural crop species. Many of the most widely cultivated vegetable crops belong to the genus Brassica: these include plants grown for their root (turnip, rutabaga/swede), their swollen stem base (kohlrabi), their leaves (cabbage, kale, pak choi) and their inflorescence (cauliflower, broccoli). Characterization at the genome, transcript, protein and metabolite levels has illustrated the complexity of the cellular response to a whole series of environmental stresses, including nutrient deficiency, pathogen attack, heavy metal toxicity, cold acclimation, and excessive and sub-optimal irradiation. This review covers recent applications of 'omics technologies to the brassicaceous vegetables, and discusses future scenarios in achieving improvements in crop end-use quality.