Long-read sequencing reveals widespread intragenic structural variants in a recent allopolyploid crop plant.

Genome structural variation (SV) contributes strongly to trait variation in eukaryotic species and may have an even higher functional significance than single-nucleotide polymorphism (SNP). In recent years, there have been a number of studies associating large chromosomal scale SV ranging from hundreds of kilobases all the way up to a few megabases to key agronomic traits in plant genomes. However, there have been little or no efforts towards cataloguing small- (30-10 000 bp) to mid-scale (10 000-30 000 bp) SV and their impact on evolution and adaptation-related traits in plants. This might be attributed to complex and highly duplicated nature of plant genomes, which makes them difficult to assess using high-throughput genome screening methods. Here, we describe how long-read sequencing technologies can overcome this problem, revealing a surprisingly high level of widespread, small- to mid-scale SV in a major allopolyploid crop species, Brassica napus. We found that up to 10% of all genes were affected by small- to mid-scale SV events. Nearly half of these SV events ranged between 100 bp and 1000 bp, which makes them challenging to detect using short-read Illumina sequencing. Examples demonstrating the contribution of such SV towards eco-geographical adaptation and disease resistance in oilseed rape suggest that revisiting complex plant genomes using medium-coverage long-read sequencing might reveal unexpected levels of functional gene variation, with major implications for trait regulation and crop improvement.

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

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

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

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

"Doubled-haploid" allohexaploid Brassica lines lose fertility and viability and accumulate genetic variation due to genomic instability.

Microspore culture stimulates immature pollen grains to develop into plants via tissue culture and is used routinely in many crop species to produce "doubled haploids": homozygous, true-breeding lines. However, microspore culture is also often used on material that does not have stable meiosis, such as interspecific hybrids. In this case, the resulting progeny may lose their "doubled haploid" homozygous status as a result of chromosome missegregation and homoeologous exchanges. However, little is known about the frequency of these effects. We assessed fertility, meiosis and genetic variability in self-pollinated progeny sets (the MDL2 population) resulting from first-generation plants (the MDL1 population) derived from microspores of a near-allohexaploid interspecific hybrid from the cross (Brassica napus × B. carinata) × B. juncea. Allelic inheritance and copy number variation were predicted using single nucleotide polymorphism marker data from the Illumina Infinium 60K Brassica array. Seed fertility and viability decreased substantially from the MDL1 to the MDL2 generation. In the MDL2 population, 87% of individuals differed genetically from their MDL1 parent. These genetic differences resulted from novel homoeologous exchanges between chromosomes, chromosome loss and gain, and segregation and instability of pre-existing karyotype abnormalities. Novel karyotype change was extremely common, with 2.2 new variants observed per MDL2 individual. Significant differences between progeny sets in the number of novel genetic variants were also observed. Meiotic instability clearly has the potential to dramatically change karyotypes (often without detectable effects on the presence or absence of alleles) in putatively homozygous, microspore-derived lines, resulting in loss of fertility and viability.

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.

Inherited allelic variants and novel karyotype changes influence fertility and genome stability in Brassica allohexaploids.

Synthetic allohexaploid Brassica hybrids (2n = AABBCC) do not exist naturally, but can be synthesized by crosses between diploid and/or allotetraploid Brassica species. Using these hybrids, we aimed to identify how novel allohexaploids restore fertility and normal meiosis after formation. Chromosome inheritance, genome structure, fertility and meiotic behaviour were assessed in three segregating allohexaploid populations derived from the cross (B. napus × B. carinata) × B. juncea using a combination of molecular marker genotyping, phenotyping and cytogenetics. Plants with unbalanced A-C translocations in one direction (where a C-genome chromosome fragment replaces an A-genome fragment) but not the other (where an A-genome fragment replaces a C-genome fragment) showed significantly reduced fertility across all populations. Genomic regions associated with fertility contained several meiosis genes with putatively causal mutations inherited from the parents (copies of SCC2 in the A genome, PAIR1/PRD3, PRD1 and ATK1/KATA/KIN14a in the B genome, and MSH2 and SMC1/TITAN8 in the C genome). Reduced seed fertility associated with the loss of chromosome fragments from only one subgenome following homoeologous exchanges could comprise a mechanism for biased genome fractionation in allopolyploids. Pre-existing meiosis gene variants present in allotetraploid parents may help to stabilize meiosis in novel allohexaploids.

Drought stress has transgenerational effects on seeds and seedlings in winter oilseed rape (Brassica napus L.).

BACKGROUND: Drought stress has a negative effect on both seed yield and seed quality in Brassica napus (oilseed rape, canola). Here we show that while drought impairs the maternal plant performance, it also increases the vigour of progeny of stressed maternal plants. We investigated the transgenerational influence of abiotic stress by detailed analysis of yield, seed quality, and seedling performance on a growth-related and metabolic level. Seeds of eight diverse winter oilseed rape genotypes were generated under well-watered and drought stress conditions under controlled-environment conditions in large plant containers. RESULTS: We found a decrease in seed quality in seeds derived from mother plants that were exposed to drought stress. At the same time, the seeds that developed under stress conditions showed higher seedling vigour compared to non-stressed controls.This effect on seed quality and seedling vigour was found to be independent of maternal plant yield performance. CONCLUSIONS: Drought stress has a positive transgenerational effect on seedling vigour. Three potential causes for stress-induced improvement of seedling vigour are discussed: (1) Heterotic effects caused by a tendency towards a higher outcrossing rate in response to stress; (2) an altered reservoir of seed storage metabolites to which the seedling resorts during early growth, and (3) inter-generational stress memory, formed by stress-induced changes in the epigenome of the seedling.

Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus.

Homoeologous exchanges (HEs) have been shown to generate novel gene combinations and phenotypes in a range of polyploid species. Gene presence/absence variation (PAV) is also a major contributor to genetic diversity. In this study, we show that there is an association between these two events, particularly in recent Brassica napus synthetic accessions, and that these represent a novel source of genetic diversity, which can be captured for the improvement of this important crop species. By assembling the pangenome of B. napus, we show that 38% of the genes display PAV behaviour, with some of these variable genes predicted to be involved in important agronomic traits including flowering time, disease resistance, acyl lipid metabolism and glucosinolate metabolism. This study is a first and provides a detailed characterization of the association between HEs and PAVs in B. napus at the pangenome level.

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.

Mapping of homoeologous chromosome exchanges influencing quantitative trait variation in Brassica napus.

Genomic rearrangements arising during polyploidization are an important source of genetic and phenotypic variation in the recent allopolyploid crop Brassica napus. Exchanges among homoeologous chromosomes, due to interhomoeologue pairing, and deletions without compensating homoeologous duplications are observed in both natural B. napus and synthetic B. napus. Rearrangements of large or small chromosome segments induce gene copy number variation (CNV) and can potentially cause phenotypic changes. Unfortunately, complex genome restructuring is difficult to deal with in linkage mapping studies. Here, we demonstrate how high-density genetic mapping with codominant, physically anchored SNP markers can detect segmental homoeologous exchanges (HE) as well as deletions and accurately link these to QTL. We validated rearrangements detected in genetic mapping data by whole-genome resequencing of parental lines along with cytogenetic analysis using fluorescence in situ hybridization with bacterial artificial chromosome probes (BAC-FISH) coupled with PCR using primers specific to the rearranged region. Using a well-known QTL region influencing seed quality traits as an example, we confirmed that HE underlies the trait variation in a DH population involving a synthetic B. napus trait donor, and succeeded in narrowing the QTL to a small defined interval that enables delineation of key candidate genes.

Growth factor receptor-bound protein 2 contributes to (hem)immunoreceptor tyrosine-based activation motif-mediated signaling in platelets.

RATIONALE: Platelets are anuclear cell fragments derived from bone marrow megakaryocytes (MKs) that safeguard vascular integrity but may also cause pathological vessel occlusion. One major pathway of platelet activation is triggered by 2 receptors that signal through an (hem)immunoreceptor tyrosine-based activation motif (ITAM), the activating collagen receptor glycoprotein (GP) VI and the C-type lectin-like receptor 2 (CLEC-2). Growth factor receptor-bound protein 2 (Grb2) is a ubiquitously expressed adapter molecule involved in signaling processes of numerous receptors in different cell types, but its function in platelets and MKs is unknown. OBJECTIVE: We tested the hypothesis that Grb2 is a crucial adapter protein in (hem)immunoreceptor tyrosine-based activation motif signaling in platelets. METHODS AND RESULTS: Here, we show that genetic ablation of Grb2 in MKs and platelets did not interfere with MK differentiation or platelet production. However, Grb2-deficiency severely impaired glycoprotein VI-mediated platelet activation because of defective stabilization of the linker of activated T-cell (LAT) signalosome and activation of downstream signaling proteins that resulted in reduced adhesion, aggregation, and coagulant activity on collagen in vitro. Similarly, CLEC-2-mediated signaling was impaired in Grb2-deficient platelets, whereas the cells responded normally to stimulation of G protein-coupled receptors. In vivo, this selective (hem)immunoreceptor tyrosine-based activation motif signaling defect resulted in prolonged bleeding times but affected arterial thrombus formation only after concomitant treatment with acetylsalicylic acid, indicating that defective glycoprotein VI signaling in the absence of Grb2 can be compensated through thromboxane A2-induced G protein-coupled receptor signaling pathways. CONCLUSIONS: These results reveal an important contribution of Grb2 in (hem)immunoreceptor tyrosine-based activation motif signaling in platelets in hemostasis and thrombosis by stabilizing the LAT signalosome.

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.

Characterizing root response phenotypes by neural network analysis.

Roots play an immediate role as the interface for water acquisition. To improve sustainability in low-water environments, breeders of major crops must therefore pay closer attention to advantageous root phenotypes; however, the complexity of root architecture in response to stress can be difficult to quantify. Here, the Sholl method, an established technique from neurobiology used for the characterization of neural network anatomy, was adapted to more adequately describe root responses to osmotic stress. This method was used to investigate the influence of in vitro osmotic stress on early root architecture and distribution in drought-resistant and -susceptible genotypes of winter oilseed rape. Interactive changes in root architecture can be easily captured by individual intersection profiles generated by Sholl analysis. Validation using manual measurements confirmed that the number of lateral roots decreased, while mean lateral root length was enhanced, under osmotic stress conditions. Both genotypes reacted to osmotic stress with a shift in their intersection patterns measured with Sholl analysis. Changes in interactive root architecture and distribution under stress were more pronounced in the drought-resistant genotype, indicating that these changes may contribute to drought resistance under mild osmotic stress conditions. The Sholl methodology is presented as a promising tool for selection of cultivars with advantageous root phenotypes under osmotic stress conditions.

Genetic factors inherited from both diploid parents interact to affect genome stability and fertility in resynthesized allotetraploid Brassica napus.

Established allopolyploids are known to be genomically stable and fertile. However, in contrast, most newly resynthesized allopolyploids are infertile and meiotically unstable. Identifying the genetic factors responsible for genome stability in newly formed allopolyploid is key to understanding how 2 genomes come together to form a species. One hypothesis is that established allopolyploids may have inherited specific alleles from their diploid progenitors which conferred meiotic stability. Resynthesized Brassica napus lines are often unstable and infertile, unlike B. napus cultivars. We tested this hypothesis by characterizing 41 resynthesized B. napus lines produced by crosses between 8 Brassica rapa and 8 Brassica oleracea lines for copy number variation resulting from nonhomologous recombination events and fertility. We resequenced 8 B. rapa and 5 B. oleracea parent accessions and analyzed 19 resynthesized lines for allelic variation in a list of meiosis gene homologs. SNP genotyping was performed using the Illumina Infinium Brassica 60K array for 3 individuals per line. Self-pollinated seed set and genome stability (number of copy number variants) were significantly affected by the interaction between both B. rapa and B. oleracea parental genotypes. We identified 13 putative meiosis gene candidates which were significantly associated with frequency of copy number variants and which contained putatively harmful mutations in meiosis gene haplotypes for further investigation. Our results support the hypothesis that allelic variants inherited from parental genotypes affect genome stability and fertility in resynthesized rapeseed.

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.

Genomic insights into the origin, domestication and diversification of Brassica juncea.

Despite early domestication around 3000 BC, the evolutionary history of the ancient allotetraploid species Brassica juncea (L.) Czern & Coss remains uncertain. Here, we report a chromosome-scale de novo assembly of a yellow-seeded B. juncea genome by integrating long-read and short-read sequencing, optical mapping and Hi-C technologies. Nuclear and organelle phylogenies of 480 accessions worldwide supported that B. juncea is most likely a single origin in West Asia, 8,000-14,000 years ago, via natural interspecific hybridization. Subsequently, new crop types evolved through spontaneous gene mutations and introgressions along three independent routes of eastward expansion. Selective sweeps, genome-wide trait associations and tissue-specific RNA-sequencing analysis shed light on the domestication history of flowering time and seed weight, and on human selection for morphological diversification in this versatile species. Our data provide a comprehensive insight into the origin and domestication and a foundation for genomics-based breeding of B. juncea.

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.

Transcriptomics reveal high regulatory diversity of drought tolerance strategies in a biennial oil crop.

Spring droughts are expected to become more frequent in Central Europe as a result of climate change. Their coincidence with flowering of biennial crops like winter oilseed rape (Brassica napus) can cause major impact for yield development. However, no data is available on the diversity of genetic regulation of drought tolerance during this stage under realistic conditions. Here, we assessed the phenotypic plasticity of drought response for eight diverse B. napus accessions under field-like conditions and linked their stress response to gene and miRNA expression during early and late stress. We observed highly diverse responses, both on the phenotypic and on the gene expression level. Our data suggest that drought tolerant accessions have more effective molecular protection mechanisms like ROS scavenging, source/sink ratio and regulation of developmental timing, compared to otherwise phenotypically similar accessions. Bna.MAP3K13.C05 expression was found to be protective independently of the tolerance mechanism, indicating cross-talk to nitrogen signaling. Moreover, we identified putative miRNA genes in the B. napus genome which respond to stress and may also be involved in protective mechanisms, representing possible breeding targets.

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.

Flowering Time Gene Variation in Brassica Species Shows Evolutionary Principles.

Flowering time genes have a strong influence on successful reproduction and life cycle adaptation. However, their regulation is highly complex and only well understood in diploid model systems. For crops with a polyploid background from the genus Brassica, data on flowering time gene variation are scarce, although indispensable for modern breeding techniques like marker-assisted breeding. We have deep-sequenced all paralogs of 35 Arabidopsis thaliana flowering regulators using Sequence Capture followed by Illumina sequencing in two selected accessions of the vegetable species Brassica rapa and Brassica oleracea, respectively. Using these data, we were able to call SNPs, InDels and copy number variations (CNVs) for genes from the total flowering time network including central flowering regulators, but also genes from the vernalisation pathway, the photoperiod pathway, temperature regulation, the circadian clock and the downstream effectors. Comparing the results to a complementary data set from the allotetraploid species Brassica napus, we detected rearrangements in B. napus which probably occurred early after the allopolyploidisation event. Those data are both a valuable resource for flowering time research in those vegetable species, as well as a contribution to speciation genetics.