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

A locus-dependent mixed inheritance in the segmental allohexaploid sweetpotato (Ipomoea batatas [L.] Lam)

Two interrelated aspects of the sweetpotato genome, its polyploid origin and inheritance type, remain uncertain. We recently proposed a segmental allohexaploid sweetpotato and thus sought to clarify its inheritance type by direct analyses of homoeolog segregations at selected single-copy loci. For such analyses, we developed a digital quantitative PCR genotyping method using one nondiscriminatory and three discriminatory probes for each selected locus to discriminate and quantify three homoeolog-differentiating variation types (homoeolog-types) in genomic DNA samples for genotype fitting and constructed a F2 population for segregation analyses. We confirmed inter-subgenomic distinctions of three identified homoeolog-types at each of five selected loci by their interspecific differentiations among 14 species in Ipomoea section batatas and genotyped the loci in 549 F2 lines, selected F1 progenies, and their founding parents. Segregation and genotype analyses revealed a locus-dependent mixed inheritance (disomic, polysomic, and intermediate types) of the homoeolog-types at 4 loci in the F2 population, displaying estimated disomic-inheritance frequencies of 0, 2.72%, 14.52%, and 36.92%, and probably in the F1 population too. There were also low-frequency non-hexaploid F1 and F2 genotypes that were probably derived from double-reduction recombination or partially unreduced gametes, and F2 genotypes of apparent aneuploids/dysploids with neopolyploid-like frequencies. Additional analyses of homoeolog-type genotypes at the 5 loci in 46 lines from various regions revealed locus-dependent selection biases, favoring genotypes having more of one homoeolog-type, i.e. more of di- or homogenized homoeolog-type composition, and one-direction ploidy trending among apparent aneuploids/dysploids. These inheritance features pointed to an evolving segmental allohexaploid sweetpotato impacted by selection biases.

A self-compatible pear mutant derived from γ-irradiated pollen carries an 11-Mb duplication in chromosome 17

Self-compatibility is a highly desirable trait for pear breeding programs. Our breeding program previously developed a novel self-compatible pollen-part Japanese pear mutant (Pyrus pyrifolia Nakai), '415-1', by using γ-irradiated pollen. '415-1' carries the S-genotype S4dS5S5, with "d" indicating a duplication of S 5 responsible for breakdown of self-incompatibility. Until now, the size and inheritance of the duplicated segment was undetermined, and a reliable detection method was lacking. Here, we examined genome duplications and their inheritance in 140 F1 seedlings resulting from a cross between '515-20' (S1S3) and '415-1'. Amplicon sequencing of S-RNase and SFBB18 clearly detected S-haplotype duplications in the seedlings. Intriguingly, 30 partially triploid seedlings including genotypes S1S4dS5, S3S4dS5, S1S5dS5, S3S5dS5, and S3S4dS4 were detected among the 140 seedlings. Depth-of-coverage analysis using ddRAD-seq showed that the duplications in those individuals were limited to chromosome 17. Further analysis through resequencing confirmed an 11-Mb chromosome duplication spanning the middle to the end of chromosome 17. The duplicated segment remained consistent in size across generations. The presence of an S3S4dS4 seedling provided evidence for recombination between the duplicated S5 segment and the original S4haplotype, suggesting that the duplicated segment can pair with other parts of chromosome 17. This research provides valuable insights for improving pear breeding programs using partially triploid individuals.

Construction of relatedness matrices in autopolyploid populations using low-depth high-throughput sequencing data

KEY MESSAGE

An improved estimator of genomic relatedness using low-depth high-throughput sequencing data for autopolyploids is developed. Its outputs strongly correlate with SNP array-based estimates and are available in the package GUSrelate. High-throughput sequencing (HTS) methods have reduced sequencing costs and resources compared to array-based tools, facilitating the investigation of many non-model polyploid species. One important quantity that can be computed from HTS data is the genetic relatedness between all individuals in a population. However, HTS data are often messy, with multiple sources of errors (i.e. sequencing errors or missing parental alleles) which, if not accounted for, can lead to bias in genomic relatedness estimates. We derive a new estimator for constructing a genomic relationship matrix (GRM) from HTS data for autopolyploid species that accounts for errors associated with low sequencing depths, implemented in the R package GUSrelate. Simulations revealed that GUSrelate performed similarly to existing GRM methods at high depth but reduced bias in self-relatedness estimates when the sequencing depth was low. Using a panel consisting of 351 tetraploid potato genotypes, we found that GUSrelate produced GRMs from genotyping-by-sequencing (GBS) data that were highly correlated with a GRM computed from SNP array data, and less biased than existing methods when benchmarking against the array-based GRM estimates. GUSrelate provides researchers with a tool to reliably construct GRMs from low-depth HTS data.

Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids

Discovery and analysis of genetic variants underlying agriculturally important traits are key to molecular breeding of crops. Reduced representation approaches have provided cost-efficient genotyping using next-generation sequencing. However, accurate genotype calling from next-generation sequencing data is challenging, particularly in polyploid species due to their genome complexity. Recently developed Bayesian statistical methods implemented in available software packages, polyRAD, EBG, and updog, incorporate error rates and population parameters to accurately estimate allelic dosage across any ploidy. We used empirical and simulated data to evaluate the three Bayesian algorithms and demonstrated their impact on the power of genome-wide association study (GWAS) analysis and the accuracy of genomic prediction. We further incorporated uncertainty in allelic dosage estimation by testing continuous genotype calls and comparing their performance to discrete genotypes in GWAS and genomic prediction. We tested the genotype-calling methods using data from two autotetraploid species, Miscanthus sacchariflorus and Vaccinium corymbosum, and performed GWAS and genomic prediction. In the empirical study, the tested Bayesian genotype-calling algorithms differed in their downstream effects on GWAS and genomic prediction, with some showing advantages over others. Through subsequent simulation studies, we observed that at low read depth, polyRAD was advantageous in its effect on GWAS power and limit of false positives. Additionally, we found that continuous genotypes increased the accuracy of genomic prediction, by reducing genotyping error, particularly at low sequencing depth. Our results indicate that by using the Bayesian algorithm implemented in polyRAD and continuous genotypes, we can accurately and cost-efficiently implement GWAS and genomic prediction in polyploid crops.

Bayesian tests for random mating in polyploids

Hardy-Weinberg proportions (HWP) are often explored to evaluate the assumption of random mating. However, in autopolyploids, organisms with more than two sets of homologous chromosomes, HWP and random mating are different hypotheses that require different statistical testing approaches. Currently, the only available methods to test for random mating in autopolyploids (i) heavily rely on asymptotic approximations and (ii) assume genotypes are known, ignoring genotype uncertainty. Furthermore, these approaches are all frequentist, and so do not carry the benefits of Bayesian analysis, including ease of interpretability, incorporation of prior information, and consistency under the null. Here, we present Bayesian approaches to test for random mating, bringing the benefits of Bayesian analysis to this problem. Our Bayesian methods also (i) do not rely on asymptotic approximations, being appropriate for small sample sizes, and (ii) optionally account for genotype uncertainty via genotype likelihoods. We validate our methods in simulations and demonstrate on two real datasets how testing for random mating is more useful for detecting genotyping errors than testing for HWP (in a natural population) and testing for Mendelian segregation (in an experimental S1 population). Our methods are implemented in Version 2.0.2 of the hwep R package on the Comprehensive R Archive Network https://cran.r-project.org/package=hwep.

Double reduction estimation and equilibrium tests in natural autopolyploid populations

Many bioinformatics pipelines include tests for equilibrium. Tests for diploids are well studied and widely available, but extending these approaches to autopolyploids is hampered by the presence of double reduction, the comigration of sister chromatid segments into the same gamete during meiosis. Though a hindrance for equilibrium tests, double reduction rates are quantities of interest in their own right, as they provide insights about the meiotic behavior of autopolyploid organisms. Here, we develop procedures to (i) test for equilibrium while accounting for double reduction, and (ii) estimate the double reduction rate given equilibrium. To do so, we take two approaches: a likelihood approach, and a novel U-statistic minimization approach that we show generalizes the classical equilibrium χ2 test in diploids. For small sample sizes and uncertain genotypes, we further develop a bootstrap procedure based on our U-statistic to test for equilibrium. We validate our methods on both simulated and real data.

Polyploid GWAS reveals the basis of molecular marker development for complex breeding traits including starch content in the storage roots of sweet potato

Given the importance of prioritizing genome-based breeding of sweet potato to enable the promotion of food and nutritional security for future human societies, here, we aimed to dissect the genetic basis of storage root starch content (SC) when associated with a complex set of breeding traits including dry matter (DM) rate, storage root fresh weight (SRFW), and anthocyanin (AN) content in a mapping population containing purple-fleshed sweet potato. A polyploid genome-wide association study (GWAS) was extensively exploited using 90,222 single-nucleotide polymorphisms (SNPs) obtained from a bi-parental 204 F₁ population between 'Konaishin' (having high SC but no AN) and 'Akemurasaki' (having high AN content but moderate SC). Through the comparison of polyploid GWAS on the whole set of the 204 F₁, 93 high-AN-containing F₁, and 111 low-AN-containing F₁ populations, a total of two (consists of six SNPs), two (14 SNPs), four (eight SNPs), and nine (214 SNPs) significantly associated signals were identified for the variations of SC, DM, SRFW, and the relative AN content, respectively. Of them, a novel signal associated with SC, which was most consistent in 2019 and 2020 in both the 204 F₁ and 111 low-AN-containing F₁ populations, was identified in homologous group 15. The five SNP markers associated with homologous group 15 could affect SC improvement with a degree of positive effect (~4.33) and screen high-starch-containing lines with higher efficiency (~68%). In a database search of 62 genes involved in starch metabolism, five genes including enzyme genes granule-bound starch synthase I (IbGBSSI), α-amylase 1D, α-amylase 1E, and α-amylase 3, and one transporter gene ATP/ADP-transporter were located on homologous group 15. In an extensive qRT-PCR of these genes using the storage roots harvested at 2, 3, and 4 months after field transplantation in 2022, IbGBSSI, which encodes the starch synthase isozyme that catalyzes the biosynthesis of amylose molecule, was most consistently elevated during starch accumulation in sweet potato. These results would enhance our understanding of the underlying genetic basis of a complex set of breeding traits in the starchy roots of sweet potato, and the molecular information, particularly for SC, would be a potential platform for molecular marker development for this trait.

A genome-wide association and fine-mapping study of white rust resistance in hexaploid chrysanthemum cultivars with a wild diploid reference genome

White rust caused by Puccinia horiana is one of the most serious diseases of chrysanthemum (Chrysanthemum × morifolium). In this study, we report the DNA markers associated with resistance against P. horiana via a simple approach using the genome of a wild diploid relative, Chrysanthemum seticuspe. First, we identified the important region of the genome in the resistant cultivar "Ariesu" via a genome-wide association study. Simplex single nucleotide polymorphism (SNP) markers mined from ddRAD-Seq were used in a biparental population originating from crosses between resistant "Ariesu" and susceptible "Yellow Queen". The C. seticuspe genome was used as a reference. For the fine mapping of P. horiana resistance locus 2 (Phr2), a comparative whole genome sequencing study was conducted. Although the genome sequences of chrysanthemum cultivars assembled via the short-read approach were fragmented, reliable genome alignments were reconstructed by mapping onto the chromosome level of the C. seticuspe pseudomolecule. Base variants were then identified by comparing the assembled genome sequences of resistant "Ariesu" and susceptible "Yellow Queen". Consequently, SNP markers that were closer to Phr2 compared with ddRAD-Seq markers were obtained. These SNP markers co-segregated with resistance in F₁ progenies originating from resistant "Ariesu" and showed robust transferability for detecting Phr2-conferring resistance among chrysanthemum genetic resources. The wild C. seticuspe pseudomolecule, a de facto monoploid genome used for ddRAD-Seq analysis and assembled genome sequence comparison, demonstrated this method's utility as a model for developing DNA markers in hexaploid chrysanthemum cultivars.

Genetic mapping of loci involved in oil tocopherol composition control in Russian sunflower (Helianthus annuus L.) lines

Tocopherols are antioxidants that preserve oil lipids against oxidation and serve as a natural source of vitamin E in the human diet. Compared with other major oilseeds like rapeseed and soybean, sunflower (Helianthus annuus L.) exhibits low phenotypic diversity of tocopherol composition, both in wild and cultivated accessions from germplasm collections. Two major mutations that alter tocopherol composition were identified in genetic collections, and several studies suggested additional loci controlling tocopherol composition, with their expression possibly depending on the genetic background. In the present study, we performed QTL mapping of tocopherol composition in two independent F2 crosses between lines with contrasting tocopherol composition from the Pustovoit All-Russia Research Institute of Oil Crops (VNIIMK) collection. We used genotyping-bysequencing (GBS) to construct single nucleotide polymorphism-based genetic maps, and performed QTL mapping using quantitative and qualitative encoding for phenotypic traits. Our results support the notion that the tocopherol composition in the assessed crosses is controlled by two loci. We additionally selected and validated two single nucleotide polymorphism markers for each cross which could be used for marker-assisted selection.

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.

High-Density Single Nucleotide Polymorphisms Genetic Map Construction and Quantitative Trait Locus Mapping of Color-Related Traits of Purple Sweet Potato [Ipomoea batatas (L.) Lam.]

Flesh color (FC), skin color (SC), and anthocyanin content (AC) are three important traits being used for commodity evaluation in purple-fleshed sweet potato. However, to date, only a few reports are available on the inheritance of these traits. In this study, we used a biparental mapping population of 274 F1 progeny generated from a cross between a dark purple-fleshed (Xuzishu8) and white-fleshed (Meiguohong) sweet potato variety for genetic analyses. Correlation analysis showed a significant positive correlation among AC, SC, and FC. Medium-to-high heritability was observed for these traits. We detected single nucleotide polymorphisms (SNPs) by specific length amplified fragment sequencing (SLAF-seq) with the average sequencing depth of 51.72 and 25.76 for parents and progeny, respectively. Then we constructed an integrated genetic map consisting of 15 linkage groups (LGS) of sweet potato spanning on 2,233.66 cm with an average map distance of 0.71 cm between adjacent markers. Based on the linkage map, ten major quantitative trait loci (QTLs) associated to FC, SC, and AC were identified on LG12 between 0 and 64.97 cm distance, such as one QTL for SC and FC, respectively, which explained 36.3 and 45.9% of phenotypic variation; eight QTLs for AC, which explained 10.5-28.5% of the variation. These major QTLs were highly consistent and co-localized on LG12. Positive correlation, high heritability, and co-localization of QTLs on the same LG group confirm the significance of this study to establish a marker-assisted breeding program for sweet potato improvement.

Meiosis in Polyploids and Implications for Genetic Mapping: A Review

Plant cytogenetic studies have provided essential knowledge on chromosome behavior during meiosis, contributing to our understanding of this complex process. In this review, we describe in detail the meiotic process in auto- and allopolyploids from the onset of prophase I through pairing, recombination, and bivalent formation, highlighting recent findings on the genetic control and mode of action of specific proteins that lead to diploid-like meiosis behavior in polyploid species. During the meiosis of newly formed polyploids, related chromosomes (homologous in autopolyploids; homologous and homoeologous in allopolyploids) can combine in complex structures called multivalents. These structures occur when multiple chromosomes simultaneously pair, synapse, and recombine. We discuss the effectiveness of crossover frequency in preventing multivalent formation and favoring regular meiosis. Homoeologous recombination in particular can generate new gene (locus) combinations and phenotypes, but it may destabilize the karyotype and lead to aberrant meiotic behavior, reducing fertility. In crop species, understanding the factors that control pairing and recombination has the potential to provide plant breeders with resources to make fuller use of available chromosome variations in number and structure. We focused on wheat and oilseed rape, since there is an abundance of elucidating studies on this subject, including the molecular characterization of the Ph1 (wheat) and PrBn (oilseed rape) loci, which are known to play a crucial role in regulating meiosis. Finally, we exploited the consequences of chromosome pairing and recombination for genetic map construction in polyploids, highlighting two case studies of complex genomes: (i) modern sugarcane, which has a man-made genome harboring two subgenomes with some recombinant chromosomes; and (ii) hexaploid sweet potato, a naturally occurring polyploid. The recent inclusion of allelic dosage information has improved linkage estimation in polyploids, allowing multilocus genetic maps to be constructed.

Genome-wide study on the polysomic genetic factors conferring plasticity of flower sexuality in hexaploid persimmon

Sexuality is one of the fundamental mechanisms that work towards maintaining genetic diversity within a species. In diploid persimmons (Diospyros spp.), separated sexuality, the presence of separate male and female individuals (dioecy), is controlled by the Y chromosome-encoded small-RNA gene, OGI. On the other hand, sexuality in hexaploid Oriental persimmon (Diospyros kaki) is more plastic, with OGI-bearing genetically male individuals, able to produce both male and female flowers (monoecy). This is thought to be linked to the partial inactivation of OGI by a retrotransposon insertion, resulting in DNA methylation of the OGI promoter region. To identify the genetic factors regulating branch sexual conversion, genome-wide correlation/association analyses were conducted using ddRAD-Seq data from an F₁ segregating population, and using both quantitative and diploidized genotypes, respectively. We found that allelic ratio at the Y-chromosomal region, including OGI, was correlated with male conversion based on quantitative genotypes, suggesting that OGI can be activated in cis in a dosage-dependent manner. Genome-wide association analysis based on diploidized genotypes, normalized for the effect of OGI allele dosage, detected three fundamental loci associated with male conversion. These loci underlie candidate genes, which could potentially act epigenetically for the activation of OGI expression.

DNA marker for resistance to Puccinia horiana in chrysanthemum (Chrysanthemum morifolium Ramat.) "Southern Pegasus"

White rust caused by Puccinia horiana Henn. adversely affects chrysanthemum (Chrysanthemum morifolium Ramat.) production. The breeding of resistant varieties is effective in controlling the disease. Here we aimed to develop DNA markers for the strong resistance to P. horiana. We conducted a linkage analysis based on the genome-wide association study (GWAS) method. We employed a biparental population for the GWAS, wherein the single nucleotide polymorphism (SNP) allele frequency could be predicted. The population was derived from crosses between a strong resistant "Southern Pegasus" and a susceptible line. The GWAS used simplex and double-simplex SNP markers selected out of SNP candidates mined from ddRAD-Seq data of an F₁ biparental population. These F₁ individuals segregated in a 1:1 ratio of resistant to susceptible. Twenty-one simplex SNPs were significantly associated with P. horiana resistance in "Southern Pegasus" and generated one linkage group. These results show the presence of a single resistance gene in "Southern Pegasus". We identified the nearest SNP marker located 2.2 cM from P. horiana resistance locus and demonstrated this SNP marker-resistance link using an independent population. This is the first report of an effective DNA marker linked to a gene for P. horiana resistance in chrysanthemum.

Identification of genome-wide single-nucleotide polymorphisms among geographically diverse radish accessions

Radish (Raphanus sativus L.) is cultivated around the world as a vegetable crop and exhibits diverse morphological and physiological features. DNA polymorphisms are responsible for differences in traits among cultivars. In this study, we determined genome-wide single-nucleotide polymorphisms (SNPs) among geographically diverse radish accessions using the double-digest restriction site-associated DNA sequencing (ddRAD-Seq) method. A total of 52,559 SNPs was identified in a collection of over 500 radish accessions (cultivated and wild) from East Asia, South and Southeast Asia, and the Occident and Near East. In addition, 2,624 SNP sites without missing data (referred to as common SNP sites) were identified among 510 accessions. Genetic diversity analyses, based on the common SNP sites, divided the cultivated radish accessions into four main groups, each derived from four geographical areas (Japan, East Asia, South and Southeast Asia, and the Occident and Near East). Furthermore, we discuss the origin of cultivated radish and its migration from the West to East Asia. SNP data generated in this work will facilitate further genetic studies on the radish breeding and production of DNA markers.

Genetic Mapping in Autohexaploid Sweet Potato with Low-Coverage NGS-Based Genotyping Data

Next-generation sequencing (NGS)-based genotyping methods can generate numerous genetic markers in a single experiment and have contributed to plant genetic mapping. However, for high precision genetic analysis, the complicated genetic segregation mode in polyploid organisms requires high-coverage NGS data and elaborate analytical algorithms. In the present study, we propose a simple strategy for the genetic mapping of polyploids using low-coverage NGS data. The validity of the strategy was investigated using simulated data. Previous studies indicated that accurate allele dosage estimation from low-coverage NGS data (read depth < 40) is difficult. Therefore, we used allele dosage probabilities calculated from read counts in association analyses to detect loci associated with phenotypic variations. The allele dosage probabilities showed significant detection power, although higher allele dosage estimation accuracy resulted in higher detection power. On the contrary, differences in the segregation patterns between the marker and causal genes resulted in a drastic decrease in detection power even if the marker and casual genes were in complete linkage and the allele dosage estimation was accurate. These results indicated that the use of a larger number of markers is advantageous, even if the accuracy of allele dosage estimation is low. Finally, we applied the strategy for the genetic mapping of autohexaploid sweet potato (Ipomoea batatas) populations to detect loci associated with agronomic traits. Our strategy could constitute a cost-effective approach for preliminary experiments done performed to large-scale studies.

ddRAD-seq derived genome-wide SNPs, high density linkage map and QTLs for fruit quality traits in strawberry (Fragaria x ananassa)

Understanding the genetic determinants are essential for improving the fruit quality traits of strawberry. In this study, we focused on mapping the loci for fruit-length (FL), -diameter (FD), -weight (FW) and -soluble solid content (SSC) using the genome-wide single nucleotide polymorphisms (SNPs) identified via ddRAD-sequencing of the F₁ population raised from Maehyang (♀) X Festival (♂). A total of 12,698 high quality SNPs were identified of which 1554 SNPs that showed significant Mendelian segregation (p < 0.05) were mapped to 53 linkage groups (LG) spanning a total of 2937.93 cM with an average marker density of 2.14 cM/locus. Six QTLs for FL and four QTLs for each of FD, FW and SSC were identified that explained 24–35%, 21–42%, 24–54% and 23–50% of overall phenotypic variations, respectively. The genes that lie within these QTL regions were extracted and discussed thoroughly. In addition, a high resolution melting marker (MF154) were designed based on the SNP A1⁷²³G of the UDP-glucose 4-epimerase GEPI48-like gene FAN_iscf00021287. The marker detected the high vs low sugar containing F₁ plants and commercial cultivars with 81.39% and 86.95% detection accuracy, respectively. These SNPs, linkage map, QTLs and candidate genes will be helpful in understanding and improving the fruit quality traits of strawberry.

QTL analysis and GWAS of agronomic traits in sweetpotato (Ipomoea batatas L.) using genome wide SNPs

While sweetpotato (Ipomoea batatas L.) improvement has generally been done by field-based selection, molecular genetic studies on traits of interest, i.e., molecular markers are needed for enhancing the breeding program of this world's 7^th most important crop, as such markers facilitate marker-assisted selection. Here, we performed a combined approach of QTLs analyses and GWAS of storage root β-carotene content (BC), dry-matter (DM) and starch content (SC) using the genetic linkage maps constructed with 5,952 and 5,640 SNPs obtained from F₁ progenies between cultivars 'J-Red' and 'Choshu'. BC was negatively correlated with DM (r = -0.45) and SC (r = -0.51), while DM was positively correlated with SC (r = 0.94). In both parental maps, a total of five, two and five QTL regions on linkage groups 7 and 8 were associated with BC, DM and SC, respectively. In GWAS of BC, one strong signal (P = 1.04 × 10^-9) was observed on linkage group 8, which co-located with one of the above QTL regions. The SNPs markers found here, particularly for β-carotene, would be useful base resources for future marker-assisted selection program with this trait.

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.

Priors for genotyping polyploids

MOTIVATION

Empirical Bayes techniques to genotype polyploid organisms usually either (i) assume technical artifacts are known a priori or (ii) estimate technical artifacts simultaneously with the prior genotype distribution. Case (i) is unappealing as it places the onus on the researcher to estimate these artifacts, or to ensure that there are no systematic biases in the data. However, as we demonstrate with a few empirical examples, case (ii) makes choosing the class of prior genotype distributions extremely important. Choosing a class is either too flexible or too restrictive results in poor genotyping performance.

RESULTS

We propose two classes of prior genotype distributions that are of intermediate levels of flexibility: the class of proportional normal distributions and the class of unimodal distributions. We provide a complete characterization of and optimization details for the class of unimodal distributions. We demonstrate, using both simulated and real data that using these classes results in superior genotyping performance.

AVAILABILITY AND IMPLEMENTATION

Genotyping methods that use these priors are implemented in the updog R package available on the Comprehensive R Archive Network: https://cran.r-project.org/package=updog. All code needed to reproduce the results of this article is available on GitHub: https://github.com/dcgerard/reproduce_prior_sims.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Unraveling the Hexaploid Sweetpotato Inheritance Using Ultra-Dense Multilocus Mapping

The hexaploid sweetpotato (Ipomoea batatas (L.) Lam., 2n = 6x = 90) is an important staple food crop worldwide and plays a vital role in alleviating famine in developing countries. Due to its high ploidy level, genetic studies in sweetpotato lag behind major diploid crops significantly. We built an ultra-dense multilocus integrated genetic map and characterized the inheritance system in a sweetpotato full-sib family using our newly developed software, MAPpoly. The resulting genetic map revealed 96.5% collinearity between I. batatas and its diploid relative I. trifida We computed the genotypic probabilities across the whole genome for all individuals in the mapping population and inferred their complete hexaploid haplotypes. We provide evidence that most of the meiotic configurations (73.3%) were resolved in bivalents, although a small portion of multivalent signatures (15.7%), among other inconclusive configurations (11.0%), were also observed. Except for low levels of preferential pairing in linkage group 2, we observed a hexasomic inheritance mechanism in all linkage groups. We propose that the hexasomic-bivalent inheritance promotes stability to the allelic transmission in sweetpotato.

Unraveling the Hexaploid Sweetpotato Inheritance Using Ultra-Dense Multilocus Mapping

The hexaploid sweetpotato (Ipomoea batatas (L.) Lam., 2n = 6x = 90) is an important staple food crop worldwide and plays a vital role in alleviating famine in developing countries. Due to its high ploidy level, genetic studies in sweetpotato lag behind major diploid crops significantly. We built an ultra-dense multilocus integrated genetic map and characterized the inheritance system in a sweetpotato full-sib family using our newly developed software, MAPpoly. The resulting genetic map revealed 96.5% collinearity between I. batatas and its diploid relative I. trifida We computed the genotypic probabilities across the whole genome for all individuals in the mapping population and inferred their complete hexaploid haplotypes. We provide evidence that most of the meiotic configurations (73.3%) were resolved in bivalents, although a small portion of multivalent signatures (15.7%), among other inconclusive configurations (11.0%), were also observed. Except for low levels of preferential pairing in linkage group 2, we observed a hexasomic inheritance mechanism in all linkage groups. We propose that the hexasomic-bivalent inheritance promotes stability to the allelic transmission in sweetpotato.

Genome-wide genetic diversity detection and population structure analysis in sweetpotato (Ipomoea batatas) using RAD-seq

Sweetpotato (Ipomoea batatas L.) is one of the most important food and grain-forage crops globally. It has been planted in >100 countries. Due to the complexity of the sweetpotato genome, its research is far behind other major food crops. At present, limited information about the sweetpotato genome is available. Thus, it is central to find an efficient approach for the investigation of sweetpotato genome. In this study, RAD-seq (Restriction site-associated DNA sequencing) was used to evaluate sweetpotato genetic structure diversity and to develop relevant SSR markers. The study yielded >128 Gb reliable sequence data from 81 sweetpotato accessions. By analyzing polymorphic tags from each accession, a total of 55,622 restriction-site associated DNA sequencing tags (RAD-seq) were found, containing 907,010 SNP. Genetic analysis divided 81 accessions into five major clusters based on their SNP genotype, which matches the results of genetic analysis and the genetic family tree. In addition, 18,320 SSRs loci were detected and 9336 SSR primer pairs were developed. Eighty-three primer pairs were amplified in different sweetpotato genotypes, 76 of which successfully amplified polymorphism bands. These results provide significant information about sweetpotato genome, which can be used to identify novel gene and to further develop the gene chip. And more significant, clustering results based on the SNP genotype provide an essential reference for breeders to match parent plants in breeding program. Additionally, SSR markers developed in this study will supply a wealth of markers for marker-assisted selection in sweetpotato 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.

Current status in whole genome sequencing and analysis of Ipomoea spp

The recent advances of next-generation sequencing have made it possible to construct reference genome sequences in divergent species. However, de novo assembly at the chromosome level remains challenging in polyploid species, due to the existence of more than two pairs of homoeologous chromosomes in one nucleus. Cultivated sweet potato (Ipomoea batatas (L.) Lam) is a hexaploid species with 90 chromosomes (2n = 6X = 90). Although the origin of sweet potato is also still under discussion, diploid relative species, I. trifida and I. triloba have been considered as one of the most possible progenitors. In this manuscript, we review the recent results and activities of whole-genome sequencing in the genus Ipomoea series Batatas, I. trifida, I. triloba and sweet potato (I. batatas). Most of the results of genome assembly suggest that the genomes of sweet potato consist of two pairs and four pairs of subgenomes, i.e., B1B1B2B2B2B2. The results also revealed the relation between sweet potato and other Ipomoea species. Together with the development of bioinformatics approaches, the large-scale publicly available genome and transcript sequence resources and international genome sequencing streams are expected to promote the genome sequence dissection in sweet potato.

Linkage Analysis and Haplotype Phasing in Experimental Autopolyploid Populations with High Ploidy Level Using Hidden Markov Models

Modern SNP genotyping technologies allow measurement of the relative abundance of different alleles for a given locus and consequently estimation of their allele dosage, opening a new road for genetic studies in autopolyploids. Despite advances in genetic linkage analysis in autotetraploids, there is a lack of statistical models to perform linkage analysis in organisms with higher ploidy levels. In this paper, we present a statistical method to estimate recombination fractions and infer linkage phases in full-sib populations of autopolyploid species with even ploidy levels for a set of SNP markers using hidden Markov models. Our method uses efficient two-point procedures to reduce the search space for the best linkage phase configuration and reestimate the final parameters by maximizing the likelihood of the Markov chain. To evaluate the method, and demonstrate its properties, we rely on simulations of autotetraploid, autohexaploid and autooctaploid populations and on a real tetraploid potato data set. The results show the reliability of our approach, including situations with complex linkage phase scenarios in hexaploid and octaploid populations.

Linkage Analysis and Haplotype Phasing in Experimental Autopolyploid Populations with High Ploidy Level Using Hidden Markov Models

Modern SNP genotyping technologies allow measurement of the relative abundance of different alleles for a given locus and consequently estimation of their allele dosage, opening a new road for genetic studies in autopolyploids. Despite advances in genetic linkage analysis in autotetraploids, there is a lack of statistical models to perform linkage analysis in organisms with higher ploidy levels. In this paper, we present a statistical method to estimate recombination fractions and infer linkage phases in full-sib populations of autopolyploid species with even ploidy levels for a set of SNP markers using hidden Markov models. Our method uses efficient two-point procedures to reduce the search space for the best linkage phase configuration and reestimate the final parameters by maximizing the likelihood of the Markov chain. To evaluate the method, and demonstrate its properties, we rely on simulations of autotetraploid, autohexaploid and autooctaploid populations and on a real tetraploid potato data set. The results show the reliability of our approach, including situations with complex linkage phase scenarios in hexaploid and octaploid populations.

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.

Analyses of single nucleotide polymorphisms identified by ddRAD-seq reveal genetic structure of tea germplasm and Japanese landraces for tea breeding

To obtain genetic information about the germplasm of tea (Camellia sinensis L.) in Japan, 167 accessions including 138 var. sinensis (96 Japanese var. sinensis and 42 exotic var. sinensis) and 29 Assam hybrids were analyzed using single nucleotide polymorphisms (SNPs) markers identified by double-digest restriction-site-associated DNA sequencing (ddRAD-seq) analysis. Approximately 10,000 SNPs were identified by ddRAD-seq and were mapped across the whole genome. The 167 tea accessions were classified into three genetic subgroups: (1) Japanese var. sinensis; (2) Japanese and exotic var. sinensis; (3) Assam hybrids and exotic var. sinensis. Leaf morphology varied widely within each genetic subgroups. The 96 Japanese var. sinensis were classified into four genetic subgroups as follows; two subgroups of Shizuoka (the largest tea production region) landraces, Uji (most ancient tea production region) landraces, and the pedigree of 'Yabukita', the leading green tea cultivar in Japan. These results indicated that the SNP markers obtained from ddRAD-seq are a useful tool to investigate the geographical background and breeding history of Japanese tea. This genetic information revealed the ancestral admixture situation of the 'Yabukita' pedigree, and showed that the genome structure of 'Yabukita' is clearly different from those of other Japanese accessions.

De novo whole-genome assembly in Chrysanthemum seticuspe, a model species of Chrysanthemums, and its application to genetic and gene discovery analysis

Cultivated chrysanthemum (Chrysanthemum morifolium Ramat.) is one of the most economically important ornamental crops grown worldwide. It has a complex hexaploid genome (2n = 6x = 54) and large genome size. The diploid Chrysanthemum seticuspe is often used as a model of cultivated chrysanthemum, since the two species are closely related. To expand our knowledge of the cultivated chrysanthemum, we here performed de novo whole-genome assembly in C. seticuspe using the Illumina sequencing platform. XMRS10, a C. seticuspe accession developed by five generations of self-crossing from a self-compatible strain, AEV2, was used for genome sequencing. The 2.72 Gb of assembled sequences (CSE_r1.0), consisting of 354,212 scaffolds, covered 89.0% of the 3.06 Gb C. seticuspe genome estimated by k-mer analysis. The N50 length of scaffolds was 44,741 bp. For protein-encoding genes, 71,057 annotated genes were deduced (CSE_r1.1_cds). Next, based on the assembled genome sequences, we performed linkage map construction, gene discovery and comparative analyses for C. seticuspe and cultivated chrysanthemum. The generated C. seticuspe linkage map revealed skewed regions in segregation on the AEV2 genome. In gene discovery analysis, candidate flowering-related genes were newly found in CSE_r1.1_cds. Moreover, single nucleotide polymorphism identification and annotation on the C. × morifolium genome showed that the C. seticuspe genome was applicable to genetic analysis in cultivated chrysanthemums. The genome sequences assembled herein are expected to contribute to future chrysanthemum studies. In addition, our approach demonstrated the usefulness of short-read genome assembly and the importance of choosing an appropriate next genome sequencing technology based on the purpose of the post-genome analysis.

High density linkage map construction and QTL mapping for runner production in allo-octoploid strawberry Fragaria × ananassa based on ddRAD-seq derived SNPs

Recent advances in high-throughput genome sequencing technologies are now making the genetic dissection of the complex genome of cultivated strawberry easier. We sequenced Maehyang (short-day cultivar) × Albion (day-neutral cultivar) crossing populations using double digest restriction-associated DNA (ddRAD) sequencing technique that yielded 978,968 reads, 80.2% of which were aligned to strawberry genome allowing the identification of 13,181 high quality single nucleotide polymorphisms (SNPs). Total 3051 SNPs showed Mendelian segregation in F₁, of which 1268 were successfully mapped to 46 linkage groups (LG) spanning a total of 2581.57 cM with an average interval genetic distance of 2.22 cM. The LGs were assigned to the 28 chromosomes of Fragaria × ananassa as determined by positioning the sequence tags on F. vesca genome. In addition, seven QTLs namely, qRU-5D, qRU-3D1, qRU-1D2, qRU-4D, qRU-4C, qRU-5C and qRU-2D2 were identified for runner production with LOD value ranging from 3.5–7.24 that explained 22–38% of phenotypic variation. The key candidate genes having putative roles in meristem differentiation for runnering and flowering within these QTL regions were identified. These will enhance our understanding of the vegetative vs sexual reproductive behavior in strawberry and will aid in setting breeding targets for developing perpetual flowering and profuse runnering cultivar.

Genetic Mapping With Allele Dosage Information in Tetraploid Urochloa decumbens (Stapf) R. D. Webster Reveals Insights Into Spittlebug (Notozulia entreriana Berg) Resistance

Urochloa decumbens (Stapf) R. D. Webster is one of the most important African forage grasses in Brazilian beef production. Currently available genetic-genomic resources for this species are restricted mainly due to polyploidy and apomixis. Therefore, crucial genomic-molecular studies such as the construction of genetic maps and the mapping of quantitative trait loci (QTLs) are very challenging and consequently affect the advancement of molecular breeding. The objectives of this work were to (i) construct an integrated U. decumbens genetic map for a full-sibling progeny using GBS-based markers with allele dosage information, (ii) detect QTLs for spittlebug (Notozulia entreriana) resistance, and (iii) seek putative candidate genes involved in defense against biotic stresses. We used the Setaria viridis genome a reference to align GBS reads and selected 4,240 high-quality SNP markers with allele dosage information. Of these markers, 1,000 were distributed throughout nine homologous groups with a cumulative map length of 1,335.09 cM and an average marker density of 1.33 cM. We detected QTLs for resistance to spittlebug, an important pasture insect pest, that explained between 4.66 and 6.24% of the phenotypic variation. These QTLs are in regions containing putative candidate genes related to defense against biotic stresses. Because this is the first genetic map with SNP autotetraploid dosage data and QTL detection in U. decumbens, it will be useful for future evolutionary studies, genome assembly, and other QTL analyses in Urochloa spp. Moreover, the results might facilitate the isolation of spittlebug-related candidate genes and help clarify the mechanism of spittlebug resistance. These approaches will improve selection efficiency and accuracy in U. decumbens molecular breeding and shorten the breeding cycle.

Genome-wide association studies on the phyllosphere microbiome: Embracing complexity in host-microbe interactions

Environmental sequencing shows that plants harbor complex communities of microbes that vary across environments. However, many approaches for mapping plant genetic variation to microbe-related traits were developed in the relatively simple context of binary host-microbe interactions under controlled conditions. Recent advances in sequencing and statistics make genome-wide association studies (GWAS) an increasingly promising approach for identifying the plant genetic variation associated with microbes in a community context. This review discusses early efforts on GWAS of the plant phyllosphere microbiome and the outlook for future studies based on human microbiome GWAS. A workflow for GWAS of the phyllosphere microbiome is then presented, with particular attention to how perspectives on the mechanisms, evolution and environmental dependence of plant-microbe interactions will influence the choice of traits to be mapped.

Genotyping Polyploids from Messy Sequencing Data

Detecting and quantifying the differences in individual genomes (i.e., genotyping), plays a fundamental role in most modern bioinformatics pipelines. Many scientists now use reduced representation next-generation sequencing (NGS) approaches for genotyping. Genotyping diploid individuals using NGS is a well-studied field, and similar methods for polyploid individuals are just emerging. However, there are many aspects of NGS data, particularly in polyploids, that remain unexplored by most methods. Our contributions in this paper are fourfold: (i) We draw attention to, and then model, common aspects of NGS data: sequencing error, allelic bias, overdispersion, and outlying observations. (ii) Many datasets feature related individuals, and so we use the structure of Mendelian segregation to build an empirical Bayes approach for genotyping polyploid individuals. (iii) We develop novel models to account for preferential pairing of chromosomes, and harness these for genotyping. (iv) We derive oracle genotyping error rates that may be used for read depth suggestions. We assess the accuracy of our method in simulations, and apply it to a dataset of hexaploid sweet potato (Ipomoea batatas). An R package implementing our method is available at https://cran.r-project.org/package=updog.

Tools for Genetic Studies in Experimental Populations of Polyploids

Polyploid organisms carry more than two copies of each chromosome, a condition rarely tolerated in animals but which occurs relatively frequently in the plant kingdom. One of the principal challenges faced by polyploid organisms is to evolve stable meiotic mechanisms to faithfully transmit genetic information to the next generation upon which the study of inheritance is based. In this review we look at the tools available to the research community to better understand polyploid inheritance, many of which have only recently been developed. Most of these tools are intended for experimental populations (rather than natural populations), facilitating genomics-assisted crop improvement and plant breeding. This is hardly surprising given that a large proportion of domesticated plant species are polyploid. We focus on three main areas: (1) polyploid genotyping; (2) genetic and physical mapping; and (3) quantitative trait analysis and genomic selection. We also briefly review some miscellaneous topics such as the mode of inheritance and the availability of polyploid simulation software. The current polyploid analytic toolbox includes software for assigning marker genotypes (and in particular, estimating the dosage of marker alleles in the heterozygous condition), establishing chromosome-scale linkage phase among marker alleles, constructing (short-range) haplotypes, generating linkage maps, performing genome-wide association studies (GWAS) and quantitative trait locus (QTL) analyses, and simulating polyploid populations. These tools can also help elucidate the mode of inheritance (disomic, polysomic or a mixture of both as in segmental allopolyploids) or reveal whether double reduction and multivalent chromosomal pairing occur. An increasing number of polyploids (or associated diploids) are being sequenced, leading to publicly available reference genome assemblies. Much work remains in order to keep pace with developments in genomic technologies. However, such technologies also offer the promise of understanding polyploid genomes at a level which hitherto has remained elusive.

Rapid identification of candidate genes for resistance to tomato late blight disease using next-generation sequencing technologies

Tomato late blight caused by Phytophthora infestans (Mont.) de Bary, also known as the Irish famine pathogen, is one of the most destructive plant diseases. Wild relatives of tomato possess useful resistance genes against this disease, and could therefore be used in breeding to improve cultivated varieties. In the genome of a wild relative of tomato, Solanum habrochaites accession LA1777, we identified a new quantitative trait locus for resistance against blight caused by an aggressive Egyptian isolate of P. infestans. Using double-digest restriction site-associated DNA sequencing (ddRAD-Seq) technology, we determined 6,514 genome-wide SNP genotypes of an F2 population derived from an interspecific cross. Subsequent association analysis of genotypes and phenotypes of the mapping population revealed that a 6.8 Mb genome region on chromosome 6 was a candidate locus for disease resistance. Whole-genome resequencing analysis revealed that 298 genes in this region potentially had functional differences between the parental lines. Among of them, two genes with missense mutations, Solyc06g071810.1 and Solyc06g083640.3, were considered to be potential candidates for disease resistance. SNP and SSR markers linking to this region can be used in marker-assisted selection in future breeding programs for late blight disease, including introgression of new genetic loci from wild species. In addition, the approach developed in this study provides a model for identification of other genes for attractive agronomical traits.

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.