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

Identification and characterization of the plasma membrane H+-ATPase genes in Brassica napus and functional analysis of BnHA9 in salt tolerance

As a primary proton pump, plasma membrane (PM) H+-ATPase plays critical roles in regulating plant growth, development, and stress responses. PM H+-ATPases have been well characterized in many plant species. However, no comprehensive study of PM H+-ATPase genes has been performed in Brassica napus (rapeseed). In this study, we identified 32 PM H+-ATPase genes (BnHAs) in the rapeseed genome, and they were distributed on 16 chromosomes. Phylogenetical and gene duplication analyses showed that the BnHA genes were classified into five subfamilies, and the segmental duplication mainly contributed to the expansion of the rapeseed PM H+-ATPase gene family. The conserved domain and subcellular analyses indicated that BnHAs encoded canonical PM H+-ATPase proteins with 14 highly conserved domains and localized on PM. Cis-acting regulatory element and expression pattern analyses indicated that the expression of BnHAs possessed tissue developmental stage specificity. The 25 upstream open reading frames with the canonical initiation codon ATG were predicted in the 5' untranslated regions of 11 BnHA genes and could be used as potential target sites for improving rapeseed traits. Protein interaction analysis showed that BnBRI1.c associated with BnHA2 and BnHA17, indicating that the conserved activity regulation mechanism of BnHAs may be present in rapeseed. BnHA9 overexpression in Arabidopsis enhanced the salt tolerance of the transgenic plants. Thus, our results lay a foundation for further research exploring the biological functions of PM H+-ATPases in rapeseed.

WRKY22 Transcription Factor from Iris laevigata Regulates Flowering Time and Resistance to Salt and Drought

Iris laevigata Fisch. is an excellent ornamental plant in cold regions due to its unique ornamental ability and strong cold resistance. However, the flowering period of the population is only about 20 days, greatly limiting its potential uses in landscaping and the cutting flower industry. In addition, I. laevigata is often challenged with various abiotic stresses including high salinity and drought in its native habitats. Thus, breeding novel cultivars with delayed flowering time and higher resistance to abiotic stress is of high importance. In this study, we utilized sequencing data from the I. laevigata transcriptome to identify WRKYs and characterized IlWRKY22, a key transcription factor that modulates flowering time and abiotic stress responses. IlWRKY22 is induced by salt and drought stress. We cloned IlWRKY22 and found that it is a Group IIe WRKY localized in the nucleus. Overexpressing IlWRKY22 in Arabidopsis thaliana (L.) Heynh. and Nicotiana tabacum L. resulted in a delayed flowering time in the transgenic plants. We created transgenic N. tabacum overexpressing IlWRKY22, which showed significantly improved resistance to both salt and drought compared to the control plants. Thus, our study revealed a unique dual function of IlWRKY22, an excellent candidate gene for breeding novel Iris cultivars of desirable traits.

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.

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

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

Transcriptome Analysis Reveals Key Genes and Pathways Associated with the Regulation of Flowering Time in Cabbage (Brassica oleracea L. var. capitata)

Flowering time is an important agronomic trait in cabbage (Brassica oleracea L. var. capitata), but the molecular regulatory mechanism underlying flowering time regulation in cabbage remains unclear. In this study, transcriptome analysis was performed using two sets of cabbage materials: (1) the early-flowering inbred line C491 (P1) and late-flowering inbred line B602 (P2), (2) the early-flowering individuals F2-B and late-flowering individuals F2-NB from the F2 population. The analysis revealed 9508 differentially expressed genes (DEGs) common to both C491_VS_ B602 and F2-B_VS_F2-NB. The Kyoto Encyclopedia of Genes and Genomes (KEGGs) analysis showed that plant hormone signal transduction and the MAPK signaling pathway were mainly enriched in up-regulated genes, and ribosome and DNA replication were mainly enriched in down-regulated genes. We identified 321 homologues of Arabidopsis flowering time genes (Ft) in cabbage. Among them, 25 DEGs (11 up-regulated and 14 down-regulated genes) were detected in the two comparison groups, and 12 gene expression patterns closely corresponded with the different flowering times in the two sets of materials. Two genes encoding MADS-box proteins, Bo1g157450 (BoSEP2-1) and Bo5g152700 (BoSEP2-2), showed significantly reduced expression in the late-flowering parent B602 compared with the early-flowering parent C491 via qRT-PCR analysis, which was consistent with the RNA-seq data. Next, the expression levels of Bo1g157450 (BoSEP2-1) and Bo5g152700 (BoSEP2-2) were analyzed in two other groups of early-flowering and late-flowering inbred lines, which showed that their expression patterns were consistent with those in the parents. Sequence analysis revealed that three and one SNPs between B602 and C491 were identified in Bo1g157450 (BoSEP2-1) and Bo5g152700 (BoSEP2-2), respectively. Therefore, BoSEP2-1 and BoSEP2-2 were designated as candidates for flowering time regulation through a potential new regulatory pathway. These results provide new insights into the molecular mechanisms underlying flowering time regulation in cabbage.

Brassica napus Haploid and Double Haploid Production and Its Latest Applications

Rapeseed is one of the most important oil crops in the world. Increasing demand for oil and limited agronomic capabilities of present-day rapeseed result in the need for rapid development of new, superior cultivars. Double haploid (DH) technology is a fast and convenient approach in plant breeding as well as genetic research. Brassica napus is considered a model species for DH production based on microspore embryogenesis; however, the molecular mechanisms underlying microspore reprogramming are still vague. It is known that morphological changes are accompanied by gene and protein expression patterns, alongside carbohydrate and lipid metabolism. Novel, more efficient methods for DH rapeseed production have been reported. This review covers new findings and advances in Brassica napus DH production as well as the latest reports related to agronomically important traits in molecular studies employing the double haploid rapeseed lines.

Comparative transcriptome profiling reveals the multiple levels of crosstalk in phytohormone networks in Brassica napus

Plant hormones are the intrinsic factors that control plant development. The integration of different phytohormone pathways in a complex network of synergistic, antagonistic and additive interactions has been elucidated in model plants. However, the systemic level of transcriptional responses to hormone crosstalk in Brassica napus is largely unknown. Here, we present an in-depth temporal-resolution study of the transcriptomes of the seven hormones in B. napus seedlings. Differentially expressed gene analysis revealed few common target genes that co-regulated (up- and down-regulated) by seven hormones; instead, different hormones appear to regulate distinct members of protein families. We then constructed the regulatory networks between the seven hormones side by side, which allowed us to identify key genes and transcription factors that regulate the hormone crosstalk in B. napus. Using this dataset, we uncovered a novel crosstalk between gibberellin and cytokinin in which cytokinin homeostasis was mediated by RGA-related CKXs expression. Moreover, the modulation of gibberellin metabolism by the identified key transcription factors was confirmed in B. napus. Furthermore, all data were available online from http://yanglab.hzau.edu.cn/BnTIR/hormone. Our study reveals an integrated hormone crosstalk network in Brassica napus, which also provides a versatile resource for future hormone studies in plant species.

Identification, evolution and expression analyses of the whole genome-wide PEBP gene family in Brassica napus L

BACKGROUND

With the release of genomic data for B.rapa, B.oleracea, and B.napus, research on the genetic and molecular functions of Brassica spp. has entered a new stage. PEBP genes in plants play an important role in the transition to flowering as well as seed development and germination. Molecular evolutionary and functional analyses of the PEBP gene family in B.napus based on molecular biology methods can provide a theoretical basis for subsequent investigations of related regulators.

RESULTS

In this paper, we identified a total of 29 PEBP genes from B.napus that were located on 14 chromosomes and 3 random locations. Most members contained 4 exons and 3 introns; motif 1 and motif 2 were the characteristic motifs of PEBP members. On the basis of intraspecific and interspecific collinearity analyses, it is speculated that fragment replication and genomic replication are the main drivers of for the amplification and evolution of the PEBP gene in the B.napus genome. The results of promoter cis-elements prediction suggest that BnPEBP family genes are inducible promoters, which may directly or indirectly participate in multiple regulatory pathways of plant growth cycle. Furthermore, the tissue-specific expression results show that the expression levels of BnPEBP family genes in different tissues were quite different, but the gene expression organization and patterns of the same subgroup were basically the same. qRT‒PCR revealed certain spatiotemporal patterns in the expression of the PEBP subgroups in roots, stems, leaves, buds, and siliques, was tissue-specific, and related to function.

CONCLUSIONS

A systematic comparative analysis of the B.napus PEBP gene family was carried out at here. The results of gene identification, phylogenetic tree construction, structural analysis, gene duplication analysis, prediction of promoter cis-elements and interacting proteins, and expression analysis provide a reference for exploring the molecular mechanisms of BnPEBP family genes in future research.

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.

Investigation of Brassica and its relative genomes in the post-genomics era

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae 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.

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.

Transcriptome Characterization and Identification of Molecular Markers (SNP, SSR, and Indels) in the Medicinal Plant Sarcandra glabra spp

Sarcandra glabra has significant metabolically active bioingredients of pharmaceutical importance. The deficiency of molecular markers for S. glabra is a hindrance in molecular breeding for genetic improvement. In this study, 57.756 million pair-end reads were generated by transcriptome sequencing in S. glabra (Thunb.) Nakai and its subspecies S. glabra ssp. brachystachys. A total of 141,954 unigenes with 646.63 bp average length were assembled. A total of 25,620 simple sequence repeats, 726,476 single nucleotide polymorphisms, and 42,939 insertions and deletions were identified, and the associated unigenes and differentially expressed genes were characterized. This work enhanced the molecular marker resources and will facilitate molecular breeding and gene mining in S. glabra spp.

Germinating seed can sense low temperature for the floral transition and vernalization of winter rapeseed (Brassica rapa)

The hybrid production of winter rapeseed is limited by the difficult vernalization processes. Thus, floral regulation of winter rapeseed parental lines cannot be executed through selection of sowing time during hybrid production. Therefore, in this study, strong winter rapeseed was used as the material to analyse the floral transition mechanism of germinating seed vernalization. Results demonstrated that germinating seeds could sense low temperatures and complete vernalization following a low temperature treatment for 56.5 d with a 100 % vernalization rate. The regression equation between vernalization rate (y) and vernalization treatment days (x) was determined as y = 0.019x - 0.0765 (R² = 0.8529). When the vernalization treatment time was prolonged, the vernalization rate and fruiting ability increased rapidly, and variations were observed in the membrane lipid oxidation and physiological characteristics. Furthermore, at the prolonged treatment time of 10-50 d, the salicylic acid (SA) content continued to decrease, with values significantly lower than those of the control. SA content is significantly positively correlated with the level of BrFLC transcription and a significantly negatively correlated with the vernalization rate of germinating seeds. Moreover, the expressions of genes associated with SA biosynthesis, SA signal transduction, the flowering key negative regulators were suppressed and that of positive regulators were promoted during vernalization. These results suggest that SA as a floral repressor is involved in the regulation of the vernalization process of winter rapeseed germination seeds. In addition, SA may be related to the counting dosage of vernalization.

Nitrogen Signaling Genes and SOC1 Determine the Flowering Time in a Reciprocal Negative Feedback Loop in Chinese Cabbage (Brassica rapa L.) Based on CRISPR/Cas9-Mediated Mutagenesis of Multiple BrSOC1 Homologs

Precise flowering timing is critical for the plant life cycle. Here, we examined the molecular mechanisms and regulatory network associated with flowering in Chinese cabbage (Brassica rapa L.) by comparative transcriptome profiling of two Chinese cabbage inbred lines, “4004” (early bolting) and “50” (late bolting). RNA-Seq and quantitative reverse transcription PCR (qPCR) analyses showed that two positive nitric oxide (NO) signaling regulator genes, nitrite reductase (BrNIR) and nitrate reductase (BrNIA), were up-regulated in line “50” with or without vernalization. In agreement with the transcription analysis, the shoots in line “50” had substantially higher nitrogen levels than those in “4004”. Upon vernalization, the flowering repressor gene Circadian 1 (BrCIR1) was significantly up-regulated in line “50”, whereas the flowering enhancer genes named SUPPRESSOR OF OVEREXPRESSION OF CONSTANCE 1 homologs (BrSOC1s) were substantially up-regulated in line “4004”. CRISPR/Cas9-mediated mutagenesis in Chinese cabbage demonstrated that the BrSOC1-1/1-2/1-3 genes were involved in late flowering, and their expression was mutually exclusive with that of the nitrogen signaling genes. Thus, we identified two flowering mechanisms in Chinese cabbage: a reciprocal negative feedback loop between nitrogen signaling genes (BrNIA1 and BrNIR1) and BrSOC1s to control flowering time and positive feedback control of the expression of BrSOC1s.

The transition to flowering in winter rapeseed during vernalization

Flowering time is a major determinant of adaptation, fitness and yield in the allopolyploid species rapeseed (Brassica napus). Despite being a close relative to Arabidopsis thaliana, little is known about the timing of floral transition and the genes that govern this process. Winter, semi-winter and spring type plants have important life history characteristics that differ in vernalization requirements for flowering and are important for growing rapeseed in different regions of the world. In this study, we investigated the timing of vernalization-driven floral transition in winter rapeseed and the effect of photoperiod and developmental age on flowering time and vernalization responsiveness. Microscopy and whole transcriptome analyses at the shoot apical meristems of plants grown under controlled conditions showed that floral transition is initiated within few weeks of vernalization. Certain Bna.SOC1 and Bna.SPL5 homeologs were among the induced genes, suggesting that they are regulating the timing of cold-induced floral transition. Moreover, the flowering response of plants with shorter pre-vernalization period correlated with a delayed expression of Bna.SOC1 and Bna.SPL5 genes. In essence, this study presents a detailed analysis of vernalization-driven floral transition and the aspects of juvenility and dormancy and their effect on flowering time in rapeseed.

Regulation of flowering under short photoperiods based on transcriptomic and metabolomic analysis in Phaseolus vulgaris L

Common bean (Phaseolus vulgaris L.) is a short-day plant and its flowering time, and consequently, pod yield and quality is influenced by photoperiod. In this study, the photoperiodic-sensitive variety 'Hong jin gou', which flowers 31 days (d) earlier in short-day than in long-day, was used as the experimental material. Samples were collected to determine the growth and photosynthetic parameters in each daylength treatment, and transcriptome and metabolome data were conducted. We identified eight genes related to flowering by further screening for differentially expressed genes. These genes function to regulate the biological clock. The combination of differentially expressed genes and metabolites, together with the known regulation network of flowering time and the day-night expression pattern of related genes allow us to speculate on the regulation of flowering time in the common bean and conclude that TIMING OF CAB EXPRESSION1 (TOC1) plays a pivotal role in the network and its upregulation or downregulation causes corresponding changes in the expression of downstream genes. The regulatory network is also influenced by gibberellic acid (GA) and jasmonic acid (JA). These regulatory pathways jointly comprise the flowering regulatory network in common bean.

High-resolution temporal dynamic transcriptome landscape reveals a GhCAL-mediated flowering regulatory pathway in cotton (Gossypium hirsutum L.)

The transition from vegetative to reproductive growth is very important for early maturity in cotton. However, the genetic control of this highly dynamic and complex developmental process remains unclear. A high-resolution tissue- and stage-specific transcriptome profile was generated from six developmental stages using 72 samples of two early-maturing and two late-maturing cotton varieties. The results of histological analysis of paraffin sections showed that flower bud differentiation occurred at the third true leaf stage (3TLS) in early-maturing varieties, but at the fifth true leaf stage (5TLS) in late-maturing varieties. Using pairwise comparison and weighted gene co-expression network analysis, 5312 differentially expressed genes were obtained, which were divided into 10 gene co-expression modules. In the MElightcyan module, 46 candidate genes regulating cotton flower bud differentiation were identified and expressed at the flower bud differentiation stage. A novel key regulatory gene related to flower bud differentiation, GhCAL, was identified in the MElightcyan module. Anti-GhCAL transgenic cotton plants exhibited late flower bud differentiation and flowering time. GhCAL formed heterodimers with GhAP1-A04/GhAGL6-D09 and regulated the expression of GhAP1-A04 and GhAGL6-D09. GhAP1-A04- and GhAGL6-D09-silenced plants also showed significant late flowering. Finally, we propose a new flowering regulatory pathway mediated by GhCAL. This study elucidated the molecular mechanism of cotton flowering regulation and provides good genetic resources for cotton early-maturing breeding.

Gene editing of three BnITPK genes in tetraploid oilseed rape leads to significant reduction of phytic acid in seeds

Commercialization of Brassica napus. L (oilseed rape) meal as protein diet is gaining more attention due to its well-balanced amino acid and protein contents. Phytic acid (PA) is a major source of phosphorus in plants but is considered as anti-nutritive for monogastric animals including humans due to its adverse effects on essential mineral absorption. The undigested PA causes eutrophication, which potentially threatens aquatic life. PA accounts to 2-5% in mature seeds of oilseed rape and is synthesized by complex pathways involving multiple enzymes. Breeding polyploids for recessive traits is challenging as gene functions are encoded by several paralogs. Gene redundancy often requires to knock out several gene copies to study their underlying effects. Therefore, we adopted CRISPR-Cas9 mutagenesis to knock out three functional paralogs of BnITPK. We obtained low PA mutants with an increase of free phosphorus in the canola grade spring cultivar Haydn. These mutants could mark an important milestone in rapeseed breeding with an increase in protein value and no adverse effects on oil contents.

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.

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

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

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Comparative transcriptome and metabolite profiling of four tissues from Alisma orientale (Sam.) Juzep reveals its inflorescence developmental and medicinal characteristics

Alisma orientale (Sam.) Juzep (A. orientale) is an important medicinal plant in traditional Chinese medicine. In this study, de novo RNA-seq of A. orientale was performed based on the cDNA libraries from four different tissues, roots, leaves, scapes and inflorescences. A total of 41,685 unigenes were assembled, 25,024 unigene functional annotations were obtained by searching against the five public sequence databases, and 3,411 simple sequence repeats in A. orientale were reported for the first time. 15,402 differentially expressed genes were analysed. The morphological characteristics showed that compared to the other tissues, the leaves had more chlorophyll, the scapes had more vascular bundles, and the inflorescences contained more starch granules and protein. In addition, the metabolic profiles of eight kinds of alisols metabolite profiling, which were measured by ultra-Performance liquid chromatography-triple quadrupole-mass spectrometry showed that alisol B 23-acetate and alisol B were the major components of the four tissues at amounts of 0.068~0.350 mg/g and 0.046~0.587 mg/g, respectively. In addition, qRT-PCR validated that farnesyl pyrophosphate synthase and 3-hydroxy-3-methylglutaryl-CoA reductase should be considered the critical candidate genes involved in alisol biosynthesis. These transcriptome and metabolic profiles of A. orientale may help clarify the molecular mechanisms underlying the medicinal characteristics of A. orientale.

A Novel Chimeric Mitochondrial Gene Confers Cytoplasmic Effects on Seed Oil Content in Polyploid Rapeseed (Brassica napus)

Cytoplasmic effects (CEs) have been discovered to influence a diverse array of agronomic traits in crops, and understanding the underlying mechanisms can help accelerate breeding programs. Seed oil content (SOC) is of great agricultural, nutritional, and economic importance. However, the genetic basis of CEs on SOC (CE-SOC) remains enigmatic. In this study, we use an optimized approach to sequence the cytoplasmic (plastid and mitochondrial) genomes of allotetraploid oilseed rape (Brassica napus) cultivars, 51218 and 56366, that bear contrasting CE-SOC. By combining comparative genomics and genome-wide transcriptome analysis, we identify mitochondria-encoded orf188 as a potential CE-SOC determinant gene. Functional analyses in the model system Arabidopsis thaliana and rapeseed demonstrated that orf188 governs CE-SOC and could significantly increase SOC, strikingly, through promoting the yield of ATP. Consistent with this finding, transcriptional profiling with microarray and RNA sequencing revealed that orf188 affects transcriptional reprogramming of mitochondrial energy metabolism to facilitate ATP production. Intriguingly, orf188 is a previously uncharacterized chimeric gene, and the presence of this genetic novelty endows rapeseed with positive CE-SOC. Our results shed light on the molecular basis of CEs on a key quantitative trait in polyploid crops and enrich the theory of maternal control of oil content, providing new scientific guidance for breeding high-oil rapeseed germplasms.

Complex Horticultural Quality Traits in Broccoli Are Illuminated by Evaluation of the Immortal BolTBDH Mapping Population

Improving horticultural quality in regionally adapted broccoli (Brassica oleracea var. italica) and other B. oleracea crops is challenging due to complex genetic control of traits affecting morphology, development, and yield. Mapping horticultural quality traits to genomic loci is an essential step in these improvement efforts. Understanding the mechanisms underlying horticultural quality enables multi-trait marker-assisted selection for improved, resilient, and regionally adapted B. oleracea germplasm. The publicly-available biparental double-haploid BolTBDH mapping population (Chinese kale × broccoli; N = 175) was evaluated for 25 horticultural traits in six trait classes (architecture, biomass, phenology, leaf morphology, floral morphology, and head quality) by multiple quantitative trait loci mapping using 1,881 genotype-by-sequencing derived single nucleotide polymorphisms. The physical locations of 56 single and 41 epistatic quantitative trait locus (QTL) were identified. Four head quality QTL (OQ_C03@57.0, OQ_C04@33.3, OQ_CC08@25.5, and OQ_C09@49.7) explain a cumulative 81.9% of phenotypic variance in the broccoli heading phenotype, contain the FLOWERING LOCUS C (FLC) homologs Bo9g173400 and Bo9g173370, and exhibit epistatic effects. Three key genomic hotspots associated with pleiotropic control of the broccoli heading phenotype were identified. One phenology hotspot reduces days to flowering by 7.0 days and includes an additional FLC homolog Bo3g024250 that does not exhibit epistatic effects with the three horticultural quality hotspots. Strong candidates for other horticultural traits were identified: BoLMI1 (Bo3g002560) associated with serrated leaf margins and leaf apex shape, BoCCD4 (Bo3g158650) implicated in flower color, and BoAP2 (Bo1g004960) implicated in the hooked sepal horticultural trait. The BolTBDH population provides a framework for B. oleracea improvement by targeting key genomic loci contributing to high horticultural quality broccoli and enabling de novo mapping of currently unexplored traits.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Translating Flowering Time From Arabidopsis thaliana to Brassicaceae and Asteraceae Crop Species

Flowering and seed set are essential for plant species to survive, hence plants need to adapt to highly variable environments to flower in the most favorable conditions. Endogenous cues such as plant age and hormones coordinate with the environmental cues like temperature and day length to determine optimal time for the transition from vegetative to reproductive growth. In a breeding context, controlling flowering time would help to speed up the production of new hybrids and produce high yield throughout the year. The flowering time genetic network is extensively studied in the plant model species Arabidopsis thaliana, however this knowledge is still limited in most crops. This article reviews evidence of conservation and divergence of flowering time regulation in A. thaliana with its related crop species in the Brassicaceae and with more distant vegetable crops within the Asteraceae family. Despite the overall conservation of most flowering time pathways in these families, many genes controlling this trait remain elusive, and the function of most Arabidopsis homologs in these crops are yet to be determined. However, the knowledge gathered so far in both model and crop species can be already exploited in vegetable crop breeding for flowering time control.

Genome-Wide Identification of Flowering-Time Genes in Brassica Species and Reveals a Correlation between Selective Pressure and Expression Patterns of Vernalization-Pathway Genes in Brassica napus

Flowering time is a key agronomic trait, directly influencing crop yield and quality. Many flowering-time genes have been identified and characterized in the model plant Arabidopsis thaliana; however, these genes remain uncharacterized in many agronomically important Brassica crops. In this study, we identified 1064, 510, and 524 putative orthologs of A. thaliana flowering-time genes from Brassica napus, Brassica rapa, and Brassica oleracea, respectively, and found that genes involved in the aging and ambient temperature pathways were fewer than those in other flowering pathways. Flowering-time genes were distributed mostly on chromosome C03 in B. napus and B. oleracea, and on chromosome A09 in B. rapa. Calculation of non-synonymous (Ka)/synonymous substitution (Ks) ratios suggested that flowering-time genes in vernalization pathways experienced higher selection pressure than those in other pathways. Expression analysis showed that most vernalization-pathway genes were expressed in flowering organs. Approximately 40% of these genes were highly expressed in the anther, whereas flowering-time integrator genes were expressed in a highly organ-specific manner. Evolutionary selection pressures were negatively correlated with the breadth and expression levels of vernalization-pathway genes. These findings provide an integrated framework of flowering-time genes in these three Brassica crops and provide a foundation for deciphering the relationship between gene expression patterns and their evolutionary selection pressures in Brassica napus.

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.