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

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.

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.

Identification of vernalization-related genes and cold memory element (CME) required for vernalization response in radish (Raphanus sativus L.)

Floral transition is accelerated by exposure to long-term cold like winter in plants, which is called as vernalization. Acceleration of floral transition by vernalization is observed in a diversity of biennial and perennial plants including Brassicaceae family plants. Scientific efforts to understand molecular mechanism underlying vernalization-mediated floral transition have been intensively focused in model plant Arabidopsis thaliana. To get a better understanding on floral transition by vernalization in radish (Raphanus sativus L.), we investigated transcriptomic changes taking place during vernalization in radish. Thousands of genes were differentially regulated along time course of vernalization compared to non-vernalization (NV) sample. Twelve major clusters of DEGs were identified based on distinctive expression profiles during vernalization. Radish FLC homologs were shown to exert an inhibition of floral transition when transformed into Arabidopsis plants. In addition, DNA region containing RY motifs located within a Raphanus sativus FLC homolog, RsFLC1 was found to be required for repression of RsFLC1 by vernalization. Transgenic plants harboring disrupted RY motifs were impaired in the enrichment of H3K27me3 on RsFLC1 chromatin, thus resulting in the delayed flowering in Arabidopsis. Taken together, we report transcriptomic profiles of radish during vernalization and demonstrate the requirement of RY motif for vernalization-mediated repression of RsFLC homologs in radish (Raphanus sativus L.).

Genome-wide characterization and expression analysis of MADS-box transcription factor gene family in Perilla frutescens

MADS-box transcription factors are widely involved in the regulation of plant growth, developmental processes, and response to abiotic stresses. Perilla frutescens, a versatile plant, is not only used for food and medicine but also serves as an economical oil crop. However, the MADS-box transcription factor family in P. frutescens is still largely unexplored. In this study, a total of 93 PfMADS genes were identified in P. frutescens genome. These genes, including 37 Type I and 56 Type II members, were randomly distributed across 20 chromosomes and 2 scaffold regions. Type II PfMADS proteins were found to contain a greater number of motifs, indicating more complex structures and diverse functions. Expression analysis revealed that most PfMADS genes (more than 76 members) exhibited widely expression model in almost all tissues. The further analysis indicated that there was strong correlation between some MIKCC-type PfMADS genes and key genes involved in lipid synthesis and flavonoid metabolism, which implied that these PfMADS genes might play important regulatory role in the above two pathways. It was further verified that PfMADS47 can effectively mediate the regulation of lipid synthesis in Chlamydomonas reinhardtii transformants. Using cis-acting element analysis and qRT-PCR technology, the potential functions of six MIKCC-type PfMADS genes in response to abiotic stresses, especially cold and drought, were studied. Altogether, this study is the first genome-wide analysis of PfMADS. This result further supports functional and evolutionary studies of PfMADS gene family and serves as a benchmark for related P. frutescens breeding studies.

Characterization of FLOWERING LOCUS C 5 in Brassica rapa L

Brassica rapa L., which includes Chinese cabbage, turnip, and pak choi, has more complex flowering time regulation than does Arabidopsis thaliana due to the presence of multiple paralogous flowering time genes. FLOWERING LOCUS C (FLC) is one of the key genes regulating the flowering time, and B. rapa has four FLC paralogs. BrFLC5 on the reference genome is deemed a pseudogene because of a mutation (from G to A) in the splice site of the third intron, but there are some accessions with a G nucleotide in the splice site. In this study, we genotyped 310 B. rapa accessions and found that 19 had homozygous and 81 had heterozygous putative functional BrFLC5 alleles. Accessions of turnip showed the highest proportion with a functional BrFLC5 allele. BrFLC5 acts as a floral repressor when overexpressed in A. thaliana. The BrFLC5 expression level varied in pre-vernalized plants, and this transcriptional variation was not associated with the G/A polymorphism in the third intron. Three accessions having a higher BrFLC5 expression in pre-vernalized plants had a 584-bp insertion in the promoter region. Many regions homologous to this 584-bp sequence are present in the B. rapa genome, and this 584-bp inserted region has tandem duplications of an AT-rich sequence in its central region. The possibility that a high expression of a functional BrFLC5 could contribute to producing premature bolting-resistant lines in B. rapa vegetables is discussed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01405-0.

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.

Genome-wide identification and functional analysis of Dof transcription factor family in Camelina sativa

BACKGROUND

Dof transcription factors (TFs) containing C₂-C₂ zinc finger domains are plant-specific regulatory proteins, playing crucial roles in a variety of biological processes. However, little is known about Dof in Camelina sativa, an important oil crop worldwide, with high stress tolerance. In this study, a genome-wide characterization of Dof proteins is performed to examine their basic structural characteristics, phylogenetics, expression patterns, and functions to identify the regulatory mechanism underlying lipid/oil accumulation and the candidate Dofs mediating stress resistance regulation in C. sativa.

RESULTS

Total of 103 CsDof genes unevenly distributed on 20 chromosomes were identified from the C. sativa genome, and they were classified into four groups (A, B, C and D) based on the classification of Arabidopsis Dof gene family. All of the CsDof proteins contained the highly-conserved typic CX₂C-X₂₁-CX₂C structure. Segmental duplication and purifying selection were detected for CsDof genes. 61 CsDof genes were expressed in multiple tissues, and 20 of them showed tissue-specific expression patterns, suggesting that CsDof genes functioned differentially in different tissues of C. sativa. Remarkably, a set of CsDof members were detected to be possible involved in regulation of oil/lipid biosynthesis in C. sativa. Six CsDof genes exhibited significant expression changes in seedlings under salt stress treatment.

CONCLUSIONS

The present data reveals that segmental duplication is the key force responsible for the expansion of CsDof gene family, and a strong purifying pressure plays a crucial role in CsDofs' evolution. Several CsDof TFs may mediate lipid metabolism and stress responses in C. sativa. Several CsDof TFs may mediate lipid metabolism and stress responses in C. sativa. Collectively, our findings provide a foundation for deep understanding the roles of CsDofs and genetic improvements of oil yield and salt stress tolerance in this species and the related crops.

Whole-Genome DNA Methylation Analysis in Brassica rapa subsp. perviridis in Response to Albugo candida Infection

DNA methylation is an epigenetic mark associated with several mechanisms in plants including immunity mechanisms. However, little is known about the regulatory role of DNA methylation in the resistance response of Brassica species against fungal diseases. White rust, caused by the fungus Albugo candida, is one of the most widespread and destructive diseases of all the cultivated Brassica species, particularly Brassica rapa L. and Brassica juncea (L.) Czern and Coss. Here, we investigate whole-genome DNA methylation modifications of B. rapa subsp. perviridis in response to white rust. As a result, 233 and 275 differentially methylated regions (DMRs) in the susceptible cultivar "Misugi" and the resistant cultivar "Nanane" were identified, respectively. In both cultivars, more than half of the DMRs were associated with genes (DMR-genes). Gene expression analysis showed that 13 of these genes were also differentially expressed between control and infected samples. Gene ontology enrichment analysis of DMR genes revealed their involvement in various biological processes including defense mechanisms. DMRs were unevenly distributed around genes in susceptible and resistant cultivars. In "Misugi," DMRs tended to be located within genes, while in "Nanane," DMRs tended to be located up and downstream of the genes. However, CG DMRs were predominantly located within genes in both cultivars. Transposable elements also showed association with all three sequence contexts of DMRs but predominantly with CHG and CHH DMRs in both cultivars. Our findings indicate the occurrence of DNA methylation modifications in B. rapa in response to white rust infection and suggest a potential regulatory role of DNA methylation modification in defense mechanisms which could be exploited to improve disease resistance.

Conserved and distinct roles of H3K27me3 demethylases regulating flowering time in Brassica rapa

Epigenetic regulation is necessary for optimal organism development and preservation of gene expression profiles in the cell. In plants, the trimethylation of histone H3 lysine 27 (H3K27me3) is a silencing epigenetic mark relevant for developmental transitions like flowering. The floral transition is a key agronomic trait; however, the epigenetic mechanisms of flowering time regulation in crops remain poorly understood. Here we study the Jumonji H3K27me3 demethylases BraA.REF6 and BraA.ELF6 in Brassica rapa. Phenotypic characterization of novel mutant lines and genome-wide H3K27me3 chromatin immunoprecipitation and transcriptomic analyses indicated that BraA.REF6 plays a greater role than BraA.ELF6 in fine-tuning H3K27me3 levels. In addition, we found that braA.elf6 mutants were early flowering due to high H3K27me3 levels at B. rapa homologs of the floral repressor FLC. Unlike mutations in Arabidopsis thaliana, braA.ref6 mutants were late flowering without altering the expression of B. rapa FLC genes. Remarkably, we found that BraA.REF6 regulated a number of gibberellic acid (GA) biosynthetic genes, including a homolog of GA1, and that GA-treatment complemented the late flowering mutant phenotype. This study increases our understanding of the epigenetic regulation of flowering time in B. rapa, highlighting conserved and distinct regulatory mechanisms between model and crop species.

Transcriptome analysis reveals anthocyanin regulation in Chinese cabbage (Brassica rapa L.) at low temperatures

Chinese cabbage that prefers cold conditions is also affected by low-temperature stress, such as the accumulation of leaf anthocyanins. Research on anthocyanin biosynthesis and regulation mechanisms has made great progress. However, research on anthocyanin accumulation for resistance to biological and non-biological stress is still lacking. To study the relationship between anthocyanin accumulation of Chinese cabbage and resistance under low-temperature conditions, RNA sequencing (RNA-seq) was performed on Chinese cabbage ‘Xiao Baojian’ grown at a low temperature for four time periods and at a control temperature for five time periods. In Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, 7954 differentially expressed genes (DEGs) were enriched, of which 587 DEGs belonged to "biosynthesis of other secondary metabolites." Gene temporal expression patterns were used to discover enriched genes related to phenylpropanoid biosynthesis; flavonoid biosynthesis and anthocyanin biosynthesis pathways were found in cluster 1. The interaction networks were constructed, and hub genes were selected, showing that flavonoid biosynthesis pathway genes (DFR, ANS, F3H, FLS1, CHS1, CHS3, and TT8) and defense mechanisms-related genes (DFR, SNL6, and TKPR1) interact with each other. Anthocyanin biosynthesis DEGs in Chinese cabbage were evaluated under low-temperature conditions to map the relevant pathways, and expression maps of transcription factors in the flavonoid pathway were created at various periods. Low temperature upregulated the expression of genes related to anthocyanin biosynthesis. Taken together, our results provide further analysis of the relationship between plant anthocyanin synthesis and stress resistance and may also provide further insights for the future development of high-quality color and cold-tolerant Chinese cabbage germplasm resources.

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

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

Development of a New DNA Marker for Fusarium Yellows Resistance in Brassica rapa Vegetables

In vegetables of Brassica rapa L., Fusarium oxysporum f. sp. rapae (For) or F. oxysporum f. sp. conglutinans (Foc) cause Fusarium yellows. A resistance gene against Foc (FocBr1) has been identified, and deletion of this gene results in susceptibility (focbr1-1). In contrast, a resistance gene against For has not been identified. Inoculation tests showed that lines resistant to Foc were also resistant to For, and lines susceptible to Foc were susceptible to For. However, prediction of disease resistance by a dominant DNA marker on FocBr1 (Bra012688m) was not associated with disease resistance of For in some komatsuna lines using an inoculation test. QTL-seq using four F₂ populations derived from For susceptible and resistant lines showed one causative locus on chromosome A03, which covers FocBr1. Comparison of the amino acid sequence of FocBr1 between susceptible and resistant alleles (FocBr1 and FocBo1) showed that six amino acid differences were specific to susceptible lines. The presence and absence of FocBr1 is consistent with For resistance in F₂ populations. These results indicate that FocBr1 is essential for For resistance, and changed amino acid sequences result in susceptibility to For. This susceptible allele is termed focbr1-2, and a new DNA marker (focbr1-2m) for detection of the focbr1-2 allele was developed.