Flowering time variation in oilseed rape (Brassica napus L.) is associated with allelic variation in the FRIGIDA homologue BnaA.FRI.a

Molecular Evaluation of the Effects of FLC Homologs and Coordinating Regulators on the Flowering Responses to Vernalization in Cabbage (Brassica oleracea var. capitata) Genotypes

The flowering loci of cabbage must be understood to boost their productivity. In this study, to clarify the flowering mechanisms of cabbage, we examined the three flowering repressors BoFLC1, 2 and 3, and the flowering regulators BoGI, BoCOOLAIR, and BoVIN3 of early (CAB1), middle (CAB3), and late (CAB5) flowering cabbage genotypes. Analysis of allele-specifically amplified genomic DNA and various sequence alignments demonstrated that maximal insertions and deletions influenced cabbage flowering behavior, notably in CAB3 and CAB5. Phylogenetic studies showed that BoFLC1, 2, and 3 in the CAB1, 3, and 5 genotypes had the highest homologies to other Brassica species, with CAB3 and 5 the most similar. Although CAB3 and CAB5 have comparable genetic patterns, flowering repressors and flowering regulators were investigated individually with and without vernalization to determine their minor flowering differences. The expression investigation revealed that vernalized CAB5 downregulated all BoFLC genes compared to CAB3 and, in contrast, CAB3 exhibited upregulated BoCOOLAIR. We hypothesized that the CAB3 BoFLC locus' additional insertions may have led to BoCOOLAIR overexpression and BoFLC downregulation. This study sheds light on cabbage genotypes-particularly those of CAB1 and CAB5-and suggests that structural variations in BoFLC2 and 3 bind flowering regulators, such as COOLAIR, which may affect cabbage flowering time.

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

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.

Genetic dissection of branch architecture in oilseed rape (Brassica napus L.) germplasm

Branch architecture is an important factor influencing rapeseed planting density, mechanized harvest, and yield. However, its related genes and regulatory mechanisms remain largely unknown. In this study, branch angle (BA) and branch dispersion degree (BD) were used to evaluate branch architecture. Branch angle exhibited a dynamic change from an increase in the early stage to a gradual decrease until reaching a stable state. Cytological analysis showed that BA variation was mainly due to xylem size differences in the vascular bundle of the branch junction. The phenotypic analysis of 327 natural accessions revealed that BA in six environments ranged from 24.3° to 67.9°, and that BD in three environments varied from 4.20 cm to 21.4 cm, respectively. A total of 115 significant loci were detected through association mapping in three models (MLM, mrMLM, and FarmCPU), which explained 0.53%-19.4% of the phenotypic variations. Of them, 10 loci were repeatedly detected in different environments and models, one of which qBAD.A03-2 was verified as a stable QTL using a secondary segregation population. Totally, 1066 differentially expressed genes (DEGs) were identified between branch adaxial- and abaxial- sides from four extremely large or small BA/BD accessions through RNA sequencing. These DEGs were significantly enriched in the pathways related to auxin biosynthesis and transport as well as cell extension such as indole alkaloid biosynthesis, other glycan degradation, and fatty acid elongation. Four known candidate genes BnaA02g16500D (PIN1), BnaA03g10430D (PIN2), BnaC03g06250D (LAZY1), and BnaC06g20640D (ARF17) were identified by both GWAS and RNA-seq, all of which were involved in regulating the asymmetric distribution of auxins. Our identified association loci and candidate genes provide a theoretical basis for further study of gene cloning and genetic improvement of branch architecture.

Genetic mapping and genomic prediction of sclerotinia stem rot resistance to rapeseed/canola (Brassica napus L.) at seedling stage

GWAS detected ninety-eight significant SNPs associated with Sclerotinia sclerotiorum resistance. Six statistical models resulted in medium to high predictive ability, depending on trait, indicating potential of genomic prediction for disease resistance breeding. The lack of complete host resistance and a complex resistance inheritance nature between rapeseed/canola and Sclerotinia sclerotiorum often limits the development of functional molecular markers that enable breeding for sclerotinia stem rot (SSR) resistance. However, genomics-assisted selection has the potential to accelerate the breeding for SSR resistance. Therefore, genome-wide association (GWA) mapping and genomic prediction (GP) were performed using a diverse panel of 337 rapeseed/canola genotypes. Three-week-old seedlings were screened using the petiole inoculation technique (PIT). Days to wilt (DW) up to 2 weeks and lesion phenotypes (LP) at 3, 4, and 7 days post-inoculation (dpi) were recorded. A strong correlation (r = - 0.90) between DW and LP_4dpi implied that a single time point scoring at four days could be used as a proxy trait. GWA analyses using single-locus (SL) and multi-locus (ML) models identified a total of 41, and 208 significantly associated SNPs, respectively. Out of these, ninety-eight SNPs were identified by a combination of the SL model and any of the ML models, at least two ML models, or two traits. These SNPs explained 1.25-12.22% of the phenotypic variance and considered as significant, could be associated with SSR resistance. Eighty-three candidate genes with a function in disease resistance were associated with the significant SNPs. Six GP models resulted in moderate to high (0.42-0.67) predictive ability depending on SSR resistance traits. The resistant genotypes and significant SNPs will serve as valuable resources for future SSR resistance breeding. Our results also highlight the potential of genomic selection to improve rapeseed/canola breeding for SSR resistance.

Linkage disequilibrium and population structure in a core collection of Brassica napus (L.)

Estimation of genetic diversity in rapeseed is important for sustainable breeding program to provide an option for the development of new breeding lines. The objective of this study was to elucidate the patterns of genetic diversity within and among different structural groups, and measure the extent of linkage disequilibrium (LD) of 383 globally distributed rapeseed germplasm using 8,502 single nucleotide polymorphism (SNP) markers. We divided the germplasm collection into five subpopulations (P1 to P5) according to geographic and growth habit-related patterns. All subpopulations showed moderate genetic diversity (average H = 0.22 and I = 0.34). The pairwise Fst comparison revealed a great degree of divergence (Fst > 0.24) between most of the combinations. The rutabaga type showed highest divergence with spring and winter types. Higher divergence was also found between winter and spring types. Admixture model based structure analysis, principal component and neighbor-joining tree analysis placed all subpopulations into three distinct clusters. Admixed genotype constituted 29.24% of total genotypes, while remaining 70.76% belongs to identified clusters. Overall, mean linkage disequilibrium was 0.03 and it decayed to its half maximum within < 45 kb distance for whole genome. The LD decay was slower in C genome (< 93 kb); relative to the A genome (< 21 kb) which was confirmed by availability of larger haplotype blocks in C genome than A genome. The findings regarding LD pattern and population structure will help to utilize the collection as an important resource for association mapping efforts to identify genes useful in crop improvement as well as for selection of parents for hybrid breeding.

Heterosis Derived From Nonadditive Effects of the BnFLC Homologs Coordinates Early Flowering and High Yield in Rapeseed (Brassica napus L.)

Early flowering facilitates crops to adapt multiple cropping systems or growing regions with a short frost-free season; however, it usually brings an obvious yield loss. In this study, we identified that the three genes, namely, BnFLC.A2, BnFLC.C2, and BnFLC.A3b, are the major determinants for the flowering time (FT) variation of two elite rapeseed (Brassica napus L.) accessions, i.e., 616A and R11. The early-flowering alleles (i.e., Bnflc.a2 and Bnflc.c2) and late-flowering allele (i.e., BnFLC.A3b) from R11 were introgressed into the recipient parent 616A through a breeding strategy of marker-assisted backcross, giving rise to eight homozygous near-isogenic lines (NILs) associated with these three loci and 19 NIL hybrids produced by the mutual crossing of these NILs. Phenotypic investigations showed that NILs displayed significant variations in both FT and plant yield (PY). Notably, genetic analysis indicated that BnFLC.A2, BnFLC.C2, and BnFLC.A3b have additive effects of 1.446, 1.365, and 1.361 g on PY, respectively, while their dominant effects reached 3.504, 2.991, and 3.284 g, respectively, indicating that the yield loss caused by early flowering can be successfully compensated by exploring the heterosis of FT genes in the hybrid NILs. Moreover, we further validated that the heterosis of FT genes in PY was also effective in non-NIL hybrids. The results demonstrate that the exploration of the potential heterosis underlying the FT genes can coordinate early flowering (maturation) and high yield in rapeseed (B. napus L.), providing an effective strategy for early flowering breeding in crops.

Identification of WRKY transcription factors responding to abiotic stresses in Brassica napus L

A total of 278 BnWRKYs were identified and analyzed. Ectopic expression of BnWRKY149 and BnWRKY217 suggests that they function in the ABA signaling pathway. WRKY transcription factors play an important role in plant development, however, their function in Brassica napus L. abiotic stress response is still unclear. In this study, a total of 278 BnWRKY transcription factors were identified from the B. napus genome data, and they were subsequently distributed in three main groups. The protein motifs and classification of BnWRKY transcription factors were analyzed, and the locations of their corresponding encoding genes were mapped on the chromosomes of B. napus. Transcriptome analysis of rapeseed seedlings exposed to drought, salt, heat, cold and abscisic acid treatment revealed that 99 BnWRKYs responded to at least one of these stresses. The expression profiles of 12 BnWRKYs were examined with qPCR and the result coincided with RNA-seq analysis. Two genes of interest, BnWRKY149 and BnWRKY217 (homologs of AtWRKY40), were overexpressed in Arabidopsis, and the corresponding proteins were located to the nucleus. Transgene plants of BnWRKY149 and BnWRKY217 were less sensitive to ABA than Arabidopsis Col-0 plants, suggesting they might play important roles in the responses of rapeseed to abiotic stress.

SNP- and Haplotype-Based GWAS of Flowering-Related Traits in Brassica napus

Traits related to flowering time are the most promising agronomic traits that directly impact the seed yield and oil quality of rapeseed (Brassica napus L.). Developing early flowering and maturity rapeseed varieties is an important breeding objective in B. napus. Many studies have reported on days to flowering, but few have reported on budding, bolting, and the interval between bolting and DTF. Therefore, elucidating the genetic architecture of QTLs and genes regulating flowering time, we presented an integrated investigation on SNP and haplotype-based genome-wide association study of 373 diverse B. napus germplasm, which were genotyped by the 60K SNP array and were phenotyped in the four environments. The results showed that a total of 15 and 37 QTLs were detected from SNP and haplotype-based GWAS, respectively. Among them, seven QTL clusters were identified by haplotype-based GWAS. Moreover, three and eight environmentally stable QTLs were detected by SNP-GWAS and haplotype-based GWAS, respectively. By integrating the above two approaches and by co-localizing the four traits, ten (10) genomic regions were under selection on chromosomes A03, A07, A08, A10, C06, C07, and C08. Interestingly, the genomic regions FT.A07.1, FT.A08, FT.C06, and FT.C07 were identified as novel. In these ten regions, a total of 197 genes controlling FT were detected, of which 14 highly expressed DEGs were orthologous to 13 Arabidopsis thaliana genes after integration with transcriptome results. In a nutshell, the above results uncovered the genetic architecture of important agronomic traits related to flowering time and provided a basis for multiple molecular marker-trait associations in B. napus.

A 598-bp InDel Variation in the Promoter Region of Bna.SOC1.A05 Is Predominantly Present in Winter Type Rapeseeds

During rapeseed domestication and breeding, genetic diversity allowed to adapt it to different eco-geographical regions and to shape its useful traits. Structural variations (SVs), including presence/absence variations (PAVs), are thought to play a major role in the genetic diversity and phenotypic plasticity of rapeseed. In this study, we detected a 598-bp PAV within the promoter region of an Arabidopsis ortholog of a major flowering time gene and a downstream target of FLC, SOC1, which is one of the first genes that are upregulated in rapeseed during vernalization. Further analysis showed that the insertion is present predominantly in winter types while absent in spring types. The 589-bp sequence is present only in the A sub-genome indicating that it originated from Brassica rapa. Since the genomic region around Bna.SOC1.A05 showed a strong reduction in nucleotide diversity, the insertion might represent a larger selected sweep for rapeseed adaptation. Cis-element analysis showed that the insertion contains an ACGTG box, which is the strongest binding motif for the HY5 transcription factor in Arabidopsis. In addition, expression analyses showed that mRNA levels of Bna.SOC1.A05 were lower in accessions carrying the insertion compared to the ones that had no insertion.

Two Brassica napus cultivars differ in gene expression, but not in their response to submergence

Heavy rainfall causes flooding of natural ecosystems as well as farmland, negatively affecting plant performance. While the responses of the wild model organism Arabidopsis thaliana to such stress conditions is well understood, little is known about the responses of its relative, the important oil crop plant Brassica napus. For the first time, we analyzed the molecular response of Brassica napus seedlings to full submergence in a natural light-dark cycle. We used two cultivars in this study, a European hybrid cultivar and an Asian flood-tolerant cultivar. Despite their genomic differences, those genotypes showed no major differences in their responses to submergence. The molecular responses to submergence included the induction of defense- and hormone-related pathways and the repression of biosynthetic processes. Furthermore, RNAseq revealed a strong carbohydrate-starvation response under submergence in daylight, which corresponded with a fast depletion of sugars. Consequently, both B. napus cultivars exhibited a strong growth repression under water, but there was no indication of a low-oxygen response. The ability of the European hybrid cultivar to form a short-lived leaf gas film neither increased underwater net photosynthesis, underwater dark respiration nor growth during submergence. Due to the high sensitivity of both cultivars, the analysis of other cultivars or related species with higher submergence tolerance is required in order to improve flood tolerance of this crop species. One major target could be the improvement of underwater photosynthesis efficiency in order to enhance submergence survival.

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.

The ubiquitin system affects agronomic plant traits

In a single vascular plant species, the ubiquitin system consists of thousands of different proteins involved in attaching ubiquitin to substrates, recognizing or processing ubiquitinated proteins, or constituting or regulating the 26S proteasome. The ubiquitin system affects plant health, reproduction, and responses to the environment, processes that impact important agronomic traits. Here we summarize three agronomic traits influenced by ubiquitination: induction of flowering, seed size, and pathogen responses. Specifically, we review how the ubiquitin system affects expression of genes or abundance of proteins important for determining when a plant flowers (focusing on FLOWERING LOCUS C, FRIGIDA, and CONSTANS), highlight some recent studies on how seed size is affected by the ubiquitin system, and discuss how the ubiquitin system affects proteins involved in pathogen or effector recognition with details of recent studies on FLAGELLIN SENSING 2 and SUPPRESSOR OF NPR CONSTITUTIVE 1, respectively, as examples. Finally, we discuss the effects of pathogen-derived proteins on plant host ubiquitin system proteins. Further understanding of the molecular basis of the above processes could identify possible genes for modification or selection for crop improvement.

Transposon insertions within alleles of BnaFLC.A10 and BnaFLC.A2 are associated with seasonal crop type in rapeseed

In Brassicaceae, the requirement for vernalization is conferred by high expression of FLOWERING LOCUS C (FLC). The expression of FLC is known to be repressed by prolonged exposure to cold. Rapeseed (Brassica napus L.) cultivars can be classified into spring, winter, and semi-winter crop types, depending on their respective vernalization requirements. In addition to two known distinct transposon insertion events, here we identified a 4.422 kb hAT and a 5.625 kb long interspersed nuclear element transposon insertion within BnaFLC.A10, and a 810 bp miniature inverted-repeat transposable element (MITE) in BnaFLC.A2. Quantitative PCR demonstrated that these insertions lead to distinct gene expression patterns and contribute differentially to the vernalization response. Transgenic and haplotype analysis indicated that the known 621 bp MITE in the promoter region of BnaFLC.A10 is a transcriptional enhancer that appears to be the main determinant of rapeseed vernalization, and has contributed to the adaptation of rapeseed in winter cultivation environments. In the absence of this transposon insertion, the functional allele of BnaFLC.A2 is a major determinant of vernalization demand. Thus, the combination of BnaFLC.A10 carrying the 621 bp MITE insertion and a functional BnaFLC.A2 appears necessary to establish the winter rapeseed crop phenotype.

Combining QTL Mapping with Genome Resequencing Identifies an Indel in an R Gene that is Associated with Variation in Leaf Rust Disease Resistance in Poplar

Poplar trees (Populus spp.) are important and are widely grown worldwide. However, the extensive occurrence of leaf rust disease caused by Melampsora spp. seriously inhibits their growth and reduces their biomass. In our previous study, a high-quality genetic map was constructed for the poplar F1 population I-69 × XYY by using next-generation sequencing-based genotyping-by-sequencing. Here, we collected phenotypic data on leaf rust disease resistance on three different dates for all 300 progenies of the F1 population. Combining a high-quality genetic map and phenotypic data, we were able to detect 11 major quantitative trait loci (QTLs) for leaf rust disease resistance. Among these 11 QTLs, two pairs were detected on at least two dates. In the corresponding genomic sequence, we found that resistance (R) gene clusters were located in these two QTL regions. By using genome resequencing, PCR confirmation and statistical analysis, a 611-bp deletion within an R gene in one QTL region was found to be associated with variation in leaf rust disease resistance. A PCR-based examination of this 611-bp deletion was performed. This 611-bp deletion was also found to affect mRNA splicing and form a new protein with the loss of some key protein domains. Based on this study, we were able to determine the genetic architecture of variation in poplar leaf rust disease resistance, and the 611-bp deletion in the R gene could be used as a diagnostic marker for future poplar molecular breeding.

Equivocal Evidence for Colony Level Stress Effects on Bumble Bee Pollination Services

Climate change poses a threat to global food security with extreme heat events causing drought and direct damage to crop plants. However, by altering behavioural or physiological responses of insects, extreme heat events may also affect pollination services on which many crops are dependent. Such effects may potentially be exacerbated by other environmental stresses, such as exposure to widely used agro-chemicals. To determine whether environmental stressors interact to affect pollination services, we carried out field cage experiments on the buff-tailed bumble bee (Bombus terrestris). Using a Bayesian approach, we assessed whether heat stress (colonies maintained at an ambient temperature of 25 °C or 31 °C) and insecticide exposure (5 ng g-1 of the neonicotinoid insecticide clothianidin) could induce behavioural changes that affected pollination of faba bean (Vicia faba). Only the bumble bee colonies and not the plants were exposed to the environmental stress treatments. Bean plants exposed to heat-stressed bumble bee colonies (31 °C) had a lower proportional pod set compared to colonies maintained at 25 °C. There was also weak evidence that heat stressed colonies caused lower total bean weight. Bee exposure to clothianidin was found to have no clear effect on plant yields, either individually or as part of an interaction. We identified no effect of either colony stressor on bumble bee foraging behaviours. Our results suggest that extreme heat stress at the colony level may impact on pollination services. However, as the effect for other key yield parameters was weaker (e.g. bean yields), our results are not conclusive. Overall, our study highlights the need for further research on how environmental stress affects behavioural interactions in plant-pollinator systems that could impact on crop yields.

EjFRI, FRIGIDA (FRI) Ortholog from Eriobotrya japonica, Delays Flowering in Arabidopsis

In the model species Arabidopsis thaliana, FRIGIDA (FRI) is a key regulator of flowering time and can inhibit flowering without vernalization. However, little information is available on the function in the Rosaceae family. Loquat (Eriobotrya japonica) belongs to the family Rosaceae and is a distinctive species, in which flowering can be induced without vernalization, followed by blooming in late-autumn or winter. To investigate the functional roles of FRI orthologs in this non-vernalization species, we isolated an FRI ortholog, dubbed as EjFRI, from loquat. Analyses of the phylogenetic tree and protein sequence alignment showed that EjFRI is assigned to eurosids I FRI lineage. Expression analysis revealed that the highest expression level of EjFRI was after flower initiation. Meanwhile, EjFRI was widely expressed in different tissues. Subcellular localization of EjFRI was only detected to be in the nucleus. Ectopic expression of EjFRI in wild-type Arabidopsis delayed flowering time. The expression levels of EjFRI in transgenic wild-type Arabidopsis were significantly higher than those of nontransgenic wild-type lines. However, the expression levels of AtFRI showed no significant difference between transgenic and nontransgenic wild-type lines. Furthermore, the upregulated AtFLC expression in the transgenic lines indicated that EjFRI functioned similarly to the AtFRI of the model plant Arabidopsis. Our study provides a foundation to further explore the characterization of EjFRI, and also contributes to illuminating the molecular mechanism about flowering in loquat.

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.

Quantitative Trait Locus Mapping and Identification of Candidate Genes Controlling Flowering Time in Brassica napus L

Flowering time plays a vital role in determining the life-cycle period, yield, and seed quality of rapeseed (Brassica napus L.) in certain environments. Quantitative trait locus (QTL) mapping to identify the genetic architecture of genes controlling flowering time helps accelerate the early maturity breeding process. In this study, simple sequence repeats (SSR) and specific-locus amplified fragment sequencing (SLAF-seq) technologies were adopted to map the QTLs for flowering time in four environments. As a result, three target intervals, FTA09, FTA10, and FTC05 were identified. Among this, FTA09 was considered as a novel interval, FTA10 and FTC05 as stable regions. Based on the parental re-sequencing data, 7,022 single nucleotide polymorphisms (SNPs) and 2,195 insertion-deletions (InDels) between the two parents were identified in these three target regions. A total of 186 genes possessed genetic variations in these intervals, 14 of which were related to flowering time involved in photoperiod, circadian clock, vernalization, and gibberellin pathways. Six InDel markers linked to flowering time were developed in the three target intervals, indicating that the results were credible in this study. These results laid a good foundation for further genetic studies on flowering-time regulation in B. napus L.

The vernalisation regulator FLOWERING LOCUS C is differentially expressed in biennial and annual Brassica napus

Plants in temperate areas evolved vernalisation requirement to avoid pre-winter flowering. In Brassicaceae, a period of extended cold reduces the expression of the flowering inhibitor FLOWERING LOCUS C (FLC) and paves the way for the expression of downstream flowering regulators. As with all polyploid species of the Brassicaceae, the model allotetraploid Brassica napus (rapeseed, canola) is highly duplicated and carries 9 annotated copies of Bna.FLC. To investigate whether these multiple homeologs and paralogs have retained their original function in vernalisation or undergone subfunctionalisation, we compared the expression patterns of all 9 copies between vernalisation-dependent (biennial, winter type) and vernalisation-independent (annual, spring type) accessions, using RT-qPCR with copy-specific primers and RNAseq data from a diversity set. Our results show that only 3 copies – Bna.FLC.A03b, Bna.FLC.A10 and to some extent Bna.FLC.C02 – are differentially expressed between the two growth types, showing that expression of the other 6 copies does not correlate with growth type. One of those 6 copies, Bna.FLC.C03b, was not expressed at all, indicating a pseudogene, while three further copies, Bna.FLC.C03a and Bna.FLC.C09ab, did not respond to cold treatment. Sequence variation at the COOLAIR binding site of Bna.FLC.A10 was found to explain most of the variation in gene expression. However, we also found that Bna.FLC.A10 expression is not fully predictive of growth type.

The role of FRIGIDA and FLOWERING LOCUS C genes in flowering time of Brassica rapa leafy vegetables

There is a wide variation of flowering time among lines of Brassica rapa L. Most B. rapa leafy (Chinese cabbage etc.) or root (turnip) vegetables require prolonged cold exposure for flowering, known as vernalization. Premature bolting caused by low temperature leads to a reduction in the yield/quality of these B. rapa vegetables. Therefore, high bolting resistance is an important breeding trait, and understanding the molecular mechanism of vernalization is necessary to achieve this goal. In this study, we demonstrated that BrFRIb functions as an activator of BrFLC in B. rapa. We showed a positive correlation between the steady state expression levels of the sum of the BrFLC paralogs and the days to flowering after four weeks of cold treatment, suggesting that this is an indicator of the vernalization requirement. We indicate that BrFLCs are repressed by the accumulation of H3K27me3 and that the spreading of H3K27me3 promotes stable FLC repression. However, there was no clear relationship between the level of H3K27me3 in the BrFLC and the vernalization requirement. We also showed that if there was a high vernalization requirement, the rate of repression of BrFLC1 expression following prolonged cold treatments was lower.

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.

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.

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.

A Genomic Variation Map Provides Insights into the Genetic Basis of Spring Chinese Cabbage (Brassica rapa ssp. pekinensis) Selection

Chinese cabbage is the most consumed leafy crop in East Asian countries. However, premature bolting induced by continuous low temperatures severely decreases the yield and quality of the Chinese cabbage, and therefore restricts its planting season and geographic distribution. In the past 40 years, spring Chinese cabbage with strong winterness has been selected to meet the market demand. Here, we report a genome variation map of Chinese cabbage generated from the resequencing data of 194 geographically diverse accessions of three ecotypes. In-depth analyses of the selection sweeps and genome-wide patterns revealed that spring Chinese cabbage was selected from a specific population of autumn Chinese cabbage around the area of Shandong peninsula in northern China. We identified 23 genomic loci that underwent intensive selection, and further demonstrated by gene expression and haplotype analyses that the incorporation of elite alleles of VERNALISATION INSENTIVE 3.1 (BrVIN3.1) and FLOWER LOCUS C 1 (BrFLC1) is a determinant genetic source of variation during selection. Moreover, we showed that the quantitative response of BrVIN3.1 to cold due to the sequence variations in the cis elements of the BrVIN3.1 promoter significantly contributes to bolting-time variation in Chinese cabbage. Collectively, our study provides valuable insights into the genetic basis of spring Chinese cabbage selection and will facilitate the breeding of bolting-resistant varieties by molecular-marker-assisted selection, transgenic or gene editing approaches.

Construction of a high-density genetic map and its application for leaf shape QTL mapping in poplar

High-quality and dense genetic maps were constructed, and leaf shape variation was dissected by QTL mapping in poplar. Species in the genus Populus, also known as poplars, are important woody species and considered model plants for perennial trees. High-density genetic maps are valuable genomic resources for population genetics. Here, we generated a high-quality and dense genetic map for an F1 poplar population using high-throughput NGS-based genotyping. A total of 92,097 high-quality SNP markers were developed by stringent filtering and identification. In total, 889 and 1650 SNPs formed the female and male genetic maps, respectively. To test the application of the genetic maps, QTL mapping of leaf shape was conducted for this F1 population. A total of nine parameters were scored for leaf shape variation in three different environments. Combining genetic maps and measurements of the nine leaf shape parameters, we mapped a total of 42 significant QTLs. The highest LOD score of all QTLs was 9.2, and that QTL explained the most (15.13%) trait variation. A total of nine QTLs could be detected in at least two environments, and they were located in two genomic regions. Within these two QTL regions, some candidate genes for regulating leaf shape were predicted through functional annotation. The successful mapping of leaf shape QTLs demonstrated the utility of our genetic maps. According to the performance of this study, we were able to provide high-quality and dense genetic maps and dissect the leaf shape variation in poplar.

Improving Flooding Tolerance of Crop Plants

A major problem of climate change is the increasing duration and frequency of heavy rainfall events. This leads to soil flooding that negatively affects plant growth, eventually leading to death of plants if the flooding persists for several days. Most crop plants are very sensitive to flooding, and dramatic yield losses occur due to flooding each year. This review summarizes recent progress and approaches to enhance crop resistance to flooding. Most experiments have been done on maize, barley, and soybean. Work on other crops such as wheat and rape has only started. The most promising traits that might enhance crop flooding tolerance are anatomical adaptations such as aerenchyma formation, the formation of a barrier against radial oxygen loss, and the growth of adventitious roots. Metabolic adaptations might be able to improve waterlogging tolerance as well, but more studies are needed in this direction. Reasonable approaches for future studies are quantitative trait locus (QTL) analyses or genome-wide association (GWA) studies in combination with specific tolerance traits that can be easily assessed. The usage of flooding-tolerant relatives or ancestral cultivars of the crop of interest in these experiments might enhance the chances of finding useful tolerance traits to be used in breeding.

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.

Histone lysine methyltransferases BnaSDG8.A and BnaSDG8.C are involved in the floral transition in Brassica napus

Although increasing experimental evidence demonstrates that histone methylations play important roles in Arabidopsis plant growth and development, little information is available regarding Brassica napus. In this study, we characterized two genes encoding homologues of the Arabidopsis histone 3 lysine 36 (H3K36) methyltransferase SDG8, namely, BnaSDG8.A and BnaSDG8.C. Although no duplication of SDG8 homologous genes had been previously reported to occur during the evolution of any sequenced species, a domain-duplication was uncovered in BnaSDG8.C. This duplication led to the identification of a previously unknown NNH domain in the SDG8 homologues, providing a useful reference for future studies and revealing the finer mechanism of SDG8 function. One NNH domain is present in BnaSDG8.A, while two adjacent NNH domains are present in BnaSDG8.C. Reverse transcriptase-quantitative polymerase chain reaction analysis revealed similar patterns but with varied levels of expression of BnaSDG8.A/C in different plant organs/tissues. To directly investigate their function, BnaSDG8.A/C cDNA was ectopically expressed to complement the Arabidopsis mutant. We observed that the expression of either BnaSDG8.A or BnaSDG8.C could rescue the Arabidopsis sdg8 mutant to the wild-type phenotype. Using RNAi and CRISPR/Cas9-mediated gene editing, we obtained BnaSDG8.A/C knockdown and knockout mutants with the early flowering phenotype as compared with the control. Further analysis of two types of the mutants revealed that BnaSDG8.A/C are required for H3K36 m2/3 deposition and prevent the floral transition of B. napus by directly enhancing the H3K36 m2/3 levels at the BnaFLC chromatin loci. This observation on the floral transition by epigenetic modification in B. napus provides useful information for breeding early-flowering varieties.

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.

Diversity and Genome Analysis of Australian and Global Oilseed Brassica napus L. Germplasm Using Transcriptomics and Whole Genome Re-sequencing

Intensive breeding of Brassica napus has resulted in relatively low diversity, such that B. napus would benefit from germplasm improvement schemes that sustain diversity. As such, samples representative of global germplasm pools need to be assessed for existing population structure, diversity and linkage disequilibrium (LD). Complexity reduction genotyping-by-sequencing (GBS) methods, including GBS-transcriptomics (GBS-t), enable cost-effective screening of a large number of samples, while whole genome re-sequencing (WGR) delivers the ability to generate large numbers of unbiased genomic single nucleotide polymorphisms (SNPs), and identify structural variants (SVs). Furthermore, the development of genomic tools based on whole genomes representative of global oilseed diversity and orientated by the reference genome has substantial industry relevance and will be highly beneficial for canola breeding. As recent studies have focused on European and Chinese varieties, a global diversity panel as well as a substantial number of Australian spring types were included in this study. Focusing on industry relevance, 633 varieties were initially genotyped using GBS-t to examine population structure using 61,037 SNPs. Subsequently, 149 samples representative of global diversity were selected for WGR and both data sets used for a side-by-side evaluation of diversity and LD. The WGR data was further used to develop genomic resources consisting of a list of 4,029,750 high-confidence SNPs annotated using SnpEff, and SVs in the form of 10,976 deletions and 2,556 insertions. These resources form the basis of a reliable and repeatable system allowing greater integration between canola genomics studies, with a strong focus on breeding germplasm and industry applicability.

Novel and major QTL for branch angle detected by using DH population from an exotic introgression in rapeseed (Brassica napus L.)

A high-density SNP map was constructed and several novel QTL for branch angle across six environments in Brassica napus were identified. Branch angle is a major determinant for the ideotype of a plant, while the mechanisms underlying this trait in Brassica napus remain elusive. Herein, we developed one doubled haploid population from a cross involving one Capsella bursa-pastoris derived B. napus intertribal introgression line with the compressed branches and wooden stems, and constructed a high-density SNP map covering the genetic distance of 2242.14 cM, with an average marker interval of 0.73 cM. After phenotypic measurements across six environments, the inclusive composite interval mapping algorithm was conducted to analyze the QTL associated with branch angle. In single-environment analysis, a total of 17 QTL were detected and mainly distributed on chromosomes A01, A03, A09 and C03. Of these, three major QTL, qBA.A03-2, qBA.C03-3 and qBA.C03-4 were steadily expressed, each explaining more than 10% of the phenotypic variation in at least two environments. Compared with other results on rapeseed branch angle, these major QTL were newly detected. In QTL by environment interactions (QEI) mapping, 10 QTL were identified, and the QTL average effect and QEI effect were estimated. Of these, 7 QTL were detected in both single-environment analysis and QEI mapping. Based on the physical positions of SNPs and the functional annotation of the Arabidopsis thaliana genome, 27 genes within the QTL regions were selected as candidate genes, including early auxin-responsive genes, small auxin-up RNA, auxin/indoleacetic acid and gretchenhagen-3. These results may pave the way for deciphering the genetic control of branch angle in B. napus.

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.

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

QTL regions on chromosomes C06 and C09 are involved in temperature dependent time to curd induction in cauliflower. Temperature is the main environmental factor influencing curding time of cauliflower (Brassica oleracea var. botrytis). Temperatures above 20-22 °C inhibit development towards curding even in many summer cultivars. To identify quantitative trait loci (QTL) controlling curding time and its related traits in a wide range of different temperature regimes from 12 to 27 °C, a doubled haploid (DH) mapping population segregating for curding time was developed and days to curd initiation (DCI), leaf appearance rate (LAR), and final leaf number (FLN) were measured. The population was genotyped with 176 single nucleotide polymorphism (SNP) markers. Composite interval mapping (CIM) revealed repeatedly detected QTL for DCI on C06 and C09. The estimated additive effect increased at high temperatures. Significant QTL × environment interactions (Q × E) for FLN and DCI on C06 and C09 suggest that these hotspot regions have major influences on temperature mediated curd induction. 25 % of the DH lines did not induce curds at temperatures higher than 22 °C. Applying a binary model revealed a QTL with LOD >15 on C06. Nearly all lines carrying the allele of the reliable early maturing parental line (PL) on that locus induced curds at high temperatures while only half of the DH lines carrying the allele of the unreliable PL reached the generative phase during the experiment. Large variation in LAR was observed. QTL for LAR were detected repeatedly in several environments on C01, C04 and C06. Negative correlations between LAR and DCI and QTL co-localizations on C04 and C06 suggest that LAR has also effects on development towards curd induction.

Genome-wide association analyses reveal complex genetic architecture underlying natural variation for flowering time in canola

Optimum flowering time is the key to maximize canola production in order to meet global demand of vegetable oil, biodiesel and canola-meal. We reveal extensive variation in flowering time across diverse genotypes of canola under field, glasshouse and controlled environmental conditions. We conduct a genome-wide association study and identify 69 single nucleotide polymorphism (SNP) markers associated with flowering time, which are repeatedly detected across experiments. Several associated SNPs occur in clusters across the canola genome; seven of them were detected within 20 Kb regions of a priori candidate genes; FLOWERING LOCUS T, FRUITFUL, FLOWERING LOCUS C, CONSTANS, FRIGIDA, PHYTOCHROME B and an additional five SNPs were localized within 14 Kb of a previously identified quantitative trait loci for flowering time. Expression analyses showed that among FLC paralogs, BnFLC.A2 accounts for ~23% of natural variation in diverse accessions. Genome-wide association analysis for FLC expression levels mapped not only BnFLC.C2 but also other loci that contribute to variation in FLC expression. In addition to revealing the complex genetic architecture of flowering time variation, we demonstrate that the identified SNPs can be modelled to predict flowering time in diverse canola germplasm accurately and hence are suitable for genomic selection of adaptative traits in canola improvement programmes.

Comparison of statistical models for nested association mapping in rapeseed (Brassica napus L.) through computer simulations

BACKGROUND

Rapeseed (Brassica napus L.) is an important oilseed crop throughout the world, serving as source for edible oil and renewable energy. Development of nested association mapping (NAM) population and methods is of importance for quantitative trait locus (QTL) mapping in rapeseed. The objectives of the research were to compare the power of QTL detection 1- β(∗) (β(∗) is the empirical type II error rate) (i) of two mating designs, double haploid (DH-NAM) and backcross (BC-NAM), (ii) of different statistical models, and (iii) for different genetic situations.

RESULTS

The computer simulations were based on the empirical data of a single nucleotide polymorphism (SNP) set of 790 SNPs from 30 sequenced conserved genes of 51 accessions of world-wide diverse B. napus germplasm. The results showed that a joint composite interval mapping (JCIM) model had significantly higher power of QTL detection than a single marker model. The DH-NAM mating design showed a slightly higher power of QTL detection than the BC-NAM mating design. The JCIM model considering QTL effects nested within subpopulations showed higher power of QTL detection than the JCIM model considering QTL effects across subpopulations, when examing a scenario in which there were interaction effects by a few QTLs interacting with a few background markers as well as a scenario in which there were interaction effects by many QTLs (≥ 25) each with more than 10 background markers and the proportion of total variance explained by the interactions was higher than 75 %.

CONCLUSIONS

The results of our study support the optimal design as well as analysis of NAM populations, especially in rapeseed.

Allelic Variation of BnaC.TT2.a and Its Association with Seed Coat Color and Fatty Acids in Rapeseed (Brassica napus L.)

Efficient molecular markers for the selection of rapeseed genetic materials with high seed oil content and ideal fatty acid (FA) composition are preferred by rapeseed breeders. Recently, we reported the molecular mechanism of TRANSPARENT TESTA 2 (TT2) in inhibiting seed FA biosynthesis in Arabidopsis. However, evidence showing the association of rapeseed TT2 homologs and seed FA production are still insufficient. In this study, we collected 83 rapeseed (Brassica napus L.) landraces from different geographical backgrounds to conduct association mapping of BnaC.TT2.a in relation to seed coat color and FA biosynthesis. Population background was corrected by 84 pairs of SSR markers that were uniformly distributed among the linkage groups of the Tapidor-Ningyou-7 DH population. A single copy of BnaC.TT2.a for single nucleotide polymorphism (SNP) assay was cloned by a pair of previously reported specific primers. From the analysis of BnaC.TT2.a allelic variations using GLM+Q model, four SNPs on intron 1 of BnaC.TT2.a that were associated with seed FA were discovered. Moreover, an InDel at position 738 on exon 3 of BnaC.TT2.a indicated a change of protein function that was significantly associated with seed coat color, linoleic acid (C18:2), and total FA content. These findings revealed the role of BnaC.TT2.a in regulating the seed color formation and seed FA biosynthesis in rapeseed, thereby suggesting effective molecular markers for rapeseed breeding.

Agronomic and Seed Quality Traits Dissected by Genome-Wide Association Mapping in Brassica napus

In Brassica napus breeding, traits related to commercial success are of highest importance for plant breeders. However, such traits can only be assessed in an advanced developmental stage. Molecular markers genetically linked to such traits have the potential to accelerate the breeding process of B. napus by marker-assisted selection. Therefore, the objectives of this study were to identify (i) genome regions associated with the examined agronomic and seed quality traits, (ii) the interrelationship of population structure and the detected associations, and (iii) candidate genes for the revealed associations. The diversity set used in this study consisted of 405 B. napus inbred lines which were genotyped using a 6K single nucleotide polymorphism (SNP) array and phenotyped for agronomic and seed quality traits in field trials. In a genome-wide association study, we detected a total of 112 associations between SNPs and the seed quality traits as well as 46 SNP-trait associations for the agronomic traits with a P < 1.28e-05 (Bonferroni correction of α = 0.05) for the inbreds of the spring and winter trial. For the seed quality traits, a single SNP-sulfur concentration in seeds (SUL) association explained up to 67.3% of the phenotypic variance, whereas for the agronomic traits, a single SNP-blossom color (BLC) association explained up to 30.2% of the phenotypic variance. In a basic local alignment search tool (BLAST) search within a distance of 2.5 Mbp around these SNP-trait associations, 62 hits of potential candidate genes with a BLAST-score of ≥100 and a sequence identity of ≥70% to A. thaliana or B. rapa could be found for the agronomic SNP-trait associations and 187 hits of potential candidate genes for the seed quality SNP-trait associations.

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.

Identification of Flowering-Related Genes Responsible for Differences in Bolting Time between Two Radish Inbred Lines

Late bolting after cold exposure is an economically important characteristic of radish (Raphanus sativus L.), an important Brassicaceae root vegetable crop. However, little information is available regarding the genes and pathways that govern flowering time in this species. We performed high-throughput RNA sequencing analysis to elucidate the molecular mechanisms that determine the differences in flowering times between two radish lines, NH-JS1 (late bolting) and NH-JS2 (early bolting). In total, 71,188 unigenes were identified by reference-guided assembly, of which 309, 788, and 980 genes were differentially expressed between the two inbred lines after 0, 15, and 35 days of vernalization, respectively. Among these genes, 218 homologs of Arabidopsis flowering-time (Ft) genes were identified in the radish, and 49 of these genes were differentially expressed between the two radish lines in the presence or absence of vernalization treatment. Most of the Ft genes up-regulated in NH-JS1 vs. NH-JS2 were repressors of flowering, such as RsFLC, consistent with the late-bolting phenotype of NH-JS1. Although, the functions of genes down-regulated in NH-JS1 were less consistent with late-bolting characteristics than the up-regulated Ft genes, several Ft enhancer genes, including RsSOC1, a key floral integrator, showed an appropriate expression to the late-bolting phenotype. In addition, the patterns of gene expression related to the vernalization pathway closely corresponded with the different bolting times of the two inbred lines. These results suggest that the vernalization pathway is conserved between radish and Arabidopsis.

Genome-Wide Analysis of Seed Acid Detergent Lignin (ADL) and Hull Content in Rapeseed (Brassica napus L.)

A stable yellow-seeded variety is the breeding goal for obtaining the ideal rapeseed (Brassica napus L.) plant, and the amount of acid detergent lignin (ADL) in the seeds and the hull content (HC) are often used as yellow-seeded rapeseed screening indices. In this study, a genome-wide association analysis of 520 accessions was performed using the Q + K model with a total of 31,839 single-nucleotide polymorphism (SNP) sites. As a result, three significant associations on the B. napus chromosomes A05, A09, and C05 were detected for seed ADL content. The peak SNPs were within 9.27, 14.22, and 20.86 kb of the key genes BnaA.PAL4, BnaA.CAD2/BnaA.CAD3, and BnaC.CCR1, respectively. Further analyses were performed on the major locus of A05, which was also detected in the seed HC examination. A comparison of our genome-wide association study (GWAS) results and previous linkage mappings revealed a common chromosomal region on A09, which indicates that GWAS can be used as a powerful complementary strategy for dissecting complex traits in B. napus. Genomic selection (GS) utilizing the significant SNP markers based on the GWAS results exhibited increased predictive ability, indicating that the predictive ability of a given model can be substantially improved by using GWAS and GS.

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

BACKGROUND

Flowering time, plant height and seed yield are strongly influenced by climatic and day-length adaptation in crop plants. To investigate these traits under highly diverse field conditions in the important oilseed crop Brassica napus, we performed a genome-wide association study using data from diverse agroecological environments spanning three continents.

METHODS

A total of 158 European winter-type B.napus inbred lines were genotyped with 21,623 unique, single-locus single-nucleotide polymorphism (SNP) markers using the Brassica 60 K-SNP Illumina® Infinium consortium array. Phenotypic associations were calculated in the panel over the years 2010-2012 for flowering time, plant height and seed yield in 5 highly diverse locations in Germany, China and Chile, adding up to 11 diverse environments in total.

RESULTS

We identified 101 genome regions associating with the onset of flowering, 69 with plant height, 36 with seed yield and 68 cross-trait regions with potential adaptive value. Within these regions, B.napus orthologs for a number of candidate adaptation genes were detected, including central circadian clock components like CIRCADIAN CLOCK- ASSOCIATED 1 (Bna.CCA1) and the important flowering-time regulators FLOWERING LOCUS T (Bna.FT) and FRUITFUL (Bna.FUL).

DISCUSSION

Gene ontology (GO) enrichment analysis of candidate regions suggested that selection of genes involved in post-transcriptional and epigenetic regulation of flowering time may play a potential role in adaptation of B. napus to highly divergent environments. The classical flowering time regulators Bna.FLC and Bna.CO were not found among the candidate regions, although both show functional variation. Allelic effects were additive for plant height and yield, but not for flowering time. The scarcity of positive minor alleles for yield in this breeding pool points to a lack of diversity for adaptation that could restrict yield gain in the face of environmental change.

CONCLUSIONS

Our study provides a valuable framework to further improve the adaptability and yield stability of this recent allopolyploid crop under changing environments. The results suggest that flowering time regulation within an adapted B. napus breeding pool is driven by a high number of small modulating processes rather than major transcription factors like Bna.CO. In contrast, yield regulation appears highly parallel, therefore yield could be increased by pyramiding positively associated haplotypes.

How can an understanding of plant-pollinator interactions contribute to global food security?

Pollination of crops by animals is an essential part of global food production, but evidence suggests that wild pollinator populations may be declining while a number of problems are besetting managed honey bee colonies. Animal-pollinated crops grown today, bred in an environment where pollination was less likely to limit fruit set, are often suboptimal in attracting and sustaining their pollinator populations. Research into plant-pollinator interactions is often conducted in a curiosity-driven, ecological framework, but may inform breeding and biotechnological approaches to enhance pollinator attraction and crop yield. In this article we review key topics in current plant-pollinator research that have potential roles in future crop breeding for enhanced global food security.

Weight Watching and the Effect of Landscape on Honeybee Colony Productivity: Investigating the Value of Colony Weight Monitoring for the Beekeeping Industry

Over the last few decades, a gradual departure away from traditional agricultural practices has resulted in alterations to the composition of the countryside and landscapes across Europe. In the face of such changes, monitoring the development and productivity of honey bee colonies from different sites can give valuable insight on the influence of landscape on their productivity and might point towards future directions for modernized beekeeping practices. Using data on honeybee colony weights provided by electronic scales spread across Denmark, we investigated the effect of the immediate landscape on colony productivity. In order to extract meaningful information, data manipulation was necessary prior to analysis as a result of different management regimes or scales malfunction. Once this was carried out, we were able to show that colonies situated in landscapes composed of more than 50% urban areas were significantly more productive than colonies situated in those with more than 50% agricultural areas or those in mixed areas. As well as exploring some of the potential reasons for the observed differences, we discuss the value of weight monitoring of colonies on a large scale.

Genome-wide gene expression perturbation induced by loss of C2 chromosome in allotetraploid Brassica napus L

Aneuploidy with loss of entire chromosomes from normal complement disrupts the balanced genome and is tolerable only by polyploidy plants. In this study, the monosomic and nullisomic plants losing one or two copies of C2 chromosome from allotetraploid Brassica napus L. (2n = 38, AACC) were produced and compared for their phenotype and transcriptome. The monosomics gave a plant phenotype very similar to the original donor, but the nullisomics had much smaller stature and also shorter growth period. By the comparative analyses on the global transcript profiles with the euploid donor, genome-wide alterations in gene expression were revealed in two aneuploids, and their majority of differentially expressed genes (DEGs) resulted from the trans-acting effects of the zero and one copy of C2 chromosome. The higher number of up-regulated genes than down-regulated genes on other chromosomes suggested that the genome responded to the C2 loss via enhancing the expression of certain genes. Particularly, more DEGs were detected in the monosomics than nullisomics, contrasting with their phenotypes. The gene expression of the other chromosomes was differently affected, and several dysregulated domains in which up- or downregulated genes obviously clustered were identifiable. But the mean gene expression (MGE) for homoeologous chromosome A2 reduced with the C2 loss. Some genes and their expressions on C2 were correlated with the phenotype deviations in the aneuploids. These results provided new insights into the transcriptomic perturbation of the allopolyploid genome elicited by the loss of individual chromosome.

SNP markers-based map construction and genome-wide linkage analysis in Brassica napus

An Illumina Infinium array comprising 5306 single nucleotide polymorphism (SNP) markers was used to genotype 175 individuals of a doubled haploid population derived from a cross between Skipton and Ag-Spectrum, two Australian cultivars of rapeseed (Brassica napus L.). A genetic linkage map based on 613 SNP and 228 non-SNP (DArT, SSR, SRAP and candidate gene markers) covering 2514.8 cM was constructed and further utilized to identify loci associated with flowering time and resistance to blackleg, a disease caused by the fungus Leptosphaeria maculans. Comparison between genetic map positions of SNP markers and the sequenced Brassica rapa (A) and Brassica oleracea (C) genome scaffolds showed several genomic rearrangements in the B. napus genome. A major locus controlling resistance to L. maculans was identified at both seedling and adult plant stages on chromosome A07. QTL analyses revealed that up to 40.2% of genetic variation for flowering time was accounted for by loci having quantitative effects. Comparative mapping showed Arabidopsis and Brassica flowering genes such as Phytochrome A/D, Flowering Locus C and agamous-Like MADS box gene AGL1 map within marker intervals associated with flowering time in a DH population from Skipton/Ag-Spectrum. Genomic regions associated with flowering time and resistance to L. maculans had several SNP markers mapped within 10 cM. Our results suggest that SNP markers will be suitable for various applications such as trait introgression, comparative mapping and high-resolution mapping of loci in B. napus.

Arabidopsis genes, AtNPR1, AtTGA2 and AtPR-5, confer partial resistance to soybean cyst nematode (Heterodera glycines) when overexpressed in transgenic soybean roots

BACKGROUND

Extensive studies using the model system Arabidopsis thaliana to elucidate plant defense signaling and pathway networks indicate that salicylic acid (SA) is the key hormone triggering the plant defense response against biotrophic and hemi-biotrophic pathogens, while jasmonic acid (JA) and derivatives are critical to the defense response against necrotrophic pathogens. Several reports demonstrate that SA limits nematode reproduction.

RESULTS

Here we translate knowledge gained from studies using Arabidopsis to soybean. The ability of thirty-one Arabidopsis genes encoding important components of SA and JA synthesis and signaling in conferring resistance to soybean cyst nematode (SCN: Heterodera glycines) are investigated. We demonstrate that overexpression of three of thirty-one Arabidoposis genes in transgenic soybean roots of composite plants decreased the number of cysts formed by SCN to less than 50% of those found on control roots, namely AtNPR1(33%), AtTGA2 (38%), and AtPR-5 (38%). Three additional Arabidopsis genes decreased the number of SCN cysts by 40% or more: AtACBP3 (53% of the control value), AtACD2 (55%), and AtCM-3 (57%). Other genes having less or no effect included AtEDS5 (77%), AtNDR1 (82%), AtEDS1 (107%), and AtPR-1 (80%), as compared to control. Overexpression of AtDND1 greatly increased susceptibility as indicated by a large increase in the number of SCN cysts (175% of control).

CONCLUSIONS

Knowledge of the pathogen defense system gained from studies of the model system, Arabidopsis, can be directly translated to soybean through direct overexpression of Arabidopsis genes. When the genes, AtNPR1, AtGA2, and AtPR-5, encoding specific components involved in SA regulation, synthesis, and signaling, are overexpressed in soybean roots, resistance to SCN is enhanced. This demonstrates functional compatibility of some Arabidopsis genes with soybean and identifies genes that may be used to engineer resistance to nematodes.