A Tourist-like MITE insertion in the upstream region of the BnFLC.A10 gene is associated with vernalization requirement in rapeseed (Brassica napus L.)

Adaptation and the Geographic Spread of Crop Species

Crops are plant species that were domesticated starting about 11,000 years ago from several centers of origin, most prominently the Fertile Crescent, East Asia, and Mesoamerica. From their domestication centers, these crops spread across the globe and had to adapt to differing environments as a result of this dispersal. We discuss broad patterns of crop spread, including the early diffusion of crops associated with the rise and spread of agriculture, the later movement via ancient trading networks, and the exchange between the Old and New Worlds over the last ∼550 years after the European colonization of the Americas. We also examine the various genetic mechanisms associated with the evolutionary adaptation of crops to their new environments after dispersal, most prominently seasonal adaptation associated with movement across latitudes, as well as altitudinal, temperature, and other environmental factors.

The story of a decade: Genomics, functional genomics, and molecular breeding in Brassica napus

Rapeseed (Brassica napus L.) is one of the major global sources of edible vegetable oil and is also used as a feed and pioneer crop and for sightseeing and industrial purposes. Improvements in genome sequencing and molecular marker technology have fueled a boom in functional genomic studies of major agronomic characters such as yield, quality, flowering time, and stress resistance. Moreover, introgression and pyramiding of key functional genes have greatly accelerated the genetic improvement of important traits. Here we summarize recent progress in rapeseed genomics and genetics, and we discuss effective molecular breeding strategies by exploring these findings in rapeseed. These insights will extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture throughout the world.

Insertion of a miniature inverted-repeat transposable element into the promoter of OsTCP4 results in more tillers and a lower grain size in rice

A class I PCF type protein, TCP4, was identified as a transcription factor associated with both grain size and tillering through a DNA pull-down-MS assay combined with a genome-wide association study. This transcription factor was found to have a significant role in the variations among the 533 rice accessions, dividing them into two main subspecies. A Tourist-like miniature inverted-repeat transposable element (MITE) was discovered in the promoter of TCP4 in japonica/geng accessions (TCP4M+), which was found to suppress the expression of TCP4 at the transcriptional level. The MITE-deleted haplotype (TCP4M-) was mainly found in indica/xian accessions. ChIP-qPCR and EMSA demonstrated the binding of TCP4 to promoters of grain reservoir genes such as SSIIa and Amy3D in vivo and in vitro, respectively. The introduction of the genomic sequence of TCP4M+ into different TCP4M- cultivars was found to affect the expression of TCP4 in the transgenic rice, resulting in decreased expression of its downstream target gene SSIIa, increased tiller number, and decreased seed length. This study revealed that a Tourist-like MITE contributes to subspecies divergence by regulating the expression of TCP4 in response to environmental pressure, thus influencing source-sink balance by regulating starch biosynthesis in rice.

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.

BnPLP1 Positively Regulates Flowering Time, Plant Height, and Main Inflorescence Length in Brassica napus

Protein prenylation mediated by the Arabidopsis thaliana PLURIPETALA (AtPLP) gene plays a crucial role in plant growth, development, and environmental response by adding a 15-carbon farnesyl group or one to two 20-carbon geranylgeranyl groups onto one to two cysteine residues at the C-terminus of the target protein. However, the homologous genes and their functions of AtPLP in rapeseed are unclear. In this study, bioinformatics analysis and gene cloning demonstrated the existence of two homologous genes of AtPLP in the Brassica napus L. genome, namely, BnPLP1 and BnPLP2. Evolutionary analysis revealed that BnPLP1 originated from the B. rapa L. genome, while BnPLP2 originated from the B. oleracea L. genome. Genetic transformation analysis revealed that the overexpression of BnPLP1 in Arabidopsis plants exhibited earlier flowering initiation, a prolonged flowering period, increased plant height, and longer main inflorescence length compared to the wild type. Contrarily, the downregulation of BnPLP1 expression in B. napus plants led to delayed flowering initiation, shortened flowering period, decreased plant height, and reduced main inflorescence length compared to the wild type. These findings indicate that the BnPLP1 gene positively regulates flowering time, plant height, and main inflorescence length. This provides a new gene for the genetic improvement of flowering time and plant architecture in rapeseed.

Non-vernalization requirement for flowering in Brassica rapa conferred by a dominant allele of FLOWERING LOCUS T

KEY MESSAGE

We identified and characterized a dominant FT allele for flowering without vernalization in Brassica rapa, while demonstrating its potential for deployment in breeding to accelerate flowering in various Brassicaceae crops. Controlling the timing of flowering is key to improving yield and quality of several agricultural crops including the Brassicas. Many Brassicaceae crops possess a conserved flowering mechanism in which FLOWERING LOCUS C (FLC) represses the transcription of flowering activators such as FLOWERING LOCUS T (FT) during vernalization. Here, we employed genetic analysis based on next-generation sequencing to identify a dominant FT allele, BraA.FT.2-C, for flowering in the absence of vernalization in the Brassica rapa cultivar 'CHOY SUM EX CHINA 3'. BraA.FT.2-C harbors two large insertions upstream of its coding region and is expressed without vernalization, despite FLC expression. We show that BraA.FT.2-C offers an opportunity to introduce flowering without vernalization requirement into winter-type brassica crops, including B. napus, which have many functional FLC paralogs. Furthermore, we demonstrated the feasibility of using B. rapa harboring BraA.FT.2-C as rootstock for grafting to induce flowering in radish (Raphanus sativus), which requires vernalization for flowering. We believe that the ability of BraA.FT.2-C to overcome repression by FLC can have significant applications in brassica crops breeding to increase yields by accelerating or delaying flowering.

The genome of Orychophragmus violaceus provides genomic insights into the evolution of Brassicaceae polyploidization and its distinct traits

Orychophragmus violaceus, referred to as "eryuelan" (February orchid) in China, is an early-flowering ornamental plant. The high oil content and abundance of unsaturated fatty acids in O. violaceus seeds make it a potential high-quality oilseed crop. Here, we generated a whole-genome assembly for O. violaceus using Nanopore and Hi-C sequencing technologies. The assembled genome of O. violaceus was ∼1.3 Gb in size, with 12 pairs of chromosomes. Through investigation of ancestral genome evolution, we determined that the genome of O. violaceus experienced a tetraploidization event from a diploid progenitor with the translocated proto-Calepineae karyotype. Comparisons between the reconstructed subgenomes of O. violaceus identified indicators of subgenome dominance, indicating that subgenomes likely originated via allotetraploidy. O. violaceus was phylogenetically close to the Brassica genus, and tetraploidy in O. violaceus occurred approximately 8.57 million years ago, close in time to the whole-genome triplication of Brassica that likely arose via an intermediate tetraploid lineage. However, the tetraploidization in Orychophragmus was independent of the hexaploidization in Brassica, as evidenced by the results from detailed phylogenetic analyses and comparisons of the break and fusion points of ancestral genomic blocks. Moreover, identification of multi-copy genes regulating the production of high-quality oil highlighted the contributions of both tetraploidization and tandem duplication to functional innovation in O. violaceus. These findings provide novel insights into the polyploidization evolution of plant species and will promote both functional genomic studies and domestication/breeding efforts in O. violaceus.

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

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

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

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

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

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

Non-vernalization requirement in Chinese kale caused by loss of BoFLC and low expressions of its paralogs

We identified the loss of BoFLC gene as the cause of non-vernalization requirement in B. oleracea. Our developed codominant marker of BoFLC gene can be used for breeding program of B. oleracea crops. Many species of the Brassicaceae family, including some Brassica crops, require vernalization to avoid pre-winter flowering. Vernalization is an unfavorable trait for Chinese kale (Brassica oleracea var. chinensis Lei), a stem vegetable, and therefore it has been lost during its domestication/breeding process. To reveal the genetics of vernalization variation, we constructed an F₂ population through crossing a Chinese kale (a non-vernalization crop) with a kale (a vernalization crop). Using bulked segregant analysis (BSA) and RNA-seq, we identified one major quantitative trait locus (QTL) controlling vernalization and fine-mapped it to a region spanning 80 kb. Synteny analysis and PCR-based sequencing results revealed that compared to that of the kale parent, the candidate region of the Chinese kale parent lost a 9,325-bp fragment containing FLC homolog (BoFLC). In addition to the BoFLC gene, there are four other FLC homologs in the genome of B. oleracea, including Bo3g005470, Bo3g024250, Bo9g173370, and Bo9g173400. The qPCR analysis showed that the BoFLC had the highest expression among the five members of the FLC family. Considering the low expression levels of the four paralogs of BoFLC, we speculate that its paralogs cannot compensate the function of the lost BoFLC, therefore the presence/absence (PA) polymorphism of BoFLC determines the vernalization variation. Based on the PA polymorphism of BoFLC, we designed a codominant marker for the vernalization trait, which can be used for breeding programs of B. oleracea crops.

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.

Functional homoeologous alleles of CONSTANS contribute to seasonal crop type in rapeseed

Two CO paralogs in Brassica napus were confirmed and shown distinct expression pattern and function in promoting flowering and allelic variation s within BnaCO.A10 were found closely associated with ecotype divergence. CONSTANS (CO) is a key gene that responds to photoperiod and in Arabidopsis can promote flowering under long-day (LD) conditions. Brassica napus L. is a major oil crop and close relative of Arabidopsis, and arose via allopolyploidization from the diploids B. rapa (A genome) and B. oleracea (C genome). In this study, we confirmed that B. napus has two CO genes located on the A10 (BnaCO.A10) and C9 (BnaCO.C9) chromosomes. Significant differences in level and temporal pattern of transcription, as well as in protein function, of these homoeologous may have resulted from sequence variation in the promoter as well as in the coding region. Apart from two insertions of 527 bp and 2002 bp in the promoter of BnaCO.C9 that function as transcriptional enhancers, this gene is otherwise highly conserved in both promoter and coding region. However, BnaCO.A10 was classified into two haplotypes and transgene analysis in Arabidopsis and backcross analysis in rapeseed indicated that the winter-type haplotype had a greater effect in promoting flowering than the spring type. We discuss the contribution of CO alleles to species evolution, and for eco-geographic radiation following crop domestication, alongside scope for managing this locus in future breeding.

Exploring the gene pool of Brassica napus by genomics-based approaches

De novo allopolyploidization in Brassica provides a very successful model for reconstructing polyploid genomes using progenitor species and relatives to broaden crop gene pools and understand genome evolution after polyploidy, interspecific hybridization and exotic introgression. B. napus (AACC), the major cultivated rapeseed species and the third largest oilseed crop in the world, is a young Brassica species with a limited genetic base resulting from its short history of domestication, cultivation, and intensive selection during breeding for target economic traits. However, the gene pool of B. napus has been significantly enriched in recent decades that has been benefit from worldwide effects by the successful introduction of abundant subgenomic variation and novel genomic variation via intraspecific, interspecific and intergeneric crosses. An important question in this respect is how to utilize such variation to breed crops adapted to the changing global climate. Here, we review the genetic diversity, genome structure, and population-level differentiation of the B. napus gene pool in relation to known exotic introgressions from various species of the Brassicaceae, especially those elucidated by recent genome-sequencing projects. We also summarize progress in gene cloning, trait-marker associations, gene editing, molecular marker-assisted selection and genome-wide prediction, and describe the challenges and opportunities of these techniques as molecular platforms to exploit novel genomic variation and their value in the rapeseed gene pool. Future progress will accelerate the creation and manipulation of genetic diversity with genomic-based improvement, as well as provide novel insights into the neo-domestication of polyploid crops with novel genetic diversity from reconstructed genomes.

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.

Genome-wide identification of loci modifying spike-branching in tetraploid wheat

Genetic modification of spike architecture is essential for improving wheat yield. Newly identified loci for the 'Miracle wheat' phenotype on chromosomes 1AS and 2BS have significant effects on spike traits. The wheat (Triticum ssp.) inflorescence, also known as a spike, forms an unbranched inflorescence in which the inflorescence meristem generates axillary spikelet meristems (SMs) destined to become sessile spikelets. Previously, we identified the putatively causative mutation in the branched head^t (bh^t) gene (TtBH-A1) of tetraploid wheat (T. turgidum convar. compositum (L.f.) Filat.) responsible for the loss of SM identity, converting the non-branching spike to a branched wheat spike. In the current study, we performed whole-genome quantitative trait loci (QTL) analysis using 146 recombinant inbred lines (RILs) derived from a cross between spike-branching wheat ('Miracle wheat') and an elite durum wheat cultivar showing broad phenotypic variation for spike architecture. Besides the previously found gene at the bh^t-A1 locus on the short arm of chromosome 2A, we also mapped two new modifier QTL for spike-branching on the short arm of chromosome 1A, termed bh^t-A2, and 2BS. Using biparental mapping population and GWAS in 302 diverse accessions, the 2BS locus was highly associated with coding sequence variation found at the homoeo-allele of TtBH-B1 (bh^t-B1). Thus, RILs that combined both bh^t-A1 and bh^t-B1 alleles showed additive genetic effects leading to increased penetrance and expressivity of the supernumerary spikelet and/or mini-spike formation.

Distribution of MITE family Monkey King in rapeseed (Brassica napus L) and its influence on gene expression

Miniature inverted-repeat transposable elements (MITEs) are a group of class II transposable elements. The MITE Monkey King (MK) was first discovered upstream of BnFLC.A10. In this study, genome resequencing of four selected B. napus accessions, revealed more than 4000 distributed copies of MKs constituting ~2.4 Mb of the B. napus genomic sequence and caused 677 polymorphisms among the four accessions. MK -polymorphism-related markers across 128 natural and 58 synthetic accessions revealed more polymorphic MKs in natural than synthetic accessions. Ten MK -induced indels significantly affected the expression levels of the nearest gene based on RNAseq analysis, six of these effects were subsequently confirmed using qRT-PCR. Decreased expression pattern of MK -derived miRNA-bna-miR6031 was also observed under various stress treatments. Further research focused on the MITE families should promote not only our understanding of gene regulatory networks but also inform crop improvement efforts.

INDITTO2 transposon conveys auxin-mediated DRO1 transcription for rice drought avoidance

Transposable elements exist widely throughout plant genomes and play important roles in plant evolution. Auxin is an important regulator that is traditionally associated with root development and drought stress adaptation. The DEEPER ROOTING 1 (DRO1) gene is a key component of rice drought avoidance. Here, we identified a transposon that acts as an autonomous auxin-responsive promoter and its presence at specific genome positions conveys physiological adaptations related to drought avoidance. Rice varieties with a high and auxin-mediated transcription of DRO1 in the root tip show deeper and longer root phenotypes and are thus better adapted to drought. The INDITTO2 transposon contains an auxin response element and displays auxin-responsive promoter activity; it is thus able to convey auxin regulation of transcription to genes in its proximity. In the rice Acuce, which displays DRO1-mediated drought adaptation, the INDITTO2 transposon was found to be inserted at the promoter region of the DRO1 locus. Transgenesis-based insertion of the INDITTO2 transposon into the DRO1 promoter of the non-adapted rice variety Nipponbare was sufficient to promote its drought avoidance. Our data identify an example of how transposons can act as promoters and convey hormonal regulation to nearby loci, improving plant fitness in response to different abiotic stresses.

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.

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

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

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.

A copia-like retrotransposon insertion in the upstream region of the SHATTERPROOF1 gene, BnSHP1.A9, is associated with quantitative variation in pod shattering resistance in oilseed rape

Seed loss resulting from pod shattering is a major constraint in production of oilseed rape (Brassica napus L.). However, the molecular mechanisms underlying pod shatter resistance are not well understood. Here, we show that the pod shatter resistance at quantitative trait locus qSRI.A9.1 is controlled by one of the B. napus SHATTERPROOF1 homologs, BnSHP1.A9, in a doubled haploid population generated from parents designated R1 and R2 as well as in a diverse panel of oilseed rape. The R1 maternal parental line of the doubled haploid population carried the allele for shattering at qSRI.A9.1, while the R2 parental line carried the allele for shattering resistance. Quantitative RT-PCR showed that BnSHP1.A9 was expressed specifically in flower buds, flowers, and developing siliques in R1, while it was not expressed in any tissue of R2. Transgenic plants constitutively expressing either of the BnSHP1.A9 alleles from the R1 and R2 parental lines showed that both alleles are responsible for pod shattering, via a mechanism that promotes lignification of the enb layer. These findings indicated that the allelic differences in the BnSHP1.A9 gene per se are not the causal factor for quantitative variation in shattering resistance at qSRI.A9.1. Instead, a highly methylated copia-like long terminal repeat retrotransposon insertion (4803 bp) in the promotor region of the R2 allele of BnSHP1.A9 repressed the expression of BnSHP1.A9, and thus contributed to pod shatter resistance. Finally, we showed a copia-like retrotransposon-based marker, BnSHP1.A9R2, can be used for marker-assisted breeding targeting the pod shatter resistance trait in oilseed rape.

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.

Perspectives of CRISPR/Cas-mediated cis-engineering in horticulture: unlocking the neglected potential for crop improvement

Directed breeding of horticultural crops is essential for increasing yield, nutritional content, and consumer-valued characteristics such as shape and color of the produce. However, limited genetic diversity restricts the amount of crop improvement that can be achieved through conventional breeding approaches. Natural genetic changes in cis-regulatory regions of genes play important roles in shaping phenotypic diversity by altering their expression. Utilization of CRISPR/Cas editing in crop species can accelerate crop improvement through the introduction of genetic variation in a targeted manner. The advent of CRISPR/Cas-mediated cis-regulatory region engineering (cis-engineering) provides a more refined method for modulating gene expression and creating phenotypic diversity to benefit crop improvement. Here, we focus on the current applications of CRISPR/Cas-mediated cis-engineering in horticultural crops. We describe strategies and limitations for its use in crop improvement, including de novo cis-regulatory element (CRE) discovery, precise genome editing, and transgene-free genome editing. In addition, we discuss the challenges and prospects regarding current technologies and achievements. CRISPR/Cas-mediated cis-engineering is a critical tool for generating horticultural crops that are better able to adapt to climate change and providing food for an increasing world population.

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

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

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

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

Eight high-quality genomes reveal pan-genome architecture and ecotype differentiation of Brassica napus

Rapeseed (Brassica napus) is the second most important oilseed crop in the world but the genetic diversity underlying its massive phenotypic variations remains largely unexplored. Here, we report the sequencing, de novo assembly and annotation of eight B. napus accessions. Using pan-genome comparative analysis, millions of small variations and 77.2-149.6 megabase presence and absence variations (PAVs) were identified. More than 9.4% of the genes contained large-effect mutations or structural variations. PAV-based genome-wide association study (PAV-GWAS) directly identified causal structural variations for silique length, seed weight and flowering time in a nested association mapping population with ZS11 (reference line) as the donor, which were not detected by single-nucleotide polymorphisms-based GWAS (SNP-GWAS), demonstrating that PAV-GWAS was complementary to SNP-GWAS in identifying associations to traits. Further analysis showed that PAVs in three FLOWERING LOCUS C genes were closely related to flowering time and ecotype differentiation. This study provides resources to support a better understanding of the genome architecture and acceleration of the genetic improvement of B. napus.

Transcriptomic Identification of Floral Transition and Development-Associated Genes in Styrax japonicus

Styrax japonicus (S. japonicus) is an important flowering tree species in temperate regions, and it is regarded as a nectariferous plant. However, there have been few studies to date analyzing floral development in this species. In order to understand gene expression dynamics during S. japonicus flower development, we; therefore, prepared cDNA libraries from three distinct stages of S. japonicus. Illumina sequencing generated 31,471 differentially expressed unigenes during flower development. We additionally conducted pathway enrichment analyses using the GO and KEGG database in order to assess the functions of genes differentially expressed during different stages of the floral development process, revealing these genes to be associated with pathways including phytohormone signaling, Transcription factor, protein kinase, and circadian rhythms. In total, 4828 TF genes, 8402 protein kinase genes, and 78 DEGs related to hormone pathways were identified in flower development stages. Six genes were selected for confirmation of expression levels using quantitative real-time PCR. The gene expression data presented herein represent the most comprehensive dataset available regarding the flowering of S. japonicus, thus offering a reference for future studies of the flowering of this and other Styracaceae species.

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.

Fine mapping of a silique length- and seed weight-related gene in Brassica napus

Using microarray analysis combined with map-based cloning, a major locus positively regulating SL and SW was mapped to a 98.47 kb interval on A09 in rapeseed. In rapeseed, seed yield is closely associated with silique-related traits such as silique length (SL) and seed weight (SW). Previously identified quantitative trait loci (QTLs) revealed that SL and SW are complex traits and many QTLs overlap. However, the genetic characterization of the association between SL and SW is poorly understood. In the present study, a BC₃F₃ near isogenic line developed from a short silique plant and the long silique cultivar 'ZS11' was analyzed to identify the locus related to SL. Map-based cloning indicated that a major locus acting as a single Mendelian factor was mapped to a 98.47 kb region on chromosome A09. BLAST analysis and DNA sequencing showed SNP variations and a fragment replacement in the upstream region of the candidate gene BnaA09g55530D may alter gene expression and influence SL. The results showed that this SL locus may also positively affect SW as well as in the 186 rapeseed accessions identified by the associated markers. Therefore, selecting plants with appropriate SL and developing functional markers for the associated gene could play important roles in the molecular breeding of high-yield rapeseed varieties.

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.

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.

Non-coding RNAs and transposable elements in plant genomes: emergence, regulatory mechanisms and roles in plant development and stress responses

This review will provide evidence for the indispensable function of these elements in regulating plant development and resistance to biotic and abiotic stresses, as well as their evolutionary role in facilitating plant adaptation. Over millions of years of evolution, plant genomes have acquired a complex constitution. Plant genomes consist not only of protein coding sequences, but also contain large proportions of non-coding sequences. These include introns of protein-coding genes, and intergenic sequences such as non-coding RNA, repeat sequences and transposable elements. These non-coding sequences help to regulate gene expression, and are increasingly being recognized as playing an important role in genome organization and function. In this review, we summarize the known molecular mechanisms by which gene expression is regulated by several species of non-coding RNAs (microRNAs, long non-coding RNAs, and circular RNAs) and by transposable elements. We further discuss how these non-coding RNAs and transposable elements evolve and emerge in the genome, and the potential influence and importance of these non-coding RNAs and transposable elements in plant development and in stress responses.

Molecular markers based on sequence variation in BoFLC1.C9 for characterizing early- and late-flowering cabbage genotypes

BACKGROUND

Cabbage (Brassica oleracea var. capitata) is popular worldwide for consumption as a leafy vegetable. Premature flowering is triggered by low temperature, and deteriorates quality of cabbage as vegetable. In general, growers prefer late-flowering varieties to assure good quality compact head. Here, we report BoFLC1.C9 as a gene with clear sequence variation between cabbage lines with different flowering times, and proposed as molecular marker to characterize early- and late-flowering cabbage lines.

RESULTS

We identified sequence variation of 67 bp insertions in intron 2, which were contributed in flowering time variation between two inbred lines through rapid down-regulation of the BoFLC1.C9 gene in early-flowering line compared to late-flowering one upon vernalization. One set of primer 'F7R7' proposed as marker, of which was explained with 83 and 80% of flowering time variation in 141 F₂ individuals and 20 commercial lines, respectively.

CONCLUSIONS

This F7R7 marker could be used as genetic tools to characterize flowering time variation and to select as well to develop early- and late-flowering cabbage cultivars.

Joint QTL mapping and transcriptome sequencing analysis reveal candidate flowering time genes in Brassica napus L

Background

Optimum flowering time is a key agronomic trait in Brassica napus. To investigate the genetic architecture and genetic regulation of flowering time in this important crop, we conducted quantitative trait loci (QTL) analysis of flowering time in a recombinant inbred line (RIL) population, including lines with extreme differences in flowering time, in six environments, along with RNA-Seq analysis.

Results

We detected 27 QTLs distributed on eight chromosomes among six environments, including one major QTL on chromosome C02 that explained 11–25% of the phenotypic variation and was stably detected in all six environments. RNA-Seq analysis revealed 105 flowering time-related differentially expressed genes (DEGs) that play roles in the circadian clock/photoperiod, autonomous pathway, and hormone and vernalization pathways. We focused on DEGs related to the regulation of flowering time, especially DEGs in QTL regions.

Conclusions

We identified 45 flowering time-related genes in these QTL regions, eight of which are DEGs, including key flowering time genes PSEUDO RESPONSE REGULATOR 7 (PRR7) and FY (located in a major QTL region on C02). These findings provide insights into the genetic architecture of flowering time in B. napus.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5356-8) 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.

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.

A comprehensive and precise set of intervarietal substitution lines to identify candidate genes and quantitative trait loci in oilseed rape (Brassica napus L.)

A set of intervarietal substitution lines were developed in rapeseed by recurrent backcrossing and marker-assisted selection and employed for mapping both qualitative and quantitative traits. Intervarietal substitution lines (ISLs) may be assembled into advanced secondary mapping populations that have remarkable potential for resolving trait loci and mapping candidate genes. To facilitate the identification of important genes in oilseed rape (canola, Brassica napus), we developed 89 ISLs using an elite cultivar 'Zhongyou 821' (ZY821) as the recipient and a re-synthesized line 'No.2127' as the donor. In the whole process of ISLs development, the target chromosome segments were selected based on the genotypes of 300 microsatellite markers evenly distributed across the genome. Eighty-nine ISLs fixed at BC₅F₄ were genotyped by sequencing using double digestion to survey the lengths of target substitution segments from the donor parent and the background segments from the recurrent parent. The total length of the substituted chromosome segments was 3030.27 Mb, representing 3.56 × of the Darmor-bzh reference genome sequence (version 4.1). Gene mapping was conducted for two qualitative traits, flower colour and seed-coat colour, and nine quantitative traits including yield- and quality-related traits, with 19 QTLs identified for the latter. Overlapping substitution segments were identified for flower colour and seed-coat colour loci, as well as for QTLs consistently detected in 2 or 3 years. These results demonstrate the value of these ISLs for locus resolution and subsequent cloning, targeted mutation or editing of genes controlling important traits in oilseed rape.

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.

Effect of vernalization on tuberization and flowering in the Tibetan turnip is associated with changes in the expression of FLC homologues

The turnip (Brassica rapa var. rapa) is a biennial crop that is planted in late summer/early fall and forms fleshy tubers for food in temperate regions. The harvested tubers then overwinter and are planted again the next spring for flowering and seeds. FLOWERING LOCUS C (FLC) is a MADS-box transcription factor that acts as a major repressor of floral transition by suppressing the flowering promoters FT and SOC1. Here we show that vernalization effectively represses tuber formation and promotes flowering in Tibetan turnip. We functionally characterized four FLC homologues (BrrFLC1, FLC2, FLC3, and FLC5), and found that BrrFLC2 and BrrFLC1 play a major role in repressing flowering in turnip and in transgenic Arabidopsis. In contrast, tuber formation was correlated with BrrFLC1 expression in the hypocotyl and was repressed under cold treatment following the quantitative downregulation of BrrFLC1. Grafting experiments of non-vernalized and vernalized turnips revealed that vernalization independently suppressed tuberization in the tuber or hypocotyl of the rootstock or scion, which occurred in parallel with the reduction in BrrFLC1 activity. Together, our results demonstrate that the Tibetan turnip is highly responsive to cold exposure, which is associated with the expression levels of BrrFLC genes.

MADS-box genes and crop domestication: the jack of all traits

MADS-box genes are key regulators of virtually every aspect of plant reproductive development. They play especially prominent roles in flowering time control, inflorescence architecture, floral organ identity determination, and seed development. The developmental and evolutionary importance of MADS-box genes is widely acknowledged. However, their role during flowering plant domestication is less well recognized. Here, we provide an overview illustrating that MADS-box genes have been important targets of selection during crop domestication and improvement. Numerous examples from a diversity of crop plants show that various developmental processes have been shaped by allelic variations in MADS-box genes. We propose that new genomic and genome editing resources provide an excellent starting point for further harnessing the potential of MADS-box genes to improve a variety of reproductive traits in crops. We also suggest that the biophysics of MADS-domain protein-protein and protein-DNA interactions, which is becoming increasingly well characterized, makes them especially suited to exploit coding sequence variations for targeted breeding approaches.

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.

Integrated QTL and eQTL Mapping Provides Insights and Candidate Genes for Fatty Acid Composition, Flowering Time, and Growth Traits in a F2 Population of a Novel Synthetic Allopolyploid Brassica napus

Brassica napus (B. napus, AACC), is an economically important allotetraploid crop species that resulted from hybridization between two diploid species, Brassica rapa (AA) and Brassica olereacea (CC). We have created one new synthetic B. napus genotype Da-Ae (AACC) and one introgression line Da-Ol-1 (AACC), which were used to generate an F₂ mapping population. Plants in this F₂ mapping population varied in fatty acid content, flowering time, and growth-related traits. Using quantitative trait locus (QTL) mapping, we aimed to determine if Da-Ae and Da-Ol-1 provided novel genetic variation beyond what has already been found in B. napus. Making use of the genotyping information generated from RNA-seq data of these two lines and their F₂ mapping population of 166 plants, we constructed a genetic map consisting of 2,021 single nucleotide polymorphism markers that spans 2,929 cM across 19 linkage groups. Besides the known major QTL identified, our high resolution genetic map facilitated the identification of several new QTL contributing to the different fatty acid levels, flowering time, and growth-related trait values. These new QTL probably represent novel genetic variation that existed in our new synthetic B. napus strain. By conducting genome-wide expression variation analysis in our F₂ mapping population, genetic regions that potentially regulate many genes across the genome were revealed. A FLOWERING LOCUS C gene homolog, which was identified as a candidate regulating flowering time and multiple growth-related traits, was found underlying one of these regions. Integrated QTL and expression QTL analyses also helped us identified candidate causative genes associated with various biological traits through expression level change and/or possible protein function modification.

Fine-tuning timing: natural variation informs the mechanistic basis of the switch to flowering in Arabidopsis thaliana

The evolution of diverse life history strategies has allowed Arabidopsis thaliana to adapt to worldwide locations, spanning a range of latitudinal and environmental conditions. Arabidopsis thaliana accessions are either vernalization-requiring winter annuals or rapid cyclers, with extensive natural variation in vernalization requirement and response. Genetic and molecular analysis of this variation has enhanced our understanding of the mechanisms involved in life history determination, with translation to both natural and crop systems in the Brassicaceae and beyond.

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.

Changing Responses to Changing Seasons: Natural Variation in the Plasticity of Flowering Time

The mechanisms by which the environment regulates flowering time have evolved as crops and wild populations have adapted to diverse climates, and the specific variants involved are increasingly known.

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.

Lessons from Domestication: Targeting Cis-Regulatory Elements for Crop Improvement

Domestication of wild plant species has provided us with crops that serve our human nutritional needs. Advanced DNA sequencing has propelled the unveiling of underlying genetic changes associated with domestication. Interestingly, many changes reside in cis-regulatory elements (CREs) that control the expression of an unmodified coding sequence. Sequence variation in CREs can impact gene expression levels, but also developmental timing and tissue specificity of expression. When genes are involved in multiple pathways or active in several organs and developmental stages CRE modifications are favored in contrast to mutations in coding regions, due to the lack of detrimental pleiotropic effects. Therefore, learning from domestication, we propose that CREs are interesting targets for genome editing to create new alleles for plant breeding.

Identification of Polymorphisms Associated with Drought Adaptation QTL in Brassica napus by Resequencing

Brassica napus is a globally important oilseed for which little is known about the genetics of drought adaptation. We previously mapped twelve quantitative trait loci (QTL) underlying drought-related traits in a biparental mapping population created from a cross between winter and spring B. napus cultivars. Here we resequence the genomes of the mapping population parents to identify genetic diversity across the genome and within QTL regions. We sequenced each parental cultivar on the Illumina HiSeq platform to a minimum depth of 23 × and performed a reference based assembly in order to describe the molecular variation differentiating them at the scale of the genome, QTL and gene. Genome-wide patterns of variation were characterized by an overall higher single nucleotide polymorphism (SNP) density in the A genome and a higher ratio of nonsynonymous to synonymous substitutions in the C genome. Nonsynonymous substitutions were used to categorize gene ontology terms differentiating the parent genomes along with a list of putative functional variants contained within each QTL. Marker assays were developed for several of the discovered polymorphisms within a pleiotropic QTL on chromosome A10. QTL analysis with the new, denser map showed the most associated marker to be that developed from an insertion/deletion polymorphism located in the candidate gene Bna.FLC.A10, and it was the only candidate within the QTL interval with observed polymorphism. Together, these results provide a glimpse of genome-wide variation differentiating annual and biennial B. napus ecotypes as well as a better understanding of the genetic basis of root and drought phenotypes.

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.