QTL-seq identifies BnaFT.A02 and BnaFLC.A02 as candidates for variation in vernalization requirement and response in winter oilseed rape (Brassica napus)

Transcriptional Reprogramming and Allelic Variation in Pleiotropic QTL Regulates Days to Flowering and Growth Habit in Pigeonpea

The present study investigated the linkage between days to flowering (DTF) and growth habit (GH) in pigeonpea using QTL mapping, QTL-seq, and GWAS approaches. The linkage map developed here is the largest to date, spanning 1825.56 cM with 7987 SNP markers. In total, eight and four QTLs were mapped for DTF and GH, respectively, harbouring 78 pigeonpea orthologs of Arabidopsis flowering time genes. Corroboratively, QTL-seq analysis identified a single linked QTL for both traits on chromosome 3, possessing 15 genes bearing genic variants. Together, these 91 genes were clustered primarily into autonomous, photoperiod, and epigenetic pathways. Further, we identified 39 associations for DTF and 111 associations for GH through GWAS in the QTL regions. Of these, nine associations were consistent and constituted nine haplotypes (five late, two early, one each for super-early and medium duration). The involvement of multiple genes explained the range of allelic effects and the presence of multiple LD blocks. Further, the linked QTL on chromosome 3 was fine-mapped to the 0.24-Mb region with an LOD score of 8.56, explaining 36.47% of the phenotypic variance. We identified a 10-bp deletion in the first exon of TFL1 gene of the ICPL 20338 variety, which may affect its interaction with the Apetala1 and Leafy genes, resulting in determinate GH and early flowering. Further, the genic marker developed for the deletion in the TFL1 gene could be utilized as a foreground marker in marker-assisted breeding programmes to develop early-flowering pigeonpea varieties.

Allelic Variation of BnaFTA2 and BnaFTC6 Is Associated With Flowering Time and Seasonal Crop Type in Rapeseed (Brassica napus L.)

Different ecological types of rapeseed (Brassica napus L.), including winter, spring, and semi-winter cultivars, exhibit varying flowering times and cannot be planted in the same cultivation areas. FLOWERING LOCUS T (FT) plays a key role in regulating flowering. In allotetraploid B. napus six copies of FT (BnaFT) have been reported. However, there is uncertainty about how the translated products of each paralog, as well as cis-allelic variations at each locus, contribute functionally to flowering time and define specific crop types. In this study, we confirm that BnaFT exhibit distinct expression patterns in different crop types of rapeseed. Using the CRISPR/Cas9 gene editing system, we provide functional evidence that the mutants between Bnaft paralogues affects the regulation of flowering time. Furthermore, we identify a new haplotype of BnaFT.A2 that is associated with early flowering time, although this appears necessary but not sufficient to confer a spring type phenotype. Three haplotypes of BnaFT.C6 were further identified and associated with both flowering time and crop types. We speculate that variations in both BnaFT.A2 and BnaFT.C6 may have undergone diversifying selection during the divergence of seasonal crop types in rapeseed.

A complex receptor locus confers responsiveness to necrosis and ethylene-inducing like peptides in Brassica napus

Brassica crops are susceptible to diseases which can be mitigated by breeding for resistance. MAMPs (microbe-associated molecular patterns) are conserved molecules of pathogens that elicit host defences known as pattern-triggered immunity (PTI). Necrosis and Ethylene-inducing peptide 1-like proteins (NLPs) are MAMPs found in a wide range of phytopathogens. We studied the response to BcNEP2, a representative NLP from Botrytis cinerea, and showed that it contributes to disease resistance in Brassica napus. To map regions conferring NLP response, we used the production of reactive oxygen species (ROS) induced during PTI across a population of diverse B. napus accessions for associative transcriptomics (AT), and bulk segregant analysis (BSA) on DNA pools created from a cross of NLP-responsive and non-responsive lines. In silico mapping with AT identified two peaks for NLP responsiveness on chromosomes A04 and C05 whereas the BSA identified one peak on A04. BSA delimited the region for NLP-responsiveness to 3 Mbp, containing ~245 genes on the Darmor-bzh reference genome and four co-segregating KASP markers were identified. The same pipeline with the ZS11 genome confirmed the highest-associated region on chromosome A04. Comparative BLAST analysis revealed unannotated clusters of receptor-like protein (RLP) homologues on ZS11 chromosome A04. However, no specific RLP homologue conferring NLP response could be identified. Our results also suggest that BR-SIGNALLING KINASE1 may be involved with modulating the NLP response. Overall, we demonstrate that responsiveness to NLP contributes to disease resistance in B. napus and define the associated genomic location. These results can have practical application in crop improvement.

Genome-wide association study and transcriptome analysis reveal natural variation of key genes regulation flowering time in rapeseed

Flowering time is a crucial determinant of both the yield and oil quality of rapeseed as well as a key indicator of plant maturity. We performed a genome-wide association study and transcriptome analysis to identify key genes/loci regulating flowering time in Brassica napus L. Forty-six haplotype regions harboring candidate genes were determined to be significantly associated with flowering time, and 28 of these haplotype regions overlapped with previously reported quantitative trait loci. A further investigation of these haplotype regions revealed nucleotide variations in the genes BnaFT-A02, BnaFRI-A10, and BnaFPA-A09 that correlated with phenotypic variations in flowering time. Furthermore, the co-expression network analysis indicated that BnaFT-A02 is directly linked to BnaFRI-A10 and BnaFPA-A09, in a subnetwork and also associated with 13 vernalization, 31 photoperiod, 33 autonomous pathway, and 10 gibberellin pathway genes, forming a potential network regulating flowering time in rapeseed. These results provide valuable haplotype markers for the breeding of early maturing rapeseed varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01479-4.

Identification candidate genes for salt resistance through quantitative trait loci-sequencing in Brassica napus L

Rapeseed (Brassica napus L.) is one of the most important oil crops worldwide. However, its yield is greatly limited by salt stress, one of the primary abiotic stresses. Identification of salt-tolerance genes and breeding salt-tolerant varieties is an effective approach to address this issue. Unfortunately, little is known about the salt-tolerance quantitative trait locus (QTL) and the identification of salt tolerance genes in rapeseed. In this study, high-throughput quantitative trait locus sequencing (QTL-seq) was applied to identifying salt-tolerant major QTLs based on two DNA pools from an F2:3 population of a cross between rapeseed line 2205 (salt tolerant) and 1423 (salt sensitive). A total of twelve major QTLs related to the salt tolerance rating (STR) were detected on chromosomes A03, A08, C02, C03, C04, C06, C07 and C09. To further enhance our understanding, we integrated QTL-seq data with transcriptome analysis of the two parental rapeseed plants subjected to salt stress, through which ten candidate genes for salt tolerance were identified within the major QTLs by gene differential expression, variation and annotated functions analysis. The marker SNP820 linked to salt tolerance was successfully validated and would be useful as a diagnostic marker in marker-assisted breeding. These findings provide valuable insights for future breeding programs aimed at developing rapeseed cultivars resistant to salt stresses.

A transcriptomic time-series reveals differing trajectories during pre-floral development in the apex and leaf in winter and spring varieties of Brassica napus

Oilseed rape (Brassica napus) is an important global oil crop, with spring and winter varieties grown commercially. To understand the transcriptomic differences between these varieties, we collected transcriptomes from apex and leaf tissue from a spring variety, Westar, and a winter variety, Tapidor, before, during, and after vernalisation treatment, until the plants flowered. Large transcriptomic differences were noted in both varieties during the vernalisation treatment because of temperature and day length changes. Transcriptomic alignment revealed that the apex transcriptome reflects developmental state, whereas the leaf transcriptome is more closely aligned to the age of the plant. Similar numbers of copies of genes were expressed in both varieties during the time series, although key flowering time genes exhibited expression pattern differences. BnaFLC copies on A2 and A10 are the best candidates for the increased vernalisation requirement of Tapidor. Other BnaFLC copies show tissue-dependent reactivation of expression post-cold, with these dynamics suggesting some copies have retained or acquired a perennial nature. BnaSOC1 genes, also related to the vernalisation pathway, have expression profiles which suggest tissue subfunctionalisation. This understanding may help to breed varieties with more consistent or robust vernalisation responses, of special importance due to the milder winters resulting from climate change.

Fine mapping and analysis of candidate genes for qBT2 and qBT7.2 locus controlling bolting time in radish (Raphanus sativus L.)

KEY MESSAGE

Two major QTLs for bolting time in radish were mapped to chromosome 02 and 07 in a 0.37 Mb and 0. 52 Mb interval, RsFLC1 and RsFLC2 is the critical genes. Radish (Raphanus sativus L.) is an important vegetable crop of Cruciferae. The premature bolting and flowering reduces the yield and quality of the fleshy root of radish. However, the molecular mechanism underlying bolting and flowering in radish remains unknown. In YZH (early bolting) × XHT (late bolting) F2 population, a high-density genetic linkage map was constructed with genetic distance of 2497.74 cM and an average interval of 2.31 cM. A total of nine QTLs for bolting time and two QTLs for flowering time were detected. Three QTLs associated with bolting time in radish were identified by QTL-seq using radish GDE (early bolting) × GDL (late bolting) F2 population. Fine mapping narrowed down qBT2 and qBT7.2 to an 0.37 Mb and 0.52 Mb region on chromosome 02 and 07, respectively. RNA-seq and qRT-PCR analysis showed that RsFLC1 and RsFLC2 were the candidate gene for qBT7.2 and qBT2 locus, respectively. Subcellular localization exhibited that RsFLC1 and RsFLC2 were mainly expressed in the nucleus. A 1856-bp insertion in the first intron of RsFLC1 was responsible for bolting time. Overexpression of RsFLC2 in Arabidopsis was significantly delayed flowering. These findings will provide new insights into the exploring the molecular mechanism of late bolting and promote the marker-assisted selection for breeding late-bolting varieties in radish.

High ambient temperature impacts on flowering time in Brassica napus through both H2A.Z-dependent and independent mechanisms

Knowledge concerning the integration of genetic pathways mediating the responses to environmental cues controlling flowering initiation in crops is scarce. Here, we reveal the diversity in oilseed rape (OSR) flowering response to high ambient temperature. Using a set of different spring OSR varieties, we found a consistent flowering delay at elevated temperatures. Remarkably, one of the varieties assayed exhibited the opposite behaviour. Several FT-like paralogs are plausible candidates to be part of the florigen in OSR. We revealed that BnaFTA2 plays a major role in temperature-dependent flowering initiation. Analysis of the H2A.Z histone variant occupancy at this locus in different Brassica napus varieties produced contrasting results, suggesting the involvement of additional molecular mechanisms in BnaFTA2 repression at high ambient temperature. Moreover, BnARP6 RNAi plants showed little accumulation of H2A.Z at high temperature while maintaining temperature sensitivity and delayed flowering. Furthermore, we found that H3K4me3 present in BnaFTA2 under inductive flowering conditions is reduced at high temperature, suggesting a role for this hallmark of transcriptionally active chromatin in the OSR flowering response to warming. Our work emphasises the plasticity of flowering responses in B. napus and offers venues to optimise this process in crop species grown under suboptimal environmental conditions.

Rapid identification of a major locus qPRL-C06 affecting primary root length in Brassica napus by QTL-seq

BACKGROUND AND AIMS

Brassica napus is one of the most important oilseed crops worldwide. Seed yield of B. napus significantly correlates with the primary root length (PRL). The aims of this study were to identify quantitative trait loci (QTLs) for PRL in B. napus.

METHODS

QTL-seq and conventional QTL mapping were jointly used to detect QTLs associated with PRL in a B. napus double haploid (DH) population derived from a cross between 'Tapidor' and 'Ningyou 7'. The identified major locus was confirmed and resolved by an association panel of B. napus and an advanced backcross population. RNA-seq analysis of two long-PRL lines (Tapidor and TN20) and two short-PRL lines (Ningyou 7 and TN77) was performed to identify differentially expressed genes in the primary root underlying the target QTLs.

KEY RESULTS

A total of 20 QTLs impacting PRL in B. napus grown at a low phosphorus (P) supply were found by QTL-seq. Eight out of ten QTLs affecting PRL at a low P supply discovered by conventional QTL mapping could be detected by QTL-seq. The locus qPRL-C06 identified by QTL-seq was repeatedly detected at both an optimal P supply and a low P supply by conventional QTL mapping. This major constitutive QTL was further confirmed by regional association mapping. qPRL-C06 was delimited to a 0.77 Mb genomic region on chromosome C06 using an advanced backcross population. A total of 36 candidate genes within qPRL-C06 were identified that showed variations in coding sequences and/or exhibited significant differences in mRNA abundances in primary root between the long-PRL and short-PRL lines, including five genes involved in phytohormone biosynthesis and signaling.

CONCLUSIONS

These results both demonstrate the power of the QTL-seq in rapid QTL detection for root traits and will contribute to marker-assisted selective breeding of B. napus cultivars with increased PRL.

Fine mapping and causal gene identification of a novel QTL for early flowering by QTL-seq, Target-seq and RNA-seq in spring oilseed rape

A novel quantitative trait locus for early flowering in spring oilseed rape, BnaC08cqDTF, was mapped to an 86-kb region on chromosome C08, and its causal gene, CRY2, was uncovered. Days to flowering is a very important agronomic and adaptive trait of Brassica napus oilseed rape (AACC, 2n = 38). We previously identified BnaC08cqDTF as a novel candidate quantitative trait locus (QTL) for early flowering in spring oilseed rape. Here, we present fine mapping of the locus and a study of its causal gene. Initial mapping was performed by QTL sequencing of DNA pools of BC₃F₂ plants with extreme flowering times derived from crosses between the spring-type cv. No. 4512 (early flowering) and cv. No. 5246 (late flowering), along with fine mapping by target sequencing of the BC₃F₂ and BC₄F₂ populations. Fine mapping narrowed down BnaC08cqDTF to an 86-kb region on chromosome C08. The region harbored fifteen genes. After comparative analyses of the DNA sequences for mutation between A and C syntenic regions and detected by RNA-seq and qRT-PCR between the two parents, we found that BnaC08G0010400ZS harbors an A/G nonsynonymous mutation in exon 3. This single nucleotide polymorphism (SNP) haplotype was also correlated with early flowering in a 256 accession panel. BnaC08G0010400ZS is a homolog of the AT1G04400 gene (CRY2) in Arabidopsis. The analyses of transgenic Arabidopsis verified that BnaC08G0010400ZS is responsible for early flowering. Our results contribute to a better understanding of the genetic control mechanism of early flowering in spring Brassica napus and will promote the breeding for early mature varieties.

QTL mapping for seed density per silique in Brassica napus

Seed density per silique (SDPS) and valid silique length (VSL) are two important yield-influencing traits in rapeseed. SDPS has a direct or indirect effect on rapeseed yield through its effect on seed per silique. In this study, a quantitative trait locus (QTL) for SDPS was detected on chromosome A09 using the QTL-seq approach and confirmed via linkage analysis in the mapping population obtained from 4263 × 3001 cross. Furthermore, one major QTL for SDPS (qSD.A9-1) was mapped to a 401.8 kb genomic interval between SSR markers Nys9A190 and Nys9A531. In the same genomic region, a QTL (qSL.A9) linked to VSL was also detected. The phenotypic variation of qSD.A9-1 and qSL.A9 was 53.1% and 47.6%, respectively. Results of the additive and dominant effects demonstrated that the expression of genes controlling SDPS and VSL were derived from a different parent in this population. Subsequently, we identified 56 genes that included 45 specific genes with exonic (splicing) variants. Further analysis identified specific genes containing mutations that may be related to seed density as well as silique length. These genes could be used for further studies to understand the details of these traits of rapeseed.

Whole genome resequencing identifies candidate genes and allelic diagnostic markers for resistance to Ralstonia solanacearum infection in cultivated peanut (Arachis hypogaea L.)

Bacterial wilt disease (BWD), caused by Ralstonia solanacearum is a major challenge for peanut production in China and significantly affects global peanut field productivity. It is imperative to identify genetic loci and putative genes controlling resistance to R. solanacearum (RRS). Therefore, a sequencing-based trait mapping approach termed "QTL-seq" was applied to a recombination inbred line population of 581 individuals from the cross of Yueyou 92 (resistant) and Xinhuixiaoli (susceptible). A total of 381,642 homozygous single nucleotide polymorphisms (SNPs) and 98,918 InDels were identified through whole genome resequencing of resistant and susceptible parents for RRS. Using QTL-seq analysis, a candidate genomic region comprising of 7.2 Mb (1.8-9.0 Mb) was identified on chromosome 12 which was found to be significantly associated with RRS based on combined Euclidean Distance (ED) and SNP-index methods. This candidate genomic region had 180 nonsynonymous SNPs and 14 InDels that affected 75 and 11 putative candidate genes, respectively. Finally, eight nucleotide binding site leucine rich repeat (NBS-LRR) putative resistant genes were identified as the important candidate genes with high confidence. Two diagnostic SNP markers were validated and revealed high phenotypic variation in the different resistant and susceptible RIL lines. These findings advocate the expediency of the QTL-seq approach for precise and rapid identification of candidate genomic regions, and the development of diagnostic markers that are applicable in breeding disease-resistant peanut varieties.

Necrosis and ethylene-inducing-like peptide patterns from crop pathogens induce differential responses within seven brassicaceous species

Translational research is required to advance fundamental knowledge on plant immunity towards application in crop improvement. Recognition of microbe/pathogen-associated molecular patterns (MAMPs/PAMPs) triggers a first layer of immunity in plants. The broadly occurring family of necrosis- and ethylene-inducing peptide 1 (NEP1)-like proteins (NLPs) contains immunogenic peptide patterns that are recognized by a number of plant species. Arabidopsis can recognize NLPs by the pattern recognition receptor AtRLP23 and its co-receptors SOBIR1, BAK1, and BKK1, leading to induction of defence responses including the production of reactive oxygen species (ROS) and elevation of intracellular [Ca2+]. However, little is known about NLP perception in Brassica crop species. Within 12 diverse accessions for each of six Brassica crop species, we demonstrate variation in response to Botrytis cinerea NLP BcNEP2, with Brassica oleracea (CC genome) being nonresponsive and only two Brassica napus cultivars responding to BcNEP2. Peptides derived from four fungal pathogens of these crop species elicited responses similar to BcNEP2 in B. napus and Arabidopsis. Induction of ROS by NLP peptides was strongly reduced in Atrlp23, Atsobir1 and Atbak1-5 Atbkk1-1 mutants, confirming that recognition of Brassica pathogen NLPs occurs in a similar manner to that of HaNLP3 from Hyaloperonospora arabidopsidis in Arabidopsis. In silico analysis of the genomes of two B. napus accessions showed similar presence of homologues for AtBAK1, AtBKK1 and AtSOBIR1 but variation in the organization of AtRLP23 homologues. We could not detect a strong correlation between the ability to respond to NLP peptides and resistance to B. cinerea.

Systematic trait dissection in oilseed rape provides a comprehensive view, further insight, and exact roadmap for yield determination

Background

Yield is the most important and complex trait that is influenced by numerous relevant traits with very complicated interrelations. While there are a large number of studies on the phenotypic relationship and genetic basis of yield traits, systematic studies with further dissection focusing on yield are limited. Therefore, there is still lack of a comprehensive and in-depth understanding of the determination of yield.

Results

In this study, yield was systematically dissected at the phenotypic, genetic to molecular levels in oilseed rape (Brassica napus L.). The analysis of correlation, network, and principal component for 21 traits in BnaZN-RIL population showed that yield was determined by a complex trait network with key contributors. The analysis of the constructed high-density single nucleotide polymorphism (SNP) linkage map revealed the concentrated distribution of distorted and heterozygous markers, likely due to selection on genes controlling the growth period and yield heterosis. A total of 134 consensus quantitative trait loci (QTL) were identified for 21 traits, of which all were incorporated into an interconnecting QTL network with dozens of hub-QTL. Four representative hub-QTL were further dissected to the target or candidate genes that governed the causal relationships between the relevant traits.

Conclusions

The highly consistent results at the phenotypic, genetic, and molecular dissecting demonstrated that yield was determined by a multilayer composite network that involved numerous traits and genes showing complex up/down-stream and positive/negative regulation. This provides a systematic view, further insight, and exact roadmap for yield determination, which represents a significant advance toward the understanding and dissection of complex traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02134-w.

Mapping of Two Major QTLs Controlling Flowering Time in Brassica napus Using a High-Density Genetic Map

Research on the flowering habit of rapeseed is important for the selection of varieties adapted to specific ecological environments. Here, quantitative trait loci (QTL) for the days-to-flowering trait were identified using a doubled haploid population of 178 lines derived from a cross between the winter type SGDH284 and the semi-winter type 158A. A linkage map encompassing 3268.01 cM was constructed using 2777 bin markers obtained from next-generation sequencing. The preliminary mapping results revealed 56 QTLs for the days to flowering in the six replicates in the three environments. Twelve consensus QTLs were identified by a QTL meta-analysis, two of which (cqDTF-C02 and cqDTF-C06) were designated as major QTLs. Based on the micro-collinearity of the target regions between B. napus and Arabidopsis, four genes possibly related to flowering time were identified in the cqDTF-C02 interval, and only one gene possibly related to flowering time was identified in the cqDTF-C06 interval. A tightly linked insertion-deletion marker for the cqFT-C02 locus was developed. These findings will aid the breeding of early maturing B. napus varieties.

Transposon insertions within alleles of BnaFT.A2 are associated with seasonal crop type in rapeseed

We identified two new transposon insertions within the promoter of BnaFT.A2 in addition to an existing 288 bp MITE within the second intron. Each insertion event corresponds to a distinct BnaFT.A2 haplotype and is closely associated with established crop seasonal ecotypes. Florigen, encoded by FLOWERING LOCUS T (FT), plays key roles not only as a flowering hormone, but also a universal growth factor affecting several aspects of plant architecture. In rapeseed, BnaFT.A2 has been revealed as one of the major loci associated with flowering time and different ecotypes. However, it is unclear how allelic variations of BnaFT.A2 affect its function in flowering time regulation and beyond. In this study, we confirmed an existing 288 bp miniature inverted-repeat transposable element (MITE) insertion within the second intron and identified two new insertions within the promoter of BnaFT.A2-a 3971 bp CACTA and a 1079 bp Helitron. Each insertion event corresponds to a distinct BnaFT.A2 haplotype and is closely associated with established crop seasonal ecotypes. These alleles have similar tissue-specific expression patterns but discrete transcriptional patterns tightly associated with rapeseed flowering time and ecotype. RNAi lines and mutants of BnaFT.A2 flowered significantly later than controls. Differentially expressed genes (DEGs), identified in transcriptomic profiling of seedling leaves from two loss-of-function mutants (Bnaft.a2-L1 and Bnaft.a2-L2) compared with controls, indicated significant enrichment for hormone metabolic genes and roles related to plant cell wall synthesis and photosynthesis. Plants with loss-of-function BnaFT.A2 had smaller leaves and lower net photosynthetic rate compared to controls. These findings not only further clarify the genetic basis of flowering time variation and ecotype formation in B. napus, but also provide an additional toolbox for genetic improvement of seasonal adaptation and production.

BnGF14-2c Positively Regulates Flowering via the Vernalization Pathway in Semi-Winter Rapeseed

14-3-3s are general regulatory factors (GF14s or GRFs) involved in a variety of physiological regulations in plants, including the control of flowering time. However, there are poorly relevant reports in rapeseed so far. In this study, we identified a homologous 14-3-3 gene BnGF14-2c (AtGRF2_Like in Brassica napus) in rapeseed based on bioinformatic analysis by using the sequences of the flowering-related 14-3-3s in other plant species. Then, we found that overexpression of BnGF14-2c in the semi-winter rapeseed “93275” promoted flowering without vernalization. Moreover, both yeast two-hybrid and bimolecular fluorescence complementation analysis indicated that BnGF14-2c may interact with two vernalization-related flowering regulators BnFT.A02 and BnFLC.A10., respectively. qPCR analysis showed that the expression of BnFT (AtFT_Like) was increased and the expression of two selected vernalization-related genes were reduced in the overexpression transgenic plants. Further investigation on subcellular localization demonstrated that BnGF14-2c localized in the nucleus and cytoplasm. The results of RNA-seq analysis and GUS staining indicated that BnGF14-2c is ubiquitously expressed except for mature seed coat. In general, the interaction of 14-3-3 and FLC was firstly documented in this study, indicating BnGF14-2c may act as a positive regulator of flowering in rapeseed, which is worthy for more in-depth exploration.

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

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

QTL mapping for low temperature germination in rapeseed

Rapeseed, a major oil crop in the world, is easily affected by low-temperature stress. A low temperature delays seed germination and increases seedling mortality, adversely affecting rapeseed growth and production. In the present study, a tolerant cultivar (Huyou21) was crossed with a susceptible genotype (3429) to develop a mapping population consisting of 574 F₂ progenies and elucidate the genetic mechanisms of seed germination under low temperatures. Two quantitative trait loci (QTL) for low-temperature germination (LTG) were detected, one on chromosome A09 (named qLTGA9-1) and the other on chromosome C01 (named qLTGC1-1), using the QTL-seq approach and confirmed via linkage analysis in the mapping population. Further, qLTGA9-1 was mapped to a 341.86 kb interval between the SSR markers Nys9A212 and Nys9A215. In this region, 69 genes including six specific genes with moderate or high effect function variants were identified based on the Ningyou7 genome sequence. Meanwhile, qLTGC1-1 was mapped onto a 1.31 Mb interval between SSR markers Nys1C96 and Nys1C117. In this region, 133 genes including five specific genes with moderate effect function variants were identified. These specific genes within the two QTL could be used for further studies on cold tolerance and as targets in rapeseed breeding programs.

Integrative analysis of GWAS and transcriptome to reveal novel loci regulation flowering time in semi-winter rapeseed

Flowering is an important turning point from vegetative growth to reproductive growth, and vernalization is an essential condition for the flowering of annual winter plants. To investigate the genetic architecture of flowering time in rapeseed, we used the 60 K Brassica Infinium SNP array to perform a genome-wide analysis of haplotype blocks associated with flowering time in 203 Chinese semi-winter rapeseed inbred lines. Twenty-one haplotype regions carrying one or more candidate genes showed a significant association with flowering time. Interestingly, we detected a SNP (Bn-scaff_22728_1-p285715) located in exon 3 of the BnVIN3-C03 gene that showed a significant association with flowering time on chromosome C03. Based on the SNP alleles A and G, two groups of accessions with early and late flowering time phenotypes were selected, respectively, and PCR amplification and gene expression analysis were combined to reveal the structural variation of the BnVIN3-C03 gene that affected flowering time. Moreover, we found that BnVIN3-C03 inhibited the expression of BnFLC-A02, BnFLC-A03.1, BnFLC-A10 and BnFLC-C03.1, thus modulating the flowering time of Brassica napus. This result provides insight into the genetic improvement of flowering time in B. napus.

Understanding Omics Driven Plant Improvement and de novo Crop Domestication: Some Examples

In the current era, one of biggest challenges is to shorten the breeding cycle for rapid generation of a new crop variety having high yield capacity, disease resistance, high nutrient content, etc. Advances in the "-omics" technology have revolutionized the discovery of genes and bio-molecules with remarkable precision, resulting in significant development of plant-focused metabolic databases and resources. Metabolomics has been widely used in several model plants and crop species to examine metabolic drift and changes in metabolic composition during various developmental stages and in response to stimuli. Over the last few decades, these efforts have resulted in a significantly improved understanding of the metabolic pathways of plants through identification of several unknown intermediates. This has assisted in developing several new metabolically engineered important crops with desirable agronomic traits, and has facilitated the de novo domestication of new crops for sustainable agriculture and food security. In this review, we discuss how "omics" technologies, particularly metabolomics, has enhanced our understanding of important traits and allowed speedy domestication of novel crop plants.

A novel deletion in FLOWERING LOCUS T modulates flowering time in winter oilseed rape

A novel structural variant was discovered in the FLOWERING LOCUS T orthologue BnaFT.A02 by long-read sequencing. Nested association mapping in an elite winter oilseed rape population revealed that this 288 bp deletion associates with early flowering, putatively by modification of binding-sites for important flowering regulation genes. Perfect timing of flowering is crucial for optimal pollination and high seed yield. Extensive previous studies of flowering behavior in Brassica napus (canola, rapeseed) identified mutations in key flowering regulators which differentiate winter, semi-winter and spring ecotypes. However, because these are generally fixed in locally adapted genotypes, they have only limited relevance for fine adjustment of flowering time in elite cultivar gene pools. In crosses between ecotypes, the ecotype-specific major-effect mutations mask minor-effect loci of interest for breeding. Here, we investigated flowering time in a multiparental mapping population derived from seven elite winter oilseed rape cultivars which are fixed for major-effect mutations separating winter-type rapeseed from other ecotypes. Association mapping revealed eight genomic regions on chromosomes A02, C02 and C03 associating with fine modulation of flowering time. Long-read genomic resequencing of the seven parental lines identified seven structural variants coinciding with candidate genes for flowering time within chromosome regions associated with flowering time. Segregation patterns for these variants in the elite multiparental population and a diversity set of winter types using locus-specific assays revealed significant associations with flowering time for three deletions on chromosome A02. One of these was a previously undescribed 288 bp deletion within the second intron of FLOWERING LOCUS T on chromosome A02, emphasizing the advantage of long-read sequencing for detection of structural variants in this size range. Detailed analysis revealed the impact of this specific deletion on flowering-time modulation under extreme environments and varying day lengths in elite, winter-type oilseed rape.

A Vernalization Response in a Winter Safflower (Carthamus tinctorius) Involves the Upregulation of Homologs of FT, FUL, and MAF

Safflower (Carthamus tinctorius) is a member of the Asteraceae family that is grown in temperate climates as an oil seed crop. Most commercially grown safflower varieties can be sown in late winter or early spring and flower rapidly in the absence of overwintering. There are winter-hardy safflower accessions that can be sown in autumn and survive over-wintering. Here, we show that a winter-hardy safflower possesses a vernalization response, whereby flowering is accelerated by exposing germinating seeds to prolonged cold. The impact of vernalization was quantitative, such that increasing the duration of cold treatment accelerated flowering to a greater extent, until the response was saturated after 2 weeks exposure to low-temperatures. To investigate the molecular-basis of the vernalization-response in safflower, transcriptome activity was compared and contrasted between vernalized versus non-vernalized plants, in both 'winter hardy' and 'spring' cultivars. These genome-wide expression analyses identified a small set of transcripts that are both differentially expressed following vernalization and that also have different expression levels in the spring versus winter safflowers. Four of these transcripts were quantitatively induced by vernalization in a winter hardy safflower but show high basal levels in spring safflower. Phylogenetic analyses confidently assigned that the nucleotide sequences of the four differentially expressed transcripts are related to FLOWERING LOCUS T (FT), FRUITFUL (FUL), and two genes within the MADS-like clade genes. Gene models were built for each of these sequences by assembling an improved safflower reference genome using PacBio-based long-read sequencing, covering 85% of the genome, with N50 at 594,000 bp in 3000 contigs. Possible evolutionary relationships between the vernalization response of safflower and those of other plants are discussed.

Total FLC transcript dynamics from divergent paralogue expression explains flowering diversity in Brassica napus

Flowering time is a key adaptive and agronomic trait. In Arabidopsis, natural variation in expression levels of the floral repressor FLOWERING LOCUS C (FLC) leads to differences in vernalization. In Brassica napus there are nine copies of FLC. Here, we study how these multiple FLC paralogues determine vernalization requirement as a system. We collected transcriptome time series for Brassica napus spring, winter, semi-winter, and Siberian kale crop types. Modelling was used to link FLC expression dynamics to floral response following vernalization. We show that relaxed selection pressure has allowed expression of FLC paralogues to diverge, resulting in variation of FLC expression during cold treatment between paralogues and accessions. We find that total FLC expression dynamics best explains differences in cold requirement between cultivars, rather than expression of specific FLC paralogues. The combination of multiple FLC paralogues with different expression dynamics leads to rich behaviour in response to cold and a wide range of vernalization requirements in B. napus. We find evidence for different strategies to determine the response to cold in existing winter rapeseed accessions.

Regulation and Subfunctionalization of Flowering Time Genes in the Allotetraploid Oil Crop Brassica napus

Flowering is a vulnerable, but crucial phase in building crop yield. Proper timing of this period is therefore decisive in obtaining optimal yields. However, genetic regulation of flowering integrates many different environmental signals and is therefore extremely complex. This complexity increases in polyploid crops which carry two or more chromosome sets, like wheat, potato or rapeseed. Here, I summarize the current state of knowledge about flowering time gene copies in rapeseed (Brassica napus), an important oil crop with a complex polyploid history and a close relationship to Arabidopsis thaliana. The current data show a high demand for more targeted studies on flowering time genes in crops rather than in models, allowing better breeding designs and a deeper understanding of evolutionary principles. Over evolutionary time, some copies of rapeseed flowering time genes changed or lost their original role, resulting in subfunctionalization of the respective homologs. For useful applications in breeding, such patterns of subfunctionalization need to be identified and better understood.

Genome Triplication Leads to Transcriptional Divergence of FLOWERING LOCUS C Genes During Vernalization in the Genus Brassica

The genus Brassica includes oil crops, vegetables, condiments, fodder crops, and ornamental plants. Brassica species underwent a whole genome triplication event after speciation between ancestral species of Brassica and closely related genera including Arabidopsis thaliana. Diploid species such as Brassica rapa and Brassica oleracea have three copies of genes orthologous to each A. thaliana gene, although deletion in one or two of the three homologs has occurred in some genes. The floral transition is one of the crucial events in a plant's life history, and time of flowering is an important agricultural trait. There is a variation in flowering time within species of the genus Brassica, and this variation is largely dependent on a difference in vernalization requirements. In Brassica, like in A. thaliana, the key gene of vernalization is FLOWERING LOCUS C (FLC). In Brassica species, the vernalization response including the repression of FLC expression by cold treatment and the enrichment of the repressive histone modification tri-methylated histone H3 lysine 27 (H3K27me3) at the FLC locus is similar to A. thaliana. B. rapa and B. oleracea each have four paralogs of FLC, and the allotetraploid species, Brassica napus, has nine paralogs. The increased number of paralogs makes the role of FLC in vernalization more complicated; in a single plant, paralogs vary in the expression level of FLC before and after vernalization. There is also variation in FLC expression levels between accessions. In this review, we focus on the regulatory circuits of the vernalization response of FLC expression in the genus Brassica.