A high-density genetic map of cucumber derived from Specific Length Amplified Fragment sequencing (SLAF-seq)

Construction of the first high-density genetic linkage map and QTL mapping of flavonoid and leaf-size related traits in Epimedium

BACKGROUND

Leaves are the main medicinal organ in Epimedium herbs, and leaf flavonoid content is an important criterion of Epimedium herbs. However, the underlying genes that regulate leaf size and flavonoid content are unclear, which limits the use of breeding for Epimedium development. This study focuses on QTL mapping of flavonoid and leaf-size related traits in Epimedium.

RESULTS

We constructed the first high-density genetic map (HDGM) using 109 F1 hybrids of Epimedium leptorrhizum and Epimedium sagittatum over three years (2019-2021). Using 5,271 single nucleotide polymorphism (SNP) markers, an HDGM with an overall distance of 2,366.07 cM and a mean gap of 0.612 cM was generated by utilizing genotyping by sequencing (GBS) technology. Every year for three years, 46 stable quantitative trait loci (QTLs) for leaf size and flavonoid contents were discovered, including 31 stable loci for Epimedin C (EC), one stable locus for total flavone content (TFC), 12 stable loci for leaf length (LL), and two stable loci for leaf area (LA). For flavonoid content and leaf size, the phenotypic variance explained for these loci varied between 4.00 and 16.80% and 14.95 and 17.34%, respectively.

CONCLUSIONS

Forty-six stable QTLs for leaf size and flavonoid content traits were repeatedly detected over three years. The HDGM and stable QTLs are laying the basis for breeding and gene investigation in Epimedium and will contribute to accelerating the identification of desirable genotypes for Epimedium breeding.

Construction of a High-Density Recombination Bin-Based Genetic Map Facilitates High-Resolution Mapping of a Major QTL Underlying Anthocyanin Pigmentation in Eggplant

High-density genetic maps can significantly improve the resolution of QTL mapping. We constructed a high-density recombination bin-based genetic map of eggplant based on 200 F₂ plants from an interspecific cross (Solanum melongena × S. incanum) using the whole genome resequencing strategy. The map was 2022.8 cM long, covering near 99% of the eggplant genome. The map contained 3776 bins, with 3644 (96.5%) being effective (position non-redundant) ones, giving a nominal average distance of 0.54 cM and an effective average distance of 0.56 cM between adjacent bins, respectively. Using this map and 172 F_2:3 lines, a major QTL with pleiotropic effects on two anthocyanin pigmentation-related traits, leaf vein color (LVC) and fruit pericarp color (FPC), was steadily detected in a bin interval of 2.28 cM (or 1.68 Mb) on chromosome E10 in two cropping seasons, explaining ~65% and 55% of the phenotypic variation in LVC and FPC, respectively. Genome-wide association analysis in this population validated the QTL and demonstrated the correctness of mapping two bins of chromosome E02 onto E10. Bioinformatics analysis suggested that a WDR protein gene inside the bin interval with reliable effective variation between the two parents could be a possible candidate gene of the QTL.

Construction of High-Density Genetic Map and Identification of a Bruchid Resistance Locus in Mung Bean (Vigna radiata L.)

Mung bean (Vigna radiata L.) is an economically important grain legume cultivated in Asian countries. High-density genetic linkage is a valuable and effective tool for mapping quantitative trait loci (QTL). In the current study, a high-resolution genetic map containing 4,180 single-nucleotide polymorphisms (SNPs) was assigned to 11 linkage groups (LGs) and spanning 1,751.39 cM in length was constructed for mung bean, and the average distance between adjacent markers was 0.42 cM. Bruchids (Callosobruchus spp.) cause significant damage to and loss of legume seeds. A locus for bruchid resistance was detected. The gene Vradi05g03810, encoding a probable resistance-specific protein, was found to be the most likely key candidate gene in mung beans. A 69-bp sequence deletion was identified in the coding region by comparing the cDNA sequences of bruchid-resistant and bruchid-susceptible lines. This SNP-based high-density linkage map is one of the first to be constructed across the mung bean genome. This map will not only facilitate the genetic mapping of genes or complex loci that control important agronomic traits but also offer a tool for promoting future genetics and comparative genomic studies in Vigna.

Construction of a high-density genetic map based on specific-locus amplified fragment sequencing and identification of loci controlling anthocyanin pigmentation in Yunnan red radish

Radish (Raphanus sativus L.) belongs to the family Brassicaceae. The Yunnan red radish variety contains fairly relatively large amounts of anthocyanins, making them important raw materials for producing edible red pigment. However, the genetic mechanism underlying this pigmentation has not been fully characterized. Herein, the radish inbred line YAAS-WR1 (white root-skin and white root-flesh) was crossed with the inbred line YAAS-RR1 (red root-skin and red root-flesh) to produce F1, F2, BC1P1, and BC1P2 populations. Genetic analyses revealed that the pigmented/non-pigmented (PiN) and purple/red (PR) traits were controlled by two genetic loci. The F2 population and the specific-locus amplified fragment sequencing (SLAF-seq) technique were used to construct a high-density genetic map (1230.16 cM), which contained 4032 markers distributed in nine linkage groups, with a mean distance between markers of 0.31 cM. Additionally, two QTL (QAC1 and QAC2) considerably affecting radish pigmentation were detected. A bioinformatics analysis of the QAC1 region identified 58 predicted protein-coding genes. Of these genes, RsF3'H, which is related to anthocyanin biosynthesis, was revealed as a likely candidate gene responsible for the PR trait. The results were further verified by analyzing gene structure and expression. Regarding QAC2, RsMYB1.3 was determined to be a likely candidate gene important for the PiN trait, with a 4-bp insertion in the first exon that introduced a premature termination codon in the YAAS-WR1 sequence. Assays demonstrated that RsMYB1.3 interacted with RsTT8 and activates RsTT8 and RsUFGT expression. These findings may help clarify the complex regulatory mechanism underlying radish anthocyanin synthesis. Furthermore, this study's results may be relevant for the molecular breeding of radish to improve the anthocyanin content and appearance of the taproots.

Identification of quantitative trait loci for powdery mildew resistance in highly resistant cucumber (Cucumis sativus L.) using ddRAD-seq analysis

Powdery mildew, caused by Podosphaera xanthii (syn. Sphaerotheca fuliginea ex Fr. Poll.), is one of the most economically important foliar diseases in cucumber (Cucumis sativus L.). Cucumber parental line 'Kyuri Chukanbohon Nou 5 Go', developed from weedy cucumber line CS-PMR1, is highly resistant to powdery mildew and is promising breeding material. We performed quantitative trait locus (QTL) analysis using double-digest restriction-site-associated DNA sequencing (ddRAD-Seq) in a population from a cross between 'Kyuri Chukanbohon Nou 5 Go' and the Japanese native cultivar 'Kaga-aonaga-fushinari', which is susceptible to powdery mildew. The resistance of the population and its parents was evaluated using leaf disc assays and image analysis. We detected one major QTL on Chr. 5 that was effective at both 20°C and 25°C and one minor QTL on Chr. 1 effective at 20°C. We detected two additional QTLs in subpopulation: one on Chr. 3 effective at 20°C and one on Chr. 5 effective at both 20°C and 25°C in a position different from the major QTL. The resistance alleles at all four QTLs were contributed by 'Kyuri Chukanbohon Nou 5 Go'. The results of this study can be used to develop practical DNA markers tightly linked to genes for powdery mildew resistance.

Genome-wide association study of body size traits in Wenshang Barred chickens based on the specific-locus amplified fragment sequencing technology

Chicken body size (BS) is an economically important trait, which has been assessed in many studies for genetic selection. However, previous reports detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) of purebred Wenshang Barred chickens. A total of 250 one-day-old male chickens were assessed in this study. Body size in individual birds was measured at 56 days. SLAF-seq was used to genotype and GWAS analysis was carried out using the general linear model (GLM) of the TASSEL program. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, of which 175,211 were tested as candidate SNPs for genome-wide association analysis using the TASSEL general linear model. Three SNPs markers reached genome-wide significance. Of these, chrZ:81729634, chrZ:81841715, and chrZ:81954149 at 81,729,634, 81,841,715, and 81,954,149 bp of GGA Z were significantly associated with body diagonal length at 56 days (BDL56); and tibia length at 56 days (TL56). These SNPs were close to three genes, including ZCCHC7, PAX5, and MELK. These results open new horizons for studies on BS and should promote the use of Chinese chickens, especially Wenshang Barred chickens.

Genome-Wide Association Study of Root and Shoot Related Traits in Spring Soybean (Glycine max L.) at Seedling Stages Using SLAF-Seq

Root systems can display variable genetic architectures leading to nutrient foraging or improving abiotic stress tolerance. Breeding for new soybean varieties with efficient root systems has tremendous potential in enhancing resource use efficiency and plant adaptation for challenging climates. In this study, root related traits were analyzed in a panel of 260 spring soybean with genome-wide association study (GWAS). Genotyping was done with specific locus amplified fragment sequencing (SLAF-seq), and five GWAS models (GLM, MLM, CMLM, FaST-LMM, and EMMAX) were used for analysis. A total of 179,960 highly consistent SNP markers distributed over the entire genome with an inter-marker distance of 2.36 kb was used for GWAS analysis. Overall, 27 significant SNPs with a phenotypic contribution ranging from 20 to 72% and distributed on chromosomes 2, 6, 8, 9, 13, 16 and 18 were identified and two of them were found to be associated with multiple root-related traits. Based on the linkage disequilibrium (LD) distance of 9.5 kb for the different chromosomes, 11 root and shoot regulating genes were detected based on LD region of a maximum 55-bp and phenotypic contribution greater than 22%. Expression analysis revealed an association between expression levels of those genes and the degree of root branching number. The current study provides new insights into the genetic architecture of soybean roots, and the underlying SNPs/genes could be critical for future breeding of high-efficient root system in soybean.

QTL Mapping for Leaf Area of Tea Plants (Camellia sinensis) Based on a High-Quality Genetic Map Constructed by Whole Genome Resequencing

High-quality genetic maps play important roles in QTL mapping and molecular marker-assisted breeding. Tea leaves are not only important vegetative organs but are also the organ for harvest with important economic value. However, the key genes and genetic mechanism of regulating leaf area have not been clarified. In this study, we performed whole-genome resequencing on "Jinxuan," "Yuncha 1" and their 96 F1 hybrid offspring. From the 1.84 Tb of original sequencing data, abundant genetic variation loci were identified, including 28,144,625 SNPs and 2,780,380 indels. By integrating the markers of a previously reported genetic map, a high-density genetic map consisting of 15 linkage groups including 8,956 high-quality SNPs was constructed. The total length of the genetic map is 1,490.81 cM, which shows good collinearity with the genome. A total of 25 representative markers (potential QTLs) related to leaf area were identified, and there were genes differentially expressed in large and small leaf samples near these markers. GWAS analysis further verified the reliability of QTL mapping. Thirty-one pairs of newly developed indel markers located near these potential QTLs showed high polymorphism and had good discrimination between large and small leaf tea plant samples. Our research will provide necessary support and new insights for tea plant genetic breeding, quantitative trait mapping and yield improvement.

Genetic diversity and population structure of eddoe taro in China using genome-wide SNP markers

Taro (Colocasia esculenta) is an important root and tuber crop cultivated worldwide. There are two main types of taro that vary in morphology of corm and cormel, ‘dasheen’ and ‘eddoe’. The eddoe type (Colocasia esculenta var. antiquorium) is predominantly distributed throughout China. Characterizing the genetic diversity present in the germplasm bank of taro is fundamental to better manage, conserve and utilize the genetic resources of this species. In this study, the genetic diversity of 234 taro accessions from 16 provinces of China was assessed using 132,869 single nucleotide polymorphism (SNP) markers identified by specific length amplified fragment-sequencing (SLAF-seq). Population structure and principal component analysis permitted the accessions to be categorized into eight groups. The genetic diversity and population differentiation of the eight groups were evaluated using the characterized SNPs. Analysis of molecular variance showed that the variation among eight inferred groups was higher than that within groups, while a relatively small variance was found among the two morphological types and 16 collection regions. Further, a core germplasm set comprising 41 taro accessions that maintained the genetic diversity of the entire collection was developed based on the genotype. This research is expected to be valuable for genetic characterization, germplasm conservation, and breeding of taro.

Novel stable QTLs identification for berry quality traits based on high-density genetic linkage map construction in table grape

BACKGROUND

Aroma, berry firmness and berry shape are three main quality traits in table grape production, and also the important target traits in grapevine breeding. However, the information about their genetic mechanisms is limited, which results in low accuracy and efficiency of quality breeding in grapevine. Mapping and isolation of quantitative trait locus (QTLs) based on the construction of genetic linkage map is a powerful approach to decipher the genetic determinants of complex quantitative traits.

RESULTS

In the present work, a final integrated map consisting of 3411 SLAF markers on 19 linkage groups (LGs) with an average distance of 0.98 cM between adjacent markers was generated using the specific length amplified fragment sequencing (SLAF-seq) technique. A total of 9 significant stable QTLs for Muscat flavor, berry firmness and berry shape were identified on two linkage groups among the hybrids analyzed over three consecutive years from 2016 to 2018. Notably, new stable QTLs for berry firmness and berry shape were found on LG 8 respectively for the first time. Based on biological function and expression profiles of candidate genes in the major QTL regions, 3 genes (VIT_08s0007g00440, VIT_08s0040g02740 and VIT_08s0040g02350) related to berry firmness and 3 genes (VIT_08s0032g01110, VIT_08s0032g01150 and VIT_08s0105g00200) linked to berry shape were highlighted. Overexpression of VIT_08s0032g01110 in transgenic Arabidopsis plants caused the change of pod shape.

CONCLUSIONS

A new high-density genetic map with total 3411 markers was constructed with SLAF-seq technique, and thus enabled the detection of narrow interval QTLs for relevant traits in grapevine. VIT_08s0007g00440, VIT_08s0040g02740 and VIT_08s0040g02350 were found to be related to berry firmness, while VIT_08s0032g01110, VIT_08s0032g01150 and VIT_08s0105g00200 were linked to berry shape.

Construction of a high-density genetic linkage map and identification of gene controlling resistance to cucumber mosaic virus in Luffa cylindrica (L.) Roem. based on specific length amplified fragment sequencing

Luffa cylindrica L. is a cash crop which has important health, medicinal and industrial value, but no high saturation genetic map has been constructed owing to a lack of efficient markers. Furthermore, no genes were reportedly responsible for CMV resistance in Luffa spp. Specific length amplified fragment sequencing (SLAF-seq) is a valuable tool for large-scale discovery of markers and genetic mapping. The present study reported the construction of a high-density genetic map and the mapping of CMV resistant genes by using an F₂ population of 130 individuals and their two inbred line parents. A total of 271.01 Mb pair-end reads were generated. 100,077 high-quality SLAFs were detected, and 7404 of them were polymorphic. Finally, 3701 of the polymorphic markers were selected for genetic map construction, and 13 linkage groups were generated. The map spanned 1518.56 cM with an average distance of 0.41 cM between adjacent markers. Based on the newly constructed high-density map, one gene located on chromosome 1 (100.072-100.457 cM) was identified to regulate CMV resistance in L. cylindrica. A gag-polypeptide of LTR copia-type retrotransposon was predicted as the candidate gene responsible for CMV resistance in L. cylindrica. The high-density genetic map and the CMV resistant gene mapped and predicted in this study will be useful in future research.

Genome-wide assessment of population structure and genetic diversity of Chinese Lou onion using specific length amplified fragment (SLAF) sequencing

Lou onion (Allium fistulosum L. var. viviparum) is an abundant source of flavonols which provides additional health benefits to diseases. Genome-wide specific length amplified fragment (SLAF) sequencing method is a rapidly developed deep sequencing technologies used for selection and identification of genetic loci or markers. This study aimed to elucidate the genetic diversity of 122 onion accessions in China using the SLAF-seq method. A set of 122 onion accessions including 107 A.fistulosum L. var. viviparum Makino, 3 A.fistulosum L. var. gigantum Makino, 3 A.mongolicum Regel and 9 A.cepa L. accessions (3 whites, 3 reds and 3 yellows) from different regions in China were enrolled. Genomic DNA was isolated from young leaves and prepared for the SLAF-seq, which generated a total of 1,387.55 M reads and 162,321 high quality SNPs (integrity >0.5 and MAF >0.05). These SNPs were used for the construction of neighbor-joining phylogenetic tree, in which 10 A.fistulosum L. var. viviparum Makino accessions from Yinchuan (Ningxia province) and Datong (Qinghai province) had close genetic relationship. The 3 A.cepa L. clusters (red, white and yellow) had close genetic relationship especially with the 97 A.fistulosum L. var. viviparum Makino accessions. Population structure analysis suggested entire population could be clustered into 3 groups, while principal component analysis (PCA) showed there were 4 genetic groups. We confirmed the SLAF-seq approach was effective in genetic diversity analysis in red onion accessions. The key findings would provide a reference to the Lou onion germplasm in China.

Identification of candidate genes for key fibre-related QTLs and derivation of favourable alleles in Gossypium hirsutum recombinant inbred lines with G. barbadense introgressions

Fine mapping QTLs and identifying candidate genes for cotton fibre-quality and yield traits would be beneficial to cotton breeding. Here, we constructed a high-density genetic map by specific-locus amplified fragment sequencing (SLAF-seq) to identify QTLs associated with fibre-quality and yield traits using 239 recombinant inbred lines (RILs), which was developed from LMY22 (a high-yield Gossypium hirsutumL. cultivar) × LY343 (a superior fibre-quality germplasm with G. barbadenseL. introgressions). The genetic map spanned 3426.57 cM, including 3556 SLAF-based SNPs and 199 SSR marker loci. A total of 104 QTLs, including 67 QTLs for fibre quality and 37 QTLs for yield traits, were identified with phenotypic data collected from 7 environments. Among these, 66 QTLs were co-located in 19 QTL clusters on 12 chromosomes, and 24 QTLs were detected in three or more environments and determined to be stable. We also investigated the genomic components of LY343 and their contributions to fibre-related traits by deep sequencing the whole genome of LY343, and we found that genomic components from G. hirsutum races (which entered LY343 via its G. barbadense parent) contributed more favourable alleles than those from G. barbadense. We further identified six putative candidate genes for stable QTLs, including Gh_A03G1147 (GhPEL6), Gh_D07G1598 (GhCSLC6) and Gh_D13G1921 (GhTBL5) for fibre-length QTLs and Gh_D03G0919 (GhCOBL4), Gh_D09G1659 (GhMYB4) and Gh_D09G1690 (GhMYB85) for lint-percentage QTLs. Our results provide comprehensive insight into the genetic basis of the formation of fibre-related traits and would be helpful for cloning fibre-development-related genes as well as for marker-assisted genetic improvement in cotton.

Transcriptome profiling reveals key genes related to astringency during cucumber fruit development

The goal of this study was to provide quantitative data on the catechin contents and underlying molecular regulatory mechanisms in cucumber during fruit development. The dynamic changes in the total catechin contents and RNA-seq-based transcriptome profiling of the flesh and peel of the cucumber cultivar 'YanBai', which is strongly astringent, were examined at three key developmental stages 3, 6 and 9 days post-pollination. The total catechin content decreased as cucumber fruit developed and was significantly lower in the flesh than in the peel. In total, 5092 and 4004 genes were found to be differently expressed in the peel and flesh, respectively. Based on a functional annotation, eight structural genes encode enzymes involved in the catechin biosynthesis pathway. Three genes encoding 4-coumarate-CoA ligases, two genes encoding chalcone isomerases, two genes encoding dihydroflavonol-4-reductase and one gene each encoding a phenylalanine ammonia-lyase, flavanone 3-hydroxylase and cinnamate 4-hydroxylase were identified as affecting the catechin content of cucumber. The transcriptome data also revealed the significance of transcription factors, including WD40-repeat proteins, MYB and bHLH, in regulating catechin biosynthesis. These findings help increase our understanding of the molecular mechanisms controlling catechin biosynthesis and astringency development in cucumber fruit.

Identification of QTLs for Stripe Rust Resistance in a Recombinant Inbred Line Population

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most devastating fungal diseases of wheat worldwide. It is essential to discover more sources of stripe rust resistance genes for wheat breeding programs. Specific locus amplified fragment sequencing (SLAF-seq) is a powerful tool for the construction of high-density genetic maps. In this study, a set of 200 recombinant inbred lines (RILs) derived from a cross between wheat cultivars Chuanmai 42 (CH42) and Chuanmai 55 (CH55) was used to construct a high-density genetic map and to identify quantitative trait loci (QTLs) for stripe rust resistance using SLAF-seq technology. A genetic map of 2828.51 cM, including 21 linkage groups, contained 6732 single nucleotide polymorphism markers (SNP). Resistance QTLs were identified on chromosomes 1B, 2A, and 7B; Qyr.saas-7B was derived from CH42, whereas Qyr.saas-1B and Qyr.saas-2A were from CH55. The physical location of Qyr.saas-1B, which explained 6.24–34.22% of the phenotypic variation, overlapped with the resistance gene Yr29. Qyr.saas-7B accounted for up to 20.64% of the phenotypic variation. Qyr.saas-2A, a minor QTL, was found to be a likely new stripe rust resistance locus. A significant additive effect was observed when all three QTLs were combined. The combined resistance genes could be of value in breeding wheat for stripe rust resistance.

A high-density genetic map developed by specific-locus amplified fragment (SLAF) sequencing and identification of a locus controlling anthocyanin pigmentation in stalk of Zicaitai (Brassica rapa L. ssp. chinensis var. purpurea)

BACKGROUND

Caixin and Zicaitai (Brassica rapa) belong to Southern and Central China respectively. Zicaitai contains high amount of anthocyanin in leaf and stalk resulting to the purple color. Stalk is the major edible part and stalk color is an economically important trait for the two vegetables. The aim of this study is to construct a high density genetic map using the specific length amplified fragment sequencing (SLAF-seq) technique to explore genetic basis for anthocyanin pigmentation traits via quantitative trait loci (QTL) mapping.

RESULTS

We constructed a high generation linkage map with a mapping panel of F2 populations derived from 150 individuals of parental lines "Xianghongtai 01" and "Yinong 50D" with purple and green stalk respectively. The map was constructed containing 4253 loci, representing 10,940 single nucleotide polymorphism (SNP) markers spanning 1030.04 centiMorgans (cM) over 10 linkage groups (LGs), with an average distance between markers of 0.27 cM. Quantitative trait loci (QTL) analysis revealed that a major locus on chromosome 7 and 4 minor QTLs explaining 2.69-61.21% of phenotypic variation (PVE) were strongly responsible for variation in stalk color trait. Bioinformatics analysis of the major locus identified 62 protein-coding genes. Among the major locus, there were no biosynthetic genes related to anthocyanin. However, there were several transcription factors like helix-loop-helix (bHLH) bHLH, MYB in the locus. Seven predicted candidate genes were selected for the transcription level analysis. Only bHLH49 transcription factor, was significantly higher expressed in both stalks and young leaves of Xianghongtai01 than Yinong50D. An insertion and deletion (InDel) marker developed from deletion/insertion in the promoter region of bHLH49 showed significant correlation with the stalk color trait in the F2 population.

CONCLUSION

Using the constructed high-qualified linkage map, this study successfully identified QTLs for stalk color trait. The identified valuable markers and candidate genes for anthocyanin accumulation in stalk will provide useful information for molecular regulation of anthocyanin biosynthesis. Overall our findings will lay a foundation for functional gene cloning, marker-assisted selection (MAS) and molecular breeding of important economic traits in B. rapa.

Molecular Analysis of Evolution and Origins of Cultivated Hawthorn (Crataegus spp.) and Related Species in China

Hawthorn is of high economic value owing to its medicinal properties and health benefits. Crataegus is a member of the Rosaceae family; the genus has a complicated taxonomic history, and several theories on its origin have been proposed. In this study, 53 accessions from seven Crataegus taxa native to China and accessions of exotic Crataegus species (two from Europe and one from North America) were analyzed by specific locus amplified fragment sequencing (SLAF-seq). In total, 933,450 single-nucleotide polymorphisms were identified after filtering and used to investigate the species' genomic evolution. Phylogenetic trees derived from nuclear simple sequence repeats (SSRs) and SLAF-seq data showed the same topology, in which Crataegus maximowiczii and Crataegus sanguineae formed a closely related cluster that was clearly separated from the cluster composed of Crataegus hupehensis, Crataegus pinnatifida, Crataegus pinnatifida var. major, Crataegus bretschneideri and Crataegus scabrifolia. Phylogenetic and structure analysis indicated that the seven Chinese Crataegus taxa had two separate speciation events. Plants that evolved the southwestern route shared the genepool with the European species, whereas plants along the northeastern route shared the genepool with the North American species. TreeMix genetic analysis revealed that C. bretschneideri may have a hybrid origin. This study provides valuable information on the origins of Chinese Crataegus and suggests an evolutionary model for the main Crataegus species that native to China.

High-Density Genetic Linkage Maps Provide Novel Insights Into ZW/ZZ Sex Determination System and Growth Performance in Mud Crab (Scylla paramamosain)

Mud crab, Scylla paramamosain is one of the most important crustacean species in global aquaculture. To determine the genetic basis of sex and growth-related traits in S. paramamosain, a high-density genetic linkage map with 16,701 single nucleotide polymorphisms (SNPs) was constructed using SLAF-seq and a full-sib family. The consensus map has 49 linkage groups, spanning 5,996.66 cM with an average marker-interval of 0.81 cM. A total of 516 SNP markers, including 8 female-specific SNPs segregated in two quantitative trait loci (QTLs) for phenotypic sex were located on LG32. The presence of female-specific SNP markers only on female linkage map, their segregation patterns and lower female: male recombination rate strongly suggest the conformation of a ZW/ZZ sex determination system in S. paramamosain. The QTLs of most (90%) growth-related traits were found within a small interval (25.18–33.74 cM) on LG46, highlighting the potential involvement of LG46 in growth. Four markers on LG46 were significantly associated with 10–16 growth-related traits. BW was only associated with marker 3846. Based on the annotation of transcriptome data, 11 and 2 candidate genes were identified within the QTL regions of sex and growth-related traits, respectively. The newly constructed high-density genetic linkage map with sex-specific SNPs, and the identified QTLs of sex- and growth-related traits serve as a valuable genetic resource and solid foundation for marker-assisted selection and genetic improvement of crustaceans.

Construction of a High-Density Genetic Map and Analysis of Seed-Related Traits Using Specific Length Amplified Fragment Sequencing for Cucurbita maxima

Seed traits are agronomically important for Cucurbita breeding, but the genes controlling seed size, seed weight and seed number have not been mapped in Cucurbita maxima (C. maxima). In this study, 100 F₂ individual derived from two parental lines, "2013-12" and "9-6", were applied to construct a 3,376.87-cM genetic map containing 20 linkage groups (LGs) with an average genetic distance of 0.47 cM using a total of 8,406 specific length amplified fragment (SLAF) markers in C. maxima. Ten quantitative trait loci (QTLs) of seed width (SW), seed length (SL) and hundred-seed weight (HSW) were identified using the composite interval mapping (CIM) method. The QTLs affecting SW, SL and HSW explained a maximum of 38.6%, 28.9% and 17.2% of the phenotypic variation and were detected in LG6, LG6 and LG17, respectively. To validate these results, an additional 150 F₂ individuals were used for QTL mapping of SW and SL with cleaved amplified polymorphic sequence (CAPS) markers. We found that two major QTLs, SL6-1 and SW6-1, could be detected in both SLAF-seq and CAPS markers in an overlapped region. Based on gene annotation and non-synonymous single-nucleotide polymorphisms (SNPs) in the major SWand SL-associated regions, we found that two genes encoding a VQ motif and an E3 ubiquitin-protein ligase may be candidate genes influencing SL, while an F-box and leucinerich repeat (LRR) domain-containing protein is the potential regulator for SW in C. maxima. This study provides the first high-density linkage map of C. maxima using SNPs developed by SLAF-seq technology, which is a powerful tool for associated mapping of important agronomic traits, map-based gene cloning and marker-assisted selection (MAS)-based breeding in C. maxima.

Construction of High-Density Genetic Map and Mapping Quantitative Trait Loci for Growth Habit-Related Traits of Peanut (Arachis hypogaea L.)

Plant growth habit is an important and complex agronomic trait and is associated with yield, disease resistance, and mechanized harvesting in peanuts. There are at least two distinct growth habits (erect and prostrate) and several intermediate forms existing in the peanut germplasm. A recombinant inbred line population containing 188 individuals was developed from a cross of "Jihua 5" and "M130" for genetically dissecting the architecture of the growth habit. A new high-density genetic linkage map was constructed by using specific locus amplified fragment sequencing technology. The map contains 2,808 single-nucleotide polymorphism markers distributed on 20 linkage groups with a total length of 1,308.20 cM and an average inter-marker distance of 0.47 cM. The quantitative trait locus (QTL) analysis of the growth habit-related traits was conducted based on phenotyping data from seven environments. A total of 39 QTLs for growth habit-related traits was detected on 10 chromosomes explaining 4.55-27.74% of the phenotypic variance, in which 6 QTLs were for lateral branch angle, 8 QTLs were for extent radius, 7 QTLs were for the index of plant type, 11 QTLs were for main stem height, and 7 QTLs were for lateral branch length. Among these QTLs, 12 were co-localized on chromosome B05 spanning an approximately 0.17 Mb physical interval in comparison with the allotetraploid reference genome of "Tifrunner." Analysis of the co-localized genome region has shown that the putative genes are involved in light and hormones and will facilitate peanut growth habit molecular breeding and study of peanut domestication.

Construction of a high-density genetic linkage map and QTL mapping of oleic acid content and three agronomic traits in sunflower (Helianthus annuus L.) using specific-locus amplified fragment sequencing (SLAF-seq)

High-density genetic linkage maps are particularly important for quantitative trait loci (QTL) mapping, genome assembly, and marker-assisted selection (MAS) in plants. In this study, a high-density genetic linkage map of sunflower (Helianthus annuus L.) was constructed using an F₂ population generated from a cross between Helianthus annuus L. '86-1' and 'L-1-OL-1' via specific-locus amplified fragment sequencing (SLAF-seq). After sequence preprocessing, 530.50 M reads (105.60 Gb) were obtained that contained a total of 343,197 SLAFs, of which 39,589 were polymorphic. Of the polymorphic SLAFs, 6,136 were organized into a linkage map consisting of 17 linkage groups (LGs) spanning 2,221.86 cM, with an average genetic distance of 0.36 cM between SLAFs. Based on this high-density genetic map, QTL analysis was performed that focused on four sunflower phenotypic traits: oleic acid content (OAC), plant height (PH), head diameter (HD), and stem diameter (SD). Subsequently, for these four traits eight QTLs were detected that will likely be useful for increasing our understanding of genetic factors underlying these traits and for use in marker-assisted selection (MAS) for future sunflower breeding.

Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached "suggestive" genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.

A high-density genetic map of extra-long staple cotton (Gossypium barbadense) constructed using genotyping-by-sequencing based single nucleotide polymorphic markers and identification of fiber traits-related QTL in a recombinant inbred line population

Background

Gossypium barbadense (Sea Island, Egyptian or Pima cotton) cotton has high fiber quality, however, few studies have investigated the genetic basis of its traits using molecular markers. Genome complexity reduction approaches such as genotyping-by-sequencing have been utilized to develop abundant markers for the construction of high-density genetic maps to locate quantitative trait loci (QTLs).

Results

The Chinese G. barbadense cultivar 5917 and American Pima S-7 were used to develop a recombinant inbred line (RIL) population with 143 lines. The 143 RILs together with their parents were tested in three replicated field tests for lint yield traits (boll weight and lint percentage) and fiber quality traits (fiber length, fiber elongation, fiber strength, fiber uniformity and micronaire) and then genotyped using GBS to develop single-nucleotide polymorphism (SNP) markers. A high-density genetic map with 26 linkage groups (LGs) was constructed using 3557 GBS SNPs spanning a total genetic distance of 3076.23 cM at an average density of 1.09 cM between adjacent markers. A total of 42 QTLs were identified, including 24 QTLs on 12 LGs for fiber quality and 18 QTLs on 7 LGs for lint yield traits, with LG1 (9 QTLs), LG10 (7 QTLs) and LG14 (6 QTLs) carrying more QTLs. Common QTLs for the same traits and overlapping QTLs for different traits were detected. Each individual QTLs explained 0.97 to 20.7% of the phenotypic variation.

Conclusions

This study represents one of the first genetic mapping studies on the fiber quality and lint yield traits in a RIL population of G. barbadense using GBS-SNPs. The results provide important information for the subsequent fine mapping of QTLs and the prediction of candidate genes towards map-based cloning and marker-assisted selection in cotton.

High-Density Genetic Map Construction and Identification of QTLs Controlling Oleic and Linoleic Acid in Peanut using SLAF-seq and SSRs

The cultivated peanut, A. hypogaea L., is an important oil and food crop globally.High-density genetic linkage mapping is a valuable and effective method for exploring complex quantitative traits. In this context, a recombinant inbred line (RIL) of 146 lines was developed by crossing Huayu28 and P76. We developed 433,679 high-quality SLAFs, of which 29,075 were polymorphic. 4,817 SLAFs were encoded and grouped into different segregation patterns. A high-resolution genetic map containing 2,334 markers (68 SSRs and 2,266 SNPs) on 20 linkage groups (LGs) spanning 2586.37 cM was constructed for peanut. The average distance between adjacent markers was 2.25 cM. Based on phenotyping in seven environments, QTLs for oleic acid (C18:1), linoleic acid (C18:2) and the ratio of oleic acid to linoleic acid (O/L) were identified and positioned on linkage groups A03, A04, A09, B09 and B10. Marker2575339 and Marker2379598 in B09 were associated with C18:1, C18:2 and O/L in seven environments, Marker4391589 and Marker4463600 in A09 were associated with C18:1, C18:2 and O/L in six environments. This map exhibits high resolution and accuracy, which will facilitate QTL discovery for essential agronomic traits in peanut.

QTL mapping of downy and powdery mildew resistances in PI 197088 cucumber with genotyping-by-sequencing in RIL population

Host resistances in PI 197088 cucumber to downy and powdery mildew pathogens are conferred by 11 (3 with major effect) and 4 (1 major effect) QTL, respectively, and three of which are co-localized. The downy mildew (DM) and powdery mildew (PM) are the two most important foliar diseases of cucurbit crops worldwide. The cucumber accession PI 197088 exhibits high-level resistances to both pathogens. Here, we reported QTL mapping results for DM and PM resistances with 148 recombinant inbred lines from a cross between PI 197088 and the susceptible line 'Coolgreen'. Phenotypic data on responses to natural DM and PM infection were collected in multi-year and multi-location replicated field trials. A high-density genetic map with 2780 single nucleotide polymorphisms (SNPs) from genotyping-by-sequencing and 55 microsatellite markers was developed, which revealed genomic regions with segregation distortion and mis-assemblies in the '9930' cucumber draft genome. QTL analysis identified 11 and 4 QTL for DM and PM resistances accounting for more than 73.5 and 63.0% total phenotypic variance, respectively. Among the 11 DM resistance QTL, dm5.1, dm5.2, and dm5.3 were major-effect contributing QTL, whereas dm1.1, dm2.1, and dm6.2 conferred susceptibility. Of the 4 QTL for PM resistance, pm5.1 was the major-effect QTL explaining 32.4% phenotypic variance and the minor-effect QTL pm6.1 contributed to disease susceptibility. Three PM QTL, pm2.1, pm5.1, and pm6.1, were co-localized with DM QTL dm2.1, dm5.2, and dm6.1, respectively, which was consistent with the observed linkage of PM and DM resistances in PI 197088. The genetic architecture of DM resistance in PI 197088 and another resistant line WI7120 (PI 330628) was compared, and the potential of using PI 197088 in cucumber breeding for downy and powdery mildew resistances is discussed.

Construction of a high-density genetic map and identification of QTLs for cucumber mosaic virus resistance in pepper (Capsicum annuum L.) using specific length amplified fragment sequencing (SLAF-seq)

Pepper (Capsicum) is one of the most important vegetable and spice crops. Aphid-transmitted cucumber mosaic virus (CMV) causes significant damage to pepper crops across the world. The genetic basis of CMV resistance in pepper is complex, and the mechanisms underlying resistance remain largely unknown. Here, we employed a SLAF-seq approach to generate a high-density genetic map of pepper. The map spanned 1,785.46 cM, containing 12,727 markers on 12 chromosomes, with a mean marker distance of 0.16 cM between adjacent markers. We used this map and the interval mapping (IM) and multiple QTL mapping (MQM) procedures to detect genetic regions associated with quantitative trait for CMV resistance. Three QTLs, qcmv11.1, qcmv11.2 and qcmv12.1, conferred resistance to CMV and showed trait variation of 10.2%, 19.2% and 7.3% respectively. Our results will help to develop markers linked to CMV-resistant QTLs to improve pepper resistance to CMV.

QTL Mapping of Fiber-Related Traits Based on a High-Density Genetic Map in Flax (Linum usitatissimum L.)

A genetic map is an important and valuable tool for quantitative trait locus (QTL) mapping, marker-assisted selection (MAS)-based breeding, and reference-assisted chromosome assembly. In this study, 112 F₂ plants from a cross between Linum usitatissimum L. "DIANE" and "NY17" and parent plants were subjected to high-throughput sequencing and specific-locus amplified fragment (SLAF) library construction. After preprocessing, 61.64 Gb of raw data containing 253.71 Mb paired-end reads, each 101 bp in length, were obtained. A total of 192,797 SLAFs were identified, of which 23,115 were polymorphic, with a polymorphism rate of 11.99%. Finally, 2,339 SLAFs were organized into a linkage map consisting of 15 linkage groups (LGs). The total length of the genetic map was 1483.25 centimorgans (cM) and the average distance between adjacent markers was 0.63 cM. Combined with flax chromosome-scale pseudomolecules, 12 QTLs associating with 6 flax fiber-related traits were mapped on the chromosomal scaffolds. This high-density genetic map of flax should serve as a foundation for flax fine QTL mapping, draft genome assembly, and MAS-guided breeding. Ultimately, the genomic regions identified in this research could potentially be valuable for improving flax fiber cultivars, as well as for identification of candidate genes involved in flax fiber formation processes.

SIGNIFICANCE STATEMENT

A high-density genetic map of flax was constructed, and QTLs were identified on the sequence scaffolds to be interrelated with fiber-related traits. The results of this study will not only provide a platform for gene/QTL fine mapping, map-based gene isolation, and molecular breeding for flax, but also provide a reference to help position sequence scaffolds on the physical map and assist in the process of assembling the flax genome sequence.

Construction of an SNP-based high-density linkage map for flax (Linum usitatissimum L.) using specific length amplified fragment sequencing (SLAF-seq) technology

Flax is an important crop for oil and fiber, however, no high-density genetic maps have been reported for this species. Specific length amplified fragment sequencing (SLAF-seq) is a high-resolution strategy for large scale de novo discovery and genotyping of single nucleotide polymorphisms. In this study, SLAF-seq was employed to develop SNP markers in an F2 population to construct a high-density genetic map for flax. In total, 196.29 million paired-end reads were obtained. The average sequencing depth was 25.08 in male parent, 32.17 in the female parent, and 9.64 in each F2 progeny. In total, 389,288 polymorphic SLAFs were detected, from which 260,380 polymorphic SNPs were developed. After filtering, 4,638 SNPs were found suitable for genetic map construction. The final genetic map included 4,145 SNP markers on 15 linkage groups and was 2,632.94 cM in length, with an average distance of 0.64 cM between adjacent markers. To our knowledge, this map is the densest SNP-based genetic map for flax. The SNP markers and genetic map reported in here will serve as a foundation for the fine mapping of quantitative trait loci (QTLs), map-based gene cloning and marker assisted selection (MAS) for flax.

Construction of the First High-Density Genetic Linkage Map and Analysis of Quantitative Trait Loci for Growth-Related Traits in Sinonovacula constricta

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.

Construction of a high-density genetic map for grape using specific length amplified fragment (SLAF) sequencing

Genetic maps are important tools in plant genomics and breeding. We report a large-scale discovery of single nucleotide polymorphisms (SNPs) using the specific length amplified fragment sequencing (SLAF-seq) technique for the construction of high-density genetic maps for two elite wine grape cultivars, ‘Chardonnay’ and ‘Beibinghong’, and their 130 F₁ plants. A total of 372.53 M paired-end reads were obtained after preprocessing. The average sequencing depth was 33.81 for ‘Chardonnay’ (the female parent), 48.20 for ‘Beibinghong’ (the male parent), and 12.66 for the F₁ offspring. We detected 202,349 high-quality SLAFs of which 144,972 were polymorphic; 10,042 SNPs were used to construct a genetic map that spanned 1,969.95 cM, with an average genetic distance of 0.23 cM between adjacent markers. This genetic map contains the largest molecular marker number of the grape maps so far reported. We thus demonstrate that SLAF-seq is a promising strategy for the construction of high-density genetic maps; the map that we report here is a good potential resource for QTL mapping of genes linked to major economic and agronomic traits, map-based cloning, and marker-assisted selection of grape.

Genome-wide association study of rust traits in orchardgrass using SLAF-seq technology

BACKGROUND

While orchardgrass (Dactylis glomerata L.) is a well-known perennial forage species, rust diseases cause serious reductions in the yield and quality of orchardgrass; however, genetic mechanisms of rust resistance are not well understood in orchardgrass.

RESULTS

In this study, a genome-wide association study (GWAS) was performed using specific-locus amplified fragment sequencing (SLAF-seq) technology in orchardgrass. A total of 2,334,889 SLAF tags were generated to produce 2,309,777 SNPs. ADMIXTURE analysis revealed unstructured subpopulations for 33 accessions, indicating that this orchardgrass population could be used for association analysis. Linkage disequilibrium (LD) analysis revealed an average r2 of 0.4 across all SNP pairs, indicating a high extent of LD in these samples. Through GWAS, a total of 4,604 SNPs were found to be significantly (P < 0.01) associated with the rust trait. The bulk analysis discovered a number of 5,211 SNPs related to rust trait. Two candidate genes, including cytochrome P450, and prolamin were implicated in disease resistance through prediction of functional genes surrounding each high-quality SNP (P < 0.01) associated with rust traits based on GWAS analysis and bulk analysis.

CONCLUSIONS

The large number of SNPs associated with rust traits and these two candidate genes may provide the basis for further research on rust resistance mechanisms and marker-assisted selection (MAS) for rust-resistant lineages.

Genome-wide assessment of population structure and genetic diversity and development of a core germplasm set for sweet potato based on specific length amplified fragment (SLAF) sequencing

Sweet potato, Ipomoea batatas (L.) Lam., is an important food crop that is cultivated worldwide. However, no genome-wide assessment of the genetic diversity of sweet potato has been reported to date. In the present study, the population structure and genetic diversity of 197 sweet potato accessions most of which were from China were assessed using 62,363 SNPs. A model-based structure analysis divided the accessions into three groups: group 1, group 2 and group 3. The genetic relationships among the accessions were evaluated using a phylogenetic tree, which clustered all the accessions into three major groups. A principal component analysis (PCA) showed that the accessions were distributed according to their population structure. The mean genetic distance among accessions ranged from 0.290 for group 1 to 0.311 for group 3, and the mean polymorphic information content (PIC) ranged from 0.232 for group 1 to 0.251 for group 3. The mean minor allele frequency (MAF) ranged from 0.207 for group 1 to 0.222 for group 3. Analysis of molecular variance (AMOVA) showed that the maximum diversity was within accessions (89.569%). Using CoreHunter software, a core set of 39 accessions was obtained, which accounted for approximately 19.8% of the total collection. The core germplasm set of sweet potato developed will be a valuable resource for future sweet potato improvement strategies.

High-Density Genetic Map Construction and Gene Mapping of Basal Branching Habit and Flowers per Leaf Axil in Sesame

A good genetic map can provide the framework for quantitative trait loci (QTL) analysis, map-based gene cloning, and genome sequence assembling. The main objectives of this study were to develop a high-density genetic linkage map using specific length amplified fragment sequencing (SLAF-seq) in sesame. In the result, a high-resolution genetic map with 9,378 SLAF markers and 13 linkage groups (LGs) was constructed. The map spanned a total genetic distance of 1,974.23 cM, and the mean LG length was 151.86 cM, with an average genetic distance of 0.22 cM between adjacent markers. Based on the newly constructed genetic map, genes for basal branching habit (SiBH) and flowers per leaf axil (SiFA) were mapped to LG5 and LG11, respectively.

Genome-Wide Association Study Identifying Candidate Genes Influencing Important Agronomic Traits of Flax (Linum usitatissimum L.) Using SLAF-seq

Flax (Linum usitatissimum L.) is an important cash crop, and its agronomic traits directly affect yield and quality. Molecular studies on flax remain inadequate because relatively few flax genes have been associated with agronomic traits or have been identified as having potential applications. To identify markers and candidate genes that can potentially be used for genetic improvement of crucial agronomic traits, we examined 224 specimens of core flax germplasm; specifically, phenotypic data for key traits, including plant height, technical length, number of branches, number of fruits, and 1000-grain weight were investigated under three environmental conditions before specific-locus amplified fragment sequencing (SLAF-seq) was employed to perform a genome-wide association study (GWAS) for these five agronomic traits. Subsequently, the results were used to screen single nucleotide polymorphism (SNP) loci and candidate genes that exhibited a significant correlation with the important agronomic traits. Our analyses identified a total of 42 SNP loci that showed significant correlations with the five important agronomic flax traits. Next, candidate genes were screened in the 10 kb zone of each of the 42 SNP loci. These SNP loci were then analyzed by a more stringent screening via co-identification using both a general linear model (GLM) and a mixed linear model (MLM) as well as co-occurrences in at least two of the three environments, whereby 15 final candidate genes were obtained. Based on these results, we determined that UGT and PL are candidate genes for plant height, GRAS and XTH are candidate genes for the number of branches, Contig1437 and LU0019C12 are candidate genes for the number of fruits, and PHO1 is a candidate gene for the 1000-seed weight. We propose that the identified SNP loci and corresponding candidate genes might serve as a biological basis for improving crucial agronomic flax traits.

Construction of a high-density genetic map and QTL mapping for pearl quality-related traits in Hyriopsis cumingii

A high-density genetic map is essential for quantitative trait locus (QTL) fine mapping. In this study, 4,508 effective single nucleotide polymorphism markers (detected using specific-locus amplified fragment sequencing) and 475 microsatellites were mapped to 19 linkage groups (LGs) using a family with 157 individuals. The map spanned 2,713 cM, with an average of 259 markers and 79 loci per LG and an average inter-marker distance of 1.81 cM. To identify QTLs for pearl quality traits, 26 putatively significant QTLs were detected for 10 traits, including, three for shell width, seven for body weight, two for shell weight, two for margin mantle weight, five for inner mantle weight, and seven for shell nacre colour. Among them, five QTLs associated with shell nacre colour were mapped to LG17 and explained 19.7% to 22.8% of the trait variation; this suggests that some important genes or loci determine shell nacre colour in LG17. The linkage map and mapped QTLs for shell nacre colour would be useful for improving the quality of Hyriopsis cumingii via marker-assisted selection.

Identification and Validation of SNP Markers Linked to Dwarf Traits Using SLAF-Seq Technology in Lagerstroemia

The genetic control of plant architecture is a promising approach to breed desirable cultivars, particularly in ornamental flowers. In this study, the F1 population (142 seedlings) derived from Lagerstroemia fauriei (non-dwarf) × L. indica 'Pocomoke' (dwarf) was phenotyped for six traits (plant height (PH), internode length (IL), internode number, primary lateral branch height (PLBH), secondary lateral branch height and primary branch number), and the IL and PLBH traits were positively correlated with the PH trait and considered representative indexes of PH. Fifty non-dwarf and dwarf seedlings were pooled and subjected to a specific-locus amplified fragment sequencing (SLAF-seq) method, which screened 1221 polymorphic markers. A total of 3 markers segregating between bulks were validated in the F1 population, with the M16337 and M38412 markers highly correlated with the IL trait and the M25207 marker highly correlated with the PLBH trait. These markers provide a predictability of approximately 80% using a single marker (M25207) and a predictability of 90% using marker combinations (M16337 + M25207) in the F1 population, which revealed that the IL and the PLBH traits, especially the PLBH, were the decisive elements for PH in terms of molecular regulation. Further validation was performed in the BC1 population and a set of 28 Lagerstroemia stocks using allele-specific PCR (AS-PCR) technology, and the results showed the stability and reliability of the SNP markers and the co-determination of PH by multiple genes. Our findings provide an important theoretical and practical basis for the early prediction and indirect selection of PH using the IL and the PLBH, and the detected SNPs may be useful for marker-assisted selection (MAS) in crape myrtle.

Construction of high-density genetic linkage map and identification of flowering-time QTLs in orchardgrass using SSRs and SLAF-seq

Orchardgrass (Dactylis glomerata L.) is one of the most economically important perennial, cool-season forage species grown and pastured worldwide. High-density genetic linkage mapping is a valuable and effective method for exploring complex quantitative traits. In this study, we developed 447,177 markers based on SLAF-seq and used them to perform a comparative genomics analysis. Perennial ryegrass sequences were the most similar (5.02%) to orchardgrass sequences. A high-density linkage map of orchardgrass was constructed using 2,467 SLAF markers and 43 SSRs, which were distributed on seven linkage groups spanning 715.77 cM. The average distance between adjacent markers was 0.37 cM. Based on phenotyping in four environments, 11 potentially significant quantitative trait loci (QTLs) for two target traits–heading date (HD) and flowering time (FT)–were identified and positioned on linkage groups LG1, LG3, and LG5. Significant QTLs explained 8.20–27.00% of the total phenotypic variation, with the LOD ranging from 3.85–12.21. Marker167780 and Marker139469 were associated with FT and HD at the same location (Ya’an) over two different years. The utility of SLAF markers for rapid generation of genetic maps and QTL analysis has been demonstrated for heading date and flowering time in a global forage grass.

A High-Density Genetic Map of Tetraploid Salix matsudana Using Specific Length Amplified Fragment Sequencing (SLAF-seq)

As a salt-tolerant arbor tree species, Salix matsudana plays an important role in afforestation and greening in the coastal areas of China. To select superior Salix varieties that adapt to wide saline areas, it is of paramount importance to understand and identify the mechanisms of salt-tolerance at the level of the whole genome. Here, we describe a high-density genetic linkage map of S. matsudana that represents a good coverage of the Salix genome. An intraspecific F₁ hybrid population was established by crossing the salt-sensitive “Yanjiang” variety as the female parent with the salt-tolerant “9901” variety as the male parent. This population, along with its parents, was genotyped by specific length amplified fragment sequencing (SLAF-seq), leading to 277,333 high-quality SLAF markers. By marker analysis, we found that both the parents and offspring were tetraploid. The mean sequencing depth was 53.20-fold for “Yanjiang”, 47.41-fold for “9901”, and 11.02-fold for the offspring. Of the SLAF markers detected, 42,321 are polymorphic with sufficient quality for map construction. The final genetic map was constructed using 6,737 SLAF markers, covering 38 linkage groups (LGs). The genetic map spanned 5,497.45 cM in length, with an average distance of 0.82 cM. As a first high-density genetic map of S. matsudana constructed from salt tolerance-varying varieties, this study will provide a foundation for mapping quantitative trait loci that modulate salt tolerance and resistance in Salix and provide important references for molecular breeding of this important forest tree.

Construction of a high-density genetic map by specific locus amplified fragment sequencing (SLAF-seq) and its application to Quantitative Trait Loci (QTL) analysis for boll weight in upland cotton (Gossypium hirsutum.)

BACKGROUND

Upland Cotton (Gossypium hirsutum) is one of the most important worldwide crops it provides natural high-quality fiber for the industrial production and everyday use. Next-generation sequencing is a powerful method to identify single nucleotide polymorphism markers on a large scale for the construction of a high-density genetic map for quantitative trait loci mapping.

RESULTS

In this research, a recombinant inbred lines population developed from two upland cotton cultivars 0-153 and sGK9708 was used to construct a high-density genetic map through the specific locus amplified fragment sequencing method. The high-density genetic map harbored 5521 single nucleotide polymorphism markers which covered a total distance of 3259.37 cM with an average marker interval of 0.78 cM without gaps larger than 10 cM. In total 18 quantitative trait loci of boll weight were identified as stable quantitative trait loci and were detected in at least three out of 11 environments and explained 4.15-16.70 % of the observed phenotypic variation. In total, 344 candidate genes were identified within the confidence intervals of these stable quantitative trait loci based on the cotton genome sequence. These genes were categorized based on their function through gene ontology analysis, Kyoto Encyclopedia of Genes and Genomes analysis and eukaryotic orthologous groups analysis.

CONCLUSIONS

This research reported the first high-density genetic map for Upland Cotton (Gossypium hirsutum) with a recombinant inbred line population using single nucleotide polymorphism markers developed by specific locus amplified fragment sequencing. We also identified quantitative trait loci of boll weight across 11 environments and identified candidate genes within the quantitative trait loci confidence intervals. The results of this research would provide useful information for the next-step work including fine mapping, gene functional analysis, pyramiding breeding of functional genes as well as marker-assisted selection.

Genome-Wide Single-Nucleotide Polymorphisms Discovery and High-Density Genetic Map Construction in Cauliflower Using Specific-Locus Amplified Fragment Sequencing

Molecular markers and genetic maps play an important role in plant genomics and breeding studies. Cauliflower is an important and distinctive vegetable; however, very few molecular resources have been reported for this species. In this study, a novel, specific-locus amplified fragment (SLAF) sequencing strategy was employed for large-scale single nucleotide polymorphism (SNP) discovery and high-density genetic map construction in a double-haploid, segregating population of cauliflower. A total of 12.47 Gb raw data containing 77.92 M pair-end reads were obtained after processing and 6815 polymorphic SLAFs between the two parents were detected. The average sequencing depths reached 52.66-fold for the female parent and 49.35-fold for the male parent. Subsequently, these polymorphic SLAFs were used to genotype the population and further filtered based on several criteria to construct a genetic linkage map of cauliflower. Finally, 1776 high-quality SLAF markers, including 2741 SNPs, constituted the linkage map with average data integrity of 95.68%. The final map spanned a total genetic length of 890.01 cM with an average marker interval of 0.50 cM, and covered 364.9 Mb of the reference genome. The markers and genetic map developed in this study could provide an important foundation not only for comparative genomics studies within Brassica oleracea species but also for quantitative trait loci identification and molecular breeding of cauliflower.

Mapping of a Novel Race Specific Resistance Gene to Phytophthora Root Rot of Pepper (Capsicum annuum) Using Bulked Segregant Analysis Combined with Specific Length Amplified Fragment Sequencing Strategy

Phytophthora root rot caused by Phytophthora capsici (P. capsici) is a serious limitation to pepper production in Southern China, with high temperature and humidity. Mapping PRR resistance genes can provide linked DNA markers for breeding PRR resistant varieties by molecular marker-assisted selection (MAS). Two BC1 populations and an F2 population derived from a cross between P. capsici-resistant accession, Criollo de Morelos 334 (CM334) and P. capsici-susceptible accession, New Mexico Capsicum Accession 10399 (NMCA10399) were used to investigate the genetic characteristics of PRR resistance. PRR resistance to isolate Byl4 (race 3) was controlled by a single dominant gene, PhR10, that was mapped to an interval of 16.39Mb at the end of the long arm of chromosome 10. Integration of bulked segregant analysis (BSA) and Specific Length Amplified Fragment sequencing (SLAF-seq) provided an efficient genetic mapping strategy. Ten polymorphic Simple Sequence Repeat (SSR) markers were found within this region and used to screen the genotypes of 636 BC1 plants, delimiting PhR10 to a 2.57 Mb interval between markers P52-11-21 (1.5 cM away) and P52-11-41 (1.1 cM). A total of 163 genes were annotated within this region and 31 were predicted to be associated with disease resistance. PhR10 is a novel race specific gene for PRR, and this paper describes linked SSR markers suitable for marker-assisted selection of PRR resistant varieties, also laying a foundation for cloning the resistance gene.

High-Density Genetic Mapping Identifies New Major Loci for Tolerance to Low-Phosphorus Stress in Soybean

Phosphorus (P) is essential for all living cells and organisms, and low-P stress represents a major constraint on plant growth and yield worldwide. Soybean is an important economical resource of protein and oil for human and animals, and soybean is also a high-P demand species that is sensitive to low-P stress, which is considered a major constraint on soybean production. However, P efficiency is an important complex quantitative trait involving multiple genes, and the mechanisms underlying soybean P efficiency are largely unknown. Here, we reported the construction of a high-density genetic map using a specific-locus amplified fragment sequencing (SLAF-seq) strategy in soybean. This map, spanning 3020.59 cM in length, contained 6159 markers on 20 chromosomes, with an average distance of 0.49 cM between adjacent markers. Based on this map, 20 loci, including eight novel loci, associated with P efficiency-related traits were identified across multiple years and treatments. The confidence intervals of almost all QTLs were refined significantly, and the accuracy of this map was evidenced by coincident detections of the previously identified P efficiency-related genes GmACP1 and GmPT1. Notably, a highly significant novel QTL located on chromosome 4, q4-2, was identified across traits, years and treatments. Several candidate genes, such as a pectin methylesterase-encoding gene (Glyma.04G214000) and a protein kinase gene (Glyma.13G161900), with significantly differential expression upon low-P stress were considered as promising candidates involved in regulating soybean P efficiency. Markers that tightly associated with P efficiency could be used for marker-assisted selection in a soybean P efficient breeding program. Further, dissection of these QTLs will facilitate gene cloning underlying P efficiency in soybean, and increase our understanding of efficient use of P in enhancing crop yield.

Cloning and Characterization of TaTGW-7A Gene Associated with Grain Weight in Wheat via SLAF-seq-BSA

Thousand-grain weight (TGW) of wheat (Triticum aestivum L.) contributes significantly to grain yield. In the present study, a candidate gene associated with TGW was identified through specific-locus amplified fragment sequencing (SLAF-seq) of DNA bulks of recombinant inbred lines (RIL) derived from the cross between Jing 411 and Hongmangchun 21. The gene was located on chromosome 7A, designated as TaTGW-7A with a complete genome sequence and an open reading frame (ORF). A single nucleotide polymorphism (SNP) was present in the first exon between two alleles at TaTGW-7A locus, resulting in a Val to Ala substitution, corresponding to a change from higher to lower TGW. Cleaved amplified polymorphic sequence (CAPS) (TGW7A) and InDel (TG9) markers were developed to discriminate the two alleles TaTGW-7Aa and TaTGW-7Ab for higher and lower TGW, respectively. A major QTL co-segregating with TaTGW-7A explained 21.7-27.1% of phenotypic variance for TGW in the RIL population across five environments. The association of TaTGW-7A with TGW was further validated in a natural population and Chinese mini-core collections. Quantitative real-time PCR revealed higher transcript levels of TaTGW-7Aa than those of TaTGW-7Ab during grain development. High frequencies of the superior allele TaTGW-7Aa for higher TGW in Chinese mini-core collections (65.0%) and 501 wheat varieties (86.0%) indicated a strong and positive selection of this allele in wheat breeding. The molecular markers TGW7A and TG9 can be used for improvement of TGW in breeding programs.

Construction of a High-Density Genetic Map Based on Large-Scale Marker Development in Mango Using Specific-Locus Amplified Fragment Sequencing (SLAF-seq)

Genetic maps are particularly important and valuable tools for quantitative trait locus (QTL) mapping and marker assisted selection (MAS) of plant with desirable traits. In this study, 173 F1 plants from a cross between Mangifera indica L. "Jin-Hwang" and M. indica L. "Irwin" and their parent plants were subjected to high-throughput sequencing and specific-locus amplified fragment (SLAF) library construction. After preprocessing, 66.02 Gb of raw data containing 330.64 M reads were obtained. A total of 318,414 SLAFs were detected, of which 156,368 were polymorphic. Finally, 6594 SLAFs were organized into a linkage map consisting of 20 linkage groups (LGs). The total length of the map was 3148.28 cM and the average distance between adjacent markers was 0.48 cM. This map could be considered, to our knowledge, the first high-density genetic map of mango, and might form the basis for fine QTL mapping and MAS of mango.

A High-Density Genetic Linkage Map for Cucumber (Cucumis sativus L.): Based on Specific Length Amplified Fragment (SLAF) Sequencing and QTL Analysis of Fruit Traits in Cucumber

High-density genetic linkage map plays an important role in genome assembly and quantitative trait loci (QTL) fine mapping. Since the coming of next-generation sequencing, makes the structure of high-density linkage maps much more convenient and practical, which simplifies SNP discovery and high-throughput genotyping. In this research, a high-density linkage map of cucumber was structured using specific length amplified fragment sequencing, using 153 F2 populations of S1000 × S1002. The high-density genetic map composed 3,057 SLAFs, including 4,475 SNP markers on seven chromosomes, and spanned 1061.19 cM. The average genetic distance is 0.35 cM. Based on this high-density genome map, QTL analysis was performed on two cucumber fruit traits, fruit length and fruit diameter. There are 15 QTLs for the two fruit traits were detected.

Next generation sequencing and omics in cucumber (Cucumis sativus L.) breeding directed research

In the post-genomic era the availability of genomic tools and resources is leading us to novel generation methods in plant breeding, as they facilitate the study of the genotype and its relationship with the phenotype, in particular for complex traits. In this study we have mainly concentrated on the Cucumis sativus and (but much less) Cucurbitaceae family several important vegetable crops. There are many reports on research conducted in Cucurbitaceae plant breeding programs on the ripening process, phloem transport, disease resistance, cold tolerance and fruit quality traits. This paper presents the role played by new omic technologies in the creation of knowledge on the mechanisms of the formation of the breeding features. The analysis of NGS (NGS-next generation sequencing) data allows the discovery of new genes and regulatory sequences, their positions, and makes available large collections of molecular markers. Genome-wide expression studies provide breeders with an understanding of the molecular basis of complex traits. Firstly a high density map should be created for the reference genome, then each re-sequencing data could be mapped and new markers brought out into breeding populations. The paper also presents methods that could be used in the future for the creation of variability and genomic modification of the species in question. It has been shown also the state and usefulness in breeding the chloroplastomic and mitochondriomic study.

High-density genetic map construction and gene mapping of pericarp color in wax gourd using specific-locus amplified fragment (SLAF) sequencing

Background

High-density map is a valuable tool for genetic and genomic analysis. Although wax gourd is a widely distributed vegetable of Cucurbitaceae and has important medicinal and health value, no genetic map has been constructed because of the lack of efficient markers. Specific-locus amplified fragment sequencing (SLAF-seq) is a newly developed high-throughput strategy for large-scale single nucleotide polymorphism (SNP) discovery and genotyping.

Results

In our present study, we constructed a high-density genetic map by using SLAF-seq and identified a locus controlling pericarp color in wax gourd. An F2 population of 140 individuals and their two parents were subjected to SLAF-seq. A total of 143.38 M pair-end reads were generated. The average sequencing depth was 26.51 in the maternal line (B214), 27.01 in the parental line (B227), and 5.11 in each F2 individual. When filtering low-depth SLAF tags, a total of 142,653 high-quality SLAFs were detected, and 22,151 of them were polymorphic, with a polymorphism rate of 15.42 %. And finally, 4,607 of the polymorphic markers were selected for genetic map construction, and 12 linkage groups (LGs) were generated. The map spanned 2,172.86 cM with an average distance between adjacent markers for 0.49 cM. The inheritance of pericarp color was also studied, which showed that the pericarp color was controlled by one single gene. And based on the newly constructed high-density map, a single locus locating on chromosome 5 was identified for controlling the pericarp color of wax gourd.

Conclusions

This is the first report of high-density genetic map construction and gene mapping in wax gourd, which will be served as an invaluable tool for gene mapping, marker assisted breeding, map-based gene cloning, comparative mapping and draft genome assembling of wax gourd.

Electronic supplementary material

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

Fine genetic mapping of target leaf spot resistance gene cca-3 in cucumber, Cucumis sativus L

The cucumber target leaf spot resistance gene cca - 3 was fine mapped in a 79-kb region harboring a CC-NB-ARC type R gene that may be responsible for the hypersensitive responses to infection of the target leaf spot pathogen in cucumber. The target leaf spot (TLS) is one of the most important foliar diseases in cucumber (Cucumis sativus L.). In this study, we conducted fine genetic mapping of a simply inherited recessive resistance gene, cca-3 against TLS with 193 F2:3 families and 890 F2 plants derived from the resistant cucumber inbred line D31 and the susceptible line D5. Initial mapping with microsatellite markers and bulked segregant analysis placed cca-3 in a 2.5-Mbp region of cucumber chromosome 6. The D5 and D31 lines were re-sequenced at 10× genome coverage to explore new markers in the target region. Genetic mapping in the large F2 population delimited the cca-3 locus in a 79-kb region with flanking markers Indel16874230 and Indel16953846. Additional fine mapping and gene annotation in this region revealed that a CC-NB-ARC type resistance gene analog, Csa6M375730, seems to be the candidate gene for cca-3. One single nucleotide polymorphism (SNP) was found in the NB-ARC domain of this candidate gene sequence between D31 and D5 that may lead to amino acid change, thus altering the function of the conserved NB-ARC motif. This SNP was validated in the segregating population as well as 24 independent cucumber lines. There was significantly higher level of cca-3 expression in the leaves of D5 (susceptible) than in D31 (resistant), and the expression level was positively correlated with the areas of necrotic spots on leaves after inoculation. It seems the cca-3 resistance gene was able to induce hypersensitive responses to the infection by TLS pathogen.

Construction of a high-density genetic map and mapping of a sex-linked locus for the brown alga Undaria pinnatifida (Phaeophyceae) based on large scale marker development by specific length amplified fragment (SLAF) sequencing

BACKGROUND

Undaria pinnatifida is an important economic brown alga in East Asian countries. However, its genetic and genomic information is very scarce, which hinders further research in this species. A high-density genetic map is a basic tool for fundamental and applied research such as discovery of functional genes and mapping of quantitative trait loci (QTL). In this study the recently developed specific length amplified fragment sequencing (SLAF-seq) technology was employed to construct a high-density genetic linkage map and locate a sex determining locus for U. pinnatifida.

RESULTS

A total of 28.06 Gb data including 140.31 M pair-end reads was obtained. After linkage analysis 4626 SLAF markers were mapped onto the genetic map. After adding the sex linked simple sequence repeat (SSR) marker [GenBank:AY738602.1], the final genetic map was 1816.28 cM long, consisting of 30 linkage groups with an average distance of 0.39 cM between adjacent markers. The length of LGs ranged from 20.12 to 106.95 cM. A major sex associated QTL was mapped to LG22 within a window starting at 29.01 cM and ending at 68.81 cM with a total of 68 SLAF markers. The SSR marker and five SLAF markers (Marker6556, 19020, 43089, 60771 and 26359) were identified as tightly sex-linked markers, as indicated by the absence of recombination between them and the sex phenotype. These markers were located at the position of 59.50 cM, which was supposed to be the sex determining region.

CONCLUSIONS

A high-density genetic linkage map was constructed using SLAF-seq technique and F1 gametophyte population for the first time in the economically important brown alga U. pinnatifida. For the first time, a major sex associated QTL suggesting a sex determining region was mapped to a single LG. This map will facilitate the further fundamental and applied research such as QTL mapping and map-based gene clone in U. pinnatifida and provide a reference for studies in other kelp species.

QTL mapping of cucumber fruit flesh thickness by SLAF-seq

Cucumber is an agriculturally and economically important vegetable crop worldwide. Fruit flesh thickness is an important trait for cucumber and also a central determinant of yield, yet little is known about the underlying mechanism of this trait. In this study, bulked segregant analysis (BSA) combined with specific length amplified fragment sequencing (SLAF-seq) was applied to finely map the gene that underlies fruit flesh thickness in cucumber. A 0.19-Mb-long quantitative trait locus on chromosome 2 controlling fruit flesh thickness (QTL fft2.1) was identified and further confirmed by simple sequence repeat (SSR) marker-based classical QTL mapping in 138 F₂ individuals. Gene prediction of this 0.19-Mb region identified 20 genes. Quantitative RT-PCR revealed higher expression levels of Csa2 M058670.1 (SET domain protein-lysine methyltransferase) in D8 (thick fruit flesh parent) compared with that in XUE1 (thin fruit flesh parent) during fruit development. Sequence alignment analysis of Csa2M058670.1 from thick and thin fruit flesh cucumber lines revealed a 4-bp deletion mutation in the promoter region of this candidate gene, which may result in the loss of Csa2M058670.1 activation in thin fruit flesh lines. The data presented herein suggest that Csa2M058670.1 is a possible candidate gene for controlling flesh thickness in cucumber.