Characterization of drought stress-responsive root transcriptome of faba bean (Vicia faba L.) using RNA sequencing

Functional Study of Amorpha fruticosa WRKY20 Gene in Response to Drought Stress

The WRKY gene family in plants regulates the plant's response to drought through regulatory networks and hormone signaling. AfWRKY20 (MT859405) was cloned from Amorpha fruticosa (A. fruticosa) seedlings using RT-PCR. The binding properties of the AfWRKY20 protein and the W-box (a DNA cis-acting element) were verified both in vivo and in vitro using EMSA and Dual-Luciferase activity assays. RT-qPCR detected that the total expression level of AfWRKY20 in leaves and roots was 22 times higher in the 30% PEG6000 simulated drought treatment compared to the untreated group. Under the simulated drought stress treatments of sorbitol and abscisic acid (ABA), the transgenic tobacco with the AfWRKY20 gene showed enhanced drought resistance at the germination stage, with significantly increased germination rate, green leaf rate, fresh weight, and root length compared to the wild-type (WT) tobacco. In addition, the superoxide dismutase (SOD) activity, chlorophyll content, and Fv/Fm ratio of AfWRKY20 transgenic tobacco were significantly higher than those of the WT tobacco under natural drought stress, while the malondialdehyde (MDA) content and 3,3'-diaminobenzidine (DAB) and nitroblue tetrazolium (NBT) staining levels were lower. The expression levels of oxidation kinase genes (NbSOD, NbPOD, and NbCAT) in transgenic tobacco under drought stress were significantly higher than those in WT tobacco. This enhancement in gene expression improved the ability of transgenic tobacco to detoxify reactive oxygen species (ROS). The survival rate of transgenic tobacco after natural drought rehydration was four times higher than that of WT tobacco. In summary, this study revealed the regulatory mechanism of AfWRKY20 in response to drought stress-induced ABA signaling, particularly in relation to ROS. This finding provides a theoretical basis for understanding the pathways of WRKY20 involved in drought stress, and offers genetic resources for molecular plant breeding aimed at enhancing drought resistance.

Identifying physiological and genetic determinants of faba bean transpiration response to evaporative demand

BACKGROUND AND AIMS

Limiting maximum transpiration rate (TR) under high vapour pressure deficit (VPD) works as a water conservation strategy. While some breeding programmes have incorporated this trait into some crops to boost yields in water-limited environments, its underlying physiological mechanisms and genetic regulation remain unknown for faba bean (Vicia faba). Thus, we aimed to identify genetic variation in the TR response to VPD in a population of faba bean recombinant inbred lines (RILs) derived from two parental lines with contrasting water use (Mélodie/2 and ILB 938/2).

METHODS

Plants were grown in well-watered soil in a climate-controlled glasshouse with diurnally fluctuating VPD and light conditions. Whole plant transpiration was measured in a gas exchange chamber that tightly regulated VPD around the shoot under constant light, while whole-plant hydraulic conductance and its components (root and stem hydraulic conductance) were calculated from dividing TR by water potential gradients measured with a pressure chamber.

KEY RESULTS

Although TR of Mélodie/2 increased linearly with VPD, ILB 938/2 limited its TR above 2.0 kPa. Nevertheless, Mélodie/2 had a higher leaf water potential than ILB 938/2 at both low (1.0 kPa) and high (3.2 kPa) VPD. Almost 90 % of the RILs limited their TR at high VPD with a break-point (BP) range of 1.5-3.0 kPa and about 10 % had a linear TR response to VPD. Thirteen genomic regions contributing to minimum and maximum transpiration, and whole-plant and root hydraulic conductance, were identified on chromosomes 1 and 3, while one locus associated with BP transpiration was identified on chromosome 5.

CONCLUSIONS

This study provides insight into the physiological and genetic control of transpiration in faba bean and opportunities for marker-assisted selection to improve its performance in water-limited environments.

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.

Resistance screening and in silico characterization of cloned novel RGA from multi-race resistant lentil germplasm against Fusarium wilt (Fusarium oxysporum f. sp. lentis)

Fusarium wilt caused by Fusarium oxysporum f. sp. lentis (Fol) is the most devastating disease of lentil present worldwide. Identification of multi-race fusarium wilt resistance genes and their incorporation into existing cultivars will help to reduce yield losses. In the present study, 100 lentil germplasms belonging to seven lentil species were screened against seven prevalent races of Fol, and accessions IC201561 (Lens culinaris subsp. culinaris), EC714243 (L. c. subsp. odemensis), and EC718238 (L. nigricans) were identified as resistant. The typical R gene codes for the nucleotide-binding site and leucine-rich repeats (NBS-LRR) at the C terminal are linked to either the Toll/interleukin 1-like receptor (TIR) or coiled coil (CC) at the N terminal. In the present study, degenerate primers, designed from the NBS region amplifying the P-loop to the GLPLA motif, isolated forty-five resistance gene analogues (RGAs) from identified resistant accessions. The sequence alignment identified both classes of RGAs, TIR and non-TIR, based on the presence of aspartate (D) and tryptophan (W) at the end of the kinase motif, respectively. The phylogenetic analysis grouped the RGAs into six classes, from LRGA1 to LRGA6, which determined the diversity of the RGAs present in the host. Grouping of the RGAs identified from Lens nigricans, LnRGA 2, 9, 13 with I2 revealed the structural similarity with the fusarium resistance gene. The similarity index ranged from 27.85% to 86.98% among the RGAs and from 26.83% to 49.41% among the known R genes, I2, Gpa2, M, and L6. The active binding sites present along the conserved motifs grouped the RGAs into 13 groups. ADP/ATP, being the potential ligand, determines the ATP binding and ATP hydrolysis activity of the RGAs. The isolated RGAs can be used to develop markers linked to the functional R gene. Furthermore, expression analysis and full-length gene isolation pave the path to identifying the molecular mechanism involved in resistance.

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.

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.

Developing a new genic SSR primer database in faba bean (Vicia faba L.)

The goal of this research was to develop a new genetic database of simple sequence repetition (SSR) primers for faba and classify them according to their target genes and respective biological processes. Approximately 75,605 and 148,196 previously published genomic and transcriptomic faba sequences, respectively, have been used to detect possible SSRs. The numbers of identified SSRs from each dataset were 25,502 and 12,319, respectively. The distribution of different repeat classes indicated that trinucleotides represent the largest number of repeat counts, followed by dinucleotides. The extracted genic SSR sequences were used to design 1091 polymerase chain reaction (PCR) primers, of which only 238 (21.8%) primers target genomic sequences and the other 853 PCR primers targeted transcriptomic sequences. The annotation of gene-targeted SSRs showed that approximately 897 genes were targeted by our SSR primers. Approximately 1890 gene ontology (GO) identification codes have been obtained. The GO keywords were distributed among distinct molecular cell features. The highest redundancies involved 554 technical words, 196 domains, and 160 molecular feature phrases. These GO codes belonged to the general level of GO and included molecular function, cellular component, and biological process (544, 670, and 676 GOs, respectively). Twenty-seven SSR PCR primers were synthesized to 12 Egyptian faba bean genotypes. Approximately 11 SSR provided one to two PCR bands, whereas other SSRs provided only one sharp band with polymorphic band size. There were 13 polymorphic primers. The polymorphism information content was 0.3, which implied moderate informativeness.

The influence of soil drought stress on the leaf transcriptome of faba bean (Vicia faba L.) in the Qinghai-Tibet Plateau

Water deficit has a significant impact on growth, development and yield of fava bean (Vicia fava L.) in arid and semi-arid climates. The aim of this study was to identify differentially expressed genes in the Qinghai 13 genotype under soil drought through leaf transcriptome analysis. A total of 256.95 M clean reads were obtained and assembled into 176334 unigenes, with an average length of 766 bp. A total of 9126 (4439 upregulated and 4687 downregulated) differentially expressed genes (DEGs) were identified in faba bean leaves under soil drought. In total, 324 putative transcription factors were identified and classified as belonging to different transcription factor families. According to GO and KEGG analysis, the soil drought stress-inducible DEGs encoded proteins mainly involved in regulating photosynthesis, osmotic adjustment, detoxification, autophagy and other functions. In addition, a large portion of DEGs appeared to be novel because they could not be annotated in any functional databases, therefore, suggesting a specific response to soil drought in faba bean. Finally, RNA-seq analysis was validated by quantitative reverse-transcription PCR analysis. This work provides comprehensive and valuable information for understanding the molecular mechanisms which faba bean uses to respond to soil drought.

Characterization of Drought-Responsive Transcriptome During Seed Germination in Adzuki Bean (Vigna angularis L.) by PacBio SMRT and Illumina Sequencing

The full-length single-molecular sequencing and short reads Illumina sequencing were combined to generate the transcripts of adzuki bean with high-quality. A total of 17,636 loci and 60,454 transcripts were detected in this study. To characterize the drought-responsive genes during seed germination in adzuki bean, two varieties, i.e., tolerant and sensitive to drought stress, were selected to conduct analysis of alternative splicing dynamics (AS) and differentially expressed genes (DEGs) by combining the newly assembled draft genome and public adzuki bean reference genome. AS analysis indicated that both the two varieties underwent a little more AS events under control conditions than under drought stress. Among the AS events, IR (intron retention) predominately accounted for 34.3%, whereas AD (alternative donor site) was the least frequent with 15.8%. Meanwhile, 562 long non-coding RNAs, 409 fusion genes and 1208 transcription factors were identified. Moreover, a total of 5,337 DEGs were identified in comparison of the two varieties with drought or control treatments. Notably, 82 DEGs were discovered in the two varieties under drought stress, which might be the candidate in regulation of seed germination to answer for different drought tolerance. The DEGs encoded proteins involved in primary or second metabolism, plant hormone signal transduction, transcript or translation processes, ubiquitin proteasome system, transcription factor, transporters, and so on. The results facilitate to increase the knowledge about the mechanism of drought tolerance during crop seed germination, and provide reference for the breeding of drought-tolerant adzuki bean.

Legume genomics and transcriptomics: From classic breeding to modern technologies

Legumes are essential and play a significant role in maintaining food standards and augmenting physiochemical soil properties through the biological nitrogen fixation process. Biotic and abiotic factors are the main factors limiting legume production. Classical breeding methodologies have been explored extensively about the problem of truncated yield in legumes but have not succeeded at the desired rate. Conventional breeding improved legume genotypes but with more resources and time. Recently, the invention of next-generation sequencing (NGS) and high-throughput methods for genotyping have opened new avenues for research and developments in legume studies. During the last decade, genome sequencing for many legume crops documented. Sequencing and re-sequencing of important legume species have made structural variation and functional genomics conceivable. NGS and other molecular techniques such as the development of markers; genotyping; high density genetic linkage maps; quantitative trait loci (QTLs) identification, expressed sequence tags (ESTs), single nucleotide polymorphisms (SNPs); and transcription factors incorporated into existing breeding technologies have made possible the accurate and accelerated delivery of information for researchers. The application of genome sequencing, RNA sequencing (transcriptome sequencing), and DNA sequencing (re-sequencing) provide considerable insights for legume development and improvement programs.

Moreover, RNA-Seq helps to characterize genes, including differentially expressed genes, and can be applied for functional genomics studies, especially when there is limited information available for the studied genomes. Genome-based crop development studies and the availability of genomics data as well as decision-making gears look be specific for breeding programs. This review mainly presents an overview of the path from classical breeding to new emerging genomics tools, which will trigger and accelerate genomics-assisted breeding for recognition of novel genes for yield and quality characters for sustainable legume crop production.

Transcriptome-based analysis of carotenoid accumulation-related gene expression in petals of Chinese cabbage (Brassica rapa L.)

To identify genes associated with carotenoid accumulation in petals of Chinese cabbage, the composition and content of carotenoids were analyzed, and comparative transcriptome sequencing was performed between the yellow flower line, 92S105, and the orange flower line, 94C9. High-performance liquid chromatography (HPLC) revealed that petals of 92S105 were high in violaxanthin as well as lutein, whereas petals of 94C9 showed considerable levels of lutein and β-carotene. Transcriptome analysis showed that 3534 and 3833 genes were up- and down-regulated in 94C9, respectively. Among these differentially expressed genes (DEGs), many related to carotenoid accumulation were identified, including 12 carotenoid biosynthesis pathway genes, 4 transcription factor genes, and 1028 specifically expressed genes. β-carotene hydroxylase 1 (BrBCH1), BrBCH2, zeaxanthin epoxidase (BrZEP), and MYB transcription factor gene (BrGAMYB) were down-regulated in petals of 94C9 when compared with petals of 92S105, which caused β-carotene accumulation and may lead to orange petal color in 94C9. Expression levels of 20 DEGs were verified by qPCR and the results were highly consistent with those of transcriptome sequencing. Moreover, Gene Ontology (GO) enrichment analysis revealed that membrane, binding, and metabolic processes were the most significantly enriched GO terms in cellular component, molecular function, and biological process ontologies, respectively. In conclusion, our study analyzed the differences in composition and content of carotenoids between 92S105 and 94C9 and identified potential candidate genes related to carotenoid accumulation in petals, thereby creating a solid foundation for future studies on the mechanism regulating carotenoid accumulation in petals of Chinese cabbage.