Transcriptome profile analysis of two Vicia faba cultivars with contrasting salinity tolerance during seed germination

Accurate Long-Read RNA Sequencing Analysis Reveals the Key Pathways and Candidate Genes under Drought Stress in the Seed Germination Stage in Faba Bean

Faba bean is an important pulse. It provides proteins for the human diet and is used in industrial foodstuffs, such as flours. Drought stress severely reduces the yield of faba bean, and this can be efficiently overcome through the identification and application of key genes in response to drought. In this study, PacBio and Illumina RNA sequencing techniques were used to identify the key pathways and candidate genes involved in drought stress response. During seed germination, a total of 17,927 full-length transcripts and 12,760 protein-coding genes were obtained. There were 1676 and 811 differentially expressed genes (DEGs) between the varieties E1 and C105 at 16 h and 64 h under drought stress, respectively. Six and nine KEGG pathways were significantly enriched at 16 h and 64 h under drought stress, which produced 40 and 184 nodes through protein-protein interaction (PPI) analysis, respectively. The DEGs of the PPI nodes were involved in the ABA (abscisic acid) and MAPK (mitogen-activated protein kinase) pathways, N-glycosylation, sulfur metabolism, and sugar metabolism. Furthermore, the ectopic overexpression of a key gene, AAT, encoding aspartate aminotransferase (AAT), in tobacco, enhanced drought tolerance. The activities of AAT and peroxidase (POD), the contents of cysteine and isoleucine, were increased, and the contents of malonaldehyde (MDA) and water loss decreased in the overexpressed plants. This study provides a novel insight into genetic response to drought stress and some candidate genes for drought tolerance genetic improvements in this plant.

Transcriptome Analysis of Stigmas of Vicia faba L. Flowers

Pollination in angiosperms depends on complex communication between pollen grains and stigmas, classified as wet or dry, depending on the presence or absence of secretions at the stigma surface, respectively. In species with wet stigma, the cuticle is disrupted and the presence of exudates is indicative of their receptivity. Most stigma studies are focused on a few species and families, many of them with self-incompatibility systems. However, there is scarce knowledge about the stigma composition in Fabaceae, the third angiosperm family, whose stigmas have been classified as semidry. Here we report the first transcriptome profiling and DEGs of Vicia faba L. styles and stigmas from autofertile (flowers able to self-fertilize in the absence of manipulation, whose exudate is released spontaneously) and autosterile (flowers that need to be manipulated to break the cuticle and release the exudates to be receptive) inbred lines. From the 76,269 contigs obtained from the de novo assembly, only 45.1% of the sequences were annotated with at least one GO term. A total of 115,920, 75,489, and 70,801 annotations were assigned to Biological Process (BP), Cellular Component (CC), and Molecular Function (MF) categories, respectively, and 5918 differentially expressed genes (DEGs) were identified between the autofertile and the autosterile lines. Among the most enriched metabolic pathways in the DEGs subset were those related with amino acid biosynthesis, terpenoid metabolism, or signal transduction. Some DEGs have been related with previous QTLs identified for autofertility traits, and their putative functions are discussed. The results derived from this work provide an important transcriptomic reference for style-stigma processes to aid our understanding of the molecular mechanisms involved in faba bean fertilization.

Thaumatin-like protein family genes VfTLP4-3 and VfTLP5 are critical for faba bean's response to drought stress at the seedling stage

Thaumatin-like proteins (TLPs) are a diverse family of pathogenesis-related proteins (PR-5) found in various plant species. Faba bean is an economically important crop known for its nutritional value and resilience to harsh environmental conditions, including drought. In this study, we conducted a comprehensive analysis of the gene structure, phylogenetics, and expression patterns of TLP genes in faba bean, with a specific focus on their response to drought stress. A total of 10 TLP genes were identified and characterized from the faba bean transcriptome, which could be classified into four distinct groups based on their evolutionary relationships. Conserved cysteine residues and REDDD motifs, which are characteristic features of TLPs, were found in most of the identified VfTLP members, and these proteins were likely to reside in the cytoplasm. Two genes, VfTLP4-3 and VfTLP5, exhibited significant upregulation under drought conditions. Additionally, ectopically expressing VfTLP4-3 and VfTLP5 in tobacco leaves resulted in enhanced drought tolerance and increased peroxidase (POD) activity. Moreover, the protein VfTLP4-3 was hypothesized to interact with glycoside hydrolase family 18 (GH18), endochitinase, dehydrin, Barwin, and aldolase, all of which are implicated in chitin metabolism. Conversely, VfTLP5 was anticipated to associate with peptidyl-prolyl cis-trans isomerase-like 3, a molecule linked to the synthesis of proline. These findings suggest that these genes may play crucial roles in mediating the drought response in faba bean through the regulation of these metabolic pathways, and serve as a foundation for future genetic improvement strategies targeting enhanced drought resilience in this economically important crop.

Novel Tri-Segmented Rhabdoviruses: A Data Mining Expedition Unveils the Cryptic Diversity of Cytorhabdoviruses

Cytorhabdoviruses (genus Cytorhabdovirus, family Rhabdoviridae) are plant-infecting viruses with enveloped, bacilliform virions. Established members of the genus Cytorhabdovirus have unsegmented single-stranded negative-sense RNA genomes (ca. 10-16 kb) which encode four to ten proteins. Here, by exploring large publicly available metatranscriptomics datasets, we report the identification and genomic characterization of 93 novel viruses with genetic and evolutionary cues of cytorhabdoviruses. Strikingly, five unprecedented viruses with tri-segmented genomes were also identified. This finding represents the first tri-segmented viruses in the family Rhabdoviridae, and they should be classified in a novel genus within this family for which we suggest the name "Trirhavirus". Interestingly, the nucleocapsid and polymerase were the only typical rhabdoviral proteins encoded by those tri-segmented viruses, whereas in three of them, a protein similar to the emaravirus (family Fimoviridae) silencing suppressor was found, while the other predicted proteins had no matches in any sequence databases. Genetic distance and evolutionary insights suggest that all these novel viruses may represent members of novel species. Phylogenetic analyses, of both novel and previously classified plant rhabdoviruses, provide compelling support for the division of the genus Cytorhabdovirus into three distinct genera. This proposed reclassification not only enhances our understanding of the evolutionary dynamics within this group of plant rhabdoviruses but also illuminates the remarkable genomic diversity they encompass. This study not only represents a significant expansion of the genomics of cytorhabdoviruses that will enable future research on the evolutionary peculiarity of this genus but also shows the plasticity in the rhabdovirus genome organization with the discovery of tri-segmented members with a unique evolutionary trajectory.

Genome-Wide Identification of bZIP Transcription Factors in Faba Bean Based on Transcriptome Analysis and Investigation of Their Function in Drought Response

Faba bean is an important cool-season edible legume crop that is constantly threatened by abiotic stresses such as drought. The basic leucine zipper (bZIP) gene family is one of the most abundant and diverse families of transcription factors in plants. It regulates plant growth and development and plays an important role in the response to biotic and abiotic stresses. In this study, we identified 18 members of the faba bean bZIP transcription factor family at the genome-wide level based on previous faba bean drought stress transcriptome sequencing data. A phylogenetic tree was constructed to group the 18 VfbZIP proteins into eight clades. Analysis of cis-acting elements in the promoter region suggested that these 18 VfbZIPs may be involved in regulating abiotic stress responses such as drought. Transcriptome data showed high expression of seven genes (VfbZIP1, VfbZIP2, VfbZIP5, VfbZIP7, VfbZIP15, VfbZIP17, and VfbZIP18) in the drought-tolerant cultivar under drought stress, in which VfbZIP1, VfbZIP2, and VfbZIP5 were consistently expressed as detected by quantitative real-time polymerase chain reaction (qRT-PCR) compared to the transcriptome data. Ectopic overexpression of the three VfbZIPs in tobacco, based on the potato Virus X (PVX) vector, revealed that VfbZIP5 enhanced the drought tolerance. Overexpressed VfbZIP5 in plants showed lower levels of proline (PRO), malondialdehyde (MDA), and peroxidase (POD) compared to those overexpressing an empty vector under 10 days of drought stress. Protein-protein interaction (PPI) analysis showed that VfbZIP5 interacted with seven proteins in faba bean, including VfbZIP7 and VfbZIP10. The results depict the importance of VfbZIPs in response to drought stress, and they would be useful for the improvement of drought tolerance.

RNA-Seq and genetic diversity analysis of faba bean (Vicia faba L.) varieties in China

Background

Faba bean (Vicia faba L) is one of the most important legumes in the world. However, there is relatively little genomic information available for this species owing to its large genome. The lack of data impedes the discovery of molecular markers and subsequent genetic research in faba bean. The objective of this study was to analyze the faba bean transcriptome, and to develop simple sequence repeat (SSR) markers to determine the genetic diversity of 226 faba bean varieties derived from different regions in China.

Methods

Faba bean varieties with different phenotype were used in transcriptome analysis. The functions of the unigenes were analyzed using various database. SSR markers were developed and the polymorphic markers were selected to conduct genetic diversity analysis.

Results

A total of 92.43 Gb of sequencing data was obtained in this study, and 133,487 unigene sequences with a total length of 178,152,541 bp were assembled. A total of 5,200 SSR markers were developed on the basis of RNA-Seq analysis. Then, 200 SSR markers were used to evaluate polymorphisms. In total, 103 (51.5%) SSR markers showed significant and repeatable bands between different faba bean varieties. Clustering analysis revealed that 226 faba bean materials were divided into five groups. Genetic diversity analysis revealed that the relationship between different faba beans in China was related, especially in the same region. These results provided a valuable data resource for annotating genes to different categories and developing SSR markers.

Recent advancement in OMICS approaches to enhance abiotic stress tolerance in legumes

The world is facing rapid climate change and a fast-growing global population. It is believed that the world population will be 9.7 billion in 2050. However, recent agriculture production is not enough to feed the current population of 7.9 billion people, which is causing a huge hunger problem. Therefore, feeding the 9.7 billion population in 2050 will be a huge target. Climate change is becoming a huge threat to global agricultural production, and it is expected to become the worst threat to it in the upcoming years. Keeping this in view, it is very important to breed climate-resilient plants. Legumes are considered an important pillar of the agriculture production system and a great source of high-quality protein, minerals, and vitamins. During the last two decades, advancements in OMICs technology revolutionized plant breeding and emerged as a crop-saving tool in wake of the climate change. Various OMICs approaches like Next-Generation sequencing (NGS), Transcriptomics, Proteomics, and Metabolomics have been used in legumes under abiotic stresses. The scientific community successfully utilized these platforms and investigated the Quantitative Trait Loci (QTL), linked markers through genome-wide association studies, and developed KASP markers that can be helpful for the marker-assisted breeding of legumes. Gene-editing techniques have been successfully proven for soybean, cowpea, chickpea, and model legumes such as Medicago truncatula and Lotus japonicus. A number of efforts have been made to perform gene editing in legumes. Moreover, the scientific community did a great job of identifying various genes involved in the metabolic pathways and utilizing the resulted information in the development of climate-resilient legume cultivars at a rapid pace. Keeping in view, this review highlights the contribution of OMICs approaches to abiotic stresses in legumes. We envisage that the presented information will be helpful for the scientific community to develop climate-resilient legume cultivars.

Combined Proteomic and Metabolomic Analysis of the Molecular Mechanism Underlying the Response to Salt Stress during Seed Germination in Barley

Salt stress is a major abiotic stress factor affecting crop production, and understanding of the response mechanisms of seed germination to salt stress can help to improve crop tolerance and yield. The differences in regulatory pathways during germination in different salt-tolerant barley seeds are not clear. Therefore, this study investigated the responses of different salt-tolerant barley seeds during germination to salt stress at the proteomic and metabolic levels. To do so, the proteomics and metabolomics of two barley seeds with different salt tolerances were comprehensively examined. Through comparative proteomic analysis, 778 differentially expressed proteins were identified, of which 335 were upregulated and 443 were downregulated. These proteins, were mainly involved in signal transduction, propanoate metabolism, phenylpropanoid biosynthesis, plant hormones and cell wall stress. In addition, a total of 187 salt-regulated metabolites were identified in this research, which were mainly related to ABC transporters, amino acid metabolism, carbohydrate metabolism and lipid metabolism; 72 were increased and 112 were decreased. Compared with salt-sensitive materials, salt-tolerant materials responded more positively to salt stress at the protein and metabolic levels. Taken together, these results suggest that salt-tolerant germplasm may enhance resilience by repairing intracellular structures, promoting lipid metabolism and increasing osmotic metabolites. These data not only provide new ideas for how seeds respond to salt stress but also provide new directions for studying the molecular mechanisms and the metabolic homeostasis of seeds in the early stages of germination under abiotic stresses.

Overexpression of DsEXLA2 gene from Dendrocalamus sinicus accelerates the plant growth rate of Arabidopsis

Expansins play crucial roles in cell wall loosening and a range of life activities involving cell wall modification. Nevertheless, the biological functions of expansin genes during fast growth of bamboo remain unclear. In this study, Dendrocalamus sinicus, the largest and fastest growing bamboo species in the world, was used as the research material, and the full length of DsEXLA2 was cloned. Bioinformatics analysis revealed that DsEXLA2 contained expansin family typical domains (DPBB_1 and Pollen_allerg_1, CDRC motif) and amino acid sequence was highly conserved among different species. The expression level of DsEXLA2 increased from top section to basal section in different internodes. Subcellular localization verified that DsEXLA2 protein was located in the cell wall. Further genetic transformation studies in Arabidopsis indicated that compared with the wild type, DsEXLA2 overexpressed transgenic plants exhibited higher plant height, thicker stem, larger leaf, and less epidermal hair number and smaller stomatal aperture in the prophase and metaphase of growth. In addition, the cellulose content in the stem of transgenic plants was increased, and cell wall was thickened significantly. Moreover, a total of 1656 differentially expressed genes (DEGs) were identified by RNA-seq. The upregulated genes were predominantly enriched in the plant-pathogen interaction, MAPK signaling pathway-plant, plant hormone signal transduction, lipid metabolism and amino acid metabolism, while the downregulated genes were mainly enriched in energy metabolism, carbohydrate metabolism, plant hormone signal transduction and ribosome. These data implied that overexpression of DsEXLA2 gene accelerates the plant growth rate of Arabidopsis. This study is helpful to reveal the molecular mechanism of DsEXLA2 in culm growth and development of D. sinicus, and to understand the rapid growth of bamboos.

iTRAQ based protein profile analysis revealed key proteins involved in regulation of drought-tolerance during seed germination in Adzuki bean

Adzuki bean is an important legume crop due to its high-quality protein, fiber, vitamins, minerals as well as rich bioactive substances. However, it is vulnerable to drought at the germination stage. However, little information is available about the genetic control of drought tolerance during seed germination in adzuki bean. In this study, some differential expression proteins (DEPs) were identified during seed germination between the drought-tolerant variety 17235 and drought-sensitive variety 17033 in adzuki bean using iTRAQ method. A total of 2834 proteins were identified in the germinating seeds of these two adzuki beans. Compared with the variety 17033, 87 and 80 DEPs were increased and decreased accumulation in variety 17235 under drought, respectively. Meanwhile, in the control group, a few DEPs, including 9 up-regulated and 21 down-regulated proteins, were detected in variety 17235, respectively. GO, KEGG, and PPI analysis revealed that the DEPs related to carbohydrate metabolism and energy production were significantly increased in response to drought stresses. To validate the proteomic function, the ectopic overexpression of V-ATPase in tobacco was performed and the result showed that V-ATPase upregulation could enhance the drought tolerance of tobacco. The results provide valuable insights into genetic response to drought stress in adzuki bean, and the DEPs could be applied to develop biomarkers related to drought tolerant in adzuki bean breeding projects.

The Adaptation and Tolerance of Major Cereals and Legumes to Important Abiotic Stresses

Abiotic stresses, including drought, extreme temperatures, salinity, and waterlogging, are the major constraints in crop production. These abiotic stresses are likely to be amplified by climate change with varying temporal and spatial dimensions across the globe. The knowledge about the effects of abiotic stressors on major cereal and legume crops is essential for effective management in unfavorable agro-ecologies. These crops are critical components of cropping systems and the daily diets of millions across the globe. Major cereals like rice, wheat, and maize are highly vulnerable to abiotic stresses, while many grain legumes are grown in abiotic stress-prone areas. Despite extensive investigations, abiotic stress tolerance in crop plants is not fully understood. Current insights into the abiotic stress responses of plants have shown the potential to improve crop tolerance to abiotic stresses. Studies aimed at stress tolerance mechanisms have resulted in the elucidation of traits associated with tolerance in plants, in addition to the molecular control of stress-responsive genes. Some of these studies have paved the way for new opportunities to address the molecular basis of stress responses in plants and identify novel traits and associated genes for the genetic improvement of crop plants. The present review examines the responses of crops under abiotic stresses in terms of changes in morphology, physiology, and biochemistry, focusing on major cereals and legume crops. It also explores emerging opportunities to accelerate our efforts to identify desired traits and genes associated with stress tolerance.

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.

Fab Advances in Fabaceae for Abiotic Stress Resilience: From 'Omics' to Artificial Intelligence

Legumes are a better source of proteins and are richer in diverse micronutrients over the nutritional profile of widely consumed cereals. However, when exposed to a diverse range of abiotic stresses, their overall productivity and quality are hugely impacted. Our limited understanding of genetic determinants and novel variants associated with the abiotic stress response in food legume crops restricts its amelioration. Therefore, it is imperative to understand different molecular approaches in food legume crops that can be utilized in crop improvement programs to minimize the economic loss. 'Omics'-based molecular breeding provides better opportunities over conventional breeding for diversifying the natural germplasm together with improving yield and quality parameters. Due to molecular advancements, the technique is now equipped with novel 'omics' approaches such as ionomics, epigenomics, fluxomics, RNomics, glycomics, glycoproteomics, phosphoproteomics, lipidomics, regulomics, and secretomics. Pan-omics-which utilizes the molecular bases of the stress response to identify genes (genomics), mRNAs (transcriptomics), proteins (proteomics), and biomolecules (metabolomics) associated with stress regulation-has been widely used for abiotic stress amelioration in food legume crops. Integration of pan-omics with novel omics approaches will fast-track legume breeding programs. Moreover, artificial intelligence (AI)-based algorithms can be utilized for simulating crop yield under changing environments, which can help in predicting the genetic gain beforehand. Application of machine learning (ML) in quantitative trait loci (QTL) mining will further help in determining the genetic determinants of abiotic stress tolerance in pulses.

Recent advances in faba bean genetic and genomic tools for crop improvement

Faba bean (Vicia faba L.), a member of the Fabaceae family, is one of the important food legumes cultivated in cool temperate regions. It holds great importance for human consumption and livestock feed because of its high protein content, dietary fibre, and nutritional value. Major faba bean breeding challenges include its mixed breeding system, unknown wild progenitor, and genome size of ~13 Gb, which is the largest among diploid field crops. The key breeding objectives in faba bean include improved resistance to biotic and abiotic stress and enhanced seed quality traits. Regarding quality traits, major progress on reduction of vicine-convicine and seed coat tannins, the main anti-nutritional factors limiting faba bean seed usage, have been recently achieved through gene discovery. Genomic resources are relatively less advanced compared with other grain legume species, but significant improvements are underway due to a recent increase in research activities. A number of bi-parental populations have been constructed and mapped for targeted traits in the last decade. Faba bean now benefits from saturated synteny-based genetic maps, along with next-generation sequencing and high-throughput genotyping technologies that are paving the way for marker-assisted selection. Developing a reference genome, and ultimately a pan-genome, will provide a foundational resource for molecular breeding. In this review, we cover the recent development and deployment of genomic tools for faba bean breeding.

Defining the Cell Wall, Cell Cycle and Chromatin Landmarks in the Responses of Brachypodium distachyon to Salinity

Excess salinity is a major stress that limits crop yields. Here, we used the model grass Brachypodium distachyon (Brachypodium) reference line Bd21 in order to define the key molecular events in the responses to salt during germination. Salt was applied either throughout the germination period ("salt stress") or only after root emergence ("salt shock"). Germination was affected at ≥100 mM and root elongation at ≥75 mM NaCl. The expression of arabinogalactan proteins (AGPs), FLA1, FLA10, FLA11, AGP20 and AGP26, which regulate cell wall expansion (especially FLA11), were mostly induced by the "salt stress" but to a lesser extent by "salt shock". Cytological assessment using two AGP epitopes, JIM8 and JIM13 indicated that "salt stress" increases the fluorescence signals in rhizodermal and exodermal cell wall. Cell division was suppressed at >75 mM NaCl. The cell cycle genes (CDKB1, CDKB2, CYCA3, CYCB1, WEE1) were induced by "salt stress" in a concentration-dependent manner but not CDKA, CYCA and CYCLIN-D4-1-RELATED. Under "salt shock", the cell cycle genes were optimally expressed at 100 mM NaCl. These changes were consistent with the cell cycle arrest, possibly at the G1 phase. The salt-induced genomic damage was linked with the oxidative events via an increased glutathione accumulation. Histone acetylation and methylation and DNA methylation were visualized by immunofluorescence. Histone H4 acetylation at lysine 5 increased strongly whereas DNA methylation decreased with the application of salt. Taken together, we suggest that salt-induced oxidative stress causes genomic damage but that it also has epigenetic effects, which might modulate the cell cycle and AGP expression gene. Based on these landmarks, we aim to encourage functional genomics studies on the responses of Brachypodium to salt.

Ethylene: A Master Regulator of Salinity Stress Tolerance in Plants

Salinity stress is one of the major threats to agricultural productivity across the globe. Research in the past three decades, therefore, has focused on analyzing the effects of salinity stress on the plants. Evidence gathered over the years supports the role of ethylene as a key regulator of salinity stress tolerance in plants. This gaseous plant hormone regulates many vital cellular processes starting from seed germination to photosynthesis for maintaining the plants' growth and yield under salinity stress. Ethylene modulates salinity stress responses largely via maintaining the homeostasis of Na+/K+, nutrients, and reactive oxygen species (ROS) by inducing antioxidant defense in addition to elevating the assimilation of nitrates and sulfates. Moreover, a cross-talk of ethylene signaling with other phytohormones has also been observed, which collectively regulate the salinity stress responses in plants. The present review provides a comprehensive update on the prospects of ethylene signaling and its cross-talk with other phytohormones to regulate salinity stress tolerance in plants.