Next Article in Journal
Genetic Characterization of Spring Wheat Germplasm for Macro-, Microelements and Trace Metals
Previous Article in Journal
Climatic Factors Determine the Distribution Patterns of Leaf Nutrient Traits at Large Scales
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 11
Issue 16
10.3390/plants11162172
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Construction of A GBS-Based High-Density Genetic Map and Flower Color-Related Loci Mapping in Grasspea (Lathyrus sativus L.)

by
Xiaopeng Hao
1,2,3,4,
Tao Yang
5,
Yan Wang
2,3,4,
Rong Liu
5,
Xue Dong
2,3,4,
Jiandong Zhao
2,3,4,
Jucai Han
1,
Xuxiao Zong
5,
Jianwu Chang
2,3,4 and
Huiping Liu
1,*
1
College of Plant Protection, Shanxi Agricultural University, Taigu 030801, China
2
Center for Agricultural Genetic Resources Research, Shanxi Agricultural University, Taiyuan 030031, China
3
Key Laboratory of Crop Gene Resources and Germplasm Enhancement on Loess Plateau, Ministry of Agriculture, Taiyuan 030031, China
4
Shanxi Key Laboratory of Genetic Resources and Genetic Improvement of Minor Crops, Taiyuan 030031, China
5
Institute of Crop Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100089, China
*
Author to whom correspondence should be addressed.
Plants 2022, 11(16), 2172; https://doi.org/10.3390/plants11162172
Submission received: 8 July 2022 / Revised: 17 August 2022 / Accepted: 18 August 2022 / Published: 21 August 2022
(This article belongs to the Section Crop Physiology and Crop Production)
Browse Figures
Versions Notes

Abstract

:
Grasspea (Lathyrus sativus L.), a legume crop with excellent resistance to a broad array of environmental stressors, has, to this point, been poorly genetically characterized. High-density genetic linkage maps are critical for draft genome assembly, quantitative trait loci (QTLs) analysis, and gene mining. The lack of a high-density genetic linkage map has limited both genomic studies and selective breeding in grasspea. Here, we developed a high-density genetic linkage map of grasspea using genotyping-by-sequencing (GBS) to sequence 154 grasspea plants, comprising 2 parents and 152 F2 progeny. In all, 307.74 Gb of data was produced, including 2,108,910,938 paired-end reads, as well as 3536 SNPs mapped to seven linkage groups (LG1–LG7). With an average length of 996.52 cM per LG, the overall genetic distance was 6975.68 cM. Both the χ2 test and QTL analysis, based on the Kruskal–Wallis (KW) test and interval mapping (IM) analysis, revealed the monogenic inheritance of flower color in grasspea, with the responsible QTL located between 308.437 cM and 311.346 cM in LG4. The results can aid grasspea genome assembly and accelerate the selective breeding of new grasspea germplasm resources.

1. Introduction

The grasspea (Lathyrus sativus L.) is an annual, cool-season legume with an ~8.2 Gb diploid genome (2n = 14) [1,2,3,4]. Grasspea originated in the Eastern Mediterranean and was domesticated in the Balkans about 8000 BC [5]. Today, grasspea is considered an orphan crop due to its limited cultivation area. However, it is a hardy and climate-smart crop that can be cultivated in environments which are unsuitable for the cultivation of other crops [6,7,8,9]. Grasspea is highly tolerant of an array of abiotic stressors, including drought, flooding, nutrient scarcity, and high soil salinity and alkalinity, and shows good resistance to several plant pathogens, including rust and powdery mildew [10,11,12,13,14,15]. Like other legumes, grasspea has the ability to fix nitrogen, thereby potentially reducing the need for exogenous fertilizer application. Because of these characteristics, grasspea has the potential to meet the food demands of low-income countries and others with marginal land [16]. Unfortunately, grasspea contains the neurotoxin N-oxalyl-L-α,β-diaminopropionic acid (β-ODAP), which causes neurolathyrism when consumed in large quantities [17,18,19]. Despite this limitation, as of 2014, the grasspea planting area reached 1,500,000 hectares, with an annual production of 1,200,000 tons [20]. Even with moderate acreage and productivity, grasspea is gaining increasing attention due to its resilience under unfavorable environmental conditions [21].
Genetic linkage mapping is vital for genome assembly and trait loci mapping. Since the first genetic linkage map was published in 1913, a plethora of linkage maps based on phenotypic, protein, and DNA markers have been produced [22]. Due to the limitations of phenotypic and protein markers, DNA markers are increasingly utilized for the generation of genetic linkage maps [23]. In 1980, restriction fragment length polymorphism (RFLP) markers were used to construct the first DNA-based human genetic map [24]. More recently, the advancement of next-generation sequencing (NGS) technology has accelerated the development and application of molecular markers for use in constructing genetic maps, including genotyping-by-sequencing (GBS), restriction-site-associated DNA sequencing (RAD-seq), specific-locus amplified fragment sequencing (SLAF-seq), and re-sequencing [25,26,27,28]. Of these, GBS is considered one of the most efficient and cost-effective technologies for creating linkage maps. In comparison to the other genotypeing methods, GBS is more precise, efficient, and cost-effective [29]. However, despite its potential to address food insecurity, no major genomic evaluation has been conducted in grasspea. Because of its poor yield, large genome, high outcrossing rate, and neurotoxic disadvantage, grasspea has received little attention from researchers and, as a result, grasspea molecular biology and genomic research has fallen significantly behind that of more prominent legumes [16,18,30].
In recent years, smart breeding methods based on marker-assisted selection (MAS) have been used for crop breeding, and a genetic linkage map is a prerequisite for this technique. It is now possible to create a high-density genetic map for the grasspea [31]. Chowdhury (1999) developed the first grasspea genetic linkage map based on 69 markers produced from an F2 population of 102 parents and progeny, including 65 random amplified polymorphic DNA (RAPD) markers, 3 isozyme markers, and 1 phenotypic marker, mapped to 14 linkage groups (LGs) [32]. Next, Skiba (2004) developed a genetic map using backcross individuals produced from parents either resistant or susceptible to Ascochyta Blight induced by Mycosphaerella pinodes [33]. Using 47 RAPD markers, 7 sequence-tagged microsatellite sites (STMS) markers, and 13 sequence-tagged sites (STS)/cleaved amplified polymorphic sequence (CAPS) markers, the map covered a genetic distance of 803.1 cM, with an average of 15.8 cM per marker, mapped to nine LGs [33]. Most recently, Diversity Array Technology sequencing (DArTseq) has been used to create a true high-density genetic map using 105 recombinant inbred lines (RILs) and 2149 markers, including 1472 SilicoDArT markers and 677 single-nucleotide polymorphisms (SNPs) [34]. Based on RIL screening, LsMLO1, a key gene responsible for powdery mildew susceptibility in both L. sativus and L. cicera, was identified, described, and mapped to the upper portion of LG1, at position 18.246 cM, in the L. sativus genetic map [34]. As the closest wild relative of L. sativus, it is perhaps equally important to examine the genome of L. cicera. To that end, Santos (2018, 2020) created two genetic maps for L. cicera, including 307 and 1468 markers, respectively, mapped to nine LGs [35,36].
Phenotypic markers, such as flower, seed, and fruit color and leaf form, are critical for differentiating between various germplasm resources [37,38,39]. In grasspea, flower color shows the greatest phenotypic diversity, with accessions exhibiting blue, white, or pink flowers. This diversity may be related to functional differences in pollinator attraction and success [40]. In this work, artificial hybridization between two grasspea accessions, IF1347 and K714, was undertaken in order to develop new, more diverse germplasm resources. An F2 population of 152 individuals was obtained and used to produce a high-density genetic map, resulting in the discovery of several flower-color-associated loci. It is our hope that these findings can aid in the assembly of a draft grasspea genome and enable future selective breeding efforts for this promising legume crop.

2. Results

2.1. GBS Data Analysis and Assessment

A total of 154 samples, including 2 parents and 152 F2 progeny, were sequenced, yielding 307.77 Gb of raw data. After quality control, 307.74 Gb of high-quality reads, including 2,108,910,938 paired-end reads, were available for SNP calling. The raw data contained an average of 99.97% clean reads, with a Q20 of 95.24%, a Q30 of 87.66%, and GC of 38.6%. The male (P1, 2,272,102,272 bp) and female (P2, 2,300,245,632 bp) parents were sequenced and the raw data contained an average of 100% clean reads (Table S1). For the 152 F2 progeny, clean reads varied from 1,301,804,352 to 2,670,197,184 bp, with an average of 1,994,496,457 bp.

2.2. SNP Calling and Genotyping

Clean reads were aligned using the ‘Rbp’ reference genome, based on the European grasspea accession ‘LS007’, supplied by JIC [7] (Table S2). A total of 15,617,673 clean reads were mapped for the male parent and 15,808,872 clean reads were mapped for the female parent. An average of 13,667,660 clean reads were mapped to the draft genome for the 152 F2 progeny. In all, approximately 98.68% of the data could be mapped to the ‘Rbp’ genome. Only the clean reads mapped to the draft genome were used for SNP calling and genotyping.
The clean reads of both parents and F2 progeny were aligned to the ‘Rbp’ draft reference genome using the Burrows-Wheeler Aligner (BWA) software with the setting “mem -t 4 -k 32 -M -R” [41]. SAMtools was used to call SNPs in all samples once the alignment files (SAM) were transformed into BAM files [42]. For the male parent (P1), 551,846 SNPs were called, with 39.06% (215,557) homozygous and 60.94% heterozygous SNPs (Table S3). For the female parent (P2), 522,270 SNPs were called, with 61.19% (319,583) homozygous and 38.81% heterozygous SNPs (Table S3). An average of 244,236 homozygous SNPs were called across F2 progeny, ranging from 38.42% to 55.95% in each F2 individual. Polymorphic SNPs in P1 and P2 were categorized into seven segregation patterns, including ‘ef × eg’, ‘hk × hk’, ‘lm × ll’, ‘nn × np’, ‘aa × bb’, ‘ab × cc’, and ‘cc × ab‘. However, only markers with ‘aa × bb’-type segregation, encompassing 77,346 SNPs between the parents, were chosen for genetic mapping of the F2 population (Figure 1).

2.3. Genetic Linkage Map Constructing and Verification

For the F2 progeny, abnormal markers with <70% accuracy were eliminated from further analyses. To identify segregation distortion, chi-square (χ2) tests were performed on all SNPs. Markers exhibiting a statistically significant χ2 result (p < 0.001) or aberrant bases were eliminated. Finally, 4428 markers with ‘aa × bb’-type segregation were kept. A high-resolution genetic linkage map was created using JoinMap4.1, based on the Kosambi genetic distance calculation and regression ordering method. A total of 3536 SNP markers were mapped to seven LGs, with a total genetic length of 6975.68 cM and an average LG length of 996.52 cM, ranging from 504.11 cM (LG1) to 1367.47 cM (LG4). The average length between two neighboring markers was 2.21 cM. The number of anchoring markers per LG ranged from 290 (LG4) to 1240 (LG1), with an average of 505 markers per LG (Supplementary file S1). LG1 had the highest marker density (1.1 cM per interval), while LG7 had the lowest marker density (2.98 cM per interval). The majority of the gaps between markers were less than 5 cM in length, with the largest gap consisting of 31.47 cM between mk1930 and mk4366 in LG4. In total, 3205 SNPs (<5 cM), accounting for 90.64% of 3536, were mapped to LG1-LG7 (Table 1 and Table 2, and Figure S1).
Both haplotype map and heat map analysis are crucial for evaluating the quality of high-resolution genetic maps. Haplotype mapping can identify markers with genotyping errors in the segregating population and heat mapping can visualize linkage connections per marker in each LG. These two methods were employed to discover any possible ordering errors or linkage connections between markers within a single LG across the seven L. sativus LGs. Haplotype maps were constructed using genotypes of SNP markers from both parents and F2 progeny. The recombination connection between markers in the haplotype map was utilized to detect marker ordering mistakes, which were indicated by double crossings (Figure S2). Between both parents and F2 progeny, the majority of areas across the seven haplotype maps exhibited similar color, indicating clear recombination blocks and ensuring the creation of a high-quality genetic map. Based on the heat maps of the seven LGs, we observed a significant connection between two nearby markers in seven LGs. Specifically, as the genetic distance between two markers increased, the link between them decreased, revealing the correct sequence of markers in the seven LGs (Figure S3).

2.4. QTL Analysis of Flower Color-Related Traits

Individuals with blue flowers (male) outnumbered those with white flowers (female) by a ratio of 115:37, in agreement with Mendel’s idealized 3:1 segregation ratio, demonstrating monogenic inheritance (χ2 = 0.0175, p > 0.05). Kruskal–Wallis (KW) analysis and interval mapping (IM) were employed to identify the significant loci affecting flower color in the F2 population. According to the KW analysis, the QTL regulating flower color was found in LG4 at location 308.437 cM (mk3948) and had a KW value of 117.99 (Figure 2). IM similarly identified the genes in LG4, but with a location of 311.346 cM (mk2691), a LOD value of 103.98, and a percentage of the variance explained (PVE) of 96.5% (Figure 3).

3. Discussion

Grasspea, an orphan food crop, shows strong resilience in the face of difficult environmental conditions, including drought, flooding, high salinity, and cold temperatures. Although grasspea has great potential to be developed into a robust crop for challenging environments, genetic research has so far been lacking [1,7,16,30]. Genetic linkage mapping is an essential tool for assembling genomes and gene mining. Previously, RAPD, simple sequence repeat (SSR), SNP, and DArT markers have been used to create grasspea genetic maps using PCR, transcriptomic, and GBS technologies [32,33,34,35,36]. Based on NGS techniques, GBS is an efficient and powerful method for identifying enrichment variation [43]. Through DNA extraction, restriction enzyme digestion, fragment selection, library assembly, and bioinformatics analysis, many variants, including SNPs and insertions/deletions (InDels), may be identified [44,45]. After the most highly significant genes are revealed, GBS may be used for high-density genetic mapping, genetic diversity and population structure analysis, QTL mapping, and genome-wide association studies (GWAS) [46,47,48]. For cool-season legumes, particularly peas, lentils, chickpeas, fava beans, and grasspeas, with large genomes (~1.0–13.0 Gb), GBS has significant benefits over de novo sequencing and re-sequencing in terms of speed, efficiency, and labor- and cost-effectiveness. In this study, GBS captured 307.74 Gb of clean data comprising 2,108,910,938 high-quality paired-end reads, with a Q30 of 87.66%. Clean reads were mapped to the ‘Rbp’ draft grasspea genome, ultimately resulting in 77,346 SNPs with an ‘aa × bb’ segregation pattern (other patterns were excluded from further analysis). Of these, 3536 SNPs were mapped to the seven LGs, resulting in a total genomic length of 6975.68 cM with an average genomic marker interval of 2.21 cM. This genetic map can be used to assist in locating continuous sequences (contigs) on the physical map for L. sativus draft genome assembly.
High-resolution genetic maps are critical for modern selective breeding. The use of such maps has led to the identification of a large number of agronomically relevant QTLs in a wide array of food crops, including maize, foxtail millet, and jujube [49,50,51]. Recently, high-resolution genetic maps have been created for both L. sativus and the closely related L. cicera [34,36]. In the most recent of these maps, 1468 markers were identified across nine LGs (seven main LGs and two minor LGs), with a total length of 712.35 cM and an average interval distance of 0.65 cM. These markers included 730 silicoDArT markers, 623 SNPs, 110 expressed sequence tag-derived SSR (ESSR) markers, and 5 intron targeted amplified polymorphism (ITAP) markers. Three Erysiphe pisi–responsive and one Erysiphe trifolii-responsive QTL were found to explain 13.0%, 10.6%, 8.2%, and 16% of the phenotypic variation, respectively [36]. More recently, utilizing on DArT-Seq, 2149 markers were mapped to 10 LGs (7 main LGs and 3 minor LGs), with a total genetic distance of 674.4 cM and an average interval distance of 0.41 cM [34]. The SNP related to LsMLO1, controlling powdery mildew susceptibility, was found on LG1 at position 18.246 cM, and a transversion (G > T) is likely responsible for phenotypic variation in powdery mildew resistance [33]. In this work, the constructed genetic map covered a total genetic distance of 6975.68 cM with an average interval distance of 2.21 cM, which is three–five times larger than previous reported for L. sativus or L. cicera [27,28,31,34,36,50]. Due to a high outcrossing rate, our parental accessions exhibited high heterogeneity, only containing between 38.81% and 60.94% heterozygous SNPs, resulting in longer marker interval distances. The purification of parental germplasm may be necessary to develop a high-quality high-resolution genetic linkage map with increased SNP density and decreased marker interval distances. Additionally, we found that more than 1/3 of all SNPs (1240) were mapped to LG1, which contained 2.6–4.3 times more SNPs than the other LGs. As mentioned above, previous high-resolution genetic maps of L. sativus and L. cicera contained two or three minor LGs alongside to the seven primary LGs. This may be due to chromosomal rearrangement, as chromosomal instability has been reported in grasspea since 1982 [52].
Flower color is an important reproductive characteristic in diverse plant species, with both monogenic or polygenic inheritance [53,54]. We found that, in grasspea, flower color is controlled by a single gene through monogenic inheritance with a segregation ratio of 3:1 (χ2 = 0.0175, p > 0.05). This is similar to both cowpea and soybean, in that QTL analysis of these species revealed the presence of a single QTL controlling flower color [55,56]. Specifically, the gene Vigun07g110700, responsible for flavonoid biosynthesis, was identified as controlling flower color, and is similar to the Arabidopsis TT8 gene, which is responsible for the biosynthesis of anthocyanins [55]. These results are in contrast to the unrelated species Salvia miltiorrhiza, in which two QTLs appear to be responsible for flower color (qfcRGB4 and qfcRGB5) [38]. These flower color QTLs were identified on LG4 with a KW value of 117.99 at position 308.437 cM (mk3948) and a LOD value of 103.98 at position 311.346 cM (mk2691, PVE = 96.5%), based on both KW and IM analysis, demonstrating the consistency of these two methods. Only flower size is more important than flower color in attracting bees and ensuring successful pollination [57,58].

4. Materials and Methods

4.1. Plant Materials

The F1 population was created through the artificial hybridization of two grasspea cultivars: the blue-flowered ICARDA accession ‘IF1347’ was used as the male parent (P1) and the white-flowered Polish accession ‘k714’ was used as the female parent (P2). The P1, P2, and F2 population from 2016 were planted at the experimental location (38.14° N, 113.04° E, altitude 1215 m), Nantou Village, Yu County, Shanxi Province, China.

4.2. DNA Sample Preparation

Young leaves from both parents and progeny were harvested four weeks after seeding in May 2018. The leaves were immediately frozen in liquid nitrogen and held at −80 °C in the laboratory. The genomic DNA (gDNA) from 154 grasspea plants (2 parents and 152 F2 progeny) was extracted using a DNAsecure Plant Kit (Tiangen Biotec Co., Ltd., Beijing, China). The quality of the gDNA was determined using 1% agarose gel electrophoresis. The purity and concentration of the gDNA were determined using a Nanophotometer spectrophotometer (Implen, Westlake, CA, USA) and a Fluorometer 2.0 (Life Technologies, Carlsbad, CA, USA).

4.3. Library Construction and GBS Sequencing

GBS sequencing was conducted according to the method of Elshire et al., 2011 [29], with minor improvements. Briefly, the gDNA from 154 grasspea plants (2 parents and 152 F2 progeny) was first digested by MseI (5’-T!TAA-3’) (New England Biolabs, ‘NEB’, Ipswich, MA, USA) at 37 °C, then subjected to restriction-ligation by T4 DNA ligase (NEB) and ATP (NEB) with a MseI Y-adapter N-containing barcode at 65 °C, followed by a second digestion with HaeIII (5’-GG!CC-3’) (NEB) at 37 °C. The digested fragments were purified using Agencourt AMPure XP (Beckman Coulter, ‘Beckman’, Brea, CA, USA) and used for PCR amplification, utilizing Phusion Master Mix (NEB) to add universal and index primers, as well as i5 and i7 index sequences, to the digested fragments. The purified PCR fragments (425–490 bp including indexes and adaptors) were screened and extracted using a 2% agarose gel with the Gel Extraction Kit (Qiagen, Valencia, CA, USA). The amplification products were purified using Agencourt AMPure XP (Beckman) and diluted before sequencing. Finally, an Illumina NovaSeq 6000 platform (Illumina, San Diego, CA, USA) was used to perform 150-bp paired-end sequencing by Tianjin Novogene Technology Co., Ltd. (Tianjin, China).

4.4. Sequencing Data Analysis and Quality Control

The barcodes in the raw data (reads) were used to classify the sequencing data from both parents and progeny. Because low-quality paired-end reads, caused primarily by base-calling duplicates and adapter contamination, can lead to unreliable raw data, the raw data (in FASTQ format) was subjected to quality control (QC) using C and Perl scripts. The QC criteria were as follows: (1) filtering read pairs with sequencing adapter; (2) filtering read pairs with one read comprising N > 10%; (3) filtering read pairs with >50% low-quality bases (Q ≤ 5).

4.5. SNP Identification and Genotyping

The BWA software was used to map the clean data from each sample to the reference draft genome (settings: mem -t 4 -k 32 -M -R) [41]. The SAMtools program was used to convert alignment data to BAM files [41]. Polymorphic markers between P1 and P2 were identified and classified into 7 segregation patterns (ef × eg, hk × hk, lm × ll, nn × np, aa × bb, ab × cc, and cc × ab) using the JoinMap 4.1 software [59]. Only markers with ‘aa × bb’-type segregation were chosen for genetic mapping of the F2 population.

4.6. Genetic Linkage Map Construction and Evaluation

All SNPs with discordant segregation ratios (p < 0.001), lacking ≥ 30% of the genotype data, or containing aberrant bases were removed using an in-house Perl script. JoinMap 4.1 software was used to calculate genetic distance using the Kosambi function, and the regression mapping approach, with a LOD value of 2–30, was used to arrange the markers into 7 LGs [60]. Finally, JoinMap 4.1 was used to identify marker order in each LG, and a Perl script SVG was used to visualize output maps. Both haplotype and heat maps were generated to analyze the genetic maps using R 4.1.0.

4.7. Flower Color Data Collection and QTL Analysis

At flowering time, the flower color of both grasspea parents and F2 progeny was examined. All progeny exhibited either blue or white, in accordance with the flower color of the male (blue) and female (white) parent. MapQTL 5.0 and CMplot by R 4.1.0 were used for the qualitative trait analysis and Manhattan figure drawing. χ2 tests, KW tests, and IM analyses were used to identify monogenic and oligogenic features, with a LOD > 3.0 indicating statistical significance.

5. Conclusions

An F2 population consisting of 152 offspring, along with the 2 parents, was sequenced by GBS and used to create a high-density genetic linkage map for grasspea. In all, a total of 3536 SNP markers were mapped to seven LGs, with a total genetic distance of 6975.68 cM and an average genetic distance of 2.21 cM. The phenotypic segregation ratio analysis (χ2 test) of F2 flower color revealed that this trait is governed by monogenic inheritance. The grasspea flower-color-related QTLs were located on LG4 by both KW and IM analysis. These results provide an important resource for the accelerated development of improved grasspea germplasm.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants11162172/s1, Table S1. Statistic of sequencing data by GBS in male (P1), female (P2), and F2 individuals of grasspea. Table S2. Statistic of clean reads, mapped reads, and sequencing depth in male (P1), female (P2), and F2 individuals of grasspea. Table S3. Statistic of homozygous SNP, heterozygous SNP, and Homozygous SNP rate in male (P1), female (P2), and F2 individuals of grasspea. Figure S1. Genetic linkage map and SNP distribution across 7 LGs. The x-axis indicates the linkage group, and the y-axis indicates genetic distance of markers. Figure S2 (1–7). Similar source distribution of each marker from 1 to 7 LGs. The x-axis indicates the F2 individuals and the y-axis indicates the mapped SNPs. The ‘green (a)’ indicates genotype ‘aa’, ‘blue green (b)’ indicates ‘bb’, ‘purple (h)’ indicates ‘ab’, and ‘orange (-)’ indicates SNP absent. Figure S3 (1–7). Heat map of linkage relationships between adjacent markers from 1 to 7 LGs. The x-axis and y-axis indicate mapped SNPs, and ‘red’ to ‘yellow’ showed that linkage relationships (genetic distance) were from far to near. Supplementary file S1. Markers name, genetic distances, positions and genotyping of 3536 SNPs in 7 LGs.

Author Contributions

Conceptualization, T.Y., X.H. and X.Z.; methodology, T.Y., X.H. and H.L., formal analysis, X.H. and Y.W.; data curation, X.H. and H.L.; writing—original draft preparation, X.H.; writing—review and editing, T.Y., Y.W., X.D., R.L., X.Z. and H.L.; supervision, J.Z., J.H. and J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Fundamental Research Program of Shanxi Province, under grant no. 20210302123417, and the China Agriculture Research System of MOF and MARA- Food Legumes, under grant no. CARS-08, and the Agricultural Science and Technology Innovation Program (ASTIP) of the Chinese Academy of Agricultural Sciences, and the National Key Research and Development Program of China, under grant no. 2021YFD1600600.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in http://db.cngb.org/cnsa/project/CNP0003158_ec6a3961/reviewlink/ (Project accession: CNP0003158), which can be accessed on 1 January 2023.

Acknowledgments

We thank Peter M. F. Emmrich (John Innes Centre, JIC, Norwich), who kindly provided the grasspea draft genome (Rbp.fa.gz) for the SNP calling, and thank the Noel Ellis’ discussion and advice.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Narayan, R.K.J. Discontinuous DNA variation in the evolution of plant species—The genus Lathyrus. Evol. 1982, 36, 877–891. [Google Scholar] [CrossRef] [PubMed]
  2. Roy, B. On the somatic chromosomes in Lathyrus. Cytologia 1936, 7, 424–430. [Google Scholar] [CrossRef] [Green Version]
  3. Talukdar, D. Cytogenetic characterization of induced autotetraploids in grass pea (Lathyrus sativus L.). Caryologia 2010, 63, 62–72. [Google Scholar]
  4. Verma, S.C.; Ohri, D. Chromosome and nucclear phenotype in the legume Lathyrus sativus L. Cytologia 1979, 44, 77–90. [Google Scholar] [CrossRef] [Green Version]
  5. Kislev, M.E. Origins of the cultivation of Lathyrus sativus and L. cicera (fabaceae). Econ. Bot. 1989, 43, 262–270. [Google Scholar] [CrossRef]
  6. Cullis, C.; Kunert, K.J. Unlocking the potential of orphan legumes. J. Exp. Bot. 2017, 68, 1895–1903. [Google Scholar] [CrossRef] [Green Version]
  7. Emmrich, P.M.F.; Sarkar, A.; Njaci, I.; Kaithakottil, G.G.; Ellis, N.; Moore, C.; Edwards, A.; Heavens, D.; Waite, D.; Cheema, J.; et al. A draft genome of grass pea (Lathyrus sativus), a resilient diploid legume. bioRxiv 2020, 1–33. [Google Scholar] [CrossRef]
  8. Ramya, K.R.; Tripathi, K.; Pandey, A.; Barpete, S.; Gore, P.G.; Raina, A.P.; Khawar, K.M.; Swain, N.; Sarker, A. Rediscovering the potential of multifaceted orphan legume grasspea—A sustainable resource with high nutritional values. Front. Nutr. 2022, 8, 826208. [Google Scholar] [CrossRef]
  9. Sarkar, A.; Emmrich, P.M.F.; Sarker, A.; Zong, X.; Martin, C.; Wang, T.L. Grass pea: Remodeling an ancient insurance crop for climate resilience. In Genomic Designing of Climate-Smart Pulse Crops; Kole, C., Ramanna, R., Eds.; Springer Nature: Cham, Switzerland, 2019; pp. 426–434. [Google Scholar]
  10. Patto, M.C.V.; Fernández-Aparicio, M.; Moral, A.; Rubiales, D. Characterization of resistance to powdery mildew (Erysiphe pisi) in a germplasm collection of Lathyrus sativus. Plant Breed. 2006, 125, 308–310. [Google Scholar] [CrossRef]
  11. Patto, M.C.V.; Rubiales, D. Resistance to rust and powdery mildew in Lathyrus crops. Czech J. Genet. Plant 2014, 50, 116–122. [Google Scholar] [CrossRef] [Green Version]
  12. Tokarz, B.; Wójtowicz, T.; Makowski, W.; Jędrzejczyk, R.J.; Tokarz, K.M. What is the difference between the response of grass pea (Lathyrus sativus L.) to salinity and drought stress?—A physiological study. Agronomy 2020, 10, 833. [Google Scholar] [CrossRef]
  13. Wiraguna, E.; Malik, A.I.; Colmer, T.D.; Erskine, W. Tolerance of four grain legume species to waterlogging, hypoxia and anoxia at germination and recovery. AoB Plants 2021, 13, plab052. [Google Scholar] [CrossRef] [PubMed]
  14. Wiraguna, E.; Malik, A.I.C.; Timothy, D.; Erskine, W. Waterlogging tolerance of grass pea (Lathyrus sativus L.) at germination related to country of origin. Exp. Agr. 2020, 56, 837–850. [Google Scholar] [CrossRef]
  15. Zhou, L.; Cheng, W.; Hou, H.; Peng, R.; Hai, N.; Bian, Z.; Jiao, C.; Wang, C. Antioxidative responses and morpho-anatomical alterations for coping with flood-induced hypoxic stress in grass pea (Lathyrus sativus L.) in comparison with pea (Pisum sativum). J. Plant Growth Regul. 2016, 35, 690–700. [Google Scholar] [CrossRef]
  16. Lambein, F.; Travella, S.; Kuo, Y.-H.; Van Montagu, M.; Heijde, M. Grass pea (Lathyrus sativus L.): Orphan crop, nutraceutical or just plain food? Planta 2019, 250, 821–838. [Google Scholar] [CrossRef] [Green Version]
  17. Das, A.; Parihar, A.K.; Barpete, S.; Kumar, S.; Gupta, S. Current perspectives on reducing the β-ODAP content and improving potential agronomic traits in grass pea (Lathyrus sativus L.). Front. Plant Sci. 2021, 12, 703275. [Google Scholar] [CrossRef]
  18. Singh, S.S.; Rao, S.L.N. Lessons from neurolathyrism: A disease of the past & the future of Lathyrus sativus (Khesari dal). Indian J. Med. Res. 2013, 138, 32–37. [Google Scholar]
  19. Yan, Z.Y.; Spencer, P.S.; Li, Z.X.; Liang, Y.M.; Wang, Y.F.; Wang, C.Y.; Li, F.M. Lathyrus sativus (grass pea) and its neurotoxin ODAP. Phytochemistry 2006, 67, 107–121. [Google Scholar] [CrossRef]
  20. Sammour, R.H. Genetic diversity in Lathyrus sativus L. germplasm. Res. Rev. BioSciences 2014, 8, 325–336. [Google Scholar]
  21. Tadele, Z. African orphan crops under abiotic stresses: Challenges and opportunities. Hindawi Sci. 2018, 2018, 1451894. [Google Scholar] [CrossRef] [Green Version]
  22. Crow, J.F. A diamond anniversary: The first chromosomal map. Genetics 1988, 118, 1–3. [Google Scholar] [CrossRef] [PubMed]
  23. Singh, B.D.; Singh, A.K. Marker-Assisted Plant Breeding: Principles and Practices; Springer: New Delhi, India, 2015; pp. 19–25. [Google Scholar]
  24. Botstein, D.; White, R.L.; Skolnick, M.; Davis, R.W. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum. Genet. 1980, 32, 314–331. [Google Scholar] [PubMed]
  25. Qi, Z.M.; Huang, L.; Zhu, R.S.; Xin, D.W.; Liu, C.Y.; Han, X.; Jiang, H.W.; Hong, W.G.; Hu, G.H.; Zheng, H.K.; et al. A high-density genetic map for soybean based on specific length amplified fragment sequencing. PLoS ONE 2014, 9, e104871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Siddique, M.I.; Lee, J.H.; Ahn, J.H.; Kusumawardhani, M.K.; Safitri, R.; Harpenas, A.; Kwon, J.K.; Kang, B.C. Genotyping-by-sequencing-based QTL mapping reveals novel loci for Pepper yellow leaf curl virus (PepYLCV) resistance in Capsicum annuum. PLoS ONE 2022, 17, e0264026. [Google Scholar] [CrossRef]
  27. Wu, K.; Liu, H.; Yang, M.; Tao, Y.; Ma, H.; Wu, W.; Zuo, Y.; Zhao, Y. High-density genetic map construction and QTLs analysis of grain yield-related traits in sesame (Sesamum indicum L.) based on RAD-Seq techonology. BMC Plant Biol. 2014, 14, 274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Zhang, S.; Yu, Z.; Qi, X.; Wang, Z.; Zheng, Y.; Ren, H.; Liang, S.; Zheng, X. Construction of a high-density genetic map and identification of leaf trait-related QTLs in chinese bayberry (Myrica rubra). Front. Plant Sci. 2021, 12, 675855. [Google Scholar] [CrossRef] [PubMed]
  29. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef] [Green Version]
  30. Chowdry, M.A.; Slinkard, A.E. Natural outcrossing in grasspea. J. Hered. 1997, 88, 154–156. [Google Scholar] [CrossRef] [Green Version]
  31. Li, B.; Tian, L.; Zhang, J.; Huang, L.; Han, F.; Yan, S.; Wang, L.; Zheng, H.; Sun, J. Construction of a high-density genetic map based on large-scale markers developed by specific length amplified fragment sequencing (SLAF-seq) and its application to QTL analysis for isoflavone content in Glycine max. BMC Genom. 2014, 15, 1086. [Google Scholar] [CrossRef] [Green Version]
  32. Chowdhury, M.A.; Slinkard, A.E. Linkage of random amplified polymorphic DNA, isozyme and morphological markers in grasspea (Lathyrus sativus). J. Agr. Sci. 1999, 133, 389–395. [Google Scholar] [CrossRef]
  33. Skiba, B.; Ford, R.; Pang, E.C.K. Construction of a linkage map based on a Lathyrus sativus backcross population and preliminary investigation of QTLs associated with resistance to ascochyta blight. Theor. Appl. Genet. 2004, 109, 1726–1735. [Google Scholar] [CrossRef] [PubMed]
  34. Santos, C.; Polanco, C.; Rubiales, D.; Vaz Patto, M.C. The MLO1 powdery mildew susceptibility gene in Lathyrus species: The power of high-density linkage maps in comparative mapping and synteny analysis. Plant Genome 2021, 14, e20090. [Google Scholar] [CrossRef] [PubMed]
  35. Santos, C.; Almeida, N.F.; Alves, M.L.; Horres, R.; Krezdorn, N.; Leitão, S.T.; Aznar-Fernández, T.; Rotter, B.; Winter, P.; Rubiales, D.; et al. First genetic linkage map of Lathyrus cicera based on RNA sequencing-derived markers: Key tool for genetic mapping of disease resistance. Hortic. Res. 2018, 5, 45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Santos, C.; Martins, D.; Rubiales, D.; Vaz Patto, M.C. Partial resistance against Erysiphe pisi and E. trifolii under different genetic control in Lathyrus cicera: Outcomes from a linkage mapping approach. Plant Dis. 2020, 104, 2875–2884. [Google Scholar] [CrossRef] [PubMed]
  37. Pereira, L.; Ruggieri, V.; Pérez, S.; Alexiou, K.G.; Fernández, M.; Jahrmann, T.; Pujol, M.; Garcia-Mas, J. QTL mapping of melon fruit quality traits using a high-density GBS-based genetic map. BMC Plant Biol. 2018, 18, 324. [Google Scholar] [CrossRef] [PubMed]
  38. Stanton, M.L. Reproductive biology of petal color variants in wild populations of raphanus sativus: I. Pollinator response to color morphs. Am. J. Bot. 1987, 74, 178–187. [Google Scholar] [CrossRef]
  39. Yu, Y.; Yang, Z.; Jiang, Y.; Wang, L.; Wu, Y.; Liao, J.; Yang, R.; Zhang, L. Inheritance and QTL mapping for flower color in salvia miltiorrhiza bunge. J. Hered. 2022, 113, esac012. [Google Scholar] [CrossRef]
  40. Rahman, M.M.; Kumar, J.; Rahman, M.A.; Afzal, M.A. Natural outcrossing in Lathyrus satlvus L. Indian J. Genet. Plant Breed. 1995, 55, 204–207. [Google Scholar]
  41. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  42. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [Green Version]
  43. Torkamaneh, D.; Laroche, J.; Belzile, F. Genome-Wide SNP Calling from Genotyping by Sequencing (GBS) Data: A Comparison of Seven Pipelines and Two Sequencing Technologies. PLoS ONE 2016, 11, e0161333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kaur, G.; Pathak, M.; Singla, D.; Sharma, A.; Chhuneja, P.; Sarao, N.K. High-density gbs-based genetic linkage map construction and QTL identification associated with yellow mosaic disease resistance in bitter gourd (Momordica charantia L.). Front. Plant Sci. 2021, 12, 671620. [Google Scholar] [CrossRef] [PubMed]
  45. Lee, Y.R.; Kim, C.W.; Han, J.W.; Choi, H.J.; Han, K.; Lee, E.S.; Kim, D.S.; Lee, J.; Siddique, M.I.; Lee, H.E. Genotyping-by-sequencing derived genetic linkage map and quantitative trait loci for sugar content in onion (Allium cepa L.). Plants 2021, 10, 2267. [Google Scholar] [CrossRef] [PubMed]
  46. Brhane, H.; Haileselassie, T.; Tesfaye, K.; Ortiz, R.; Hammenhag, C.; Abreha, K.B.; Geleta, M. Novel gbs-based snp markers for finger millet and their use in genetic diversity analyses. Front. Genet. 2022, 13, 848627. [Google Scholar] [CrossRef]
  47. Li, L.; Yang, X.; Wang, Z.; Ren, M.; An, C.; Zhu, S.; Xu, R. Genetic mapping of powdery mildew resistance genes in wheat landrace Guizi 1 via genotyping by sequencing. Mol. Biol. Rep. 2022, 49, 4461–4468. [Google Scholar] [CrossRef]
  48. Nyirahabimana, F.; Shimira, F.; Zahid, G.; Solmaz, I. Recent status of genotyping by sequencing (GBS) technology in cucumber (Cucumis sativus L.): A review. Mol. Biol. Rep. 2022, 49, 5547–5554. [Google Scholar] [CrossRef]
  49. Jia, G.; Wang, H.; Tang, S.; Zhi, H.; Liu, S.; Wen, Q.; Qiao, Z.; Diao, X. Detection of genomic loci associated with chromosomal recombination using high-density linkage mapping in Setaria. Sci. Rep. 2017, 7, 15180. [Google Scholar] [CrossRef] [Green Version]
  50. Su, C.; Wang, W.; Gong, S.; Zuo, J.; Li, S.; Xu, S. High density linkage map construction and mapping of yield trait QTLs in maize (Zea mays) using the genotyping-by-sequencing (GBS) technology. Front. Plant Sci. 2017, 8, 706. [Google Scholar] [CrossRef] [Green Version]
  51. Zhang, Z.; Wei, T.; Zhong, Y.; Li, X.; Huang, J. Construction of a high-density genetic map of Ziziphus jujuba Mill. using genotyping by sequencing technology. Tree Genet. Genomes 2016, 12, 76. [Google Scholar] [CrossRef]
  52. Lavania, U.C. Chromosomal instability in Lathyrus sativus L. Theor. Appl. Genet. 1982, 62, 135–138. [Google Scholar] [CrossRef]
  53. Mollel, M.; Haileselassie, T. Flower colour inheritance in Nicotiana alata (Solanaceae) and its use as a genetic marker for gene flow studies. Huria. J. Open Univ. Tanzan. 2012, 12, 94–101. [Google Scholar]
  54. Nieuwhof, M.; Eijk, J.P.; Garretsen, G.; Eikelboom, W. Inheritance of flower colour in relation to inheritance of flower pigments in tulip (Tulipa L.). J. Genet. Breed. 1990, 44, 277–280. [Google Scholar]
  55. Lo, S.; Muñoz-Amatriaín, M.; Boukar, O.; Herniter, I.; Cisse, N.; Guo, Y.; Roberts, P.A.; Xu, S.; Fatokun, C.; Close, T.J. Identifcation of QTL controlling domestication-related traits in cowpea (Vigna unguiculata L. Walp). Sci. Rep. 2018, 8, 6216. [Google Scholar] [CrossRef] [Green Version]
  56. Nesi, N.; Debeaujon, I.; Jond, C.; Pelletier, G.; Caboche, M.; Lepiniec, L. The TT8 gene encodes a basic helix-loop-helix domain protein required for expression of DFR and BAN genes in Arabidopsis siliques. Plant Cell 2000, 12, 1863–1878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Bauer, A.A.; Clayton, M.K.; Brunet, J. Floral traits infl uencing plant attractiveness to three bee species: Consequences for plant reproductive success. Am. J. Bot. 2017, 104, 772–781. [Google Scholar] [CrossRef] [Green Version]
  58. Campbell, D.R.; Bischoff, M.; Lord, J.M.; Robertson, A.W. Flower color influences insect visitation in alpine New Zealand. Ecology 2010, 91, 2638–2649. [Google Scholar] [CrossRef]
  59. Ooijen, J. Multipoint maximum likelihood mapping in a full-sib family of an outbreeding species. Genet. Res. 2011, 93, 343–349. [Google Scholar] [CrossRef]
  60. Narasimha, R. Kosambi and proper orthogonal decomposition. Resonance 2011, 16, 574–581. [Google Scholar] [CrossRef]
Plants 11 02172 g001 550
Figure 1. Frequency of polymorphism segregation types between K714 and IF1347. Different colors (from left to right) represent the 7 segregation types.
Figure 1. Frequency of polymorphism segregation types between K714 and IF1347. Different colors (from left to right) represent the 7 segregation types.
Plants 11 02172 g001
Plants 11 02172 g002 550
Figure 2. Kruskal–Wallis (KW) analysis (p < 0.01). The grey arrow and blue point indicate the most significant value. Chromosomes (left to right, from 1 to 7) are represented with different colors along the x-axis.
Figure 2. Kruskal–Wallis (KW) analysis (p < 0.01). The grey arrow and blue point indicate the most significant value. Chromosomes (left to right, from 1 to 7) are represented with different colors along the x-axis.
Plants 11 02172 g002
Plants 11 02172 g003 550
Figure 3. Interval mapping (IM) test (LOD > 3.0). The grey arrow and green point indicate the most significant value. Chromosomes (left to right) are represented from 1 to 7 along the x-axis.
Figure 3. Interval mapping (IM) test (LOD > 3.0). The grey arrow and green point indicate the most significant value. Chromosomes (left to right) are represented from 1 to 7 along the x-axis.
Plants 11 02172 g003
Table
Table 1. Grasspea genetic linkage map statistical information.
Table 1. Grasspea genetic linkage map statistical information.
LG aNumber of SNPsGenetic Distance (cM)/LGAvg. Genetic Distance (cM)/LGMax. Gap (cM)/LG
LG112401367.471.110.41
LG24761328.022.7929.15
LG3354751.702.1218.29
LG4290504.111.7431.47
LG5367924.762.5229.16
LG6427959.592.2518.42
LG73821140.032.9825.22
Total35366975.682.2131.47
a: grasspea linkage group (LG).
Table
Table 2. Statistics regarding the gap between markers across the grasspea genetic linkage map.
Table 2. Statistics regarding the gap between markers across the grasspea genetic linkage map.
LG aNumber of SNPs in <5 cMNumber of SNPs in 5–10 cMNumber of SNPs in 10–20 cMNumber of SNPs in >20 cMNumber of SNPs in <5 cM Ratio (%)
LG11228111099.03
LG23995122483.82
LG33162513089.27
LG4272115293.79
LG53063920283.38
LG63803512088.99
LG73046017179.58
Total320523290990.64
a: grasspea linkage group (LG).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hao, X.; Yang, T.; Wang, Y.; Liu, R.; Dong, X.; Zhao, J.; Han, J.; Zong, X.; Chang, J.; Liu, H. Construction of A GBS-Based High-Density Genetic Map and Flower Color-Related Loci Mapping in Grasspea (Lathyrus sativus L.). Plants 2022, 11, 2172. https://doi.org/10.3390/plants11162172

AMA Style

Hao X, Yang T, Wang Y, Liu R, Dong X, Zhao J, Han J, Zong X, Chang J, Liu H. Construction of A GBS-Based High-Density Genetic Map and Flower Color-Related Loci Mapping in Grasspea (Lathyrus sativus L.). Plants. 2022; 11(16):2172. https://doi.org/10.3390/plants11162172

Chicago/Turabian Style

Hao, Xiaopeng, Tao Yang, Yan Wang, Rong Liu, Xue Dong, Jiandong Zhao, Jucai Han, Xuxiao Zong, Jianwu Chang, and Huiping Liu. 2022. "Construction of A GBS-Based High-Density Genetic Map and Flower Color-Related Loci Mapping in Grasspea (Lathyrus sativus L.)" Plants 11, no. 16: 2172. https://doi.org/10.3390/plants11162172

APA Style

Hao, X., Yang, T., Wang, Y., Liu, R., Dong, X., Zhao, J., Han, J., Zong, X., Chang, J., & Liu, H. (2022). Construction of A GBS-Based High-Density Genetic Map and Flower Color-Related Loci Mapping in Grasspea (Lathyrus sativus L.). Plants, 11(16), 2172. https://doi.org/10.3390/plants11162172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop