Construction of a chromosome-scale long-read reference genome assembly for potato

Transcriptome, hormonal, and secondary metabolite changes in leaves of DEFENSE NO DEATH 1 (DND1) silenced potato plants

DEFENSE NO DEATH 1 (DND1) is a cyclic nucleotide-gated ion channel protein. Earlier, it was shown that the silencing of DND1 in the potato (Solanum tuberosum L.) leads to resistance to late blight, powdery mildew, and gray mold diseases. At the same time, however, it can reduce plant growth and cause leaf necrosis. To obtain knowledge of the molecular events behind the pleiotropic effect of DND1 downregulation in the potato, metabolite and transcriptome analyses were performed on three DND1 silenced lines of the cultivar 'Désirée.' A massive increase in the salicylic acid content of leaves was detected. Concentrations of jasmonic acid and chlorogenic acid and their derivatives were also elevated. Expression of 1866 genes was altered in the same way in all three DND1 silenced lines, including those related to the synthesis of secondary metabolites. The activation of several alleles of leaf rust, late blight, and other disease resistance genes, as well as the induction of pathogenesis-related genes, was detected. WRKY and NAC transcription factor families were upregulated, whereas bHLHs were downregulated, indicating their central role in transcriptome changes. These results suggest that the maintenance of the constitutive defense state leads to the reduced growth of DND1 silenced potato plants.

Chromosome-scale genome assembly and annotation of the tetraploid potato cultivar Diacol Capiro adapted to the Andean region

Potato (Solanum tuberosum) is an essential crop for food security and is ranked as the third most important crop worldwide for human consumption. The Diacol Capiro cultivar holds the dominant position in Colombian cultivation, primarily catering to the food processing industry. This highly heterozygous, autotetraploid cultivar belongs to the Andigenum group and it stands out for its adaptation to a wide variety of environments spanning altitudes from 1,800 to 3,200 meters above sea level. Here, a chromosome-scale assembly, referred to as DC, is presented for this cultivar. The assembly was generated by combining circular consensus sequencing with proximity ligation Hi-C for the scaffolding and represents 2.369 Gb with 48 pseudochromosomes covering 2,091 Gb and an anchor rate of 88.26%. The reference genome metrics, including an N50 of 50.5 Mb, a BUSCO (Benchmarking Universal Single-Copy Orthologue) score of 99.38%, and an Long Terminal Repeat Assembly Index score of 13.53, collectively signal the achieved high assembly quality. A comprehensive annotation yielded a total of 154,114 genes, and the associated BUSCO score of 95.78% for the annotated sequences attests to their completeness. The number of predicted NLR (Nucleotide-Binding and Leucine-Rich-Repeat genes) was 2107 with a large representation of NBARC (for nucleotide binding domain shared by Apaf-1, certain R gene products, and CED-4) containing domains (99.85%). Further comparative analysis of the proposed annotation-based assembly with high-quality known potato genomes, showed a similar genome metrics with differences in total gene numbers related to the ploidy status. The genome assembly and annotation of DC presented in this study represent a valuable asset for comprehending potato genetics. This resource aids in targeted breeding initiatives and contributes to the creation of enhanced, resilient, and more productive potato varieties, particularly beneficial for countries in Latin America.

Alternative splicing regulation in plants by SP7-like effectors from symbiotic arbuscular mycorrhizal fungi

Most plants in natural ecosystems associate with arbuscular mycorrhizal (AM) fungi to survive soil nutrient limitations. To engage in symbiosis, AM fungi secrete effector molecules that, similar to pathogenic effectors, reprogram plant cells. Here we show that the Glomeromycotina-specific SP7 effector family impacts on the alternative splicing program of their hosts. SP7-like effectors localize at nuclear condensates and interact with the plant mRNA processing machinery, most prominently with the splicing factor SR45 and the core splicing proteins U1-70K and U2AF35. Ectopic expression of these effectors in the crop plant potato and in Arabidopsis induced developmental changes that paralleled to the alternative splicing modulation of a specific subset of genes. We propose that SP7-like proteins act as negative regulators of SR45 to modulate the fate of specific mRNAs in arbuscule-containing cells. Unraveling the communication mechanisms between symbiotic fungi and their host plants will help to identify targets to improve plant nutrition.

Tuber transcriptome analysis reveals a novel WRKY transcription factor StWRKY70 potentially involved in potato pigmentation

Tuber flesh pigmentation, conferred by the presence of secondary metabolite anthocyanins, is one of many key agronomic traits for potato tubers. Although several genes of potato anthocyanin biosynthesis have been reported, transcription factors (TFs) contributing to tuber flesh pigmentation are still not fully understood. In this study, transcriptomic profiling of diploid potato accessions with or without tuber flesh pigmentation was conducted and genes of the anthocyanin biosynthesis pathway were found significantly enriched within the 1435 differentially expressed genes (DEGs). Weighted Gene Co-expression Network Analysis (WGCNA) and connectivity analysis pinpointed a subset of 173 genes closely related to the key biosynthetic gene StDFR. Of the eight transcription factors in the subset, group III WRKY StWRKY70, was chosen for showing high connectivity to StDFR and ten other anthocyanin biosynthetic genes and homology to known WRKYs of anthocyanin pathway. The transient activation assay showed StWRKY70 predominantly stimulated the expression of StDFR and StANS as well as the accumulation of anthocyanins by enhancing the function of the MYB transcription factor StAN1. Furthermore, the interaction between StWRKY70 and StAN1 was verified by Y2H and BiFC. Our analysis discovered a new transcriptional activator StWRKY70 which potentially involved in tuber flesh pigmentation, thus may lay the foundation for deciphering how the WRKY-MYB-bHLH-WD40 (WRKY-MBW) complex regulate the accumulation of anthocyanins and provide new strategies to breed for more nutritious potato varieties with enhanced tuber flesh anthocyanins.

Genome-wide identification and characterization of SNW/SKIP domain-containing proteins in plants

Sessile organisms, such as plants, developed various ways to sense and respond to external and internal stimuli to maximize their fitness through evolutionary time. Transcripts and protein regulation are, among many, the main mechanisms that plants use to respond to environmental changes. SKIP protein is one such, presenting an SNKW interacting domain, which is highly conserved among eukaryotes, where SKI interacting protein acts in regulating key processes. In the present work, many bioinformatics tools, such as phylogenetic relationships, gene structure, physical-chemical properties, conserved motifs, prediction of regulatory cis-elements, chromosomal localization, and protein-protein interaction network, were used to better understand the genome-wide SNW/SKIP domain-containing proteins. In total, 28 proteins containing the SNW/SKIP domain were identified in different plant species, including plants of agronomic interest. Two main protein clusters were formed in phylogenetic analysis, and gene structure analysis revealed that, in general, the coding region had no introns. Also, expression of these genes is possibly induced by abiotic stress stimuli. Primary structure analysis of the proteins revealed the existence of an evolutionarily conserved functional unit. But physicochemical properties show that proteins containing the SNW/SKIP domain are commonly unstable under in vivo conditions. In addition, the protein network, demonstrated that SKIP homologues could act by modulating plant fitness through gene expression regulation at the transcriptional and post-transcriptional levels. This could be corroborated by the expression number of gene copies of SKIP proteins in many species, highlighting it's crucial role in plant development and tolerance through the course of evolution.

Genotyping-by-sequencing targets genic regions and improves resolution of genome-wide association studies in autotetraploid potato

KEY MESSAGE

De novo genotyping in potato using methylation-sensitive GBS discovers SNPs largely confined to genic or gene-associated regions and displays enhanced effectiveness in estimating LD decay rates, population structure and detecting GWAS associations over 'fixed' SNP genotyping platform. Study also reports the genetic architectures including robust sequence-tagged marker-trait associations for sixteen important potato traits potentially carrying higher transferability across a wider range of germplasm. This study deploys recent advancements in polyploid analytical approaches to perform complex trait analyses in cultivated tetraploid potato. The study employs a 'fixed' SNP Infinium array platform and a 'flexible and open' genome complexity reduction-based sequencing method (GBS, genotyping-by-sequencing) to perform genome-wide association studies (GWAS) for several key potato traits including the assessment of population structure and linkage disequilibrium (LD) in the studied population. GBS SNPs discovered here were largely confined (~ 90%) to genic or gene-associated regions of the genome demonstrating the utility of using a methylation-sensitive restriction enzyme (PstI) for library construction. As compared to Infinium array SNPs, GBS SNPs displayed enhanced effectiveness in estimating LD decay rates and discriminating population subgroups. GWAS using a combined set of 30,363 SNPs identified 189 unique QTL marker-trait associations (QTL-MTAs) covering all studied traits. The majority of the QTL-MTAs were from GBS SNPs potentially illustrating the effectiveness of marker-dense de novo genotyping platforms in overcoming ascertainment bias and providing a more accurate correction for different levels of relatedness in GWAS models. GWAS also detected QTL 'hotspots' for several traits at previously known as well as newly identified genomic locations. Due to the current study exploiting genome-wide genotyping and de novo SNP discovery simultaneously on a large tetraploid panel representing a greater diversity of the cultivated potato gene pool, the reported sequence-tagged MTAs are likely to have higher transferability across a wider range of potato germplasm and increased utility for expediting genomics-assisted breeding for the several complex traits studied.

Modern Breeding Strategies and Tools for Durable Late Blight Resistance in Potato

Cultivated potato (Solanum tuberosum) is a major crop worldwide. It occupies the second place after cereals (corn, rice, and wheat). This important crop is threatened by the Oomycete Phytophthora infestans, the agent of late blight disease. This pathogen was first encountered during the Irish famine during the 1840s and is a reemerging threat to potatoes. It is mainly controlled chemically by using fungicides, but due to health and environmental concerns, the best alternative is resistance. When there is no disease, no treatment is required. In this study, we present a summary of the ongoing efforts concerning resistance breeding of potato against this devastating pathogen, P. infestans. This work begins with the search for and selection of resistance genes, whether they are from within or from outside the species. The genetic methods developed to date for gene mining, such as effectoromics and GWAS, provide researchers with the ability to identify genes of interest more efficiently. Once identified, these genes are cloned using molecular markers (MAS or QRL) and can then be introduced into different cultivars using somatic hybridization or recombinant DNA technology. More innovative technologies have been developed lately, such as gene editing using the CRISPR system or gene silencing, by exploiting iRNA strategies that have emerged as promising tools for managing Phytophthora infestans, which can be employed. Also, gene pyramiding or gene stacking, which involves the accumulation of two or more R genes on the same individual plant, is an innovative method that has yielded many promising results. All these advances related to the development of molecular techniques for obtaining new potato cultivars resistant to P. infestans can contribute not only to reducing losses in agriculture but especially to ensuring food security and safety.

Identification, Elucidation and Deployment of a Cytoplasmic Male Sterility System for Hybrid Potato

Recent advances in diploid F1 hybrid potato breeding rely on the production of inbred lines using the S-locus inhibitor (Sli) gene. As a result of this method, female parent lines are self-fertile and require emasculation before hybrid seed production. The resulting F1 hybrids are self-fertile as well and produce many undesirable berries in the field. Utilization of cytoplasmic male sterility would eliminate the need for emasculation, resulting in more efficient hybrid seed production and male sterile F1 hybrids. We observed plants that completely lacked anthers in an F2 population derived from an interspecific cross between diploid S. tuberosum and S. microdontum. We studied the antherless trait to determine its suitability for use in hybrid potato breeding. We mapped the causal locus to the short arm of Chromosome 6, developed KASP markers for the antherless (al) locus and introduced it into lines with T and A cytoplasm. We found that antherless type male sterility is not expressed in T and A cytoplasm, proving that it is a form of CMS. We hybridized male sterile al/al plants with P cytoplasm with pollen from al/al plants with T and A cytoplasm and we show that the resulting hybrids set significantly fewer berries in the field. Here, we show that the antherless CMS system can be readily deployed in diploid F1 hybrid potato breeding to improve hybridization efficiency and reduce berry set in the field.

Targeted genotyping-by-sequencing of potato and data analysis with R/polyBreedR

Mid-density targeted genotyping-by-sequencing (GBS) combines trait-specific markers with thousands of genomic markers at an attractive price for linkage mapping and genomic selection. A 2.5K targeted GBS assay for potato (Solanum tuberosum L.) was developed using the DArTag technology and later expanded to 4K targets. Genomic markers were selected from the potato Infinium single nucleotide polymorphism (SNP) array to maximize genome coverage and polymorphism rates. The DArTag and SNP array platforms produced equivalent dendrograms in a test set of 298 tetraploid samples, and 83% of the common markers showed good quantitative agreement, with RMSE (root mean squared error) <0.5. DArTag is suited for genomic selection candidates in the clonal evaluation trial, coupled with imputation to a higher density platform for the training population. Using the software polyBreedR, an R package for the manipulation and analysis of polyploid marker data, the RMSE for imputation by linkage analysis was 0.15 in a small half-diallel population (N = 85), which was significantly lower than the RMSE of 0.42 with the random forest method. Regarding high-value traits, the DArTag markers for resistance to potato virus Y, golden cyst nematode, and potato wart appeared to track their targets successfully, as did multi-allelic markers for maturity and tuber shape. In summary, the potato DArTag assay is a transformative and publicly available technology for potato breeding and genetics.

Genomics and resequencing of Fagopyrum dibotrys from different geographic regions reveals species evolution and genetic diversity

Fagopyrum dibotrys, belonging to the family Polygonaceae and genus Fagopyrum, is used in traditional Chinese medicine and is rich in beneficial components, such as flavonoids. As its abundant medicinal value has become increasingly recognized, its excessive development poses a considerable challenge to wild germplasm resources, necessitating artificial cultivation and domestication. Considering these factors, a high-quality genome of F. dibotrys was assembled and the evolutionary relationships within Caryophyllales were compared, based on which 58 individual samples of F. dibotrys were re-sequenced. We found that the samples could be categorized into three purebred populations and regions distributed at distinct elevations. Our varieties were cultivated from the parental populations of the subpopulation in central Yunnan. F. dibotrys is speculated to have originated in the high-altitude Tibetan Plateau region, and that its combination with flavonoids can protect plants against ultraviolet radiation; this infers a subpopulation with a high accumulation of flavonoids. This study assembled a high-quality genome and provided a theoretical foundation for the future introduction, domestication, and development of cultivated varieties of F. dibotrys.

Evidence-based unification of potato gene models with the UniTato collaborative genome browser

Potato (Solanum tuberosum) is the most popular tuber crop and a model organism. A variety of gene models for potato exist, and despite frequent updates, they are not unified. This hinders the comparison of gene models across versions, limits the ability to reuse experimental data without significant re-analysis, and leads to missing or wrongly annotated genes. Here, we unify the recent potato double monoploid v4 and v6 gene models by developing an automated merging protocol, resulting in a Unified poTato genome model (UniTato). We subsequently established an Apollo genome browser (unitato.nib.si) that enables public access to UniTato and further community-based curation. We demonstrate how the UniTato resource can help resolve problems with missing or misplaced genes and can be used to update or consolidate a wider set of gene models or genome information. The automated protocol, genome annotation files, and a comprehensive translation table are provided at github.com/NIB-SI/unitato.

Identification of the potato (Solanum tuberosum L.) P-type ATPase gene family and investigating the role of PHA2 in response to Pep13

P-type ATPase family members play important roles in plant growth and development and are involved in plant resistance to various biotic and abiotic factors. Extensive studies have been conducted on the P-type ATPase gene families in Arabidopsis thaliana and rice but our understanding in potato remains relatively limited. Therefore, this study aimed to screen and analyze 48 P-type ATPase genes from the potato (Solanum tuberosum L.) genome database at the genome-wide level. Potato P-type ATPase genes were categorized into five subgroups based on the phylogenetic classification of the reported species. Additionally, several bioinformatic analyses, including gene structure analysis, chromosomal position analysis, and identification of conserved motifs and promoter cis-acting elements, were performed. Interestingly, the plasma membrane H+-ATPase (PM H+-ATPase) genes of one of the P3 subgroups showed differential expression in different tissues of potato. Specifically, PHA2, PHA3, and PHA7 were highly expressed in the roots, whereas PHA8 was expressed in potatoes only under stress. Furthermore, the small peptide Pep13 inhibited the expression of PHA1, PHA2, PHA3, and PHA7 in potato roots. Transgenic plants heterologously overexpressing PHA2 displayed a growth phenotype sensitive to Pep13 compared with wild-type plants. Further analysis revealed that reducing potato PM H+-ATPase enzyme activity enhanced resistance to Pep13, indicating the involvement of PM H+-ATPase in the physiological process of potato late blight and the enhancement of plant disease resistance. This study confirms the critical role of potato PHA2 in resistance to Pep13.

A de novo genome assembly of Solanum bulbocastanum Dun., a Mexican diploid species reproductively isolated from the A-genome species, including cultivated potatoes

Potato and its wild relatives are distributed mainly in the Mexican highlands and central Andes of South America. The South American A-genome species, including cultivated potatoes, are reproductively isolated from Mexican diploid species. Whole-genome sequencing has disclosed genome structure and similarity, mostly in cultivated potatoes and their closely related species. In this study, we generated a chromosome-scale assembly of the genome of a Mexican diploid species, Solanum bulbocastanum Dun., using PacBio long-read sequencing, optical mapping, and Hi-C scaffolding technologies. The final sequence assembly consisted of 737.9 Mb, among which 647.0 Mb were anchored to the 12 chromosomes. Compared with chromosome-scale assemblies of S. lycopersicum (tomato), S. etuberosum (non-tuber-bearing species with E-genome), S. verrucosum, S. chacoense, S. multidissectum, and S. phureja (all four are A-genome species), the S. bulbocastnum genome was the shortest. It contained fewer transposable elements (56.2%) than A-genome species. A cluster analysis was performed based on pairwise ratios of syntenic regions among the seven chromosome-scale assemblies, showing that the A-genome species were first clustered as a distinct group. Then, this group was clustered with S. bulbocastanum. Sequence similarity in 1,624 single-copy orthologous gene groups among 36 Solanum species and clones separated S. bulbocastanum as a specific group, including other Mexican diploid species, from the A-genome species. Therefore, the S. bulbocastanum genome differs in genome structure and gene sequences from the A-genome species. These findings provide important insights into understanding and utilizing the genetic diversity of S. bulbocastanum and the other Mexican diploid species in potato breeding.

Construction of a Full-Length Transcriptome of Western Honeybee Midgut Tissue and Improved Genome Annotation

Honeybees are an indispensable pollinator in nature with pivotal ecological, economic, and scientific value. However, a full-length transcriptome for Apis mellifera, assembled with the advanced third-generation nanopore sequencing technology, has yet to be reported. Here, nanopore sequencing of the midgut tissues of uninoculated and Nosema ceranae-inoculated A. mellifera workers was conducted, and the full-length transcriptome was then constructed and annotated based on high-quality long reads. Next followed improvement of sequences and annotations of the current reference genome of A. mellifera. A total of 5,942,745 and 6,664,923 raw reads were produced from midguts of workers at 7 days post-inoculation (dpi) with N. ceranae and 10 dpi, while 7,100,161 and 6,506,665 raw reads were generated from the midguts of corresponding uninoculated workers. After strict quality control, 6,928,170, 6,353,066, 5,745,048, and 6,416,987 clean reads were obtained, with a length distribution ranging from 1 kb to 10 kb. Additionally, 16,824, 17,708, 15,744, and 18,246 full-length transcripts were respectively detected, including 28,019 nonredundant ones. Among these, 43,666, 30,945, 41,771, 26,442, and 24,532 full-length transcripts could be annotated to the Nr, KOG, eggNOG, GO, and KEGG databases, respectively. Additionally, 501 novel genes (20,326 novel transcripts) were identified for the first time, among which 401 (20,255), 193 (13,365), 414 (19,186), 228 (12,093), and 202 (11,703) were respectively annotated to each of the aforementioned five databases. The expression and sequences of three randomly selected novel transcripts were confirmed by RT-PCR and Sanger sequencing. The 5' UTR of 2082 genes, the 3' UTR of 2029 genes, and both the 5' and 3' UTRs of 730 genes were extended. Moreover, 17,345 SSRs, 14,789 complete ORFs, 1224 long non-coding RNAs (lncRNAs), and 650 transcription factors (TFs) from 37 families were detected. Findings from this work not only refine the annotation of the A. mellifera reference genome, but also provide a valuable resource and basis for relevant molecular and -omics studies.

The genomic and enzymatic basis for iridoid biosynthesis in cat thyme (Teucrium marum)

Iridoids are non-canonical monoterpenoids produced by both insects and plants. An example is the cat-attracting and insect-repelling volatile iridoid nepetalactone, produced by Nepeta sp. (catmint) and aphids. Recently, both nepetalactone biosynthetic pathways were elucidated, showing a remarkable convergent evolution. The iridoid, dolichodial, produced by Teucrium marum (cat thyme) and multiple insect species, has highly similar properties to nepetalactone but its biosynthetic origin remains unknown. We set out to determine the genomic, enzymatic, and evolutionary basis of iridoid biosynthesis in T. marum. First, we generated a de novo chromosome-scale genome assembly for T. marum using Oxford Nanopore Technologies long reads and proximity-by-ligation Hi-C reads. The 610.3 Mb assembly spans 15 pseudomolecules with a 32.9 Mb N50 scaffold size. This enabled identification of iridoid biosynthetic genes, whose roles were verified via activity assays. Phylogenomic analysis revealed that the evolutionary history of T. marum iridoid synthase, the iridoid scaffold-forming enzyme, is not orthologous to typical iridoid synthases but is derived from its conserved paralog. We discovered an enzymatic route from nepetalactol to diverse iridoids through the coupled activity of an iridoid oxidase cytochrome P450 and acetyltransferases, via an inferred acylated intermediate. This work provides a genomic resource for specialized metabolite research in mints and demonstration of the role of acetylation in T. marum iridoid diversity. This work will enable future biocatalytic or biosynthetic production of potent insect repellents, as well as comparative studies into iridoid biosynthesis in insects.

A chromosome-level genome of mango exclusively from long-read sequence data

Improvements in long-read sequencing techniques have greatly accelerated plant genome sequencing. Current de novo assemblies are routinely achieved by assembling long-read sequence data into contigs that are assembled to chromosome level by chromatin conformation capture. We report here a chromosome-level mango genome using only PacBio high-fidelity (HiFi) long reads. HiFi reads at high coverage (204x) resulted in the assembly of 17 chromosomes, each as a single contig with telomeres at both ends. The remaining three chromosomes were represented each by two contigs, with telomeres at one end and ribosomal repeats at the other end. Analyzing contig ends allowed them to be paired and linked to generate the remaining three complete chromosomes, telomere-to-telomere but with ribosomal repeats of uncertain length. The assembled genome was 365 Mb with 100% completeness as assessed by Benchmarking Universal Single-Copy Orthologs analysis. The haplotypes assembled demonstrated extensive structural differences. This approach using very high genome coverage may be useful for assembling high-quality genomes for many other plants.

The phased Solanum okadae genome and Petota pangenome analysis of 23 other potato wild relatives and hybrids

Potato is an important crop in the genus Solanum section Petota. Potatoes are susceptible to multiple abiotic and biotic stresses and have undergone constant improvement through breeding programs worldwide. Introgression of wild relatives from section Petota with potato is used as a strategy to enhance the diversity of potato germplasm. The current dataset contributes a phased genome assembly for diploid S. okadae, and short read sequences and de novo assemblies for the genomes of 16 additional wild diploid species in section Petota that were noted for stress resistance and were of interest to potato breeders. Genome sequence data for three additional genomes representing polyploid hybrids with cultivated potato, and an additional genome from non-tuberizing S. etuberosum, which is outside of section Petota, were also included. High quality short reads assemblies were achieved with genome sizes ranging from 575 to 795 Mbp and annotations were performed utilizing transcriptome sequence data. Genomes were compared for presence/absence of genes and phylogenetic analyses were carried out using plastome and nuclear sequences.

Integrated analysis of ATAC-seq and transcriptomic reveals the ScDof3-ScproC molecular module regulating the cold acclimation capacity of potato

Low temperature severely affects the geographical distribution and production of potato, which may incur cold damage in early spring or winter. Cultivated potatoes, mainly derived from Solanum tuberosum, are sensitive to freezing stress, but wild species of potato such as S. commersonii exhibit both constitutive freezing tolerance and/or cold acclimation tolerance. Hence, such wild species could assist in cold hardiness breeding. Yet the key transcription factors and their downstream functional genes that confer freezing tolerance are far from clear, hindering the breeding process. Here, we used ATAC-seq (Assay for Transposase-Accessible Chromatin with high-throughput sequencing) alongside RNA-seq to investigate the variation in chromatin accessibility and patterns of gene expression in freezing-tolerant CMM5 (S. commersonii), before and after its cold treatment. Our results suggest that after exposure to cold, transcription factors including Dof3, ABF2, PIF4, and MYB4 were predicted to further control the genes active in the synthetic/metabolic pathways of plant hormones, namely abscisic acid, polyamine, and reductive glutathione (among others). This suggests these transcription factors could regulate freezing tolerance of CMM5 leaves. In particular, ScDof3 was proven to regulate the expression of ScproC (pyrroline-5-carboxylate reductase, P5CR) according to dual-LUC assays. Overexpressing ScDof3 in Nicotiana benthamiana leaves led to an increase in both the proline content and expression level of NbproC (homolog of ScproC). These results demonstrate the ScDof3-ScproC module regulates the proline content and thus promotes freezing tolerance in potato. Our research provides valuable genetic resources to further study the molecular mechanisms underpinning cold tolerance in potato.

Molecular dissection of an intronic enhancer governing cold-induced expression of the vacuolar invertase gene in potato

Potato (Solanum tuberosum) is the third most important food crop in the world. Potato tubers must be stored at cold temperatures to minimize sprouting and losses due to disease. However, cold temperatures strongly induce the expression of the potato vacuolar invertase gene (VInv) and cause reducing sugar accumulation. This process, referred to as "cold-induced sweetening," is a major postharvest problem for the potato industry. We discovered that the cold-induced expression of VInv is controlled by a 200 bp enhancer, VInvIn2En, located in its second intron. We identified several DNA motifs in VInvIn2En that bind transcription factors involved in the plant cold stress response. Mutation of these DNA motifs abolished VInvIn2En function as a transcriptional enhancer. We developed VInvIn2En deletion lines in both diploid and tetraploid potato using clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease 9 (Cas9)-mediated gene editing. VInv transcription in cold-stored tubers was significantly reduced in the deletion lines. Interestingly, the VInvIn2En sequence is highly conserved among distantly related Solanum species, including tomato (Solanum lycopersicum) and other non-tuber-bearing species. We conclude that the VInv gene and the VInvIn2En enhancer have adopted distinct roles in the cold stress response in tubers of tuber-bearing Solanum species.

Tomato root specialized metabolites evolved through gene duplication and regulatory divergence within a biosynthetic gene cluster

Tremendous plant metabolic diversity arises from phylogenetically restricted specialized metabolic pathways. Specialized metabolites are synthesized in dedicated cells or tissues, with pathway genes sometimes colocalizing in biosynthetic gene clusters (BGCs). However, the mechanisms by which spatial expression patterns arise and the role of BGCs in pathway evolution remain underappreciated. In this study, we investigated the mechanisms driving acylsugar evolution in the Solanaceae. Previously thought to be restricted to glandular trichomes, acylsugars were recently found in cultivated tomato roots. We demonstrated that acylsugars in cultivated tomato roots and trichomes have different sugar cores, identified root-enriched paralogs of trichome acylsugar pathway genes, and characterized a key paralog required for root acylsugar biosynthesis, SlASAT1-LIKE (SlASAT1-L), which is nested within a previously reported trichome acylsugar BGC. Last, we provided evidence that ASAT1-L arose through duplication of its paralog, ASAT1, and was trichome-expressed before acquiring root-specific expression in the Solanum genus. Our results illuminate the genomic context and molecular mechanisms underpinning metabolic diversity in plants.

A coordinated switch in sucrose and callose metabolism enables enhanced symplastic unloading in potato tubers

One of the early changes upon tuber induction is the switch from apoplastic to symplastic unloading. Whether and how this change in unloading mode contributes to sink strength has remained unclear. In addition, developing tubers also change from energy to storage-based sucrose metabolism. Here, we investigated the coordination between changes in unloading mode and sucrose metabolism and their relative role in tuber sink strength by looking into callose and sucrose metabolism gene expression combined with a model of apoplastic and symplastic unloading. Gene expression analysis suggests that callose deposition in tubers is decreased by lower callose synthase expression. Furthermore, changes in callose and sucrose metabolism are strongly correlated, indicating a well-coordinated developmental switch. Modelling indicates that symplastic unloading is not the most efficient unloading mode per se. Instead, it is the concurrent metabolic switch that provides the physiological conditions necessary to potentiate symplastic transport and thereby enhance tuber sink strength .

StCoExpNet: a global co-expression network analysis facilitates identifying genes underlying agronomic traits in potatoes

KEY MESSAGE

We constructed a gene expression atlas and co-expression network for potatoes and identified several novel genes associated with various agronomic traits. This resource will accelerate potato genetics and genomics research. Potato (Solanum tuberosum L.) is the world's most crucial non-cereal food crop and ranks third in food production after wheat and rice. Despite the availability of several potato transcriptome datasets at public databases like NCBI SRA, an effort has yet to be put into developing a global transcriptome atlas and a co-expression network for potatoes. The objectives of our study were to construct a global expression atlas for potatoes using publicly available transcriptome datasets, identify housekeeping and tissue-specific genes, construct a global co-expression network and identify co-expression clusters, investigate the transcriptional complexity of genes involved in various essential biological processes related to agronomic traits, and provide a web server (StCoExpNet) to easily access the newly constructed expression atlas and co-expression network to investigate the expression and co-expression of genes of interest. In this study, we used data from 2299 publicly available potato transcriptome samples obtained from 15 different tissues to construct a global transcriptome atlas. We found that roughly 87% of the annotated genes exhibited detectable expression in at least one sample. Among these, we identified 281 genes with consistent and stable expression levels, indicating their role as housekeeping genes. Conversely, 308 genes exhibited marked tissue-specific expression patterns. We exemplarily linked some co-expression clusters to important agronomic traits of potatoes, such as self-incompatibility, anthocyanin biosynthesis, tuberization, and defense responses against multiple pathogens. The dataset compiled here constitutes a new resource (StCoExpNet), which can be accessed at https://stcoexpnet.julius-kuehn.de . This transcriptome atlas and the co-expression network will accelerate potato genetics and genomics research.

Characterization of FBA genes in potato (Solanum tuberosum L.) and expression patterns in response to light spectrum and abiotic stress

Fructose-1, 6-bisphosphate aldolase (FBA) plays vital roles in plant growth, development, and response to abiotic stress. However, genome-wide identification and structural characterization of the potato (Solanum tuberosum L.) FBA gene family has not been systematically analyzed. In this study, we identified nine StFBA gene members in potato, with six StFBA genes localized in the chloroplast and three in the cytoplasm. The analysis of gene structures, protein structures, and phylogenetic relationships indicated that StFBA genes were divided into Class I and II, which exhibited significant differences in structure and function. Synteny analysis revealed that segmental duplication events promoted the expansion of the StFBA gene family. Promoter analysis showed that most StFBA genes contained cis-regulatory elements associated with light and stress responses. Expression analysis showed that StFBA3, StFBA8, and StFBA9 showing significantly higher expression levels in leaf, stolon, and tuber under blue light, indicating that these genes may improve photosynthesis and play an important function in regulating the induction and expansion of microtubers. Expression levels of the StFBA genes were influenced by drought and salt stress, indicating that they played important roles in abiotic stress. This work offers a theoretical foundation for in-depth understanding of the evolution and function of StFBA genes, as well as providing the basis for the genetic improvement of potatoes.

Genotype-by-environment interactions and local adaptation shape selection in the US National Chip Processing Trial

KEY MESSAGE

We find evidence of selection for local adaptation and extensive genotype-by-environment interaction in the potato National Chip Processing Trial (NCPT). We present a novel method for dissecting the interplay between selection, local adaptation and environmental response in plant breeding schemes. Balancing local adaptation and the desire for widely adapted cultivars is challenging for plant breeders and makes genotype-by-environment interactions (GxE) an important target of selection. Selecting for GxE requires plant breeders to evaluate plants across multiple environments. One way breeders have accomplished this is to test advanced materials across many locations. Public potato breeders test advanced breeding material in the National Chip Processing Trial (NCPT), a public-private partnership where breeders from ten institutions submit advanced chip lines to be evaluated in up to ten locations across the country. These clones are genotyped and phenotyped for important agronomic traits. We used these data to interrogate the NCPT for GxE. Further, because breeders submitting clones to the NCPT select in a relatively small geographic range for the first 3 years of selection, we examined these data for evidence of incidental selection for local adaptation, and the alleles underlying it, using an environmental genome-wide association study (envGWAS). We found genomic regions associated with continuous environmental variables and discrete breeding programs, as well as regions of the genome potentially underlying GxE for yield.

A k-mer-based bulked segregant analysis approach to map seed traits in unphased heterozygous potato genomes

High-throughput sequencing-based methods for bulked segregant analysis (BSA) allow for the rapid identification of genetic markers associated with traits of interest. BSA studies have successfully identified qualitative (binary) and quantitative trait loci (QTLs) using QTL mapping. However, most require population structures that fit the models available and a reference genome. Instead, high-throughput short-read sequencing can be combined with BSA of k-mers (BSA-k-mer) to map traits that appear refractory to standard approaches. This method can be applied to any organism and is particularly useful for species with genomes diverged from the closest sequenced genome. It is also instrumental when dealing with highly heterozygous and potentially polyploid genomes without phased haplotype assemblies and for which a single haplotype can control a trait. Finally, it is flexible in terms of population structure. Here, we apply the BSA-k-mer method for the rapid identification of candidate regions related to seed spot and seed size in diploid potato. Using a mixture of F1 and F2 individuals from a cross between 2 highly heterozygous parents, candidate sequences were identified for each trait using the BSA-k-mer approach. Using parental reads, we were able to determine the parental origin of the loci. Finally, we mapped the identified k-mers to a closely related potato genome to validate the method and determine the genomic loci underlying these sequences. The location identified for the seed spot matches with previously identified loci associated with pigmentation in potato. The loci associated with seed size are novel. Both loci are relevant in future breeding toward true seeds in potato.

Chromosome-scale assembly of the streamlined picoeukaryote Picochlorum sp. SENEW3 genome reveals Rabl-like chromatin structure and potential for C4 photosynthesis

Genome sequencing and assembly of the photosynthetic picoeukaryotic Picochlorum sp. SENEW3 revealed a compact genome with a reduced gene set, few repetitive sequences, and an organized Rabl-like chromatin structure. Hi-C chromosome conformation capture revealed evidence of possible chromosomal translocations, as well as putative centromere locations. Maintenance of a relatively few selenoproteins, as compared to similarly sized marine picoprasinophytes Mamiellales, and broad halotolerance compared to others in Trebouxiophyceae, suggests evolutionary adaptation to variable salinity environments. Such adaptation may have driven size and genome minimization and have been enabled by the retention of a high number of membrane transporters. Identification of required pathway genes for both CAM and C4 photosynthetic carbon fixation, known to exist in the marine mamiellale pico-prasinophytes and seaweed Ulva, but few other chlorophyte species, further highlights the unique adaptations of this robust alga. This high-quality assembly provides a significant advance in the resources available for genomic investigations of this and other photosynthetic picoeukaryotes.

Identification and expression analyses of B3 genes reveal lineage-specific evolution and potential roles of REM genes in pepper

BACKGROUND

The B3 gene family, one of the largest plant-specific transcription factors, plays important roles in plant growth, seed development, and hormones. However, the B3 gene family, especially the REM subfamily, has not been systematically and functionally studied.

RESULTS

In this study, we performed genome-wide re-annotation of B3 genes in five Solanaceae plants, Arabidopsis thaliana, and Oryza sativa, and finally predicted 1,039 B3 genes, including 231 (22.2%) newly annotated genes. We found a striking abundance of REM genes in pepper species (Capsicum annuum, Capsicum baccatum, and Capsicum chinense). Comparative motif analysis revealed that REM and other subfamilies (ABI3/VP1, ARF, RAV, and HSI) consist of different amino acids. We verified that the large number of REM genes in pepper were included in the specific subgroup (G8) through the phylogenetic analysis. Chromosome location and evolutionary analyses suggested that the G8 subgroup genes evolved mainly via a pepper-specific recent tandem duplication on chromosomes 1 and 3 after speciation between pepper and other Solanaceae. RNA-seq analyses suggested the potential functions of REM genes under salt, heat, cold, and mannitol stress conditions in pepper (C. annuum).

CONCLUSIONS

Our study provides evolutionary and functional insights into the REM gene family in pepper.

Identification of two mutant JASON-RELATED genes associated with unreduced pollen production in potato

KEY MESSAGE

Multiple QTLs control unreduced pollen production in potato. Two major-effect QTLs co-locate with mutant alleles of genes with homology to AtJAS, a known regulator of meiotic spindle orientation. In diploid potato the production of unreduced gametes with a diploid (2n) rather than a haploid (n) number of chromosomes has been widely reported. Besides their evolutionary important role in sexual polyploidisation, unreduced gametes also have a practical value for potato breeding as a bridge between diploid and tetraploid germplasm. Although early articles argued for a monogenic recessive inheritance, the genetic basis of unreduced pollen production in potato has remained elusive. Here, three diploid full-sib populations were genotyped with an amplicon sequencing approach and phenotyped for unreduced pollen production across two growing seasons. We identified two minor-effect and three major-effect QTLs regulating this trait. The two QTLs with the largest effect displayed a recessive inheritance and an additive interaction. Both QTLs co-localised with genes encoding for putative AtJAS homologs, a key regulator of meiosis II spindle orientation in Arabidopsis thaliana. The function of these candidate genes is consistent with the cytological phenotype of mis-oriented metaphase II plates observed in the parental clones. The alleles associated with elevated levels of unreduced pollen showed deleterious mutation events: an exonic transposon insert causing a premature stop, and an amino acid change within a highly conserved domain. Taken together, our findings shed light on the natural variation underlying unreduced pollen production in potato and will facilitate interploidy breeding by enabling marker-assisted selection for this trait.

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.

Plant necrotrophic bacterial disease resistance phenotypes, QTL, and metabolites identified through integrated genetic mapping and metabolomics in Solanum species

Most food crops are susceptible to necrotrophic bacteria that cause rotting and wilting diseases in fleshy organs and foods. All varieties of cultivated potato (Solanum tuberosum L.) are susceptible to diseases caused by Pectobacterium species, but resistance has been demonstrated in wild potato relatives including S. chacoense. Previous studies demonstrated that resistance is in part mediated by antivirulence activity of phytochemicals in stems and tubers. Little is known about the genetic basis of antivirulence traits, and the potential for inheritance and introgression into cultivated potato is unclear. Here, the metabolites and genetic loci associated with antivirulence traits in S. chacoense were elucidated by screening a sequenced S. tuberosum x S. chacoense recombinant inbred line (RIL) population for antivirulence traits of its metabolite extracts. Metabolite extracts from the RILs exhibited a quantitative distribution for two antivirulence traits that were positively correlated: quorum sensing inhibition and exo-protease inhibition, with some evidence of transgressive segregation, supporting the role of multiple loci and metabolites regulating these resistance-associated systems. Metabolomics was performed on the highly resistant and susceptible RILs that revealed 30 metabolites associated with resistance, including several alkaloids and terpenes. Specifically, several prenylated metabolites were more abundant in resistant RILs. We constructed a high-density linkage map with 795 SNPs mapped to 12 linkage groups, spanning a length of 1,507 cM and a density of 1 marker per 1.89 cM. Genetic mapping of the antivirulence and metabolite data identified five quantitative trait loci (QTLs) related to quorum sensing inhibition that explained 8-28% of the phenotypic variation and two QTLs for protease activity inhibition that explained 14-19% of the phenotypic variation. Several candidate genes including alkaloid, and secondary metabolite biosynthesis that are related to disease resistance were identified within these QTLs. Taken together, these data support that quorum sensing inhibition and exo-protease inhibition assays may serve as breeding targets to improve resistance to nectrotrophic bacterial pathogens in potato and other plants. The identified candidate genes and metabolites can be utilized in marker assisted selection and genomic selection to improve soft- rot and blackleg disease resistance.

Whole-genome sequencing of tetraploid potato varieties reveals different strategies for drought tolerance

Climate changes leading to increasingly longer seasonal drought periods in large parts of the world increase the necessity for breeding drought-tolerant crops. Cultivated potato (Solanum tuberosum), the third most important vegetable crop worldwide, is regarded as drought-sensitive due to its shallow root architecture. Two German tetraploid potato cultivars differing in drought tolerance and their F1-progeny were evaluated under various drought scenarios. Bulked segregant analyses were combined with whole-genome sequencing (BSA-Seq) using contrasting bulks of drought-tolerant and drought-sensitive F1-clones. Applying QTLseqr, 15 QTLs comprising 588,983 single nucleotide polymorphisms (SNPs) in 2325 genes associated with drought stress tolerance were identified. SeqSNP analyses in an association panel of 34 mostly starch potato varieties using 1-8 SNPs for each of 188 selected genes narrowed the number of candidate genes down to 10. In addition, ent-kaurene synthase B was the only gene present under QTL 10. Eight of the identified genes (StABP1, StBRI1, StKS, StLEA, StPKSP1, StPKSP2, StYAB5, and StZOG1) address plant development, the other three genes (StFATA, StHGD and StSYP) contribute to plant protection under drought stress. Allelic variation in these genes might be explored in future breeding for drought-tolerant potato varieties.

Reproducible Quantitative Trait Loci for Resistance to Soft Rot Caused by Dickeya dianthicola Derived from the Wild Potato Solanum microdontum (PI 458355) Are Located on Chromosomes 1, 3, and 5

The potato wild relative Solanum microdontum is a breeder-friendly source of genetic resistance to soft rot. Our objectives were to (i) identify loci associated with soft rot resistance in S. microdontum germplasm and (ii) develop bi-parental populations in a self-compatible S. tuberosum genetic background to recover segregating F2 progenies, construct a linkage map, and identify quantitative trait loci (QTLs). Under objective (i), tubers from 103 S. microdontum genotypes from the United States Potato Genebank were inoculated with a high virulence strain of Dickeya dianthicola, and lesion size was measured after a 24-h incubation period at 30°C. Association analysis using 3,490 polymorphic Infinium array SNP markers identified soft rot resistance loci on chromosomes 1, 2, 3, 5, 7, 8, 11, and 12. Under objective (ii), a resistant S. microdontum accession PI 458355 was crossed with a highly fertile, self-compatible, diploid S. tuberosum pollen parent (PI 654351) to generate segregating F2 populations. Composite interval mapping was conducted using a genetic linkage map with 970 GBS-based SNP markers. Reproducible QTLs were detected on chromosomes 1, 3, and 5, explaining 11, 13, and 23% of the phenotypic variation, respectively. Homozygous S. microdontum alleles at the QTL on chromosome 5 and heterozygous or homozygous S. microdontum alleles at QTLs on chromosomes 1 and 3 significantly decrease lesion size compared with the homozygous S. tuberosum parent. The germplasm created in these studies provides a resource for studying traits from S. microdontum, and we can use the advanced F2 selections for future potato improvement. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Genome-wide genetic architecture for plant maturity and drought tolerance in diploid potatoes

Cultivated potato (Solanum tuberosum) is known to be highly susceptible to drought. With climate change and its frequent episodes of drought, potato growers will face increased challenges to achieving their yield goals. Currently, a high proportion of untapped potato germplasm remains within the diploid potato relatives, and the genetic architecture of the drought tolerance and maturity traits of diploid potatoes is still unknown. As such, a panel of 384 ethyl methanesulfonate-mutagenized diploid potato clones were evaluated for drought tolerance and plant maturity under field conditions. Genome-wide association studies (GWAS) were conducted to dissect the genetic architecture of the traits. The results obtained from the genetic structure analysis of the panel showed five main groups and seven subgroups. Using the Genome Association and Prediction Integrated Tool-mixed linear model GWAS statistical model, 34 and 17 significant quantitative trait nucleotides (QTNs) were found associated with maturity and drought traits, respectively. Chromosome 5 carried most of the QTNs, some of which were also detected by using the restricted two-stage multi-locus multi-allele-GWAS haploblock-based model, and two QTNs were found to be pleiotropic for both maturity and drought traits. Using the non-parametric U-test, one and three QTNs, with 5.13%-7.4% phenotypic variations explained, showed favorable allelic effects that increase the maturity and drought trait values. The quantitaive trait loci (QTLs)/QTNs associated with maturity and drought trait were found co-located in narrow (0.5-1 kb) genomic regions with 56 candidate genes playing roles in plant development and senescence and in abiotic stress responses. A total of 127 potato clones were found to be late maturing and tolerant to drought, while nine were early to moderate-late maturing and tolerant to drought. Taken together, the data show that the studied germplasm panel and the identified candidate genes are prime genetic resources for breeders and biologists in conventional breeding and targeted gene editing as climate adaptation tools.

Identification and Characterization of the BZR Transcription Factor Genes Family in Potato (Solanum tuberosum L.) and Their Expression Profiles in Response to Abiotic Stresses

Brassinazole resistant (BZR) genes act downstream of the brassinosteroid signaling pathway regulating plant growth and development and participating in plant stress responses. However, the BZR gene family has not systematically been characterized in potato. We identified eight BZR genes in Solanum tuberosum, which were distributed among seven chromosomes unequally and were classified into three subgroups. Potato and tomato BZR proteins were shown to be closely related with high levels of similarity. The BZR gene family members in each subgroup contained similar conserved motifs. StBZR genes exhibited tissue-specific expression patterns, suggesting their functional differentiation during evolution. StBZR4, StBZR7, and StBZR8 were highly expressed under white light in microtubers. StBZR1 showed a progressive up-regulation from 0 to 6 h and a progressive down-regulation from 6 to 24 h after drought and salt stress. StBZR1, StBZR2, StBZR4, StBZR5, StBZR6, StBZR7 and StBZR8 were significantly induced from 0 to 3 h under BR treatment. This implied StBZR genes are involved in phytohormone and stress response signaling pathways. Our results provide a theoretical basis for understanding the functional mechanisms of BZR genes in potato.

Phylogenetic and Expression Analysis of the Sucrose Synthase and Sucrose Phosphate Synthase Gene Family in Potatoes

Sucrose synthase (SUS) and sucrose phosphate synthase (SPS) are essential in plant sucrose metabolism. The potato is an important crop worldwide, but systematic analyses of the StSUS and StSPS gene families in potatoes are still lacking. Ten sucrose metabolism-related genes were identified in this study. The SUSs and SPSs could each be split into three subgroups through phylogenetic analysis. StSUSIc was the most highly expressed gene in different developmental tissues. Ka/Ks analysis showed that StSUSIb and StSUSIc were subjected to more-significant homozygous selection pressure. Our cis-acting element analysis of the StSUS and StSPS promoter sequences showed four elements: defense- and stress-responsive, hormone-responsive, light-responsive, and transcription factor elements. The expression of StSUS and StSPS genes was found to be regulated by circadian rhythm. In the treatments of 1% to 5% sucrose, glucose, and fructose, the expression of StSUS and StSPS family genes was enhanced by sucrose, but inhibited at high-glucose and fructose concentrations. This study identified six StSUS and four StSPS genes and analyzed their gene structure, conserved motifs, chromosome position, promoter elements, phylogenetic tree, and tissue-specific expression patterns. Our results will motivate more research into the biological process underlying the genes of sucrose metabolism in potatoes.

Haplotype-resolved assembly of a tetraploid potato genome using long reads and low-depth offspring data

Potato is one of the world's major staple crops, and like many important crop plants, it has a polyploid genome. Polyploid haplotype assembly poses a major computational challenge. We introduce a novel strategy for the assembly of polyploid genomes and present an assembly of the autotetraploid potato cultivar Altus. Our method uses low-depth sequencing data from an offspring population to achieve chromosomal clustering and haplotype phasing on the assembly graph. Our approach generates high-quality assemblies of individual chromosomes with haplotype-specific sequence resolution of whole chromosome arms and can be applied in common breeding scenarios where collections of offspring are available.

Crossover shortage in potato is caused by StMSH4 mutant alleles and leads to either highly uniform unreduced pollen or sterility

The balanced segregation of homologous chromosomes during meiosis is essential for fertility and is mediated by crossovers (COs). A strong reduction of CO number leads to the unpairing of homologous chromosomes after the withdrawal of the synaptonemal complex. This results in the random segregation of univalents during meiosis I and ultimately to the production of unbalanced and sterile gametes. However, if CO shortage is combined with another meiotic alteration that restitutes the first meiotic division, then uniform and balanced unreduced male gametes, essentially composed of nonrecombinant homologs, are produced. This mitosis-like division is of interest to breeders because it transmits most of the parental heterozygosity to the gametes. In potato, CO shortage, a recessive trait previously referred to as desynapsis, was tentatively mapped to chromosome 8. In this article, we have fine-mapped the position of the CO shortage locus and identified StMSH4, an essential component of the class I CO pathway, as the most likely candidate gene. A 7 base-pair insertion in the second exon of StMSH4 was found to be associated with CO shortage in our mapping population. We also identified a second allele with a 3,820 base-pair insertion and confirmed that both alleles cannot complement each other. Such nonfunctional alleles appear to be common in potato cultivars. More than half of the varieties we tested are carriers of mutational load at the StMSH4 locus. With this new information, breeders can choose to remove alleles associated with CO shortage from their germplasm to improve fertility or to use them to produce highly uniform unreduced male gametes in alternative breeding schemes.

Chromosome-scale genome assembly of the 'Munstead' cultivar of Lavandula angustifolia

OBJECTIVES

Lavandula angustifolia (English lavender) is commercially important not only as an ornamental species but also as a major source of fragrances. To better understand the genomic basis of chemical diversity in lavender, we sequenced, assembled, and annotated the 'Munstead' cultivar of L. angustifolia.

DATA DESCRIPTION

A total of 80 Gb of Oxford Nanopore Technologies reads was used to assemble the 'Munstead' genome using the Canu genome assembler software. Following multiple rounds of error correction and scaffolding using Hi-C data, the final chromosome-scale assembly represents 795,075,733 bp across 25 chromosomes with an N50 scaffold length of 31,371,815 bp. Benchmarking Universal Single Copy Orthologs analysis revealed 98.0% complete orthologs, indicative of a high-quality assembly representative of genic space. Annotation of protein-coding sequences revealed 58,702 high-confidence genes encoding 88,528 gene models. Access to the 'Munstead' genome will permit comparative analyses within and among lavender accessions and provides a pivotal species for comparative analyses within Lamiaceae.

Karyotyping of aneuploid and polyploid plants from low coverage whole-genome resequencing

BACKGROUND

Karyotype, as a basic characteristic of species, provides valuable information for fundamental theoretical research and germplasm resource innovation. However, traditional karyotyping techniques, including fluorescence in situ hybridization (FISH), are challenging and low in efficiency, especially when karyotyping aneuploid and polyploid plants. The use of low coverage whole-genome resequencing (lcWGR) data for karyotyping was explored, but existing methods are complicated and require control samples.

RESULTS

In this study, a new protocol for molecular karyotype analysis was provided, which proved to be a simpler, faster, and more accurate method, requiring no control. Notably, our method not only provided the copy number of each chromosome of an individual but also an accurate evaluation of the genomic contribution from its parents. Moreover, we verified the method through FISH and published resequencing data.

CONCLUSIONS

This method is of great significance for species evolution analysis, chromosome engineering, crop improvement, and breeding.

StHAB1, a negative regulatory factor in abscisic acid signaling, plays crucial roles in potato drought tolerance and shoot branching

Abscisic acid (ABA) is critical in drought tolerance and plant growth. Group A protein type 2C phosphatases (PP2Cs) are negative regulators of ABA signaling and plant adaptation to stress. Knowledge about the functions of potato group A PP2Cs is limited. Here, we report that the potato group A PP2C StHAB1 is broadly expressed in potato plants and strongly induced by ABA and drought. Suppression of StHAB1 enhanced potato ABA sensitivity and drought tolerance, whereas overexpression of the dominant mutant StHAB1G276D compromised ABA sensitivity and drought tolerance. StHAB1 interacts with almost all ABA receptors and the Snf1-Related Kinase OST1. Suppressing StHAB1 and overexpressing StHAB1G276D alter potato growth morphology; notably, overexpression of StHAB1G276D causes excessive shoot branching. RNA-sequencing analyses identified that the auxin efflux carrier genes StPIN3, StPIN5, and StPIN8 were up-regulated in StHAB1G276D-overexpressing axillary buds. Correspondingly, the auxin concentration was reduced in StHAB1G276D-overexpressing axillary buds, consistent with the role of auxin in repressing lateral branch outgrowth. The expression of BRANCHED1s (StBRC1a and StBRC1b) was unchanged in StHAB1G276D-overexpressing axillary buds, suggesting that StHAB1G276D overexpression does not cause axillary bud outgrowth via regulation of BRC1 expression. Our findings demonstrate that StHAB1 is vital in potato drought tolerance and shoot branching.

Genome-wide identification of 2-oxoglutarate and Fe (II)-dependent dioxygenase family genes and their expression profiling under drought and salt stress in potato

The 2-Oxoglutatrate-dependent dioxygenases (2OGDs) comprise the 2-Oxoglutatrate and Fe(II)-dependent dioxygenases (2ODD) enzyme families that facilitate the biosynthesis of various compounds like gibberellin, ethylene, etc. The 2OGDs are also involved in various catabolism pathways, such as auxin and salicylic acid catabolism. Despite their important roles, 2ODDs have not been studied in potato, which is the third most important crop globally. In this study, a comprehensive genome wide analysis was done to identify all 2ODDs in potatoes, and the putative genes were analysed for the presence of the signature 2OG-FeII_Oxy (PF03171) domain and the conserved DIOX_N (PF14226) domain. A total of 205 St2ODDs were identified and classified into eight groups based on their function. The physiochemical properties, gene structures, and motifs were analysed, and gene duplication events were also searched for St2ODDs. The active amino acid residues responsible for binding with 2-oxoglutarate and Fe (II) were conserved throughout the St2ODDs. The three-dimensional (3D) structures of the representative members of flavanol synthase (FNS), 1-aminocyclopropane-1-carboxylic acid oxidases (ACOs), and gibberellin oxidases (GAOXs) were made and docked with their respective substrates, and the potential interactions were visualised. The expression patterns of the St2ODDs under abiotic stressors such as heat, salt, and drought were also analysed. We found altered expression levels of St2ODDs under abiotic stress conditions, which was further confirmed for drought and salt stress using qRT-PCR. The expression levels of St2ODD115, St2ODD34, and St2ODD99 were found to be upregulated in drought stress with 2.2, 1.8, and 2.6 fold changes, respectively. After rewatering, the expression levels were normal. In salt stress, the expression levels of St2ODD151, St2ODD76, St2ODD91, and St2ODD34 were found to be upregulated after 24 hours (h), 48 hours (h), 72 hours (h), and 96 hours (h). Altogether, the elevated expression levels suggest the importance of St2ODDs under abiotic stresses, i.e., drought and salt. Overall, our study provided a knowledge base for the 2ODD gene family in potato, which can be used further to study the important roles of 2ODDs in potato plants.

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

Haplotype-resolved genome assembly of Populus tremula × P. alba reveals aspen-specific megabase satellite DNA

Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717-1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence-sensitive analyses. Here we report a telomere-to-telomere genome for this hybrid aspen with two chromosome-scale, haplotype-resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array-based and array-less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus) but PtaM147-like sequences occur in LTR-retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis-driven research to probe into the function of the abundant and aspen-specific PtaM147 satellite DNA.

SMRT-AgRenSeq-d in potato (Solanum tuberosum) as a method to identify candidates for the nematode resistance Gpa5

Potato is the third most important food crop in the world. Diverse pathogens threaten sustainable crop production but can be controlled, in many cases, through the deployment of disease resistance genes belonging to the family of nucleotide-binding, leucine-rich-repeat (NLR) genes. To identify effective disease resistance genes in established varieties, we have successfully established SMRT-AgRenSeq in tetraploid potatoes and have further enhanced the methodology by including dRenSeq in an approach that we term SMR-AgRenSeq-d. The inclusion of dRenSeq enables the filtering of candidates after the association analysis by establishing a presence/absence matrix across resistant and susceptible varieties that is translated into an F1 score. Using a SMRT-RenSeq-based sequence representation of the NLRome from the cultivar Innovator, SMRT-AgRenSeq-d analyses reliably identified the late blight resistance benchmark genes Rpi-R1, Rpi-R2-like, Rpi-R3a, and Rpi-R3b in a panel of 117 varieties with variable phenotype penetrations. All benchmark genes were identified with an F1 score of 1, which indicates absolute linkage in the panel. This method also identified nine strong candidates for Gpa5 that controls the potato cyst nematode (PCN) species Globodera pallida (pathotypes Pa2/3). Assuming that NLRs are involved in controlling many types of resistances, SMRT-AgRenSeq-d can readily be applied to diverse crops and pathogen systems.

MFPPDB: a comprehensive multi-functional plant peptide database

Plants produce a wide range of bioactive peptides as part of their innate defense mechanisms. With the explosive growth of plant-derived peptides, verifying the therapeutic function using traditional experimental methods are resources and time consuming. Therefore, it is necessary to predict the therapeutic function of plant-derived peptides more effectively and accurately with reduced waste of resources and thus expedite the development of plant peptides. We herein developed a repository of plant peptides predicted to have multiple therapeutic functions, named as MFPPDB (multi-functional plant peptide database). MFPPDB including 1,482,409 single or multiple functional plant origin therapeutic peptides derived from 121 fundamental plant species. The functional categories of these therapeutic peptides include 41 different features such as anti-bacterial, anti-fungal, anti-HIV, anti-viral, and anti-cancer. The detailed physicochemical information of these peptides was presented in functional search and physicochemical property search module, which can help users easily access the peptide information by the plant peptide species, ID, and functions, or by their peptide ID, isoelectric point, peptide sequence, and molecular weight through web-friendly interface. We further matched the predicted peptides to nine state-of-the-art curated functional peptide databases and found that at least 293,408 of the peptides possess functional potentials. Overall, MFPPDB integrated a massive number of plant peptides have single or multiple therapeutic functions, which will facilitate the comprehensive research in plant peptidomics. MFPPDB can be freely accessed through http://124.223.195.214:9188/mfppdb/index.

Genome-wide association mapping and genomic prediction for late blight and potato cyst nematode resistance in potato (Solanum tuberosum L.)

Potatoes are an important source of food for millions of people worldwide. Biotic stresses, notably late blight and potato cyst nematodes (PCN) pose a major threat to potato production worldwide, and knowledge of genes controlling these traits is limited. A genome-wide association mapping study was conducted to identify the genomic regulators controlling these biotic stresses, and the genomic prediction accuracy was worked out using the GBLUP model of genomic selection (GS) in a panel of 222 diverse potato accessions. The phenotype data on resistance to late blight and two PCN species (Globodera pallida and G. rostochiensis) were recorded for three and two consecutive years, respectively. The potato panel was genotyped using genotyping by sequencing (GBS), and 1,20,622 SNP markers were identified. A total of 7 SNP associations for late blight resistance, 9 and 11 for G. pallida and G. rostochiensis, respectively, were detected by additive and simplex dominance models of GWAS. The associated SNPs were distributed across the chromosomes, but most of the associations were found on chromosomes 5, 10 and 11, which have been earlier reported as the hotspots of disease-resistance genes. The GS prediction accuracy estimates were low to moderate for resistance to G. pallida (0.04-0.14) and G. rostochiensis (0.14-0.21), while late blight resistance showed a high prediction accuracy of 0.42-0.51. This study provides information on the complex genetic nature of these biotic stress traits in potatoes and putative SNP markers for resistance breeding.

Deciphering the underlying immune network of the potato defense response inhibition by Phytophthora infestans nuclear effector Pi07586 through transcriptome analysis

Phytophthora infestans, a highly destructive plant oomycete pathogen, is responsible for causing late blight in potatoes worldwide. To successfully infect host cells and evade immunity, P. infestans secretes various effectors into host cells and exclusively targets the host nucleus. However, the precise mechanisms by which these effectors manipulate host gene expression and reprogram defenses remain poorly understood. In this study, we focused on a nuclear-targeted effector, Pi07586, which has been implicated in immune suppression. Quantitative real-time PCR (qRT-PCR) analysis showed Pi07586 was significant up-regulation during the early stages of infection. Agrobacterium-induced transient expression revealed that Pi07586 localized in the nucleus of leaf cells. Overexpression of Pi07586 resulted in increased leaf colonization by P. infestans. RNA-seq analysis revealed that Pi07586 effectively suppressed the expression of PR-1C-like and photosynthetic antenna protein genes. Furthermore, high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS) analysis indicated that Pi07586 overexpression led to a substantial decrease in abscisic acid (ABA), jasmonic acid (JA), and jasmonoyl-isoleucine (JA-Ile) levels, while not affecting salicylic acid (SA) and indole-3-acetic acid (IAA) production. These findings shed new light on the modulation of plant immunity by Pi07586 and enhance our understanding of the intricate relationship between P. infestans and host plants.