Classification and characterization of species within the genus lens using genotyping-by-sequencing (GBS)

Genetic structure and ecological niche space of lentil's closest wild relative, Lens orientalis (Boiss.) Schmalh

Crops arose from wild ancestors and to understand their domestication it is essential to compare the cultivated species with their crop wild relatives. These represent an important source of further crop improvement, in particular in relation to climate change. Although there are about 58,000 Lens accessions held in genebanks, only 1% are wild. We examined the geographic distribution and genetic diversity of the lentil's immediate progenitor L. orientalis. We used Genotyping by Sequencing (GBS) to identify and characterize differentiation among accessions held at germplasm collections. We then determined whether genetically distinct clusters of accessions had been collected from climatically distinct locations. Of the 195 genotyped accessions, 124 were genuine L. orientalis with four identified genetic groups. Although an environmental distance matrix was significantly correlated with geographic distance in a Mantel test, the four identified genetic clusters were not found to occupy significantly different environmental space. Maxent modelling gave a distinct predicted distribution pattern centred in the Fertile Crescent, with intermediate probabilities of occurrence in parts of Turkey, Greece, Cyprus, Morocco, and the south of the Iberian Peninsula with NW Africa. Future projections did not show any dramatic alterations in the distribution according to the climate change scenarios tested. We have found considerable diversity in L. orientalis, some of which track climatic variability. The results of the study showed the genetic diversity of wild lentil and indicate the importance of ongoing collections and in situ conservation for our future capacity to harness the genetic variation of the lentil progenitor.

Drought stress in Lens culinaris: effects, tolerance mechanism, and its smart reprogramming by using modern biotechnological approaches

Among legumes, lentil serves as an imperative source of dietary proteins and are considered an important pillar of global food and nutritional security. The crop is majorly cultivated in arid and semi-arid regions and exposed to different abiotic stresses. Drought stress is a polygenic stress that poses a major threat to the crop productivity of lentils. It negatively influenced the seed emergence, water relations traits, photosynthetic machinery, metabolites, seed development, quality, and yield in lentil. Plants develop several complex physiological and molecular protective mechanisms for tolerance against drought stress. These complicated networks are enabled to enhance the cellular potential to survive under extreme water-scarce conditions. As a result, proper drought stress-mitigating novel and modern approaches are required to improve lentil productivity. The currently existing biotechnological techniques such as transcriptomics, genomics, proteomics, metabolomics, CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/cas9), and detection of QTLs (quantitative trait loci), proteins, and genes responsible for drought tolerance have gained appreciation among plant breeders for developing climate-resilient lentil varieties. In this review, we critically elaborate the impact of drought on lentil, mechanisms employed by plants to tolerate drought, and the contribution of omics approaches in lentils for dealing with drought, providing deep insights to enhance lentil productivity and improve resistance against abiotic stresses. We hope this updated review will directly help the lentil breeders to develop resistance against drought stress.

Graphical Abstract

Limited-transpiration trait in response to high vapor pressure deficit from wild to cultivated species: study of the Lens genus

Lentil (Lens culinaris Medik.) is commonly grown in drought-prone areas where terminal heat and drought are frequent. The limited-transpiration (TRlim) trait under high vapor pressure deficit (VPD) could be a way to conserve water and increase yield under water deficit conditions. The TRlim trait was examined in cultivated and wild lentil species together with its evolution throughout the breeding pipeline. Sixty-one accessions representing the six wild lentil species (L. orientalis, L. tomentosus, L. odemensis, L. lamottei, L. ervoides, and L. nigricans) and 13 interspecific advanced lines were evaluated in their transpiration response to high VPD. A large variation in transpiration rate (TR) response to increased VPD was recorded among wild lentil accessions, with 43 accessions exhibiting a breakpoint (BP) in their TR response to increasing VPD, with values ranging from 0.92 kPa to 3.38 kPa under greenhouse conditions. Ten genotypes for the interspecific advanced lines displayed a BP with an average of 1.95 kPa, much lower than previously reported for cultivated lentil. Results from field experiments suggest that the TRlim trait (BP=0.97 kPa) positively affected yield and yield-related parameters during the years with late-season water stress. The selection of TRlim genotypes for high VPD environments could improve lentil productivity in drought-prone areas.

Genome Size Variation across a Cypriot Fabeae Tribe Germplasm Collection

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

Assessment of the Genetic Diversity and Population Structure of the Peruvian Andean Legume, Tarwi (Lupinus mutabilis), with High Quality SNPs

Lupinus mutabilis Sweet (Fabaceae), “tarwi” or “chocho”, is an important grain legume in the Andean region. In Peru, studies on tarwi have mainly focused on morphological features; however, they have not been molecularly characterized. Currently, it is possible to explore the genetic parameters of plants with reliable and modern methods such as genotyping by sequencing (GBS). Here, for the first time, we used single nucleotide polymorphism (SNP) markers to infer the genetic diversity and population structure of 89 accessions of tarwi from nine Andean regions of Peru. A total of 5922 SNPs distributed along all chromosomes of tarwi were identified. STRUCTURE analysis revealed that this crop is grouped into two clusters. A dendrogram was generated using the UPGMA clustering algorithm and, like the principal coordinate analysis (PCoA), it showed two groups that correspond to the geographic origin of the tarwi samples. AMOVA showed a reduced variation between clusters (7.59%) and indicated that variability within populations is 92.41%. Population divergence (Fst) between clusters 1 and 2 revealed low genetic difference (0.019). We also detected a negative Fis for both populations, demonstrating that, like other Lupinus species, tarwi also depends on cross-pollination. SNP markers were powerful and effective for the genotyping process in this germplasm. We hope that this information is the beginning of the path towards a modern genetic improvement and conservation strategies of this important Andean legume.

The Prospects of gene introgression from crop wild relatives into cultivated lentil for climate change mitigation

Crop wild relatives (CWRs), landraces and exotic germplasm are important sources of genetic variability, alien alleles, and useful crop traits that can help mitigate a plethora of abiotic and biotic stresses and crop yield reduction arising due to global climatic changes. In the pulse crop genus Lens, the cultivated varieties have a narrow genetic base due to recurrent selections, genetic bottleneck and linkage drag. The collection and characterization of wild Lens germplasm resources have offered new avenues for the genetic improvement and development of stress-tolerant, climate-resilient lentil varieties with sustainable yield gains to meet future food and nutritional requirements. Most of the lentil breeding traits such as high-yield, adaptation to abiotic stresses and resistance to diseases are quantitative and require the identification of quantitative trait loci (QTLs) for marker assisted selection and breeding. Advances in genetic diversity studies, genome mapping and advanced high-throughput sequencing technologies have helped identify many stress-responsive adaptive genes, quantitative trait loci (QTLs) and other useful crop traits in the CWRs. The recent integration of genomics technologies with plant breeding has resulted in the generation of dense genomic linkage maps, massive global genotyping, large transcriptomic datasets, single nucleotide polymorphisms (SNPs), expressed sequence tags (ESTs) that have advanced lentil genomic research substantially and allowed for the identification of QTLs for marker-assisted selection (MAS) and breeding. Assembly of lentil and its wild species genomes (~4Gbp) opens up newer possibilities for understanding genomic architecture and evolution of this important legume crop. This review highlights the recent strides in the characterization of wild genetic resources for useful alleles, development of high-density genetic maps, high-resolution QTL mapping, genome-wide studies, MAS, genomic selections, new databases and genome assemblies in traditionally bred genus Lens for future crop improvement amidst the impending global climate change.

Conventional and new-breeding technologies for improving disease resistance in lentil (Lens culinaris Medik)

Lentil, an important cool season food legume, is a rich source of easily digestible protein, folic acid, bio-available iron, and zinc nutrients. Lentil grows mainly as a sole crop in the winter after harvesting rice in South Asia. However, the annual productivity is low due to its slow growth during the early phase, competitive weed infestation, and disease outbreaks during the crop growth period. Disease resistance breeding has been practiced for a long time to enhance resistance to various diseases. Often the sources of resistance are available in wild crop relatives. Thus, wide hybridization and the ovule rescue technique have helped to introgress the resistance trait into cultivated lentils. Besides hybridization, induced mutagenesis contributed immensely in creating variability for disease tolerance, and several disease-resistant mutant lines have been developed. However, to overcome the limitations of traditional breeding approaches, advancement in molecular marker technologies, and genomics has helped to develop disease-resistant and climate-resilient lentil varieties with more precision and efficiency. This review describes types of diseases, disease screening methods, the role of conventional and new breeding technologies in alleviating disease-incurred damage and progress toward making lentil varieties more resilient to disease outbreaks under the shadow of climate change.

Unlocking the hidden variation from wild repository for accelerating genetic gain in legumes

The fluctuating climates, rising human population, and deteriorating arable lands necessitate sustainable crops to fulfil global food requirements. In the countryside, legumes with intriguing but enigmatic nitrogen-fixing abilities and thriving in harsh climatic conditions promise future food security. However, breaking the yield plateau and achieving higher genetic gain are the unsolved problems of legume improvement. Present study gives emphasis on 15 important legume crops, i.e., chickpea, pigeonpea, soybean, groundnut, lentil, common bean, faba bean, cowpea, lupin, pea, green gram, back gram, horse gram, moth bean, rice bean, and some forage legumes. We have given an overview of the world and India's area, production, and productivity trends for all legume crops from 1961 to 2020. Our review article investigates the importance of gene pools and wild relatives in broadening the genetic base of legumes through pre-breeding and alien gene introgression. We have also discussed the importance of integrating genomics, phenomics, speed breeding, genetic engineering and genome editing tools in legume improvement programmes. Overall, legume breeding may undergo a paradigm shift once genomics and conventional breeding are integrated in the near future.

Diploid chromosome-level reference genome and population genomic analyses provide insights into Gypenoside biosynthesis and demographic evolution of Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum (Thunb.) Makino is a perennial creeping herbaceous plant in the family Cucurbitaceae, which has great medicinal value and commercial potential, but urgent conservation efforts are needed due to the gradual decreases and fragmented distribution of its wild populations. Here, we report the high-quality diploid chromosome-level genome of G. pentaphyllum obtained using a combination of next-generation sequencing short reads, Nanopore long reads, and Hi-C sequencing technologies. The genome is anchored to 11 pseudo-chromosomes with a total size of 608.95 Mb and 26 588 predicted genes. Comparative genomic analyses indicate that G. pentaphyllum is estimated to have diverged from Momordica charantia 60.7 million years ago, with no recent whole-genome duplication event. Genomic population analyses based on genotyping-by-sequencing and ecological niche analyses indicated low genetic diversity but a strong population structure within the species, which could classify 32 G. pentaphyllum populations into three geographical groups shaped jointly by geographic and climate factors. Furthermore, comparative transcriptome analyses showed that the genes encoding enzyme involved in gypenoside biosynthesis had higher expression levels in the leaves and tendrils. Overall, the findings obtained in this study provide an effective molecular basis for further studies of demographic genetics, ecological adaption, and systematic evolution in Cucurbitaceae species, as well as contributing to molecular breeding, and the biosynthesis and biotransformation of gypenoside.

Complete chloroplast genome sequence of Lens ervoides and comparison to Lens culinaris

Lens is a member of the Papilionoideae subfamily of Fabaceae and is generally used as a source of vegetable protein as part of human diets in many regions worldwide. Chloroplast (cp) genomes are highly active genetic components of plants and can be utilized as molecular markers for various purposes. As one of the wild lentil species, the Lens ervoides cp genome has been sequenced for the first time in this study using next-generation sequencing. The de novo assembly of the cp genome resulted in a single 122,722 bp sequence as two separate coexisting structural haplotypes with similar lengths. Results indicated that the cp genome of L. ervoides belongs to the inverted repeat lacking clade. Several noteworthy divergences within the coding regions were observed in ndhB, ndhF, rbcL, rpoC2, and ycf2 genes. Analysis of relative synonymous codon usage showed that certain genes, psbN, psaI, psbI, psbE, psbK, petD, and ndhC, preferred using biased codons more often and therefore might have elevated expression and translation efficiencies. Overall, this study exhibited the divergence level between the wild-type and cultured lentil cp genomes and pointed to certain regions that can be utilized as distinction markers for various goals.

RNA-Seq and Gene Ontology Analysis Reveal Differences Associated With Low R/FR-Induced Shade Responses in Cultivated Lentil and a Wild Relative

Lentil is an important pulse crop not only because of its high nutrient value but also because of its ecological advantage in a sustainable agricultural system. Our previous work showed that the cultivated lentil and wild lentil germplasm respond differently to light environments, especially to low R/FR-induced shade conditions. Little is known about how cultivated and wild lentils respond to shade at the level of gene expression and function. In this study, transcriptomic profiling of a cultivated lentil (Lupa, L. culinaris) and a wild lentil (BGE 016880, L. orientalis) at several growth stages is presented. De novo transcriptomes were assembled for both genotypes, and differential gene expression analysis and gene ontology enrichment analysis were performed. The transcriptomic resources generated in this study provide fundamental information regarding biological processes and genes associated with shade responses in lentils. BGE 016880 and Lupa shared a high similarity in their transcriptomes; however, differential gene expression profiles were not consistent between these two genotypes. The wild lentil BGE 016880 had more differentially expressed genes than the cultivated lentil Lupa. Upregulation of genes involved in gibberellin, brassinosteroid, and auxin synthesis and signaling pathways, as well as cell wall modification, in both genotypes explains their similarity in stem elongation response under the shade. Genes involved in jasmonic acid and flavonoid biosynthesis pathways were downregulated in BGE 016880 only, and biological processes involved in defense responses were significantly enriched in the wild lentil BGE 016880 only. Downregulation of WRKY and MYB transcription factors could contribute to the reduced defense response in BGE 016880 but not in Lupa under shade conditions. A better understanding of shade responses of pulse crop species and their wild relatives will play an important role in developing genetic strategies for crop improvement in response to changes in light environments.

Supporting in situ conservation of the genetic diversity of crop wild relatives using genomic technologies

The last decade has witnessed huge technological advances in genomics, particularly in DNA sequencing. Here, we review the actual and potential application of genomics in supporting in situ conservation of crop wild relatives (CWRs). In addition to helping in prioritization of protection of CWR taxa and in situ conservation sites, genome analysis is allowing the identification of novel alleles that need to be prioritized for conservation. Genomics is enabling the identification of potential sources of important adaptive traits that can guide the establishment or enrichment of in situ genetic reserves. Genomic tools also have the potential for developing a robust framework for monitoring and reporting genome‐based indicators of genetic diversity changes associated with factors such as land use or climate change. These tools have been demonstrated to have an important role in managing the conservation of populations, supporting sustainable access and utilization of CWR diversity, enhancing accelerated domestication of new crops and forensic genomics thus preventing misappropriation of genetic resources. Despite this great potential, many policy makers and conservation managers have failed to recognize and appreciate the need to accelerate the application of genomics to support the conservation and management of biodiversity in CWRs to underpin global food security. Funding and inadequate genomic expertise among conservation practitioners also remain major hindrances to the widespread application of genomics in conservation.

Protein Biofortification in Lentils (Lens culinaris Medik.) Toward Human Health

Lentil (Lens culinaris Medik.) is a nutritionally dense crop with significant quantities of protein, low-digestible carbohydrates, minerals, and vitamins. The amino acid composition of lentil protein can impact human health by maintaining amino acid balance for physiological functions and preventing protein-energy malnutrition and non-communicable diseases (NCDs). Thus, enhancing lentil protein quality through genetic biofortification, i.e., conventional plant breeding and molecular technologies, is vital for the nutritional improvement of lentil crops across the globe. This review highlights variation in protein concentration and quality across Lens species, genetic mechanisms controlling amino acid synthesis in plants, functions of amino acids, and the effect of antinutrients on the absorption of amino acids into the human body. Successful breeding strategies in lentils and other pulses are reviewed to demonstrate robust breeding approaches for protein biofortification. Future lentil breeding approaches will include rapid germplasm selection, phenotypic evaluation, genome-wide association studies, genetic engineering, and genome editing to select sequences that improve protein concentration and quality.

Integrating Genomic and Phenomic Approaches to Support Plant Genetic Resources Conservation and Use

Plant genebanks provide genetic resources for breeding and research programs worldwide. These programs benefit from having access to high-quality, standardized phenotypic and genotypic data. Technological advances have made it possible to collect phenomic and genomic data for genebank collections, which, with the appropriate analytical tools, can directly inform breeding programs. We discuss the importance of considering genebank accession homogeneity and heterogeneity in data collection and documentation. Citing specific examples, we describe how well-documented genomic and phenomic data have met or could meet the needs of plant genetic resource managers and users. We explore future opportunities that may emerge from improved documentation and data integration among plant genetic resource information systems.

Genetic diversity and main functional composition of Lingzhi strains from main producing areas in China

In this study, we used genotyping by sequencing (GBS) to examine the genetic diversity of 22 strains of Lingzhi and the quality differences in 15 fruit bodies of Lingzhi from different Chinese regions. The phylogenetic trees of 22 strains were constructed based on ITS (Internal transcribed spacer) and SNP (single nucleotide polymorphism). Moisture, ash, water-soluble extracts, alcohol-soluble extracts, polysaccharides, and triterpenoids from 15 fruit bodies of Lingzhi were detected and analyzed based on Chinese Pharmacopoeia and the US Pharmacopoeia references. Moreover, the monosaccharide composition of polysaccharides was studied using PMP-HPLC, and the effect of polysaccharides on the proliferation rate of splenocytes was investigated in vitro. The identification results of these strains by the phylogenetic trees which were constructed based on ITS sequences and SNPs showed that most of the strains applied in the main producing areas of Lingzhi in China were accurate except for a few inaccurate strains. The moisture, ash, water and alcohol soluble extractive, polysaccharide and triterpenoid content of all samples were meet the requirements of the Chinese Pharmacopoeia, while the polysaccharide and triterpenoid content of less than half of the samples meet the requirements of the U.S. Pharmacopoeia. The polysaccharide extracted from these samples have different effects on the proliferation rate of spleen cells. To sum up, this is the first study that reported on the differences in Lingzhi strains from the main producing areas in China. The quality of some fruit bodies did not meet the pharmacopeia requirements, and wrong strains were used in some production areas; thus, strains should be given special attention before legal processing.

Diversity in Surface Microstructures of Trichomes, Epidermal Cells, and Stomata in Lentil Germplasm

To develop crops capable of withstanding challenges posed by climate change, breeding strategies must focus on addressing multiple stresses occurring concurrently in plants. Leaf epidermal structures such as trichomes, stomata, and epidermal cells play an important role in mediating plant defense and could be essential traits that impart wide-ranging tolerance to biotic and abiotic stresses. Consequently, it is important to inform on the underlying diversity in these traits in lentil germplasm (Lens spp.). In this study, we characterized foliar microstructures of 12 genotypes belonging to seven wild and cultivated Lens species. We performed scanning electron microscopy on leaflet and pod surfaces for their qualitative characterization. For quantitative characterization, we observed surface imprints via light microscopy and quantified trichome density (TD), trichome length (TL), stomatal density (SD), epidermal cell density (ECD), and stomatal index (SI) on adaxial and abaxial leaflet surfaces for each genotype. We also assessed the heritability of trichome traits by evaluating interspecific recombinant inbred lines (RILs) derived from the cross Lens culinaris CDC Redberry × Lens tomentosus IG 72805. Comparing foliar microstructures, we found that TD and TL varied widely among cultivated and wild lentil genotypes. However, in most lentil genotypes, the adaxial leaflet surface had lower TD and longer trichomes compared to the abaxial surface. Pubescence on pods comprised five major phenotypes: no trichomes or glabrous pods, very short trichomes at low density, short trichomes at high density, medium-length trichomes at high density, and long trichomes at high density. Leaves of all species were amphistomatous, and SI, SD, and ECD were all higher on the adaxial compared to the abaxial surface. Adaxial surfaces had slightly sunken stomata, which might be an adaptive trait to conserve water. Quantifying TD and TL on the leaflets of interspecific RILs revealed transgressive segregation of these traits, suggesting that TD and TL are quantitative in nature. While taxonomic implications of this study are limited, a detailed description of agronomically relevant morphophysiological traits presented in this paper along with the mode of inheritance of trichomes may serve as a resource for scientists developing lentil adapted to concurrent biotic and abiotic stresses of the future.

Intelligent Characterization of Lentil Genetic Resources: Evolutionary History, Genetic Diversity of Germplasm, and the Need for Well-Represented Collections

The genetic and phenotypic characterization of crops allows us to elucidate their evolutionary and domestication history, the genetic basis of important traits, and the use of variation present in landraces and wild relatives to enhance resilience. In this context, we aim to provide an overview of the main genetic resources developed for lentil and their main outcomes, and to suggest protocols for continued work on this important crop. Lens culinaris is the third-most-important cool-season grain and its use is increasing as a quick-cooking, nutritious, plant-based source of protein. L. culinaris was domesticated in the Fertile Crescent, and six additional wild taxa (L. orientalis, L. tomentosus, L. odemensis, L. lamottei, L. ervoides, and L. nigricans) are recognized. Numerous genetic diversity studies have shown that wild relatives present high levels of genetic variation and provide a reservoir of alleles that can be used for breeding programs. Furthermore, the integration of genetics/genomics and breeding techniques has resulted in identification of quantitative trait loci and genes related to attributes of interest. Genetic maps, massive genotyping, marker-assisted selection, and genomic selection are some of the genetic resources generated and applied in lentil. In addition, despite its size (∼4 Gbp) and complexity, the L. culinaris genome has been assembled, allowing a deeper understanding of its architecture. Still, major knowledge gaps exist in lentil, and a deeper understanding and characterization of germplasm resources, including wild relatives, is critical to lentil breeding and improvement. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Recording of lentil seed descriptors Basic Protocol 2: Lentil seed imaging Basic Protocol 3: Lentil seed increase Basic Protocol 4: Recording of primary lentil seed INCREASE descriptors.

The complete chloroplast genome of Abeliophyllum distichum f. lilacinum Nakai (Oleaceae) from the Chungbuk Province, Korea

The chloroplast genome of Abeliophyllum distichum f. lilacinum Nakai, classified to a monotypic in this genus, and an endemic species in Korea, was sequenced to understand the genetic differences among intraspecies and cultivars of A. distichum. The chloroplast genome length is 156,015 bp (GC ratio is 37.8%) and has a typical quadripartite structure: 86,779 bp large single copy (35.8%) and 17,828 bp small single copy (31.9%) regions separated by two 25,704 bp inverted repeat (43.2%) regions. The genome encodes for 133 genes (88 protein-coding genes, eight rRNAs, and 37 tRNAs). Six to 99 SNPs and seven to 18 INDEL regions (19 bp to 72 bp) were identified against available chloroplast genomes of A. distichum. Phylogenetic trees show that A. distichum f. lilacinum is clustered with the Dae Ryun cultivar which has a larger fruit body. Our analyses suggest additional research, such as Genotyping-By-Sequencing, for understanding relationship between morphology and genotype of A. distichum.

Application of Genomics to Understand Salt Tolerance in Lentil

Soil salinity is a major abiotic stress, limiting lentil productivity worldwide. Understanding the genetic basis of salt tolerance is vital to develop tolerant varieties. A diversity panel consisting of 276 lentil accessions was screened in a previous study through traditional and image-based approaches to quantify growth under salt stress. Genotyping was performed using two contrasting methods, targeted (tGBS) and transcriptome (GBS-t) genotyping-by-sequencing, to evaluate the most appropriate methodology. tGBS revealed the highest number of single-base variants (SNPs) (c. 56,349), and markers were more evenly distributed across the genome compared to GBS-t. A genome-wide association study (GWAS) was conducted using a mixed linear model. Significant marker-trait associations were observed on Chromosome 2 as well as Chromosome 4, and a range of candidate genes was identified from the reference genome, the most plausible being potassium transporters, which are known to be involved in salt tolerance in related species. Detailed mineral composition performed on salt-treated and control plant tissues revealed the salt tolerance mechanism in lentil, in which tolerant accessions do not transport Na⁺ ions around the plant instead localize within the root tissues. The pedigree analysis identified two parental accessions that could have been the key sources of tolerance in this dataset.

The History of Lentil (Lens culinaris subsp. culinaris) Domestication and Spread as Revealed by Genotyping-by-Sequencing of Wild and Landrace Accessions

Protein-rich legumes accompanied carbohydrate-rich cereals since the beginning of agriculture and yet their domestication history is not as well understood. Lentil (Lens culinaris Medik. subsp. culinaris) was first cultivated in Southwest Asia (SWA) 8000-10,000 years ago but archeological evidence is unclear as to how many times it may have been independently domesticated, in which SWA region(s) this may have happened, and whether wild species within the Lens genus have contributed to the cultivated gene pool. In this study, we combined genotyping-by-sequencing (GBS) of 190 accessions from wild (67) and domesticated (123) lentils from the Old World with archeological information to explore the evolutionary history, domestication, and diffusion of lentils to different environments. GBS led to the discovery of 87,647 single-nucleotide polymorphisms (SNPs), which allowed us to infer the phylogeny of genus Lens. We confirmed previous studies proposing four groups within it. The only gene flow detected was between cultivated varieties and their progenitor (L. culinaris subsp. orientalis) albeit at very low levels. Nevertheless, a few putative hybrids or naturalized cultivars were identified. Within cultivated lentil, we found three geographic groups. Phylogenetics, population structure, and archeological data coincide in a scenario of protracted domestication of lentils, with two domesticated gene pools emerging in SWA. Admixed varieties are found throughout their range, suggesting a relaxed selection process. A small number of alleles involved in domestication and adaptation to climatic variables were identified. Both novel mutation and selection on standing variation are presumed to have played a role in adaptation of lentils to different environments. The results presented have implications for understanding the process of plant domestication (past), the distribution of genetic diversity in germplasm collections (present), and targeting genes in breeding programs (future).

Effects of Lentil Genotype on the Colonization of Beneficial Trichoderma Species and Biocontrol of Aphanomyces Root Rot

Trichoderma species are opportunistic plant symbionts that are common in the root and rhizosphere ecosystems. Many Trichoderma species may enhance plant growth, nutrient acquisition, and disease resistance, and for these reasons, they are widely used in agriculture as biofertilizers or biocontrol agents. Host plant genotype and other microorganisms, such as root pathogens, may influence the efficacy of Trichoderma inoculants. Aphanomyces euteiches is an important soil-borne oomycete in western Canada that causes root rot in legume crops such as lentil and pea, and there is not yet any significantly resistant varieties or effective treatments available to control the disease. In this study, the composition of root-associated fungal communities and the abundance of Trichoderma species, T. harzianum strain T-22 and T. virens strain G41, was determined in the roots of eight Lens genotypes based on internal transcribed spacer (ITS) Illumina MiSeq paired-end sequencing, both in the presence and the absence of the root rot pathogen Aphanomyces euteiches. Biocontrol effects of T. harzianum on A. euteiches was also examined. Significant genotypic variations were observed in the composition of root-associated fungal communities and the abundance of the different Trichoderma species in the lentil roots. The presence of A. euteiches altered the composition of Trichoderma found associated to the lentil genotypes. Biocontrol of A. euteiches by T. harzianum T22 species was observed in vitro and positive correlations between the abundance of Trichoderma and plant root and shoot biomass were observed in vivo. These findings revealed that lentil genotype and infection by the phytopathogen A. euteiches greatly influenced the colonization of root-associated fungi and the abundance of the Trichoderma species, as well as the effect on plant growth promotion. The multipartite interactions observed among lentil genotypes, Trichoderma species and A. euteiches suggest possibilities to select compatible host-beneficial microbe combinations in lentil breeding programs and to develop application strategies to harness the beneficial effects of Trichoderma inoculants in sustainable crop production systems.

Characterization of Genetic and Allelic Diversity Amongst Cultivated and Wild Lentil Accessions for Germplasm Enhancement

Intensive breeding of cultivated lentil has resulted in a relatively narrow genetic base, which limits the options to increase crop productivity through selection. Assessment of genetic diversity in the wild gene pool of lentil, as well as characterization of useful and novel alleles/genes that can be introgressed into elite germplasm, presents new opportunities and pathways for germplasm enhancement, followed by successful crop improvement. In the current study, a lentil collection consisting of 467 wild and cultivated accessions that originated from 10 diverse geographical regions was assessed, to understand genetic relationships among different lentil species/subspecies. A total of 422,101 high-confidence SNP markers were identified against the reference lentil genome (cv. CDC Redberry). Phylogenetic analysis clustered the germplasm collection into four groups, namely, Lens culinaris/Lens orientalis, Lens lamottei/Lens odemensis, Lens ervoides, and Lens nigricans. A weak correlation was observed between geographical origin and genetic relationship, except for some accessions of L. culinaris and L. ervoides. Genetic distance matrices revealed a comparable level of variation within the gene pools of L. culinaris (Nei’s coefficient 0.01468–0.71163), L. ervoides (Nei’s coefficient 0.01807–0.71877), and L. nigricans (Nei’s coefficient 0.02188–1.2219). In order to understand any genic differences at species/subspecies level, allele frequencies were calculated from a subset of 263 lentil accessions. Among all cultivated and wild lentil species, L. nigricans exhibited the greatest allelic differentiation across the genome compared to all other species/subspecies. Major differences were observed on six genomic regions with the largest being on Chromosome 1 (c. 1 Mbp). These results indicate that L. nigricans is the most distantly related to L. culinaris and additional structural variations are likely to be identified from genome sequencing studies. This would provide further insights into evolutionary relationships between cultivated and wild lentil germplasm, for germplasm improvement and introgression.

Dissecting the Genetic Architecture of Aphanomyces Root Rot Resistance in Lentil by QTL Mapping and Genome-Wide Association Study

Lentil (Lens culinaris Medikus) is an important source of protein for people in developing countries. Aphanomyces root rot (ARR) has emerged as one of the most devastating diseases affecting lentil production. In this study, we applied two complementary quantitative trait loci (QTL) analysis approaches to unravel the genetic architecture underlying this complex trait. A recombinant inbred line (RIL) population and an association mapping population were genotyped using genotyping by sequencing (GBS) to discover novel single nucleotide polymorphisms (SNPs). QTL mapping identified 19 QTL associated with ARR resistance, while association mapping detected 38 QTL and highlighted accumulation of favorable haplotypes in most of the resistant accessions. Seven QTL clusters were discovered on six chromosomes, and 15 putative genes were identified within the QTL clusters. To validate QTL mapping and genome-wide association study (GWAS) results, expression analysis of five selected genes was conducted on partially resistant and susceptible accessions. Three of the genes were differentially expressed at early stages of infection, two of which may be associated with ARR resistance. Our findings provide valuable insight into the genetic control of ARR, and genetic and genomic resources developed here can be used to accelerate development of lentil cultivars with high levels of partial resistance to ARR.

Evaluation and identification of wild lentil accessions for enhancing genetic gains of cultivated varieties

Domesticated lentil has a relatively narrow genetic base globally and most released varieties are susceptible to severe biotic and abiotic stresses. The crop wild relatives could provide new traits of interest for tailoring novel germplasm and cultivated lentil improvement. The primary objective of this study was to evaluate wild lentil accessions for identification of economically viable agro-morphological traits and resistance against major biotic stresses. The study has revealed substantial variations in seed yield and its important component characters. Further, the diversity analysis of wild accessions showed two major clusters which were bifurcated into sub-clusters, thereby suggesting their wider genetic divergence. However, principal component analysis exhibited that seed yield plant-1, number of seeds plant-1, number of pods plant-1, harvest index and biological yield plant-1 contributed significantly to the total genetic variation assessed in wild lentil taxa. Moreover, some of the wild accessions collected from Syria and Turkey regions showed resistance against more than one disease indicating rich diversity of lentil genetic resources. The identification of most promising genotypes carrying resistance against major biotic stresses could be utilized in the cultivated or susceptible varieties of lentil for enhancing genetic gains. The study has also identified some trait specific accessions, which could also be taken into the consideration while planning distant hybridization in lentil.

Spatiotemporal Differentiation of Alpine Butterfly Parnassius glacialis (Papilionidae: Parnassiinae) in China: Evidence from Mitochondrial DNA and Nuclear Single Nucleotide Polymorphisms

The Apollo butterfly, Parnassius glacialis, is one of the most charming members of its genus and includes two subspecies locally distributed in montane areas of south-central China and Japan. In this study, we investigated the genetic structure and demographic history of P. glacialis by analyzing partial sequences of four mitochondrial genes and nuclear single nucleotide polymorphisms (SNPs) via genotyping-by-sequencing (GBS) of samples from nearly the entire known distributional range in China. The mitochondrial DNA (mtDNA) data demonstrated that a total of 39 haplotypes were present, and the species was estimated to have diverged about 0.95 million years ago during the middle Pleistocene transition into two main clades that likely formed during the Kunlun-Huanghe tectonic movement. The two clades then dispersed independently in distinct geographic areas alongside the mountainous routes in central and southern China, most likely driven by the Pleistocene glacial-interglacial cycles. Nuclear SNP analysis was generally congruent with mtDNA results at the individual level. A minor incongruence of genetic structures that was detected between mtDNA and nuclear SNP data from the Laojunshan and Tiantangzhai populations was likely due to secondary contact and male-biased dispersal. Our work demonstrates that complicated dispersal-vicariance evolutionary processes likely led to the current geographic distribution of P. glacialis in China, particularly the uplift of the Qinghai-Tibet Plateau and related climatic oscillations during the Quaternary period.

Genotyping-by-Sequencing Derived Single Nucleotide Polymorphisms Provide the First Well-Resolved Phylogeny for the Genus Triticum (Poaceae)

Wheat (Triticum spp.) has been an important staple food crop for mankind since the beginning of agriculture. The genus Triticum L. is composed of diploid, tetraploid, and hexaploid species, majority of which have not yet been discriminated clearly, and hence their phylogeny and classification remain unresolved. Genotyping-by-sequencing (GBS) is an easy and affordable method that allows us to generate genome-wide single nucleotide polymorphism (SNP) markers. In this study, we used GBS to obtain SNPs covering all seven chromosomes from 283 accessions of Triticum-related genera. After filtering low-quality and redundant SNPs based on haplotype information, the GBS assay provided 14,188 high-quality SNPs that were distributed across the A (71%), B (26%), and D (2.4%) genomes. Cluster analysis and discriminant analysis of principal components (DAPC) allowed us to distinguish six distinct groups that matched well with Triticum species complexity. We constructed a Bayesian phylogenetic tree using 14,188 SNPs, in which 17 Triticum species and subspecies were discriminated. Dendrogram analysis revealed that the polyploid wheat species could be divided into groups according to the presence of A, B, D, and G genomes with strong nodal support and provided new insight into the evolution of spelt wheat. A total of 2,692 species-specific SNPs were identified to discriminate the common (T. aestivum) and durum (T. turgidum) wheat cultivar and landraces. In principal component analysis grouping, the two wheat species formed individual clusters and the SNPs were able to distinguish up to nine groups of 10 subspecies. This study demonstrated that GBS-derived SNPs could be used efficiently in genebank management to classify Triticum species and subspecies that are very difficult to distinguish by their morphological characters.

Genomic Analyses Provide Insights Into the Evolutionary History and Genetic Diversity of Auricularia Species

Species in the genus Auricularia play important roles for people’s food and nutrition especially Auricularia cornea and A. heimuer. To understand their evolutionary history, genome structure, and population-level genetic variation, we performed a high-quality genome sequencing of Auricularia cornea and the corresponding comparative genomic analysis. The genome size of A. cornea was similar to Auricularia subglabra, but 1.5 times larger than that of A. heimuer. Several factors were responsible for genome size variation including gene numbers, repetitive elements, and gene lengths. Phylogenomic analysis revealed that the estimated divergence time between A. heimuer and other Auricularia is ∼79.1 million years ago (Mya), while the divergence between A. cornea and A. subglabra occurred in ∼54.8 Mya. Population genomic analysis also provided insight into the demographic history of A. cornea and A. heimuer, indicating that their populations fluctuated over time with global climate change during Marine Isotope Stage 5-2. Moreover, despite the highly similar external morphologies of A. cornea and A. heimuer, their genomic properties were remarkably different. The A. cornea genome only shared 14% homologous syntenic blocks with A. heimuer and possessed more genes encoding carbohydrate-active enzymes and secondary metabolite biosynthesis proteins. The cross-taxa transferability rates of simple sequence repeat (SSR) and insertion or deletion (InDel) markers within the genus Auricularia were also lower than that previously observed for species within the same genus. Taken together, these results indicate a high level of genetic differentiation between these two Auricularia species. Consequently, our study provides new insights into the genomic evolution and genetic differentiation of Auricularia species that will facilitate future genetic breeding.

Genotyping by Sequencing of Cultivated Lentil (Lens culinaris Medik.) Highlights Population Structure in the Mediterranean Gene Pool Associated With Geographic Patterns and Phenotypic Variables

Cultivated lentil (Lens culinaris Medik.) is one of the oldest domesticated crops and one of the most important grain legumes worldwide. The Mediterranean Basin holds large part of lentil biodiversity; however, no genetic structure was defined within the Mediterranean gene pool. In this study, we used high-throughput genotyping by sequencing to resolve the genetic structure of the Mediterranean ex situ lentil collection held at the Italian National Research Council. Sequencing of a 188-plex genotyping-by-sequencing library and bioinformatics treatment of data yielded 6,693 single nucleotide polymorphisms. Analysis of nonredundant genotypes with nonparametric and parametric methods highlighted the occurrence of five highly differentiated genetic clusters. Clustering could be related to geographic patterns and phenotypic traits, indicating that post-domestication routes introducing cultivation in Mediterranean countries and selection were major forces shaping lentil population structure. The estimation of the fixation index F_ST at individual single nucleotide polymorphism loci allowed the identification of distinctive alleles across clusters, suggesting the possibility to set up molecular keys for the assignment of lentil germplasm to specific genetic groups. Finally, significant associations between markers and phenotypic data were identified. Overall, the results of this study are of major importance for lentil conservation genetics and breeding and provide insights on the lentil evolutionary history.

Seed Protein of Lentils: Current Status, Progress, and Food Applications

Grain legumes are widely recognized as staple sources of dietary protein worldwide. Lentil seeds are an excellent source of plant-based proteins and represent a viable alternative to animal and soybean proteins for food processing formulations. Lentil proteins provide not only dietary amino acids but are also a source of bioactive peptides that provide health benefits. This review focuses on the current knowledge of seed protein, extraction and isolation methods, bioactive peptides, and food applications of lentil protein. Lentil is the most rapidly expanding crop for direct human consumption, and has potential for greater impact as a protein source for food processing applications. Improvements in lentil protein quality, amino acid composition, and processing fractions will enhance the nutritional quality of this rapidly expanding crop globally.

Genetic diversity and population structure of Ethiopian Capsicum germplasms

We established a collection of 142 Capsicum genotypes from different geographical areas of Ethiopia with the aim of capturing genetic diversity. Morphological traits and high-resolution melting analysis distinguished one Capsicum baccatum, nine Capsicum frutescens and 132 Capsicum annuum accessions in the collection. Measurement of plant growth parameters revealed variation between germplasms in parameters including plant height, stem thickness, internode length, number of side branches, fruit width, and fruit length. Broad-sense heritability was maximum for fruit weight, followed by length and width of leaves. We used genotyping by sequencing (GBS) to identify single-nucleotide polymorphisms (SNPs) in the panel of 142 Capsicum germplasms and found 2,831,791 genome-wide SNP markers. Among these, we selected 53,284 high-quality SNPs and used them to estimate the level of genetic diversity, population structure, and phylogenetic relationships. From model-based ancestry analysis, the phylogenetic tree and principal-coordinate analysis (PCoA), we identified two distinct genetic populations: one comprising 132 C. annuum accessions and the other comprising the nine C. frutescens accessions. GWAS analysis detected 509 SNP markers that were significantly associated with fruit-, stem- and leaf-related traits. This is the first comprehensive report of the analysis of genetic variation in Ethiopian Capsicum species involving a large number of accessions. The results will help breeders utilize the germplasm collection to improve existing commercial cultivars.

Construction of a high-density interspecific (Lens culinaris x L. odemensis) genetic map based on functional markers for mapping morphological and agronomical traits, and QTLs affecting resistance to Ascochyta in lentil

Usage of high-throughput sequencing approaches allow for the generation and characterization of reference transcriptome datasets that support gene-based marker discovery, which in turn can be used to build genetic maps among other purposes. We have obtained a transcriptome assembly including 49,453 genes for the lentil (Lens culinaris Medik.) cultivar Alpo using RNAseq methodology. This transcriptome was used as reference to obtain 6,306 quality polymorphic markers (SNPs and short indels) analyzing genotype data from a RIL population at F7 generation derived from the interspecific cross between L. culinaris cv. Alpo and L. odemensis accession ILWL235. L. odemensis is a wild species included in the secondary gene pool and can be used as a source for gene introgression in lentil breeding programs. Marker data were used to construct the first genetic interspecific map between these two species. This linkage map has been used to precisely identify regions of the CDC-Redberry lentil draft genome in which the candidate genes for some qualitative traits (seed coat spotting pattern, flower color, and stem pigmentation) could be located. The genome regions corresponding to a significant single quantitative trait locus (QTL) controlling "time to flowering" located in chromosome 6 and three QTLs regulating seed size and positioned in chromosomes 1 and 5 (two QTLs) were also identified. Significant QTLs for Ascochyta blight resistance in lentil were mapped to chromosome 6 in the genome region or close to it where QTLs for Ascochyta blight resistance have previously been reported.

Polyphenolic Composition of Lentil Roots in Response to Infection by Aphanomyces euteiches

Polyphenols comprise the largest group of plant secondary metabolites and have critical roles in plant physiology and response to the biotic and abiotic environment. Changes in the content of polyphenols in the root extracts and root tissues of wild (Lens ervoides) and cultivated (Lens culinaris) lentil genotypes were examined in response to infection by Aphanomyces euteiches using liquid chromatography mass spectrometry (LC-MS). Genotype, infection and their interaction determined the composition of polyphenols in lentil roots. The levels of several polyphenols were lower in the root extract of the low-tannin genotype L. culinaris ZT-4 compared to L. ervoides L01-827A. Kaempferol derivatives including kaempferol dirutinoside and kaempferol 3-robinoside 7-rhamnoside were more concentrated in the healthy root tissues of L. ervoides L01-827A than in L. culinaris genotypes. Infection increased the concentration of kaempferol, apigenin, and naringenin in the root tissues of all genotypes, but had no effect on some polyphenols in the low-tannin genotype L. culinaris ZT-4. The concentrations of apigenin, naringenin, apigenin 4-glucoside, naringenin7-rutinoside, diosmetin, and hesperetin 7-rutinoside were higher in the infected root tissues of L. ervoides L01-827A compared with the L. culinaris genotypes. Organic acids including coumaric acid, vanillic acid, 4-aminosalicylic acid, 4-hydroxybenzoic acid, and 3,4-dihydroxybenzoic acid effectively suppressed the in-vitro hyphal growth of A. euteiches. Some of these bioactive polyphenols were more concentrated in roots of L. ervoides L01-827A but were low to undetectable in ZT-4. This study shows that genotypic differences exist in the composition of root polyphenols in lentil, and is related to the response to infection caused by A. euteiches. Polyphenols, particularly the organic acid content could be useful for selection and breeding of lentil genotypes that are resistant to Aphanomyces root rot (ARR) disease.

Capturing variation in Lens (Fabaceae): Development and utility of an exome capture array for lentil

PREMISE OF THE STUDY

Lentil is an important legume crop with reduced genetic diversity caused by domestication bottlenecks. Due to its large and complex genome, tools for reduced representation sequencing are needed. We developed an exome capture array for use in various genetic diversity studies.

METHODS

Based on the CDC Redberry draft genome, we developed an exome capture array using multiple sources of transcript resources. The probes were designed to target not only the cultivated lentil, but also wild species. We assessed the utility of the developed method by applying the generated data set to population structure and phylogenetic analyses.

RESULTS

The data set includes 16 wild lentils and 22 cultivar accessions of lentil. Alignment rates were over 90%, and the genic regions were well represented in the capture array. After stringent filtering, 6.5 million high-quality variants were called, and the data set was used to assess the interspecific relationships within the genus Lens.

DISCUSSION

The developed exome capture array provides large amounts of genomic data to be used in many downstream analyses. The method will have useful applications in marker-assisted breeding programs aiming to improve the quality of cultivated lentil.

Naturally occurring variations in the nod-independent model legume Aeschynomene evenia and relatives: a resource for nodulation genetics

BACKGROUND

Among semi-aquatic species of the legume genus Aeschynomene, some have the unique property of being root and stem-nodulated by photosynthetic Bradyrhizobium lacking the nodABC genes necessary for the production of Nod factors. These species provide an excellent biological system with which to explore the evolution of nodulation in legumes. Among them, Aeschynomene evenia has emerged as a model legume to undertake the genetic dissection of the so-called Nod-independent symbiosis. In addition to the genetic analysis of nodulation on a reference line, natural variation in a germplasm collection could also be surveyed to uncover genetic determinants of nodulation. To this aim, we investigated the patterns of genetic diversity in a collection of 226 Nod-independent Aeschynomene accessions.

RESULTS

A combination of phylogenetic analyses, comprising ITS and low-copy nuclear genes, along with cytogenetic experiments and artificial hybridizations revealed the richness of the Nod-independent Aeschynomene group with the identification of 13 diploid and 6 polyploid well-differentiated taxa. A set of 54 SSRs was used to further delineate taxon boundaries and to identify different genotypes. Patterns of microsatellite diversity also illuminated the genetic basis of the Aeschynomene taxa that were all found to be predominantly autogamous and with a predicted simple disomic inheritance, two attributes favorable for genetics. In addition, taxa displaying a pronounced genetic diversity, notably A. evenia, A. indica and A. sensitiva, were characterized by a clear geographically-based genetic structure and variations in root and stem nodulation.

CONCLUSION

A well-characterized germplasm collection now exists as a major genetic resource to thoroughly explore the natural variation of nodulation in response to different bradyrhizobial strains. Symbiotic polymorphisms are expected to be found notably in the induction of nodulation, in nitrogen fixation and also in stem nodulation. Subsequent genetic analysis and locus mapping will pave the way for the identification of the underlying genes through forward or reverse genetics. Such discoveries will significantly contribute to our understanding of the molecular mechanisms underpinning how some Aeschynomene species can be efficiently nodulated in a Nod-independent fashion.

Root Traits, Nodulation and Root Distribution in Soil for Five Wild Lentil Species and Lens culinaris (Medik.) Grown under Well-Watered Conditions

The efficient use of resources such as water and nutrients by plants is increasingly important as the world population food demand continues to grow. With the increased production of lentil in the temperate zones of North America, improvement in yield needs to be maintained. The use of wild lentil genotypes as sources of genetic diversity for introgression into cultivated lentil is an important breeding strategy, but little is known about their root systems. We evaluated the root systems of five wild lentil species and Lens culinaris under fully watered conditions. Plants were grown in 60 cm tubes containing equal volumes of soil collected from the reconstructed A, B, and C horizons. Significant differences were observed for root traits and fine root distribution between and within species and the proportion of root biomass partitioned into each soil layer was unique for each genotype. We also observed variability in nodule number and nodule shape within and between genotypes. Some genotypes more efficiently used water for either biomass or seed production. The allocation of resources to seed production also varied between genotypes. These observations could have impact on the design of future lentil breeding in the context of strategies for managing changes in rainfall amount and distribution for lentil production ecosystems.

QTL mapping reveals genetic determinants of fungal disease resistance in the wild lentil species Lens ervoides

Lens ervoides, a wild relative of lentil is an important source of allelic diversity for enhancing the genetic resistance of the cultivated species against economically important fungal diseases, such as anthracnose and Stemphylium blight caused by Colletotrichum lentis and Stemphylium botryosum, respectively. To unravel the genetic control underlying resistance to these fungal diseases, a recombinant inbred line (RIL) population (n = 94, F9) originating from a cross between two L. ervoides accessions, L01-827A and IG 72815, was genotyped on the Illumina HiSeq 2500 platform. A total of 289.07 million 100 bp paired-end reads were generated, giving an average 7.53-fold genomic coverage to the RILs and identifying 2,180 high-quality SNPs that assembled in 543 unique haplotypes. Seven linkage groups were resolved among haplotypes, equal to the haploid chromosome number in L. ervoides. The genetic map spanned a cumulative distance of 740.94 cM. Composite interval mapping revealed five QTLs with a significant association with resistance to C. lentis race 0, six QTLs for C. lentis race 1 resistance, and three QTLs for S. botryosum resistance. Taken together, the data obtained in the study reveal that the expression of resistance to fungal diseases in L. ervoides is a result of rearrangement of resistant alleles contributed by both parental accessions.

Analysis of genetic structure and interrelationships in lentil species using morphological and SSR markers

Genetic structure and relationships of 130 lentil accessions belonging to six taxa were analysed. For this purpose, seven morphological traits and 31 polymorphic simple sequence repeat (SSR) primers were used for this purpose. Morphological traits grouped lentil accessions into five main clusters. SSR primers collectively amplified 139 polymorphic alleles in a range of 2-10 with an average of 4.48 alleles. The size of amplified alleles varied from 50 to 650 bp. Polymorphism information content (PIC) ranged from 0.02 to 0.85 with an average of 0.46. Neighbour-joining tree grouped accessions broadly according to their taxonomic ranks, except L. culinaris ssp. odemensis. Analysis of molecular variance (AMOVA) revealed that a major portion (82.0%) of genetic variance resided within species, while only 18% resided among species. Bayesian model-based STRUCTURE analysis assigned all accessions into five clusters and showed some admixture within individuals. Cluster analysis showed that cultivated Lens accessions of Ethiopian origin clustered separately, from other cultivated accessions indicating its distinct lineage. Among the analysed lentil species, L. culinaris ssp. odemensis seemed to have conserved genetic background and needs revision of its taxonomic status. Results of present study provide important information on genetic diversity and relationships among different wild and cultivated taxa of lentil. Thus, these results can be useful in designing breeding strategies for future improvement and taxonomic implications in lentil.

Genetic Diversity, Population Structure, and Linkage Disequilibrium of a Core Collection of Ziziphus jujuba Assessed with Genome-wide SNPs Developed by Genotyping-by-sequencing and SSR Markers

Chinese jujube (Ziziphus jujuba Mill) is an economically important fruit species native to China with high nutritious and medicinal value. Genotyping-by-sequencing was used to detect and genotype single nucleotide polymorphisms (SNPs) in a core collection of 150 Chinese jujube accessions and further to characterize their genetic diversity, population structure, and linkage disequilibrium (LD). A total of 4,680 high-quality SNPs were identified, of which 38 sets of tri-allelic SNPs were detected. The average polymorphism information content (PIC) values based on bi-allelic SNPs and tri-allelic SNPs were 0.27 and 0.38, respectively. STRUCTURE and principal coordinate analyses based on SNPs revealed that the 150 accessions could be clustered into two groups. However, neighbor-joining trees indicated the accessions should be grouped into three major clusters. Our data confirm that the resolving power for genetic diversity was similar for the SSRs and SNPs. In contrast, regarding population structure, the resolving power was higher for SSRs than for SNPs. The LD pattern in Chinese jujube was investigated for the first time. We observed a relatively rapid LD decay with a short range (∼10 kb) for all pseudo-chromosomes and for individual pseudo-chromosomes. Our findings provide important information for future genome-wide association analyses and marker-assisted selective breeding of Chinese jujube.

Low-depth genotyping-by-sequencing (GBS) in a bovine population: strategies to maximize the selection of high quality genotypes and the accuracy of imputation

Background

Genotyping-by-sequencing (GBS) has emerged as a powerful and cost-effective approach for discovering and genotyping single-nucleotide polymorphisms. The GBS technique was largely used in crop species where its low sequence coverage is not a drawback for calling genotypes because inbred lines are almost homozygous. In contrast, only a few studies used the GBS technique in animal populations (with sizeable heterozygosity rates) and many of those that have been published did not consider the quality of the genotypes produced by the bioinformatic pipelines. To improve the sequence coverage of the fragments, an alternative GBS preparation protocol that includes selective primers during the PCR amplification step has been recently proposed. In this study, we compared this modified protocol with the conventional two-enzyme GBS protocol. We also described various procedures to maximize the selection of high quality genotypes and to increase the accuracy of imputation.

Results

The in silico digestions of the bovine genome showed that the combination of PstI and MspI is more suitable for sequencing bovine GBS libraries than the use of single digestions with PstI or ApeKI. The sequencing output of the GBS libraries generated a total of 123,666 variants with the selective-primer approach and 272,103 variants with the conventional approach. Validating our data with genotypes obtained from mass spectrometry and Illumina’s bovine SNP50 array, we found that the genotypes produced by the conventional GBS method were concordant with those produced by these alternative genotyping methods, whereas the selective-primer method failed to call heterozygotes with confidence. Our results indicate that high accuracy in genotype calling (>97%) can be obtained using low read-depth thresholds (3 to 5 reads) provided that markers are simultaneously filtered for genotype quality scores. We also show that factors such as the minimum call rate and the minor allele frequency positively influence the accuracy of imputation of missing GBS data. The highest accuracies (around 85%) of imputed GBS markers were obtained with the FIMPUTE program when GBS and SNP50 array genotypes were combined (80,190 to 100,297 markers) before imputation.

Conclusions

We discovered that the conventional two-enzyme GBS protocol could produce a large number of high-quality genotypes provided that appropriate filtration criteria were used. In contrast, the selective-primer approach resulted in a substantial proportion of miscalled genotypes and should be avoided for livestock genotyping studies. Overall, our study demonstrates that carefully adjusting the different filtering parameters applied to the GBS data is critical to maximize the selection of high quality genotypes and to increase the accuracy of imputation of missing data. The strategies and results presented here provide a framework to maximize the output of the GBS technique in animal populations and qualified the PstI/MspI GBS assay as a low-cost high-density genotyping platform. The conclusions reported here regarding read-depth and genotype quality filtering could benefit many GBS applications, notably genome-wide association studies, where there is a need to increase the density of markers genotyped across the target population while preserving the quality of genotypes.

Flowering and Growth Responses of Cultivated Lentil and Wild Lens Germplasm toward the Differences in Red to Far-Red Ratio and Photosynthetically Active Radiation

Understanding environmental responses of pulse crop species and their wild relatives will play an important role in developing genetic strategies for crop improvement in response to changes in climate. This study examined how cultivated lentil and wild Lens germplasm responded to different light environments, specifically differences in red/far-red ratio (R/FR) and photosynthetically active radiation (PAR). Three genotypes of each the seven Lens species were grown in environmentally controlled growth chambers equipped to provide light treatments consisting of different R/FR ratios and PAR values. Our results showed that overall, days to flower of Lens genotypes were mainly influenced by the R/FR induced light quality change but not by the PAR related light intensity change. The cultivated lentil (L. culinaris) showed consistent, accelerated flowering in response to the low R/FR light environment together with three wild lentil genotypes (L. orientalis IG 72611, L. tomentosus IG 72830, and L. ervoides IG 72815) while most wild lentil genotypes had reduced responses and flowering time was not significantly affected. The longest shoot length, longest internode length, and largest leaflet area were observed under the low R/FR low PAR environment for both cultivated and wild lentils. The distinctly different responses between flowering time and elongation under low R/FR conditions among wild Lens genotypes suggests discrete pathways controlling flowering and elongation, which are both components of shade avoidance responses. The yield and above-ground biomass of Lens genotypes were the highest under high R/FR high PAR conditions, intermediate under low R/FR low PAR conditions, and lowest under high R/FR low PAR light conditions. Three L. lamottei genotypes (IG 110809, IG 110810, and IG 110813) and one L. ervoides genotype (IG 72646) were less sensitive in their time to flower responses while maintaining similar yield, biomass, and harvest index across all three light environments; these are indications of better adaptability toward changes in light environment.

Genotyping-by-sequencing provides the first well-resolved phylogeny for coffee (Coffea) and insights into the evolution of caffeine content in its species: GBS coffee phylogeny and the evolution of caffeine content

A comprehensive and meaningful phylogenetic hypothesis for the commercially important coffee genus (Coffea) has long been a key objective for coffee researchers. For molecular studies, progress has been limited by low levels of sequence divergence, leading to insufficient topological resolution and statistical support in phylogenetic trees, particularly for the major lineages and for the numerous species occurring in Madagascar. We report here the first almost fully resolved, broadly sampled phylogenetic hypothesis for coffee, the result of combining genotyping-by-sequencing (GBS) technology with a newly developed, lab-based workflow to integrate short read next-generation sequencing for low numbers of additional samples. Biogeographic patterns indicate either Africa or Asia (or possibly the Arabian Peninsula) as the most likely ancestral locality for the origin of the coffee genus, with independent radiations across Africa, Asia, and the Western Indian Ocean Islands (including Madagascar and Mauritius). The evolution of caffeine, an important trait for commerce and society, was evaluated in light of our phylogeny. High and consistent caffeine content is found only in species from the equatorial, fully humid environments of West and Central Africa, possibly as an adaptive response to increased levels of pest predation. Moderate caffeine production, however, evolved at least one additional time recently (between 2 and 4Mya) in a Madagascan lineage, which suggests that either the biosynthetic pathway was already in place during the early evolutionary history of coffee, or that caffeine synthesis within the genus is subject to convergent evolution, as is also the case for caffeine synthesis in coffee versus tea and chocolate.

A Novel Lens orientalis Resistance Source to the Recently Evolved Highly Aggressive Australian Ascochyta lentis Isolates

Substantial yield losses and poor seed quality are frequently associated with Ascochyta blight infection of lentil caused by Ascochyta lentis. Recently reported changes in aggressiveness of A. lentis have led to decreased resistance within cultivars, such as Northfield and Nipper in Australia. Furthermore, the narrow genetic base of the current breeding program remains a risk for further selective pathogen evolution to overcome other currently used resistances. Therefore, incorporation of potentially novel and diverse resistance genes into the advanced lines will aid to improve cultivar stability. To identify these, 30 genotypes sourced from five wild species (Lens orientalis, L. odomensis, L. ervoides, L. nigricans and L. lamottei), including eight previously reported resistance sources, were screened for disease reaction to two recently isolated and highly aggressive isolates. Subsequently, two L. orientalis accessions were found highly resistant and a further six L. nigricans, one L. odomensis, one L. ervoides, one L. lamottei, and one L. orientalis accessions were moderately resistant. Several of these were more resistant than the currently deployed resistance source, ILL 7537. Furthermore, L. orientalis accession ILWL 180 was consistently resistant against other highly aggressive isolates recovered from diverse geographical lentil growing regions and host genotypes, suggesting stability and potential for future use of this accession in the Australian lentil breeding program.

Evaluation of Wild Lentil Species as Genetic Resources to Improve Drought Tolerance in Cultivated Lentil

Increasingly unpredictable annual rainfall amounts and distribution patterns have far reaching implications for pulse crop biology. Seedling and whole plant survival will be affected given that water is a key factor in plant photosynthesis and also influences the evolving disease spectrum that affects crops. The wild relatives of cultivated lentil are native to drought prone areas, making them good candidates for the evaluation of drought tolerance traits. We evaluated root and shoot traits of genotypes of cultivated lentil and five wild species grown under two water deficit regimes as well as fully watered conditions over a 13 week period indoors. Plants were grown in sectioned polyvinyl chloride (PVC) tubes containing field soil from the A, B, and C horizons. We found that root distribution into different soil horizons varied among wild lentil genotypes. Secondly, wild lentil genotypes employed diverse strategies such as delayed flowering, reduced transpiration rates, reduced plant height, and deep root systems to either escape, evade or tolerate drought conditions. In some cases, more than one drought strategy was observed within the same genotype. Sequence based classification of wild and cultivated genotypes did not explain patterns of drought response. The environmental conditions at their centers of origin may explain the patterns of drought strategies observed in wild lentils. The production of numerous small seeds by wild lentil genotypes may have implications for yield improvement in lentil breeding programs.

Genotyping-by-sequencing provides the discriminating power to investigate the subspecies of Daucus carota (Apiaceae)

BACKGROUND

The majority of the subspecies of Daucus carota have not yet been discriminated clearly by various molecular or morphological methods and hence their phylogeny and classification remains unresolved. Recent studies using 94 nuclear orthologs and morphological characters, and studies employing other molecular approaches were unable to distinguish clearly many of the subspecies. Fertile intercrosses among traditionally recognized subspecies are well documented. We here explore the utility of single nucleotide polymorphisms (SNPs) generated by genotyping-by-sequencing (GBS) to serve as an effective molecular method to discriminate the subspecies of the D. carota complex.

RESULTS

We used GBS to obtain SNPs covering all nine Daucus carota chromosomes from 162 accessions of Daucus and two related genera. To study Daucus phylogeny, we scored a total of 10,814 or 38,920 SNPs with a maximum of 10 or 30 % missing data, respectively. To investigate the subspecies of D. carota, we employed two data sets including 150 accessions: (i) rate of missing data 10 % with a total of 18,565 SNPs, and (ii) rate of missing data 30 %, totaling 43,713 SNPs. Consistent with prior results, the topology of both data sets separated species with 2n = 18 chromosome from all other species. Our results place all cultivated carrots (D. carota subsp. sativus) in a single clade. The wild members of D. carota from central Asia were on a clade with eastern members of subsp. sativus. The other subspecies of D. carota were in four clades associated with geographic groups: (1) the Balkan Peninsula and the Middle East, (2) North America and Europe, (3) North Africa exclusive of Morocco, and (4) the Iberian Peninsula and Morocco. Daucus carota subsp. maximus was discriminated, but neither it, nor subsp. gummifer (defined in a broad sense) are monophyletic.

CONCLUSIONS

Our study suggests that (1) the morphotypes identified as D. carota subspecies gummifer (as currently broadly circumscribed), all confined to areas near the Atlantic Ocean and the western Mediterranean Sea, have separate origins from sympatric members of other subspecies of D. carota, (2) D. carota subsp. maximus, on two clades with some accessions of subsp. carota, can be distinguished from each other but only with poor morphological support, (3) D. carota subsp. capillifolius, well distinguished morphologically, is an apospecies relative to North African populations of D. carota subsp. carota, (4) the eastern cultivated carrots have origins closer to wild carrots from central Asia than to western cultivated carrots, and (5) large SNP data sets are suitable for species-level phylogenetic studies in Daucus.

Genetic diversity analysis of Gossypium arboreum germplasm accessions using genotyping-by-sequencing

The diploid cotton species Gossypium arboreum possesses many favorable agronomic traits such as drought tolerance and disease resistance, which can be utilized in the development of improved upland cotton cultivars. The USDA National Plant Germplasm System maintains more than 1600 G. arboreum accessions. Little information is available on the genetic diversity of the collection thereby limiting the utilization of this cotton species. The genetic diversity and population structure of the G. arboreum germplasm collection were assessed by genotyping-by-sequencing of 375 accessions. Using genome-wide single nucleotide polymorphism sequence data, two major clusters were inferred with 302 accessions in Cluster 1, 64 accessions in Cluster 2, and nine accessions unassigned due to their nearly equal membership to each cluster. These two clusters were further evaluated independently resulting in the identification of two sub-clusters for the 302 Cluster 1 accessions and three sub-clusters for the 64 Cluster 2 accessions. Low to moderate genetic diversity between clusters and sub-clusters were observed indicating a narrow genetic base. Cluster 2 accessions were more genetically diverse and the majority of the accessions in this cluster were landraces. In contrast, Cluster 1 is composed of varieties or breeding lines more recently added to the collection. The majority of the accessions had kinship values ranging from 0.6 to 0.8. Eight pairs of accessions were identified as potential redundancies due to their high kinship relatedness. The genetic diversity and genotype data from this study are essential to enhance germplasm utilization to identify genetically diverse accessions for the detection of quantitative trait loci associated with important traits that would benefit upland cotton improvement.

GBS-SNP-CROP: a reference-optional pipeline for SNP discovery and plant germplasm characterization using variable length, paired-end genotyping-by-sequencing data

BACKGROUND

With its simple library preparation and robust approach to genome reduction, genotyping-by-sequencing (GBS) is a flexible and cost-effective strategy for SNP discovery and genotyping, provided an appropriate reference genome is available. For resource-limited curation, research, and breeding programs of underutilized plant genetic resources, however, even low-depth references may not be within reach, despite declining sequencing costs. Such programs would find value in an open-source bioinformatics pipeline that can maximize GBS data usage and perform high-density SNP genotyping in the absence of a reference.

RESULTS

The GBS SNP-Calling Reference Optional Pipeline (GBS-SNP-CROP) developed and presented here adopts a clustering strategy to build a population-tailored "Mock Reference" from the same GBS data used for downstream SNP calling and genotyping. Designed for libraries of paired-end (PE) reads, GBS-SNP-CROP maximizes data usage by eliminating unnecessary data culling due to imposed read-length uniformity requirements. Using 150 bp PE reads from a GBS library of 48 accessions of tetraploid kiwiberry (Actinidia arguta), GBS-SNP-CROP yielded on average three times as many SNPs as TASSEL-GBS analyses (32 and 64 bp tag lengths) and over 18 times as many as TASSEL-UNEAK, with fewer genotyping errors in all cases, as evidenced by comparing the genotypic characterizations of biological replicates. Using the published reference genome of a related diploid species (A. chinensis), the reference-based version of GBS-SNP-CROP behaved similarly to TASSEL-GBS in terms of the number of SNPs called but had an improved read depth distribution and fewer genotyping errors. Our results also indicate that the sets of SNPs detected by the different pipelines above are largely orthogonal to one another; thus GBS-SNP-CROP may be used to augment the results of alternative analyses, whether or not a reference is available.

CONCLUSIONS

By achieving high-density SNP genotyping in populations for which no reference genome is available, GBS-SNP-CROP is worth consideration by curators, researchers, and breeders of under-researched plant genetic resources. In cases where a reference is available, especially if from a related species or when the target population is particularly diverse, GBS-SNP-CROP may complement other reference-based pipelines by extracting more information per sequencing dollar spent. The current version of GBS-SNP-CROP is available at https://github.com/halelab/GBS-SNP-CROP.git.

Identification of Quantitative Trait Loci Controlling Root and Shoot Traits Associated with Drought Tolerance in a Lentil (Lens culinaris Medik.) Recombinant Inbred Line Population

Drought is one of the major abiotic stresses limiting lentil productivity in rainfed production systems. Specific rooting patterns can be associated with drought avoidance mechanisms that can be used in lentil breeding programs. In all, 252 co-dominant and dominant markers were used for Quantitative Trait Loci (QTL) analysis on 132 lentil recombinant inbred lines based on greenhouse experiments for root and shoot traits during two seasons under progressive drought-stressed conditions. Eighteen QTLs controlling a total of 14 root and shoot traits were identified. A QTL-hotspot genomic region related to a number of root and shoot characteristics associated with drought tolerance such as dry root biomass, root surface area, lateral root number, dry shoot biomass and shoot length was identified. Interestingly, a QTL (QRSratioIX-2.30) related to root-shoot ratio, an important trait for drought avoidance, explaining the highest phenotypic variance of 27.6 and 28.9% for the two consecutive seasons, respectively, was detected. This QTL was closed to the co-dominant SNP marker TP6337 and also flanked by the two SNP TP518 and TP1280. An important QTL (QLRNIII-98.64) related to lateral root number was found close to TP3371 and flanked by TP5093 and TP6072 SNP markers. Also, a QTL (QSRLIV-61.63) associated with specific root length was identified close to TP1873 and flanked by F7XEM6b SRAP marker and TP1035 SNP marker. These two QTLs were detected in both seasons. Our results could be used for marker-assisted selection in lentil breeding programs targeting root and shoot characteristics conferring drought avoidance as an efficient alternative to slow and labor-intensive conventional breeding methods.

Genetic Diversity in Lens Species Revealed by EST and Genomic Simple Sequence Repeat Analysis

Low productivity of pilosae type lentils grown in South Asia is attributed to narrow genetic base of the released cultivars which results in susceptibility to biotic and abiotic stresses. For enhancement of productivity and production, broadening of genetic base is essentially required. The genetic base of released cultivars can be broadened by using diverse types including bold seeded and early maturing lentils from Mediterranean region and related wild species. Genetic diversity in eighty six accessions of three species of genus Lens was assessed based on twelve genomic and thirty one EST-SSR markers. The evaluated set of genotypes included diverse lentil varieties and advanced breeding lines from Indian programme, two early maturing ICARDA lines and five related wild subspecies/species endemic to the Mediterranean region. Genomic SSRs exhibited higher polymorphism in comparison to EST SSRs. GLLC 598 produced 5 alleles with highest gene diversity value of 0.80. Among the studied subspecies/species 43 SSRs detected maximum number of alleles in L. orientalis. Based on Nei's genetic distance cultivated lentil L. culinaris subsp. culinaris was found to be close to its wild progenitor L. culinaris subsp. orientalis. The Prichard's structure of 86 genotypes distinguished different subspecies/species. Higher variability was recorded among individuals within population than among populations.