Bridging model and crop legumes through comparative genomics

Wild Lathyrus-A Treasure of Novel Diversity

Grasspea (Lathyrus sativus L.) is a climate-smart legume crop with adaptation to fragile agroecosystems. The genus Lathyrus is recognized for its vast genetic diversity, encompassing over 160 species, many of which are cultivated for various purposes across different regions of the world. Among these, Lathyrus sativus is widely cultivated as food, feed, and fodder in South Asia, Sub-Saharan Africa, and the Central and West Asia and North Africa (CWANA) regions. Its global cultivation has declined substantially due to the stigma posed by the presence of neurotoxin β-N-oxalyl-L-α, β-diaminopropionic acid (β-ODAP) in its seeds and green foliage. Overconsumption for a longer period of grasspea seeds harvested from landraces may lead to a neurological disorder called neurolathyrism in humans. ODAP is an obstacle for grasspea expansion, but crop wild relatives (CWRs) have been found to offer a solution. The incorporation of CWRs, particularly Lathyrus cicera, and landraces into breeding programs may reduce the ODAP content in grasspea varieties to a safer level. Recent advances in genomics-assisted breeding have expanded the potential for utilizing challenging CWRs to develop grasspea varieties that combine ultra-low ODAP levels with improved yield, stability, and adaptability. Further progress in omics technologies-such as transcriptomics, proteomics, and metabolomics-along with genome sequencing and editing, has greatly accelerated the development of grasspea varieties with reduced or zero ODAP content, while also enhancing the plant's agronomic value. This review highlights the significance of utilizing CWRs in pre-breeding programs, and harnessing advanced tools and technologies to enhance the performance, adaptability, and resilience of grasspea in response to changing environmental conditions.

A comparative plastome approach enhances the assessment of genetic variation in the Melilotus genus

BACKGROUND

Melilotus, a member of the Fabaceae family, is a pivotal forage crop that is extensively cultivated in livestock regions globally due to its notable productivity and ability to withstand abiotic stress. However, the genetic attributes of the chloroplast genome and the evolutionary connections among different Melilotus species remain unresolved.

RESULTS

In this study, we compiled the chloroplast genomes of 18 Melilotus species and performed a comprehensive comparative analysis. Through the examination of protein-coding genes, we successfully established a robust phylogenetic tree for these species. This conclusion is further supported by the phylogeny derived from single-nucleotide polymorphisms (SNPs) across the entire chloroplast genome. Notably, our findings revealed that M. infestus, M. siculus, M. sulcatus, and M. speciosus formed a distinct subgroup within the phylogenetic tree. Additionally, the chloroplast genomes of these four species exhibit two shared inversions. Moreover, inverted repeats were observed to have reemerged in six species within the IRLC. The distribution patterns of single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) within protein-coding genes indicated that ycf1 and ycf2 accumulated nonconservative alterations during evolutionary development. Furthermore, an examination of the evolutionary rate of protein-coding genes revealed that rps18, rps7, and rpl16 underwent positive selection specifically in Melilotus.

CONCLUSIONS

We present a comparative analysis of the complete chloroplast genomes of Melilotus species. This study represents the most thorough and detailed exploration of the evolution and variability within the genus Melilotus to date. Our study provides valuable chloroplast genomic information for improving phylogenetic reconstructions and making biogeographic inferences about Melilotus and other Papilionoideae species.

Medicago truncatula quantitative resistance to a new strain of Verticillium alfalfae from Iran revealed by a genome-wide association study

Verticillium wilt is a major threat to many crops, among them alfalfa (Medicago sativa). The model plant Medicago truncatula, a close relative of alfalfa was used to study the genetic control of resistance towards a new Verticillium alfalfae isolate. The accidental introduction of pathogen strains through global trade is a threat to crop production and such new strains might also be better adapted to global warming. Isolates of V. alfalfae were obtained from alfalfa fields in Iran and characterized. The Iranian isolate AF1 was used in a genome-wide association study (GWAS) involving 242 accessions from the Mediterranean region. Root inoculations were performed with conidia at 25°C and symptoms were scored regularly. Maximum Symptom Score and Area under Disease Progess Curve were computed as phenotypic traits to be used in GWAS and for comparison to a previous study with French isolate V31.2 at 20°C. This comparison showed high correlation with a shift to higher susceptibility, and similar geographical distribution of resistant and susceptible accessions to AF1 at 25°C, with resistant accessions mainly in the western part. GWAS revealed 30 significant SNPs linked to resistance towards isolate AF1. None of them were common to the previous study with isolate V31.2 at 20°C. To confirm these loci, the expression of nine underlying genes was studied. All genes were induced in roots following inoculation, in susceptible and resistant plants. However, in resistant plants induction was higher and lasted longer. Taken together, the use of a new pathogen strain and a shift in temperature revealed a completely different genetic control compared to a previous study that demonstrated the existence of two major QTLs. These results can be useful for Medicago breeding programs to obtain varieties better adapted to future conditions.

Iron accumulation and partitioning in hydroponically grown wild and cultivated chickpea (Cicer arietinum L)

Chickpea (Cicer arietinum L.) is a staple food in many developing countries where iron (Fe) deficiency often occurs in their population. The crop is a good source of protein, vitamins, and micronutrients. Fe biofortification in chickpea can be part of long-term strategy to enhance Fe intake in human diet to help to alleviate Fe deficiency. To develop cultivars with high Fe concentration in seeds, understanding the mechanisms of absorption and translocation of Fe into the seeds is critical. An experiment was conducted using a hydroponic system to evaluate Fe accumulation in seeds and other organs at different growth stages of selected genotypes of cultivated and wild relatives of chickpea. Plants were grown in media with Fe zero and Fe added conditions. Six chickpea genotypes were grown and harvested at six different growth stages: V3, V10, R2, R5, R6, and RH for analysis of Fe concentration in roots, stems, leaves, and seeds. The relative expression of genes related to Fe-metabolism including FRO2, IRT1, NRAMP3, V1T1, YSL1, FER3, GCN2, and WEE1 was analyzed. The results showed that the highest and lowest accumulation of Fe throughout the plant growth stages were found in the roots and stems, respectively. Results of gene expression analysis confirmed that the FRO2 and IRT1 were involved in Fe uptake in chickpeas and expressed more in roots under Fe added condition. All transporter genes: NRAMP3, V1T1, YSL1 along with storage gene FER3 showed higher expression in leaves. In contrast, candidate gene WEE1 for Fe metabolism expressed more in roots under Fe affluent condition; however, GCN2 showed over-expression in roots under Fe zero condition. Current finding will contribute to better understanding of Fe translocation and metabolism in chickpea. This knowledge can further be used to develop chickpea varieties with high Fe in seeds.

Genome-wide characterization of AINTEGUMENTA-LIKE family in Medicago truncatula reveals the significant roles of AINTEGUMENTAs in leaf growth

AINTEGUMENTA-LIKE (AIL) transcription factors are widely studied and play crucial roles in plant growth and development. However, the functions of the AIL family in legume species are largely unknown. In this study, 11 MtAIL genes were identified in the model legume Medicago truncatula, of which four of them are MtANTs. In situ analysis showed that MtANT1 was highly expressed in the shoot apical meristem (SAM) and leaf primordium. Characterization of mtant1 mtant2 mtant3 mtant4 quadruple mutants and MtANT1-overexpressing plants revealed that MtANTs were not only necessary but also sufficient for the regulation of leaf size, and indicated that they mainly function in the regulation of cell proliferation during secondary morphogenesis of leaves in M. truncatula. This study systematically analyzed the MtAIL family at the genome-wide level and revealed the functions of MtANTs in leaf growth. Thus, these genes may provide a potential application for promoting the biomass of legume forages.

Molecular Characterization of the miR156/MsSPL Model in Regulating the Compound Leaf Development and Abiotic Stress Response in Alfalfa

Plant leaf patterns and shapes are spectacularly diverse. Changing the complexity of leaflet numbers is a valuable approach to increase its nutrition and photosynthesis. Alfalfa (Medicago sativa) is the most important forage legume species and has diversified compound leaf patterns, which makes it a model species for studying compound leaf development. However, transcriptomic information from alfalfa remains limited. In this study, RNA-Seq technology was used to identify 3746 differentially expressed genes (DEGs) between multifoliate and trifoliate alfalfa. Through an analysis of annotation information and expression data, SPL, one of the key regulators in modifiable plant development and abiotic stress response, was further analyzed. Here, thirty MsSPL genes were obtained from the alfalfa genome, of which 16 had the putative miR156 binding site. A tissue expression pattern analysis showed that the miR156-targeted MsSPLs were divided into two classes, namely, either tissue-specific or widely expressed in all tissues. All miR156-targeted SPLs strongly showed diversification and positive roles under drought and salt conditions. Importantly, miR156/MsSPL08 was significantly suppressed in multifoliate alfalfa. Furthermore, in the paralogous mutant of MsSPL08 isolated from Medicago truncatula, the phenotypes of mutant plants reveal that miR156/MsSPL08 is involved not only involved the branches but also especially regulates the number of leaflets. The legume is a typical compound leaf plant; the ratio of the leaflet often affects the quality of the forage. This study sheds light on new functions of SPL genes that regulate leaflet number development.

Comprehensive RNAseq analysis for identification of genes expressed under heat stress in lentil

Lentils are highly sensitive to abrupt increases in temperature during the mid to late reproductive stages, leading to severe biomass and seed yield reduction. Therefore, we carried out an RNAseq analysis between IG4258 (heat tolerant) and IG3973 (heat sensitive) lentil genotypes at the reproductive stage under both normal and heat stress conditions in the field. It resulted in 209,549 assembled transcripts and among these 161,809 transcripts had coding regions, of which 94,437 transcripts were annotated. The differential gene expression analysis showed upregulation of 678 transcripts and downregulation of 680 transcripts between the tolerant and sensitive genotypes at the early reproductive stage. While 76 transcripts were upregulated and 47 transcripts were downregulated at the late reproductive stage under heat stress conditions. The validation of 12 up-or downregulated transcripts through RT-PCR corresponded well with the expression analysis data of RNAseq, with a correlation of R²  = 0.89. Among these transcripts, the DN364_c1_g1_i9 and DN2218_c0_g1_i5 transcripts encoded enzymes involved in the tryptophan pathway, indicating that tryptophan biosynthesis plays a role under heat stress in lentil. Moreover, KEGG pathways enrichment analysis identified transcripts associated with genes encoding proteins/regulating factors related to different metabolic pathways including signal transduction, fatty acid biosynthesis, rRNA processing, ribosome biogenesis, gibberellin (GA) biosynthesis, and riboflavin biosynthesis. This analysis also identified 6852 genic-SSRs leading to the development of 4968 SSR primers that are potential genomic resources for molecular mapping of heat-tolerant genes in lentil.

Integrated Transcriptomic and Bioinformatics Analyses Reveal the Molecular Mechanisms for the Differences in Seed Oil and Starch Content Between Glycine max and Cicer arietinum

The seed oil and starch content of soybean are significantly different from that of chickpea. However, there are limited studies on its molecular mechanisms. To address this issue, we conducted integrated transcriptomic and bioinformatics analyses for species-specific genes and acyl-lipid-, starch-, and carbon metabolism-related genes. Among seven expressional patterns of soybean-specific genes, four were highly expressed at the middle- and late oil accumulation stages; these genes significantly enriched fatty acid synthesis and carbon metabolism, and along with common acetyl CoA carboxylase (ACCase) highly expressed at soybean middle seed development stage, common starch-degrading enzyme beta-amylase-5 (BAM5) was highly expressed at soybean early seed development stage and oil synthesis-related genes ACCase, KAS, KAR, ACP, and long-chain acyl-CoA synthetase (LACS) were co-expressed with WRI1, which may result in high seed oil content and low seed starch content in soybean. The common ADP-glucose pyrophosphorylase (AGPase) was highly expressed at chickpea middle seed development stage, along with more starch biosynthesis genes co-expressed with four-transcription-factor homologous genes in chickpea than in soybean, and the common WRI1 was not co-expressed with oil synthesis genes in chickpea, which may result in high seed starch content and low seed oil content in chickpea. The above results may be used to improve chickpea seed oil content in two ways. One is to edit CaWRI1 to co-express with oil synthesis-related genes, which may increase carbon metabolites flowing to oil synthesis, and another is to increase the expression levels of miRNA159 and miRNA319 to inhibit the expression of MYB33, which may downregulate starch synthesis-related genes, making more carbon metabolites flow into oil synthesis. Our study will provide a basis for future breeding efforts to increase the oil content of chickpea seeds.

Comparative analysis of alfalfa (Medicago sativa L.) seedling transcriptomes reveals genotype-specific drought tolerance mechanisms

Drought is one of the main abiotic factors that affect alfalfa yield. The identification of genes that control this complex trait can provide important insights for alfalfa breeding. However, little is known about how alfalfa responds and adapts to drought stress, particularly in cultivars of differing drought tolerance. In this study, the drought-tolerant cultivar Dryland 'DT' and the drought-sensitive cultivar WL343HQ 'DS' were used to characterize leaf and root physiological responses and transcriptional changes in response to water deficit. Under drought stress, Dryland roots (DTR) showed more differentially expressed genes than WL343HQ roots (DSR), whereas WL343HQ leaves (DSL) showed more differentially expressed genes than Dryland leaves (DTL). Many of these genes were involved in stress-related pathways, carbohydrate metabolism, and lignin and wax biosynthesis, which may have improved the drought tolerance of alfalfa. We also observed that several genes related to ABA metabolism, root elongation, peroxidase activity, cell membrane stability, ubiquitination, and genetic processing responded to drought stress in alfalfa. We highlighted several candidate genes, including sucrose synthase, xylan 1,4-beta-xylosidase, primary-amine oxidase, and alcohol-forming fatty acyl-CoA reductase, for future studies on drought stress resistance in alfalfa and other plant species. In summary, our results reveal the unique drought adaptation and resistance characteristics of two alfalfa genotypes. These findings, which may be valuable for drought resistance breeding, warrant further gene functional analysis to augment currently available information and to clarify the drought stress regulatory mechanisms of alfalfa and other plants.

Comprehensive analysis and identification of drought-responsive candidate NAC genes in three semi-arid tropics (SAT) legume crops

Background

Chickpea, pigeonpea, and groundnut are the primary legume crops of semi-arid tropics (SAT) and their global productivity is severely affected by drought stress. The plant-specific NAC (NAM - no apical meristem, ATAF - Arabidopsis transcription activation factor, and CUC - cup-shaped cotyledon) transcription factor family is known to be involved in majority of abiotic stresses, especially in the drought stress tolerance mechanism. Despite the knowledge available regarding NAC function, not much information is available on NAC genes in SAT legume crops.

Results

In this study, genome-wide NAC proteins – 72, 96, and 166 have been identified from the genomes of chickpea, pigeonpea, and groundnut, respectively, and later grouped into 10 clusters in chickpea and pigeonpea, while 12 clusters in groundnut. Phylogeny with well-known stress-responsive NACs in Arabidopsis thaliana, Oryza sativa (rice), Medicago truncatula, and Glycine max (soybean) enabled prediction of putative stress-responsive NACs in chickpea (22), pigeonpea (31), and groundnut (33). Transcriptome data revealed putative stress-responsive NACs at various developmental stages that showed differential expression patterns in the different tissues studied. Quantitative real-time PCR (qRT-PCR) was performed to validate the expression patterns of selected stress-responsive, Ca_NAC (Cicer arietinum - 14), Cc_NAC (Cajanus cajan - 15), and Ah_NAC (Arachis hypogaea - 14) genes using drought-stressed and well-watered root tissues from two contrasting drought-responsive genotypes of each of the three legumes. Based on expression analysis, Ca_06899, Ca_18090, Ca_22941, Ca_04337, Ca_04069, Ca_04233, Ca_12660, Ca_16379, Ca_16946, and Ca_21186; Cc_26125, Cc_43030, Cc_43785, Cc_43786, Cc_22429, and Cc_22430; Ah_ann1.G1V3KR.2, Ah_ann1.MI72XM.2, Ah_ann1.V0X4SV.1, Ah_ann1.FU1JML.2, and Ah_ann1.8AKD3R.1 were identified as potential drought stress-responsive candidate genes.

Conclusion

As NAC genes are known to play role in several physiological and biological activities, a more comprehensive study on genome-wide identification and expression analyses of the NAC proteins have been carried out in chickpea, pigeonpea and groundnut. We have identified a total of 21 potential drought-responsive NAC genes in these legumes. These genes displayed correlation between gene expression, transcriptional regulation, and better tolerance against drought. The identified candidate genes, after validation, may serve as a useful resource for molecular breeding for drought tolerance in the SAT legume crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07602-5.

Agronomical parameters of host and non-host legumes inoculated with Melilotus indicus-isolated rhizobial strains in desert unreclaimed soil

In a search for identification of rhizobial strains with superior N₂-fixation efficiency and improved plant agronomic characteristics upon inoculation, four strains, 4.21, 9.17, 11.2 and 14.1, isolated from root nodules of wild-grown Melilotus indicus have been used to inoculate field-grown common bean, pea, cowpea and fenugreek plants. Uninoculated plants and those inoculated with host-specific commercial inoculants were used as a control. The root length, shoot height, shoot dry weight and root dry weight and the grain yield of the plants were determined after harvest. The content of N, organic C and carbohydrates content of the grain were also recorded. The inoculation with the strains 4.21 and 14.1 increased the grain yield of the fenugreek compared both with the uninoculated plants and those inoculated with the commercial strain ARC-1. The grain yield of the common bean treated with the strains 9.17 and 14.1 was also higher than that of the uninoculated and the commercial strains ARC-301. In contrast, none of the strains increased the grain yield of the pea and cowpea plants compared to the commercial strains ARC-201 and ARC-169, respectively. Significant increases of some agronomical parameters were observed in some plant-bacterium couples, albeit nodulation was not observed. It is possible that the positive effects of rhizobial inoculation on the agronomical parameters of the non-nodule forming legumes could be due to plant growth promotion characteristic of the strains used for inoculation. Analysis of the phylogeny of the almost complete 16S rRNA sequence of the rhizobial inoculants revealed that the strains 4.21 and 9.17 clustered together with R. skierniewicense and R. rosettiformans, respectively, and that the strains 11.2 and 14.1 grouped with E. meliloti. All the four strains produced IAA, and showed biocontrol activity against Rhizotocnia solani, Fusarium oxysporum, Pythium ultimum, Alternaria alternata and Sclerotonia rolsfi, albeit to a different extent.

Genetically diverse lentil- and faba bean-nodulating rhizobia are present in soils across Central and Southern Ethiopia

In total 196 bacterial isolates were obtained from root nodules of lentil (Lens culinaris) and faba bean (Vicia faba) grown on soil samples collected from 10 different sites in central and southern parts of Ethiopia. All isolates were identified as members of the genus Rhizobium by using recA gene sequence analysis. In the recA phylogenetic tree 195 rhizobial strains were classified into nine genospecies. The phylogeny of symbiotic genes nodC and nifH revealed five and six distinct groups respectively, largely dominated by symbiovar viciae. A multivariate analysis showed that environmental variables of the sampling sites considered in this study had more effect on the distribution and composition of the genospecies than the host legumes of the strains. Twenty representative strains, selected based on their isolation site, host plant and nodC group, were able to nodulate all lentil, faba bean, field pea (Pisum abyssinicum) and grass pea (Lathyrus sativus) plants in a greenhouse test in axenic conditions. The majority of the rhizobial strains were effective nitrogen-fixing symbionts for all tested legumes, indicating their potential to serve as broad host-range inoculants in agriculture. The present work suggests the presence of taxonomically and symbiotically diverse rhizobial species for legumes in the Viciae tribe in Ethiopia.

Genomics of Plant Disease Resistance in Legumes

The constant interactions between plants and pathogens in the environment and the resulting outcomes are of significant importance for agriculture and agricultural scientists. Disease resistance genes in plant cultivars can break down in the field due to the evolution of pathogens under high selection pressure. Thus, the protection of crop plants against pathogens is a continuous arms race. Like any other type of crop plant, legumes are susceptible to many pathogens. The dawn of the genomic era, in which high-throughput and cost-effective genomic tools have become available, has revolutionized our understanding of the complex interactions between legumes and pathogens. Genomic tools have enabled a global view of transcriptome changes during these interactions, from which several key players in both the resistant and susceptible interactions have been identified. This review summarizes some of the large-scale genomic studies that have clarified the host transcriptional changes during interactions between legumes and their plant pathogens while highlighting some of the molecular breeding tools that are available to introgress the traits into breeding programs. These studies provide valuable insights into the molecular basis of different levels of host defenses in resistant and susceptible interactions.

Salinity stress response and 'omics' approaches for improving salinity stress tolerance in major grain legumes

Sustaining yield gains of grain legume crops under growing salt-stressed conditions demands a thorough understanding of plant salinity response and more efficient breeding techniques that effectively integrate modern omics knowledge. Grain legume crops are important to global food security being an affordable source of dietary protein and essential mineral nutrients to human population, especially in the developing countries. The global productivity of grain legume crops is severely challenged by the salinity stress particularly in the face of changing climates coupled with injudicious use of irrigation water and improper agricultural land management. Plants adapt to sustain under salinity-challenged conditions through evoking complex molecular mechanisms. Elucidating the underlying complex mechanisms remains pivotal to our knowledge about plant salinity response. Improving salinity tolerance of plants demand enriching cultivated gene pool of grain legume crops through capitalizing on 'adaptive traits' that contribute to salinity stress tolerance. Here, we review the current progress in understanding the genetic makeup of salinity tolerance and highlight the role of germplasm resources and omics advances in improving salt tolerance of grain legumes. In parallel, scope of next generation phenotyping platforms that efficiently bridge the phenotyping-genotyping gap and latest research advances including epigenetics is also discussed in context to salt stress tolerance. Breeding salt-tolerant cultivars of grain legumes will require an integrated "omics-assisted" approach enabling accelerated improvement of salt-tolerance traits in crop breeding programs.

Adapting legume crops to climate change using genomic approaches

Our agricultural system and hence food security is threatened by combination of events, such as increasing population, the impacts of climate change, and the need to a more sustainable development. Evolutionary adaptation may help some species to overcome environmental changes through new selection pressures driven by climate change. However, success of evolutionary adaptation is dependent on various factors, one of which is the extent of genetic variation available within species. Genomic approaches provide an exceptional opportunity to identify genetic variation that can be employed in crop improvement programs. In this review, we illustrate some of the routinely used genomics-based methods as well as recent breakthroughs, which facilitate assessment of genetic variation and discovery of adaptive genes in legumes. Although additional information is needed, the current utility of selection tools indicate a robust ability to utilize existing variation among legumes to address the challenges of climate uncertainty.

QTL Mapping Combined With Bulked Segregant Analysis Identify SNP Markers Linked to Leaf Shape Traits in Pisum sativum Using SLAF Sequencing

Leaf shape is an important trait that influences the utilization rate of light, and affects quality and yield of pea (Pisum sativum). In the present study, a joint method of high-density genetic mapping using specific locus amplified fragment sequencing (SLAF-seq) and bulked segregant analysis (BSA) was applied to rapidly detect loci with leaf shape traits. A total of 7,146 polymorphic SLAFs containing 12,213 SNP markers were employed to construct a high-density genetic map for pea. We conducted quantitative trait locus (QTL) mapping on an F2 population to identify QTLs associated with leaf shape traits. Moreover, SLAF-BSA was conducted on the same F2 population to identify the single nucleotide polymorphism (SNP) markers linked to leaf shape in pea. Two QTLs (qLeaf_or-1, qLeaf_or-2) were mapped on linkage group 7 (LG7) for pea leaf shape. Through alignment of SLAF markers with Cicer arietinum, Medicago truncatula, and Glycine max, the pea LGs were assigned to their corresponding homologous chromosomal groups. The comparative genetic analysis showed that pea is more closely related to M. truncatula. Based on the sequencing results of two pools with different leaf shape, 179 associated markers were obtained after association analysis. The joint analysis of SLAF-seq and BSA showed that the QTLs obtained from mapping on a high-density genetic map are convincing due to the closely associated map region with the BSA results, which provided more potential markers related to leaf shape. Thus, the identified QTLs could be used in marker-assisted selection for pea breeding in the future. Our study revealed that joint analysis of QTL mapping on a high-density genetic map and BSA-seq is a cost-effective and accurate method to reveal genetic architecture of target traits in plant species without a reference genome.

Mining alfalfa (Medicago sativa L.) nodules for salinity tolerant non-rhizobial bacteria to improve growth of alfalfa under salinity stress

There are fewer reports on plant growth promoting (PGP) bacteria living in nodules as helper to tolerance to abiotic stress such as salinity and drought. The study was conducted to isolate rhizobial and non-rhizobial drought and salinity tolerant bacteria from the surface sterilized root nodules of alfalfa, grown in saline soils, and evaluate the effects of effective isolates on plant growth under salt stress. Based on drought and salinity tolerance of bacterial isolates and having multiple PGP traits, two non-rhizobial endophytic isolates and one rhizobial endophytic isolate were selected for further identification and characterization. Based on partial sequences of 16 S rRNA genes, non-rhizobial isolates and rhizobial isolate were closely related to Klebsiella sp., Kosakonia cowanii, and Sinorhizobium meliloti, respectively. None of the two non-rhizobial strains were able to form nodules on alfalfa roots under greenhouse and in vitro conditions. Co-inoculation of alfalfa plant with Klebsiella sp. A36, K. cowanii A37, and rhizobial strain S. meliloti ARh29 had a positive effect on plant growth indices under salinity stress. In addition, the single inoculation of non-rhizobial strains without rhizobial strain resulted in an increase in alfalfa growth indices compared to the plants non-inoculated and the ones inoculated with S. meliloti ARh29 alone under salinity stress, indicating that nodule non-rhizobial strains have PGP potentials and may be a promising way for improving effectiveness of Rhizobium bio-fertilizers in salt-affected soils.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

A transcriptome-SNP-derived linkage map of Apios americana (potato bean) provides insights about genome re-organization and synteny conservation in the phaseoloid legumes

KEY MESSAGE

We report a linkage map for Apios americana and describe synteny with selected warm-season legumes. A translocation event in common bean and soybean is confirmed against Apios and Vigna species. Apios (Apios americana; "apios"), a tuberous perennial legume in the Phaseoleae tribe, was widely used as a food by Native Americans. Work in the last 40 years has led to several improved breeding lines. Aspects of the pollination biology (complex floral structure and tripping mechanism) have made controlled crosses difficult, and the previous reports indicated that the plant is likely primarily an outcrosser. We used a pseudo-testcross strategy to construct a genetic map specific to the maternal parent. The map was built using single-nucleotide polymorphism markers identified by comparing the expressed sequences of individuals in the mapping population against a de novo maternal reference transcriptome assembly. The apios map consists of 11 linkage groups and 1121 recombinationally distinct loci, covering ~ 938.6 cM. By sequencing the transcriptomes of all potential pollen parents, we were able to identify the probable pollen donors and to discover new aspects of the pollination biology in apios. No selfing was observed, but multiple pollen parents were seen within individual pods. Comparisons with genome sequences in other species in the Phaseoleae showed extended synteny for most apios linkage groups. This synteny supports the robustness of the map, and also sheds light on the history of the Phaseoleae, as apios is relatively early diverging in this tribe. We detected a translocation event that separates apios and two Vigna species from Phaseolus vulgaris and Glycine max. This apios mapping work provides a general protocol for sequencing-based construction of high-density linkage maps in outcrossing species with heterogeneous pollen parents.

Genome-Wide Development of MicroRNA-Based SSR Markers in Medicago truncatula with Their Transferability Analysis and Utilization in Related Legume Species

Microsatellite (simple sequence repeats, SSRs) marker is one of the most widely used markers in marker-assisted breeding. As one type of functional markers, MicroRNA-based SSR (miRNA-SSR) markers have been exploited mainly in animals, but the development and characterization of miRNA-SSR markers in plants are still limited. In the present study, miRNA-SSR markers for Medicago truncatula (M. truncatula) were developed and their cross-species transferability in six leguminous species was evaluated. A total of 169 primer pairs were successfully designed from 130 M. truncatula miRNA genes, the majority of which were mononucleotide repeats (70.41%), followed by dinucleotide repeats (14.20%), compound repeats (11.24%) and trinucleotide repeats (4.14%). Functional classification of SSR-containing miRNA genes showed that all targets could be grouped into three Gene Ontology (GO) categories: 17 in biological process, 11 in molecular function, and 14 in cellular component. The miRNA-SSR markers showed high transferability in other six leguminous species, ranged from 74.56% to 90.53%. Furthermore, 25 Mt-miRNA-SSR markers were used to evaluate polymorphisms in 20 alfalfa accessions, and the polymorphism information content (PIC) values ranged from 0.39 to 0.89 with an average of 0.71, the allele number per marker varied from 3 to 18 with an average of 7.88, indicating a high level of informativeness. The present study is the first time developed and characterized of M. truncatula miRNA-SSRs and demonstrated their utility in transferability, these novel markers will be valuable for genetic diversity analysis, marker-assisted selection and genotyping in leguminous species.

Natural rhizobial diversity helps to reveal genes and QTLs associated with biological nitrogen fixation in common bean

Common bean is one of the most important crops for human feed, and the most important legume for direct consumption by millions of people, especially in developing countries. It is a promiscuous host legume in terms of nodulation, able to associate with a broad and diverse range of rhizobia, although the competitiveness for nodulation and the nitrogen fixation capacity of most of these strains is generally low. As a result, common bean is very inefficient for symbiotic nitrogen fixation, and nitrogen has to be supplied with chemical fertilizers. In the last years, symbiotic nitrogen fixation has received increasing attention as a sustainable alternative to nitrogen fertilizers, and also as a more economic and available one in poor countries. Therefore, optimization of nitrogen fixation of bean-rhizobia symbioses and selection of efficient rhizobial strains should be a priority, which begins with the study of the natural diversity of the symbioses and the rhizobial populations associated. Natural rhizobia biodiversity that nodulates common bean may be a source of adaptive alleles acting through phenotypic plasticity. Crosses between accessions differing for nitrogen fixation may combine alleles that never meet in nature. Another way to discover adaptive genes is to use association genetics to identify loci that common bean plants use for enhanced biological nitrogen fixation and, in consequence, for marker assisted selection for genetic improvement of symbiotic nitrogen fixation. In this review, rhizobial biodiversity resources will be discussed, together with what is known about the loci that underlie such genetic variation, and the potential candidate genes that may influence the symbiosis' fitness benefits, thus achieving an optimal nitrogen fixation capacity in order to help reduce reliance on nitrogen fertilizers in common bean.

Hierarchically Aligning 10 Legume Genomes Establishes a Family-Level Genomics Platform

Mainly due to their economic importance, genomes of 10 legumes, including soybean (Glycine max), wild peanut (Arachis duranensis and Arachis ipaensis), and barrel medic (Medicago truncatula), have been sequenced. However, a family-level comparative genomics analysis has been unavailable. With grape (Vitis vinifera) and selected legume genomes as outgroups, we managed to perform a hierarchical and event-related alignment of these genomes and deconvoluted layers of homologous regions produced by ancestral polyploidizations or speciations. Consequently, we illustrated genomic fractionation characterized by widespread gene losses after the polyploidizations. Notably, high similarity in gene retention between recently duplicated chromosomes in soybean supported the likely autopolyploidy nature of its tetraploid ancestor. Moreover, although most gene losses were nearly random, largely but not fully described by geometric distribution, we showed that polyploidization contributed divergently to the copy number variation of important gene families. Besides, we showed significantly divergent evolutionary levels among legumes and, by performing synonymous nucleotide substitutions at synonymous sites correction, redated major evolutionary events during their expansion. This effort laid a solid foundation for further genomics exploration in the legume research community and beyond. We describe only a tiny fraction of legume comparative genomics analysis that we performed; more information was stored in the newly constructed Legume Comparative Genomics Research Platform (www.legumegrp.org).

Microsatellite markers: what they mean and why they are so useful

Microsatellites or Single Sequence Repeats (SSRs) are extensively employed in plant genetics studies, using both low and high throughput genotyping approaches. Motivated by the importance of these sequences over the last decades this review aims to address some theoretical aspects of SSRs, including definition, characterization and biological function. The methodologies for the development of SSR loci, genotyping and their applications as molecular markers are also reviewed. Finally, two data surveys are presented. The first was conducted using the main database of Web of Science, prospecting for articles published over the period from 2010 to 2015, resulting in approximately 930 records. The second survey was focused on papers that aimed at SSR marker development, published in the American Journal of Botany's Primer Notes and Protocols in Plant Sciences (over 2013 up to 2015), resulting in a total of 87 publications. This scenario confirms the current relevance of SSRs and indicates their continuous utilization in plant science.

Understanding the Impact of Drought on Foliar and Xylem Invading Bacterial Pathogen Stress in Chickpea

In field conditions, plants are concurrently exposed to multiple stresses, where one stressor impacts the plant's response to another stressor, and the resultant net effect of these stresses differs from individual stress response. The present study investigated the effect of drought stress on interaction of chickpea with Pseudomonas syringae pv. phaseolicola (Psp; foliar pathogen) and Ralstonia solanacearum (Rs; xylem inhabiting wilt causing pathogen), respectively, and the net-effect of combined stress on chlorophyll content and cell death. Two type of stress treatments were used to study the influence of each stress factor during combined stress, viz., imposition of drought stress followed by pathogen challenge (DP), and pathogen inoculated plants imposed with drought in course of pathogen infection (PD). Drought stress was imposed at different levels with pathogen inoculum to understand the influence of different stress intensities on stress interaction and their net impact. Drought stressed chickpea plants challenged with Psp infection (DPsp) showed reduced in planta bacterial number compared to Psp infection alone. Similarly, Rs infection of chickpea plants showed reduced in planta bacterial number under severe drought stress. Combined drought and Psp (DPsp) infected plants showed decreased cell death compared to plants infected only with Psp but the extent of cell death was similar to drought stressed plants. Similarly, chlorophyll content in plants under combined stress was similar to the individual drought stressed plants; however, the chlorophyll content was more compared to pathogen only infected plants. Under combined drought and Rs infection (DRs), cell death was similar to individual drought stress but significantly less compared to only Rs infected plants. Altogether, the study proposes that both stress interaction and net effect of combined stress could be majorly influenced by first occurring stress, for example, drought stress in DP treatment. In addition, our results indicate that the outcome of the two stress interaction in plant depends on timing of stress occurrence and nature of infecting pathogen.

A gene-based map of the Nod factor-independent Aeschynomene evenia sheds new light on the evolution of nodulation and legume genomes

Aeschynomene evenia has emerged as a new model legume for the deciphering of the molecular mechanisms of an alternative symbiotic process that is independent of the Nod factors. Whereas most of the research on nitrogen-fixing symbiosis, legume genetics and genomics has so far focused on Galegoid and Phaseolid legumes, A. evenia falls in the more basal and understudied Dalbergioid clade along with peanut (Arachis hypogaea). To provide insights into the symbiotic genes content and the structure of the A. evenia genome, we established a gene-based genetic map for this species. Firstly, an RNAseq analysis was performed on the two parental lines selected to generate a F₂ mapping population. The transcriptomic data were used to develop molecular markers and they allowed the identification of most symbiotic genes. The resulting map comprised 364 markers arranged in 10 linkage groups (2n = 20). A comparative analysis with the sequenced genomes of Arachis duranensis and A. ipaensis, the diploid ancestors of peanut, indicated blocks of conserved macrosynteny. Altogether, these results provided important clues regarding the evolution of symbiotic genes in a Nod factor-independent context. They provide a basis for a genome sequencing project and pave the way for forward genetic analysis of symbiosis in A. evenia.

An update on hypoallergenicity of peanut and soybean: where are we now?

Legumes are one of the major sources of proteins and positively correlate with the development of modern society. At the same time, unfortunately, they significantly contribute to the rising prevalence of food allergy.

CNMS: The preferred genic markers for comparative genomic, molecular phylogenetic, functional genetic diversity and differential gene regulatory expression analyses in chickpea

The intra/inter-genomic comparative mapping-based phylogenetic footprinting identified 5 paralogous and 656 orthologous genome-wide CNMS markers in the upstream sequences of chickpea genes. These CNMS markers revealed a high-degree of gene-based syntenic relationship between chickpea and Medicago genomes while minimum between chickpea and Vitis genomes. The time of divergence and duplication estimated using CNMS markers highlight the expected phylogenetic relationships between chickpea and six dicot (legume) species as well as occurrence of ancient genome (approximately 53 Mya) with small-scale recent segmental (approximately 10 Mya) duplication events in chickpea. A wider level of functional molecular diversity (14 to 88 percent) and admixed population genetic structure was detected among desi, kabuli and wild genotypes by genic CNMS markers at a genome-wide scale suggesting their utility in large-scale genetic analysis in chickpea. The subfunctionalization at the cis-regulatory element region and TFBS (transcription factor binding site) motif levels in the upstream sequences of CNMS marker-associated orthologous genes than the paralogues was predominant. Functional constraint might have considerable effect on these CNMScontaining regulatory elements controlling consistent orthologous gene expression in dicots. A rapid subfunctionalization based on diverge differential expression of paralogous CNMS marker-associated genes particularly those that underwent recent small-scale segmental duplication events in chickpea was apparent. The differential regulation of expression and subfunctionalization potential of Ultra CNMS marker-associated genes suggest their utility in deciphering the complex gene regulatory function as well as identification and targeted mapping of potential genes/QTLs governing vital agronomic traits in chickpea. The gene-based CNMS markers with desirable inherent genetic attributes like higher degree of comparative genome mapping, functional genetic diversity and differential gene regulatory expression potential can significantly propel the genomics-assisted chickpea crop improvement.

Molecular characterization of a mutation affecting abscisic acid biosynthesis and consequently stomatal responses to humidity in an agriculturally important species

Mutants deficient in the phytohormone abscisic acid (ABA) have been instrumental in determining not only the biosynthetic pathway for this hormone, but also its physiological role in land plants. The wilty mutant of Pisum sativum is one of the classical, well-studied ABA-deficient mutants; however, this mutant remains uncharacterized at a molecular level. Using a candidate gene approach, we show that the wilty mutation affects the xanthoxin dehydrogenase step in ABA biosynthesis. To date, this step has only been represented by mutants in the ABA2 gene of Arabidopsis thaliana. Functional ABA biosynthesis appears to be critical for normal stomatal responses to changes in humidity in angiosperms, with wilty mutant plants having no increase in foliar ABA levels in response to a doubling in vapour pressure deficit, and no closure of stomata. Phylogenetic analysis of the ABA2 gene family from diverse land plants indicates that an ABA-biosynthesis-specific short-chain dehydrogenase (ABA2) evolved in the earliest angiosperms. The relatively recent origin of specificity in this step has important implications for both the evolution of ABA biosynthesis and action in land plants.

Enhanced salt tolerance of alfalfa (Medicago sativa) by rstB gene transformation

Generating salt tolerance forage plant is essential for use of the land affected by high salinity. A salt tolerance gene rstB was used as a selectable marker gene in Agrobacterium-mediated transformation of tobacco under a selective regime of 170mM NaCl. The transgenic plants showed clear improvement in salt tolerance. To improve salt tolerance of alfalfa (Medicago sativa L.), rstB gene was introduced into alfalfa genome by Agrobacterium-mediated transformation. No abnormal phenotype was observed among the transgenic plants when compared with wild type (wt) plants. Significant enhancement of resistance to salt-shock treatment was noted on the rstB transgenic (T0) plants. Transgenic second-generation (T1) seeds showed improved germination rate and seedling growth under salt-stress condition. Hindered Na(+) accumulation, but enhanced Ca(2+) accumulation was observed on the rstB T1 plants when subjected to salt-stresses. Enhanced calcium accumulation in transgenic plants was also verified by cytohistochemical localization of calcium. Under salt-stress of 50mM NaCl, about 15% of the transgenic plants finished their life-cycle but the wt plants had no flower formation. The results demonstrated that the expression of rstB gene improved salt tolerance in transgenic alfalfa.

Naturally occurring diversity helps to reveal genes of adaptive importance in legumes

Environmental changes challenge plants and drive adaptation to new conditions, suggesting that natural biodiversity may be a source of adaptive alleles acting through phenotypic plasticity and/or micro-evolution. Crosses between accessions differing for a given trait have been the most common way to disentangle genetic and environmental components. Interestingly, such man-made crosses may combine alleles that never meet in nature. Another way to discover adaptive alleles, inspired by evolution, is to survey large ecotype collections and to use association genetics to identify loci of interest. Both of these two genetic approaches are based on the use of biodiversity and may eventually help us in identifying the genes that plants use to respond to challenges such as short-term stresses or those due to global climate change. In legumes, two wild species, Medicago truncatula and Lotus japonicus, plus the cultivated soybean (Glycine max) have been adopted as models for genomic studies. In this review, we will discuss the resources, limitations and future plans for a systematic use of biodiversity resources in model legumes to pinpoint genes of adaptive importance in legumes, and their application in breeding.

A combinatorial approach of comprehensive QTL-based comparative genome mapping and transcript profiling identified a seed weight-regulating candidate gene in chickpea

High experimental validation/genotyping success rate (94-96%) and intra-specific polymorphic potential (82-96%) of 1536 SNP and 472 SSR markers showing in silico polymorphism between desi ICC 4958 and kabuli ICC 12968 chickpea was obtained in a 190 mapping population (ICC 4958 × ICC 12968) and 92 diverse desi and kabuli genotypes. A high-density 2001 marker-based intra-specific genetic linkage map comprising of eight LGs constructed is comparatively much saturated (mean map-density: 0.94 cM) in contrast to existing intra-specific genetic maps in chickpea. Fifteen robust QTLs (PVE: 8.8-25.8% with LOD: 7.0-13.8) associated with pod and seed number/plant (PN and SN) and 100 seed weight (SW) were identified and mapped on 10 major genomic regions of eight LGs. One of 126.8 kb major genomic region harbouring a strong SW-associated robust QTL (Caq'SW1.1: 169.1-171.3 cM) has been delineated by integrating high-resolution QTL mapping with comprehensive marker-based comparative genome mapping and differential expression profiling. This identified one potential regulatory SNP (G/A) in the cis-acting element of candidate ERF (ethylene responsive factor) TF (transcription factor) gene governing seed weight in chickpea. The functionally relevant molecular tags identified have potential to be utilized for marker-assisted genetic improvement of chickpea.

Current knowledge in lentil genomics and its application for crop improvement

Most of the lentil growing countries face a certain set of abiotic and biotic stresses causing substantial reduction in crop growth, yield, and production. Until-to date, lentil breeders have used conventional plant breeding techniques of selection-recombination-selection cycle to develop improved cultivars.These techniques have been successful in mainstreaming some of the easy-to-manage monogenic traits. However, in case of complex quantitative traits, these conventional techniques are less precise. As most of the economic traits are complex, quantitative, and often influenced by environments and genotype-environment interaction, the genetic improvement of these traits becomes difficult. Genomics assisted breeding is relatively powerful and fast approach to develop high yielding varieties more suitable to adverse environmental conditions. New tools such as molecular markers and bioinformatics are expected to generate new knowledge and improve our understanding on the genetics of complex traits. In the past, the limited availability of genomic resources in lentil could not allow breeders to employ these tools in mainstream breeding program.The recent application of the next generation sequencing and genotyping by sequencing technologies has facilitated to speed up the lentil genome sequencing project and large discovery of genome-wide single nucleotide polymorphism (SNP) markers. Currently, several linkage maps have been developed in lentil through the use of expressed sequenced tag (EST) derived simple sequence repeat (SSR) and SNP markers.These maps have emerged as useful genomic resources to identify quantitative trait loci imparting tolerance to biotic and abiotic stresses in lentil. In this review, the current knowledge on available genomic resources and its application in lentil breeding program are discussed.

Identification of molecular markers associated with Verticillium wilt resistance in alfalfa (Medicago sativa L.) using high-resolution melting

Verticillium wilt, caused by the soilborne fungus, Verticillium alfalfae, is one of the most serious diseases of alfalfa (Medicago sativa L.) worldwide. To identify loci associated with resistance to Verticillium wilt, a bulk segregant analysis was conducted in susceptible or resistant pools constructed from 13 synthetic alfalfa populations, followed by association mapping in two F1 populations consisted of 352 individuals. Simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were used for genotyping. Phenotyping was done by manual inoculation of the pathogen to replicated cloned plants of each individual and disease severity was scored using a standard scale. Marker-trait association was analyzed by TASSEL. Seventeen SNP markers significantly associated with Verticillium wilt resistance were identified and they were located on chromosomes 1, 2, 4, 7 and 8. SNP markers identified on chromosomes 2, 4 and 7 co-locate with regions of Verticillium wilt resistance loci reported in M. truncatula. Additional markers identified on chromosomes 1 and 8 located the regions where no Verticillium resistance locus has been reported. This study highlights the value of SNP genotyping by high resolution melting to identify the disease resistance loci in tetraploid alfalfa. With further validation, the markers identified in this study could be used for improving resistance to Verticillium wilt in alfalfa breeding programs.

Genetic mapping of legume orthologs reveals high conservation of synteny between lentil species and the sequenced genomes of Medicago and chickpea

Lentil (Lens culinaris Medik.) is a global food crop with increasing importance for food security in south Asia and other regions. Lens ervoides, a wild relative of cultivated lentil, is an important source of agronomic trait variation. Lens is a member of the galegoid clade of the Papilionoideae family, which includes other important dietary legumes such as chickpea (Cicer arietinum) and pea (Pisum sativum), and the sequenced model legume Medicago truncatula. Understanding the genetic structure of Lens spp. in relation to more fully sequenced legumes would allow leveraging of genomic resources. A set of 1107 TOG-based amplicons were identified in L. ervoides and a subset thereof used to design SNP markers for mapping. A map of L. ervoides consisting of 377 SNP markers spread across seven linkage groups was developed using a GoldenGate genotyping array and single SNP marker assays. Comparison with maps of M. truncatula and L. culinaris documented considerable shared synteny and led to the identification of a few major translocations and a major inversion that distinguish Lens from M. truncatula, as well as a translocation that distinguishes L. culinaris from L. ervoides. The identification of chromosome-level differences among Lens spp. will aid in the understanding of introgression of genes from L. ervoides into cultivated L. culinaris, furthering genetic research and breeding applications in lentil.

Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean (Vicia faba L.)

We have identified QTLs for stomatal characteristics on chromosome II of faba bean by applying SNPs derived from M. truncatula , and have identified candidate genes within these QTLs using synteny between the two species. Faba bean (Vicia faba L.) is a valuable food and feed crop worldwide, but drought often limits its production, and its genome is large and poorly mapped. No information is available on the effects of genomic regions and genes on drought adaptation characters such as stomatal characteristics in this species, but the synteny between the sequenced model legume, Medicago truncatula, and faba bean can be used to identify candidate genes. A mapping population of 211 F5 recombinant inbred lines (Mélodie/2 × ILB 938/2) were phenotyped to identify quantitative trait loci (QTL) affecting stomatal morphology and function, along with seed weight, under well-watered conditions in a climate-controlled glasshouse in 2013 and 2014. Canopy temperature (CT) was evaluated in 2013 under water-deficit (CTd). In total, 188 polymorphic single nucleotide polymorphisms (SNPs), developed from M. truncatula genome data, were assigned to nine linkage groups that covered ~928 cM of the faba bean genome with an average inter-marker distance of 5.8 cM. 15 putative QTLs were detected, of which eight (affecting stomatal density, length and conductance and CT) co-located on chromosome II, in the vicinity of a possible candidate gene-a receptor-like protein kinase found in the syntenic interval of M. truncatula chromosome IV. A ribose-phosphate pyrophosphokinase from M. truncatula chromosome V, postulated as a possible candidate gene for the QTL for CTd, was found some distance away in the same chromosome. These results demonstrate that genomic information from M. truncatula can successfully be translated to the faba bean genome.

Comparative and evolutionary analysis of major peanut allergen gene families

Peanut (Arachis hypogaea L.) causes one of the most serious food allergies. Peanut seed proteins, Arah1, Arah2, and Arah3, are considered to be among the most important peanut allergens. To gain insights into genome organization and evolution of allergen-encoding genes, approximately 617 kb from the genome of cultivated peanut and 215 kb from a wild relative were sequenced including three Arah1, one Arah2, eight Arah3, and two Arah6 gene family members. To assign polarity to differences between homoeologous regions in peanut, we used as outgroups the single orthologous regions in Medicago, Lotus, common bean, chickpea, and pigeonpea, which diverged from peanut about 50 Ma and have not undergone subsequent polyploidy. These regions were also compared with orthologs in many additional dicot plant species to help clarify the timing of evolutionary events. The lack of conservation of allergenic epitopes between species, and the fact that many different proteins can be allergenic, makes the identification of allergens across species by comparative studies difficult. The peanut allergen genes are interspersed with low-copy genes and transposable elements. Phylogenetic analyses revealed lineage-specific expansion and loss of low-copy genes between species and homoeologs. Arah1 syntenic regions are conserved in soybean, pigeonpea, tomato, grape, Lotus, and Arabidopsis, whereas Arah3 syntenic regions show genome rearrangements. We infer that tandem and segmental duplications led to the establishment of the Arah3 gene family. Our analysis indicates differences in conserved motifs in allergen proteins and in the promoter regions of the allergen-encoding genes. Phylogenetic analysis and genomic organization studies provide new insights into the evolution of the major peanut allergen-encoding genes.

Development of EST-SSR markers and construction of a linkage map in faba bean (Vicia faba)

To develop a high density linkage map in faba bean, a total of 1,363 FBES (Faba bean expressed sequence tag [EST]-derived simple sequence repeat [SSR]) markers were designed based on 5,090 non-redundant ESTs developed in this study. A total of 109 plants of a 'Nubaria 2' × 'Misr 3' F2 mapping population were used for map construction. Because the parents were not pure homozygous lines, the 109 F2 plants were divided into three subpopulations according to the original F1 plants. Linkage groups (LGs) generated in each subpopulation were integrated by commonly mapped markers. The integrated 'Nubaria 2' × 'Misr 3' map consisted of six LGs, representing a total length of 684.7 cM, with 552 loci. Of the mapped loci, 47% were generated from multi-loci diagnostic (MLD) markers. Alignment of homologous sequence pairs along each linkage group revealed obvious syntenic relationships between LGs in faba bean and the genomes of two model legumes, Lotus japonicus and Medicago truncatula. In a polymorphic analysis with ten Egyptian faba bean varieties, 78.9% (384/487) of the FBES markers showed polymorphisms. Along with the EST-SSR markers, the dense map developed in this study is expected to accelerate marker assisted breeding in faba bean.

Characterization of five subgroups of the sieve element occlusion gene family in Glycine max reveals genes encoding non-forisome P-proteins, forisomes and forisome tails

P-proteins are structural phloem proteins discussed to be involved in the rapid sealing of injured sieve elements. P-proteins are found in all dicotyledonous and some monocotyledonous plants, but additional crystalloid P-proteins, known as forisomes, have evolved solely in the Fabaceae. Both types are encoded by members of the sieve element occlusion (SEO) gene family, which comprises seven phylogenetic subgroups. The Fabaceae-specific subgroup 1 contains genes encoding forisome subunits in e.g. Medicago truncatula, Vicia faba, Dipteryx panamensis and Canavalia gladiata whereas basal subgroup 5 encodes P-proteins in Nicotiana tabacum (tobacco) and Arabidopsis thaliana. The function of remaining subgroups is still unknown. We chose Glycine max (soybean) as a model to investigate SEO proteins representing different subgroups in one species. We isolated native P-proteins to determine the SEO protein composition and analyzed the expression pattern, localization and structure of the G. max SEO proteins representing five of the subgroups. We found that subgroup 1 GmSEO genes encode forisome subunits, a member of subgroup 5 encodes a non-forisome P-protein and subgroup 2 GmSEO genes encode the components of forisome tails, which are present in a restricted selection of Fabaceaen species. We therefore present the first molecular characterization of a Fabaceae non-forisome P-protein and the first evidence that forisome tails are encoded by a phylogenetically-distinct branch of the SEO gene family.

Genome duplication and gene loss affect the evolution of heat shock transcription factor genes in legumes

Whole-genome duplication events (polyploidy events) and gene loss events have played important roles in the evolution of legumes. Here we show that the vast majority of Hsf gene duplications resulted from whole genome duplication events rather than tandem duplication, and significant differences in gene retention exist between species. By searching for intraspecies gene colinearity (microsynteny) and dating the age distributions of duplicated genes, we found that genome duplications accounted for 42 of 46 Hsf-containing segments in Glycine max, while paired segments were rarely identified in Lotus japonicas, Medicago truncatula and Cajanus cajan. However, by comparing interspecies microsynteny, we determined that the great majority of Hsf-containing segments in Lotus japonicas, Medicago truncatula and Cajanus cajan show extensive conservation with the duplicated regions of Glycine max. These segments formed 17 groups of orthologous segments. These results suggest that these regions shared ancient genome duplication with Hsf genes in Glycine max, but more than half of the copies of these genes were lost. On the other hand, the Glycine max Hsf gene family retained approximately 75% and 84% of duplicated genes produced from the ancient genome duplication and recent Glycine-specific genome duplication, respectively. Continuous purifying selection has played a key role in the maintenance of Hsf genes in Glycine max. Expression analysis of the Hsf genes in Lotus japonicus revealed their putative involvement in multiple tissue-/developmental stages and responses to various abiotic stimuli. This study traces the evolution of Hsf genes in legume species and demonstrates that the rates of gene gain and loss are far from equilibrium in different species.

Dt2 is a gain-of-function MADS-domain factor gene that specifies semideterminacy in soybean

Similar to Arabidopsis thaliana, the wild soybeans (Glycine soja) and many cultivars exhibit indeterminate stem growth specified by the shoot identity gene Dt1, the functional counterpart of Arabidopsis TERMINAL FLOWER1 (TFL1). Mutations in TFL1 and Dt1 both result in the shoot apical meristem (SAM) switching from vegetative to reproductive state to initiate terminal flowering and thus produce determinate stems. A second soybean gene (Dt2) regulating stem growth was identified, which, in the presence of Dt1, produces semideterminate plants with terminal racemes similar to those observed in determinate plants. Here, we report positional cloning and characterization of Dt2, a dominant MADS domain factor gene classified into the APETALA1/SQUAMOSA (AP1/SQUA) subfamily that includes floral meristem (FM) identity genes AP1, FUL, and CAL in Arabidopsis. Unlike AP1, whose expression is limited to FMs in which the expression of TFL1 is repressed, Dt2 appears to repress the expression of Dt1 in the SAMs to promote early conversion of the SAMs into reproductive inflorescences. Given that Dt2 is not the gene most closely related to AP1 and that semideterminacy is rarely seen in wild soybeans, Dt2 appears to be a recent gain-of-function mutation, which has modified the genetic pathways determining the stem growth habit in soybean.

Genomics-assisted breeding in four major pulse crops of developing countries: present status and prospects

KEY MESSAGE

Given recent advances in pulse molecular biology, genomics-driven breeding has emerged as a promising approach to address the issues of limited genetic gain and low productivity in various pulse crops. The global population is continuously increasing and is expected to reach nine billion by 2050. This huge population pressure will lead to severe shortage of food, natural resources and arable land. Such an alarming situation is most likely to arise in developing countries due to increase in the proportion of people suffering from protein and micronutrient malnutrition. Pulses being a primary and affordable source of proteins and minerals play a key role in alleviating the protein calorie malnutrition, micronutrient deficiencies and other undernourishment-related issues. Additionally, pulses are a vital source of livelihood generation for millions of resource-poor farmers practising agriculture in the semi-arid and sub-tropical regions. Limited success achieved through conventional breeding so far in most of the pulse crops will not be enough to feed the ever increasing population. In this context, genomics-assisted breeding (GAB) holds promise in enhancing the genetic gains. Though pulses have long been considered as orphan crops, recent advances in the area of pulse genomics are noteworthy, e.g. discovery of genome-wide genetic markers, high-throughput genotyping and sequencing platforms, high-density genetic linkage/QTL maps and, more importantly, the availability of whole-genome sequence. With genome sequence in hand, there is a great scope to apply genome-wide methods for trait mapping using association studies and to choose desirable genotypes via genomic selection. It is anticipated that GAB will speed up the progress of genetic improvement of pulses, leading to the rapid development of cultivars with higher yield, enhanced stress tolerance and wider adaptability.

Lathyrus diversity: available resources with relevance to crop improvement--L. sativus and L. cicera as case studies

BACKGROUND

The Lathyrus genus includes 160 species, some of which have economic importance as food, fodder and ornamental crops (mainly L. sativus, L. cicera and L. odoratus, respectively) and are cultivated in >1·5 Mha worldwide. However, in spite of their well-recognized robustness and potential as a source of calories and protein for populations in drought-prone and marginal areas, cultivation is in decline and there is a high risk of genetic erosion.

SCOPE

In this review, current and past taxonomic treatments of the Lathyrus genus are assessed and its current status is examined together with future prospects for germplasm conservation, characterization and utilization. A particular emphasis is placed on the importance of diversity analysis for breeding of L. sativus and L. cicera.

CONCLUSIONS

Efforts for improvement of L. sativus and L. cicera should concentrate on the development of publicly available joint core collections, and on high-resolution genotyping. This will be critical for permitting decentralized phenotyping. Such a co-ordinated international effort should result in more efficient and faster breeding approaches, which are particularly needed for these neglected, underutilized Lathyrus species.

New evidence of ancestral polyploidy in the Genistoid legume Lupinus angustifolius L. (narrow-leafed lupin)

This is the first clear evidence of duplication and/or triplication of large chromosomal regions in a genome of a Genistoid legume, the most basal clade of Papilionoid legumes. Lupinus angustifolius L. (narrow-leafed lupin) is the most widely cultivated species of Genistoid legume, grown for its high-protein grain. As a member of this most basal clade of Papilionoid legumes, L. angustifolius serves as a useful model for exploring legume genome evolution. Here, we report an improved reference genetic map of L. angustifolius comprising 1207 loci, including 299 newly developed Diversity Arrays Technology markers and 54 new gene-based PCR markers. A comparison between the L. angustifolius and Medicago truncatula genomes was performed using 394 sequence-tagged site markers acting as bridging points between the two genomes. The improved L. angustifolius genetic map, the updated M. truncatula genome assembly and the increased number of bridging points between the genomes together substantially enhanced the resolution of synteny and chromosomal colinearity between these genomes compared to previous reports. While a high degree of syntenic fragmentation was observed that was consistent with the large evolutionary distance between the L. angustifolius and M. truncatula genomes, there were striking examples of conserved colinearity of loci between these genomes. Compelling evidence was found of large-scale duplication and/or triplication in the L. angustifolius genome, consistent with one or more ancestral polyploidy events.

Genetic analysis of tolerance to the root lesion nematode Pratylenchus neglectus in the legume Medicago littoralis

BACKGROUND

The nematode Pratylenchus neglectus has a wide host range and is able to feed on the root systems of cereals, oilseeds, grain and pasture legumes. Under the Mediterranean low rainfall environments of Australia, annual Medicago pasture legumes are used in rotation with cereals to fix atmospheric nitrogen and improve soil parameters. Considerable efforts are being made in breeding programs to improve resistance and tolerance to Pratylenchus neglectus in the major crops wheat and barley, which makes it vital to develop appropriate selection tools in medics.

RESULTS

A strong source of tolerance to root damage by the root lesion nematode (RLN) Pratylenchus neglectus had previously been identified in line RH-1 (strand medic, M. littoralis). Using RH-1, we have developed a single seed descent (SSD) population of 138 lines by crossing it to the intolerant cultivar Herald. After inoculation, RLN-associated root damage clearly segregated in the population. Genetic analysis was performed by constructing a genetic map using simple sequence repeat (SSR) and gene-based SNP markers. A highly significant quantitative trait locus (QTL), QPnTolMl.1, was identified explaining 49% of the phenotypic variation in the SSD population. All SSRs and gene-based markers in the QTL region were derived from chromosome 1 of the sequenced genome of the closely related species M. truncatula. Gene-based markers were validated in advanced breeding lines derived from the RH-1 parent and also a second RLN tolerance source, RH-2 (M. truncatula ssp. tricycla). Comparative analysis to sequenced legume genomes showed that the physical QTL interval exists as a synteny block in Lotus japonicus, common bean, soybean and chickpea. Furthermore, using the sequenced genome information of M. truncatula, the QTL interval contains 55 genes out of which five are discussed as potential candidate genes responsible for the mapped tolerance.

CONCLUSION

The closely linked set of SNP-based PCR markers is directly applicable to select for two different sources of RLN tolerance in breeding programs. Moreover, genome sequence information has allowed proposing candidate genes for further functional analysis and nominates QPnTolMl.1 as a target locus for RLN tolerance in economically important grain legumes, e.g. chickpea.

Natural product biosynthesis in Medicago species

The genus Medicago, a member of the legume (Fabaceae) family, comprises 87 species of flowering plants, including the forage crop M. sativa (alfalfa) and the model legume M. truncatula (barrel medic). Medicago species synthesize a variety of bioactive natural products that are used to engage into symbiotic interactions but also serve to deter pathogens and herbivores. For humans, these bioactive natural products often possess promising pharmaceutical properties. In this review, we focus on the two most interesting and well characterized secondary metabolite classes found in Medicago species, the triterpene saponins and the flavonoids, with a detailed overview of their biosynthesis, regulation, and profiling methods. Furthermore, their biological role within the plant as well as their potential utility for human health or other applications is discussed. Finally, we give an overview of the advances made in metabolic engineering in Medicago species and how the development of novel molecular and omics toolkits can influence a better understanding of this genus in terms of specialized metabolism and chemistry. Throughout, we critically analyze the current bottlenecks and speculate on future directions and opportunities for research and exploitation of Medicago metabolism.

Development and characterization of simple sequence repeat (SSR) markers based on RNA-sequencing of Medicago sativa and in silico mapping onto the M. truncatula genome

Sufficient codominant genetic markers are needed for various genetic investigations in alfalfa since the species is an outcrossing autotetraploid. With the newly developed next generation sequencing technology, a large amount of transcribed sequences of alfalfa have been generated and are available for identifying SSR markers by data mining. A total of 54,278 alfalfa non-redundant unigenes were assembled through the Illumina HiSeqTM 2000 sequencing technology. Based on 3,903 unigene sequences, 4,493 SSRs were identified. Tri-nucleotide repeats (56.71%) were the most abundant motif class while AG/CT (21.7%), AGG/CCT (19.8%), AAC/GTT (10.3%), ATC/ATG (8.8%), and ACC/GGT (6.3%) were the subsequent top five nucleotide repeat motifs. Eight hundred and thirty- seven EST-SSR primer pairs were successfully designed. Of these, 527 (63%) primer pairs yielded clear and scored PCR products and 372 (70.6%) exhibited polymorphisms. High transferability was observed for ssp falcata at 99.2% (523) and 71.7% (378) in M. truncatula. In addition, 313 of 527 SSR marker sequences were in silico mapped onto the eight M. truncatula chromosomes. Thirty-six polymorphic SSR primer pairs were used in the genetic relatedness analysis of 30 Chinese alfalfa cultivated accessions generating a total of 199 scored alleles. The mean observed heterozygosity and polymorphic information content were 0.767 and 0.635, respectively. The codominant markers not only enriched the current resources of molecular markers in alfalfa, but also would facilitate targeted investigations in marker-trait association, QTL mapping, and genetic diversity analysis in alfalfa.