Single-Molecule Real-Time and Illumina-Based RNA Sequencing Data Identified Vernalization-Responsive Candidate Genes in Faba Bean (Vicia faba L.)

Modeling the Flowering Activation Motif during Vernalization in Legumes: A Case Study of M. trancatula

In many plant species, flowering is promoted by the cold treatment or vernalization. The mechanism of vernalization-induced flowering has been extensively studied in Arabidopsis but remains largely unknown in legumes. The orthologs of the FLC gene, a major regulator of vernalization response in Arabidopsis, are absent or non-functional in the vernalization-sensitive legume species. Nevertheless, the legume integrator genes FT and SOC1 are involved in the transition of the vernalization signal to meristem identity genes, including PIM (AP1 ortholog). However, the regulatory contribution of these genes to PIM activation in legumes remains elusive. Here, we presented the theoretical and data-driven analyses of a feed-forward regulatory motif that includes a vernalization-responsive FT gene and several SOC1 genes, which independently activate PIM and thereby mediate floral transition. Our theoretical model showed that the multiple regulatory branches in this regulatory motif facilitated the elimination of no-sense signals and amplified useful signals from the upstream regulator. We further developed and analyzed four data-driven models of PIM activation in Medicago trancatula in vernalized and non-vernalized conditions in wild-type and fta1-1 mutants. The model with FTa1 providing both direct activation and indirect activation via three intermediate activators, SOC1a, SOC1b, and SOC1c, resulted in the most relevant PIM dynamics. In this model, the difference between regulatory inputs of SOC1 genes was nonessential. As a result, in the M. trancatula model, the cumulative action of SOC1a, SOC1b, and SOC1c was favored. Overall, in this study, we first presented the in silico analysis of vernalization-induced flowering in legumes. The considered vernalization network motif can be supplemented with additional regulatory branches as new experimental data become available.

Genome Size Variation across a Cypriot Fabeae Tribe Germplasm Collection

DNA content is an important trait linked to the evolutionary routes of taxa and often connected to speciation. In the present study, we studied C-values variation across the Cypriot Fabeae gene pool. Several hundred plants (Vicia spp., Lens spp., Pisum spp.) were sampled across Cyprus. Accurate estimates were established by flow cytometry and propidium iodine staining for 155 discrete populations/accessions. A ten-fold variation was detected across lineages with 1C DNA content varying from 1.584 pg for V. cretica (ARI02420) to 13.983 pg for V. faba (ARI00187). In general, flow cytometry was precise for the characterization of species, even though there were instances of genome overlapping across taxa. Most analyses in the current work refer to species that have not been characterized before by flow cytometry (or any other DNA content estimation method). Still, a correlation to C-values previously reported in Kew Plant DNA C-values database was attempted. A high degree of correlation except for V. dalmatica was established. The evaluation of genome size trait in relation with the Fabeae phylogeny, revealed that Pisum and Lens genera were rather homogenous, but an astonishing fluctuation was shown for Vicia spp. Moreover, it was established that genome up- or down-scaling was not directly linked to speciation drivers. The genomic size measurements presented here could deliver extra quality control for the identification and characterization of taxa in germplasm collections, particularly in cases where species share morphological characters.

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.

Unraveling the complexity of faba bean (Vicia faba L.) transcriptome to reveal cold-stress-responsive genes using long-read isoform sequencing technology

Faba bean (Vicia faba L.), a globally important grain legume providing a stable source of dietary protein, was one of the earliest plant cytogenetic models. However, the lack of draft genome annotations and unclear structural information on mRNA transcripts have impeded its genetic improvement. To address this, we sequenced faba bean leaf transcriptome using the PacBio single-molecule long-read isoform sequencing platform. We identified 28,569 nonredundant unigenes, ranging from 108 to 9669 bp, with a total length of 94.5 Mb. Many unigenes (3597, 12.5%) had 2-20 isoforms, indicating a highly complex transcriptome. Approximately 96.5% of the unigenes matched sequences in public databases. The predicted proteins and transcription factors included NB-ARC, Myb_domain, C3H, bHLH, and heat shock proteins, implying that this genome has an abundance of stress resistance genes. To validate our results, we selected WCOR413-15785, DHN2-12403, DHN2-14197, DHN2-14797, COR15-14478, and HVA22-15 unigenes from the ICE-CBF-COR pathway to analyze their expression patterns in cold-treated samples via qRT-PCR. The expression of dehydrin-related genes was induced by cold stress. The assembled data provide the first insights into the deep sequencing of full-length RNA from faba bean at the single-molecule level. This study provides an important foundation to improve gene modeling and protein prediction.