Transcriptome Analyses Reveal Candidate Pod Shattering-Associated Genes Involved in the Pod Ventral Sutures of Common Vetch (Vicia sativa L.)

Transcriptomic insights into the potential impacts of flavonoids and nodule-specific cysteine-rich peptides on nitrogen fixation in Vicia villosa and Vicia sativa

Vicia villosa (VV) and Vicia sativa (VS) are legume forages highly valued for their excellent nitrogen fixation. However, no research has addressed the mechanisms underlying their differences in nitrogen fixation. This study employed physiological, cytological, and comparative transcriptomic approaches to elucidate the disparities in nitrogen fixation between them. Our results showed that the total amount of nitrogen fixed was 60.45% greater in VV than in VS, and the comprehensive nitrogen response performance was 94.19% greater, while the nitrogen fixation efficiency was the same. The infection zone and differentiated bacteroid proportion in mature VV root nodules were 33.76% and 19.35% greater, respectively, than those in VS. The size of the VV genome was 15.16% larger than that of the VS genome, consistent with its greater biomass. A significant enrichment of the flavonoid biosynthetic pathway was found only for VV-specific genes, among which chalcone-flavonone isomerase, caffeoyl-CoA-O-methyltransferase and stilbene synthase were extremely highly expressed. The VV-specific genes also exhibited significant enrichment in symbiotic nodulation; genes related to nodule-specific cysteine-rich peptides (NCRs) comprised 61.11% of the highly expressed genes. qRT‒PCR demonstrated that greater enrichment and expression of the dominant NCR (Unigene0004451) were associated with greater nodule bacteroid differentiation and greater nitrogen fixation in VV. Our findings suggest that the greater total nitrogen fixation of VV was attributed to its larger biomass, leading to a greater nitrogen demand and enhanced fixation physiology. This process is likely achieved by the synergistic effects of high bacteroid differentiation along with high expression of flavonoid and NCR genes.

Genome-wide identification of GH28 family and insight into its contributions to pod shattering resistance in Brassica napus L

Rapeseed (Brassica napus L.), accounts for nearly 16% of vegetable oil, is the world's second produced oilseed. However, pod shattering has caused significant yield loses in rapeseed production, particularly during mechanical harvesting. The GH28 genes can promote pod shattering by changing the structure of the pod cell wall in Arabidopsis. However, the role of the GH28 gene family in rapeseed was largely unknown. Therefore, a genome-wide comprehensive analysis was conducted to classify the role of GH28 gene family on rapeseed pod shattering. A total of 37 BnaGH28 genes in the rapeseed genome were identified. These BnaGH28s can be divided into five groups (Group A-E), based on phylogenetic and synteny analysis. Protein property, gene structure, conserved motif, cis-acting element, and gene expression profile of BnaGH28 genes in the same group were similar. Specially, the expression level of genes in group A-D was gradually decreased, but increased in group E with the development of silique. Among eleven higher expressed genes in group E, two BnaGH28 genes (BnaA07T0199500ZS and BnaC06T0206500ZS) were significantly regulated by IAA or GA treatment. And the significant effects of BnaA07T0199500ZS variation on pod shattering resistance were also demonstrated in present study. These results could open a new window for insight into the role of BnaGH28 genes on pod shattering resistance in rapeseed.

Transcriptome Profiling Reveals Molecular Responses to Salt Stress in Common Vetch (Vicia sativa L.)

Common vetch (Vicia sativa L.) is an important annual diploid leguminous forage. In the present study, transcriptomic profiling in common vetch in response to salt stress was conducted using a salt-tolerant line (460) and a salt-sensitive line (429). The common responses in common vetch and the specific responses associated with salt tolerance in 460 were analyzed. Several KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, including plant hormone and MAPK (mitogen-activated protein kinase) signaling, galactose metabolism, and phenylpropanoid phenylpropane biosynthesis, were enriched in both lines, though some differentially expressed genes (DEGs) showed distinct expression patterns. The roots in 460 showed higher levels of lignin than in 429. α-linolenic acid metabolism, carotenoid biosynthesis, the photosynthesis-antenna pathway, and starch and sucrose metabolism pathways were specifically enriched in salt-tolerant line 460, with higher levels of accumulated soluble sugars in the leaves. In addition, higher transcript levels of genes involved in ion homeostasis and reactive oxygen species (ROS) scavenging were observed in 460 than in 429 in response to salt stress. The transcriptomic analysis in common vetch in response to salt stress provides useful clues for further investigations on salt tolerance mechanism in the future.

Transcriptomic analysis revealed the candidate metabolic pathways and genes associated with cold tolerance in a mutant without anthocyanin accumulation in common vetch (Vicia sativa L.)

Common vetch (Vicia sativa L.) is a leguminous crop used to feed livestock with vegetative organs or fertilize soils by returning to the field. Survival of fall-seeded plants is often affected by freezing damage during overwintering. This study aims to investigate the transcriptomic profiling in response to cold in a mutant with reduced accumulation of anthocyanins under normal growth and low-temperature conditions for understanding the underlying mechanisms. The mutant had increased cold a tolerance with higher survival rate and biomass during overwintering compared to the wild type, which led to increased forage production. Transcriptomic analysis in combination with qRT-PCR and physiological measurements revealed that reduced anthocyanins accumulation in the mutant resulted from reduced expression of serial genes involving in anthocyanin biosynthesis, which led to the altered metabolism, with an increased accumulation of free amino acids and polyamines. The higher levels of free amino acids and proline in the mutant under low temperature were associated with improved cold tolerance. The altered expression of some genes involved in ABA and GA signaling was also associated with increased cold tolerance in the mutant.

Analysis of seed production and seed shattering in a new artificial grassland forage: pigeon pea

Pigeon pea is a perennial leguminous plant that is widely cultivated as a forage and pharmaceutical plant in subtropical and tropical areas, especially in artificial grasslands. Higher seed shattering is one of the most important factors in potentially increasing the seed yield of pigeon pea. Advance technology is necessary to increase the seed yield of pigeon pea. Through 2 consecutive years of field observations, we found that fertile tiller number was the key component of the seed yield of pigeon pea due to the direct effect of fertile tiller number per plant (0.364) on pigeon pea seed yield was the highest. Multiplex morphology, histology, and cytological and hydrolytic enzyme activity analysis showed that shatter-susceptible and shatter-resistant pigeon peas possessed an abscission layer at the same time (10 DAF); however, abscission layer cells dissolved earlier in shattering-susceptible pigeon pea (15 DAF), which led to the tearing of the abscission layer. The number of vascular bundle cells and vascular bundle area were the most significant negative factors (p< 0.01) affecting seed shattering. Cellulase and polygalacturonase were involved in the dehiscence process. In addition, we inferred that larger vascular bundle tissues and cells in the ventral suture of seed pods could effectively resist the dehiscence pressure of the abscission layer. This study provides foundation for further molecular studies to increase the seed yield of pigeon pea.

Transcript profiling of hairy vetch (Vicia villosa Roth) identified interesting genes for seed dormancy

Hairy vetch, a diploid annual legume species, has a robust growth habit, high biomass yield, and winter hardy characteristics. Seed hardness is a major constraint for growing hairy vetch commercially. Hard seeded cultivars are valuable as forages, whereas soft seeded and shatter resistant cultivars have advantages for their use as a cover crop. Transcript analysis of hairy vetch was performed to understand the genetic mechanisms associated with important hairy vetch traits. RNA was extracted from leaves, flowers, immature pods, seed coats, and cotyledons of contrasting soft and hard seeded "AU Merit" plants. A range of 31.22-79.18 Gb RNA sequence data per tissue sample were generated with estimated coverage of 1040-2639×. RNA sequence assembly and mapping of the contigs against the Medicago truncatula (V4.0) genome identified 76,422 gene transcripts. A total of 24,254 transcripts were constitutively expressed in hairy vetch tissues. Key genes, such as KNOX4 (a class II KNOTTED-like homeobox KNOXII gene), qHs1 (endo-1,4-β-glucanase), GmHs1-1 (calcineurin-like metallophosphoesterase), chitinase, shatterproof 1 and 2 (SHP1, SHP2), shatter resistant 1-5 (SHAT1-5)(NAC transcription factor), PDH1 (prephenate dehydrogenase 1), and pectin methylesterases with a potential role in seed hardness and pod shattering, were further explored based on genes involved in seed hardness from other species to query the hairy vetch transcriptome data. Identification of interesting candidate genes in hairy vetch can facilitate the development of improved cultivars with desirable seed characteristics for use as a forage and as a cover crop.

Designing Novel Strategies for Improving Old Legumes: An Overview from Common Vetch

Common vetch (Vicia sativa L.) is a grain legume used in animal feeding, rich in protein content, fatty acid, and mineral composition that makes for a very adequate component to enrich feedstuff. In addition, relevant pharmacological properties have been reported in humans. The common vetch, similar to other legumes, can fix atmospheric nitrogen, a crucial feature for sustainable agricultural systems. These properties enhance the use of vetch as a cover crop and its sowing in intercropping systems. Moreover, several studies have recently pointed out the potential of vetch in the phytoremediation of contaminated soils. These characteristics make vetch a relevant crop, which different potential improvements target. Varieties with different yields, flowering times, shattering resistance, nutritional composition, rhizobacteria associations, drought tolerance, nitrogen fixation capacity, and other agronomic-relevant traits have been identified when different vetch accessions are compared. Recently, the analysis of genomic and transcriptomic data has allowed the development of different molecular markers to be used for assisted breeding purposes, promoting crop improvement. Here, we review the potential of using the variability of V. sativa genetic resources and new biotechnological and molecular tools for selecting varieties with improved traits to be used in sustainable agriculture systems.

Screening and identification of multiple abiotic stress responsive candidate genes based on hybrid-sequencing in Vicia sativa

Common vetch is an important leguminous forage for both livestock fodder and green manure and has a tremendous latent capacity in a sustainable agroecosystem. In the present study, a comprehensive transcriptome analysis of the aboveground leaves and underground roots of common vetch under multiple abiotic stress treatments, including NaCl, drought, cold, and cold drought, was performed using hybrid-sequencing technology, i. e. single-molecule real-time sequencing technology (SMRT) and supplemented by next-generation sequencing (NGS) technology. A total of 485,038 reads of insert (ROIs) with a mean length of 2606 bp and 228,261 full-length nonchimeric (FLNC) reads were generated. After deduplication, 39,709 transcripts were generated. Of these transcripts, we identified 1059 alternative splicing (AS) events, 17,227 simple sequence repeats (SSRs), and 1647 putative transcription factors (TFs). Furthermore, 640 candidates long noncoding RNAs (lncRNAs) and 28,256 complete coding sequences (CDSs) were identified. In gene annotation analyses, a total of 38,826 transcripts (97.78%) were annotated in eight public databases. Finally, seven multiple abiotic stress-responsive candidate genes were obtained through gene expression, annotation information, and protein-protein interaction (PPI) networks. Our research not only enriched the structural information of FL transcripts in common vetch, but also provided useful information for exploring the molecular mechanism of multiple abiotic stress tolerance between aboveground and underground tissues in common vetch and related legumes.

A comparative transcriptomic analysis reveals a coordinated mechanism activated in response to cold acclimation in common vetch (Vicia sativa L.)

BACKGROUND

Due to its strong abiotic stress tolerance, common vetch is widely cultivated as a green manure and forage crop in grass and crop rotation systems. The comprehensive molecular mechanisms activated in common vetch during cold adaptation remain unknown.

RESULTS

We investigated physiological responses and transcriptome profiles of cold-sensitive (Lanjian No. 1) and cold-tolerant (Lanjian No. 3) cultivars during cold acclimation to explore the molecular mechanisms of cold acclimation. In total, 2681 and 2352 differentially expressed genes (DEGs) were identified in Lanjian No. 1 and Lanjian No. 3, respectively; 7532 DEGs were identified in both lines. DEGs involved in "plant hormone signal transduction" were significantly enriched during cold treatment, and 115 DEGs involved in cold-processed hormone signal transduction were identified. Common vetch increased the level of indoleacetic acid (IAA) by upregulating the transcriptional regulator Aux/IAA and downregulating GH3, endowing it with stronger cold tolerance. An auxin-related DEG was overexpressed in yeast and shown to possess a biological function conferring cold tolerance.

CONCLUSION

This study identifies specific genes involved in Ca²⁺ signaling, redox regulation, circadian clock, plant hormones, and transcription factors whose transcriptional differentiation during cold acclimation may improve cold tolerance and contributes to the understanding of common and unique molecular mechanisms of cold acclimation in common vetch. The candidate genes identified here also provide valuable resources for further functional genomic and breeding studies.

Novel candidate loci for morpho-agronomic and seed quality traits detected by targeted genotyping-by-sequencing in common bean

Phaseolus vulgaris L., known as common bean, is one of the most important grain legumes cultivated around the world for its immature pods and dry seeds, which are rich in protein and micronutrients. Common bean offers a cheap food and protein sources to ameliorate food shortage and malnutrition around the world. However, the genetic basis of most important traits in common bean remains unknown. This study aimed at identifying QTL and candidate gene models underlying twenty-six agronomically important traits in common bean. For this, we assembled and phenotyped a diversity panel of 200 P. vulgaris genotypes in the greenhouse, comprising determinate bushy, determinate climbing and indeterminate climbing beans. The panel included dry beans and snap beans from different breeding programmes, elite lines and landraces from around the world with a major focus on accessions of African, European and South American origin. The panel was genotyped using a cost-conscious targeted genotyping-by-sequencing (GBS) platform to take advantage of highly polymorphic SNPs detected in previous studies and in diverse germplasm. The detected single nucleotide polymorphisms (SNPs) were applied in marker-trait analysis and revealed sixty-two quantitative trait loci (QTL) significantly associated with sixteen traits. Gene model identification via a similarity-based approach implicated major candidate gene models underlying the QTL associated with ten traits including, flowering, yield, seed quality, pod and seed characteristics. Our study revealed six QTL for pod shattering including three new QTL potentially useful for breeding. However, the panel was evaluated in a single greenhouse environment and the findings should be corroborated by evaluations across different field environments. Some of the detected QTL and a number of candidate gene models only elucidate the understanding of the genetic nature of these traits and provide the basis for further studies. Finally, the study showed the possibility of using a limited number of SNPs in performing marker-trait association in common bean by applying a highly scalable targeted GBS approach. This targeted GBS approach is a cost-efficient strategy for assessment of the genetic basis of complex traits and can enable geneticists and breeders to identify novel loci and targets for marker-assisted breeding more efficiently.

Evaluation of salt tolerance in common vetch (Vicia sativa L.) germplasms and the physiological responses to salt stress

Common vetch (Vicia sativa L.) is an important leguminous crop, providing humans with starch from seeds, feeding livestock with vegetative organs, or fertilizing soils by returning to field. It is aimed to evaluate salt tolerance in common vetch collections for breeding programs and to investigate the underlined physiological mechanisms. Relative germination rate and relative seedling growth showed great difference among common vetch collections in response to salt. A lower level of Na⁺ and higher levels of K⁺ and K⁺/Na⁺ ratio were maintained in both shoots and roots in salt-tolerant collections than in salt-sensitive ones under salt stress. Expression of the genes involved in transportation and redistribution of Na⁺ and K⁺ were cooperatively responsible for salt stress. Transcript levels of NHX7, HKT1, AKT2, and HAK17 in leaves and roots were induced after salt stress, with higher transcript levels in salt-tolerant collections compared with the sensitive ones. Proline and P5CS1 transcript levels were increased after salt stress, with higher levels in salt-tolerant collection compared with salt-sensitive ones. Both O₂^- and H₂O₂ were accumulated after salt stress, and lower levels were accumulated in salt-tolerant collection compared with salt-sensitive ones. Superoxide dismutase, catalase and ascorbate peroxidase activities were altered in response to salt and higher levels were maintained in salt-tolerant collections compared with salt-sensitive ones. It is suggested that salt tolerance in common vetch is associated with maintenance of K⁺ and Na⁺ homeostasis and the associated gene expression and promoted proline accumulation and antioxidant defense system.

Histological characteristics, cell wall hydrolytic enzyme activity, and transcriptome analysis with seed shattering of Stylosanthes accessions

Stylosanthes spp. (stylo) are annual or perennial legume forages that are widely grown as forage and cover crops in tropical and subtropical regions. However, the seed yield of stylo is very low due to serious seed shattering. In the present study, we found that, although seed shattering was common among the stylo accessions, the shattering rates were genetically different. Therefore, we first synthesized the morphological, histological characteristic, physiochemical, and transcriptome analyses to determine the seed shattering mechanism in stylo. In general, the stylo germplasm with shorter lobules and thicker stems had a lower seed shattering rate and a higher seed weight. The seed and seed stalk joint is the abscission zone in stylo. Multiplex histology and hydrolytic enzyme activity analysis showed that the tearing of the abscission zone occurs due to the intense enzymatic degradation of polygalacturonase and cellulase in the seed shattering-susceptible accession TF0275. cDNA libraries from the abscission zone tissues of TF0041 and TF0275 at 14, 21, and 28 days after flowering were constructed and sequenced. A total of 47,606 unigenes were annotated and 18,606 differentially expressed genes (DEGs) were detected, including 9,140 upregulated and 9,446 downregulated unigenes. Furthermore, the 26 candidate DEGs involved in lignin biosynthesis, cellulase synthesis, and plant hormone signal transduction were found at all three developmental stages. This study provides valuable insights for future mechanistic studies of seed shattering in stylo.

Genome-Wide Identification of ERF Transcription Factor Family and Functional Analysis of the Drought Stress-Responsive Genes in Melilotus albus

As an important forage legume with high values in feed and medicine, Melilotus albus has been widely cultivated. The AP2/ERF transcription factor has been shown to play an important regulatory role in plant drought resistance, but it has not been reported in the legume forage crop M. albus. To digger the genes of M. albus in response to drought stress, we identified and analyzed the ERF gene family of M. albus at the genome-wide level. A total of 100 MaERF genes containing a single AP2 domain sequence were identified in this study, named MaERF001 to MaERF100, and bioinformatics analysis was performed. Collinearity analysis indicated that segmental duplication may play a key role in the expansion of the M. albus ERF gene family. Cis-acting element predictions suggest that MaERF genes are involved in various hormonal responses and abiotic stresses. The expression patterns indicated that MaERFs responded to drought stress to varying degrees. Furthermore, four up-regulated ERFs (MaERF008, MaERF037, MaERF054 and MaERF058) under drought stress were overexpressed in yeast and indicated their biological functions to confer the tolerance to drought. This work will advance the understanding of the molecular mechanisms underlying the drought response in M. albus. Further study of the promising potential candidate genes identified in this study will provide a valuable resource as the next step in functional genomics studies and improve the possibility of improving drought tolerance in M. albus by transgenic approaches.

Mechanism of pod shattering in the forage legume Medicago ruthenica

Pod shattering is a seed dispersal strategy and an important agronomical trait in domesticated crops. The relationship between pod shattering and pod morphology in the genus Medicago is well known; however, the detailed mechanism underlying pod dehiscence in Medicago ruthenica, a perennial legume used for forage production, is unknown. Here, the pod ventral sutures of shatter-resistant and shatter-susceptible M. ruthenica genotypes were examined at 8, 12, 16, and 20 d after flowering. The mechanism of pod shattering was analyzed through microscopic observations, polygalacturonase (PG) and cellulase (CE) activity analyses, and RNA-sequencing (RNA-Seq), and the results were verified via reverse transcriptase-quantitative polymerase chain reaction. Pod shattering at the ventral suture in M. ruthenica occurs via a combination of two mechanisms: degradation of the middle lamella at the abscission layers (ALs) and detachment of lignified cells on either side of the ALs triggered by physical forces. Increased PG and CE activities in the pod ventral suture are essential for AL cell-autolysis in the shatter-susceptible genotype. RNA-Seq revealed that 11 genes encoding PG and CE were highly expressed in the ventral sutures of the shatter-susceptible genotype. The expression levels of auxin biosynthesis-related genes decreased in the AL cells and they were negatively associated with pod dehiscence. These results enhance our understanding of the pod shattering mechanism not only in M. ruthenica but also in other leguminous plants.

Unraveling Origin, History, Genetics, and Strategies for Accelerated Domestication and Diversification of Food Legumes

Domestication is a dynamic and ongoing process of transforming wild species into cultivated species by selecting desirable agricultural plant features to meet human needs such as taste, yield, storage, and cultivation practices. Human plant domestication began in the Fertile Crescent around 12,000 years ago and spread throughout the world, including China, Mesoamerica, the Andes and Near Oceania, Sub-Saharan Africa, and eastern North America. Indus valley civilizations have played a great role in the domestication of grain legumes. Crops, such as pigeon pea, black gram, green gram, lablab bean, moth bean, and horse gram, originated in the Indian subcontinent, and Neolithic archaeological records indicate that these crops were first domesticated by early civilizations in the region. The domestication and evolution of wild ancestors into today’s elite cultivars are important contributors to global food supply and agricultural crop improvement. In addition, food legumes contribute to food security by protecting human health and minimize climate change impacts. During the domestication process, legume crop species have undergone a severe genetic diversity loss, and only a very narrow range of variability is retained in the cultivars. Further reduction in genetic diversity occurred during seed dispersal and movement across the continents. In general, only a few traits, such as shattering resistance, seed dormancy loss, stem growth behavior, flowering–maturity period, and yield traits, have prominence in the domestication process across the species. Thus, identification and knowledge of domestication responsive loci were often useful in accelerating new species’ domestication. The genes and metabolic pathways responsible for the significant alterations that occurred as an outcome of domestication might aid in the quick domestication of novel crops. Further, recent advances in “omics” sciences, gene-editing technologies, and functional analysis will accelerate the domestication and crop improvement of new crop species without losing much genetic diversity. In this review, we have discussed about the origin, center of diversity, and seed movement of major food legumes, which will be useful in the exploration and utilization of genetic diversity in crop improvement. Further, we have discussed about the major genes/QTLs associated with the domestication syndrome in pulse crops and the future strategies to improve the food legume crops.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

Development and Validation of SNP and InDel Markers for Pod-Shattering Tolerance in Soybean

Pod-shattering causes a significant yield loss in many soybean cultivars. Shattering-tolerant cultivars provide the most effective approach to minimizing this loss. We developed molecular markers for pod-shattering and validated them in soybeans with diverse genetic backgrounds. The genes Glyma.16g141200, Glyma.16g141500, and Glyma.16g076600, identified in our previous study by quantitative trait locus (QTL) mapping and whole-genome resequencing, were selected for marker development. The whole-genome resequencing of three parental lines (one shattering-tolerant and two shattering-susceptible) identified single nucleotide polymorphism (SNP) and/or insertion/deletion (InDel) regions within or near the selected genes. Two SNPs and one InDel were converted to Kompetitive Allele-Specific PCR (KASP) and InDel markers, respectively. The accuracy of the markers was examined in the two recombinant inbred line populations used for the QTL mapping, as well as the 120 varieties and elite lines, through allelic discrimination and phenotyping by the oven-drying method. Both types of markers successfully discriminated the pod shattering-tolerant and shattering-susceptible genotypes. The prediction accuracy, which was as high as 90.9% for the RILs and was 100% for the varieties and elite lines, also supported the accuracy and usefulness of these markers. Thus, the markers can be used effectively for genetic and genomic studies and the marker-assisted selection for pod-shattering tolerance in soybean.

Transcriptome and Coexpression Network Analyses Provide In-Sights into the Molecular Mechanisms of Hydrogen Cyanide Synthesis during Seed Development in Common Vetch (Vicia sativa L.)

The common vetch (Vicia sativa L.) seed is an ideal plant-based protein food for humans, but its edible value is mainly limited by the presence of cyanogenic glycosides that hydrolyze to produce toxic hydrogen cyanide (HCN), and the genes that regulate HCN synthesis in common vetch are unknown. In this study, seeds from common vetch at 5, 10, 15, 20, 25, 30, and 35 days after anthesis were sampled, and the seven stages were further divided into five developmental stages, S1, S2, S3, S4, and S5, based on morphological and transcriptome analyses. A total of 16,403 differentially expressed genes were identified in the five developmental stages. The HCN contents of seeds in these five stages were determined by alkaline titration, and weighted gene coexpression network analysis was used to explain the molecular regulatory mechanism of HCN synthesis in common vetch seeds. Eighteen key regulatory genes for HCN synthesis were identified, including the VsGT2, VsGT17 and CYP71A genes, as well as the VsGT1 gene family. VsGT1, VsGT2, VsGT17 and CYP71A jointly promoted HCN synthesis, from 5 to 25 days after anthesis, with VsGT1-1, VsGT1-4, VsGT1-11 and VsGT1-14 playing major roles. The HCN synthesis was mainly regulated by VsGT1, from 25 to 35 days after anthesis. As the expression level of VsGT1 decreased, the HCN content no longer increased. In-depth elucidation of seed HCN synthesis lays the foundations for breeding common vetch with low HCN content.

Chromosome-level assembly of the common vetch (Vicia sativa) reference genome

Vicia sativa L. (common vetch, n = 6) is an annual, herbaceous, climbing legume, originating in the Fertile Crescent of the Middle East and now widespread in the Mediterranean basin, West, Central and Eastern Asia, North and South America. V. sativa is of economic importance as a forage legume in countries such as Australia, China, and the USA, and contributes valuable nitrogen to agricultural rotation cropping systems. To accelerate precision genome breeding and genomics-based selection of this legume, we present a chromosome-level reference genome sequence for V. sativa, constructed using a combination of long-read Oxford Nanopore sequencing, short-read Illumina sequencing, and high-throughput chromosome conformation data (CHiCAGO and Hi-C) analysis. The chromosome-level assembly of six pseudo-chromosomes has a total genome length of 1.65 Gbp, with a median contig length of 684 Kbp. BUSCO analysis of the assembly demonstrated very high completeness of 98% of the dicotyledonous orthologs. RNA-seq analysis and gene modelling enabled the annotation of 53,218 protein-coding genes. This V. sativa assembly will provide insights into vetch genome evolution and be a valuable resource for genomic breeding, genetic diversity and for understanding adaption to diverse arid environments.

An Efficient Root Transformation System for Recalcitrant Vicia sativa

Common vetch (Vicia sativa) is a multi-purpose legume widely used in pasture and crop rotation systems. Vetch seeds have desirable nutritional characteristics and are often used to feed ruminant animals. Although transcriptomes are available for vetch, problems with genetic transformation and plant regeneration hinder functional gene studies in this legume species. Therefore, the aim of this study was to develop a simple, efficient and rapid hairy root transformation system for common vetch to facilitate functional gene analysis. At first, we infected the hypocotyls of 5-day-old in vitro or in vivo, soil-grown seedlings with Rhizobium rhizogenes K599 using a stabbing method and produced transgenic hairy roots after 24 days at 19 and 50% efficiency, respectively. We later improved the hairy root transformation in vitro by infecting different explants (seedling, hypocotyl-epicotyl, and shoot) with R. rhizogenes. We observed hairy root formation at the highest efficiency in shoot and hypocotyl-epicotyl explants with 100 and 93% efficiency, respectively. In both cases, an average of four hairy roots per explant were obtained, and about 73 and 91% of hairy roots from shoot and hypocotyl-epicotyl, respectively, showed stable expression of a co-transformed marker β-glucuronidase (GUS). In summary, we developed a rapid, highly efficient, hairy root transformation method by using R. rhizogenes on vetch explants, which could facilitate functional gene analysis in common vetch.

Transcriptome Analysis of Salt Stress in Hibiscus hamabo Sieb. et Zucc Based on Pacbio Full-Length Transcriptome Sequencing

Hibiscus hamabo Sieb. et Zucc is an important semi-mangrove plant with great morphological features and strong salt resistance. In this study, by combining single molecule real time and next-generation sequencing technologies, we explored the transcriptomic changes in the roots of salt stressed H. hamabo. A total of 94,562 unigenes were obtained by clustering the same isoforms using the PacBio RSII platform, and 2269 differentially expressed genes were obtained under salt stress using the Illumina platform. There were 519 differentially expressed genes co-expressed at each treatment time point under salt stress, and these genes were found to be enriched in ion signal transduction and plant hormone signal transduction. We used Arabidopsis thaliana (L.) Heynh. transformation to confirm the function of the HhWRKY79 gene and discovered that overexpression enhanced salt tolerance. The full-length transcripts generated in this study provide a full characterization of the transcriptome of H. hamabo and may be useful in mining new salt stress-related genes specific to this species, while facilitating the understanding of the salt tolerance mechanisms.

Genome survey sequencing of common vetch (Vicia sativa L.) and genetic diversity analysis of Chinese germplasm with genomic SSR markers

Background

Common vetch (Vicia sativa L.) is an annual legume with excellent suitability in cold and dry regions. Despite its great applied potential, the genomic information regarding common vetch currently remains unavailable.

Methods and results

In the present study, the whole genome survey of common vetch was performed using the next-generation sequencing (NGS). A total of 79.84 Gbp high quality sequence data were obtained and assembled into 3,754,145 scaffolds with an N50 length of 3556 bp. According to the K-mer analyses, the genome size, heterozygosity rate and GC content of common vetch genome were estimated to be 1568 Mbp, 0.4345 and 35%, respectively. In addition, a total of 76,810 putative simple sequence repeats (SSRs) were identified. Among them, dinucleotide was the most abundant SSR type (44.94%), followed by Tri- (35.82%), Tetra- (13.22%), Penta- (4.47%) and Hexanucleotide (1.54%). Furthermore, a total of 58,175 SSR primer pairs were designed and ten of them were validated in Chinese common vetch. Further analysis showed that Chinese common vetch harbored high genetic diversity and could be clustered into two main subgroups.

Conclusion

This is the first report about the genome features of common vetch, and the information will help to design whole genome sequencing strategies. The newly identified SSRs in this study provide basic molecular markers for germplasm characterization, genetic diversity and QTL mapping studies for common vetch.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-021-06875-z.

Genetic, anatomical, and environmental patterns related to pod shattering resistance in domesticated cowpea [Vigna unguiculata (L.) Walp]

Pod shattering, which causes the explosive release of seeds from the pod, is one of the main sources of yield losses in cowpea in arid and semi-arid areas. Reduction of shattering has therefore been a primary target for selection during domestication and improvement of cowpea, among other species. Using a mini-core diversity panel of 368 cowpea accessions, four regions with a statistically significant association with pod shattering were identified. Two genes (Vigun03g321100 and Vigun11g100600), involved in cell wall biosynthesis, were identified as strong candidates for pod shattering. Microscopic analysis was conducted on a subset of accessions representing the full spectrum of shattering phenotypes. This analysis indicated that the extent of wall fiber deposition was highly correlated with shattering. The results from this study also demonstrate that pod shattering in cowpea is exacerbated by arid environmental conditions. Finally, using a subset of West African landraces, patterns of historical selection for shattering resistance related to precipitation in the environment of origin were identified. Together, these results shed light on sources of resistance to pod shattering, which will, in turn, improve climate resilience of a major global nutritional staple.

Pod shattering in grain legumes: emerging genetic and environment-related patterns

A reduction in pod shattering is one of the main components of grain legume domestication. Despite this, many domesticated legumes suffer serious yield losses due to shattering, particularly under arid conditions. Mutations related to pod shattering modify the twisting force of pod walls or the structural strength of the dehiscence zone in pod sutures. At a molecular level, a growing body of evidence indicates that these changes are controlled by a relatively small number of key genes that have been selected in parallel across grain legume species, supporting partial molecular convergence. Legume homologs of Arabidopsis thaliana silique shattering genes play only minor roles in legume pod shattering. Most domesticated grain legume species contain multiple shattering-resistance genes, with mutants of each gene typically showing only partial shattering resistance. Hence, crosses between varieties with different genes lead to transgressive segregation of shattering alleles, producing plants with either enhanced shattering resistance or atavistic susceptibility to the trait. The frequency of these resistance pod-shattering alleles is often positively correlated with environmental aridity. The continued development of pod-shattering-related functional information will be vital for breeding crops that are suited to the increasingly arid conditions expected in the coming decades.

Search for Resistant Genotypes to Cuscuta campestris Infection in Two Legume Species, Vicia sativa and Vicia ervilia

The dodders (Cuscuta spp.) are parasitic plants that feed on the stems of their host plants. Cuscuta campestris is one of the most damaging parasitic plants for the worldwide agricultural production of broad-leaved crops. Its control is limited or non-existent, therefore resistance breeding is the best alternative both economically and environmentally. Common vetch (Vicia sativa) and bitter vetch (Vicia ervilia) are highly susceptible to C. campestris, but no resistant genotypes have been identified. Thus, the aim of this study was to identify in V. sativa and V.ervilia germplasm collections genotypes resistant to C. campestris infection for use in combating this parasitic plant. Three greenhouse screening were conducted to: (1) identify resistant responses in a collection of 154 accessions of bitter vetch and a collection of 135 accessions of common vetch genotypes against infection of C. campestris; (2) confirm the resistant response identified in common vetch accessions; and (3) characterize the effect of C. campestris infection on biomass of V. sativa resistant and susceptible accessions. Most common vetch and bitter vetch genotypes tested were susceptible to C. campestris. However, the V. sativa genotype Vs.1 exhibited high resistance. The resistant phenotype was characterized by a delay in the development of C. campestris posthaustorial growth and a darkening resembling a hypersensitive-like response at the penetration site. The resistant mechanism was effective in limiting the growth of C. campestris as the ratio of parasite/host shoot dry biomass was more significantly reduced than the rest of the accessions. To the best or our knowledge, this is the first identification of Cuscuta resistance in V. sativa genotypes.

First approach to pod dehiscence in faba bean: genetic and histological analyses

Pod dehiscence causes important yield losses in cultivated crops and therefore has been a key trait strongly selected against in crop domestication. In spite of the growing knowledge on the genetic basis of dehiscence in different crops, no information is available so far for faba bean. Here we conduct the first comprehensive study for faba bean pod dehiscence by combining, linkage mapping, comparative genomics, QTL analysis and histological examination of mature pods. Mapping of dehiscence-related genes revealed conservation of syntenic blocks among different legumes. Three QTLs were identified in faba bean chromosomes II, IV and VI, although none of them was stable across years. Histological analysis supports the convergent phenotypic evolution previously reported in cereals and related legume species but revealed a more complex pattern in faba bean. Contrary to common bean and soybean, the faba bean dehiscence zone appears to show functional equivalence to that described in crucifers. The lignified wall fiber layer, which is absent in the paucijuga primitive line Vf27, or less lignified and vacuolated in other dehiscent lines, appears to act as the major force triggering pod dehiscence in this species. While our findings, provide new insight into the mechanisms underlying faba bean dehiscence, full understanding of the molecular bases will require further studies combining precise phenotyping with genomic analysis.

Genome-Wide Identification and Expression Analysis of the Dof Gene Family in Medicago sativa L. Under Various Abiotic Stresses

The Dof transcription factor is a plant-specific transcriptional regulator that plays important roles in plant development and acts as a mediator in plant external stress responses. However, Dofs have previously been identified in several plants but not in alfalfa (Medicago sativa L.), one of the most widely cultivated forage legumes. In the present study, a total of 40 MsDof genes were identified, and the phylogenetic reconstruction, classification, conserved motifs, and expression patterns under abscisic acid (ABA), cold, heat, drought and salt stresses of these Dof genes were comprehensively analyzed. The Dof genes family in alfalfa could be classified into eight classes. Gene ontology (GO) and tissue-specific analysis indicated that most MsDof genes may be involved in biological functions during plant growth. Moreover, the expression profiles and quantitative real-time PCR analysis indicated that eight candidate abiotic tolerance genes were induced in response to four abiotic stresses. This study identified the possibility of abiotic tolerance candidate genes playing various roles in stress resistance at the whole genome level, which would provide new information on the Dof family in alfalfa.

Coordinated mechanisms of leaves and roots in response to drought stress underlying full-length transcriptome profiling in Vicia sativa L

BACKGROUND

Common vetch (Vicia sativa L.) is an important self-pollinating annual forage legume and is of interest for drought prone regions as a protein source to feed livestock and human consumption. However, the development and production of common vetch are negatively affected by drought stress. Plants have evolved common or distinct metabolic pathways between the aboveground and underground in response to drought stress. Little is known regarding the coordinated response of aboveground and underground tissues of common vetch to drought stress.

RESULTS

Our results showed that a total of 30,427 full-length transcripts were identified in 12 samples, with an average length of 2278.89 bp. Global transcriptional profiles of the above 12 samples were then analysed via Illumina-Seq. A total of 3464 and 3062 differentially expressed genes (DEGs) were identified in the leaves and roots, respectively. Gene Ontology (GO) enrichment analyses identified that the dehydrin genes and Δ¹-pyrroline-5-carboxylate synthase were induced for the biosynthesis of proline and water conservation. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis results indicated that the DEGs were significantly enriched in hormone signal transduction, starch and sucrose metabolism, and arginine and proline metabolism, and various drought response candidate genes were also identified. Abscisic acid (ABA; the AREB/ABF-SnRK2 pathway) regulates the activity of AMY3 and BAM1 to induce starch degradation in leaves and increase carbon export to roots, which may be associated with the drought stress responses in common vetch. Among the co-induced transcription factors (TFs), AREB/ABF, bHLH, MYB, WRKY, and AP2/ERF had divergent expression patterns and may be key in the crosstalk between leaves and roots during adaption to drought stress. In transgenic yeast, the overexpression of four TFs increased yeast tolerance to osmotic stresses.

CONCLUSION

The multipronged approach identified in the leaves and roots broadens our understanding of the coordinated mechanisms of drought response in common vetch, and further provides targets to improve drought resistance through genetic engineering.

Pod-shattering characteristics differences between two groups of soybeans are associated with specific changes in gene expression

Soybean is an economically important leguminous crop, and pod dehiscence of soybean could cause huge yield loss. In this study, we measured fruit-cracking forces and percentages of dehisced pods for ten soybean accessions, then separated them into two groups as shattering-sensitive (SS) and shattering-resistant (SR) soybeans. Pod transcriptomes from these two groups were analyzed, and 225 differentially expressed genes (DEGs) were identified between SS and SR soybeans. Some of these DEGs have been previously reported to be associated with pod dehiscence in soybean. The expression patterns of selected DEGs were validated by real-time quantitative reverse transcription PCR, which confirmed the expression changes found in RNA-seq analysis. We also de novo identified 246 soybean pod-long intergenic ncRNAs (lincRNAs), 401 intronic lncRNAs, and 23 antisense lncRNAs from these transcriptomes. Furthermore, genes and lincRNAs co-expression network analysis showed that there are distinct expression patterns between SS and SR soybeans in some co-expression modules. In conclusion, we systematically investigated potential genes and molecular pathways as candidates for differences in soybean pod dehiscence and will provide a useful resource for molecular breeding of soybeans.

Molecular bases for drought tolerance in common vetch: designing new molecular breeding tools

BACKGROUND

Common vetch (Vicia sativa L.) is a forage grain legume of high protein content and high nitrogen fixation, relevant in sustainable agriculture systems. Drought is the main limiting factor of this crop yield. Genetic resources collections are essential to provide genetic variability for breeding. The analysis of drought associated parameters has allowed us to identify drought tolerant and sensitive ecotypes in a vetch core collection.

RESULTS

To understand the mechanisms involved in drought response we analysed transcriptomic differences between tolerant and sensitive accessions. Polymorphic variants (SNPs and SSRs) in these differential expressed genes (DEGs) have also been analysed for the design of drought-associated markers. A total of 1332 transcripts were commonly deregulated in both genotypes under drought. To know the drought adaptive response, we also analysed DEGs between accessions. A total of 2646 transcripts are DEG between sensitive and tolerant ecotypes, in watered and drought conditions, including important genes involved in redox homeostasis, cell wall modifications and stress-response. The integration of this functional and genetic information will contribute to understand the molecular mechanisms of drought response and the adaptive mechanisms of drought tolerance in common vetch. The identification of polymorphic variants in these DEGs has also been screened for the design of drought-associated markers that could be used in future breeding program strategies.

CONCLUSIONS

Our studies shed light for the first time in common vetch about the genes and pathways associated with drought tolerance. In addition, we identify over 100 potential drought associated polymorphism, as SNPs or SSRs, which are differently present in drought and tolerant genotypes. The use of these molecular markers for trait prediction would enable the development of genomic tools for future engineering strategies by screening of germplasm crop collections for traits related with crop drought resilience, adaptability or yield in vetch.

Pod Dehiscence in Hairy Vetch (Vicia villosa Roth)

Hairy vetch, Vicia villosa (Roth), is a cover crop that does not exhibit a typical domestication syndrome. Pod dehiscence reduces seed yield and creates weed problems for subsequent crops. Breeding efforts aim to reduce pod dehiscence in hairy vetch. To characterize pod dehiscence in the species, we quantified visual dehiscence and force required to cause dehiscence among 606 genotypes grown among seven environments of the United States. To identify potential secondary selection traits, we correlated pod dehiscence with various morphological pod characteristics and field measurements. Genotypes of hairy vetch exhibited wide variation in pod dehiscence, from completely indehiscent to completely dehiscent ratings. Mean force to dehiscence also varied widely, from 0.279 to 8.97 N among genotypes. No morphological traits were consistently correlated with pod dehiscence among environments where plants were grown. Results indicated that visual ratings of dehiscence would efficiently screen against genotypes with high pod dehiscence early in the breeding process. Force to dehiscence may be necessary to identify the indehiscent genotypes during advanced stages of selection.

MYB transcription factors in alfalfa (Medicago sativa): genome-wide identification and expression analysis under abiotic stresses

Background

Alfalfa is the most widely cultivated forage legume and one of the most economically valuable crops in the world. Its survival and production are often hampered by environmental changes. However, there are few studies on stress-resistance genes in alfalfa because of its incomplete genomic information and rare expression profile data. The MYB proteins are characterized by a highly conserved DNA-binding domain, which is large, functionally diverse, and represented in all eukaryotes. The role of MYB proteins in plant development is essential; they function in diverse biological processes, including stress and defense responses, and seed and floral development. Studies on the MYB gene family have been reported in several species, but they have not been comprehensively analyzed in alfalfa.

Methods

To identify more comprehensive MYB transcription factor family genes, the sequences of 168 Arabidopsis thaliana, 430 Glycine max, 185 Medicago truncatula, and 130 Oryza sativa MYB proteins were downloaded from the Plant Transcription Factor Database. These sequences were used as queries in a BLAST search against the M. sativa proteome sequences provided by the Noble Research Institute.

Results

In the present study, a total of 265 MsMYB proteins were obtained, including 50 R1-MYB, 186 R2R3-MYB, 26 R1R2R3-MYB, and three atypical-MYB proteins. These predicted MsMYB proteins were divided into 12 subgroups by phylogenetic analysis, and gene ontology (GO) analysis indicated that most of the MsMYB genes are involved in various biological processes. The expression profiles and quantitative real-time PCR analysis indicated that some MsMYB genes might play a crucial role in the response to abiotic stresses. Additionally, a total of 170 and 914 predicted protein–protein and protein-DNA interactions were obtained, respectively. The interactions between MsMYB043 and MSAD320162, MsMYB253 and MSAD320162, and MsMYB253 and MSAD308489 were confirmed by a yeast two-hybrid system. This work provides information on the MYB family in alfalfa that was previously lacking and might promote the cultivation of stress-resistant alfalfa.

Morphological and Genetic Mechanisms Underlying Awn Development in Monocotyledonous Grasses

The identification of biological mechanisms underlying the development of complex quantitative traits, including those that contribute to plant architecture, yield and quality potential, and seed dispersal, is a major focus in the evolutionary biology and plant breeding. The awn, a bristle-like extension from the lemma in the floret, is one of the distinct morphological and physiological traits in grass species. Awns are taught as an evolutionary trait assisting seed dispersal and germination and increasing photosynthesis. Awn development seems to be complex process, involving dramatic phenotypic and molecular changes. Although recent advances investigated the underlying morphological and molecular genetic factors of awn development, there is little agreement about how these factors interact during awn formation and how this interaction affects variation of awn morphology. Consequently, the developmental sequence of the awn is not yet well understood. Here, we review awn morphological and histological features, awn development pathways, and molecular processes of awn development. We argue that morphological and molecular genetic mechanisms of awn development previously studied in major cereal crops, such as barley, wheat, and rice, offered intriguing insights helping to characterize this process in a comparative approach. Applying such an approach will aid to deeply understand factors involved in awn development in grass species.

Transcriptome-Wide Characterization and Functional Identification of the Aquaporin Gene Family During Drought Stress in Common Vetch

Aquaporins (AQPs) are transmembrane channels that are essential for the movement of water and other small molecules between biofilms in various physiological processes in plants. In this study, based on transcriptome-wide data, we identified and described a total of 21 AQP genes in common vetch (Vicia sativa subsp. sativa), which is an economically important pasture legume worldwide. Based on phylogenetic analyses, the VsAQPs were sorted into four subfamilies, including four plasma membrane intrinsic proteins (PIPs), six tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins, and four small basic intrinsic proteins. Furthermore, chemical and physical properties of these VsAQPs, including the isoelectric point and theoretical molecular weight, were analyzed. Analyses of the AQP signature sequences and key residues indicated the substrate specificity of each VsAQP. A set of VsAQPs was selected for gene expression analysis in a number of tissues and after drought stress treatments using real-time quantitative reverse transcription/polymerase chain reaction assays. Most of the PIPs and TIPs were proposed to have critical roles in regulating the flow of water during drought stress. Heterologous expression experiments in yeast indicated that VsPIP1;2 and VsPIP2;2 are key candidate genes for improving drought stress tolerance. The results reported in this study could be a crucial resource for further practical analyses and for genetic improvement of drought stress tolerance in common vetch.

Full-length transcript sequencing and comparative transcriptomic analysis to evaluate the contribution of osmotic and ionic stress components towards salinity tolerance in the roots of cultivated alfalfa (Medicago sativa L.)

BACKGROUND

Alfalfa is the most extensively cultivated forage legume. Salinity is a major environmental factor that impacts on alfalfa's productivity. However, little is known about the molecular mechanisms underlying alfalfa responses to salinity, especially the relative contribution of the two important components of osmotic and ionic stress.

RESULTS

In this study, we constructed the first full-length transcriptome database for alfalfa root tips under continuous NaCl and mannitol treatments for 1, 3, 6, 12, and 24 h (three biological replicates for each time points, including the control group) via PacBio Iso-Seq. This resulted in the identification of 52,787 full-length transcripts, with an average length of 2551 bp. Global transcriptional changes in the same 33 stressed samples were then analyzed via BGISEQ-500 RNA-Seq. Totals of 8861 NaCl-regulated and 8016 mannitol-regulated differentially expressed genes (DEGs) were identified. Metabolic analyses revealed that these DEGs overlapped or diverged in the cascades of molecular networks involved in signal perception, signal transduction, transcriptional regulation, and antioxidative defense. Notably, several well characterized signalling pathways, such as CDPK, MAPK, CIPK, and PYL-PP2C-SnRK2, were shown to be involved in osmotic stress, while the SOS core pathway was activated by ionic stress. Moreover, the physiological shifts of catalase and peroxidase activity, glutathione and proline content were in accordance with dynamic transcript profiles of the relevant genes, indicating that antioxidative defense system plays critical roles in response to salinity stress.

CONCLUSIONS

Overall, our study provides evidence that the response to salinity stress in alfalfa includes both osmotic and ionic components. The key osmotic and ionic stress-related genes are candidates for future studies as potential targets to improve resistance to salinity stress via genetic engineering.

Multiple Regulatory Networks Are Activated during Cold Stress in Medicago sativa L

Cultivated alfalfa (Medicago sativa L.) is one of the most important perennial legume forages in the world, and it has considerable potential as a valuable forage crop for livestock. However, the molecular mechanisms underlying alfalfa responses to cold stress are largely unknown. In this study, the transcriptome changes in alfalfa under cold stress at 4 °C for 2, 6, 24, and 48 h (three replicates for each time point) were analyzed using the high-throughput sequencing platform, BGISEQ-500, resulting in the identification of 50,809 annotated unigenes and 5283 differentially expressed genes (DEGs). Metabolic pathway enrichment analysis demonstrated that the DEGs were involved in carbohydrate metabolism, photosynthesis, plant hormone signal transduction, and the biosynthesis of amino acids. Moreover, the physiological changes of glutathione and proline content, catalase, and peroxidase activity were in accordance with dynamic transcript profiles of the relevant genes. Additionally, some transcription factors might play important roles in the alfalfa response to cold stress, as determined by the expression pattern of the related genes during 48 h of cold stress treatment. These findings provide valuable information for identifying and characterizing important components in the cold signaling network in alfalfa and enhancing the understanding of the molecular mechanisms underlying alfalfa responses to cold stress.

Pod Shattering: A Homologous Series of Variation Underlying Domestication and an Avenue for Crop Improvement

In wild habitats, fruit dehiscence is a critical strategy for seed dispersal; however, in cultivated crops it is one of the major sources of yield loss. Therefore, indehiscence of fruits, pods, etc., was likely to be one of the first traits strongly selected in crop domestication. Even with the historical selection against dehiscence in early domesticates, it is a trait still targeted in many breeding programs, particularly in minor or underutilized crops. Here, we review dehiscence in pulse (grain legume) crops, which are of growing importance as a source of protein in human and livestock diets, and which have received less attention than cereal crops and the model plant Arabidopsis thaliana. We specifically focus on the (i) history of indehiscence in domestication across legumes, (ii) structures and the mechanisms involved in shattering, (iii) the molecular pathways underlying this important trait, (iv) an overview of the extent of crop losses due to shattering, and the effects of environmental factors on shattering, and, (v) efforts to reduce shattering in crops. While our focus is mainly pulse crops, we also included comparisons to crucifers and cereals because there is extensive research on shattering in these taxa.

Potential DNA barcodes for Melilotus species based on five single loci and their combinations

Melilotus, an annual or biennial herb, belongs to the tribe Trifolieae (Leguminosae) and consists of 19 species. As an important green manure crop, diverse Melilotus species have different values as feed and medicine. To identify different Melilotus species, we examined the efficiency of five candidate regions as barcodes, including the internal transcribed spacer (ITS) and two chloroplast loci, rbcL and matK, and two non-coding loci, trnH-psbA and trnL-F. In total, 198 individuals from 98 accessions representing 18 Melilotus species were sequenced for these five potential barcodes. Based on inter-specific divergence, we analysed sequences and confirmed that each candidate barcode was able to identify some of the 18 species. The resolution of a single barcode and its combinations ranged from 33.33% to 88.89%. Analysis of pairwise distances showed that matK+rbcL+trnL-F+trnH-psbA+ITS (MRTPI) had the greatest value and rbcL the least. Barcode gap values and similarity value analyses confirmed these trends. The results indicated that an ITS region, successfully identifying 13 of 18 species, was the most appropriate single barcode and that the combination of all five potential barcodes identified 16 of the 18 species. We conclude that MRTPI is the most effective tool for Melilotus species identification. Taking full advantage of the barcode system, a clear taxonomic relationship can be applied to identify Melilotus species and enhance their practical production.