Construction of reference chromosome-scale pseudomolecules for potato: integrating the potato genome with genetic and physical maps

Characterisation and mapping of a Globodera pallida resistance derived from the wild potato species Solanum spegazzinii

KEY MESSAGE

A new resistance locus acting against the potato cyst nematode Globodera pallida was mapped to chromosome VI in the diploid wild potato species Solanum spegazzinii CPC 7195. The potato cyst nematodes (PCN) Globodera pallida and Globodera rostochiensis are economically important potato pests in almost all regions where potato is grown. One important management strategy involves deployment through introgression breeding into modern cultivars of new sources of naturally occurring resistance from wild potato species. We describe a new source of resistance to G. pallida from wild potato germplasm. The diploid species Solanum spegazzinii Bitter accession CPC 7195 shows resistance to G. pallida pathotypes Pa1 and Pa2/3. A cross and first backcross of S. spegazzinii with Solanum tuberosum Group Phureja cultivar Mayan Gold were performed, and the level of resistance to G. pallida Pa2/3 was determined in progeny clones. Bulk-segregant analysis (BSA) using generic mapping enrichment sequencing (GenSeq) and genotyping-by-sequencing were performed to identify single-nucleotide polymorphisms (SNPs) that are genetically linked to the resistance, using S. tuberosum Group Phureja clone DM1-3 516 R44 as a reference genome. These SNPs were converted into allele-specific PCR assays, and the resistance was mapped to an interval of roughly 118 kb on chromosome VI. This newly identified resistance, which we call Gpa VIlspg, can be used in future efforts to produce modern cultivars with enhanced and broad-spectrum resistances to the major pests and pathogens of potato.

QTL identified that influence tuber length-width ratio, degree of flatness, tuber size, and specific gravity in a russet-skinned, tetraploid mapping population

Potato tuber shape, size, and specific gravity are important agronomic traits in the russet market class of potatoes with an impact on quality, consistency, and product recovery of processed foods such as French fries. Therefore, identifying genetic regions associated with the three traits through quantitative trait locus/loci (QTL) analysis is a crucial process in the subsequent development of marker-assisted selection for use in potato breeding programs. QTL analysis was conducted on a tetraploid mapping population consisting of 190 individuals derived from the cross between two russet-skinned parents, Palisade Russet and the breeding clone ND028673B-2Russ. Field data collected over a 2-year period and used in the QTL analyses included tuber length-width and width-depth ratios that were obtained using a digital caliper. The width-depth ratio provided an assessment of the "flatness" of a tuber, which is of importance in potato processing. To cross-validate the accuracy and differences among tuber shape measurement methods, a trained evaluator also assessed the identical tubers based on 1-5 scale (compressed to long) visual assessment method. Furthermore, the weights of analyzed tubers and specific gravities were also collected during the phenotyping process for each mapping clone. A major tuber shape QTL was consistently observed on chromosome 10 with both the length-width ratio and visual assessments. On chromosome 4, a significant QTL for tuber shape from the visual assessment phenotypic data was also detected. Additionally, a tuber shape-related QTL on chromosome 6 was also detected from the length-width ratio data from 2020. Chromosome 2 was also identified as having a significant QTL for the width-depth ratio, which is of importance in influencing the flatness of a tuber. One significant QTL for tuber weight (i.e., tuber size) was observed on chromosome 5, and a significant QTL for specific gravity was found on chromosome 3. These significant and major QTL should be useful for developing marker-assisted selection for more efficient potato breeding.

Transcriptome Analysis Reveals Novel Genes Potentially Involved in Tuberization in Potato

The formation and development of tubers, the primary edible and economic organ of potatoes, directly affect their yield and quality. The regulatory network and mechanism of tuberization have been preliminarily revealed in recent years, but plenty of relevant genes remain to be discovered. A few candidate genes were provided due to the simplicity of sampling and result analysis of previous transcriptomes related to tuberization. We sequenced and thoroughly analyzed the transcriptomes of thirteen tissues from potato plants at the tuber proliferation phase to provide more reference information and gene resources. Among them, eight tissues were stolons and tubers at different developmental stages, which we focused on. Five critical periods of tuberization were selected to perform an analysis of differentially expressed genes (DEGs), according to the results of the tissue correlation. Compared with the unswollen stolons (Sto), 2751, 4897, 6635, and 9700 DEGs were detected in the slightly swollen stolons (Sto1), swollen stolons (Sto2), tubers of proliferation stage 1 (Tu1), and tubers of proliferation stage 4 (Tu4). A total of 854 transcription factors and 164 hormone pathway genes were identified in the DEGs. Furthermore, three co-expression networks associated with Sto-Sto1, Sto2-Tu1, and tubers of proliferation stages two to five (Tu2-Tu5) were built using the weighted gene co-expression network analysis (WGCNA). Thirty hub genes (HGs) and 30 hub transcription factors (HTFs) were screened and focalized in these networks. We found that five HGs were reported to regulate tuberization, and most of the remaining HGs and HTFs co-expressed with them. The orthologs of these HGs and HTFs were reported to regulate processes (e.g., flowering, cell division, hormone synthesis, metabolism and signal transduction, sucrose transport, and starch synthesis) that were also required for tuberization. Such results further support their potential to control tuberization. Our study provides insights and countless candidate genes of the regulatory network of tuberization, laying the foundation for further elucidating the genetic basis of tuber development.

The effect of marker types and density on genomic prediction and GWAS of key performance traits in tetraploid potato

Genomic prediction and genome-wide association studies are becoming widely employed in potato key performance trait QTL identifications and to support potato breeding using genomic selection. Elite cultivars are tetraploid and highly heterozygous but also share many common ancestors and generation-spanning inbreeding events, resulting from the clonal propagation of potatoes through seed potatoes. Consequentially, many SNP markers are not in a 1:1 relationship with a single allele variant but shared over several alleles that might exert varying effects on a given trait. The impact of such redundant "diluted" predictors on the statistical models underpinning genome-wide association studies (GWAS) and genomic prediction has scarcely been evaluated despite the potential impact on model accuracy and performance. We evaluated the impact of marker location, marker type, and marker density on the genomic prediction and GWAS of five key performance traits in tetraploid potato (chipping quality, dry matter content, length/width ratio, senescence, and yield). A 762-offspring panel of a diallel cross of 18 elite cultivars was genotyped by sequencing, and markers were annotated according to a reference genome. Genomic prediction models (GBLUP) were trained on four marker subsets [non-synonymous (29,553 SNPs), synonymous (31,229), non-coding (32,388), and a combination], and robustness to marker reduction was investigated. Single-marker regression GWAS was performed for each trait and marker subset. The best cross-validated prediction correlation coefficients of 0.54, 0.75, 0.49, 0.35, and 0.28 were obtained for chipping quality, dry matter content, length/width ratio, senescence, and yield, respectively. The trait prediction abilities were similar across all marker types, with only non-synonymous variants improving yield predictive ability by 16%. Marker reduction response did not depend on marker type but rather on trait. Traits with high predictive abilities, e.g., dry matter content, reached a plateau using fewer markers than traits with intermediate-low correlations, such as yield. The predictions were unbiased across all traits, marker types, and all marker densities >100 SNPs. Our results suggest that using non-synonymous variants does not enhance the performance of genomic prediction of most traits. The major known QTLs were identified by GWAS and were reproducible across exonic and whole-genome variant sets for dry matter content, length/width ratio, and senescence. In contrast, minor QTL detection was marker type dependent.

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.

Revitalizing agriculture: next-generation genotyping and -omics technologies enabling molecular prediction of resilient traits in the Solanaceae family

This review highlights -omics research in Solanaceae family, with a particular focus on resilient traits. Extensive research has enriched our understanding of Solanaceae genomics and genetics, with historical varietal development mainly focusing on disease resistance and cultivar improvement but shifting the emphasis towards unveiling resilience mechanisms in genebank-preserved germplasm is nowadays crucial. Collecting such information, might help researchers and breeders developing new experimental design, providing an overview of the state of the art of the most advanced approaches for the identification of the genetic elements laying behind resilience. Building this starting point, we aim at providing a useful tool for tackling the global agricultural resilience goals in these crops.

Potato and sweetpotato breeding at the international potato center: approaches, outcomes and the way forward

Root and tuber crop breeding is at the front and center of CIP's science program, which seeks to develop and disseminate sustainable agri-food technologies, information and practices to serve objectives including poverty alleviation, income generation, food security and the sustainable use of natural resources. CIP was established in 1971 in Peru, which is part of potato's center of origin and diversity, with an initial mandate on potato and expanding to include sweetpotato in 1986. Potato and sweetpotato are among the top 10 most consumed food staples globally and provide some of the most affordable sources of energy and vital nutrients. Sweetpotato plays a key role in securing food for many households in Africa and South Asia, while potato is important worldwide. Both crops grow in a range of conditions with relatively few inputs and simple agronomic techniques. Potato is adapted to the cooler environments, while sweetpotato grows well in hot climates, and hence, the two crops complement each other. Germplasm enhancement (pre-breeding), the development of new varieties and building capacity for breeding and variety testing in changing climates with emphasis on adaptation, resistance, nutritional quality and resource-use efficiency are CIP's central activities with significant benefits to the poor. Investments in potato and sweetpotato breeding and allied disciplines at CIP have resulted in the release of many varieties some of which have had documented impact in the release countries. Partnership with diverse types of organizations has been key to the centers way of working toward improving livelihoods through crop production in the global South.

Transcriptomic Analyses Reveal the Role of Cytokinin and the Nodal Stem in Microtuber Sprouting in Potato (Solanum tuberosum L.)

In potatoes, tuber secondary growth, especially sprouting, deforms the tubers and severely lowers their commercial value. Tuber sprouting is induced by signal substances, such as gibberellin (GA), which are transported to the tuber from the plant body. The molecular mechanism underlying GA-induced sprouting remains ambiguous. Here, we tried to recreate tuber secondary growth using in vitro stemmed microtubers (MTs) (with the nodal stem attached) and MT halves (with the nodal stem entirely removed). Our experiments showed that GA alone could initiate the sprouting of stemmed microtubers; however, GA failed to initiate MT halves unless 6-benzyladenine, a synthetic cytokinin CK, was co-applied. Here, we analyzed the transcriptional profiles of sprouting buds using these in vitro MTs. RNA-seq analysis revealed a downregulation of cytokinin-activated signaling but an upregulation of the "Zeatin biosynthesis" pathway, as shown by increased expression of CYP735A, CISZOG, and UGT85A1 in sprouting buds; additionally, the upregulation of genes, such as IAA15, IAA22, and SAUR50, associated with auxin-activated signaling and one abscisic acid (ABA) negative regulator, PLY4, plays a vital role during sprouting growth. Our findings indicate that the role of the nodal stem is synonymous with CK in sprouting growth, suggesting that CK signaling and homeostasis are critical to supporting GA-induced sprouting. To effectively control tuber sprouting, more effort is required to be devoted to these critical genes.

Identification of QTL associated with plant vine characteristics and infection response to late blight, early blight, and Verticillium wilt in a tetraploid potato population derived from late blight-resistant Palisade Russet

Potato late blight (causal agent Phytophthora infestans) is a disease of potatoes with economic importance worldwide. Control is primarily through field monitoring and the application of fungicides. Control of late blight with fungicides and host plant resistance is difficult, with documented cases of such control measures failing with the advent of new pathotypes of P. infestans. To better understand host plant resistance and to develop more durable late blight resistance, Quantitative Trait Locus/Loci (QTL) analysis was conducted on a tetraploid mapping population derived from late blight-resistant potato cultivar Palisade Russet. Additionally, QTL analyses for other traits such as Verticillium wilt and early blight resistance, vine size and maturity were performed to identify a potential relationship between multiple traits and prepare genetic resources for molecular markers useful in breeding programs. For this, one hundred ninety progenies from intercrossing Palisade Russet with a late blight susceptible breeding clone (ND028673B-2Russ) were assessed. Two parents and progenies were evaluated over a two-year period for response to infection by the US-8 genotype of P. infestans in inoculated field screenings in Corvallis, Oregon. In Aberdeen, Idaho, the same mapping population was also evaluated for phenotypic response to early blight and Verticillium wilt, and vine size and maturity in a field over a two-year period. After conducting QTL analyses with those collected phenotype data, it was observed that chromosome 5 has a significant QTL for all five traits. Verticillium wilt and vine maturity QTL were also observed on chromosome 1, and vine size QTL was also found on chromosomes 3 and 10. An early blight QTL was also detected on chromosome 2. The QTL identified in this study have the potential for converting into breeder-friendly molecular markers for marker-assisted selection.

Pho1a (plastid starch phosphorylase) is duplicated and essential for normal starch granule phenotype in tubers of Solanum tuberosum L

Reserve starch from seeds and tubers is a crucial plant product for human survival. Much research has been devoted to quantitative and qualitative aspects of starch synthesis and its relation to abiotic factors of importance in agriculture. Certain aspects of genetic factors and enzymes influencing carbon assimilation into starch granules remain elusive after many decades of research. Starch phosphorylase (Pho) can operate, depending on metabolic conditions, in a synthetic and degradative pathway. The plastidial form of the enzyme is one of the most highly expressed genes in potato tubers, and the encoded product is imported into starch-synthesizing amyloplasts. We identified that the genomic locus of a Pho1a-type starch phosphorylase is duplicated in potato. Our study further shows that the enzyme is of importance for a normal starch granule phenotype in tubers. Null mutants created by genome editing display rounded starch granules in an increased number that contained a reduced ratio of apparent amylose in the starch.

A physical map of traits of agronomic importance based on potato and tomato genome sequences

Potato, tomato, pepper, and eggplant are worldwide important crop and vegetable species of the Solanaceae family. Molecular linkage maps of these plants have been constructed and used to map qualitative and quantitative traits of agronomic importance. This research has been undertaken with the vision to identify the molecular basis of agronomic characters on the one hand, and on the other hand, to assist the selection of improved varieties in breeding programs by providing DNA-based markers that are diagnostic for specific agronomic characters. Since 2011, whole genome sequences of tomato and potato became available in public databases. They were used to combine the results of several hundred mapping and map-based cloning studies of phenotypic characters between 1988 and 2022 in physical maps of the twelve tomato and potato chromosomes. The traits evaluated were qualitative and quantitative resistance to pathogenic oomycetes, fungi, bacteria, viruses, nematodes, and insects. Furthermore, quantitative trait loci for yield and sugar content of tomato fruits and potato tubers and maturity or earliness were physically mapped. Cloned genes for pathogen resistance, a few genes underlying quantitative trait loci for yield, sugar content, and maturity, and several hundred candidate genes for these traits were included in the physical maps. The comparison between the physical chromosome maps revealed, in addition to known intrachromosomal inversions, several additional inversions and translocations between the otherwise highly collinear tomato and potato genomes. The integration of the positional information from independent mapping studies revealed the colocalization of qualitative and quantitative loci for resistance to different types of pathogens, called resistance hotspots, suggesting a similar molecular basis. Synteny between potato and tomato with respect to genomic positions of quantitative trait loci was frequently observed, indicating eventual similarity between the underlying genes.

Investigating the genetic components of tuber bruising in a breeding population of tetraploid potatoes

BACKGROUND

Tuber bruising in tetraploid potatoes (Solanum tuberosum) is a trait of economic importance, as it affects tubers' fitness for sale. Understanding the genetic components affecting tuber bruising is a key step in developing potato lines with increased resistance to bruising. As the tetraploid setting renders genetic analyses more complex, there is still much to learn about this complex phenotype. Here, we used capture sequencing data on a panel of half-sibling populations from a breeding programme to perform a genome-wide association analysis (GWAS) for tuber bruising. In addition, we collected transcriptomic data to enrich the GWAS results. However, there is currently no satisfactory method to represent both GWAS and transcriptomics analysis results in a single visualisation and to compare them with existing knowledge about the biological system under study.

RESULTS

When investigating population structure, we found that the STRUCTURE algorithm yielded greater insights than discriminant analysis of principal components (DAPC). Importantly, we found that markers with the highest (though non-significant) association scores were consistent with previous findings on tuber bruising. In addition, new genomic regions were found to be associated with tuber bruising. The GWAS results were backed by the transcriptomics differential expression analysis. The differential expression notably highlighted for the first time the role of two genes involved in cellular strength and mechanical force sensing in tuber resistance to bruising. We proposed a new visualisation, the HIDECAN plot, to integrate the results from the genomics and transcriptomics analyses, along with previous knowledge about genomic regions and candidate genes associated with the trait.

CONCLUSION

This study offers a unique genome-wide exploration of the genetic components of tuber bruising. The role of genetic components affecting cellular strength and resistance to physical force, as well as mechanosensing mechanisms, was highlighted for the first time in the context of tuber bruising. We showcase the usefulness of genomic data from breeding programmes in identifying genomic regions whose association with the trait of interest merit further investigation. We demonstrate how confidence in these discoveries and their biological relevance can be increased by integrating results from transcriptomics analyses. The newly proposed visualisation provides a clear framework to summarise of both genomics and transcriptomics analyses, and places them in the context of previous knowledge on the trait of interest.

Using haplotype and QTL analysis to fix favorable alleles in diploid potato breeding

At present, the potato (Solanum tuberosum L.) of international commerce is autotetraploid, and the complexity of this genetic system creates limitations for breeding. Diploid potato breeding has long been used for population improvement, and because of an improved understanding of the genetics of gametophytic self-incompatibility, there is now sustained interest in the development of uniform F₁ hybrid varieties based on inbred parents. We report here on the use of haplotype and quantitative trait locus (QTL) analysis in a modified backcrossing (BC) scheme, using primary dihaploids of S. tuberosum as the recurrent parental background. In Cycle 1, we selected XD3-36, a self-fertile F₂ individual homozygous for the self-compatibility gene Sli (S-locus inhibitor). Signatures of gametic and zygotic selection were observed at multiple loci in the F₂ generation, including Sli. In the BC₁ cycle, an F₁ population derived from XD3-36 showed a bimodal response for vine maturity, which led to the identification of late versus early alleles in XD3-36 for the gene CDF1 (Cycling DOF Factor 1). Greenhouse phenotypes and haplotype analysis were used to select a vigorous and self-fertile F₂ individual with 43% homozygosity, including for Sli and the early-maturing allele CDF1.3. Partially inbred lines from the BC₁ and BC₂ cycles have been used to initiate new cycles of selection, with the goal of reaching higher homozygosity while maintaining plant vigor, fertility, and yield.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Transcriptomic analysis of wound-healing in Solanum tuberosum (potato) tubers: Evidence for a stepwise induction of suberin-associated genes

Suberin deposition involves both phenolic and aliphatic polymer biosynthesis and deposition in the same tissue. Therefore, any consideration of exploiting suberin for crop enhancement (e.g., enhanced storage, soil borne disease resistance) requires knowledge of both phenolic and aliphatic component biosynthesis and their coordinated, temporal deposition. In the present study, we use a wound-healing potato tuber system to explore global transcriptome changes during the early stages of wound-healing. Wounding leads to initial and substantial transcriptional changes that follow distinctive temporal patterns - primary metabolic pathways were already functional, or up-regulated immediately, and maintained at levels that would allow for precursor carbon skeletons and energy to feed into downstream metabolic processes. Genes involved in pathways for phenolic production (i.e., the shikimate pathway and phenylpropanoid metabolism) were up-regulated early while those involved in aliphatic suberin production (i.e., fatty acid biosynthesis and modification) were transcribed later into the time course. The pattern of accumulation of genes associated with ABA biosynthesis and degradation steps support a role for ABA in regulating aliphatic suberin production. Evaluation of putative Casparian strip membrane-like genes pinpointed wound-responsive candidates that may mediate the suberin deposition process.

Comparative and expression analyses of AP2/ERF genes reveal copy number expansion and potential functions of ERF genes in Solanaceae

BACKGROUND

The AP2/ERF gene family is a superfamily of transcription factors that are important in the response of plants to abiotic stress and development. However, comprehensive research of the AP2/ERF genes in the Solanaceae family is lacking.

RESULTS

Here, we updated the annotation of AP2/ERF genes in the genomes of eight Solanaceae species, as well as Arabidopsis thaliana and Oryza sativa. We identified 2,195 AP2/ERF genes, of which 368 (17%) were newly identified. Based on phylogenetic analyses, we observed expansion of the copy number of these genes, especially those belonging to specific Ethylene-Responsive Factor (ERF) subgroups of the Solanaceae. From the results of chromosomal location and synteny analyses, we identified that the AP2/ERF genes of the pepper (Capsicum annuum), the tomato (Solanum lycopersicum), and the potato (Solanum tuberosum) belonging to ERF subgroups form a tandem array and most of them are species-specific without orthologs in other species, which has led to differentiation of AP2/ERF gene repertory among Solanaceae. We suggest that these genes mainly emerged through recent gene duplication after the divergence of these species. Transcriptome analyses showed that the genes have a putative function in the response of the pepper and tomato to abiotic stress, especially those in ERF subgroups.

CONCLUSIONS

Our findings will provide comprehensive information on AP2/ERF genes and insights into the structural, evolutionary, and functional understanding of the role of these genes in the Solanaceae.

Phenotypic and molecular analyses in diploid and tetraploid genotypes of Solanum tuberosum L. reveal promising genotypes and candidate genes associated with phenolic compounds, ascorbic acid contents, and antioxidant activity

Potato tubers contain biochemical compounds with antioxidant properties that benefit human health. However, the genomic basis of the production of antioxidant compounds in potatoes has largely remained unexplored. Therefore, we report the first genome-wide association study (GWAS) based on 4488 single nucleotide polymorphism (SNP) markers and the phenotypic evaluation of Total Phenols Content (TPC), Ascorbic Acid Content (AAC), and Antioxidant Activity (AA) traits in 404 diverse potato genotypes (84 diploids and 320 tetraploids) conserved at the Colombian germplasm bank that administers AGROSAVIA. The concentration of antioxidant compounds correlated to the skin tuber color and ploidy level. Especially, purple-blackish tetraploid tubers had the highest TPC (2062.41 ± 547.37 mg GAE), while diploid pink-red tubers presented the highest AA (DDPH: 14967.1 ± 4687.79 μmol TE; FRAP: 2208.63 ± 797.35 mg AAE) and AAC (4.52 mg ± 0.68 AA). The index selection allowed us to choose 20 promising genotypes with the highest values for the antioxidant compounds. Genome Association mapping identified 58 SNP-Trait Associations (STAs) with single-locus models and 28 Quantitative Trait Nucleotide (QTNs) with multi-locus models associated with the evaluated traits. Among models, eight STAs/QTNs related to TPC, AAC, and AA were detected in common, flanking seven candidate genes, from which four were pleiotropic. The combination in one population of diploid and tetraploid genotypes enabled the identification of more genetic associations. However, the GWAS analysis implemented independently in populations detected some regions in common between diploids and tetraploids not detected in the mixed population. Candidate genes have molecular functions involved in phenolic compounds, ascorbic acid biosynthesis, and antioxidant responses concerning plant abiotic stress. All candidate genes identified in this study can be used for further expression analysis validation and future implementation in marker-assisted selection pre-breeding platforms targeting fortified materials. Our study further revealed the importance of potato germplasm conserved in national genebanks, such as AGROSAVIA's, as a valuable genetic resource to improve existing potato varieties.

GWAS in tetraploid potato: identification and validation of SNP markers associated with glycoalkaloid content

Genome-wide association studies (GWAS) are a useful tool to unravel the genetic architecture of complex traits, but the results can be difficult to interpret. Population structure, genetic heterogeneity, and rare alleles easily result in false positive or false negative associations. This paper describes the analysis of a GWAS panel combined with three bi-parental mapping populations to validate GWAS results, using phenotypic data for steroidal glycoalkaloid (SGA) accumulation and the ratio (SGR) between the two major glycoalkaloids α-solanine and α-chaconine in potato tubers. SGAs are secondary metabolites in the Solanaceae family, functional as a defence against various pests and pathogens and in high quantities toxic for humans. With GWAS, we identified five quantitative trait loci (QTL) of which Sga1.1, Sgr8.1, and Sga11.1 were validated, but not Sga3.1 and Sgr7.1. In the bi-parental populations, Sga5.1 and Sga7.1 were mapped, but these were not identified with GWAS. The QTLs Sga1.1, Sga7.1, Sgr7.1, and Sgr8.1 co-localize with genes GAME9, GAME 6/GAME 11, SGT1, and SGT2, respectively. For other genes involved in SGA synthesis, no QTLs were identified. The results of this study illustrate a number of pitfalls in GWAS of which population structure seems the most important. We also show that introgression breeding for disease resistance has introduced new haplotypes to the gene pool involved in higher SGA levels in certain pedigrees. Finally, we show that high SGA levels remain unpredictable in potato but that α-solanine/α-chaconine ratio has a predictable outcome with specific SGT1 and SGT2 haplotypes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01344-2.

Identification and characterization of resistance quantitative trait loci against bacterial wilt caused by the Ralstonia solanacearum species complex in potato

Bacterial wilt (BW) caused by the Ralstonia solanacearum species complex (RSSC) represents one of the most serious diseases affecting potato cultivation. The development of BW-resistant cultivars represents the most efficient strategy to control this disease. The resistance-related quantitative trait loci (QTLs) in plants against different RSSC strains have not been studied extensively. Therefore, we performed QTL analysis for evaluating BW resistance using a diploid population derived from Solanum phureja, S. chacoense, and S. tuberosum. Plants cultivated in vitro were inoculated with different strains (phylotype I/biovar 3, phylotype I/biovar 4, and phylotype IV/biovar 2A) and incubated at 24 °C or 28 °C under controlled conditions. Composite interval mapping was performed for the disease indexes using a resistant parent-derived map and a susceptible parent-derived map consisting of single-nucleotide polymorphism markers. We identified five major and five minor resistance QTLs on potato chromosomes 1, 3, 5, 6, 7, 10, and 11. The major QTLs PBWR-3 and PBWR-7 conferred stable resistance against Ralstonia pseudosolanacearum (phylotype I) and Ralstonia syzygii (phylotype IV), whereas PBWR-6b was a strain-specific major resistance QTL against phylotype I/biovar 3 and was more effective at a lower temperature. Therefore, we suggest that broad-spectrum QTLs and strain-specific QTLs can be combined to develop the most effective BW-resistant cultivars for specific areas.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01321-9.

The genetic architectures of vine and skin maturity in tetraploid potato

The genetic architectures of potato vine and skin maturity, as well as the correlation between the traits, were investigated using multiple techniques from quantitative genetics and genomics. Potato vine and skin maturity, which refer to foliar senescence and adherence of the tuber periderm, respectively, are both important to production and therefore breeding. Our objective was to investigate the genetic architectures of these traits in a genome-wide association panel of 586 genotypes, and through joint linkage mapping in a half-diallel subset (N = 397). Skin maturity was measured by image analysis after mechanized harvest 120 days after planting. To correct for the influence of vine maturity on skin maturity under these conditions, the former was used as a covariate in the analysis. The genomic heritability based on a 10 K SNP array was 0.33 for skin maturity vs. 0.46 for vine maturity. Only minor QTLs were detected for skin maturity, the largest being on chromosome 9 and explaining 8% of the variation. As in many previous studies, S. tuberosum Cycling DOF Factor 1 (CDF1) had a large influence on vine maturity, explaining 33% of the variation in the panel as a bi-allelic SNP and 44% in the half-diallel as a multi-allelic QTL. From the estimated effects of the parental haplotypes in the half-diallel and prior knowledge of the allelic series for CDF1, the CDF1 allele for each haplotype was predicted and ultimately confirmed through whole-genome sequencing. The ability to connect statistical alleles from QTL models with biological alleles based on DNA sequencing represents a new milestone in genomics-assisted breeding for tetraploid species.

Rychc Confers Extreme Resistance to Potato virus Y in Potato

The Potato virus Y (PVY) is responsible for huge economic losses for the potato industry worldwide and is the fifth most consequential plant virus globally. The main strategies for virus control are to limit aphid vectors, produce virus-free seed potatoes, and breed virus-resistant varieties. The breeding of PVY-resistant varieties is the safest and most effective method in terms of cost and environmental protection. Rychc, a gene that confers extreme resistance to PVY, is from S. chacoense, which is a wild diploid potato species that is widely used in many PVY-resistant breeding projects. In this study, Rychc was fine mapped and successfully cloned from S. chacoense accession 40-3. We demonstrated that Rychc encodes a TIR-NLR protein by stably transforming a diploid susceptible cultivar named AC142 and a tetraploid potato variety named E3. The Rychc conferred extreme resistance to PVYO, PVYN:O and PVYNTN in both of the genotypes. To investigate the genetic events occurring during the evolution of the Rychc locus, we sequenced 160 Rychc homologs from 13 S. chacoense genotypes. Based on the pattern of sequence identities, 160 Rychc homologs were divided into 11 families. In Family 11 including Rychc, we found evidence for Type I evolutionary patterns with frequent sequence exchanges, obscured orthologous relationships and high non-synonymous to synonymous substitutions (Ka/Ks), which is consistent with rapid diversification and positive selection in response to rapid changes in the PVY genomes. Furthermore, a functional marker named MG64-17 was developed in this study that indicates the phenotype with 100% accuracy and, therefore, can be used for marker-assisted selection in breeding programs that use S. chacoense as a breeding resource.

Potato Cultivar Identification in South Africa Using a Custom SNP Panel

DNA fingerprinting is a molecular technique applied to identify genetic differences between plant cultivars or lines and is used for genetic purity testing. The suitability of single nucleotide polymorphism (SNP) panels for the fingerprinting of tetraploid potato were investigated as a new high throughput, objective, and cost-effective method instead of simple sequence repeats (SSRs) and polyacrylamide gel electrophoresis (PAGE). One-hundred and ninety (190) potato cultivars, including various cultivars currently important in South Africa, were genotyped at 500 SNP positions utilising SeqSNP by LGC Biosearch Technologies. An optimal panel of 25 SNP markers was identified that could discriminate between South African potato cultivars on genetic allele dosage. The genotypes of these SNPs were validated on selected potato genotypes using KASP (Kompetitive Allele Specific PCR) SNP assays. A database of SNP genotype profiles was compiled for all the entries of the germplasm database. The panel of 21 successful SNPs accurately identified the unique potato cultivars in the database. The KASP SNP assays of the successful SNP panel are therefore available for potato DNA fingerprinting as new germplasm, or purity test requests are submitted to ARC-VIMP. This panel provides an objective method for assigning putative cultivar identity to unknown samples submitted for fingerprinting.

Identification of novel PHD-finger genes in pepper by genomic re-annotation and comparative analyses

BACKGROUND

The plant homeodomain (PHD)-finger gene family that belongs to zinc-finger genes, plays an important role in epigenetics by regulating gene expression in eukaryotes. However, inaccurate annotation of PHD-finger genes hinders further downstream comparative, evolutionary, and functional studies.

RESULTS

We performed genome-wide re-annotation in Arabidopsis thaliana (Arabidopsis), Oryza sativa (rice), Capsicum annuum (pepper), Solanum tuberosum (potato), and Solanum lycopersicum (tomato) to better understand the role of PHD-finger genes in these species. Our investigation identified 875 PHD-finger genes, of which 225 (26% of total) were newly identified, including 57 (54%) novel PHD-finger genes in pepper. The PHD-finger genes of the five plant species have various integrated domains that may be responsible for the diversification of structures and functions of these genes. Evolutionary analyses suggest that PHD-finger genes were expanded recently by lineage-specific duplication, especially in pepper and potato, resulting in diverse repertoires of PHD-finger genes among the species. We validated the expression of six newly identified PHD-finger genes in pepper with qRT-PCR. Transcriptome analyses suggest potential functions of PHD-finger genes in response to various abiotic stresses in pepper.

CONCLUSIONS

Our data, including the updated annotation of PHD-finger genes, provide useful information for further evolutionary and functional analyses to better understand the roles of the PHD-finger gene family in pepper.

The chloroplast redox-responsive transcriptome of solanaceous plants reveals significant nuclear gene regulatory motifs associated to stress acclimation

Transcriptomes of solanaceous plants expressing a plastid-targeted antioxidant protein were analysed to identify chloroplast redox networks modulating the expression of nuclear genes associated with stress acclimation. Plastid functions depend on the coordinated expression of nuclear genes, many of them associated to developmental and stress response pathways. Plastid-generated signals mediate this coordination via retrograde signaling, which includes sensing of chloroplast redox state and levels of reactive oxygen species (ROS), although it remains a poorly understood process. Chloroplast redox poise and ROS build-up can be modified by recombinant expression of a plastid-targeted antioxidant protein, i.e., cyanobacterial flavodoxin, with the resulting plants displaying increased tolerance to multiple environmental challenges. Here we analysed the transcriptomes of these flavodoxin-expressing plants to study the coordinated transcriptional responses of the nucleus to the chloroplast redox status and ROS levels during normal growth and stress responses (drought or biotic stress) in tobacco and potato, members of the economically important Solanaceae family. We compared their transcriptomes against those from stressed and mutant plants accumulating ROS in different subcellular compartments and found distinct ROS-related imprints modulated by flavodoxin expression and/or stress. By introducing our datasets in a large-scale interaction network, we identified transcriptional factors related to ROS and stress responses potentially involved in flavodoxin-associated signaling. Finally, we discovered identical cis elements in the promoters of many genes that respond to flavodoxin in the same direction as in wild-type plants under stress, suggesting a priming effect of flavodoxin before stress manifestation. The results provide a genome-wide picture illustrating the relevance of chloroplast redox status on biotic and abiotic stress responses and suggest new cis and trans targets to generate stress-tolerant solanaceous crops.

Spatial control of potato tuberization by the TCP transcription factor BRANCHED1b

The control of carbon allocation, storage and usage is critical for plant growth and development and is exploited for both crop food production and CO₂ capture. Potato tubers are natural carbon reserves in the form of starch that have evolved to allow propagation and survival over winter. They form from stolons, below ground, where they are protected from adverse environmental conditions and animal foraging. We show that BRANCHED1b (BRC1b) acts as a tuberization repressor in aerial axillary buds, which prevents buds from competing in sink strength with stolons. BRC1b loss of function leads to ectopic production of aerial tubers and reduced underground tuberization. In aerial axillary buds, BRC1b promotes dormancy, abscisic acid responses and a reduced number of plasmodesmata. This limits sucrose accumulation and access of the tuberigen protein SP6A. BRC1b also directly interacts with SP6A and blocks its tuber-inducing activity in aerial nodes. Altogether, these actions help promote tuberization underground.

Genome-Wide Analyses of Aspartic Proteases on Potato Genome (Solanum tuberosum): Generating New Tools to Improve the Resistance of Plants to Abiotic Stress

Aspartic proteases are proteolytic enzymes widely distributed in living organisms and viruses. Although they have been extensively studied in many plant species, they are poorly described in potatoes. The present study aimed to identify and characterize S. tuberosum aspartic proteases. Gene structure, chromosome and protein domain organization, phylogeny, and subcellular predicted localization were analyzed and integrated with RNAseq data from different tissues, organs, and conditions focused on abiotic stress. Sixty-two aspartic protease genes were retrieved from the potato genome, distributed in 12 chromosomes. A high number of intronless genes and segmental and tandem duplications were detected. Phylogenetic analysis revealed eight StAP groups, named from StAPI to StAPVIII, that were differentiated into typical (StAPI), nucellin-like (StAPIIIa), and atypical aspartic proteases (StAPII, StAPIIIb to StAPVIII). RNAseq data analyses showed that gene expression was consistent with the presence of cis-acting regulatory elements on StAP promoter regions related to water deficit. The study presents the first identification and characterization of 62 aspartic protease genes and proteins on the potato genome and provides the baseline material for functional gene determinations and potato breeding programs, including gene editing mediated by CRISPR.

Unprecedented organelle genomic variations in morning glories reveal independent evolutionary scenarios of parasitic plants and the diversification of plant mitochondrial complexes

Background

The morning glories (Convolvulaceae) are distributed worldwide and produce economically important crops, medicinal herbs, and ornamentals. Members of this family are diverse in morphological characteristics and trophic modes, including the leafless parasitic Cuscuta (dodders). Organelle genomes were generally used for studying plant phylogeny and genomic variations. Notably, plastomes in parasitic plants always show non-canonical features, such as reduced size and accelerated rates. However, few organelle genomes of this group have been sequenced, hindering our understanding of their evolution, and dodder mitogenome in particular.

Results

We assembled 22 new mitogenomes and 12 new plastomes in Convolvulaceae. Alongside previously known ones, we totally analyzed organelle genomes of 23 species in the family. Our sampling includes 16 leafy autotrophic species and 7 leafless parasitic dodders, covering 8 of the 12 tribes. Both the plastid and mitochondrial genomes of these plants have encountered variations that were rarely observed in other angiosperms. All of the plastomes possessed atypical IR boundaries. Besides the gene and IR losses in dodders, some leafy species also showed gene and intron losses, duplications, structural variations, and insertions of foreign DNAs. The phylogeny reconstructed by plastid protein coding sequences confirmed the previous relationship of the tribes. However, the monophyly of ‘Merremieae’ and the sister group of Cuscuta remained uncertain. The mitogenome was significantly inflated in Cuscuta japonica, which has exceeded over 800 kb and integrated massive DNAs from other species. In other dodders, mitogenomes were maintained in small size, revealing divergent evolutionary strategies. Mutations unique to plants were detected in the mitochondrial gene ccmFc, which has broken into three fragments through gene fission and splicing shift. The unusual changes likely initially happened to the common ancestor of the family and were caused by a foreign insertion from rosids followed by double-strand breaks and imprecise DNA repairs. The coding regions of ccmFc expanded at both sides after the fission, which may have altered the protein structure.

Conclusions

Our family-scale analyses uncovered unusual scenarios for both organelle genomes in Convolvulaceae, especially in parasitic plants. The data provided valuable genetic resources for studying the evolution of Convolvulaceae and plant parasitism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01250-1.

Germplasm, Breeding, and Genomics in Potato Improvement of Biotic and Abiotic Stresses Tolerance

Potato is one of the most important food crops in the world. Late blight, viruses, soil and tuber-borne diseases, insect-pests mainly aphids, whiteflies, and potato tuber moths are the major biotic stresses affecting potato production. Potato is an irrigated and highly fertilizer-responsive crop, and therefore, heat, drought, and nutrient stresses are the key abiotic stresses. The genus Solanum is a reservoir of genetic diversity, however, a little fraction of total diversity has been utilized in potato breeding. The conventional breeding has contributed significantly to the development of potato varieties. In recent years, a tremendous progress has been achieved in the sequencing technologies from short-reads to long-reads sequence data, genomes of Solanum species (i.e., pan-genomics), bioinformatics and multi-omics platforms such as genomics, transcriptomics, proteomics, metabolomics, ionomics, and phenomics. As such, genome editing has been extensively explored as a next-generation breeding tool. With the available high-throughput genotyping facilities and tetraploid allele calling softwares, genomic selection would be a reality in potato in the near future. This mini-review covers an update on germplasm, breeding, and genomics in potato improvement for biotic and abiotic stress tolerance.

Assessment of Genetic Diversity and Relatedness in an Andean Potato Collection from Argentina by High-Density Genotyping

Native potatoes are the most diverse among cultivated potato species and thus constitute a valuable source for identifying genes for potato improvement. Nevertheless, high-density mapping, needed to reveal allelic diversity, has not been performed for native Argentinian potatoes. We present a study of the genetic variability and population structure of 96 Andigena potatoes from Northwestern Argentina performed using a subset of 5035 SNPs with no missing data and full reproducibility. These high-density markers are distributed across the genome and present a good coverage of genomic regions. A Bayesian approach revealed the presence of: (I) a major group comprised of most of the Andean accessions; (II) a smaller group containing the out-group cv. Spunta and the sequenced genotype DM; and (III) a third group containing colored flesh potatoes. This grouping was also consistent when maximum likelihood trees were constructed and further confirmed by a principal coordinate analysis. A group of 19 accessions stored as Andean varieties clustered consistently with group Tuberosum accessions. This was in agreement with previous studies and we hypothesize that they may be reintroductions of European-bred long day-adapted potatoes. The present study constitutes a valuable source for allele mining of genes of interest and thus provides a tool for association mapping studies.

QTL mapping in outbred tetraploid (and diploid) diallel populations

Over the last decade, multiparental populations have become a mainstay of genetics research in diploid species. Our goal was to extend this paradigm to autotetraploids by developing software for quantitative trait locus (QTL) mapping in connected F1 populations derived from a set of shared parents. For QTL discovery, phenotypes are regressed on the dosage of parental haplotypes to estimate additive effects. Statistical properties of the model were explored by simulating half-diallel diploid and tetraploid populations with different population sizes and numbers of parents. Across scenarios, the number of progeny per parental haplotype (pph) largely determined the statistical power for QTL detection and accuracy of the estimated haplotype effects. Multiallelic QTL with heritability 0.2 were detected with 90% probability at 25 pph and genome-wide significance level 0.05, and the additive haplotype effects were estimated with over 90% accuracy. Following QTL discovery, the software enables a comparison of models with multiple QTL and nonadditive effects. To illustrate, we analyzed potato tuber shape in a half-diallel population with three tetraploid parents. A well-known QTL on chromosome 10 was detected, for which the inclusion of digenic dominance lowered the Deviance Information Criterion (DIC) by 17 points compared to the additive model. The final model also contained a minor QTL on chromosome 1, but higher-order dominance and epistatic effects were excluded based on the DIC. In terms of practical impacts, the software is already being used to select offspring based on the effect and dosage of particular haplotypes in breeding programs.

Linkage and QTL mapping for tuber shape and specific gravity in a tetraploid mapping population of potato representing the russet market class

BACKGROUND

Tuber shape and specific gravity (dry matter) are important agronomic traits in potato processing and impact production costs, quality, and consistency of the final processed food products such as French fries and potato chips. In this study, linkage and QTL mapping were performed for these two traits to allow for the implementation of marker-assisted selection to facilitate breeding efforts in the russet market class. Two parents, Rio Grande Russet (female) and Premier Russet (male) and their 205 F1 progenies were initially phenotyped for tuber shape and specific gravity in field trials conducted in Idaho and North Carolina in 2010 and 2011, with specific gravity also being measured in Minnesota in 2011. Progenies and parents were previously genotyped using the Illumina SolCAP Infinium 8303 Potato SNP array, with ClusterCall and MAPpoly (R-packages) subsequently used for autotetraploid SNP calling and linkage mapping in this study. The 12 complete linkage groups and phenotypic data were then imported into QTLpoly, an R-package designed for polyploid QTL analyses.

RESULTS

Significant QTL for tuber shape were detected on chromosomes 4, 7, and 10, with heritability estimates ranging from 0.09 to 0.36. Significant tuber shape QTL on chromosomes 4 and 7 were specific to Idaho and North Carolina environments, respectively, whereas the QTL on chromosome 10 was significant regardless of growing environment. Single marker analyses identified alleles in the parents associated with QTL on chromosomes 4, 7, and 10 that contributed to significant differences in tuber shape among progenies. Significant QTL were also identified for specific gravity on chromosomes 1 and 5 with heritability ranging from 0.12 to 0.21 and were reflected across environments.

CONCLUSION

Fully automated linkage mapping and QTL analysis were conducted to identify significant QTL for tuber shape and dry matter in a tetraploid mapping population representing the russet market class. The findings are important for the development of molecular markers useful to potato breeders for marker-assisted selection for the long tuber shape and acceptable dry matter required by the potato industry within this important market class.

Genome-wide simple sequence repeat markers in potato: abundance, distribution, composition, and polymorphism

Simple sequence repeats (SSRs) are important sources of genetic diversity and are widely used as markers in genetics and molecular breeding. In this study, we examined four potato genomes of DM1-3 516 R44 (DM) from Solanum phureja, RH89039-16 (RH) from Solanum tuberosum, M6 from Solanum chacoense and Solanum commersonii to determine SSR abundance and distribution and develop a larger list of polymorphic markers for a potentially wide range of uses for the potato community. A total of 1,734,619 SSRs were identified across the four genomes with an average of 433,655 SSRs per genome and 2.31kb per SSR. The most abundant repeat units for mono-, di-, tri-, and tetra-nucleotide SSRs were (A/T)n, (AT/AT)n, (AAT/ATT)n, and (ATAT/ATAT)n, respectively. The SSRs were most abundant (78.79%) in intergenic regions and least abundant (3.68%) in untranslated regions. On average, 168,069 SSRs with unique flanking sequences were identified in the four genomes. Further, we identified 16,245 polymorphic SSR markers among the four genomes. Experimental validation confirmed 99.69% of tested markers could generate target bands. The high-density potato SSR markers developed in this study will undoubtedly facilitate the application of SSR markers for genetic research and marker-pyramiding in potato breeding.

Haplotype reconstruction in connected tetraploid F1 populations

In diploid species, many multiparental populations have been developed to increase genetic diversity and quantitative trait loci (QTL) mapping resolution. In these populations, haplotype reconstruction has been used as a standard practice to increase the power of QTL detection in comparison with the marker-based association analysis. However, such software tools for polyploid species are few and limited to a single biparental F1 population. In this study, a statistical framework for haplotype reconstruction has been developed and implemented in the software PolyOrigin for connected tetraploid F1 populations with shared parents, regardless of the number of parents or mating design. Given a genetic or physical map of markers, PolyOrigin first phases parental genotypes, then refines the input marker map, and finally reconstructs offspring haplotypes. PolyOrigin can utilize single nucleotide polymorphism (SNP) data coming from arrays or from sequence-based genotyping; in the latter case, bi-allelic read counts can be used (and are preferred) as input data to minimize the influence of genotype calling errors at low depth. With extensive simulation we show that PolyOrigin is robust to the errors in the input genotypic data and marker map. It works well for various population designs with ≥30 offspring per parent and for sequences with read depth as low as 10x. PolyOrigin was further evaluated using an autotetraploid potato dataset with a 3 × 3 half-diallel mating design. In conclusion, PolyOrigin opens up exciting new possibilities for haplotype analysis in tetraploid breeding populations.

Quantitative Trait Locus Mapping for Common Scab Resistance in a Tetraploid Potato Full-Sib Population

Despite the negative impact of common scab (Streptomyces spp.) on the potato industry, little is known about the genetic architecture of resistance to this bacterial disease in the crop. We evaluated a mapping population (∼150 full sibs) derived from a cross between two tetraploid potatoes ('Atlantic' × B1829-5) in three environments (MN11, PA11, ME12) under natural common scab pressure. Three measures to common scab reaction, namely percentage of scabby tubers and disease area and lesion indices, were found to be highly correlated (>0.76). Because of the large environmental effect, heritability values were zero for all three traits in MN11, but moderate to high in PA11 and ME12 (∼0.44 to 0.79). We identified a single quantitative trait locus (QTL) for lesion index in PA11, ME12, and joint analyses on linkage group 3, explaining ∼22 to 30% of the total variation. The identification of QTL haplotypes and candidate genes contributing to disease resistance can support genomics-assisted breeding approaches in the crop.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Polyploid QTL-seq towards rapid development of tightly linked DNA markers for potato and sweetpotato breeding through whole-genome resequencing

Potato (Solanum tuberosum L.) and sweetpotato (Ipomoea batatas L.), which are nutritionally and commercially important tuberous crops, possess a perplexing heredity because of their autopolyploid genomes. To reduce cross-breeding efforts for selecting superior cultivars from progenies with innumerable combinations of traits, DNA markers tightly linked to agronomical traits are required. To develop DNA markers, we developed a method for quantitative trait loci (QTL) mapping using whole-genome next-generation sequencing (NGS) in autopolyploid crops. To apply the NGS-based bulked segregant method, QTL-seq was modified. (1) Single parent-specific simplex (unique for one homologous chromosome) single-nucleotide polymorphisms (SNPs), which present a simple segregation ratio in the progenies, were exploited by filtering SNPs by SNP index (allele frequency). (2) Clusters of SNPs, which were inherited unevenly between bulked progenies with opposite phenotypes, especially those with an SNP index of 0 for the bulk that did not display the phenotypes of interest, were explored. These modifications allowed for separate tracking of alleles located on each of the multiple homologous chromosomes. By applying this method, clusters of SNPs linked to the potato cyst nematode resistance H1 gene and storage root anthocyanin (AN) content were identified in tetraploid potato and hexaploid sweetpotato, respectively, and completely linked DNA markers were developed at the site of the presented SNPs. Thus, polyploid QTL-seq is a versatile method that is free from specialized manipulation for sequencing and construction of elaborate linkage maps and facilitates rapid development of tightly linked DNA markers in autopolyploid crops, such as potato and sweetpotato.

Chromosome-scale reference genome assembly of a diploid potato clone derived from an elite variety

Potato (Solanum tuberosum L.) is one of the most important crops with a worldwide production of 370 million metric tons. The objectives of this study were (1) to create a high-quality consensus sequence across the two haplotypes of a diploid clone derived from a tetraploid elite variety and assess the sequence divergence from the available potato genome assemblies, as well as among the two haplotypes; (2) to evaluate the new assembly’s usefulness for various genomic methods; and (3) to assess the performance of phasing in diploid and tetraploid clones, using linked-read sequencing technology. We used PacBio long reads coupled with 10x Genomics reads and proximity ligation scaffolding to create the dAg1_v1.0 reference genome sequence. With a final assembly size of 812 Mb, where 750 Mb are anchored to 12 chromosomes, our assembly is larger than other available potato reference sequences and high proportions of properly paired reads were observed for clones unrelated by pedigree to dAg1. Comparisons of the new dAg1_v1.0 sequence to other potato genome sequences point out the high divergence between the different potato varieties and illustrate the potential of using dAg1_v1.0 sequence in breeding applications.

A Comprehensive Phylogenetic Analysis of the MAP4K Family in the Green Lineage

The kinase-mediated phosphorylation impacts every basic cellular process. While mitogen-activated protein kinase technology kinase kinases (MAP4Ks) are evolutionarily conserved, there is no comprehensive overview of the MAP4K family in the green lineage (Viridiplantae). In this study, we identified putative MAP4K members from representative species of the two core groups in the green lineage: Chlorophyta, which is a diverse group of green algae, and Streptophyta, which is mostly freshwater green algae and land plants. From that, we inferred the evolutionary relationships of MAP4K proteins through a phylogenetic reconstruction. Furthermore, we provided a classification of the MAP4Ks in the green lineage into three distinct.

Global multi-environment resistance QTL for foliar late blight resistance in tetraploid potato with tropical adaptation

The identification of environmentally stable and globally predictable resistance to potato late blight is challenged by the clonal and polyploid nature of the crop and the rapid evolution of the pathogen. A diversity panel of tetraploid potato germplasm bred for multiple resistance and quality traits was genotyped by genotyping by sequencing (GBS) and evaluated for late blight resistance in three countries where the International Potato Center (CIP) has established breeding work. Health-indexed, in vitro plants of 380 clones and varieties were distributed from CIP headquarters and tuber seed was produced centrally in Peru, China, and Ethiopia. Phenotypes were recorded following field exposure to local isolates of Phytophthora infestans. QTL explaining resistance in four experiments conducted across the three countries were identified in chromosome IX, and environment-specific QTL were found in chromosomes III, V, and X. Different genetic models were evaluated for prediction ability to identify best performing germplasm in each and all environments. The best prediction ability (0.868) was identified with the genomic best linear unbiased predictors (GBLUPs) when using the diploid marker data and QTL-linked markers as fixed effects. Genotypes with high levels of resistance in all environments were identified from the B3, LBHT, and B3-LTVR populations. The results show that many of the advanced clones bred in Peru for high levels of late blight resistance maintain their resistance in Ethiopia and China, suggesting that the centralized selection strategy has been largely successful.

Neofunctionalisation of the Sli gene leads to self-compatibility and facilitates precision breeding in potato

Genetic gain in potato is hampered by the heterozygous tetraploid genome of cultivated potato. Converting potato into a diploid inbred-line based F1-hybrid crop provides a promising route towards increased genetic gain. The introduction of a dominant S-locus inhibitor (Sli) gene into diploid potato germplasm allows efficient generation of self-fertilized seeds and thus the development of potato inbred lines. Little is known about the structure and function of the Sli locus. Here we describe the mapping of Sli to a 12.6 kb interval on chromosome 12 using a recombinant screen approach. One of two candidate genes present in this interval shows a unique sequence that is exclusively present in self-compatible lines. We describe an expression vector that converts self-incompatible genotypes into self-compatible and a CRISPR-Cas9 vector that converts SC genotypes into SI. The Sli gene encodes an F-box protein that is specifically expressed in pollen from self-compatible plants. A 533 bp insertion in the promotor of that gene leads to a gain of function mutation, which overcomes self-pollen rejection.

The S-locus inhibitor (Sli) gene could allow potato breeding by facilitating production of diploid inbred lines. Here the authors show that Sli encodes an F-box protein with a promoter insertion enhancing expression in pollen can overcome pollen rejection in the styles of diploid potato.

Allelic variants of the NLR protein Rpi-chc1 differentially recognize members of the Phytophthora infestans PexRD12/31 effector superfamily through the leucine-rich repeat domain

Phytophthora infestans is a pathogenic oomycete that causes the infamous potato late blight disease. Resistance (R) genes from diverse Solanum species encode intracellular receptors that trigger effective defense responses upon the recognition of cognate RXLR avirulence (Avr) effector proteins. To deploy these R genes in a durable fashion in agriculture, we need to understand the mechanism of effector recognition and the way the pathogen evades recognition. In this study, we cloned 16 allelic variants of the Rpi-chc1 gene from Solanum chacoense and other Solanum species, and identified the cognate P. infestans RXLR effectors. These tools were used to study effector recognition and co-evolution. Functional and non-functional alleles of Rpi-chc1 encode coiled-coil nucleotide-binding leucine-rich repeat (CNL) proteins, being the first described representatives of the CNL16 family. These alleles have distinct patterns of RXLR effector recognition. While Rpi-chc1.1 recognized multiple PexRD12 (Avrchc1.1) proteins, Rpi-chc1.2 recognized multiple PexRD31 (Avrchc1.2) proteins, both belonging to the PexRD12/31 effector superfamily. Domain swaps between Rpi-chc1.1 and Rpi-chc1.2 revealed that overlapping subdomains in the leucine-rich repeat (LRR) domain are responsible for the difference in effector recognition. This study showed that Rpi-chc1.1 and Rpi-chc1.2 evolved to recognize distinct members of the same PexRD12/31 effector family via the LRR domain. The biased distribution of polymorphisms suggests that exchange of LRRs during host-pathogen co-evolution can lead to novel recognition specificities. These insights will guide future strategies to breed durable resistant varieties.

Genome design of hybrid potato

Reinventing potato from a clonally propagated tetraploid into a seed-propagated diploid, hybrid potato, is an important innovation in agriculture. Due to deleterious mutations, it has remained a challenge to develop highly homozygous inbred lines, a prerequisite to breed hybrid potato. Here, we employed genome design to develop a generation of pure and fertile potato lines and thereby the uniform, vigorous F₁s. The metrics we applied in genome design included the percentage of genome homozygosity and the number of deleterious mutations in the starting material, the number of segregation distortions in the S₁ population, the haplotype information to infer the break of tight linkage between beneficial and deleterious alleles, and the genome complementarity of the parental lines. This study transforms potato breeding from a slow, non-accumulative mode into a fast-iterative one, thereby potentiating a broad spectrum of benefits to farmers and consumers.

Comparative Analysis of Re-Annotated Genes Provides Insight into Evolutionary Divergence and Expressions of Aquaporin Family in Pepper

Aquaporins (AQPs) are known to have a vital role in water transport in all living organisms including agriculturally important crops, but a comprehensive genomic study of AQPs in pepper has not been implemented. Here, we updated previous gene annotations and generated a total of 259 AQP genes from five plants, including pepper. Phylogenetic and motif analyses revealed that a large proportion of pepper AQP genes belong to the specific subgroup of tonoplast intrinsic protein (TIP) subfamily, TIP4. Chromosomal localization and estimated duplication times illustrated that genes in TIP4 formed a tandem array on the short arm of chromosome 1, resulting from pepper-specific expansion after its divergence with Solanaceae species. Transcriptome analyses under various abiotic stress conditions revealed that transport-, photosystem-, and thylakoid-related genes were generally enriched in expression clusters containing AQP genes in pepper. These results provide valuable genomic resources and insight into the evolutionary mechanism that generate genomic diversity of the AQP gene family in pepper.

The recombination landscape and multiple QTL mapping in a Solanum tuberosum cv. 'Atlantic'-derived F1 population

There are many challenges involved with the genetic analyses of autopolyploid species, such as the tetraploid potato, Solanum tuberosum (2n = 4x = 48). The development of new analytical methods has made it valuable to re-analyze an F₁ population (n = 156) derived from a cross involving 'Atlantic', a widely grown chipping variety in the USA. A fully integrated genetic map with 4285 single nucleotide polymorphisms, spanning 1630 cM, was constructed with MAPpoly software. We observed that bivalent configurations were the most abundant ones (51.0~72.4% depending on parent and linkage group), though multivalent configurations were also observed (2.2~39.2%). Seven traits were evaluated over four years (2006-8 and 2014) and quantitative trait loci (QTL) mapping was carried out using QTLpoly software. Based on a multiple-QTL model approach, we detected 21 QTL for 15 out of 27 trait-year combination phenotypes. A hotspot on linkage group 5 was identified with co-located QTL for maturity, plant yield, specific gravity, and internal heat necrosis resistance evaluated over different years. Additional QTL for specific gravity and dry matter were detected with maturity-corrected phenotypes. Among the genes around QTL peaks, we found those on chromosome 5 that have been previously implicated in maturity (StCDF1) and tuber formation (POTH1). These analyses have the potential to provide insights into the biology and breeding of tetraploid potato and other autopolyploid species.

Genetic diversity and population structure of advanced clones selected over forty years by a potato breeding program in the USA

Knowledge regarding genetic diversity and population structure of breeding materials is essential for crop improvement. The Texas A&M University Potato Breeding Program has a collection of advanced clones selected and maintained in-vitro over a 40-year period. Little is known about its genetic makeup and usefulness for the current breeding program. In this study, 214 potato clones were genotyped with the Infinium Illumina 22 K V3 Potato Array. After filtering, a total of 10,106 single nucleotide polymorphic (SNP) markers were used for analysis. Heterozygosity varied by SNP, with an overall average of 0.59. Three groups of tetraploid clones primarily based on potato market classes, were detected using STRUCTURE software and confirmed by discriminant analysis of principal components. The highest coefficient of differentiation observed between the groups was 0.14. Signatures of selection were uncovered in genes controlling potato flesh and skin color, length of plant cycle and tuberization, and carbohydrate metabolism. A core set of 43 clones was obtained using Core Hunter 3 to develop a sub-collection that retains similar genetic diversity as the whole population, minimize redundancies, and facilitates long-term conservation of genetic resources. The comprehensive molecular characterization of our breeding clone bank collection contributes to understanding the genetic diversity of existing potato resources. This analysis could be applied to other breeding programs and assist in the selection of parents, fingerprinting, protection, and management of the breeding collections.

Unravelling Differences in Candidate Genes for Drought Tolerance in Potato (Solanum tuberosum L.) by Use of New Functional Microsatellite Markers

Potato is regarded as drought sensitive and most vulnerable to climate changes. Its cultivation in drought prone regions or under conditions of more frequent drought periods, especially in subtropical areas, requires intensive research to improve drought tolerance in order to guarantee high yields under limited water supplies. A candidate gene approach was used to develop functional simple sequence repeat (SSR) markers for association studies in potato with the aim to enhance breeding for drought tolerance. SSR primer combinations, mostly surrounding interrupted complex and compound repeats, were derived from 103 candidate genes for drought tolerance. Validation of the SSRs was performed in an association panel representing 34 mainly starch potato cultivars. Seventy-five out of 154 SSR primer combinations (49%) resulted in polymorphic, highly reproducible banding patterns with polymorphic information content (PIC) values between 0.11 and 0.90. Five SSR markers identified allelic differences between the potato cultivars that showed significant associations with drought sensitivity. In all cases, the group of drought-sensitive cultivars showed predominantly an additional allele, indicating that selection against these alleles by marker-assisted breeding might confer drought tolerance. Further studies of these differences in the candidate genes will elucidate their role for an improved performance of potatoes under water-limited conditions.

Transcriptional Changes in Potato Sprouts upon Interaction with Rhizoctonia solani Indicate Pathogen-Induced Interference in the Defence Pathways of Potato

Rhizoctonia solani is the causer of black scurf disease on potatoes and is responsible for high economical losses in global agriculture. In order to increase the limited knowledge of the plants' molecular response to this pathogen, we inoculated potatoes with R. solani AG3-PT isolate Ben3 and carried out RNA sequencing with total RNA extracted from potato sprouts at three and eight days post inoculation (dpi). In this dual RNA-sequencing experiment, the necrotrophic lifestyle of R. solani AG3-PT during early phases of interaction with its host has already been characterised. Here the potato plants' comprehensive transcriptional response to inoculation with R. solani AG3 was evaluated for the first time based on significantly different expressed plant genes extracted with DESeq analysis. Overall, 1640 genes were differentially expressed, comparing control (-Rs) and with R. solani AG3-PT isolate Ben3 inoculated plants (+Rs). Genes involved in the production of anti-fungal proteins and secondary metabolites with antifungal properties were significantly up regulated upon inoculation with R. solani. Gene ontology (GO) terms involved in the regulation of hormone levels (i.e., ethylene (ET) and jasmonic acid (JA) at 3 dpi and salicylic acid (SA) and JA response pathways at 8 dpi) were significantly enriched. Contrastingly, the GO term "response to abiotic stimulus" was down regulated at both time points analysed. These results may support future breeding efforts toward the development of cultivars with higher resistance level to black scurf disease or the development of new control strategies.

Resisting Potato Cyst Nematodes With Resistance

Potato cyst nematodes (PCN) are economically important pests with a worldwide distribution in all temperate regions where potatoes are grown. Because above ground symptoms are non-specific, and detection of cysts in the soil is determined by the intensity of sampling, infestations are frequently spread before they are recognised. PCN cysts are resilient and persistent; their cargo of eggs can remain viable for over two decades, and thus once introduced PCN are very difficult to eradicate. Various control methods have been proposed, with resistant varieties being a key environmentally friendly and effective component of an integrated management programme. Wild and landrace relatives of cultivated potato have provided a source of PCN resistance genes that have been used in breeding programmes with varying levels of success. Producing a PCN resistant variety requires concerted effort over many years before it reaches what can be the biggest hurdle-commercial acceptance. Recent advances in potato genomics have provided tools to rapidly map resistance genes and to develop molecular markers to aid selection during breeding. This review will focus on the translation of these opportunities into durably PCN resistant varieties.

Identification of Solanum Immune Receptors by Bulked Segregant RNA-Seq and High-Throughput Recombinant Screening

The identification, understanding, and deployment of immune receptors are crucial to achieve high-level and durable resistance for crops against pathogens. In potato, many R genes have been identified using map-based cloning strategies. However, this is a challenging and laborious task that involves the development of a high number of molecular markers for the initial mapping, and the screening of thousands of plants for fine mapping. Bulked segregant RNA-Seq (BSR-Seq) has proven to be an efficient technique for the mapping of resistance genes. The RNA from two bulks of plants with contrasting phenotypes is sequenced and analyzed to identify single-nucleotide polymorphism (SNPs) markers linked to the target gene. Subsequently, the SNP markers that are identified can be used to delimit the mapping interval. Additionally, we designed an in vitro recombinant screening strategy that is advantageous for analyzing a large number of plants, in terms of time, space, and cost. Tips and detailed protocols, including BSR-Seq, bioinformatic analysis, and recombinant screening, are provided in this chapter.

A Bioinformatics Pipeline to Identify a Subset of SNPs for Genomics-Assisted Potato Breeding

Modern potato breeding methods following a genomic-led approach provide means for shortening breeding cycles and increasing breeding efficiency across selection cycles. Acquiring genetic data for large breeding populations remains expensive. We present a pipeline to reduce the number of single nucleotide polymorphisms (SNPs) to lower the cost of genotyping. First, we reduced the number of individuals to be genotyped with a high-throughput method according to the multi-trait variation as defined by principal component analysis of phenotypic characteristics. Next, we reduced the number of SNPs by pruning for linkage disequilibrium. By adjusting the square of the correlation coefficient between two adjacent loci, we obtained reduced subsets of SNPs. We subsequently tested these SNP subsets by two methods; (1) a genome-wide association study (GWAS) for marker identification, and (2) genomic selection (GS) to predict genomic estimated breeding values. The results indicate that both GWAS and GS can be done without loss of information after SNP reduction. The pipeline allows for creating custom SNP subsets to cover all variation found in any particular breeding population. Low-throughput genotyping will reduce the genotyping cost associated with large populations, thereby making genomic breeding methods applicable to large potato breeding populations by reducing genotyping costs.