SLAF-seq: An Efficient Method of Large-Scale De Novo SNP Discovery and Genotyping Using High-Throughput Sequencing

Large-scale genotyping plays an important role in genetic association studies. It has provided new opportunities for gene discovery, especially when combined with high-throughput sequencing technologies. Here, we report an efficient solution for large-scale genotyping. We call it specific-locus amplified fragment sequencing (SLAF-seq). SLAF-seq technology has several distinguishing characteristics: i) deep sequencing to ensure genotyping accuracy; ii) reduced representation strategy to reduce sequencing costs; iii) pre-designed reduced representation scheme to optimize marker efficiency; and iv) double barcode system for large populations. In this study, we tested the efficiency of SLAF-seq on rice and soybean data. Both sets of results showed strong consistency between predicted and practical SLAFs and considerable genotyping accuracy. We also report the highest density genetic map yet created for any organism without a reference genome sequence, common carp in this case, using SLAF-seq data. We detected 50,530 high-quality SLAFs with 13,291 SNPs genotyped in 211 individual carp. The genetic map contained 5,885 markers with 0.68 cM intervals on average. A comparative genomics study between common carp genetic map and zebrafish genome sequence map showed high-quality SLAF-seq genotyping results. SLAF-seq provides a high-resolution strategy for large-scale genotyping and can be generally applicable to various species and populations.     

High-density genetic map of <Emphasis Type="Italic">Populus deltoides</Emphasis> constructed by using specific length amplified fragment sequencing

Populus deltoides is an important industrial tree species widely planted in many areas of the world, and de novo genome sequencing of this plant has been carried out by several research groups. A dense genetic map associating genome sequences is highly desirable for reconstructing the chromosome pseudomolecules using the obtained sequence scaffold assemblies. For this purpose, an intraspecific full-sib F₁ mapping pedigree was established by using the sequenced P. deltoides as the maternal parent. With this mapping pedigree, we constructed a high-density genetic map using 92 randomly selected progenies. Single nucleotide polymorphism (SNP) markers were discovered by using specific length amplified fragment sequencing (SLAF-seq). In total, 487,038 SLAFs were generated, of which 179,360 were polymorphic. A high-density genetic map was built using HighMap software, which included 11,680 Mendelian segregation markers distributing in 2851 marker bins. The established map consisted of 19 linkage groups (LGs) that equaled to the 19 haploid chromosomes in poplar genome, and spanned a total genetic length of 3494.66 cM, with an average distance of 1.23 cM per marker bin. The map presented here will be useful for anchoring the genome sequence assemblies along chromosomes, and for many other aspects of genetic studies on P. deltoides and the closely related species.     

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

Background

Quantitative trait locus (QTL) mapping is an efficient approach to discover the genetic architecture underlying complex quantitative traits. However, the low density of molecular markers in genetic maps has limited the efficiency and accuracy of QTL mapping. In this study, specific length amplified fragment sequencing (SLAF-seq), a new high-throughput strategy for large-scale SNP discovery and genotyping based on next generation sequencing (NGS), was employed to construct a high-density soybean genetic map using recombinant inbred lines (RILs, Luheidou2 × Nanhuizao, F_5:8). With this map, the consistent QTLs for isoflavone content across various environments were identified.

Results

In total, 23 Gb of data containing 87,604,858 pair-end reads were obtained. The average coverage for each SLAF marker was 11.20-fold for the female parent, 12.51-fold for the male parent, and an average of 3.98-fold for individual RILs. Among the 116,216 high-quality SLAFs obtained, 9,948 were polymorphic. The final map consisted of 5,785 SLAFs on 20 linkage groups (LGs) and spanned 2,255.18 cM in genome size with an average distance of 0.43 cM between adjacent markers. Comparative genomic analysis revealed a relatively high collinearity of 20 LGs with the soybean reference genome. Based on this map, 41 QTLs were identified that contributed to the isoflavone content. The high efficiency and accuracy of this map were evidenced by the discovery of genes encoding isoflavone biosynthetic enzymes within these loci. Moreover, 11 of these 41 QTLs (including six novel loci) were associated with isoflavone content across multiple environments. One of them, qIF20-2, contributed to a majority of isoflavone components across various environments and explained a high amount of phenotypic variance (8.7% - 35.3%). This represents a novel major QTL underlying isoflavone content across various environments in soybean.

Conclusions

Herein, we reported a high-density genetic map for soybean. This map exhibited high resolution and accuracy. It will facilitate the identification of genes and QTLs underlying essential agronomic traits in soybean. The novel major QTL for isoflavone content is useful not only for further study on the genetic basis of isoflavone accumulation, but also for marker-assisted selection (MAS) in soybean breeding in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1086) contains supplementary material, which is available to authorized users.     

Construction of the First High-Density Genetic Linkage Map and Analysis of Quantitative Trait Loci for Growth-Related Traits in <Emphasis Type="Italic">Sinonovacula constricta</Emphasis>

The razor clam (Sinonovacula constricta) is an important aquaculture species, for which a high-density genetic linkage map would play an important role in marker-assisted selection (MAS). In this study, we constructed a high-density genetic map and detected quantitative trait loci (QTLs) for Sinonovacula constricta with an F₁ cross population by using the specific locus amplified fragment sequencing (SLAF-seq) method. A total of 315,553 SLAF markers out of 467.71 Mreads were developed. The final linkage map was composed of 7516 SLAFs (156.60-fold in the parents and 20.80-fold in each F₁ population on average). The total distance of the linkage map was 2383.85 cM, covering 19 linkage groups with an average inter-marker distance of 0.32 cM. The proportion of gaps less than 5.0 cM was on average 96.90%. A total of 16 suggestive QTLs for five growth-related traits (five QTLs for shell height, six QTLs for shell length, three QTLs for shell width, one QTL for total body weight, and one QTL for soft body weight) were identified. These QTLs were distributed on five linkage groups, and the regions showed overlapping on LG9 and LG13. In conclusion, the high-density genetic map and QTLs for S. constricta provide a valuable genetic resource and a basis for MAS.     

The First High-Density Genetic Map Construction in Tree Peony (Paeonia Sect. Moutan) using Genotyping by Specific-Locus Amplified Fragment Sequencing

Genetic linkage maps, permitting the elucidation of genome structure, are one of most powerful genomic tools to accelerate marker-assisted breeding. However, due to a lack of sufficient user-friendly molecular markers, no genetic linkage map has been developed for tree peonies (Paeonia Sect. Moutan), a group of important horticultural plants worldwide. Specific-locus amplified fragment sequencing (SLAF-seq) is a recent molecular marker development technology that enable the large-scale discovery and genotyping of sequence-based marker in genome-wide. In this study, we performed SLAF sequencing of an F1 population, derived from the cross P. ostti ‘FenDanBai’ × P. × suffruticosa ‘HongQiao’, to identify sufficient high-quality markers for the construction of high-density genetic linkage map in tree peonies. After SLAF sequencing, a total of 78 Gb sequencing data and 285,403,225 pair-end reads were generated. We detected 309,198 high-quality SLAFs from these data, of which 85,124 (27.5%) were polymorphic. Subsequently, 3518 of the polymorphic markers, which were successfully encoded in to Mendelian segregation types, and were in conformity with the criteria of high-quality markers, were defined as effective markers and used for genetic linkage mapping. Finally, we constructed an integrated genetic map, which comprised 1189 markers on the five linkage groups, and spanned 920.699 centiMorgans (cM) with an average inter-marker distance of 0.774 cM. There were 1115 ‘SNP-only’ markers, 18 ‘InDel-only’ markers, and 56 ‘SNP&InDel’ markers on the map. Among these markers, 450 (37.85%) showed significant segregation distortion (P < 0.05). In conclusion, this investigation reported the first large-scale marker development and high-density linkage map construction for tree peony. The results of this study will serve as a solid foundation not only for marker-assisted breeding, but also for genome sequence assembly for tree peony.     

Construction of a high-density genetic map using specific length amplified fragment markers and identification of a quantitative trait locus for anthracnose resistance in walnut (Juglans regia L.)

Background

Walnut (Juglans regia, 2n = 32, approximately 606 Mb per 1C genome) is an economically important tree crop. Resistance to anthracnose, caused by Colletotrichum gloeosporioides, is a major objective of walnut genetic improvement in China. The recently developed specific length amplified fragment sequencing (SLAF-seq) is an efficient strategy that can obtain large numbers of markers with sufficient sequence information to construct high-density genetic maps and permits detection of quantitative trait loci (QTLs) for molecular breeding.

Results

SLAF-seq generated 161.64 M paired-end reads. 153,820 SLAF markers were obtained, of which 49,174 were polymorphic. 13,635 polymorphic markers were sorted into five segregation types and 2,577 markers of them were used to construct genetic linkage maps: 2,395 of these fell into 16 linkage groups (LGs) for the female map, 448 markers for the male map, and 2,577 markers for the integrated map. Taking into account the size of all LGs, the marker coverage was 2,664.36 cM for the female map, 1,305.58 cM for the male map, and 2,457.82 cM for the integrated map. The average intervals between two adjacent mapped markers were 1.11 cM, 2.91 cM and 0.95 cM for three maps, respectively. ‘SNP_only’ markers accounted for 89.25 % of the markers on the integrated map. Mapping markers contained 5,043 single nucleotide polymorphisms (SNPs) loci, which corresponded to two SNP loci per SLAF marker. According to the integrated map, we used interval mapping (Logarithm of odds, LOD > 3.0) to detect our quantitative trait. One QTL was detected for anthracnose resistance. The interval of this QTL ranged from 165.51 cM to 176.33 cM on LG14, and ten markers in this interval that were above the threshold value were considered to be linked markers to the anthracnose resistance trait. The phenotypic variance explained by each marker ranged from 16.2 to 19.9 %, and their LOD scores varied from 3.22 to 4.04.

Conclusions

High-density genetic maps for walnut containing 16 LGs were constructed using the SLAF-seq method with an F1 population. One QTL for walnut anthracnose resistance was identified based on the map. The results will aid molecular marker-assisted breeding and walnut resistance genes identification.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1822-8) contains supplementary material, which is available to authorized users.     

Large-Scale SNP Discovery and Genotyping for Constructing a High-Density Genetic Map of Tea Plant Using Specific-Locus Amplified Fragment Sequencing (SLAF-seq)

Genetic maps are important tools in plant genomics and breeding. The present study reports the large-scale discovery of single nucleotide polymorphisms (SNPs) for genetic map construction in tea plant. We developed a total of 6,042 valid SNP markers using specific-locus amplified fragment sequencing (SLAF-seq), and subsequently mapped them into the previous framework map. The final map contained 6,448 molecular markers, distributing on fifteen linkage groups corresponding to the number of tea plant chromosomes. The total map length was 3,965 cM, with an average inter-locus distance of 1.0 cM. This map is the first SNP-based reference map of tea plant, as well as the most saturated one developed to date. The SNP markers and map resources generated in this study provide a wealth of genetic information that can serve as a foundation for downstream genetic analyses, such as the fine mapping of quantitative trait loci (QTL), map-based cloning, marker-assisted selection, and anchoring of scaffolds to facilitate the process of whole genome sequencing projects for tea plant.     

High-density genetic map construction and identification of a locus controlling weeping trait in an ornamental woody plant (Prunus mume Sieb. et Zucc)

High-density genetic map is a valuable tool for fine mapping locus controlling a specific trait especially for perennial woody plants. In this study, we firstly constructed a high-density genetic map of mei (Prunus mume) using SLAF markers, developed by specific locus amplified fragment sequencing (SLAF-seq). The linkage map contains 8,007 markers, with a mean marker distance of 0.195 cM, making it the densest genetic map for the genus Prunus. Though weeping trees are used worldwide as landscape plants, little is known about weeping controlling gene(s) (Pl). To test the utility of the high-density genetic map, we did fine-scale mapping of this important ornamental trait. In total, three statistic methods were performed progressively based on the result of inheritance analysis. Quantitative trait loci (QTL) analysis initially revealed that a locus on linkage group 7 was strongly responsible for weeping trait. Mutmap-like strategy and extreme linkage analysis were then applied to fine map this locus within 1.14 cM. Bioinformatics analysis of the locus identified some candidate genes. The successful localization of weeping trait strongly indicates that the high-density map constructed using SLAF markers is a worthy reference for mapping important traits for woody plants.     

Construction of a high-density SNP genetic map in flue-cured tobacco based on SLAF-seq

Tobacco (Nicotiana tabacum L., 2n = 48) is an important agronomic crop and model plant. Flue-cured tobacco is the most important type and accounts for approximately 80 % of tobacco production worldwide. The low genetic diversity of flue-cured tobacco impedes the construction of a high-density genetic linkage map using simple sequence repeat (SSR) markers and warrants the exploitation of single nucleotide polymorphic (SNP) markers from genomic regions. In this article, initially using specific locus-amplified fragment sequencing, we discovered 10,891 SNPs that were subsequently used as molecular markers for genetic map construction. Combined with SSR markers, a final high-density genetic map was generated containing 4215 SNPs and 194 SSRs distributed on 24 linkage groups (LGs). The genetic map was 2662.43 cM in length, with an average distance of 0.60 cM between adjacent markers. Furthermore, by mapping the SNP markers to the ancestral genomes of Nicotiana tomentosiformis and Nicotiana sylvestris, a large number of genome rearrangements were identified as occurring after the polyploidization event. Finally, using this novel integrated map and mapping population, two major quantitative trait loci (QTLs) were identified for flue-curing and mapped to the LG6 of tobacco. This is the first report of SNP markers and a SNP-based linkage map being developed in tobacco. The high-density genetic map and QTLs related to tobacco curing will support gene/QTL fine mapping, genome sequence assembly and molecular breeding in tobacco.     

The linkage maps of Dendrobium species based on RAPD and SRAP markers

Dendrobium plants are used commonly as tonic herbs and health food in many Asian countries,especially in China.Here we report the genetic map construction of two Dendrobium species with a double pseudo-testcross strategy using random amplified polymorphic DNA (RAPD) and sequence-related amplified polymorphism (SRAP) markers.A F1 mapping population of 90 individuals was developed from a cross between D.officinale and D.hercoglossum.A total of 307 markers,including 209 RAPD and 98 SRAP,were identified and used for genetic linkage group (LG) analysis.The D.officinale linkage map consisted of 11 major linkage groups and 3 doublets,which covered 629.4 cM by a total of 62 markers with an average locus distance of 11.2 cM between two adjacent markers.The D.hercoglossum linkage map contained 112 markers mapped on 15 major and 4 minor linkage groups,spanning a total length of 1,304.6 cM with an average distance of 11.6 cM between two adjacent markers.The maps constructed in this study covered 92.7% and 82.7% of the D.hercoglossum and D.officinale genomes respectively,providing an important basis for the mapping of horticultural and medicinal traits and for the application of marker-assisted selection in Dendrobium breeding program.     

The first genetic linkage map of Luohanguo (Siraitia grosvenorii ) based on ISSR and SRAP markers

In this study, the first genetic map of Luohanguo (Siraitia grosvenorii (Swingle) C. Jeffrey) was constructed with 150 F? population individuals using inter-simple sequence repeat (ISSR) and sequence-related amplified polymorphism (SRAP) markers. A total of 100 ISSRs and 196 SRAP primer combinations generated 51 and 222 polymorphic markers, respectively. Among the 273 markers obtained, 199 markers (29 ISSRs and 170 SRAPs) were mapped to 25 linkage groups. The map covered 1463.3 cM with a mean map distance of 7.35 cM between adjacent markers and a maximum map distance of 52.6 cM between two markers. The markers were distributed randomly in 25 groups except for minor clusters in the distal region of linkage groups. All 25 linkage groups consisted of 2-36 loci ranging in length from 19.5 to 152.6 cM and accounted for 59.8% of the total map distance. This map provides reference information for future molecular breeding work on Luohanguo.     

AFLP-based genetic linkage map for the red flour beetle (Tribolium castaneum)

The red flour beetle (Tribolium castaneum) is a major pest of stored grain and grain products and a popular model species for a variety of ecological, evolutionary, and developmental biology studies. Development of a linkage map based on reproducible and highly polymorphic molecular markers would greatly facilitate research in these disciplines. We have developed a genetic linkage map using 269 amplified fragment length polymorphism (AFLP) markers. Ten previously known random amplified polymorphic DNA (RAPD) markers were used as anchor markers for linkage group assignment. The linkage map was constructed through genotyping two independent F(2) segregating populations with 48 AFLP primer combinations. Each primer combination generated an average of 4.6 AFLP markers eligible for linkage mapping. The length of the integrated map is 573 cM, giving an average marker resolution of 2.0 cM and an average physical distance per genetic distance of 350 kb/cM. A cluster of loci on linkage group 3 exhibited significant segregation distortion. We have also identified six X-linked and two Y-linked markers. Five mapped AFLP fragments were sequenced and converted to sequence-tagged site (STS) markers.     

大豆重组自交系群体（ZKS-HX）遗传图谱的构建和形态标记的定位

以大豆品种‘合丰25’为母本,半野生大豆‘新民6号’为父本杂交得到的F2-9代122个重组自交系为试验材料,构建了含有124个SSR标记、1个EST标记、3个形态学标记的大豆遗传图谱。此图谱覆盖的基因组长度为2348．3cM．标记间平均距离为18．3cM。每个连锁群长度范围为15．1～195．9cM之间,标记数范围2—10个。本文将控制茸毛色（Pb）基因定位于LG06-C2连锁群上,与Sat_40x2的遗传距离为39．6cM;控制叶耳g（Le）、花色（4W,）基因定位于LG12-F连锁群上,它们之间的遗传距离为9．9cM,与两边的Satt348、Sat_240标记遗传距离分别为13．3cM和10．5cM。     

High-density linkage mapping revealed suppression of recombination at the sex determination locus in papaya

A high-density genetic map of papaya (Carica papaya L.) was constructed using 54 F(2) plants derived from cultivars Kapoho and SunUp with 1501 markers, including 1498 amplified fragment length polymorphism (AFLP) markers, the papaya ringspot virus coat protein marker, morphological sex type, and fruit flesh color. These markers were mapped into 12 linkage groups at a LOD score of 5.0 and recombination frequency of 0.25. The 12 major linkage groups covered a total length of 3294.2 cM, with an average distance of 2.2 cM between adjacent markers. This map revealed severe suppression of recombination around the sex determination locus with a total of 225 markers cosegregating with sex types. The cytosine bases were highly methylated in this region on the basis of the distribution of methylation-sensitive and -insensitive markers. This high-density genetic map is essential for cloning of specific genes of interest such as the sex determination gene and for the integration of genetic and physical maps of papaya.     

RAPD-based genetic linkage maps of Tribolium castaneum.

A genetic map of the red flour beetle (Tribolium castaneum) integrating molecular with morphological markers was constructed using a backcross population of 147 siblings. The map defines 10 linkage groups (LGs), presumably corresponding to the 10 chromosomes, and consists of 122 randomly amplified polymorphic DNA (RAPD) markers, six molecular markers representing identified genes, and five morphological markers. The total map length is 570 cM, giving an average marker resolution of 4.3 cM. The average physical distance per genetic distance was estimated at 350 kb/cM. A cluster of loci showing distorted segregation was detected on LG9. The process of converting RAPD markers to sequence-tagged site markers was initiated: 18 RAPD markers were cloned and sequenced, and single-strand conformational polymorphisms were identified for 4 of the 18. The map positions of all 4 coincided with those of the parent RAPD markers.     

A genetic linkage map construction for sesame (Sesamum indicum L.)

Sesame (Sesamum indicum L.) is one of the oldest oilseed crops with high seed oil quality. The first sesame genetic linkage map based on F2 segregating population of an intraspecific cross between two cultivars was constructed. Using three types of PCR-based markers, 284 polymorphic loci including 10 EST-SSR marker, 30 AFLP marker and 244 RSAMPL marker, respectively, had been screened. Subsequently, a total of 220 molecular markers were mapped in 30 linkage groups covering a genetic length of 936.72 cM, and the average distance between markers was 4.93 cM. In this map, the linkage groups contained from 2 to 33 loci each and ranged in distance from 6.44 cM to 74.52 cM. Based on map information, sesame genome length was estimated to be approximately 1,232.53 cM, and genome coverage of this map was about 76.0%. As a starting point of sesame genome study, the genetic linkage map will be hopeful to tag traits of breeding interest and further aid in the sesame molecular breeding. Furthermore, RSAMPL marker had been also appreciated in this paper, for its first usage in genetic map construction and higher utilization potential in some crop species lacking much genome information.     

Development of microsatellite markers and construction of genetic map in rice blast pathogen Magnaporthe grisea

Magnaporthe grisea is the most destructive fungal pathogen of rice and a model organism for studying plant-pathogen interaction. Molecular markers and genetic maps are useful tools for genetic studies. In this study, based on the released genome sequence data of M. grisea, we investigated 446 simple sequence repeat (SSR) loci and developed 313 SSR markers, which showed polymorphisms among nine isolates from rice (including a laboratory strain 2539). The number of alleles of each marker ranged 2-9 with an average of 3.3. The polymorphic information content (PIC) of each marker ranged 0.20-0.89 with an average of 0.53. Using a population derived from a cross between isolates Guy11 and 2539, we constructed a genetic map of M. grisea consisting of 176 SSR markers. The map covers a total length of 1247 cM, equivalent to a physical length of about 35.0 Mb or 93% of the genome, with an average distance of 7.1cM between adjacent markers. A web-based database of the SSR markers and the genetic map was established (http://ibi.zju.edu.cn/pgl/MGM/index.html).     

Construction of an SSR and RAD Marker-Based Genetic Linkage Map for Carnation (<Emphasis Type="Italic">Dianthus caryophyllus</Emphasis> L.)

A high-density genetic map, an essential tool for comparative genomic studies and quantitative trait locus fine mapping, can also facilitate genome sequence assembly. The sequence-based marker technology known as restriction site-associated DNA (RAD) enables synchronous, single nucleotide polymorphism marker discovery, and genotyping using massively parallel sequencing. We constructed a high-density linkage map for carnation (Dianthus caryophyllus L.) based on simple sequence repeat (SSR) markers in combination with RAD markers developed by double-digest RAD sequencing (ddRAD-seq). A total of 2404 (285 SSR and 2119 RAD) markers could be assigned to 15 linkage groups spanning 971.5 cM, with an average marker interval of 0.4 cM. The total length of scaffolds with identified map positions was 95.6 Mb, which is equivalent to 15.4 % of the estimated genome size. The generated map is the first SSR and RAD marker-based high-density linkage map reported for carnation. The ddRAD-seq pipeline developed in this study should also help accelerate genetic and genomics analyses and molecular breeding of carnation and other non-model crops.