Construction of a linkage map based on retrotransposon insertion polymorphisms in sweetpotato via high-throughput sequencing

Chromatographic printed array strip (C-PAS) method for cultivar-specific identification of sweetpotato cultivars 'Beniharuka' and 'Fukumurasaki'

Sweetpotato (Ipomoea batatas) cultivars grown in Japan are highly valued for their excellent sweetness, high quality, and good texture. The export volume of sweetpotato from Japan has been rising rapidly, with a 10-fold increase on a weight basis over the last 10 years. However, since sweetpotato is propagated vegetatively from storage roots, it is easy to cultivate and propagate this crop, prompting concerns that Japanese sweetpotato cultivars/lines are being exported overseas, cultivated without permission, or reimported. Therefore, a rapid and accurate cultivar identification methodology is needed. In this study, we comprehensively analyzed the insertion sites of Cl8 retrotransposon to develop a cultivar identification technique for the Japanese cultivars 'Beniharuka' and 'Fukumurasaki'. These two cultivars were successfully distinguished from other cultivars using a minimum of two marker sets. Using the chromatographic printed array strip (C-PAS) method for DNA signal detection, 'Beniharuka' and 'Fukumurasaki' can be precisely identified using a single strip of chromatographic paper based on multiplex DNA signals derived from the amplicons of the Cl8 insertion sites. Since this method can detect DNA signals in only ~15 minutes, we expect that our method will facilitate rapid, reliable, and convenient cultivar discrimination for on-site inspection of sweetpotato.

A target cultivar-specific identification system based on the chromatographic printed array strip method for eight prominent Japanese citrus cultivars

Citrus is a major cultivated crop in Japan, and new cultivars are of great interest in the Japanese and global market. Recently, the infringement of breeders' rights to citrus cultivars bred in Japan has become a problem related to the agricultural product export strategy promoted by the Japanese government. Cultivar identification systems using DNA markers are an effective tool for protecting breeders' rights. Here, a novel target cultivar-specific identification system using the chromatographic printed array strip method was developed for eight prominent Japanese citrus cultivars. A polymorphic InDel fragment specific to each cultivar was explored through the screening of published citrus InDel markers and next-generation sequencing of retrotransposon libraries. The cultivar-specific DNA marker set for each cultivar comprised 1-3 polymorphic InDel fragments in combination with a PCR-positive DNA marker for the ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit gene. The DNA markers were detected within 3 hours from DNA extraction to the detection by the C-PAS4 membrane stick following multiplex PCR. The developed system is superior as a convenient, rapid, and cost-effective DNA diagnostic method during inspection. The proposed target cultivar-specific identification system is expected to serve as an efficient tool for the injunction of suspicious registered cultivars, contributing to the protection of breeders' rights.

Exploiting the miniature inverted-repeat transposable elements insertion polymorphisms as an efficient DNA marker system for genome analysis and evolutionary studies in wheat and related species

Transposable elements (TEs) constitute ~80% of the complex bread wheat genome and contribute significantly to wheat evolution and environmental adaptation. We studied 52 TE insertion polymorphism markers to ascertain their efficiency as a robust DNA marker system for genetic studies in wheat and related species. Significant variation was found in miniature inverted-repeat transposable element (MITE) insertions in relation to ploidy with the highest number of "full site" insertions occurring in the hexaploids (32.6 ± 3.8), while the tetraploid and diploid progenitors had 22.3 ± 0.6 and 15.0 ± 3.5 "full sites," respectively, which suggested a recent rapid activation of these transposons after the formation of wheat. Constructed phylogenetic trees were consistent with the evolutionary history of these species which clustered mainly according to ploidy and genome types (SS, AA, DD, AABB, and AABBDD). The synthetic hexaploids sub-clustered near the tetraploid species from which they were re-synthesized. Preliminary genotyping in 104 recombinant inbred lines (RILs) showed predominantly 1:1 segregation for simplex markers, with four of these markers already integrated into our current DArT-and SNP-based linkage map. The MITE insertions also showed stability with no single excision observed. The MITE insertion site polymorphisms uncovered in this study are very promising as high-potential evolutionary markers for genomic studies in wheat.

Comprehensive survey of transposon mPing insertion sites and transcriptome analysis for identifying candidate genes controlling high protein content of rice

Rice is the most important crop species in the world, being staple food of more than 80% of people in Asia. About 80% of rice grain is composed of carbohydrates (starch), with its protein content as low as 7-8%. Therefore, increasing the protein content of rice offers way to create a stable protein source that contributes to improving malnutrition and health problems worldwide. We detected two rice lines harboring a significantly higher protein content (namely, HP5-7 and HP7-5) in the EG4 population. The EG4 strain of rice is a unique material in that the transposon mPing has high transpositional activity and high copy numbers under natural conditions. Other research indicated that mPing is abundant in the gene-rich euchromatic regions, suggesting that mPing amplification should create new allelic variants, novel regulatory networks, and phenotypic changes in the EG4 population. Here, we aimed to identify the candidate genes and/or mPing insertion sites causing high protein content by comprehensively identifying the mPing insertion sites and carrying out an RNA-seq-based transcriptome analysis. By utilizing the next-generation sequencing (NGS)-based methods, ca. 570 mPing insertion sites were identified per line in the EG4 population. Our results also indicated that mPing apparently has a preference for inserting itself in the region near a gene, with 38 genes in total found to contain the mPing insertion in the HP lines, of which 21 and 17 genes were specific to HP5-7 and HP7-5, respectively. Transcriptome analysis revealed that most of the genes related to protein synthesis (encoding glutelin, prolamin, and globulin) were up-regulated in HP lines relative to the control line. Interestingly, the differentially expressed gene (DEG) analysis revealed that the expression levels of many genes related to photosynthesis decreased in both HP lines; this suggests the amount of starch may have decreased, indirectly contributing to the increased protein content. The high-protein lines studied here are expected to contribute to the development of high protein-content rice by introducing valuable phenotypic traits such as high and stable yield, disease resistance, and abundant nutrients.

Mapping of Nematode Resistance in Hexaploid Sweetpotato Using an Next-Generation Sequencing-Based Association Study

The southern root-knot nematode (SRKN; Meloidogyne incognita) is a typical parasitic nematode that affects sweetpotato [Ipomoea batatas (L.) Lam.], causing a significant decrease in crop yield and commercial value. In Japan, the SRKN is classified into 10 races: SP1-SP5, SP6-1, SP6-2, and SP7-SP9, with the dominant race differing according to the cultivation area. Soil insecticides have previously been used to reduce the soil density of SRKNs; however, this practice is both costly and labor intensive. Therefore, the development of SRKN-resistant sweetpotato lines and cultivars is necessary. However, due to the complexity of polyploid inheritance and the highly heterogeneous genomic composition of sweetpotato, genetic information and research for this species are significantly lacking compared to those for other major diploid crop species. In this study, we utilized the recently developed genome-wide association approach, which uses multiple-dose markers to assess autopolyploid species. We performed an association analysis to investigate resistance toward SRKN-SP2, which is the major race in areas with high sweetpotato production in Japan. The segregation ratio of resistant and susceptible lines in the F₁ mapping population derived from the resistant "J-Red" and susceptible "Choshu" cultivars was fitted to 1: 3, suggesting that resistance to SP2 may be regulated by two loci present in the simplex. By aligning the double digest restriction-site associated DNA sequencing reads to the published Ipomoea trifida reference sequence, 46,982 single nucleotide polymorphisms (SNPs) were identified (sequencing depth > 200). The association study yielded its highest peak on chromosome 7 (Chr07) and second highest peak on chromosome 3 (Chr03), presenting as a single-dose in both loci. Selective DNA markers were developed to screen for resistant plants using the SNPs identified on Chr03 and Chr07. Our results showed that SRKN-SP2-resistant plants were selected with a probability of approximately 70% when combining the two selective DNA markers. This study serves as a model for the identification of genomic regions that control agricultural traits and the elucidation of their effects, and is expected to greatly advance marker-assisted breeding and association studies in polyploid crop species.

Assessment of genetic diversity and variety identification based on developed retrotransposon-based insertion polymorphism (RBIP) markers in sweet potato (Ipomoea batatas (L.) Lam.)

Sweet potato, a dicotyledonous and perennial plant, is the third tuber/root crop species behind potato and cassava in terms of production. Long terminal repeat (LTR) retrotransposons are highly abundant in sweet potato, contributing to genetic diversity. Retrotransposon-based insertion polymorphism (RBIP) is a high-throughput marker system to study the genetic diversity of plant species. To date, there have been no transposon marker-based genetic diversity analyses of sweet potato. Here, we reported a structure-based analysis of the sweet potato genome, a total of 21555 LTR retrotransposons, which belonged to the main LTR-retrotransposon subfamilies Ty3-gypsy and Ty1-copia were identified. After searching and selecting using Hidden Markov Models (HMMs), 1616 LTR retrotransposon sequences containing at least two models were screened. A total of 48 RBIP primers were synthesized based on the high copy numbers of conserved LTR sequences. Fifty-six amplicons with an average polymorphism of 91.07% were generated in 105 sweet potato germplasm resources based on RBIP markers. A Unweighted Pair Group Method with Arithmatic Mean (UPGMA) dendrogram, a model-based genetic structure and principal component analysis divided the sweet potato germplasms into 3 groups containing 8, 53, and 44 germplasms. All the three analyses produced significant groupwise consensus. However, almost all the germplasms contained only one primary locus. The analysis of molecular variance (AMOVA) among the groups indicated higher intergroup genetic variation (53%) than intrapopulation genetic variation. In addition, long-term self-retention may cause some germplasm resources to exhibit variable segregation. These results suggest that these sweet potato germplasms are not well evolutionarily diversified, although geographic speciation could have occurred at a limited level. This study highlights the utility of RBIP markers for determining the intraspecies variability of sweet potato and have the potential to be used as core primer pairs for variety identification, genetic diversity assessment and linkage map construction. The results could provide a good theoretical reference and guidance for germplasm research and breeding.

Identification of QTL for resistance to root rot in sweetpotato (Ipomoea batatas (L.) Lam) with SSR linkage maps

BACKGROUND

Sweetpotato root rot is a devastating disease caused by Fusarium solani that seriously endangers the yield of sweetpotato in China. Although there is currently no effective method to control the disease, breeding of resistant varieties is the most effective and economic option. Moreover, quantitative trait locus (QTL) associated with resistance to root rot have not yet been reported, and the biological mechanisms of resistance remain unclear in sweetpotato. Thus, increasing our knowledge about the mechanism of disease resistance and identifying resistance loci will assist in the development of disease resistance breeding.

RESULTS

In this study, we constructed genetic linkage maps of sweetpotato using a mapping population consisting of 300 individuals derived from a cross between Jizishu 1 and Longshu 9 by simple sequence repeat (SSR) markers, and mapped seven QTLs for resistance to root rot. In total, 484 and 573 polymorphic SSR markers were grouped into 90 linkage groups for Jizishu 1 and Longshu 9, respectively. The total map distance for Jizishu 1 was 3974.24 cM, with an average marker distance of 8.23 cM. The total map distance for Longshu 9 was 5163.35 cM, with an average marker distance of 9.01 cM. Five QTLs (qRRM_1, qRRM_2, qRRM_3, qRRM_4, and qRRM_5) were located in five linkage groups of Jizishu 1 map explaining 52.6-57.0% of the variation. Two QTLs (qRRF_1 and qRRF_2) were mapped on two linkage groups of Longshu 9 explaining 57.6 and 53.6% of the variation, respectively. Furthermore, 71.4% of the QTLs positively affected the variation. Three of the seven QTLs, qRRM_3, qRRF_1, and qRRF_2, were colocalized with markers IES43-5mt, IES68-6 fs**, and IES108-1 fs, respectively.

CONCLUSIONS

To our knowledge, this is the first report on the construction of a genetic linkage map for purple sweetpotato (Jizishu 1) and the identification of QTLs associated with resistance to root rot in sweetpotato using SSR markers. These QTLs will have practical significance for the fine mapping of root rot resistance genes and play an important role in sweetpotato marker-assisted breeding.

Multiple QTL Mapping in Autopolyploids: A Random-Effect Model Approach with Application in a Hexaploid Sweetpotato Full-Sib Population

In developing countries, the sweetpotato, Ipomoea batatas (L.) Lam. [Formula: see text], is an important autopolyploid species, both socially and economically. However, quantitative trait loci (QTL) mapping has remained limited due to its genetic complexity. Current fixed-effect models can fit only a single QTL and are generally hard to interpret. Here, we report the use of a random-effect model approach to map multiple QTL based on score statistics in a sweetpotato biparental population ('Beauregard' × 'Tanzania') with 315 full-sibs. Phenotypic data were collected for eight yield component traits in six environments in Peru, and jointly adjusted means were obtained using mixed-effect models. An integrated linkage map consisting of 30,684 markers distributed along 15 linkage groups (LGs) was used to obtain the genotype conditional probabilities of putative QTL at every centiMorgan position. Multiple interval mapping was performed using our R package QTLpoly and detected a total of 13 QTL, ranging from none to four QTL per trait, which explained up to 55% of the total variance. Some regions, such as those on LGs 3 and 15, were consistently detected among root number and yield traits, and provided a basis for candidate gene search. In addition, some QTL were found to affect commercial and noncommercial root traits distinctly. Further best linear unbiased predictions were decomposed into additive allele effects and were used to compute multiple QTL-based breeding values for selection. Together with quantitative genotyping and its appropriate usage in linkage analyses, this QTL mapping methodology will facilitate the use of genomic tools in sweetpotato breeding as well as in other autopolyploids.

Development of molecular markers associated with resistance to Meloidogyne incognita by performing quantitative trait locus analysis and genome-wide association study in sweetpotato

The southern root-knot nematode, Meloidogyne incognita, is a pest that decreases yield and the quality of sweetpotato [Ipomoea batatas (L.) Lam.]. There is a demand to produce resistant cultivars and develop DNA markers to select this trait. However, sweetpotato is hexaploid, highly heterozygous, and has an enormous genome (∼3 Gb), which makes genetic linkage analysis difficult. In this study, a high-density linkage map was constructed based on retrotransposon insertion polymorphism, simple sequence repeat, and single nucleotide polymorphism markers. The markers were developed using F₁ progeny between J-Red, which exhibits resistance to multiple races of M. incognita, and Choshu, which is susceptible to multiple races of such pest. Quantitative trait locus (QTL) analysis and a genome-wide association study detected highly effective QTLs for resistance against three races, namely, SP1, SP4, and SP6-1, in the Ib01-6 J-Red linkage group. A polymerase chain reaction marker that can identify genotypes based on single nucleotide polymorphisms located in this QTL region can discriminate resistance from susceptibility in the F₁ progeny at a rate of 70%. Thus, this marker could be helpful in selecting sweetpotato cultivars that are resistant to multiple races of M. incognita.

Developing Transposable Element Marker System for Molecular Breeding

Transposable element (TE) marker system was developed considering the useful properties of the transposable elements such as their large number in the animal and plant genomes, high rate of insertion polymorphism, and ease of detection. Various methods have been employed for developing a large number of TE markers in several crop plants for genomics studies. Here we describe some of these methods including the recent whole genome search. We also review the application of TE markers in molecular breeding.

Genome-Wide Association Studies (GWAS) for Yield and Weevil Resistance in Sweet potato (Ipomoea batatas (L.) Lam)

We apply the GWAS to sweet potato genome, and identified the SNPs associated with yield and weevil resistance. The sweet potato (Ipomoea batatas (L.) Lam) is a highly heterozygous, outcrossing, polyploid species, which presents challenges for genetic analysis. Therefore, we considered that genome-wide association studies (GWAS) may be applied to the study of the sweet potato genome. The yield of two sweet potato varieties [Purple Sweet Lord (PSL) and 90IDN-47] was assessed at two locations (Kumamoto and Okinawa prefectures) in Japan in 2013 and the yield scores were used for GWAS. The results showed that there were several single nucleotide polymorphisms (SNP) above the significance thresholds in PSL; two peaks were detected in Kumamoto and Okinawa on the Ib03-3 and Ib01-4 linkage groups of PSL, respectively. As for 90IDN-47, one relatively high peak was detected in Kumamoto on the Ib13-8 linkage group. Interestingly, although high peaks above significance thresholds were detected in Kumamoto and Okinawa in PSL, the peaks were located in different linkage groups. This result suggests that the genetic regions controlling yield may change in response to environmental conditions. Additionally, we investigated the degree of weevil damage to the plants, which is the greatest problem in sweet potato cultivation in Okinawa. In this experiment, no SNPs were identified above the significance thresholds. However, one relatively high peak was found in the 90IDN-47 genotype, which showed resistance to weevils. On the other hand, one relatively high peak was also detected in the PSL genotype, which showed susceptibility to weevils. These results suggest that two regions could affect weevil resistance and may contain the gene(s) controlling weevil resistance.

Genome-wide association study overcomes the genome complexity in autohexaploid chrysanthemum and tags SNP markers onto the flower color genes

The use of DNA markers has revolutionized selection in crop breeding by linkage mapping and QTL analysis, but major problems still remain for polyploid species where marker-assisted selection lags behind the situation in diploids because of its high genome complexity. To overcome the complex genetic mode in the polyploids, we investigated the development of a strategy of genome-wide association study (GWAS) using single-dose SNPs, which simplify the segregation patterns associated polyploids, with respect to the development of DNA markers. In addition, we employed biparental populations for the GWAS, wherein the SNP allele frequency could be predicted. The research investigated whether the method could be used to effectively develop DNA markers for petal color in autohexaploid chrysanthemum (Chrysanthemum morifolium; 2n = 6x = 54). The causal gene for this trait is already-known CmCCD4a encoding a dioxygenase which cleaves carotenoids in petals. We selected 9,219 single-dose SNPs, out of total 52,489 SNPs identified by dd-RAD-Seq, showing simplex (1 × 0) and double-simplex (1 × 1) inheritance pattern according to alternative allele frequency with respect to the SNP loci in the F1 population. GWAS, using these single-dose SNPs, discovered highly reproducible SNP markers tightly linked to the causal genes. This is the first report of a straightforward GWAS-based marker developing system for use in autohexaploid species.

An ultra-dense integrated linkage map for hexaploid chrysanthemum enables multi-allelic QTL analysis

KEY MESSAGE

We constructed the first integrated genetic linkage map in a polysomic hexaploid. This enabled us to estimate inheritance of parental haplotypes in the offspring and detect multi-allelic QTL. Construction and use of linkage maps are challenging in hexaploids with polysomic inheritance. Full map integration requires calculations of recombination frequency between markers with complex segregation types. In addition, detection of QTL in hexaploids requires information on all six alleles at one locus for each individual. We describe a method that we used to construct a fully integrated linkage map for chrysanthemum (Chrysanthemum × morifolium, 2n = 6x = 54). A bi-parental F1 population of 406 individuals was genotyped with an 183,000 SNP genotyping array. The resulting linkage map consisted of 30,312 segregating SNP markers of all possible marker dosage types, representing nine chromosomal linkage groups and 107 out of 108 expected homologues. Synteny with lettuce (Lactuca sativa) showed local colinearity. Overall, it was high enough to number the chrysanthemum chromosomal linkage groups according to those in lettuce. We used the integrated and phased linkage map to reconstruct inheritance of parental haplotypes in the F1 population. Estimated probabilities for the parental haplotypes were used for multi-allelic QTL analyses on four traits with different underlying genetic architectures. This resulted in the identification of major QTL that were affected by multiple alleles having a differential effect on the phenotype. The presented linkage map sets a standard for future genetic mapping analyses in chrysanthemum and closely related species. Moreover, the described methods are a major step forward for linkage mapping and QTL analysis in hexaploids.

AlloMap6: an R package for genetic linkage analysis in allohexaploids

Allopolyploids are a group of polyploids with more than two sets of chromosomes derived from different species. Previous linkage analysis of allopolyploids is based on the assumption that different chromosomes pair randomly during meiosis. A more sophisticated model to relax this assumption has been developed for allotetraploids by incorporating the preferential pairing behavior of homologous over homoeologous chromosomes. Here, we show that the basic principle of this model can be extended to perform linkage analysis of higher-ploidy allohexaploids, where multiple preferential pairing factors are used to characterize chromosomal-pairing meiotic features between different constituent species. We implemented the extended model into an R package, called AlloMap6, allowing the recombination fractions and preferential pairing factors to be estimated simultaneously. Allomap6 has two major functionalities, computer simulation and real-data analysis. By analyzing a real data from a full-sib family of allohexaploid persimmon, we tested and validated the usefulness and utility of this package. AlloMap6 lays a foundation for allohexaploid genetic mapping and provides a new horizon to explore the chromosomal kinship of allohexaploids.

A high-density SNP genetic map consisting of a complete set of homologous groups in autohexaploid sweetpotato (Ipomoea batatas)

Sweetpotato (Ipomoea batatas) is an autohexaploid species with 90 chromosomes (2n = 6x = 90) and a basic chromosome number of 15, and is therefore regarded as one of the most challenging species for high-density genetic map construction. Here, we used single nucleotide polymorphisms (SNPs) identified by double-digest restriction site-associated DNA sequencing based on next-generation sequencing technology to construct a map for sweetpotato. We then aligned the sequence reads onto the reference genome sequence of I. trifida, a likely diploid ancestor of sweetpotato, to detect SNPs. In addition, to simplify analysis of the complex genetic mode of autohexaploidy, we used an S1 mapping population derived from self-pollination of a single parent. As a result, 28,087 double-simplex SNPs showing a Mendelian segregation ratio in the S1 progeny could be mapped onto 96 linkage groups (LGs), covering a total distance of 33,020.4 cM. Based on the positions of the SNPs on the I. trifida genome, the LGs were classified into 15 groups, each with roughly six LGs and six small extra groups. The molecular genetic techniques used in this study are applicable to high-density mapping of other polyploid plant species, including important crops.

Challenges to genome sequence dissection in sweetpotato

The development of next generation sequencing (NGS) technologies has enabled the determination of whole genome sequences in many non-model plant species. However, genome sequencing in sweetpotato (Ipomoea batatas (L.) Lam) is still difficult because of the hexaploid genome structure. Previous studies suggested that a diploid wild relative, I. trifida (H.B.K.) Don., is the most possible ancestor of sweetpotato. Therefore, the genetic and genomic features of I. trifida have been studied as a potential reference for sweetpotato. Meanwhile, several research groups have begun the challenging task of directly sequencing the sweetpotato genome. In this manuscript, we review the recent results and activities of large-scale genome and transcriptome analysis related to genome sequence dissection in sweetpotato under the sections as follows: I. trifida genome and transcript sequencing, genome sequences of I. nil (Japanese morning glory), transcript sequences in sweetpotato, chloroplast sequences, transposable elements and transfer DNA. The recent international activities of de novo whole genome sequencing in sweetpotato are also described. The large-scale publically available genome and transcript sequence resources and the international genome sequencing streams are expected to promote the genome sequence dissection in sweetpotato.

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.

Inheritance of low pasting temperature in sweetpotato starch and the dosage effect of wild-type alleles

Sweetpotato (Ipomoea batatas (L.) Lam.), which is an outcrossing hexaploid, is one of the most important starch-producing crops in the world. During the last decade, new sweetpotato cultivars, e.g. 'Quick Sweet', which have approximately 20°C lower pasting temperature, slower retrogradation and higher digestibility of raw starch than ordinary cultivars, have been developed in Japan. Genetic analysis of these variants with low pasting temperature starch was conducted in this study. Using 8 variants and 15 normal clones, 26 families were generated. The results from analyzing these progenies suggested that this trait is a qualitative character controlled by one recessive allele (designated spt), which is inherited in a hexasomic manner. A dosage effect of the wild-type Spt allele was found for starch pasting temperature, although the effect was not linear. These results will aid breeders to develop sweetpotato cultivars with a range of starch pasting temperatures.