Diverse regulatory factors associate with flowering time and yield responses in winter-type Brassica napus

Genetic control of root morphology in response to nitrogen across rapeseed diversity

Rapeseed (Brassica napus L.) is an oil-containing crop of great economic value but with considerable nitrogen requirement. Breeding root systems that efficiently absorb nitrogen from the soil could be a driver to ensure genetic gains for more sustainable rapeseed production. The aim of this study is to identify genomic regions that regulate root morphology in response to nitrate availability. The natural variability offered by 300 inbred lines was screened at two experimental locations. Seedlings grew hydroponically with low or elevated nitrate levels. Fifteen traits related to biomass production and root morphology were measured. On average across the panel, a low nitrate level increased the root-to-shoot biomass ratio and the lateral root length. A large phenotypic variation was observed, along with important heritability values and genotypic effects, but low genotype-by-nitrogen interactions. Genome-wide association study and bulk segregant analysis were used to identify loci regulating phenotypic traits. The first approach nominated 319 SNPs that were combined into 80 QTLs. Three QTLs identified on the A07 and C07 chromosomes were stable across nitrate levels and/or experimental locations. The second approach involved genotyping two groups of individuals from an experimental F2 population created by crossing two accessions with contrasting lateral root lengths. These individuals were found in the tails of the phenotypic distribution. Co-localized QTLs found in both mapping approaches covered a chromosomal region on the A06 chromosome. The QTL regions contained some genes putatively involved in root organogenesis and represent selection targets for redesigning the root morphology of rapeseed.

Investigation of mutation load and rate in androgenic mutant lines of rapeseed in early generations evaluated by high-density SNP genotyping

Oilseed rape (Brassica napus) is an important oil crop distributed worldwide with a broad adaptation to different climate zones. The cultivation of rapeseed is one of the most commercially viable areas in crop production. Altogether 269,093 ha of rapeseed are cultivated in Kazakhstan. However, all rapeseed cultivars and lines cultivated in Kazakhstan on an industrial scale predominantly belong to the foreign breeding system. Therefore, the formation of a diverse genetic pool for breeding new, highly productive cultivars adopted to the environmental conditions of Kazakhstan is the most important goal in country selection programs. In this work, we have developed ethyl methanesulfonate (EMS) doubled haploid mutant lines from plant material of cultivars 'Galant' and 'Kris' to broad diversity of rapeseed in Kazakhstan. The development of mutant lines was performed via embryo callusogenesis or embryo secondary callusogenesis. Mutants were investigated by Brassica90k SNP array, and we were able to locate 24,657 SNPs from 26,256 SNPs filtered by quality control on the genome assembly (Bra_napus_v2.0). Only 18,831 SNPs were assigned to the available annotated genomic features. The most frequent combination of mutations according to reference controls was adenine with guanine (70%), followed by adenine with cytosine (28.8%), and only minor fractions were cytosine with guanine (0.54%) and adenine with thymine (0.59%). We revealed 5606.27 markers for 'Kris' and 4893.01 markers for 'Galant' by mutation occurrence. Most mutation occurrences were occupied by double mutations where progenitors and offspring were homozygous by different alleles, enabling the selection of appropriate genotypes in a short period of time. Regarding the biological impact of mutations, 861 variants were reported as having a low predicted impact, with 1042 as moderate and 121 as high; all others were reported as belonging to non-coding sequences, intergenic regions, and other features with the effect of modifiers. Protein encoding genes, such as wall-associated receptor kinase-like protein 5, TAO1-like disease resistance protein, receptor-like protein 12, and At5g42460-like F-box protein, contained more than two variable positions, with an impact on their biological activities. Nevertheless, the obtained mutant lines were able to survive and reproduce. Mutant lines, which include moderate and high impact mutations in encoding genes, are a perfect pool not only for MAS but also for the investigation of the fundamental basis of protein functions. For the first time, a collection of mutant lines was developed in our country to improve the selection of local rapeseed cultivars.

Genome-wide association study of 23 flowering phenology traits and 4 floral agronomic traits in tree peony (Paeonia section Moutan DC.) reveals five genes known to regulate flowering time

Tree peony is a unique traditional flower in China, with large, fragrant, and colorful flowers. However, a relatively short and concentrated flowering period limits the applications and production of tree peony. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of flowering phenology traits and ornamental phenotypes in tree peony. A diverse panel of 451 tree peony accessions was phenotyped for 23 flowering phenology traits and 4 floral agronomic traits over 3 years. Genotyping by sequencing (GBS) was used to obtain a large number of genome-wide single-nucleotide polymorphisms (SNPs) (107 050) for the panel genotypes, and 1047 candidate genes were identified by association mapping. Eighty-two related genes were observed during at least 2 years for flowering, and seven SNPs repeatedly identified for multiple flowering phenology traits over multiple years were highly significantly associated with five genes known to regulate flowering time. We validated the temporal expression profiles of these candidate genes and highlighted their possible roles in the regulation of flower bud differentiation and flowering time in tree peony. This study shows that GWAS based on GBS can be used to identify the genetic determinants of complex traits in tree peony. The results expand our understanding of flowering time control in perennial woody plants. Identification of markers closely related to these flowering phenology traits can be used in tree peony breeding programs for important agronomic traits.

Genome-wide association study identifies candidate genes and favorable haplotypes for seed yield in Brassica napus

Oilseed rape (Brassica napus L.) is one of the most essential oil crops. Genetic improvement of seed yield (SY) is a major aim of B. napus breeding. Several studies have been reported on the genetic mechanisms of SY of B. napus. Here, a genome-wide association study (GWAS) of SY was conducted using a panel of 403 natural accessions of B. napus, with more than five million high-quality single-nucleotide polymorphisms (SNPs). A total of 1773 significant SNPs were detected associated with SY, and 783 significant SNPs were co-located with previously reported QTLs. The lead SNPs chrA01__8920351 and chrA02__4555979 were jointly detected in Trial 2_2 and Trial 2_mean value, and in Trial 1_2 and Trial 1_mean value, respectively. Subsequently, two candidate genes of BnaA01g17200D and BnaA02g08680D were identified through combining transcriptome, candidate gene association analysis, and haplotype analysis. BnaA09g10430D detected through lead SNP chrA09__5160639 was associated with SY of B. napus. Our results provide valuable information for studying the genetic control of seed yield in B. napus and valuable genes, haplotypes, and cultivars resources for the breeding of high seed yield B. napus cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01332-6.

A genome-wide association study (GWAS) identifies multiple loci linked with the natural variation for Al3+ resistance in Brassica napus

Acid soils limit yields of many important crops including canola (Brassica napus ), Australia's third largest crop. Aluminium (Al3+ ) stress is the main cause of this limitation primarily because the toxic Al3+ present inhibits root growth. Breeding programmes do not target acid-soil tolerance in B. napus because genetic variation and convincing quantitative trait loci have not been reported. We conducted a genome-wide association study (GWAS) using the BnASSYST diversity panel of B. napus genotyped with 35 729 high-quality DArTseq markers. We screened 352 B. napus accessions in hydroponics with and without a toxic concentration of AlCl3 (12μM, pH 4.3) for 12days and measured shoot biomass, root biomass, and root length. By accounting for both population structure and kinship matrices, five significant quantitative trait loci for different measures of resistance were identified using incremental Al3+ resistance indices. Within these quantitative trait locus regions of B. napus , 40 Arabidopsis thaliana gene orthologues were identified, including some previously linked with Al3+ resistance. GWAS analysis indicated that multiple genes are responsible for the natural variation in Al3+ resistance in B. napus . The results provide new genetic resources and markers to enhance that Al3+ resistance of B. napus germplasm via genomic and marker-assisted selection.

Advances and Challenges for QTL Analysis and GWAS in the Plant-Breeding of High-Yielding: A Focus on Rapeseed

Yield is one of the most important agronomic traits for the breeding of rapeseed (Brassica napus L), but its genetic dissection for the formation of high yield remains enigmatic, given the rapid population growth. In the present review, we review the discovery of major loci underlying important agronomic traits and the recent advancement in the selection of complex traits. Further, we discuss the benchmark summary of high-throughput techniques for the high-resolution genetic breeding of rapeseed. Biparental linkage analysis and association mapping have become powerful strategies to comprehend the genetic architecture of complex agronomic traits in crops. The generation of improved crop varieties, especially rapeseed, is greatly urged to enhance yield productivity. In this sense, the whole-genome sequencing of rapeseed has become achievable to clone and identify quantitative trait loci (QTLs). Moreover, the generation of high-throughput sequencing and genotyping techniques has significantly enhanced the precision of QTL mapping and genome-wide association study (GWAS) methodologies. Furthermore, this study demonstrates the first attempt to identify novel QTLs of yield-related traits, specifically focusing on ovule number per pod (ON). We also highlight the recent breakthrough concerning single-locus-GWAS (SL-GWAS) and multi-locus GWAS (ML-GWAS), which aim to enhance the potential and robust control of GWAS for improved complex traits.

Uncovering Trait Associations Resulting in Maximal Seed Yield in Winter and Spring Oilseed Rape

Seed yield is a complex trait for many crop species including oilseed rape (OSR) (Brassica napus), the second most important oilseed crop worldwide. Studies have focused on the contribution of distinct factors in seed yield such as environmental cues, agronomical practices, growth conditions, or specific phenotypic traits at the whole plant level, such as number of pods in a plant. However, how female reproductive traits contribute to whole plant level traits, and hence to seed yield, has been largely ignored. Here, we describe the combined contribution of 33 phenotypic traits within a B. napus diversity set population and their trade-offs at the whole plant and organ level, along with their interaction with plant level traits. Our results revealed that both Winter OSR (WOSR) and Spring OSR (SOSR); the two more economically important OSR groups in terms of oil production; share a common dominant reproductive strategy for seed yield. In this strategy, the main inflorescence is the principal source of seed yield, producing a good number of ovules, a large number of long pods with a concomitantly high number of seeds per pod. Moreover, we observed that WOSR opted for additional reproductive strategies than SOSR, presenting more plasticity to maximise seed yield. Overall, we conclude that OSR adopts a key strategy to ensure maximal seed yield and propose an ideal ideotype highlighting crucial phenotypic traits that could be potential targets for breeding.

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.

Using wild relatives and related species to build climate resilience in Brassica crops

Climate change will have major impacts on crop production: not just increasing drought and heat stress, but also increasing insect and disease loads and the chance of extreme weather events and further adverse conditions. Often, wild relatives show increased tolerances to biotic and abiotic stresses, due to reduced stringency of selection for yield and yield-related traits under optimum conditions. One possible strategy to improve resilience in our modern-day crop cultivars is to utilize wild relative germplasm in breeding, and attempt to introgress genetic factors contributing to greater environmental tolerances from these wild relatives into elite crop types. However, this approach can be difficult, as it relies on factors such as ease of hybridization and genetic distance between the source and target, crossover frequencies and distributions in the hybrid, and ability to select for desirable introgressions while minimizing linkage drag. In this review, we outline the possible effects that climate change may have on crop production, introduce the Brassica crop species and their wild relatives, and provide an index of useful traits that are known to be present in each of these species that may be exploitable through interspecific hybridization-based approaches. Subsequently, we outline how introgression breeding works, what factors affect the success of this approach, and how this approach can be optimized so as to increase the chance of recovering the desired introgression lines. Our review provides a working guide to the use of wild relatives and related crop germplasm to improve biotic and abiotic resistances in Brassica crop species.

Using wild relatives and related species to build climate resilience in Brassica crops

Climate change will have major impacts on crop production: not just increasing drought and heat stress, but also increasing insect and disease loads and the chance of extreme weather events and further adverse conditions. Often, wild relatives show increased tolerances to biotic and abiotic stresses, due to reduced stringency of selection for yield and yield-related traits under optimum conditions. One possible strategy to improve resilience in our modern-day crop cultivars is to utilize wild relative germplasm in breeding, and attempt to introgress genetic factors contributing to greater environmental tolerances from these wild relatives into elite crop types. However, this approach can be difficult, as it relies on factors such as ease of hybridization and genetic distance between the source and target, crossover frequencies and distributions in the hybrid, and ability to select for desirable introgressions while minimizing linkage drag. In this review, we outline the possible effects that climate change may have on crop production, introduce the Brassica crop species and their wild relatives, and provide an index of useful traits that are known to be present in each of these species that may be exploitable through interspecific hybridization-based approaches. Subsequently, we outline how introgression breeding works, what factors affect the success of this approach, and how this approach can be optimized so as to increase the chance of recovering the desired introgression lines. Our review provides a working guide to the use of wild relatives and related crop germplasm to improve biotic and abiotic resistances in Brassica crop species.

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.

Genome wide association analyses to understand genetic basis of flowering and plant height under three levels of nitrogen application in Brassica juncea (L.) Czern & Coss

Timely transition to flowering, maturity and plant height are important for agronomic adaptation and productivity of Indian mustard (B. juncea), which is a major edible oilseed crop of low input ecologies in Indian subcontinent. Breeding manipulation for these traits is difficult because of the involvement of multiple interacting genetic and environmental factors. Here, we report a genetic analysis of these traits using a population comprising 92 diverse genotypes of mustard. These genotypes were evaluated under deficient (N75), normal (N100) or excess (N125) conditions of nitrogen (N) application. Lower N availability induced early flowering and maturity in most genotypes, while high N conditions delayed both. A genotyping-by-sequencing approach helped to identify 406,888 SNP markers and undertake genome wide association studies (GWAS). 282 significant marker-trait associations (MTA's) were identified. We detected strong interactions between GWAS loci and nitrogen levels. Though some trait associated SNPs were detected repeatedly across fertility gradients, majority were identified under deficient or normal levels of N applications. Annotation of the genomic region (s) within ± 50 kb of the peak SNPs facilitated prediction of 30 candidate genes belonging to light perception, circadian, floral meristem identity, flowering regulation, gibberellic acid pathways and plant development. These included over one copy each of AGL24, AP1, FVE, FRI, GID1A and GNC. FLC and CO were predicted on chromosomes A02 and B08 respectively. CDF1, CO, FLC, AGL24, GNC and FAF2 appeared to influence the variation for plant height. Our findings may help in improving phenotypic plasticity of mustard across fertility gradients through marker-assisted breeding strategies.

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

Winter, spring and biennial varieties of Brassica napus that vary in vernalization requirement are grown for vegetable and oil production. Here, we show that the obligate or facultative nature of the vernalization requirement in European winter oilseed rape is determined by allelic variation at a 10 Mbp region on chromosome A02. This region includes orthologues of the key floral regulators FLOWERING LOCUS C (BnaFLC.A02) and FLOWERING LOCUS T (BnaFT.A02). Polymorphism at BnaFLC.A02 and BnaFT.A02, mostly in cis-regulatory regions, results in distinct gene expression dynamics in response to vernalization treatment. Our data suggest allelic variation at BnaFT.A02 is associated with flowering time in the absence of vernalization, while variation at BnaFLC.A02 is associated with flowering time under vernalizing conditions. We hypothesize selection for BnaFLC.A02 and BnaFT.A02 gene expression variation has facilitated the generation of European winter oilseed rape varieties that are adapted to different winter climates. This knowledge will allow for the selection of alleles of flowering time regulators that alter the vernalization requirement of oilseed rape, informing the generation of new varieties with adapted flowering times and improved yields.

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.

Alternative splicing of flowering time gene FT is associated with halving of time to flowering in coconut

Coconut palm has two distinct types-"tall" and "dwarf"-which differ morphologically. Tall coconut varieties need 8-10 years to start flowering, while dwarf coconut varieties only require 3-5 years. We compared seedling and reproductive stage transcriptomes for both coconut types to determine potential molecular mechanisms underlying control of flowering time in coconut. Several key genes in the photoperiod pathway were differentially expressed between seedling and reproductive leaf samples in both tall and dwarf coconut. These genes included suppressor of overexpression of constans (SOC1), flowering locus T (FT), and Apetala 1 (AP1). Alternative splicing analysis of genes in the photoperiod pathway further revealed that the FT gene produces different transcripts in tall compared to dwarf coconut. The shorter alternative splice variant of FT [which included a 6 bp deletion, alternative 3' splicing sites (A3SS)] was found to be exclusively present in dwarf coconut varieties but absent in most tall coconut varieties. Our results provide a valuable information resource as well as suggesting a probable mechanism for differentiation of flowering time onset in coconut, providing a target for future breeding work in accelerating time to flowering in this crop species.

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.

Integrating GWAS, linkage mapping and gene expression analyses reveals the genetic control of growth period traits in rapeseed (Brassica napus L.)

BACKGROUND

Brassica napus is one of the most important oilseed crops, and also an important biofuel plant due to its low air pollution and renewability. Growth period are important traits that affect yield and are crucial for its adaptation to different environments in B. napus.

RESULTS

To elucidate the genetic basis of growth period traits, genome-wide association analysis (GWAS) and linkage mapping were employed to detect the quantitative trait loci (QTL) for days to initial flowering (DIF), days to final flowering (DFF), flowering period (FP), maturity time (MT), and whole growth period (GP). A total of 146 SNPs were identified by association mapping, and 83 QTLs were identified by linkage mapping using the RIL population. Among these QTLs, 19 were pleiotropic SNPs related to multiple traits, and six (q18DFF.A03-2, q18MT.A03-2, q17DFF.A05-1, q18FP.C04, q17DIF.C05 and q17GP.C09) were consistently detected using both mapping methods. Additionally, we performed RNA sequencing to analyze the differential expression of gene (DEG) transcripts between early- and late-flowering lines selected from the RIL population, and the DEGs were integrated with association mapping and linkage analysis to confirm their roles in the growth period. Consequently, 12 candidate genes associated with growth period traits were identified in B. napus. Among these genes, seven have polymorphic sites in the coding sequence and the upstream 2-kb sequence based on the resequencing data. The haplotype BnaSOC1.A05-Haplb and BnaLNK2.C06-Hapla showed more favorable phenotypic traits.

CONCLUSIONS

The candidate genes identified in this study will contribute to our genetic understanding of growth period traits and can be used as targets for target mutations or marker-assisted breeding for rapeseed adapted to different environments.

Vernalization and Floral Transition in Autumn Drive Winter Annual Life History in Oilseed Rape

Plants with winter annual life history germinate in summer or autumn and require a period of prolonged winter cold to initiate flowering, known as vernalization. In the Brassicaceae, the requirement for vernalization is conferred by high expression of orthologs of the FLOWERING LOCUS C (FLC) gene, the expression of which is known to be silenced by prolonged exposure to winter-like temperatures [1]. Based on a wealth of vernalization experiments, typically carried out in the range of 5°C-10°C, we would expect field environments during winter to induce flowering in crops with winter annual life history. Here, we show that, in the case of winter oilseed rape, expression of multiple FLC orthologs declines not during winter but predominantly during October when the average air temperature is 10°C-15°C. We further demonstrate that plants proceed through the floral transition in early November and overwinter as inflorescence meristems, which complete floral development in spring. To validate the importance of pre-winter temperatures in flowering time control, we artificially simulated climate warming in field trial plots in October. We found that increasing the temperature by 5°C in October results in raised FLC expression and delays the floral transition by 3 weeks but only has a mild effect on flowering date the following spring. Our work shows that winter annuals overwinter as a floral bud in a manner that resembles perennials and highlights the importance of studying signaling events in the field for understanding how plants transition to flowering under real environmental conditions.

QTL analysis and candidate gene identification for plant height in cotton based on an interspecific backcross inbred line population of Gossypium hirsutum × Gossypium barbadense

We constructed the first high-quality and high-density genetic linkage map for an interspecific BIL population in cotton by specific-locus amplified fragment sequencing for QTL mapping. A novel gene GhPIN3 for plant height was identified in cotton. Ideal plant height (PH) is important for improving lint yield and mechanized harvesting in cotton. Most published genetic studies on cotton have focused on fibre yield and quality traits rather than PH. To facilitate the understanding of the genetic basis in PH, an interspecific backcross inbred line (BIL) population of 250 lines derived from upland cotton (Gossypium hirsutum L.) CRI36 and Egyptian cotton (G. barbadense L.) Hai7124 was used to construct a high-density genetic linkage map for quantitative trait locus (QTL) mapping. The high-density genetic map harboured 7,709 genotyping-by-sequencing (GBS)-based single nucleotide polymorphism (SNP) markers that covered 3,433.24 cM with a mean marker interval of 0.67 cM. In total, ten PH QTLs were identified and each explained 4.27-14.92% of the phenotypic variation, four of which were stable as they were mapped in at least two tests or based on best linear unbiased prediction in seven field tests. Based on functional annotation of orthologues in Arabidopsis and transcriptome data for the genes within the stable QTL regions, GhPIN3 encoding for the hormone auxin efflux carrier protein was identified as a candidate gene located in the stable QTL qPH-Dt1-1 region. A qRT-PCR analysis showed that the expression level of GhPIN3 in apical tissues was significantly higher in four short-statured cotton genotypes than that in four tall-statured cotton genotypes. Virus-induced gene silencing cotton has significantly increased PH when the expression of the GhPIN3 gene was suppressed.

Different copies of SENSITIVITY TO RED LIGHT REDUCED 1 show strong subfunctionalization in Brassica napus

BACKGROUND

Correct timing of flowering is critical for plants to produce enough viable offspring. In Arabidopsis thaliana (Arabidopsis), flowering time is regulated by an intricate network of molecular signaling pathways. Arabidopsis srr1-1 mutants lacking SENSITIVITY TO RED LIGHT REDUCED 1 (SRR1) expression flower early, particularly under short day (SD) conditions (1). SRR1 ensures that plants do not flower prematurely in such non-inductive conditions by controlling repression of the key florigen FT. Here, we have examined the role of SRR1 in the closely related crop species Brassica napus.

RESULTS

Arabidopsis SRR1 has five homologs in Brassica napus. They can be divided into two groups, where the A02 and C02 copies show high similarity to AtSRR1 on the protein level. The other group, including the A03, A10 and C09 copies all carry a larger deletion in the amino acid sequence. Three of the homologs are expressed at detectable levels: A02, C02 and C09. Notably, the gene copies show a differential expression pattern between spring and winter type accessions of B. napus. When the three expressed gene copies were introduced into the srr1-1 background, only A02 and C02 were able to complement the srr1-1 early flowering phenotype, while C09 could not. Transcriptional analysis of known SRR1 targets in Bna.SRR1-transformed lines showed that CYCLING DOF FACTOR 1 (CDF1) expression is key for flowering time control via SRR1.

CONCLUSIONS

We observed subfunctionalization of the B. napus SRR1 gene copies, with differential expression between early and late flowering accessions of some Bna.SRR1 copies. This suggests involvement of Bna.SRR1 in regulation of seasonal flowering in B. napus. The C09 gene copy was unable to complement srr1-1 plants, but is highly expressed in B. napus, suggesting specialization of a particular function. Furthermore, the C09 protein carries a deletion which may pinpoint a key region of the SRR1 protein potentially important for its molecular function. This is important evidence of functional domain annotation in the highly conserved but unique SRR1 amino acid sequence.

GWAS hints at pleiotropic roles for FLOWERING LOCUS T in flowering time and yield-related traits in canola

Background

Transition to flowering at the right time is critical for local adaptation and to maximize grain yield in crops. Canola is an important oilseed crop with extensive variation in flowering time among varieties. However, our understanding of underlying genes and their role in canola productivity is limited.

Results

We report our analyses of a diverse GWAS panel (300–368 accessions) of canola and identify SNPs that are significantly associated with variation in flowering time and response to photoperiod across multiple locations. We show that several of these associations map in the vicinity of FLOWERING LOCUS T (FT) paralogs and its known transcriptional regulators. Complementary QTL and eQTL mapping studies, conducted in an Australian doubled haploid population, also detected consistent genomic regions close to the FT paralogs associated with flowering time and yield-related traits. FT sequences vary between accessions. Expression levels of FT in plants grown in field (or under controlled environment cabinets) correlated with flowering time. We show that markers linked to the FT paralogs display association with variation in multiple traits including flowering time, plant emergence, shoot biomass and grain yield.

Conclusions

Our findings suggest that FT paralogs not only control flowering time but also modulate yield-related productivity traits in canola.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5964-y) contains supplementary material, which is available to authorized users.

Transcriptome sequencing assisted discovery and computational analysis of novel SNPs associated with flowering in Raphanus sativus in-bred lines for marker-assisted backcross breeding

The sequencing of radish genome aids in the better understanding and tailoring of traits associated with economic importance. In order to accelerate the genomics assisted breeding and genetic selection, transcriptomes of 33 radish inbred lines with diverse traits were sequenced for the development of single nucleotide polymorphic (SNP) markers. The sequence reads ranged from 2,560,543,741 bp to 20,039,688,139 bp with the GC (%) of 47.80-49.34 and phred quality score (Q30) of 96.47-97.54%. A total of 4951 polymorphic SNPs were identified among the accessions after stringent filtering and 298 SNPs with efficient marker assisted backcross breeding (MAB) markers were generated from the polymorphic SNPs. Further, functional annotations of SNPs revealed the effects and importance of the SNPs identified in the flowering process. The SNPs were predominantly associated with the four major flowering related transcription factors such as MYB, MADS box (AG), AP2/EREB, and bHLH. In addition, SNPs in the vital flowering integrator gene (FT) and floral repressors (EMBRYONIC FLOWER 1, 2, and FRIGIDA) were identified among the radish inbred lines. Further, 50 SNPs were randomly selected from 298 SNPs and validated using Kompetitive Allele Specific PCR genotyping system (KASP) in 102 radish inbred lines. The homozygosity of the inbred lines varied from 56 to 96% and the phylogenetic analysis resulted in the clustering of inbred lines into three subgroups. Taken together, the SNP markers identified in the present study can be utilized for the discrimination, seed purity test, and adjusting parental combinations for breeding in radish.

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.

An APETALA1 ortholog affects plant architecture and seed yield component in oilseed rape (Brassica napus L.)

BACKGROUND

Increasing the productivity of rapeseed as one of the widely cultivated oil crops in the world is of upmost importance. As flowering time and plant architecture play a key role in the regulation of rapeseed yield, understanding the genetic mechanism underlying these traits can boost the rapeseed breeding. Meristem identity genes are known to have pleiotropic effects on plant architecture and seed yield in various crops. To understand the function of one of the meristem identity genes, APETALA1 (AP1) in rapeseed, we performed phenotypic analysis of TILLING mutants under greenhouse conditions. Three stop codon mutant families carrying a mutation in Bna.AP1.A02 paralog were analyzed for different plant architecture and seed yield-related traits.

RESULTS

It was evident that stop codon mutation in the K domain of Bna.AP1.A02 paralog caused significant changes in flower morphology as well as plant architecture related traits like plant height, branch height, and branch number. Furthermore, yield-related traits like seed yield per plant and number of seeds per plants were also significantly altered in the same mutant family. Apart from phenotypic changes, stop codon mutation in K domain of Bna.AP1.A02 paralog also altered the expression of putative downstream target genes like Bna.TFL1 and Bna.FUL in shoot apical meristem (SAM) of rapeseed. Mutant plants carrying stop codon mutations in the COOH domain of Bna.AP1.A02 paralog did not have a significant effect on plant architecture, yield-related traits or the expression of the downstream targets.

CONCLUSIONS

We found that Bna.AP1.A02 paralog has pleiotropic effect on plant architecture and yield-related traits in rapeseed. The allele we found in the current study with a beneficial effect on seed yield can be incorporated into rapeseed breeding pool to develop new varieties.

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.

Photoperiod response-related gene SiCOL1 contributes to flowering in sesame

BACKGROUND

Sesame is a major oilseed crop which is widely cultivated all around the world. Flowering, the timing of transition from vegetative to reproductive growth, is one of the most important events in the life cycle of sesame. Sesame is a typical short-day (SD) plant and its flowering is largely affected by photoperiod. However, the flowering mechanism in sesame at the molecular level is still not very clear. Previous studies showed that the CONSTANS (CO) gene is the crucial photoperiod response gene which plays a center role in duration of the plant vegetative growth.

RESULTS

In this study, the CO-like (COL) genes were identified and characterized in the sesame genome. Two homologs of the CO gene in the SiCOLs, SiCOL1 and SiCOL2, were recognized and comprehensively analyzed. However, sequence analysis showed that SiCOL2 lacked one of the B-box motifs. In addition, the flowering time of the transgenic Arabidopsis lines with overexpressed SiCOL2 were longer than that of SiCOL1, indicating that SiCOL1 was more likely to be the potential functional homologue of CO in sesame. Expression analysis revealed that SiCOL1 had high expressed levels before flowering in leaves and exhibited diurnal rhythmic expression in both SD and long-day (LD) conditions. In total, 16 haplotypes of SiCOL1 were discovered in the sesame collections from Asia. However, the mutated haplotypes did not express under both SD and LD conditions and was regarded as a nonfunctional allele. Notably, the sesame landraces from high-latitude regions harboring nonfunctional alleles of SiCOL1 flowered much earlier than landraces from low-latitude regions under LD condition, and adapted to the northernmost regions of sesame cultivation. The result indicated that sesame landraces from high-latitude regions might have undergone artificial selection to adapt to the LD environment.

CONCLUSIONS

Our results suggested that SiCOL1 might contribute to regulation of flowering in sesame and natural variations in SiCOL1 were probably related to the expansion of sesame cultivation to high-latitude regions. The results could be used in sesame breeding and in broadening adaptation of sesame varieties to new regions.

Finding invisible quantitative trait loci with missing data

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

Mapping QTL controlling agronomic traits in a doubled haploid population of winter oilseed rape (Brassica napus L.)

Identification of superior alleles for agronomic traits in genetic resources of oilseed rape (Brassica napus L.) would be useful for improving the performance of locally adapted cultivars in Iran. The objective of the present work was to analyse the genetic variation and inheritance of important agronomic traits in a doubled haploid population derived from a cross between two German oilseed rape cultivars, Sansibar and Oase. Field experiments were performed in 2016-2017 with 200 doubled haploid lines and the parental genotypes applying an alpha-lattice design with two replicates. Phenological traits were recorded during the cultivation period and at maturity, seed yield, yield components and seed quality traits were determined. Significant genetic variation was found in most of the traits and heritabilities ranged from medium (48.5%) for days to end of flowering to high (92.6%) for oil content. A molecular marker linkage map was used to map 36 QTL for different traits on 17 linkage groups. Between three and four QTL were identified for each seed yield, seed weight, oil and protein content. Some of the plant material and positive QTL alleles identified for agronomic traits may be useful for improving those characters in locally adapted cultivars in Iran.

Whole-transcriptome analysis reveals genetic factors underlying flowering time regulation in rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.), one of the most important sources of vegetable oil and protein-rich meals worldwide, is adapted to different geographical regions by modification of flowering time. Rapeseed cultivars have different day length and vernalization requirements, which categorize them into winter, spring, and semiwinter ecotypes. To gain a deeper insight into genetic factors controlling floral transition in B. napus, we performed RNA sequencing (RNA-seq) in the semiwinter doubled haploid line, Ningyou7, at different developmental stages and temperature regimes. The expression profiles of more than 54,000 gene models were compared between different treatments and developmental stages, and the differentially expressed genes were considered as targets for association analysis and genetic mapping to confirm their role in floral transition. Consequently, 36 genes with association to flowering time, seed yield, or both were identified. We found novel indications for neofunctionalization in homologs of known flowering time regulators like VIN3 and FUL. Our study proved the potential of RNA-seq along with association analysis and genetic mapping to identify candidate genes for floral transition in rapeseed. The candidate genes identified in this study could be subjected to genetic modification or targeted mutagenesis and genotype building to breed rapeseed adapted to certain environments.

Sequence variation and functional analysis of a FRIGIDA orthologue (BnaA3.FRI) in Brassica napus

BACKGROUND

Allelic variation at the FRIGIDA (FRI) locus is a major contributor to natural variation of flowering time and vernalization requirement in Arabidopsis thaliana. Dominant FRI inhibits flowering by activating the expression of the MADS box transcriptional repressor FLOWERING LOCUS C (FLC), which represses flowering prior to vernalization. Four FRI orthologues had been identified in the domesticated amphidiploid Brassica napus. Linkage and association studies had revealed that one of the FRI orthologues, BnaA3.FRI, contributes to flowering time variation and crop type differentiation.

RESULTS

Sequence analyses indicated that three out of the four BnaFRI paralogues, BnaA3.FRI, BnaA10.FRI and BnaC3.FRI, contained a large number of polymorphic sites. Haplotype analysis in a panel of 174 B. napus accessions using PCR markers showed that all the three paralogues had a biased distribution of haplotypes in winter type oilseed rape (P < 0.01). Association analysis indicated that only BnaA3.FRI contributes to flowering time variation in B. napus. In addition, transgenic functional complementation demonstrated that mutations in the coding sequence of BnaA3.FRI lead to weak alleles, and subsequently to flowering time variation.

CONCLUSION

This study for the first time provides a molecular basis for flowering time control by BnaA3.FRI in B. napus, and will facilitate predictive oilseed rape breeding to select varieties with favorable flowering time and better adaption to latitude and seasonal shifts due to changing climate.

Genome-wide regression models considering general and specific combining ability predict hybrid performance in oilseed rape with similar accuracy regardless of trait architecture

Genomic prediction using the Brassica 60 k genotyping array is efficient in oilseed rape hybrids. Prediction accuracy is more dependent on trait complexity than on the prediction model. In oilseed rape breeding programs, performance prediction of parental combinations is of fundamental importance. Due to the phenomenon of heterosis, per se performance is not a reliable indicator for F₁-hybrid performance, and selection of well-paired parents requires the testing of large quantities of hybrid combinations in extensive field trials. However, the number of potential hybrids, in general, dramatically exceeds breeding capacity and budget. Integration of genomic selection (GS) could substantially increase the number of potential combinations that can be evaluated. GS models can be used to predict the performance of untested individuals based only on their genotypic profiles, using marker effects previously predicted in a training population. This allows for a preselection of promising genotypes, enabling a more efficient allocation of resources. In this study, we evaluated the usefulness of the Illumina Brassica 60 k SNP array for genomic prediction and compared three alternative approaches based on a homoscedastic ridge regression BLUP and three Bayesian prediction models that considered general and specific combining ability (GCA and SCA, respectively). A total of 448 hybrids were produced in a commercial breeding program from unbalanced crosses between 220 paternal doubled haploid lines and five male-sterile testers. Predictive ability was evaluated for seven agronomic traits. We demonstrate that the Brassica 60 k genotyping array is an adequate and highly valuable platform to implement genomic prediction of hybrid performance in oilseed rape. Furthermore, we present first insights into the application of established statistical models for prediction of important agronomical traits with contrasting patterns of polygenic control.

Major Co-localized QTL for Plant Height, Branch Initiation Height, Stem Diameter, and Flowering Time in an Alien Introgression Derived Brassica napus DH Population

Plant height (PH), branch initiation height (BIH), and stem diameter (SD) are three stem-related traits that play crucial roles in plant architecture and lodging resistance. Herein, we show one doubled haploid (DH) population obtained from a cross between Y689 (one Capsella bursa-pastoris derived Brassica napus intertribal introgression) and Westar (B. napus cultivar) that these traits were significantly positively correlated with one another and with flowering time (FT). Based on a high-density SNP map, a total of 102 additive quantitative trait loci (QTL) were identified across six environments. Seventy-two consensus QTL and 49 unique QTL were identified using a two-round strategy of QTL meta-analysis. Notably, a total of 19 major QTL, including 11 novel ones, were detected for these traits, which comprised two QTL clusters on chromosomes A02 and A07. Conditional QTL mapping was performed to preliminarily evaluate the genetic basis (pleiotropy or tight linkage) of the co-localized QTL. In addition, QTL by environment interactions (QEI) mapping was performed to verify the additive QTL and estimate the QEI effect. In the genomic regions of all major QTL, orthologs of the genes involved in phytohormone biosynthesis, phytohormone signaling, flower development, and cell differentiation in Arabidopsis were proposed as candidate genes. Of these, BnaA02g02560, an ortholog of Arabidopsis GASA4, was suggested as a candidate gene for PH, SD, and FT; and BnaA02g08490, an ortholog of Arabidopsis GNL, was associated with PH, BIH and FT. These results provide useful information for further genetic studies on stem-related traits and plant growth adaptation.

Earliness traits in rapeseed (Brassica napus): SNP loci and candidate genes identified by genome-wide association analysis

Life cycle timing is critical for yield and productivity of Brassica napus (rapeseed) cultivars grown in different environments. To facilitate breeding for earliness traits in rapeseed, SNP loci and underlying candidate genes associated with the timing of initial flowering, maturity and final flowering, as well as flowering period (FP) were investigated in two environments in a diversity panel comprising 300 B. napus inbred lines. Genome-wide association studies (GWAS) using 201,817 SNP markers previously developed from SLAF-seq (specific locus amplified fragment sequencing) revealed a total of 131 SNPs strongly linked (P < 4.96E-07) to the investigated traits. Of these 131 SNPs, 40 fell into confidence intervals or were physically adjacent to previously published flowering time QTL or SNPs. Phenotypic effect analysis detected 35 elite allelic variants for early maturing, and 90 for long FP. Candidate genes present in the same linkage disequilibrium blocks (r²>0.6) or in 100 kb regions around significant trait-associated SNPs were screened, revealing 57 B. napus genes (33 SNPs) orthologous to 39 Arabidopsis thaliana flowering time genes. These results support the practical and scientific value of novel large-scale SNP data generation in uncovering the genetic control of agronomic traits in B. napus, and also provide a theoretical basis for molecular marker-assisted selection of earliness breeding in rapeseed.

Genome-wide regression models considering general and specific combining ability predict hybrid performance in oilseed rape with similar accuracy regardless of trait architecture

KEY MESSAGE

Genomic prediction using the Brassica 60 k genotyping array is efficient in oilseed rape hybrids. Prediction accuracy is more dependent on trait complexity than on the prediction model. In oilseed rape breeding programs, performance prediction of parental combinations is of fundamental importance. Due to the phenomenon of heterosis, per se performance is not a reliable indicator for F₁-hybrid performance, and selection of well-paired parents requires the testing of large quantities of hybrid combinations in extensive field trials. However, the number of potential hybrids, in general, dramatically exceeds breeding capacity and budget. Integration of genomic selection (GS) could substantially increase the number of potential combinations that can be evaluated. GS models can be used to predict the performance of untested individuals based only on their genotypic profiles, using marker effects previously predicted in a training population. This allows for a preselection of promising genotypes, enabling a more efficient allocation of resources. In this study, we evaluated the usefulness of the Illumina Brassica 60 k SNP array for genomic prediction and compared three alternative approaches based on a homoscedastic ridge regression BLUP and three Bayesian prediction models that considered general and specific combining ability (GCA and SCA, respectively). A total of 448 hybrids were produced in a commercial breeding program from unbalanced crosses between 220 paternal doubled haploid lines and five male-sterile testers. Predictive ability was evaluated for seven agronomic traits. We demonstrate that the Brassica 60 k genotyping array is an adequate and highly valuable platform to implement genomic prediction of hybrid performance in oilseed rape. Furthermore, we present first insights into the application of established statistical models for prediction of important agronomical traits with contrasting patterns of polygenic control.

A user guide to the Brassica 60K Illumina Infinium™ SNP genotyping array

The Brassica napus 60K Illumina Infinium™ SNP array has had huge international uptake in the rapeseed community due to the revolutionary speed of acquisition and ease of analysis of this high-throughput genotyping data, particularly when coupled with the newly available reference genome sequence. However, further utilization of this valuable resource can be optimized by better understanding the promises and pitfalls of SNP arrays. We outline how best to analyze Brassica SNP marker array data for diverse applications, including linkage and association mapping, genetic diversity and genomic introgression studies. We present data on which SNPs are locus-specific in winter, semi-winter and spring B. napus germplasm pools, rather than amplifying both an A-genome and a C-genome locus or multiple loci. Common issues that arise when analyzing array data will be discussed, particularly those unique to SNP markers and how to deal with these for practical applications in Brassica breeding applications.

Post-polyploidisation morphotype diversification associates with gene copy number variation

Genetic models for polyploid crop adaptation provide important information relevant for future breeding prospects. A well-suited model is Brassica napus, a recent allopolyploid closely related to Arabidopsis thaliana. Flowering time is a major adaptation trait determining life cycle synchronization with the environment. Here we unravel natural genetic variation in B. napus flowering time regulators and investigate associations with evolutionary diversification into different life cycle morphotypes. Deep sequencing of 35 flowering regulators was performed in 280 diverse B. napus genotypes. High sequencing depth enabled high-quality calling of single-nucleotide polymorphisms (SNPs), insertion-deletions (InDels) and copy number variants (CNVs). By combining these data with genotyping data from the Brassica 60 K Illumina® Infinium SNP array, we performed a genome-wide marker distribution analysis across the 4 ecogeographical morphotypes. Twelve haplotypes, including Bna.FLC.A10, Bna.VIN3.A02 and the Bna.FT promoter on C02_random, were diagnostic for the diversification of winter and spring types. The subspecies split between oilseed/kale (B. napus ssp. napus) and swedes/rutabagas (B. napus ssp. napobrassica) was defined by 13 haplotypes, including genomic rearrangements encompassing copies of Bna.FLC, Bna.PHYA and Bna.GA3ox1. De novo variation in copies of important flowering-time genes in B. napus arose during allopolyploidisation, enabling sub-functionalisation that allowed different morphotypes to appropriately fine-tune their lifecycle.

Genome-Wide Association and Transcriptome Analyses Reveal Candidate Genes Underlying Yield-determining Traits in Brassica napus

Yield is one of the most important yet complex crop traits. To improve our understanding of the genetic basis of yield establishment, and to identify candidate genes responsible for yield improvement in Brassica napus, we performed genome-wide association studies (GWAS) for seven yield-determining traits [main inflorescence pod number (MIPN), branch pod number (BPN), pod number per plant (PNP), seed number per pod (SPP), thousand seed weight, main inflorescence yield (MIY), and branch yield], using data from 520 diverse B. napus accessions from two different yield environments. In total, we detected 128 significant single nucleotide polymorphisms (SNPs), 93 of which were revealed as novel by integrative analysis. A combination of GWAS and transcriptome sequencing on 21 haplotype blocks from samples pooled by four extremely high-yielding or low-yielding accessions revealed the differential expression of 14 crucial candiate genes (such as Bna.MYB83, Bna.SPL5, and Bna.ROP3) associated with multiple traits or containing multiple SNPs associated with the same trait. Functional annotation and expression pattern analyses further demonstrated that these 14 candiate genes might be important in developmental processes and biomass accumulation, thus affecting the yield establishment of B. napus. These results provide valuable information for understanding the genetic mechanisms underlying the establishment of high yield in B. napus, and lay the foundation for developing high-yielding B. napus varieties.

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.

Genome-Wide Association Study Reveals Candidate Genes for Control of Plant Height, Branch Initiation Height and Branch Number in Rapeseed (Brassica napus L.)

Plant architecture is crucial for rapeseed yield and is determined by plant height (PH), branch initiation height (BIH), branch number (BN) and leaf and inflorescence morphology. In this study, we measured three major factors (PH, BIH, and BN) in a panel of 333 rapeseed accessions across 4 years. A genome-wide association study (GWAS) was performed via Q + K model and the panel was genotyped using the 60 k Brassica Infinium SNP array. We identified seven loci for PH, four for BIH, and five for BN. Subsequently, by determining linkage disequilibrium (LD) decay associated with 38 significant SNPs, we gained 31, 15, and 17 candidate genes for these traits, respectively. We also showed that PH is significantly correlated with BIH, while no other correlation was revealed. Notably, a GA signaling gene (BnRGA) and a flowering gene (BnFT) located on chromosome A02 were identified as the most likely candidate genes associated with PH regulation. Furthermore, a meristem initiation gene (BnLOF2) and a NAC domain transcriptional factor (BnCUC3) that may be associated with BN were identified on the chromosome A07. This study reveals novel insight into the genetic control of plant architecture and may facilitate marker-based breeding for rapeseed.

Natural variation in CIRCADIAN CLOCK ASSOCIATED 1 is associated with flowering time in Brassica rapa

Flowering time is a very important agronomic trait and the development of molecular markers associated with this trait can facilitate crop breeding. CIRCADIAN CLOCK ASSOCIATED 1 (CCA1), a core oscillator component of circadian rhythms that affect metabolic pathways in plants, has been implicated in flowering time control in species of Brassica. CCA1 gene sequences from three Brassica rapa inbred lines, showing either early flowering or late flowering phenotypes, were analyzed and a high level of sequence variation was identified, especially within the fourth intron. Using this information, three PCR primer sets were designed and tested using various inbred lines of B. rapa. The usage of InDel markers was further validated by evaluation of flowering time and high resolution melting (HRM) analysis. Both methods, PCR and HRM, validated the use of newly developed markers. Additional sequence analyses of Brassica plants with diploid (AA, BB, or CC) and allotetraploid genomes further confirmed a large number of sequence polymorphisms in the CCA1 gene, including insertions/deletions in the fourth intron. Our results demonstrated that sequence variations in CCA1 can be used to develop valuable trait-related molecular markers for Brassica crop breeding.

Oilseed rape: learning about ancient and recent polyploid evolution from a recent crop species

Oilseed rape (Brassica napus) is one of our youngest crop species, arising several times under cultivation in the last few thousand years and completely unknown in the wild. Oilseed rape originated from hybridisation events between progenitor diploid species B. rapa and B. oleracea, both important vegetable species. The diploid progenitors are also ancient polyploids, with remnants of two previous polyploidisation events evident in the triplicated genome structure. This history of polyploid evolution and human agricultural selection makes B. napus an excellent model with which to investigate processes of genomic evolution and selection in polyploid crops. The ease of de novo interspecific hybridisation, responsiveness to tissue culture, and the close relationship of oilseed rape to the model plant Arabidopsis thaliana, coupled with the recent availability of reference genome sequences and suites of molecular cytogenetic and high-throughput genotyping tools, allow detailed dissection of genetic, genomic and phenotypic interactions in this crop. In this review we discuss the past and present uses of B. napus as a model for polyploid speciation and evolution in crop species, along with current and developing analysis tools and resources. We further outline unanswered questions that may now be tractable to investigation.

Genotypes, Networks, Phenotypes: Moving Toward Plant Systems Genetics

One of the central goals in biology is to understand how and how much of the phenotype of an organism is encoded in its genome. Although many genes that are crucial for organismal processes have been identified, much less is known about the genetic bases underlying quantitative phenotypic differences in natural populations. We discuss the fundamental gap between the large body of knowledge generated over the past decades by experimental genetics in the laboratory and what is needed to understand the genotype-to-phenotype problem on a broader scale. We argue that systems genetics, a combination of systems biology and the study of natural variation using quantitative genetics, will help to address this problem. We present major advances in these two mostly disconnected areas that have increased our understanding of the developmental processes of flowering time control and root growth. We conclude by illustrating and discussing the efforts that have been made toward systems genetics specifically in plants.

Breeding signature of combining ability improvement revealed by a genomic variation map from recurrent selection population in Brassica napus

Combining ability is crucial for parent selection in crop hybrid breeding. The present investigation and results had revealed the underlying genetic factors which might contribute in adequate combining ability, further assisting in enhancing heterosis and stability. Here, we conducted a large-scale analysis of genomic variation in order to define genomic regions affecting the combining ability in recurrent selection population of rapeseed. A population of 175 individuals was genotyped with the Brassica60K SNP chip. 525 hybrids were assembled with three different testers and used to evaluate the general combining ability (GCA) in three environments. By detecting the changes of the genomic variation, we identified 376 potential genome regions, spanning 3.03% of rapeseed genome which provided QTL-level resolution on potentially selected variants. More than 96% of these regions were located in the C subgenome, indicating that C subgenome had sustained stronger selection pressure in the breeding program than the A subgenome. In addition, a high level of linkage disequilibrium in rapeseed genome was detected, suggesting that marker-assisted selection for the population improvement might be easily implemented. This study outlines the evidence for high GCA on a genomic level and provided underlying molecular mechanism for recurrent selection improvement in B. napus.

Over-Expression of GmGIa-Regulated Soybean miR172a Confers Early Flowering in Transgenic Arabidopsis thaliana

Flowering is a pivotal event in the life cycle of plants. miR172 has been widely confirmed to play critical roles in flowering time control by regulating its target gene expression in Arabidopsis. However, the role of its counterpart in soybean remains largely unclear. In the present study, we found that the gma-miR172a was regulated by a GIGANTEA ortholog, GmGIa, in soybean through miRNA metabolism. The expression analysis revealed that gma-miR172a has a pattern of diurnal rhythm expression and its abundance increased rapidly as plants grew until the initiation of flowering phase in soybean. One target gene of gma-miR172a, Glyma03g33470, was predicted and verified using a modified RLM 5'-RACE (RNA ligase-mediated rapid amplification of 5' cDNA ends) assay. Overexpression of gma-miR172a exhibited an early flowering phenotype and the expression of FT, AP1 and LFY were simultaneously increased in gma-miR172a-transgenic Arabidopsis plants, suggesting that the early flowering phenotype was associated with up-regulation of these genes. The overexpression of the gma-miR172a-resistant version of Glyma03g33470 weakened early flowering phenotype in the toe1 mutant of Arabidopsis. Taken together, our results suggested that gma-miR172a played an important role in GmGIa-mediated flowering by repressing Glyma03g33470, which in turn increased the expression of FT, AP1 and LFY to promote flowering in soybean.

Association Mapping of Flowering Time QTLs and Insight into Their Contributions to Rapeseed Growth Habits

Plants have developed sophisticated systems to adapt to local conditions during evolution, domestication and natural or artificial selection. The selective pressures of these different growing conditions have caused significant genomic divergence within species. The flowering time trait is the most crucial factor because it helps plants to maintain sustainable development. Controlling flowering at appropriate times can also prevent plants from suffering from adverse growth conditions, such as drought, winter hardness, and disease. Hence, discovering the genome-wide genetic mechanisms that influence flowering time variations and understanding their contributions to adaptation should be a central goal of plant genetics and genomics. A global core collection panel with 448 inbred rapeseed lines was first planted in four independent environments, and their flowering time traits were evaluated. We then performed a genome-wide association mapping of flowering times with a 60 K SNP array for this core collection. With quality control and filtration, 20,342 SNP markers were ultimately used for further analyses. In total, 312 SNPs showed marker-trait associations in all four environments, and they were based on a threshold p-value of 4.06 × 10(-4); the 40 QTLs showed significant association with flowering time variations. To explore flowering time QTLs and genes related to growth habits in rapeseed, selection signals related to divergent habits were screened at the genome-wide level and 117 genomic regions were found. Comparing locations of flowering time QTLs and genes with these selection regions revealed that 20 flowering time QTLs and 224 flowering time genes overlapped with 24 and 81 selected regions, respectively. Based on this study, a number of marker-trait associations and candidate genes for flowering time variations in rapeseed were revealed. Moreover, we also showed that both flowering time QTLs and genes play important roles in rapeseed growth habits. These results will be applied to rapeseed breeding programs, and they will aid in our understanding of the relation between flowering time variations and growth habits in plants.

Association Mapping of Flowering Time QTLs and Insight into Their Contributions to Rapeseed Growth Habits

Plants have developed sophisticated systems to adapt to local conditions during evolution, domestication and natural or artificial selection. The selective pressures of these different growing conditions have caused significant genomic divergence within species. The flowering time trait is the most crucial factor because it helps plants to maintain sustainable development. Controlling flowering at appropriate times can also prevent plants from suffering from adverse growth conditions, such as drought, winter hardness, and disease. Hence, discovering the genome-wide genetic mechanisms that influence flowering time variations and understanding their contributions to adaptation should be a central goal of plant genetics and genomics. A global core collection panel with 448 inbred rapeseed lines was first planted in four independent environments, and their flowering time traits were evaluated. We then performed a genome-wide association mapping of flowering times with a 60 K SNP array for this core collection. With quality control and filtration, 20,342 SNP markers were ultimately used for further analyses. In total, 312 SNPs showed marker-trait associations in all four environments, and they were based on a threshold p-value of 4.06 × 10(-4); the 40 QTLs showed significant association with flowering time variations. To explore flowering time QTLs and genes related to growth habits in rapeseed, selection signals related to divergent habits were screened at the genome-wide level and 117 genomic regions were found. Comparing locations of flowering time QTLs and genes with these selection regions revealed that 20 flowering time QTLs and 224 flowering time genes overlapped with 24 and 81 selected regions, respectively. Based on this study, a number of marker-trait associations and candidate genes for flowering time variations in rapeseed were revealed. Moreover, we also showed that both flowering time QTLs and genes play important roles in rapeseed growth habits. These results will be applied to rapeseed breeding programs, and they will aid in our understanding of the relation between flowering time variations and growth habits in plants.

Genome-Wide Association Study Provides Insight into the Genetic Control of Plant Height in Rapeseed (Brassica napus L.)

Plant height is a key morphological trait of rapeseed. In this study, we measured plant height of a rapeseed population across six environments. This population contains 476 inbred lines representing the major Chinese rapeseed genepool and 44 lines from other countries. The 60K Brassica Infinium® SNP array was utilized to genotype the association panel. A genome-wide association study (GWAS) was performed via three methods, including a robust, novel, nonparametric Anderson-Darling (A-D) test. Consequently, 68 loci were identified as significantly associated with plant height (P < 5.22 × 10(-5)), and more than 70% of the loci (48) overlapped the confidence intervals of reported QTLs from nine mapping populations. Moreover, 24 GWAS loci were detected with selective sweep signals, which reflected the signatures of historical semi-dwarf breeding. In the linkage disequilibrium (LD) decay range up-and downstream of 65 loci (r (2) > 0.1), we found plausible candidates orthologous to the documented Arabidopsis genes involved in height regulation. One significant association found by GWAS colocalized with the established height locus BnRGA in rapeseed. Our results provide insights into the genetic basis of plant height in rapeseed and may facilitate marker-based breeding.