Regulation and function of SOC1, a flowering pathway integrator

Warm temperature during floral bud transition turns off EjTFL1 gene expression and promotes flowering in Loquat (Eriobotrya japonica Lindl.)

The Rosaceae family includes several deciduous woody species whose flower development extends over two consecutive growing seasons with a winter dormant period in between. Loquat (Eriobotrya japonica Lindl.) belongs to this family, but it is an evergreen species whose flower bud initiation and flowering occur within the same growing year. Vegetative growth dominates from spring to late summer when terminal buds bloom as panicles. Thus, its floral buds do not undergo winter dormancy until flowering, but a summer heat period of dormancy is required for floral bud differentiation, and that is why we used loquat to study the mechanism by which this summer rest period contributes to floral differentiation of Rosaceae species. As for the deciduous species, the bud transition to the generative stage is initiated by the floral integrator genes. There is evidence that combinations of environmental signals and internal cues (plant hormones) control the expression of TFL1, but the mechanism by which this gene regulates its expression in loquat needs to be clarified for a better understanding of its floral initiation and seasonal growth cycles. Under high temperatures (>25ºC) after floral bud inductive period, EjTFL1 expression decreases during meristem transition to the reproductive stage, and the promoters of flowering (EjAP1 and EjLFY) increase, indicating that the floral bud differentiation is affected by high temperatures. Monitoring the apical meristem of loquat in June-August of two consecutive years under ambient and thermal controlled conditions showed that under lower temperatures (<25ºC) during the same period, shoot apex did not stop growing and a higher EjTFL1 expression was recorded, preventing the bud to flower. Likewise, temperature directly affects ABA content in the meristem paralleling EjTFL1 expression, suggesting signaling cascades could converge to refine the expression of EjTFL1 under specific conditions (Tª<25ºC) during the floral transition stage.

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.

JcSEUSS1 negatively regulates reproductive organ development in perennial woody Jatropha curcas

MAIN CONCLUSION

Overexpression of JcSEUSS1 resulted in late flowering, reduced flower number, wrinkled kernels, and decreased seed yield in Jatopha curcas, while downregulation of JcSEUSS1 increased flower number and seed production. The seed oil of Jatropha curcas is suitable as an ideal alternative for diesel fuel, yet the seed yield of Jatropha is restricted by its small number of female flowers and low seed setting rate. Therefore, it is crucial to identify genes that regulate flowering and seed set, and hence improve seed yield. In this study, overexpression of JcSEUSS1 resulted in late flowering, fewer flowers and fruits, and smaller fruits and seeds, causing reduced seed production and oil content. In contrast, the downregulation of JcSEUSS1 by RNA interference (RNAi) technology caused an increase in the flower number and seed yield. However, the flowering time, seed number per fruit, seed weight, and size exhibited no obvious changes in JcSEUSS1-RNAi plants. Moreover, the fatty acid composition also changed in JcSEUSS1 overexpression and RNAi plants, the percentage of unsaturated fatty acids (FAs) was increased in overexpression plants, and the saturated FAs were increased in RNAi plants. These results indicate that JcSEUSS1 played a negative role in regulating reproductive growth and worked redundantly with other genes in the regulation of flowering time, seed number per fruit, seed weight, and size.

Identification, characterization, and comprehensive expression profiling of floral master regulators in pigeon pea (Cajanus cajan [L.] Millspaugh)

Pigeon pea is an important protein-rich pulse crop. Identification of flowering master regulators in pigeon pea is highly imperative as indeterminacy and late flowering are impediments towards yield improvement. A genome-wide analysis was performed to explore flowering orthologous groups in pigeon pea. Among the 412 floral orthologs identified in pigeon pea, 148 genes belong to the meristem identity, photoperiod-responsive, and circadian clock-associated ortholog groups. Our comparative genomics study revealed purifying selection pressures (ka/ks) on floral orthologs, and duplication patterns and evolution through synteny with other model species. Phylogenetic analysis of floral genes substantiated a connection between pigeon pea plant architecture and flowering time as all the PEBP domain-containing genes belong to meristem identity floral networks of pigeon pea. Expression profiling of eleven major orthologs in contrasting determinate and indeterminate genotypes indicated that these orthologs might be involved in flowering regulation. Expression of floral inducer, FT, and floral repressor, TFL1, was non-comparable in indeterminate genotypes across all the developmental stages of pigeon pea. However, dynamic FT/TFL1 expression ratio detected in all tissues of both the genotypes suggested their role in floral transition. One TFL1 ortholog having high sequence conserveness across pigeon pea genotypes showed differential expression indicating genotype-dependent regulation of this ortholog. Presence of conserved 6mA-methylation patterns in light-responsive elements and in other cis-regulatory elements of FT and TFL1 across different plant genotypes indicated possible involvement of epigenetic regulation in flowering.

Gibberellin signaling modulates flowering via the DELLA-BRAHMA-NF-YC module in Arabidopsis

Gibberellin (GA) plays a key role in floral induction by activating the expression of floral integrator genes in plants, but the epigenetic regulatory mechanisms underlying this process remain unclear. Here, we show that BRAHMA (BRM), a core subunit of the chromatin-remodeling SWItch/sucrose nonfermentable (SWI/SNF) complex that functions in various biological processes by regulating gene expression, is involved in GA-signaling-mediated flowering via the formation of the DELLA-BRM-NF-YC module in Arabidopsis (Arabidopsis thaliana). DELLA, BRM, and NF-YC transcription factors interact with one another, and DELLA proteins promote the physical interaction between BRM and NF-YC proteins. This impairs the binding of NF-YCs to SOC1, a major floral integrator gene, to inhibit flowering. On the other hand, DELLA proteins also facilitate the binding of BRM to SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1). The GA-induced degradation of DELLA proteins disturbs the DELLA-BRM-NF-YC module, prevents BRM from inhibiting NF-YCs, and decreases the DNA-binding ability of BRM, which promote the deposition of H3K4me3 on SOC1 chromatin, leading to early flowering. Collectively, our findings show that BRM is a key epigenetic partner of DELLA proteins during the floral transition. Moreover, they provide molecular insights into how GA signaling coordinates an epigenetic factor with a transcription factor to regulate the expression of a flowering gene and flowering in plants.

Warmer temperature during asexual reproduction induce methylome, transcriptomic, and lasting phenotypic changes in Fragaria vesca ecotypes

Plants must adapt with increasing speed to global warming to maintain their fitness. One rapid adaptation mechanism is epigenetic memory, which may provide organisms sufficient time to adapt to climate change. We studied how the perennial Fragaria vesca adapted to warmer temperatures (28°C vs. 18°C) over three asexual generations. Differences in flowering time, stolon number, and petiole length were induced by warmer temperature in one or more ecotypes after three asexual generations and persisted in a common garden environment. Induced methylome changes differed between the four ecotypes from Norway, Iceland, Italy, and Spain, but shared methylome responses were also identified. Most differentially methylated regions (DMRs) occurred in the CHG context, and most CHG and CHH DMRs were hypermethylated at the warmer temperature. In eight CHG DMR peaks, a highly similar methylation pattern could be observed between ecotypes. On average, 13% of the differentially methylated genes between ecotypes also showed a temperature-induced change in gene expression. We observed ecotype-specific methylation and expression patterns for genes related to gibberellin metabolism, flowering time, and epigenetic mechanisms. Furthermore, we observed a negative correlation with gene expression when repetitive elements were found near (±2 kb) or inside genes. In conclusion, lasting phenotypic changes indicative of an epigenetic memory were induced by warmer temperature and were accompanied by changes in DNA methylation patterns. Both shared methylation patterns and transcriptome differences between F. vesca accessions were observed, indicating that DNA methylation may be involved in both general and ecotype-specific phenotypic variation.

Functional conservation and divergence of SEPALLATA-like genes in floral development in Cymbidium sinense

Cymbidium sinense is one of the most important traditional Chinese Orchids due to its unique and highly ornamental floral organs. Although the ABCDE model for flower development is well-established in model plant species, the precise roles of these genes in C. sinense are not yet fully understood. In this study, four SEPALLATA-like genes were isolated and identified from C. sinense. CsSEP1 and CsSEP3 were grouped into the AGL9 clade, while CsSEP2 and CsSEP4 were included in the AGL2/3/4 clade. The expression pattern of CsSEP genes showed that they were significantly accumulated in reproductive tissues and expressed during flower bud development but only mildly detected or even undetected in vegetative organs. Subcellular localization revealed that CsSEP1 and CsSEP4 were localized to the nucleus, while CsSEP2 and CsSEP3 were located at the nuclear membrane. Promoter sequence analysis predicted that CsSEP genes contained a number of hormone response elements (HREs) and MADS-box binding sites. The early flowering phenotype observed in transgenic Arabidopsis plants expressing four CsSEP genes, along with the expression profiles of endogenous genes, such as SOC1, LFY, AG, FT, SEP3 and TCPs, in both transgenic Arabidopsis and C. sinense protoplasts, suggested that the CsSEP genes played a regulatory role in the flowering transition by influencing downstream genes related to flowering. However, only transgenic plants overexpressing CsSEP3 and CsSEP4 caused abnormal phenotypes of floral organs, while CsSEP1 and CsSEP2 had no effect on floral organs. Protein-protein interaction assays indicated that CsSEPs formed a protein complex with B-class CsAP3-2 and CsSOC1 proteins, affecting downstream genes to regulate floral organs and flowering time. Our findings highlighted both the functional conservation and divergence of SEPALLATA-like genes in C. sinense floral development. These results provided a valuable foundation for future studies of the molecular network underlying floral development in C. sinense.

Genome-wide analysis of the MADS-box gene family in Lonicera japonica and a proposed floral organ identity model

BACKGROUND

Lonicera japonica Thunb. is widely used in traditional Chinese medicine. Medicinal L. japonica mainly consists of dried flower buds and partially opened flowers, thus flowers are an important quality indicator. MADS-box genes encode transcription factors that regulate flower development. However, little is known about these genes in L. japonica.

RESULTS

In this study, 48 MADS-box genes were identified in L. japonica, including 20 Type-I genes (8 Mα, 2 Mβ, and 10 Mγ) and 28 Type-II genes (26 MIKCc and 2 MIKC*). The Type-I and Type-II genes differed significantly in gene structure, conserved domains, protein structure, chromosomal distribution, phylogenesis, and expression pattern. Type-I genes had a simpler gene structure, lacked the K domain, had low protein structure conservation, were tandemly distributed on the chromosomes, had more frequent lineage-specific duplications, and were expressed at low levels. In contrast, Type-II genes had a more complex gene structure; contained conserved M, I, K, and C domains; had highly conserved protein structure; and were expressed at high levels throughout the flowering period. Eleven floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in L. japonica. By integrating expression pattern and protein interaction data for these genes, we developed a possible model for floral organ identity determination.

CONCLUSION

This study genome-widely identified and characterized the MADS-box gene family in L. japonica. Eleven floral homeotic MADS-box genes were identified and a possible model for floral organ identity determination was also developed. This study contributes to our understanding of the MADS-box gene family and its possible involvement in floral organ development in L. japonica.

miRNAs for crop improvement

Climate change significantly impacts crop production by inducing several abiotic and biotic stresses. The increasing world population, and their food and industrial demands require focused efforts to improve crop plants to ensure sustainable food production. Among various modern biotechnological tools, microRNAs (miRNAs) are one of the fascinating tools available for crop improvement. miRNAs belong to a class of small non-coding RNAs playing crucial roles in numerous biological processes. miRNAs regulate gene expression by post-transcriptional target mRNA degradation or by translation repression. Plant miRNAs have essential roles in plant development and various biotic and abiotic stress tolerance. In this review, we provide propelling evidence from previous studies conducted around miRNAs and provide a one-stop review of progress made for breeding stress-smart future crop plants. Specifically, we provide a summary of reported miRNAs and their target genes for improvement of plant growth and development, and abiotic and biotic stress tolerance. We also highlight miRNA-mediated engineering for crop improvement and sequence-based technologies available for the identification of miRNAs associated with stress tolerance and plant developmental events.

A novel MORN-motif type gene GmMRF2 controls flowering time and plant height of soybean

The flowering time of soybean is a highly important agronomic characteristic, which affects the adaptability and yield. AtMRF1, a MORN-repeat motif gene, acts as a floral promoter in Arabidopsis, its functions in soybean are not yet understood. Here, we employed qRT-PCR to analyze the tissue expression patten of MRF1 homologs in soybean and determined that the GmMRF2 gene, containing a MORN-motif, highly expressed in the shoot and responded to photoperiod. GmMRF2 overexpression soybean lines exhibited earlier flowering time under long-day (LD) conditions, and increased plant height under both LD and short-day (SD) conditions compared to wild-type (WT) plants. The expression levels of gibberellic acid (GA) pathway genes that positively regulate plant height genes and flowering-promoting genes were up-regulated in the GmMRF2 overexpression lines, were up-regulated in the GmMRF2 overexpression lines. Further study revealed that GmMRF2 interacted with GmTCP15 to co-induce the expression of GmSOC1b. Together, our results preliminarily reveal the functions and mechanisms of GmMRF2 in regulating flowering time and plant height, provide a new promising gene for soybean crop improvement.

TM3 and STM3 Promote Flowering Together with FUL2 and MBP20, but Act Antagonistically in Inflorescence Branching in Tomato

The moment at which a plant transitions to reproductive development is paramount to its life cycle and is strictly controlled by many genes. The transcription factor SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) plays a central role in this process in Arabidopsis. However, the role of SOC1 in tomato (Solanum lycopersicum) has been sparsely studied. Here, we investigated the function of four tomato SOC1 homologs in the floral transition and inflorescence development. We thoroughly characterized the SOC1-like clade throughout the Solanaceae and selected four tomato homologs that are dynamically expressed upon the floral transition. We show that of these homologs, TOMATO MADS 3 (TM3) and SISTER OF TM3 (STM3) promote the primary and sympodial transition to flowering, while MADS-BOX PROTEIN 23 (MBP23) and MBP18 hardly contribute to flowering initiation in the indeterminate cultivar Moneyberg. Protein-protein interaction assays and whole-transcriptome analysis during reproductive meristem development revealed that TM3 and STM3 interact and share many targets with FRUITFULL (FUL) homologs, including cytokinin regulators. Furthermore, we observed that mutating TM3/STM3 affects inflorescence development, but counteracts the inflorescence-branching phenotype of ful2 mbp20. Collectively, this indicates that TM3/STM3 promote the floral transition together with FUL2/MBP20, while these transcription factors have opposite functions in inflorescence development.

FLOWERING LOCUS C drives delayed flowering in Arabidopsis grown and selected at elevated CO 2

Altered flowering time at elevated [CO ₂ ] is well documented, although mechanisms are not well understood. An Arabidopsis genotype previously selected for high fitness at elevated [CO ₂ ] (SG) showed delayed flowering and larger size at flowering when grown at elevated (700 ppm) versus current (380 ppm) [CO ₂ ]. This response was correlated with prolonged expression of FLOWERING LOCUS C ( FLC ), a vernalization-responsive floral repressor gene. To determine if FLC directly delays flowering at elevated [CO ₂ ] in SG, we used vernalization (extended cold) to downregulate FLC expression. We hypothesized that vernalization would eliminate delayed flowering at elevated [CO ₂ ] through the direct reduction of FLC expression, eliminating differences in flowering time between current and elevated [CO ₂ ]. We found that with downregulation of FLC expression via vernalization, SG plants grown at elevated [CO ₂ ] no longer delayed flowering compared to current [CO ₂ ]. Thus, vernalization returned the earlier flowering phenotype, counteracting effects of elevated [CO ₂ ] on flowering. This study indicates that elevated [CO ₂ ] can delay flowering directly through FLC , and downregulation of FLC under elevated [CO ₂ ] reverses this effect. Moreover, this study demonstrates that increasing [CO ₂ ] may potentially drive major changes in development through FLC .

Molecular Cues for Phenological Events in the Flowering Cycle in Avocado

Reproductively mature horticultural trees undergo an annual flowering cycle that repeats each year of their reproductive life. This annual flowering cycle is critical for horticultural tree productivity. However, the molecular events underlying the regulation of flowering in tropical tree crops such as avocado are not fully understood or documented. In this study, we investigated the potential molecular cues regulating the yearly flowering cycle in avocado for two consecutive crop cycles. Homologues of flowering-related genes were identified and assessed for their expression profiles in various tissues throughout the year. Avocado homologues of known floral genes FT, AP1, LFY, FUL, SPL9, CO and SEP2/AGL4 were upregulated at the typical time of floral induction for avocado trees growing in Queensland, Australia. We suggest these are potential candidate markers for floral initiation in these crops. In addition, DAM and DRM1, which are associated with endodormancy, were downregulated at the time of floral bud break. In this study, a positive correlation between CO activation and FT in avocado leaves to regulate flowering was not seen. Furthermore, the SOC1-SPL4 model described in annual plants appears to be conserved in avocado. Lastly, no correlation of juvenility-related miRNAs miR156, miR172 with any phenological event was observed.

A Rapid Alkalinization Factor-like Peptide EaF82 Impairs Tapetum Degeneration during Pollen Development through Induced ATP Deficiency

In plants, the timely degeneration of tapetal cells is essential for providing nutrients and other substances to support pollen development. Rapid alkalinization factors (RALFs) are small, cysteine-rich peptides known to be involved in various aspects of plant development and growth, as well as defense against biotic and abiotic stresses. However, the functions of most of them remain unknown, while no RALF has been reported to involve tapetum degeneration. In this study, we demonstrated that a novel cysteine-rich peptide, EaF82, isolated from shy-flowering 'Golden Pothos' (Epipremnum aureum) plants, is a RALF-like peptide and displays alkalinizing activity. Its heterologous expression in Arabidopsis delayed tapetum degeneration and reduced pollen production and seed yields. RNAseq, RT-qPCR, and biochemical analyses showed that overexpression of EaF82 downregulated a group of genes involved in pH changes, cell wall modifications, tapetum degeneration, and pollen maturation, as well as seven endogenous Arabidopsis RALF genes, and decreased proteasome activity and ATP levels. Yeast two-hybrid screening identified AKIN10, a subunit of energy-sensing SnRK1 kinase, as its interacting partner. Our study reveals a possible regulatory role for RALF peptide in tapetum degeneration and suggests that EaF82 action may be mediated through AKIN10 leading to the alteration of transcriptome and energy metabolism, thereby causing ATP deficiency and impairing pollen development.

RcSPL1-RcTAF15b regulates the flowering time of rose (Rosa chinensis)

Rose (Rosa chinensis), which is an economically valuable floral species worldwide, has three types, namely once-flowering (OF), occasional or re-blooming (OR), and recurrent or continuous flowering (CF). However, the mechanism underlying the effect of the age pathway on the duration of the CF or OF juvenile phase is largely unknown. In this study, we observed that the RcSPL1 transcript levels were substantially upregulated during the floral development period in CF and OF plants. Additionally, accumulation of RcSPL1 protein was controlled by rch-miR156. The ectopic expression of RcSPL1 in Arabidopsis thaliana accelerated the vegetative phase transition and flowering. Furthermore, the transient overexpression of RcSPL1 in rose plants accelerated flowering, whereas silencing of RcSPL1 had the opposite phenotype. Accordingly, the transcription levels of floral meristem identity genes (APETALA1, FRUITFULL, and LEAFY) were significantly affected by the changes in RcSPL1 expression. RcTAF15b protein, which is an autonomous pathway protein, was revealed to interact with RcSPL1. The silencing and overexpression of RcTAF15b in rose plants led to delayed and accelerated flowering, respectively. Collectively, the study findings imply that RcSPL1-RcTAF15b modulates the flowering time of rose plants.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Antagonistic regulation of target genes by the SISTER OF TM3-JOINTLESS2 complex in tomato inflorescence branching

Inflorescence branch number is a yield-related trait controlled by cell fate determination in meristems. Two MADS-box transcription factors (TFs)-SISTER OF TM3 (STM3) and JOINTLESS 2 (J2)-have opposing regulatory roles in inflorescence branching. However, the mechanisms underlying their regulatory functions in inflorescence determinacy remain unclear. Here, we characterized the functions of these TFs in tomato (Solanum lycopersicum) floral meristem and inflorescence meristem (IM) through chromatin immunoprecipitation and sequencing analysis of their genome-wide occupancy. STM3 and J2 activate or repress the transcription of a set of common putative target genes, respectively, through recognition and binding to CArG box motifs. FRUITFULL1 (FUL1) is a shared putative target of STM3 and J2 and these TFs antagonistically regulate FUL1 in inflorescence branching. Moreover, STM3 physically interacts with J2 to mediate its cytosolic redistribution and restricts J2 repressor activity by reducing its binding to target genes. Conversely, J2 limits STM3 regulation of target genes by transcriptional repression of the STM3 promoter and reducing STM3-binding activity. Our study thus reveals an antagonistic regulatory relationship in which STM3 and J2 control tomato IM determinacy and branch number.

The BELL1-like homeobox gene MdBLH14 from apple controls flowering and plant height via repression of MdGA20ox3

Apple growth and yield are largely dependent on plant height and flowering characteristics. The BELL1-like homeobox (BLH) transcription factors regulate extensive plant biological processes. However, the BLH-mediated regulation of plant height and flowering in apple remains elusive. In the current study, 19 members of the MdBLH family were identified in the apple genome. Segmental duplication and purifying selection are the main reasons for the evolution of the MdBLH genes. A BLH1-like gene, MdBLH14, was isolated and functionally characterized. The MdBLH14 was preferentially expressed in flower buds, and downregulated during the floral induction period. The subcellular localization in tobacco leaves indicated that MdBLH14 is a nuclear protein. Overexpression of MdBLH14 in Arabidopsis led to a significant dwarfing and late-flowering phenotype by hindering active GA accumulation. Additionally, MdKNOX19, another member of the TALE superfamily, physically interacts with MdBLH14 and synergistically inhibits the expression of MdGA20ox3. This is the first report on the function of the MdBLH14 from apple, and its mechanism involving plant flower induction and growth. The data presented here provide a theoretical basis for genetically breeding new apple varieties.

The Transcription Factors WRKY41 and WRKY53 Mediate Early Flowering Induced by the Novel Plant Growth Regulator Guvermectin in Arabidopsis thaliana

Flowering is a crucial stage for plant reproductive success; therefore, the regulation of plant flowering has been widely researched. Although multiple well-defined endogenous and exogenous flowering regulators have been reported, new ones are constantly being discovered. Here, we confirm that a novel plant growth regulator guvermectin (GV) induces early flowering in Arabidopsis. Interestingly, our genetic experiments newly demonstrated that WRKY41 and its homolog WRKY53 were involved in GV-accelerated flowering as positive flowering regulators. Overexpression of WRKY41 or WRKY53 resulted in an early flowering phenotype compared to the wild type (WT). In contrast, the w41/w53 double mutants showed a delay in GV-accelerated flowering. Gene expression analysis showed that flowering regulatory genes SOC1 and LFY were upregulated in GV-treated WT, 35S:WRKY41, and 35S:WRKY53 plants, but both declined in w41/w53 mutants with or without GV treatment. Meanwhile, biochemical assays confirmed that SOC1 and LFY were both direct targets of WRKY41 and WRKY53. Furthermore, the early flowering phenotype of 35S:WRKY41 lines was abolished in the soc1 or lfy background. Together, our results suggest that GV plays a function in promoting flowering, which was co-mediated by WRKY41 and WRKY53 acting as new flowering regulators by directly activating the transcription of SOC1 and LFY in Arabidopsis.

Arabidopsis AGAMOUS-LIKE16 and SUPPRESSOR OF CONSTANS1 regulate the genome-wide expression and flowering time

Flowering transition is tightly coordinated by complex gene regulatory networks, in which AGAMOUS-LIKE 16 (AGL16) plays important roles. Here, we identified the molecular function and binding properties of AGL16 and demonstrated its partial dependency on the SUPPRESSOR OF CONSTANS 1 (SOC1) function in regulating flowering. AGL16 bound to promoters of more than 2,000 genes via CArG-box motifs with high similarity to that of SOC1 in Arabidopsis (Arabidopsis thaliana). Approximately 70 flowering genes involved in multiple pathways were potential targets of AGL16. AGL16 formed a protein complex with SOC1 and shared a common set of targets. Intriguingly, only a limited number of genes were differentially expressed in the agl16-1 loss-of-function mutant. However, in the soc1-2 knockout background, AGL16 repressed and activated the expression of 375 and 182 genes, respectively, with more than a quarter bound by AGL16. Corroborating these findings, AGL16 repressed the flowering time more strongly in soc1-2 than in the Col-0 background. These data identify a partial inter-dependency between AGL16 and SOC1 in regulating genome-wide gene expression and flowering time, while AGL16 provides a feedback regulation on SOC1 expression. Our study sheds light on the complex background dependency of AGL16 in flowering regulation, thus providing additional insights into the molecular coordination of development and environmental adaptation.

Lysine 27 of histone H3.3 is a fine modulator of developmental gene expression and stands as an epigenetic checkpoint for lignin biosynthesis in Arabidopsis

Chromatin is a dynamic platform within which gene expression is controlled by epigenetic modifications, notably targeting amino acid residues of histone H3. Among them is lysine 27 of H3 (H3K27), the trimethylation of which by the Polycomb Repressive Complex 2 (PRC2) is instrumental in regulating spatiotemporal patterns of key developmental genes. H3K27 is also subjected to acetylation and is found at sites of active transcription. Most information on the function of histone residues and their associated modifications in plants was obtained from studies of loss-of-function mutants for the complexes that modify them. To decrypt the genuine function of H3K27, we expressed a non-modifiable variant of H3 at residue K27 (H3.3^K27A ) in Arabidopsis, and developed a multi-scale approach combining in-depth phenotypical and cytological analyses, with transcriptomics and metabolomics. We uncovered that the H3.3^K27A variant causes severe developmental defects, part of them are reminiscent of PRC2 mutants, part of them are new. They include early flowering, increased callus formation and short stems with thicker xylem cell layer. This latest phenotype correlates with mis-regulation of phenylpropanoid biosynthesis. Overall, our results reveal novel roles of H3K27 in plant cell fates and metabolic pathways, and highlight an epigenetic control point for elongation and lignin composition of the stem.

Flowering in sugarcane-insights from the grasses

Flowering is a crucial phase for angiosperms to continue their species propagation and is highly regulated. In the current review, flowering in sugarcane and the associated mechanisms are elaborately presented. In sugarcane, flowering has two effects, wherein it is a beneficial factor from the breeder's perspective and crucial for crop improvement, but commercially, it depletes the sucrose reserves from the stalks; hence, less value is assigned. Different species of Saccharum genus are spread across geographical latitudes, thereby proving their ability to grow in multiple inductive daylengths of different locations according in the habituated zone. In general, sugarcane is termed an intermediate daylength plant with quantitative short-day behaviour as it requires reduction in daylength from 12 h 55 min to 12 h or 12 h 30 min. The prime concern in sugarcane flowering is its erratic flowering nature. The transition to reproductive stage which reverts to vegetative stage if there is any deviation from ambient temperature and light is also an issue. Spatial and temporal gene expression patterns during vegetative to reproductive stage transition and after reverting to vegetative state could possibly reveal how the genetic circuits are being governed. This review will also shed a light on potential roles of genes and/or miRNAs in flowering in sugarcane. Knowledge of transcriptomic background of circadian, photoperiod, and gibberellin pathways in sugarcane will enable us to better understand of variable response in floral development.

Multiple transcriptome comparisons reveal the essential roles of FLOWERING LOCUS T in floral initiation and SOC1 and SVP in floral activation in blueberry

The flowering mechanisms, especially chilling requirement-regulated flowering, in deciduous woody crops remain to be elucidated. Flower buds of northern highbush blueberry cultivar Aurora require approximately 1,000 chilling hours to bloom. Overexpression of a blueberry FLOWERING LOCUS T (VcFT) enabled precocious flowering of transgenic “Aurora” mainly in non-terminated apical buds during flower bud formation, meanwhile, most of the mature flower buds could not break until they received enough chilling hours. In this study, we highlighted two groups of differentially expressed genes (DEGs) in flower buds caused by VcFT overexpression (VcFT-OX) and full chilling. We compared the two groups of DEGs with a focus on flowering pathway genes. We found: 1) In non-chilled flower buds, VcFT-OX drove a high VcFT expression and repressed expression of a major MADS-box gene, blueberry SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (VcSOC1) resulting an increased VcFT/VcSOC1 expression ratio; 2) In fully chilled flower buds that are ready to break, the chilling upregulated VcSOC1 expression in non-transgenic “Aurora” and repressed VcFT expression in VcFT-OX “Aurora”, and each resulted in a decreased ratio of VcFT to VcSOC1; additionally, expression of a blueberry SHORT VEGETATIVE PHASE (VcSVP) was upregulated in chilled flower buds of both transgenic and non-transgenic’ “Aurora”. Together with additional analysis of VcFT and VcSOC1 in the transcriptome data of other genotypes and tissues, we provide evidence to support that VcFT expression plays a significant role in promoting floral initiation and that VcSOC1 expression is a key floral activator. We thus propose a new hypothesis on blueberry flowering mechanism, of which the ratios of VcFT-to-VcSOC1 at transcript levels in the flowering pathways determine flower bud formation and bud breaking. Generally, an increased VcFT/VcSOC1 ratio or increased VcSOC1 in leaf promotes precocious flowering and flower bud formation, and a decreased VcFT/VcSOC1 ratio with increased VcSOC1 in fully chilled flower buds contributes to flower bud breaking.

BBX7 interacts with BBX8 to accelerate flowering in chrysanthemum

The quantitative control of FLOWERING LOCUS T (FT) activation is important for the floral transition in flowering plants. However, the flowering regulation mechanisms in the day-neutral, summer-flowering chrysanthemum plant remain unclear. In this study, the chrysanthemum BBX7 homolog CmBBX7 was isolated and its flowering function was identified. The expression of CmBBX7 showed a diurnal rhythm and CmBBX7 exhibited higher expression levels than CmBBX8. Overexpression of CmBBX7 in transgenic chrysanthemum accelerated flowering, whereas lines transfected with a chimeric repressor (pSRDX-CmBBX7) exhibited delayed flowering. Yeast single hybridization, luciferase, electrophoretic mobility shift, and chromatin immunoprecipitation assays showed that CmBBX7 directly targets CmFTL1. In addition, we found that CmBBX7 and CmBBX8 interact to positively regulate the expression of CmFTL1 through binding to its promoter. Collectively, these results highlight CmBBX7 as a key cooperator in the BBX8-FT module to control chrysanthemum flowering.

Genome-wide identification of the MIKCc-type MADS-box gene family in Chrysanthemum lavandulifolium reveals their roles in the capitulum development

Chrysanthemum ×morifolium is well known throughout the world for its diverse and exquisite flower types. However, due to the complicated genetic background of C. ×morifolium, it is difficult to understand the molecular mechanism of its flower development. And it limits the molecular breeding of improving chrysanthemum flower types. C. ×morifolium has the typical radial capitulum, and many researches showed that the members of the MIKCc-type MADS box gene family play a key role in the formation and development of the capitulum. However, it has been difficult to isolate the important MIKCc and investigate their roles in this process due to the lack of genomic information in chrysanthemum. Here, we identified MIKCc-type MADS box genes at whole genome-wide level in C. lavandulifolium, a diploid species closely related to C. ×morifolium, and investigated their roles in capitulum development by gene expression pattern analysis and protein interaction analysis. A total of 40 ClMIKCc were identified and were phylogenetically grouped into 12 clades. Members of all clades showed different enriched expression patterns during capitulum formation. We speculate that the E-class genes in C. lavandulifolium underwent subfunctionalization because they have a significantly expanded, more diverse expression patterns, and specifically tissue expression than AtSEPs. Meanwhile, we detected the C-class expressed in disc floret corolla, which could be the clue to explore the morphological differences between disc and ray floret corolla. In addition, the potential roles of some MIKCcs in complex inflorescence formation were explored by comparing the number and phylogenetic relationship of MIKCc subfamily members in Asteraceae with different capitulum types. Members of the FLC branch in Asteraceae were found to be possibly related to the differentiation and development of the ray floret.

Co-regulatory effects of hormone and mRNA-miRNA module on flower bud formation of Camellia oleifera

Few flower buds in a high-yield year are the main factors restricting the yield of Camellia oleifera in the next year. However, there are no relevant reports on the regulation mechanism of flower bud formation. In this study, hormones, mRNAs, and miRNAs were tested during flower bud formation in MY3 ("Min Yu 3," with stable yield in different years) and QY2 ("Qian Yu 2," with less flower bud formation in a high-yield year) cultivars. The results showed that except for IAA, the hormone contents of GA₃, ABA, tZ, JA, and SA in the buds were higher than those in the fruit, and the contents of all hormones in the buds were higher than those in the adjacent tissues. This excluded the effect of hormones produced from the fruit on flower bud formation. The difference in hormones showed that 21-30 April was the critical period for flower bud formation in C. oleifera; the JA content in MY3 was higher than that in QY2, but a lower concentration of GA₃ contributed to the formation of the C. oleifera flower bud. JA and GA₃ might have different effects on flower bud formation. Comprehensive analysis of the RNA-seq data showed that differentially expressed genes were notably enriched in hormone signal transduction and the circadian system. Flower bud formation in MY3 was induced through the plant hormone receptor TIR1 (transport inhibitor response 1) of the IAA signaling pathway, the miR535-GID1c module of the GA signaling pathway, and the miR395-JAZ module of the JA signaling pathway. In addition, the expression of core clock components GI (GIGANTEA) and CO (CONSTANS) in MY3 increased 2.3-fold and 1.8-fold over that in QY2, respectively, indicating that the circadian system also played a role in promoting flower bud formation in MY3. Finally, the hormone signaling pathway and circadian system transmitted flowering signals to the floral meristem characteristic genes LFY (LEAFY) and AP1 (APETALA 1) via FT (FLOWERING LOCUS T) and SOC1 (SUPPRESSOR OF OVEREXPRESSION OF CO 1) to regulate flower bud formation. These data will provide the basis for understanding the mechanism of flower bud alternate formation and formulating high yield regulation measures for C. oleifera.

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.

The N-terminal region located upstream of the TCP domain is responsible for the antagonistic action of the Arabidopsis thaliana TCP8 and TCP23 transcription factors on flowering time

TCP proteins (TCPs) are plant-exclusive transcription factors that exert effects on multiple aspects of plant development, from germination to flower and fruit formation. TCPs are divided into two main classes, I and II. In this study, we found that the Arabidopsis thaliana class I TCP transcription factor TCP8 is a positive regulator of flowering time. TCP8 mutation and constitutive expression delayed and accelerated flowering, respectively. Accordingly, TCP8 mutant plants showed a delay in the maximum expression of FT and reduced SOC1 transcript levels, while plants overexpressing TCP8 presented increased transcript levels of both genes. Notably, the related class I protein TCP23 showed the opposite behavior, since TCP23 mutation and overexpression accelerated and retarded flowering, respectively. To elucidate the molecular basis of these differences, we analyzed TCP8 and TCP23 comparatively. We found that both proteins are able to physically interact and bind class I TCP motifs, but only TCP8 shows transcriptional activation activity when expressed in plants, which is negatively affected by TCP23. From the analysis of plants expressing different chimeras between the TCPs, we found that the N-terminal region located upstream of the TCP domain is responsible for the opposite effect that TCP8 and TCP23 exert over flowering time and regulation of FT and SOC1 expression. These results suggest that structural features outside the TCP domain modulate the specificity of action of class I TCPs.

Emergent molecular traits of lettuce and tomato grown under wavelength-selective solar cells

The integration of semi-transparent organic solar cells (ST-OSCs) in greenhouses offers new agrivoltaic opportunities to meet the growing demands for sustainable food production. The tailored absorption/transmission spectra of ST-OSCs impacts the power generated as well as crop growth, development and responses to the biotic and abiotic environments. To characterize crop responses to ST-OSCs, we grew lettuce and tomato, traditional greenhouse crops, under three ST-OSC filters that create different light spectra. Lettuce yield and early tomato development are not negatively affected by the modified light environment. Our genomic analysis reveals that lettuce production exhibits beneficial traits involving nutrient content and nitrogen utilization while select ST-OSCs impact regulation of flowering initiation in tomato. These results suggest that ST-OSCs integrated into greenhouses are not only a promising technology for energy-neutral, sustainable and climate-change protected crop production, but can deliver benefits beyond energy considerations.

Epistatic interactions among multiple copies of FLC genes with naturally occurring insertions correlate with flowering time variation in radish

Brassicaceae crops, which underwent whole-genome triplication during their evolution, have multiple copies of flowering-related genes. Interactions among multiple gene copies may be involved in flowering time regulation; however, this mechanism is poorly understood. In this study, we performed comprehensive, high-throughput RNA sequencing analysis to identify candidate genes involved in the extremely late-bolting (LB) trait in radish. Then, we examined the regulatory roles and interactions of radish FLOWERING LOCUS C (RsFLC) paralogs, the main flowering repressor candidates. Seven flowering integrator genes, five vernalization genes, nine photoperiodic/circadian clock genes and eight genes from other flowering pathways were differentially expressed in the early-bolting (EB) cultivar 'Aokubinagafuto' and LB radish cultivar 'Tokinashi' under different vernalization conditions. In the LB cultivar, RsFLC1 and RsFLC2 expression levels were maintained after 40 days of cold exposure. Bolting time was significantly correlated with the expression rates of RsFLC1 and RsFLC2. Using the EB × LB F2 population, we performed association analyses of genotypes with or without 1910- and 1627-bp insertions in the first introns of RsFLC1 and RsFLC2, respectively. The insertion alleles prevented the repression of their respective FLC genes under cold conditions. Interestingly, genotypes homozygous for RsFLC2 insertion alleles maintained high RsFLC1 and RsFLC3 expression levels under cold conditions, and two-way analysis of variance revealed that RsFLC1 and RsFLC3 expression was influenced by the RsFLC2 genotype. Our results indicate that insertions in the first introns of RsFLC1 and RsFLC2 contribute to the late-flowering trait in radish via different mechanisms. The RsFLC2 insertion allele conferred a strong delay in bolting by inhibiting the repression of all three RsFLC genes, suggesting that radish flowering time is determined by epistatic interactions among multiple FLC gene copies.

Their C-termini divide Brassica rapa FT-like proteins into FD-interacting and FD-independent proteins that have different effects on the floral transition

Members of the FLOWERING LOCUS T (FT)-like clade of phosphatidylethanolamine-binding proteins (PEBPs) induce flowering by associating with the basic leucine zipper (bZIP) transcription factor FD and forming regulatory complexes in angiosperm species. However, the molecular mechanism of the FT-FD heterocomplex in Chinese cabbage (Brassica rapa ssp. pekinensis) is unknown. In this study, we identified 12 BrPEBP genes and focused our functional analysis on four BrFT-like genes by overexpressing them individually in an FT loss-of-function mutant in Arabidopsis thaliana. We determined that BrFT1 and BrFT2 promote flowering by upregulating the expression of floral meristem identity genes, whereas BrTSF and BrBFT, although close in sequence to their Arabidopsis counterparts, had no clear effect on flowering in either long- or short-day photoperiods. We also simultaneously genetically inactivated BrFT1 and BrFT2 in Chinese cabbage using CRISPR/Cas9-mediated genome editing, which revealed that BrFT1 and BrFT2 may play key roles in inflorescence organogenesis as well as in the transition to flowering. We show that BrFT-like proteins, except for BrTSF, are functionally divided into FD interactors and non-interactors based on the presence of three specific amino acids in their C termini, as evidenced by the observed interconversion when these amino acids are mutated. Overall, this study reveals that although BrFT-like homologs are conserved, they may have evolved to exert functionally diverse functions in flowering via their potential to be associated with FD or independently from FD in Brassica rapa.

Systematic analysis of MADS-box gene family in the U's triangle species and targeted mutagenesis of BnaAG homologs to explore its role in floral organ identity in Brassica napus

MADS-box transcription factors play an important role in regulating floral organ development and participate in environmental responses. To date, the MADS-box gene family has been widely identified in Brassica rapa (B. rapa), Brassica oleracea (B. oleracea), and Brassica napus (B. napus); however, there are no analogous reports in Brassica nigra (B. nigra), Brassica juncea (B. juncea), and Brassica carinata (B. carinata). In this study, a whole-genome survey of the MADS-box gene family was performed for the first time in the triangle of U species, and a total of 1430 MADS-box genes were identified. Based on the phylogenetic relationship and classification of MADS-box genes in Arabidopsis thaliana (A. thaliana), 1430 MADS-box genes were categorized as M-type subfamily (627 genes), further divided into Mα, Mβ, Mγ, and Mδ subclades, and MIKC-type subfamily (803 genes), further classified into 35 subclades. Gene structure and conserved protein motifs of MIKC-type MADS-box exhibit diversity and specificity among different subclades. Comparative analysis of gene duplication events and syngenic gene pairs among different species indicated that polyploidy is beneficial for MIKC-type gene expansion. Analysis of transcriptome data within diverse tissues and stresses in B. napus showed tissue-specific expression of MIKC-type genes and a broad response to various abiotic stresses, particularly dehydration stress. In addition, four representative floral organ mutants (wtl, feml, aglf-2, and aglf-1) in the T0 generation were generated by editing four AGAMOUS (BnaAG) homoeologs in B. napus that enriched the floral organ variant phenotype. In brief, this study provides useful information for investigating the function of MADS-box genes and contributes to revealing the regulatory mechanisms of floral organ development in the genetic improvement of new varieties.

Transcriptome analysis of critical genes related to flowering in Mikania micrantha at different altitudes provides insights for a potential control

BACKGROUND

Mikania micrantha is a vine with strong invasion ability, and its strong sexual reproduction ability is not only the main factor of harm, but also a serious obstacle to control. M. micrantha spreads mainly through seed production. Therefore, inhibiting the flowering and seed production of M. micrantha is an effective strategy to prevent from continuing to spread.

RESULT

The flowering number of M. micrantha is different at different altitudes. A total of 67.01 Gb of clean data were obtained from nine cDNA libraries, and more than 83.47% of the clean reads were mapped to the reference genome. In total, 5878 and 7686 significantly differentially expressed genes (DEGs) were found in E2 vs. E9 and E13 vs. E9, respectively. Based on the background annotation and gene expression, some candidate genes related to the flowering pathway were initially screened, and their expression levels in the three different altitudes in flower bud differentiation showed the same trend. That is, at an altitude of 1300 m, the flower integration gene and flower meristem gene were downregulated (such as SOC1 and AP1), and the flowering inhibition gene was upregulated (such as FRI and SVP). Additionally, the results showed that there were many DEGs involved in the hormone signal transduction pathway in the flower bud differentiation of M. micrantha at different altitudes.

CONCLUSIONS

Our results provide abundant sequence resources for clarifying the underlying mechanisms of flower bud differentiation and mining the key factors inhibiting the flowering and seed production of M. micrantha to provide technical support for the discovery of an efficient control method.

Metabolic analysis of the regulatory mechanism of sugars on secondary flowering in Magnolia

Background

Magnolia, a traditional and important ornamental plant in urban greening, has been cultivated for about 2000 years in China for its elegant flower shape and gorgeous flower color. Most varieties of Magnolia bloom once a year in spring, whereas a few others, such as Magnolia liliiflora Desr. ‘Hongyuanbao’, also bloom for the second time in summer or early autumn. Such a twice flowering trait is desirable for its high ornamental value, while its underlying mechanism remains unclear.

Methods

Paraffin section was used to show the flowering time and phenotypic changes of M. liliiflora ‘Hongyuanbao’ during the twice flowering periods from March 28 to August 25, 2018. Gas chromatography-mass spectrometry (GC-MS) was then performed to explore the chemical metabolites through the twice flower bud differentiation process in ‘Hongyuanbao’, and the metabolites were screened and identified by orthogonal projection to latent structures discriminant analysis (OPLS-DA). Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis (KEGG) was used to reveal the relationship between the sugar metabolites and twice-flowering characteristic. To further investigate the potential role of sucrose and trehalose on flowering regulation of ‘Hongyuanbao’, the plants once finished the spring flowering were regularly sprayed with sucrose and trehalose solutions at 30 mM, 60 mM, and 90 mM concentrations from April 22, 2019. The flower bud differentiation processes of sprayed plants were observed and the expression patterns of the genes involved in sucrose and trehalose metabolic pathways were studied by quantitative reverse transcription PCR (qRT-PCR).

Results

It showed that ‘Hongyuanbao’ could complete flower bud differentiation twice in a year and flowered in both spring and summer. The metabolites of flower bud differentiation had a significant variation between the first and second flower buds. Compared to the first flower bud differentiation process, the metabolites in the sucrose and trehalose metabolic pathways were significantly up-regulated during the second flower bud differentiation process. Besides that, the expression levels of a number of trehalose-6-phosphate synthase (TPS) genes including MlTPS1, MlTPS5, MlTPS6, MlTPS7 and MlTPS9 were substantially increased in the second flower differentiation process compared with the first process. Exogenous treatments indicated that compared to the control plants (sprayed with water, CK), all three concentrations of trehalose could accelerate flowering and the effect of 60 mM concentration was the most significant. For the sucrose foliar spray, only the 60 mM concentration accelerated flowering compared with CK. It suggested that different concentration of trehalose and sucrose might have different effects. Expression analysis showed that sucrose treatment increased the transcription levels of MlTPS5 and MlTPS6, whereas trehalose treatment increased MlTPS1, showing that different MlTPS genes took part in sucrose and trehalose metabolic pathways respectively. The expression levels of a number of flowering-related genes, such as MlFT, MlLFY, and MlSPL were also increased in response to the sprays of sucrose and trehalose.

Conclusions

We provide a novel insight into the effect of sucrose and trehalose on the flowering process in Magnolia. Under the different sugar contents treatments, the time of flower bud differentiation of Magnolia was advanced. Induced and accelerated flowering in response to sucrose and trehalose foliar spray, coupled with elevated expression of trehalose regulatory and response genes, suggests that secondary flower bud formation is a promoted by altered endogenous sucrose and trehalose levels. Those results give a new understanding of sucrose and trehalose on twice-flowering in Magnolia and provide a preliminary speculation for inducing and accelerating the flowering process in Magnolia.

Drought and low temperature-induced NF-YA1 activates FT expression to promote citrus flowering

Flower induction in adult citrus is mainly regulated by drought and low temperatures. However, the mechanism of FLOWERING LOCUS T regulation of citrus flowering (CiFT) under two flower-inductive stimuli remains largely unclear. In this study, a citrus transcription factor, nuclear factor YA (CiNF-YA1), was found to specifically bind to the CiFT promoter by forming a complex with CiNF-YB2 and CiNF-YC2 to activate CiFT expression. CiNF-YA1 was induced in juvenile citrus by low temperature and drought treatments. Overexpression of CiNF-YA1 increased drought susceptibility in transgenic citrus, whereas suppression of CiNF-YA1 enhanced drought tolerance in silenced citrus plants. Furthermore, a GOLDEN₂ ^- LIKE protein (CiFE) that interacts with CiFT protein was also isolated. Further experimental evidence showed that CiFE binds to the citrus LEAFY (CiLFY) promoter and activates its expression. In addition, the expressions of CiNF-YA1 and CiFE showed a seasonal increase during the floral induction period and were induced by artificial drought and low-temperature treatments at which floral induction occurred. These results indicate that CiNF-YA1 may activate CiFT expression in response to drought and low temperatures by binding to the CiFT promoter. CiFT then forms a complex with CiFE to activate CiLFY, thereby promoting the flowering of adult citrus.

CRISPR/Cas9-Mediated Targeted Mutagenesis of BnaCOL9 Advances the Flowering Time of Brassica napus L

Rapeseed (Brassica napus L.) is one of the most important oil crops in the world. The planting area and output of rapeseed are affected by the flowering time, which is a critical agronomic feature. COL9 controls growth and development in many different plant species as a member of the zinc finger transcription factor family. However, BnaCOL9 in rapeseed has not been documented. The aim of this study was to apply CRISPR/Cas9 technology to create an early-flowering germplasm resource to provide useful material for improving the early-maturing breeding of rapeseed. We identified four COL9 homologs in rapeseed that were distributed on chromosomes A05, C05, A03, and C03. We successfully created quadruple BnaCOL9 mutations in rapeseed using the CRISPR/Cas9 platform. The quadruple mutants of BnaCOL9 flowered earlier than the wild-type. On the other hand, the flowering time of the BnaCOL9 overexpression lines was delayed. An analysis of the expression patterns revealed that these genes were substantially expressed in the leaves and flowers. A subcellular localization experiment demonstrated that BnaCOL9 was in the nucleus. Furthermore, we discovered that two key flowering-related genes, BnaCO and BnaFT, were highly elevated in the BnaCOL9 mutants, but dramatically downregulated in the BnaCOL9 overexpression lines. Our findings demonstrate that BnaCOL9 is a significant flowering inhibitor in rapeseed and may be employed as a crucial gene for early-maturing breeding.

Genome-wide characterization of PEBP gene family in Perilla frutescens and PfFT1 promotes flowering time in Arabidopsis thaliana

Phosphatidylethanolamine-binding proteins (PEBP) family plays important roles in regulating plant flowering time and morphogenesis. However, geneme-wide identification and functional analysis of PEBP genes in the rigorous short-day plant Perilla frutescens (PfPEBP) have not been studied. In this study, 10 PfPEBP were identified and divided into three subfamilies based on their phylogenetic relationships: FT-like, TFL1-like and MFT-like. Gene structure analysis showed that all PfPEBP genes contain 4 exons and 3 introns. Motifs DPDxP and GIHR essential for anion-binding activity are highly conserved in PfPEBP. A large number of light-responsive elements were detected in promoter regions of PfPEBP. Gene expression of PfFT1 exhibited a diurnal rhythm. It was highly expressed in leaves under the short-day photoperiod, but higher in flowers and seeds under the long-day photoperiod. Overexpression of PfFT1 in Arabidopsis thaliana not only promoted early flowering of Col-0 or Ler, but also rescued the late flowering phenotype of ft-1 mutant. We concluded that PfFT1 promotes early flowering by regulating the expression of flowering-related genes AtAP1, AtLFY, AtFUL and AtSOC1. In conclusion, our results provided valuable information for elucidating the functions of PfPEBP in P. frutescens and shed light on the promoting effect of PfFT1 on flowering.

Manipulating GA-Related Genes for Cereal Crop Improvement

The global population is projected to experience a rapid increase in the future, which poses a challenge to global food sustainability. The "Green Revolution" beginning in the 1960s allowed grain yield to reach two billion tons in 2000 due to the introduction of semi-dwarfing genes in cereal crops. Semi-dwarfing genes reduce the gibberellin (GA) signal, leading to short plant stature, which improves the lodging resistance and harvest index under modern fertilization practices. Here, we reviewed the literature on the function of GA in plant growth and development, and the role of GA-related genes in controlling key agronomic traits that contribute to grain yield in cereal crops. We showed that: (1) GA is a significant phytohormone in regulating plant development and reproduction; (2) GA metabolism and GA signalling pathways are two key components in GA-regulated plant growth; (3) GA interacts with other phytohormones manipulating plant development and reproduction; and (4) targeting GA signalling pathways is an effective genetic solution to improve agronomic traits in cereal crops. We suggest that the modification of GA-related genes and the identification of novel alleles without a negative impact on yield and adaptation are significant in cereal crop breeding for plant architecture improvement. We observed that an increasing number of GA-related genes and their mutants have been functionally validated, but only a limited number of GA-related genes have been genetically modified through conventional breeding tools and are widely used in crop breeding successfully. New genome editing technologies, such as the CRISPR/Cas9 system, hold the promise of validating the effectiveness of GA-related genes in crop development and opening a new venue for efficient and accelerated crop breeding.

Genome-wide identification and characterization profile of phosphatidy ethanolamine-binding protein family genes in carrot

Members of the family of Phosphatidy Ethanolamine-Binding Protein (PEBP) have been shown to be key regulators of the transition of plants from vegetative to reproductive phases. Here, a total of 12 PEBP proteins were identified in the carrot (Daucus carota L.) genome and classified into FT-like (4), TFL1-like (6), and MFT-like 2) subfamilies, that had different lengths (110-267 aa) and were distributed unevenly across seven chromosomes. Moreover, 13 and 31 PEBP proteins were identified in other two Apiaceae species, celery (Apium graveolens L.) and coriander (Coriandrum sativum L.). The phylogenetic and evolutionary results of these PEBP family proteins were obtained based on the protein sequences. In the three Apiaceae species, purifying selection was the main evolutionary force, and WGD, segmental duplication, and dispersed duplication have played key roles in the PEBP family expansion. The expression analysis showed that carrot PEBP genes exhibited relatively broad expression patterns across various tissues. In the period of bolting to flowering, the carrot FT-like subfamily genes were upregulated as positive regulators, and TFL1-like subfamily genes remained at lower expression levels as inhibitors. More interestingly, the members of carrot FT-like genes had different temporal-spatial expression characteristics, suggesting that they have different regulatory functions in the carrot reproductive phase. In summary, this study contributes to our understanding of the PEBP family proteins and provides a foundation for exploring the mechanism of carrot bolting and flowering for the breeding of cultivars with bolting resistance.

Profiling of phytohormones in apple fruit and buds regarding their role as potential regulators of flower bud formation

Apple (Malus × domestica Borkh.) cropping behavior, if not regulated, is often manifested by high yields of small-sized fruit in so called ON-years, which are usually followed by strongly reduced crop loads in OFF-years. Such cropping pattern is defined as biennial bearing and causes significant losses in apple production. The growth of apple fruit overlaps with the formation of flower buds, which remain dormant until the following spring. Earlier works proposed that some fruit-derived mobile compounds, as e.g., phytohormones, could suppress flower bud formation that thereby leads to biennial bearing. We addressed this hypothesis by analyzing 39 phytohormones in apple seeds, fruit flesh and by measuring phytohormone export from the fruits of the biennial bearing cultivar 'Fuji' and of the regular bearing cultivar 'Gala'. Moreover, we analyzed the same compounds in bourse buds from fruiting (ON-trees) and non-fruiting (OFF-trees) spurs of both apple cultivars over the period of flower bud formation. Our results showed that apple fruit exported at least 14 phytohormones including indole-3-acetic acid and gibberellin A3; however, their influence on flower bud formation was inconclusive. A gibberellin-like compound, which was detected exclusively in bourse buds, was significantly more abundant in bourse buds from ON-trees compared with OFF-trees. Cultivar differences were marked by the accumulation of trans-zeatin-O-glucoside in bourse buds of 'Gala' ON-trees, whereas the levels of this compound in 'Gala' OFF were significantly lower and comparable to those in 'Fuji' ON- and OFF-trees. Particular phytohormones including five cytokinin forms as well as abscisic acid and its degradation products had higher levels in bourse buds from OFF-trees compared with ON-trees and were therefore proposed as potential promotors of flower bud initiation. The work discusses regulatory roles of phytohormones in flower bud formation in apple based on the novel and to date most comprehensive phytohormone profiles of apple fruit and buds.

MIKC-type MADS-box transcription factor gene family in peanut: Genome-wide characterization and expression analysis under abiotic stress

Peanut (Arachis hypogaea) is one of the most important economic crops around the world, especially since it provides vegetable oil and high-quality protein for humans. Proteins encoded by MADS-box transcription factors are widely involved in regulating plant growth and development as well as responses to abiotic stresses. However, the MIKC-type MADS-box TFs in peanut remains currently unclear. Hence, in this study, 166 MIKC-type MADS-box genes were identified in both cultivated and wild-type peanut genomes, which were divided into 12 subfamilies. We found a variety of development-, hormone-, and stress-related cis-acting elements in the promoter region of peanut MIKC-type MADS-box genes. The chromosomal distribution of peanut MADS-box genes was not random, and gene duplication contributed to the expansion of the MADS-box gene family. The interaction network of the peanut AhMADS proteins was established. Expression pattern analysis showed that AhMADS genes were specifically expressed in tissues and under abiotic stresses. It was further confirmed via the qRT-PCR technique that five selected AhMADS genes could be induced by abiotic and hormone treatments and presented different expressive profiles under various stresses. Taken together, these findings provide valuable information for the exploration of candidate genes in molecular breeding and further study of AhMADS gene functions.

Divergence of flowering-related genes to control flowering in five Euphorbiaceae genomes

Reproductive growth and vegetative growth are a pair of main contradictions in the process of plant growth. Flowering, as part of reproductive growth, is a key switch in the life cycle of higher plants, which affects the yield and economic benefits of plants to a certain extent. The Euphorbiaceae species, including castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), cassava (Manihot esculenta), and rubber tree (Hevea brasiliensis), have important economic values because they are raw materials for the production of biodiesel, rubber, etc. The flowering mechanisms are still excluded in the Euphorbiaceae species. The flowering-related genes of Arabidopsis thaliana (Arabidopsis) were used as a reference to determine the orthologs of these genes in Euphorbiaceae genomes. The result showed that 146, 144, 114, 114, and 149 of 207 A. thaliana genes were respectively matched to R. communis, V. fordii, J. curcas, H. brasiliensis, and M. esculenta. These identified genes were clustered into seven pathways including gibberellins, floral meristem identity (FMI), vernalization, photoperiod, floral pathway integrators (FPIs), and autonomous pathways. Then, some key numbers of flowering-related genes are widely conserved in the Euphorbiaceae genomes including but not limited to FPI genes LFY, SOC1, FT, and FMI genes AG, CAL, and FUL. However, some genes, including FRI, FLC, and GO, were missing in several or all five Euphorbiaceae species. In this study, we proposed the putative mechanisms of flowering-related genes to control flowering and provided new candidate flowering genes for using marker-assisted breeding to improve variety quality.

Comprehensive collection of genes and comparative analysis of full-length transcriptome sequences from Japanese larch (Larix kaempferi) and Kuril larch (Larix gmelinii var. japonica)

BACKGROUND

Japanese larch (Larix kaempferi) is an economically important deciduous conifer species that grows in cool-temperate forests and is endemic to Japan. Kuril larch (L. gmelinii var. japonica) is a variety of Dahurian larch that is naturally distributed in the Kuril Islands and Sakhalin. The hybrid larch (L. gmelinii var. japonica × L. kaempferi) exhibits heterosis, which manifests as rapid juvenile growth and high resistance to vole grazing. Since these superior characteristics have been valued by forestry managers, the hybrid larch is one of the most important plantation species in Hokkaido. To accelerate molecular breeding in these species, we collected and compared full-length cDNA isoforms (Iso-Seq) and RNA-Seq short-read, and merged them to construct candidate gene as reference for both Larix species. To validate the results, candidate protein-coding genes (ORFs) related to some flowering signal-related genes ​were screened from the reference sequences, and the phylogenetic relationship with closely related species was elucidated.

RESULTS

Using the isoform sequencing of PacBio RS ll and the de novo assembly of RNA-Seq short-read sequences, we identified 50,690 and 38,684 ORFs in Japanese larch and Kuril larch, respectively. BUSCO completeness values were 90.5% and 92.1% in the Japanese and Kuril larches, respectively. After comparing the collected ORFs from the two larch species, a total of 19,813 clusters, comprising 22,571 Japanese larch ORFs and 22,667 Kuril larch ORFs, were contained in the intersection of the Venn diagram. In addition, we screened several ORFs related to flowering signals (SUPPRESSER OF OVEREXPRESSION OF CO1: SOC1, LEAFY: LFY, FLOWERING Locus T: FT, CONSTANCE: CO) from both reference sequences, and very similar found in other species.

CONCLUSIONS

The collected ORFs will be useful as reference sequences for molecular breeding of Japanese and Kuril larches, and also for clarifying the evolution of the conifer genome and investigating functional genomics.

The tomato CONSTANS-LIKE protein SlCOL1 regulates fruit yield by repressing SFT gene expression

BACKGROUND

CONSTANS (CO) and CONSTANS-LIKE (COL) transcription factors have been known to regulate a series of cellular processes including the transition from the vegetative growth to flower development in plants. However, their role in regulating fruit yield in tomato is poorly understood.

RESULT

In this study, the tomato ortholog of Arabidopsis CONSTANS, SlCOL1, was shown to play key roles in the control of flower development and fruit yield. Suppression of SlCOL1 expression in tomato was found to lead to promotion of flower and fruit development, resulting in increased tomato fruit yield. On the contrary, overexpression of SlCOL1 disturbed flower and fruit development, and significantly reduced tomato fruit yield. Genetic and biochemical evidence indicated that SlCOL1 controls inflorescence development by directly binding to the promoter region of tomato inflorescence-associated gene SINGLE-FLOWER TRUSS (SFT) and negatively regulating its expression. Additionally, we found that SlCOL1 can also negatively regulate fruit size in tomato.

CONCLUSIONS

Tomato SlCOL1 binds to the promoter of the SFT gene, down-regulates its expression, and plays a key role in reducing the fruit size.

Mechanisms of Vernalization-Induced Flowering in Legumes

Vernalization is the requirement for exposure to low temperatures to trigger flowering. The best knowledge about the mechanisms of vernalization response has been accumulated for Arabidopsis and cereals. In Arabidopsis thaliana, vernalization involves an epigenetic silencing of the MADS-box gene FLOWERING LOCUS C (FLC), which is a flowering repressor. FLC silencing releases the expression of the main flowering inductor FLOWERING LOCUS T (FT), resulting in a floral transition. Remarkably, no FLC homologues have been identified in the vernalization-responsive legumes, and the mechanisms of cold-mediated transition to flowering in these species remain elusive. Nevertheless, legume FT genes have been shown to retain the function of the main vernalization signal integrators. Unlike Arabidopsis, legumes have three subclades of FT genes, which demonstrate distinct patterns of regulation with respect to environmental cues and tissue specificity. This implies complex mechanisms of vernalization signal propagation in the flowering network, that remain largely elusive. Here, for the first time, we summarize the available information on the genetic basis of cold-induced flowering in legumes with a special focus on the role of FT genes.

Genome-wide identification of nitrate transporter genes from Spirodela polyrhiza and characterization of SpNRT1.1 function in plant development

Nitrate transporter (NRT) genes that participate in nitrate transport and distribution are indispensable for plant growth, development, and stress tolerance. Spirodela polyrhiza has the smallest genome among monocotyledon plants, and it has strong nitrate absorbance and phytoremediation abilities. However, the evolutionary history, expression patterns, and functions of the NRT gene family in S. polyrhiza are not well understood. Here, we identified 29 NRT members in the S. polyrhiza genome. Gene structure and phylogeny analyses showed that S. polyrhiza nitrate transporter (SpNRTs) genes were divided into eight clades without gene expansion compared with that in Arabidopsis. Transcriptomic analysis showed that SpNRT genes have spatiotemporal expression patterns and respond to abiotic stress. Functional analysis revealed that in S. polyrhiza, SpNRT1.1 expression was strongly induced by treatment with nitrate and ammonium. Overexpression of SpNRT1.1 significantly repressed primary root length, and the number and total length of lateral roots. This was more pronounced in high ammonium concentration medium. Overexpressed SpNRT1.1 in Arabidopsis significantly improved biomass and delayed flowering time, indicating that the nitrate transport ability of SpNRT1.1 differs from AtNRT1.1. In conclusion, our results provide valuable information about the evolution of the NRT family in higher plants and the function of SpNRT1.1.

Genome-wide identification, phylogeny, and expression profiling analysis of shattering genes in rapeseed and mustard plants

BACKGROUND

Non-synchronized pods shattering in the Brassicaceae family bring upon huge yield losses around the world. The shattering process was validated to be controlled by eight genes in Arabidopsis, including SHP1, SHP2, FUL, IND, ALC, NAC, RPL, and PG. We performed genome-wide identification, characterization, and expression analysis of shattering genes in B.napus and B. juncea to gain understanding into this gene family and to explain their expression patterns in fresh and mature siliques.

RESULTS

A comprehensive genome investigation of B.napus and B.juncea revealed 32 shattering genes, which were identified and categorized using protein motif structure, exon-intron organization, and phylogeny. The phylogenetic study revealed that these shattering genes contain little duplications, determined with a distinct chromosome number. Motifs of 32 shattering proteins were observed where motifs1 and 2 were found to be more conserved. A single motif was observed for other genes like Br-nS7, Br-nS9, Br-nS10, Br-jS21, Br-jS23, Br-jS24, Br-jS25, and Br-jS26. Synteny analysis was performed that validated a conserved pattern of blocks among these cultivars. RT-PCR based expressions profiles showed higher expression of shattering genes in B. juncea as compared to B.napus. SHP1, SHP2, and FUL gene were expressed more in mature silique. ALC gene was upregulated in fresh silique of B. napus but downregulation of ALC were observed in fresh silique of B. juncea.

CONCLUSION

This study authenticates the presence of shattering genes in the local cultivars of Brassica. It has been validated that the expression of shattering genes were more in B. juncea as compared to B.napus. The outcomes of this study contribute to the screening of more candidate genes for further investigation.

Role of FaSOC1 and FaCO in the seasonal control of reproductive and vegetative development in the perennial crop Fragaria × ananassa

The diploid woodland strawberry (F. vesca) represents an important model for the genus Fragaria. Significant advances in the understanding of the molecular mechanisms regulating seasonal alternance of flower induction and vegetative reproduction has been made in this species. However, this research area has received little attention on the cultivated octoploid strawberry (F. × ananassa) despite its enormous agronomical and economic importance. To advance in the characterization of this intricated molecular network, expression analysis of key flowering time genes was performed both in short and long days and in cultivars with seasonal and perpetual flowering. Analysis of overexpression of FaCO and FaSOC1 in the seasonal flowering 'Camarosa' allowed functional validation of a number of responses already observed in F. vesca while uncovered differences related to the regulation of FaFTs expression and gibberellins (GAs) biosynthesis. While FvCO has been shown to promote flowering and inhibit runner development in the perpetual flowering H4 accession of F. vesca, our study showed that FaCO responds to LD photoperiods as in F. vesca but delayed flowering to some extent, possibly by induction of the strong FaTFL1 repressor in crowns. A contrasting effect on runnering was observed in FaCO transgenic plants, some lines showing reduced runner number whereas in others runnering was slightly accelerated. We demonstrate that the role of the MADS-box transcription factor FaSOC1 as a strong repressor of flowering and promoter of vegetative growth is conserved in woodland and cultivated strawberry. Our study further indicates an important role of FaSOC1 in the photoperiodic repression of FLOWERING LOCUS T (FT) genes FaFT2 and FaFT3 while FaTFL1 upregulation was less prominent than that observed in F. vesca. In our experimental conditions, FaSOC1 promotion of vegetative growth do not require induction of GA biosynthesis, despite GA biosynthesis genes showed a marked photoperiodic upregulation in response to long days, supporting GA requirement for the promotion of vegetative growth. Our results also provided insights into additional factors, such as FaTEM, associated with the vegetative developmental phase that deserve further characterization in the future.

Molecular basis of high temperature-induced bolting in lettuce revealed by multi-omics analysis

BACKGROUND

High temperature induces early bolting in lettuce (Lactuca sativa L.), which affects both quality and production. However, the molecular mechanism underlying high temperature-induced bolting is still limited.

RESULTS

We performed systematical analysis of morphology, transcriptome, miRNAs and methylome in lettuce under high temperature treatment. Through a comparison of RNA-Seq data between the control and the high temperature treated lettuces at different time points totally identified 2944 up-regulated genes and 2203 down-regulated genes, which cover three floral pathways including photoperiod, age and gibberellin (GA) pathways. Genome wide analysis of miRNAs and methylome during high temperature treatment indicated miRNAs and DNA methylation might play a role controlling gene expression during high temperature-induced bolting. miRNA targets included some protein kinase family proteins, which potentially play crucial roles in this process.

CONCLUSIONS

Together, our results propose a possible regulation network involved in high temperature-induced bolting.

Transcriptomic comparison sheds new light on regulatory networks for dimorphic flower development in response to photoperiod in Viola prionantha

BACKGROUND

Chasmogamous (CH)-cleistogamous (CL) dimorphic flowers are developed in Viola prionantha. However, the environmental and genetic factors necessary for the CH-CL transition are unknown.

RESULTS

In the present work, short-day (SD) conditions induced CH flowers, whereas long days (LDs) triggered CL flowers in V. prionantha. Compared to fully developed CH flowers, CL flowers had less mature stamens, no nectar glands, and immature petals. Comparative transcriptomics revealed differentially expressed genes (DEGs) during CL and CH development. Core genes in the photoperiod pathway, such as V. prionantha orthologs of GIGANTEA (GI), CONSTANS (CO), and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1), which promote floral induction, were highly expressed in CL flowers, whereas UNUSUAL FLORAL ORGANS (UFO) and B-class MADS-box genes for floral organ identity and development showed an opposite alteration. Moreover, genes in the glycolytic process, sucrose metabolic process, and fatty acid biosynthetic process were all highly expressed in CH flowers. Interestingly, V. prionantha orthologs of the B-class MADS-box genes APETALA3 (AP3) and PISTILLATA (PI) might relate to these sugar-fatty acid processes and were co-expressed with GAIP-B-like and YABBY5 (YAB5), which regulate the development of the petal, stamen, and nectary. Compared to CH flowers, DEGs and hub genes in the most significantly correlated modules of the gene co-expression network, which are involved in abiotic and biotic responses, were upregulated in CL flowers.

CONCLUSIONS

We proposed an integrative model for transcription regulation of genes in the photoperiod pathway, floral organ development, stress response, and sugar-fatty acid processes to determine CH-CL flower development in V. prionantha. Particularly, under LDs, activated GI may induce genes involved in the stress-response pathways, and then downregulated AP3 and PI or UFO to inhibit the sugar-fatty acid metabolic processes, together forming CL flowers. In contrast, CH flowers were produced under SDs. This work provides novel insights into the developmental evolution of dimorphic flowers in Viola.