Quantitative trait loci controlling leaf appearance and curd initiation of cauliflower in relation to temperature

Genetic dissection of morphological variation between cauliflower and a rapid cycling Brassica oleracea line

To improve resolution to small genomic regions and sensitivity to small-effect loci in the identification of genetic factors conferring the enlarged inflorescence and other traits of cauliflower while also expediting further genetic dissection, 104 near-isogenic introgression lines (NIILs) covering 78.56% of the cauliflower genome, were selected from an advanced backcross population using cauliflower [Brassica oleracea var. botrytis L., mutant for Orange gene (ORG)] as the donor parent and a rapid cycling line (TO1434) as recurrent parent. Subsets of the advanced backcross population and NIILs were planted in the field for 8 seasons, finding 141 marker-trait associations for 15 leaf-, stem-, and flower-traits. Exemplifying the usefulness of these lines, we delineated the previously known flower color gene to a 4.5 MB interval on C3; a gene for small plant size to a 3.4 MB region on C8; and a gene for large plant size and flowering time to a 6.1 MB region on C9. This approach unmasked closely linked QTL alleles with opposing effects (on chr. 8) and revealed both alleles with expected phenotypic effects and effects opposite the parental phenotypes. Selected B. oleracea NIILs with short generation time add new value to widely used research and teaching materials.

Understanding population structure and detection of QTLs for curding-related traits in Indian cauliflower by genotyping by sequencing analysis

Curd initiation and development are complex traits and highly responsive for different temperature ranges in cauliflower. The present study was aimed to identify QTLs for eight traits associated with curding behaviour in diverse germplasm of Indian cauliflower. For this, 92 genotypes of cauliflower and 2 each of tropical broccoli and cabbage were genotyped through genotyping by sequencing (GBS). It generated ≈302 million reads (9.1226E + 10 bp) and identified 35,381 SNPs, maximum from chromosome 3 (4735) with a mean value of 3981.1 SNPs. Ts/Tv ratio was 1.74, suggesting transition bias. STRUCTURE analysis revealed delta value of K = 4 and four subpopulations and prominence of population admixture. In total, 121 significant SNPs were detected for eight traits, 38 for Delhi (North Indian plain) and 83 for Barapani (North-East India). Twelve QTLs were detected for traits associated with regulation of curd formation and development, five of which were for marketable curd length, curd width, days to 50% curd harvest and marketable curd weight from Delhi region and seven for curd length, curd width, days to 50% curd harvest, gross plant weight, leaf length, marketable/net curd weight and number of leaves per plant for Barapani area of North East India. The SNPs identified will be useful for development of markers for curding-related traits and their use in breeding varieties with wider curding plasticity.

Genome sequencing sheds light on the contribution of structural variants to Brassica oleracea diversification

Background

Brassica oleracea includes several morphologically diverse, economically important vegetable crops, such as the cauliflower and cabbage. However, genetic variants, especially large structural variants (SVs), that underlie the extreme morphological diversity of B. oleracea remain largely unexplored.

Results

Here we present high-quality chromosome-scale genome assemblies for two B. oleracea morphotypes, cauliflower and cabbage. Direct comparison of these two assemblies identifies ~ 120 K high-confidence SVs. Population analysis of 271 B. oleracea accessions using these SVs clearly separates different morphotypes, suggesting the association of SVs with B. oleracea intraspecific divergence. Genes affected by SVs selected between cauliflower and cabbage are enriched with functions related to response to stress and stimulus and meristem and flower development. Furthermore, genes affected by selected SVs and involved in the switch from vegetative to generative growth that defines curd initiation, inflorescence meristem proliferation for curd formation, maintenance and enlargement, are identified, providing insights into the regulatory network of curd development.

Conclusions

This study reveals the important roles of SVs in diversification of different morphotypes of B. oleracea, and the newly assembled genomes and the SVs provide rich resources for future research and breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01031-2.

From landrace to modern hybrid broccoli: the genomic and morphological domestication syndrome within a diverse B. oleracea collection

Worldwide, broccoli (Brassica oleracea var. italica) is among the most economically important, nutritionally rich, and widely-grown vegetable crops. To explore the genomic basis of the dramatic changes in broccoli morphology in the last century, we evaluated 109 broccoli or broccoli/cauliflower intermediates for 24 horticultural traits. Genotype-by-sequencing markers were used to determine four subpopulations within italica: Calabrese broccoli landraces and hybrids, sprouting broccoli, and violet cauliflower, and to evaluate between and within group relatedness and diversity. While overall horticultural quality and harvest index of improved hybrid broccoli germplasm has increased by year of cultivar release, this improvement has been accompanied by a considerable reduction in allelic diversity when compared to the larger pool of germplasm. Two landraces are the most likely founding source of modern broccoli hybrids, and within these modern hybrids, we identified 13 reduction-in-diversity genomic regions, 53 selective sweeps, and 30 (>1 Mbp) runs of homozygosity. Landrace accessions collected in southern Italy contained 4.8-fold greater unique alleles per accessions compared to modern hybrids and provide a valuable resource in subsequent improvement efforts. This work broadens the understanding of broccoli germplasm, informs conservation efforts, and enables breeding for complex quality traits and regionally adapted cultivars.

Heterosis and combining ability in cytoplasmic male sterile and doubled haploid based Brassica oleracea progenies and prediction of heterosis using microsatellites

In Brassica oleracea, heterosis is the most efficient tool providing impetus to hybrid vegetable industry. In this context, we presented the first report on identifying superior heterotic crosses for yield and commercial traits in cauliflower involving cytoplasmic male sterile (CMS) and doubled haploid (DH) lines as parents. We studied the suitability of genomic-SSRs and EST-SSRs based genetic distance (GD) and agronomic trait based phenotypic distance (PD) for predicting heterosis in F₁ hybrids using CMS and DH based parents. 120 F₁ hybrids derived from 20Ogura based CMS lines and 6 DH based testers were evaluated for 16 agronomic traits along with the 26 parental lines and 4 commercial standard checks. The genomic-SSRs and EST-SSRs based genetic structure analysis grouped the 26 parental lines into 4 distinct clusters. The CMS lines Ogu118-6A, Ogu33A, Ogu34-1A were good general combiner for developing early maturity hybrids. The SCA effects were significantly associated with heterosis suggesting non-additive gene effects for the heterotic response of hybrids. Less than unity value of σ²A/D coupled with σ²_gca/σ²_sca indicated the predominance of non-additive gene action in the expression of studied traits. The correlation analysis of genetic distance with heterosis for commercial traits suggested that microsatellites based genetic distance estimates can be helpful in heterosis prediction to some extent.

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.

Translating Flowering Time From Arabidopsis thaliana to Brassicaceae and Asteraceae Crop Species

Flowering and seed set are essential for plant species to survive, hence plants need to adapt to highly variable environments to flower in the most favorable conditions. Endogenous cues such as plant age and hormones coordinate with the environmental cues like temperature and day length to determine optimal time for the transition from vegetative to reproductive growth. In a breeding context, controlling flowering time would help to speed up the production of new hybrids and produce high yield throughout the year. The flowering time genetic network is extensively studied in the plant model species Arabidopsis thaliana, however this knowledge is still limited in most crops. This article reviews evidence of conservation and divergence of flowering time regulation in A. thaliana with its related crop species in the Brassicaceae and with more distant vegetable crops within the Asteraceae family. Despite the overall conservation of most flowering time pathways in these families, many genes controlling this trait remain elusive, and the function of most Arabidopsis homologs in these crops are yet to be determined. However, the knowledge gathered so far in both model and crop species can be already exploited in vegetable crop breeding for flowering time control.

Subtropical adaptation of a temperate plant (Brassica oleracea var. italica) utilizes non-vernalization-responsive QTLs

While many tropical plants have been adapted to temperate cultivation, few temperate plants have been adapted to the tropics. Originating in Western Europe, Brassica oleracea vernalization requires a period of low temperature and BoFLC2 regulates the transition to floral development. In B. oleracea germplasm selected in Taiwan, a non-vernalization pathway involving BoFLC3 rather than BoFLC2 regulates curd induction. In 112 subtropical breeding lines, specific haplotype combinations of BoFLC3 and PAN (involved in floral organ identity and a positional candidate for additional curd induction variation) adapt B. oleracea to high ambient temperature and short daylength. Duplicated genes permitted evolution of alternative pathways for control of flowering in temperate and tropical environments, a principle that might be utilized via natural or engineered approaches in other plants. New insight into regulation of Brassica flowering exemplifies translational agriculture, tapping knowledge of botanical models to improve food security under projected climate change scenarios.

Genome-Based Prediction of Time to Curd Induction in Cauliflower

The development of cauliflower (Brassica oleracea var. botrytis) is highly dependent on temperature due to vernalization requirements, which often causes delay and unevenness in maturity during months with warm temperatures. Integrating quantitative genetic analyses with phenology modeling was suggested to accelerate breeding strategies toward wide-adaptation cauliflower. The present study aims at establishing a genome-based model simulating the development of doubled haploid (DH) cauliflower lines to predict time to curd induction of DH lines not used for model parameterization and test hybrids derived from the bi-parental cross. Leaf appearance rate and the relation between temperature and thermal time to curd induction were examined in greenhouse trials on 180 DH lines at seven temperatures. Quantitative trait loci (QTL) analyses carried out on model parameters revealed ten QTL for leaf appearance rate (LAR), five for the slope and two for the intercept of linear temperature-response functions. Results of the QTL-based phenology model were compared to a genomic selection (GS) model. Model validation was carried out on data comprising four field trials with 72 independent DH lines, 160 hybrids derived from the parameterization set, and 34 hybrids derived from independent lines of the population. The QTL model resulted in a moderately accurate prediction of time to curd induction (R2 = 0.42-0.51) while the GS model generated slightly better results (R2 = 0.52-0.61). Predictions of time to curd induction of test hybrids from independent DH lines were less precise with R2 = 0.40 for the QTL and R2 = 0.48 for the GS model. Implementation of juvenile-to-adult phase transition is proposed for model improvement.

Genomic Prediction and Association Mapping of Curd-Related Traits in Gene Bank Accessions of Cauliflower

Genetic resources are an important source of genetic variation for plant breeding. Genome-wide association studies (GWAS) and genomic prediction greatly facilitate the analysis and utilization of useful genetic diversity for improving complex phenotypic traits in crop plants. We explored the potential of GWAS and genomic prediction for improving curd-related traits in cauliflower (Brassica oleracea var. botrytis) by combining 174 randomly selected cauliflower gene bank accessions from two different gene banks. The collection was genotyped with genotyping-by-sequencing (GBS) and phenotyped for six curd-related traits at two locations and three growing seasons. A GWAS analysis based on 120,693 single-nucleotide polymorphisms identified a total of 24 significant associations for curd-related traits. The potential for genomic prediction was assessed with a genomic best linear unbiased prediction model and BayesB. Prediction abilities ranged from 0.10 to 0.66 for different traits and did not differ between prediction methods. Imputation of missing genotypes only slightly improved prediction ability. Our results demonstrate that GWAS and genomic prediction in combination with GBS and phenotyping of highly heritable traits can be used to identify useful quantitative trait loci and genotypes among genetically diverse gene bank material for subsequent utilization as genetic resources in cauliflower breeding.