Genetic linkage analysis using DNA markers in sweetpotato.

Sweetpotato is one of the most important food crop species in the world, with more than 104,000,000 tons produced each year, and the breeding of superior cultivars with agronomically important traits, such as improved disease resistance, yield, and nutrient richness, is necessary, especially in developing countries. However, as a result of the crop's complex genomic architecture, which results from its hexaploidy (2n = 6× = 90), high heterozygosity, huge genome, and outcrossing nature, the analysis of genetic linkage in sweetpotato has been challenging. In addition, the lack of whole genome sequences or gene annotations for cultivated hexaploids has interrupted the validation of mapped QTL regions and gene cloning. In spite of these technical difficulties, linkage map construction and QTL mapping analysis have been reported. This review summarizes the results of these linkage analyses, which used SSR, AFLP, and retrotransposon-based molecular markers, and describes future directions for the genetic analysis and marker-assisted breeding of this important but genetically complex crop species.

Introduction

The sweetpotato (Ipomoea batatas (L.) Lam., Convolvulaceae) is an edible tuberous root that was originally found in Central America. Its global production volume is around 104,000,000 t (FAO 2014; http://www.fao.org/home/en/), with close to 75% of the amount produced in Asia, and China itself produces around 71,000,000 t, which accounts for 70% of the global production. Large quantities of the crop are also produced in Southeast Asian countries, such as Indonesia and Vietnam, and in African countries, such as Nigeria and Uganda. Therefore, the sweetpotato is a vital source of nutrients for human populations and is especially important for various developing countries in Africa and Southeast Asia.

In Japan, on the other hand, the annual production of sweetpotato is ~886,500 t (FAO 2014), and most of the crop is produced in Kagoshima, Ibaraki, and Chiba (JRTA 2010). The form and characteristics currently desired in Japan are not simply a matter of yield, but also run the gamut from taste, nutritional qualities, and processing. Of the cultivars cultivated for consumption in Japan, intensely sweet cultivars, such as ‘Beniharuka’, are preferred in addition to the long-cultivated ‘Kokei No. 14’, and cultivation also leans toward types that are high in sugar content and superior in flavor (Kai ). Furthermore, islands between Tanegashima in Kagoshima prefecture and the Nansei Islands of Okinawa have long been used for growing anthocyanin-containing cultivars, such as the purple-fleshed sweetpotato (JRTA 2010).

In recent years, a number of other cultivars have also been grown, in order to facilitate the expansion of uses for sweetpotatoes by producing cultivars that are well suited for the crop’s various functions. For example, purple sweetpotato cultivars, such as ‘Ayamurasaki’ and ‘Akemurasaki’, are used as raw material for producing natural anthocyanin colorants (Sakai , Yamakawa ), and other purple sweetpotato cultivars, such as the ‘Purple Sweet Lord’, are used for food (Tamiya ). In addition, several cultivars also yield tuberous roots with high levels of β-carotene, and there is currently an expectation that the demand for such functional foods will expand. Furthermore, the ‘J-Red’ cultivar is used for juice (Yamakawa ), the ‘Sunny Red’ cultivar is used for powders (Yamakawa ), and the ‘Ayakomachi’ cultivar is used in salads and other cooking materials (Kai ). In addition, in order to expand demand among regular households, various cultivars are also being grown in response to changes in consumer demand, such as the ‘Quick Sweet’ cultivar, which demonstrates an appealing flavor, even when cooked with a microwave (Katayama ).

In addition to meeting consumer demand, another important objective in the domestic breeding of sweetpotato in Japan is to confer disease resistance. Particularly because sweetpotato has a tuberous root, and its underground portion is susceptible to soil-borne diseases. For example, although many Streptomyces spp. are useful, in that they produce antibiotics, some are known to cause common scab on root crops, and the actinobacteria S. ipomoea causes sweetpotato soil rot (Person and Martin 1940). In the case of the latter, when sweetpotato is invaded by bacteria, the underground portion blackens and starts to decay, round black specks appear on the tuberous root, and the plant decomposes (JRTA 2010, Person and Martin 1940). Furthermore, the sweetpotato loses considerable commercial value if there are any black specks. Unfortunately, most domestic cultivars are susceptible to the disease. Thus, the breeding of disease-resistant cultivars is desired. Therefore, when considering both consumer demand and disease resistance, there are various plant breeding objectives in sweetpotato.

The construction of genetic linkage maps is important for gene cloning, quantitative trait locus (QTL) detection, comparative genomic research, and marker-assisted breeding. However, in sweetpotato, linkage map construction has been challenging, owing to the species’ complex genomic architecture: high heterozygosity, general self-incompatibility, and hexaploidy, with a large number of small chromosomes (2n = 6× = 90), as mentioned above. For linkage map construction in such outcrossing heterozygous species, Grattapaglia and Sederoff (1994) developed a two-way pseudo-test cross strategy, in which linkage analysis is conducted for each parent separately, since dominant markers that are heterozygous in one parent and recessive homozygous in the other parent will segregate in the F1 generation, resulting in the development of two parental linkage maps. In hexaploids like sweetpotato, only simplex, duplex, and triplex alleles can produce polymorphisms in the progeny, and alleles at higher doses produce non-polymorphism. Based on the allele dosage (simplex, duplex, or triplex), several variations in segregation ratios are detected, and these ratios are affected by cytological characteristics (autohexaploid, tetradiploid, or allohexaploid) (Table 1).

Table 1

Expected segregation ratios (presence:absence) for the inheritance of a dominant marker in a hexaloid

Marker dosage	Autohexaploid (hexasomic)	Tetradiploid (tetradisomic, tetrasomic, disomic)	Allohexaploid (disomic)
Simplex	Aaaaaa 1:1	Aaaa aa 1:1	Aa aa aa 1:1
		aaaa Aa 1:1
Duplex	AAaaaa 4:1	AAaa aa 5:1a	Aa Aa aa 3:1
		Aaaa Aa 3:1b	AA aa aa 1:0
		aaaa AA 1:0c
Triplex	AAAaaa 19:1	AAAa aa 1:0	Aa Aa Aa 7:1
		AAaa Aa 11:1	AA Aa aa 1:0
		Aaaa AA 1:0

Tetrasomic inheritance.

Tetradisomic inheritance.

Disomic inheritance.

However, the ploidy type of sweetpotato has been a mystery for many years, and it has not yet been clearly determined. A number of studies have suggested that sweetpotato is allohexaploid (Jones 1965, Magoon , Sinha and Sharma 1992, Ting and Kehr 1953), whereas Shiotani and Kawase (1989) proposed that the genome constitution of sweetpotato should be B1B1B2B2B2B2 based on the frequent occurrence of tetravalents and hexavalents. On the other hand, several relatively recent genetic linkage analyses suggest that sweetpotato is autohexaploid, with some preferential pairing (Cervantes-Flores , Kriegner , Ukoskit and Thompson 1997, Zhao ).

As mentioned above, in outcrossing heterozygous species, such as sweetpotato, a two-way pseudo-test cross strategy has been applied to construct linkage maps, in which dominant markers that are heterozygous in one parent and recessive homozygous (=null) in the other parent are analyzed. During the construction of linkage maps in hexaploids, marker types are classified according to their segregation ratios in the mapping population for three hypothetical cytological types (autohexaploid, tetradiploid, and allohexaploid) : (a) simplex markers, which are present as a single copy in one parent and absent in the other parent, segregate in a 1:1 (presence:absence) ratio, regardless of the cytological types; (b) duplex markers, which are present as two copies in one parent and absent in the other parent, segregate in a 4:1 (hexasomic), 5:1 (tetrasomic), or 3:1 (disomic or tetradisomic) ratio; and (c) triplex markers, which are present as three copies in one parent and absent in the other parent, segregate in a 19:1 (hexasomic), 11:1 (tetradisomic), or 7:1 (disomic) ratio (Table 1) (Jones 1967), whereas double simplex markers occur as single copies in both parents and segregate in a 3:1 ratio (Certavantes-Flores , Kriegner ).

Because the simple 1:1 segregation ratio is only detected in simplex markers, for any cytological characteristics (autohexaploid, tetradiploid, or allohexaploid), simplex markers are highly useful for linkage analysis, and in linkage map construction, simplex markers are used to build an initial framework map, followed by the insertion of duplex, triplex, and double simplex markers. Thus, obtaining a large number of simplex markers is critical for forming the framework map, which is a basis of the linkage groups, and duplex and triplex markers are used for identifying the corresponding homologous linkage groups in the map. Therefore, the screening and genotyping of polymorphic markers based on the calculation of segregation ratios is apparently more tedious in polyploids than in diploids. However, although genetic and linkage analysis in sweetpotato has lagged far behind that in other plant species, both linkage map construction and QTL mapping in the species have been reported. In this review, we summarize the results of previous genetic studies and describe future perspectives on genetic study in sweetpotato.

1. Linkage map construction in sweetpotato

To date, several studies have constructed genetic linkage maps in sweetpotato (Cervantes-Flores , Chang , Kriegner , Li , Monden , Ukoskit and Thompson 1997, Zhao ) (Table 2). The first one was reported by Ukoskit and Thompson (1997), who constructed low-density linkage maps using randomly amplified polymorphic DNA (RAPD) markers, with a mapping population of 76 progenies derived from ‘Vardaman’, which exhibits high early root yield, drought tolerance, and susceptibility to root-knot nematode (RKN), and ‘Regal’, which is resistant to disease, RKN, and insects. Ukoskit and Thompson analyzed the segregation ratios of 196 RAPD markers in each parent and then used the markers to construct linkage maps. Of the 196 markers, 76 (74.5%) and 67 (71.3%), 17 (16.7%) and 21 (22.3%), and four (3.9%) and three (3.2%) were simplex, duplex, and triplex markers, in ‘Vardaman’ and ‘Regal’, respectively. In order to develop linkage maps, the linkage relationships of simplex RAPD markers were first analyzed and used to construct a framework map for each parent, and then duplex and triplex markers were inserted into the framework map manually, which increased the mapping efficiency. In linkage maps that only contained simplex markers, 25 and 20 markers were included in the ‘Vardaman’ and ‘Regal’ maps, respectively, and adding duplex and triplex markers increased the total number of mapped markers to 48 and 46. In addition, the authors also reported that the main advantage of using non-simplex markers was the identification of homologous linkage groups, of which three were identified in both maps. This work was the first to report using the pseudo-test cross mapping strategy in sweetpotato, which enabled the construction of genetic linkage maps and demonstrated the high probability that sweetpotato is autohexaploid, not allohexaploid.

Table 2

Summary of linkage maps constructed in sweetpotato

Types of markers	Total number of markers	Total number of linakge groups	Totel length (cM)	Average length of markers (cM)	Parent cultivar	Populations	Software	Paper
RAPD	48	18	474.6	9.9	Vardaman	76 progenies derived from Vardaman × Regal	MAPMAKER 3.0	Ukoskit and Thompson (1997)
	46	16	489.1	10.6	Regal
AFLP	632	90	3655.6	5.8	Tanzania	94 progenies derived from Tanzania × Bikilamialiya	JoinMap 3.0	Kriegner et al. (2003)
	435	80	3011.5	6.9	Bikilamaliya
AFLP	726	90	5276	4.8	Beauregard	240 progenies derived from Tanzania × Beauregard	JoinMap 3.0 and MAPMAKER/EXP 3.0	Cervantes-Flores et al. (2008)
	947	86	5792	4.5	Tanzania
ISSR	37	9	479.8	12.9	Nancy Hall	119 progenies derived from NH × TN27	MAPMAKER 3.0	Chang et al. (2009)
	47	12	853.5	17.7	Tainung 27	112 progenies derived from TN27 × NH
SRAP	473	81	5802.46	10.2	Luoxushu 8	240 progenies derived from Luoxushu 8 × Zhengshu 20	JoinMap 3.0	Li et al. (2010)
	328	66	3967.9	12.0	Zhengshu 20
AFLP, SSR	2077	90	8184.5	3.9	Xushu 18	202 progenies derived from Xushu 18 × Xu781	JoinMap 3.0	Zhao et al. (2013)
	1954	90	8151.7	4.2	Xu 781
Retrotransposon-based	124	43	931.5	11.6	Purple Sweet Lord	98 progenies derived from PSL× 90IDN-47	JoinMap 4.1	Monden et al. (2015)
	122	47	734.3	9.8	90IDN-47

Afterward, several researchers reported the construction of genetic linkage maps in sweetpotato using molecular markers, such as amplified fragment length polymorphism (AFLP), simple sequence repeat (SSR), and sequence related amplified polymorphism (SRAP) markers. In particular, the application of AFLP markers enabled researchers to dramatically increase the number of mapped markers and linkage groups. The method used for linkage map construction was basically the same as that reported by Ukoskit and Thompson (1997): simplex markers were first used to construct a framework map, and then the duplex, triplex, and double-simplex markers were inserted, in order to detect homologous groups. Kriegner developed linkage maps using AFLPs in an F1 mapping population that was derived from a cross between ‘Tanzania’ and ‘Bikilamliya’ and obtained a total of 90 and 80 linkage groups, which included 632 and 435 markers in the ‘Tanzania’ and ‘Bikilamliya’ maps, respectively, which had total lengths of 3655.6 and 3011.5 cM (Table 2). Meanwhile, Cervantes-Flores developed another set of linkage maps using AFLPs and a mapping population that was derived from a cross between ‘Tanzania’ and ‘Beauregard’. In the ‘Tanzania’ map, 86 linkage groups, which included 947 markers, were obtained, and the entire length of the maps was estimated as 5792 cM, with an average between-marker distance of 4.5 cM (Table 2). In the ‘Beauregard’ map, 90 linkage groups, which included 726 markers, were obtained, and the entire length of the maps was 5276 cM, with an average between-marker distance of 4.8 cM (Table 2). In addition, these studies also supported the conclusion that sweetpotato is autohexaploid, not allohexaploid; however, some preferential pairing between certain chromosomes was reported, as well.

Moreover, the genetic linkage map in sweetpotato with the highest marker density was constructed by Zhao , using AFLP and SSR markers in an F1 mapping population of ‘Xushu 18’ and ‘Xu 781’. In the ‘Xushu 18’ map, 90 linkage groups were obtained, which included 2077 markers, and the total length of the maps was 8184.5 cM, with an average between-marker distance of 3.9 cM (Table 2). Meanwhile, in the ‘Xu 781’ map, 90 linkage groups were obtained, which included 1954 markers, and the total length of the maps was 8151.7 cM, with an average between-marker distance of 4.2 cM (Table 2). This was the first study to establish 90 complete linkage maps, a number which corresponds to the species’ actual chromosome number, and remains the only study to have accomplished this to date. However, the AFLP markers in Zhao et al.’s linkage map exhibited a clustered distribution, and in several linkage groups, such as XUSHU18(02.09), XUSHU18(10.51), XU781(01.03), and XU781(05.28), many markers were clustered. This clustering could have resulted from a biological basis, thus, reflecting the suppressed or reduced genetic recombination in heterochromatic regions around the centromeres and telomeric regions (Tanksley ), but other authors have suggested that AFLP markers generated using EcoRI and MseI tend to exhibit uneven distributions as a result of the frequency of EcoRI-recognized A/T-rich regions in pericentromeric heterochromatin (Isidore , Truco ).

In addition, the genotyping and screening of polymorphic classical DNA markers, such as AFLPs and SSRs, can be tedious and time-consuming, although this depends on the relatedness of the parental cultivars and the types of markers used. For example, Zhao reported that a total of 2893 AFLP primer combinations were screened using parental cultivars and a small subset of progeny, but only 304 (10.5%) produced good quality polymorphic and segregating bands. Meanwhile, in the case of SSR markers, only 68 of 173 SSR primer pairs (39.3%) exhibited polymorphic bands of acceptable quality and were used to develop SSR markers in the entire mapping population.

In recent years, efficient and genome-wide genotyping has been achieved in dramatically less time by using next-generation sequencing (NGS) technologies (Goodwin , Heard , Nielsen ), and NGS-based genetic linkage analyses have been pursued in many crop species (Baird , Baxter , Bielenberg , Chutimanitsakun , Gardner , Guajardo , İpek , Ma , Poland , Pootakham , Russell , Shirasawa , Spindel , Verma , Wang , Ward , Wu , Xiao ). However, few studies have used NGS-based approaches for the genetic analysis of sweetpotato.

In one NGS-based study in sweetpotato, Monden analyzed comprehensive retrotransposon insertion polymorphisms and developed genetic linkage maps. Because a large number of retrotransposon insertions are distributed throughput the genomes of eukaryotes and because such insertions are genetically inherited, retrotransposon insertion polymorphisms can be used to develop molecular markers (Kalendar 2011, Kalendar , Kumar and Hirochika 2001, Poczai , Schulman ). For example, in cultivated strains, Rtsp-1, which is one of the active retrotransposon families in sweetpotato, exhibits high rates of insertion polymorphisms (Monden , Tahara ), which indicates that the insertion sites are highly likely to be useful for linkage map construction. Accordingly, the authors performed high-throughput Illumina NGS of Rtsp-1 insertion sites in a number of cultivars and found that the presence or absence of insertion sites could be used as the basis for efficient genotyping. However, in this case, the number of developed markers was relatively small, and the coverage was also low since the method relied on the copy number of a targeted retrotransposon family (~250 in this case). Thus, further research is needed to increase the number of mapped markers by analyzing insertion sites in other retrotransposon families with high copy numbers or by introducing different types of molecular markers. An advantage of the research was that NGS technology allowed the authors to reduce the time and labor required for genotyping, and, interestingly, the retrotransposon-based markers also indicated an extremely high proportion of simplex markers (~90%) (Monden ), which demonstrates the utility of retrotransposon-based markers for linkage map construction in polyploid species. Therefore, more research focused on higher-density linkage map construction is needed for precise QTL detection, marker-assisted selection, and breeding.

2. QTL analysis in sweetpotato

Several QTL mapping analyses have been conducted using the constructed linkage maps (Cervantes-Flores , Zhao ) for identifying agronomically important genes or genomic regions, and so far, such analyses have targeted economically important traits, such as nematode resistance, starch content, dry-matter content, β-carotene content, and root yield. The first QTL analysis in sweetpotato was conducted by Cervantes-Flores in order to identify target genes or chromosomal regions associated with RKN resistance, by using a mapping population of ‘Tanzania’ and ‘Beauregard’. In sweetpotato, RKN is known for causing significant reductions in yield and for damaging tuberous roots, and the characteristic symptoms of RKN damage are round to spindle-shaped swellings (galls) on fibrous roots and cracks on fleshy tuberous roots (Lawrence , Nakayama ). Although nematicides provide a fairly effective means of controlling RKN, they are expensive and highly neurotoxic. Thus, it is preferable to develop resistant cultivars as more economical and sustainable strategy. Although the resistance of sweetpotato to RKN is not well understood, previous studies have suggested that several genes may be involved and that resistance to different RKN species should be conferred by different genes (Cervantes-Flores , 2002b, Sano ). In the study of Cervantes-Flores , an F1 mapping population was developed from ‘Tanzania’ and ‘Beauregard’ parental cultivars, owing to previous reports of high levels of resistance in ‘Tanzania’ to multiple RKN species and the high susceptibility of ‘Beauregard’ (Cervantes-Flores , 2002b). Genetic maps constructed by Cervantes-Flores were used for QTL mapping, and RKN resistance was assessed as the number of egg masses observed on the roots. The frequency distribution of resistance in the mapping population (n = 240) was highly skewed, with most of the progeny exhibiting moderate to high levels of resistance, which suggested that RKN resistance was a quantitative trait that was likely controlled by a few major genes. The subsequent QTL analysis identified nine chromosome regions that were associated with RKN resistance, seven of which (rkn1 to rkn7) were observed in ‘Tanzania’ and two of which (rkn8 and rkn9) were observed in ‘Beauregard’ (Table 3). All the QTLs observed in ‘Tanzania’ and their interactions explained ~40% of the variation in egg masses; however, no single major QTL explained more than 15% of the variation. One the other hand, the study also identified three unmapped duplex markers that were highly associated with RKN resistance, accounting for as much as 45% of the variation in egg masses. Therefore, the three unmapped markers might be associated with a locus or loci with major effects (data not shown); however, further research was needed to validate their QTLs, through the inclusion of additional markers. One the other hand, another study conducted QTL mapping analysis for targeting the southern root-knot nematode (SRKN) resistance, and identified genomic regions with a large effect on resistance, and successfully developed sequence-specific PCR markers associated with QTL(s) for resistance to multiple races of SRKN in sweetpotato (Nakayama ).

Table 3

Summary of QTL resultes conducted in sweetpotato

Trait	Cultivar	Name of QTL	Linkage group	Closest marker	LOD	(Additive) effect	R2 (%)	Paper
Root-knot nematode resistance	Tanzania	rkn1	T0.5.26	E42M6022	–	−0.16	4.0	Cervantes-Flores et al. (2008)
		rkn2	T07.38	E35M4414	–	−0.18	3.7
		rkn3	T07.41	E38M4512	–	−0.30	9.1
		rkn4	T08.46	E36M3811	–	−0.21	5.4
		rkn5	T01.01	E32M4920	–	0.16	3.7
		rkn6	T07.37	E46M3201	–	0.20	4.4
		rkn7	T07.39	E32M3722	–	0.32	11.5
	Beauregard	rkn8	B04.24	E40M6008	–	−0.20	4.5
		rkn9	B07.39	R42M3525	–	−0.17	2.2
Top weight	Nancy Hall	TW	NH Group 9	ISSR 807-7	4.9	198.72	20.5	Chang et al. (2009)
	Tainung 27	TW1	TN27 Group 11	ISSR 857-2	5.2	114.27	16.0
		TW2	TN27 Group 11	ISSR 824-3	5.8	140.44	20.2
		TW3	TN27 Group 12	ISSR 857-5	6.0	127.68	24.7
Root weight	Nancy Hall	RW1	NH Group 8	ISSR 822-4	3.2	121.24	25.1
		RW2	NH Group 9	ISSR 807-7	3.4	109.38	29.5
	Tainung 27	RW	TN27 Group7	ISSR 825-12	2.5	328.73	28.2
Root skin color	Nancy Hall	SC1	NH Group 1	ISSR 809-3	5.8	2.80	20.5
		SC2	NH Group 3	ISSR 826-4	4.1	2.41	14.2
		SC3	NH Group 4	ISSR 827-6	4.3	1.83	15.9
	Tainung 27	SC1	TN27 Group 1	ISSR 819-6	2.9	1.92	20.4
		SC2	TN27 Group 2	ISSR 811-5	5.0	1.18	28.3
		SC3	TN27 Group 11	ISSR 857-2	4.2	1.71	22.8
Root flesh color	Nancy Hall	FC1	NH Group 1	ISSR 809-7	2.7	1.54	10.7
		FC2	NH Group 2	ISSR 811-1	5.9	1.98	27.8
		FC3	NH Group 6	ISSR 809-14	2.5	2.54	14.1
		FC4	NH Group 7	ISSR 825-16	2.5	2.09	20.4
	Tainung 27	FC1	TN27 Group 3	ISSR 810-4	3.9	1.62	29.9
		FC2	TN27 Group 12	ISSR 857-6	3.8	1.01	22.8
Root shape	Nancy Hall	RS1	NH Group 2	ISSR 811-2	9.1	1.51	20.0
		RS2	NH Group 3	ISSR 826-2	9.4	1.80	26.1
		RS3	NH Group 6	ISSR 825-15	8.9	2.31	29.8
Root number	Nancy Hall	RN	NH Group 9	ISSR 823-7	2.7	10.62	14.8
Dry-matter content	Beauregard	dm1	B01.03	E42M3421	–	nagtive	–	Cervantes-Flores et al. (2011)
		dm2	B04.23	E43M5403	–	negative	–
		dm3	B05.26	E35M4511	–	positive	–
		dm4	B07.40	E32M3202	–	positive	–
		dm5	B11.61	E40M4010	–	positive	–
		dm6	B11.62	E36M5103	–	negative	–
		dm7	B12.70	E34M4906	–	negative	–
		dm8	B00.89	E36M5408	–	positive	–
	Tanzania	dm9	T01.05	E35M3603	–	positive	–
		dm10	T02.07	E43M3524	–	negative	–
		dm11	T05.25	E36M3808	–	positive	–
		dm12	T06.32	E31M3208	–	positive	–
		dm13	T07.37	E46M6011	–	positive	–
Starch content	Beauregard	sta1	B01.03	E42M3421	–	negative	–
		sta2	B04.23	E43M5403	–	negative	–
		sta3	B07.39	E42M3525	–	positive	–
		sta4	B07.40	E32M3202	–	positive	–
		sta5	B11.61	E40M4010	–	positive	–
		sta6	B11.62	E36M5103	–	negative	–
		sta7	B12.70	E34M4906	–	negative	–
	Tanzania	sta8	T01.04	E37M4203	–	positive	–
		sta9	T02.07	E43M3524	–	negative	–
		sta10	T05.25	E36M3808	–	positive	–
		sta11	T06.32	E31M3208	–	positive	–
		sta12	T71	E32M4908	–	negative	–
β-Carotene content	Beauregard	caro1	B04.23	E43M5403	–	positive	–
		caro2	B08.48	E38M3725	–	negative	–
		caro3	B11.62	E36M5103	–	positive	–
		caro4	B12.69	E44M4902	–	positive	–
	Tanzania	caro5	T13.74	E45M3611	–	positive	–
		caro6	T13.76	E40M3105	–	positive	–
		caro7	T78	E46M3901	–	positive	–
		caro8	T82	E36M4015	–	negative	–
Dry-matter content	Xu 781	DMM_1	XU781 (02.08)	E10M2-10s	4.7	positive	16.1	Zhao et al. (2013)
		DMM_2	XU781 (02.08)	E45M58-14s	4.0	positive	16.3
		DMM_3	XU781(03.13)	E53M27-4s	3.0	positive	11.1
		DMM_4	XU781(03.16)	E51M8-14s	3.4	positive	17.7
		DMM_5	XU781(04.19)	E5M24-19s	2.8	positive	13.0
		DMM_6	XU781(04.21)	E53M41-14s	2.5	positive	9.8
		DMM_7	XU781(05.25)	E26M1-6s	2.6	negative	12.4
		DMM_8	XU781(05.26)	E7M51-6s	3.4	positive	12.8
		DMM_9	XU781(06.31)	E42M64-4s	4.3	positive	21.6
		DMM_10	XU781(06.34)	E51M28-6d	3.4	positive	10.7
		DMM_11	XU781(06.36)	C55-5s	4.4	positive	45.1
		DMM_12	XU781(07.37)	E4M24-26s	2.8	positive	9.6
		DMM_13	XU781(07.40)	E27M62-2d/E36M63-10ds	5.3	positive	29.5
		DMM_14	XU781(08.44)	E53M7-3s	3.7	positive	13.4
		DMM_15	XU781(09.50)	E15M57-4d/C60-1s	3.4	positive	20.7
		DMM_16	XU781(09.50)	C55-11s	3.7	negative	13.2
		DMM_17	XU781(09.51)	E40M11-8s	3.4	positive	13.0
		DMM_18	XU781(11.57)	E53M16-24s	3.1	positive	11.2
		DMM_19	XU781(00.73)	E60M5-16d	2.8	positive	9.7
		DMM_20	XU781(00.74)	E1M34-7s	4.2	positive	15.2
		DMM_21	XU781(00.79)	C20-2s	2.9	positive	9.3
		DMM_22	XU781(00.80)	E51M28-6d/C24-4d	3.5	positive	31.2
	Xushu 18	DMF_1	XUSHU18(02.09)	E41M59-2s	2.8	positive	16.7
		DMF_2	XUSHU18(03.13)	E7M9-4ds	2.7	negative	12.1
		DMF_3	XUSHU18(06.35)	E55M16-14ds	2.6	negative	9.0
		DMF_4	XUSHU18(06.35)	E12M59-12d	3.0	negative	25.3
		DMF_5	XUSHU18(13.61)	E52M13-10d	2.9	negative	10.5
Starch content	Xu 781	DMFN_1	XU781(01.05)	E31M19-21s	5.0	positive	20.5	Xiao-xia et al. (2013)
		DMFN_2	XU781(02.07)	E51M28-6d	4.1	positive	14.4
		DMFN_3	XU781(03.15)	E17M55-18d	2.6	positive	10.6
		DMFN_4	XU781(07.40)	E36M63-10ds	7.2	positive	38.8
		DMFN_5	XU781(08.45)	E50M49-3d	4.6	positive	14.0
		DMFN_6	XU781(10.54)	C61_2d	2.5	positive	9.1
	Xushu 18	DMMN_1	XUSHU18(05.27)	E1M31-17ds	3.0	negative	16.1
		DMMN_2	XUSHU18(06.31)	E50M12-8s	2.5	negative	14.3
Storage root yield	Xu 781	YIEM_1	XU781(04.20)	E50M49-4d	10.1	negative	41.8	Li et al. (2014)
		YIEM_2	XU781(07.38)	IB_S10_9d	9.5	negative	37.0
		YIEM_3	XU781(07.41)	E54M5-4d	11.1	positive	51.2
		YIEM_4	XU781(07.42)	E21M8-23d	10.8	negative	46.8
		YIEM_5	XU781(12.61)	IB_S10_4d	3.0	negative	44.1
	Xushu 18	YIEF_1	XUSHU18(01.01)	E24M1-17d	5.3	negative	36.3
		YIEF_2	XUSHU18(03.16)	E53M30-12d	4.6	positive	59.3
		YIEF_3	XUSHU18(11.55)	C27_8s	9.1	positive	42.9
		YIEF_4	XUSHU18(13.62)	E63M3-10d	5.4	positive	44.0

‘–’: no information in the paper.

Subsequently, other QTL analyses have been conducted to identify genes and/or chromosomal regions that control economically important traits, such as root yield, dry-matter content, and starch content. For example, Cervantes-Flores conducted QTL mapping analysis for dry-matter, starch, and β-carotene content using a mapping population of 240 plants that were derived from a cross between two parental cultivars with significantly different traits, namely ‘Tanzania’, which is a white cream-colored African landrace that has a high dry-matter content (>30%), and ‘Beauregard’, which is an American orange-fleshed cultivar with a low dry-matter content. The authors found that the distribution of dry-matter and starch content in the mapping population was continuous, whereas that of β-carotene content was highly skewed, and transgressive segregations were also observed in all three traits, which suggested that the traits were quantitative. In addition, 13, 12, and eight QTLs were identified for dry-matter, starch, and β-carotene content, respectively (Table 3), and among the three traits, dry-matter and starch content were highly correlated (r > 0.8), whereas β-carotene content was negatively correlated with both starch (r = −0.6) and dry-matter (r = −0.3) content. Therefore, it is reasonable to assume that most of the QTL regions that significantly affected starch content also affected dry-matter content.

Zhao developed the first map that included 90 complete sweetpotato linkage groups, and subsequent analyses conducted by Xiao-xia and Li were able to identify QTLs and co-localizing markers for agronomically important traits, such as dry-matter (Zhao ), starch content (Xiao-xia ), and yield (Li ). In these three studies, QTL mapping was conducted using the linkage map based on AFLP and SSR markers, to obtain a mapping population of 202 plants that were derived from a cross between ‘Xushu 18’, which is the most widely grown cultivar in China and has high yield and moderate starch and dry-matter content, and ‘Xu 781’, which has low yield and high starch and dry-matter content. The analysis identified 27, eight, and nine QTLs for dry-matter content, starch content, and root yield (Table 3), and both dry-matter and starch content were distributed normally, whereas root yield was abnormally distributed, as reported by Brondani in rice. However, transgressive segregation was observed in all three traits, which indicates that the three traits are quantitative, and as expected, a number of QTLs were detected for all three traits. For dry-matter content, 27 QTLs (DMM-1 to DMM-22 and DMF-1 to DMF-5) were identified, in which DMM-11 QTL accounted for the largest amount of phenotypic variation in ‘Xu 781’ (45.1%) and was co-localized with a single SSR marker (C55-5s) on map XU781(06.36). For starch content, eight QTLs (DMFN-1 to DMFN-6, DMMN-1, and DMMN-2) were identified, in which DMFN-4 accounted for the largest amount of phenotypic variation in ‘Xu 781’ (38.8%) and was co-localized with the E36M63-10ds marker on map XU781(07.40). For root yield, a total of 45 QTLs were detected, using phenotypic data from two different locations and three time points, and nine of the QTLs were identified at the same genomic regions in at least two different environments and identified as stable. In particular, the YIEF-2 QTL in ‘Xushu 18’ explained the largest amount of phenotypic variation (59.3%) and was co-localized with the E53M30-12d marker on map XUSHU18(03.16). In addition, the parental cultivars also differ in their resistance to root rot, stem nematodes, and sweetpotato feathery mottle virus (SPFMV), which suggests that similar QTL analyses could be used to identify chromosome regions that affect these traits, as well.

Because of the complex structure of the sweetpotato genome, genetic linkage analysis of sweetpotato has been more difficult than that of other crop species, and although several studies have performed linkage and QTL mapping in sweetpotato, the number of such studies is quite smaller than that in other crop species. In addition, the lack of a whole genome sequence or gene annotation in sweetpotato has hindered the validation of mapped regions and then isolation of target genes. Therefore, further research, especially to establish a genome sequence and gene annotation, is needed to validate QTL mapping results and to assist in marker-assisted sweetpotato breeding programs.

3. Future research

Recently, NGS technologies have been applied to whole genome sequencing, genome-wide genetic analysis, and high-density linkage map construction, and in sweetpotato, Hirakawa reported the de novo whole-genome sequencing of Ipomoea trifida (H. B. K.) G. Don., which is considered to be the closest diploid ancestor of the hexaploid cultivated sweetpotato. In I. trifida, several ploidy levels from diploid (2n = 2× = 30) to hexaploid (2n = 6× = 90) exist (Kobayashi 1983), and it is considered a wild relative of the hexaploid sweetpotato, I. batatas, because the two species can hybridize (Komaki and Katayama 1999, Orjeda , Shiotani and Kawase 1989). In addition, since several molecular and cytogenetic studies have confirmed that the two species are closely related (Huang and Sun 2000, Jarret and Austin 1994, Komaki , Roullier , Srisuwan ), the whole genome sequence information of this species is expected to dramatically promote genomic and genetic studies in I. batatas.

Higher density linkage maps will also be constructed by applying NGS-based genotyping methods, such as restriction site-associated DNA sequencing (RAD-seq), double digest RAD-seq (ddRAD-seq), and genotyping-by-sequencing (GBS), all of which facilitate the rapid discovery of thousands of single nucleotide polymorphisms (SNPs) in large populations by sequencing only the target DNA regions that flank specific restriction enzyme sites (Baird , Elshire , Miller , Peterson , Van Tassell ). To date, these methods have already yielded a number of high-density, SNP-based linkage maps, even in non-model species (Baird , Baxter , Bielenberg , Chutimanitsakun , Gardner , Guajardo , İpek , Ma , Poland , Pootakham , Russell , Shirasawa , Spindel , Verma , Wang , Ward , Wu , Xiao ), and in sweetpotato, genotyping and linkage map construction with RAD-seq is already underway (Shirasawa ).

It is generally expected that higher density linkage maps will facilitate rapid and precise QTL detection, as well as the development of molecular markers associated with agronomically important traits, such as disease resistance, root yield, and nutrient content. On the other hand, novel NGS-based methods, such as genome wide association study, genomic selection, MutMap, and QTL-seq, have also been established and widely applied in several plant species (Abe , Das , Fekih , Huang , Illa-Berenguer , Lu , Singh , Spindel , Takagi , 2013b, Tian , Wang ), without the need to construct linkage maps. However, although these methods are extremely powerful for accelerating plant breeding programs, the application of such methods in sweetpotato are limited by several barriers related to its complex genome structure, heterozygosity, huge genome size, and polyploidy. In addition, the collection and maintenance of large mapping populations are time and labor intensive in sweetpotato, owing to crop’s size, cultivation form, and space requirements, thus further limiting the application of genetic studies. Therefore, researchers from different disciplines will need to collaborate, in order to decipher the crop’s complex genetic architecture, to facilitate genetic studies, and to develop marker-assisted breeding programs, by identifying genes and genomic regions that are involved in agronomically important traits.